the emergence of classical dynamical chaos from more profound quantum mechanics is one of the most intriguing problems in physics @xcite .
dynamical chaos in classical mechanics is a special kind of random motion in dynamical systems without any noise and random parameters that is characterized by sensitive dependence on initial conditions in a bounded phase space .
the current consensus is that _ isolated bounded quantum systems _ do not show sensitive dependence on initial conditions in the same way as classical systems because their evolution is unitary .
so the question is : what is the fundamental mechanism of arising classical chaos from quantum mechanics ?
the purpose of this lecture is to show how classical hamiltonian chaos with sensitive dependence on initial conditions , positive lyapunov exponents , fractal properties of underlying phase space , and anomalous statistical characteristics may arise from quantum electrodynamics of a single atom strongly interacting with a quantized mode of the electromagnetic radiation field in a high - quality cavity .
the study of the fundamental atom - photon interaction constitutes the rapidly growing field of cavity quantum electrodynamics ( for a review see @xcite ) .
the experimental state of the art has reached in this field the stage where the transition from classical to quantum dynamical regimes can now be probed directly .
atoms and photons , confined in a high - quality cavity , are ideal objects to study quantum - classical correspondence and quantum chaos .
real quantum systems are not isolated .
they interact with their environment and in attempts to measure their states with classical measuring devices which should , by virtue of their purpose , be in unstable states .
so we are dealing with quantum - classical hybrids .
it is not a simple question which degrees of freedom should be considered in a given physical situation as quantum ones and which ones as classical .
let us consider an excited atom in a single - mode high - quality cavity whose frequency is close to the frequency of one of the atomic electrodipole transitions .
the atom emits a photon in the cavity mode and goes to a lower - lying state . under the conditions of the strong atom - field coupling
, the photon may be reemitted and reabsorbed by the atom many times .
sooner or later , the atom will emit a photon in one of the other modes of the electromagnetic field that are not sustained by the cavity ( for example , transverse modes in an open fabry - perot cavity ) and the photon will be lost in a surrounding , the process known as spontaneous emission . in real experiments , the cavity relaxation due to the losses in the cavity walls can not be neglected .
it is usually modeled by a weak coupling of the selected cavity mode with a bath of harmonic oscillators ( spanning a wide frequency range ) in thermal equilibrium at a given temperature .
sooner or later , the atom will relax to a lower - lying state that is not resonant with the cavity mode .
it seems to be reasonable to treat two near - resonant atomic levels , strongly coupled to a cavity mode , as a quantum dynamical system weakly coupled to the other electromagnetic modes , the other atomic levels and the cavity walls .
an infinite number of the respective degrees of freedom , which are coupled to the quantum dynamical system s degrees of freedom but are not affected by them , forms an inexhaustible external reservoir ( that may be treated classically ) usually called  the environment " .
it is a kind of the external coupling which is inevitably present in reality .
there is another kind of coupling called the dynamical coupling .
for example , if the average number of photons in the cavity mode is sufficiently large , one may treat the field as a classical object dynamically coupled to quantized atoms ( the semiclassical approximation ) and take into account the feedback effect of the atoms on the radiation field .
the semiclassical approximation breaks down when we deal with single cold atoms and photons in a cavity .
in this situation we should take into account the translational ( external ) atomic degree of freedom .
when a cold atom emits and absorbs photons its momentum and position may vary significally due to the recoil effect . the external atomic degree of freedom , which may be treated classically if the values of the atomic momentum are greater than the photon momentum , is dynamically coupled to the internal atomic and field degrees of freedom which are treated quantum mechanically .
we show in this lecture that , even in the absence of any other interaction with the environment , the interaction of the quantum atom - field system with the external atomic degree of freedom provides , under appropriate conditions , the emergence of classical chaos from quantum electrodynamics .
the discovery that a single - mode laser , a symbol of coherence and stability , may exhibit deterministic chaos is especially important not only because lasers are one of the main instruments in physics but lasers , as well , provide almost ideal systems to test general ideas in quantum mechanics and statistical physics . from the standpoint of nonlinear dynamics , laser is an open dissipative system that transforms an external excitation into a coherent output in the presence of loss .
some manifestations of strange attractors and dissipative chaos have been observed with different types of lasers ( for a review , see @xcite ) .
a collection of identical two - level atoms , interacting with a single - mode electromagnetic field , provides the simplest model for laser dynamics . because of a large number of atoms ( and photons )
, laser dynamics can be adequately described in the so - called semiclassical approximation , where one treats atoms as two - level objects , which may be described by the bloch equations for two components of the collective electrodipole polarisation and collective atomic population inversion , interacting with a field mode to be governed by the classical maxwell equations for the field strengths with the right hands depending on the atomic polarisation . with a single - mode
homogeneously broadened laser , operating in resonance with the gain center , these five equations can be reduced to three real - valued equations for slowly varying amplitudes which have been shown to be equivalent to the well - known lorentz model for fluid convection @xcite .
practically in the same time , ideas of dynamical chaos have been explored with fundamental models of the matter - radiation interaction , comprising of a collection of two - level atoms _ interacting with their own radiation field _ in a perfect single - mode cavity without any loss and external excitation .
it has been shown theoretically and numerically @xcite that the semiclassical maxwell
bloch equations , following from the dicke hamiltonian @xcite , may demonstrate _
hamiltonian semiclassical chaos _
if one goes beyond the so - called rotating - wave approximation , i.  e. if one takes into account energy non - conserving terms in the hamiltonian ( for details , see the beginning of the next section ) .
this mechanism of arising hamiltonian semiclassical chaos is rather weak under realistic assumptions . on the other hand , rotating - wave approximation suppresses chaos with atoms at rest due to existence of an additional integral of motion , the conserving interaction energy .
it has been shown in a series of our papers @xcite that hamiltonian chaos may arise within the rotating - wave approximation under conditions of a modulation of the atom - field interaction . in a natural way it occurs
when atoms move through a cavity in a direction along which the cavity sustains a standing - wave field that is periodic in space @xcite .
when moving with a constant velocity , atoms `` see '' the field whose strength changes in time periodically .
it breaks down the interaction - energy integral and may cause hamiltonian semiclassical chaos of a homoclinic type @xcite .
there are another ways of modulations to be considered in @xcite .
the modulation of the detuning between the atomic transition frequency and the frequency of the field mode has been shown to produce parametric instability and hamiltonian semiclassical chaos @xcite .
structural hamiltonian chaos @xcite may arise as a result of a harmonic modulation of cavity length which causes the respective oscillations of the nodes of a standing - wave .
generally speaking , the model of the atom - field interaction should involve not only the internal atomic and field degrees of freedom but also the center - of - mass motion of the atom . when emitting and absorbing photons , atoms not only change their internal states but their position and momentum are changed as well due the photon recoil effect .
this effect may be neglected if one deals with rydberg atoms interacting with a microwave field in a cavity ( as it has been done in @xcite ) or with thermal usual atoms interacting with a visible light .
it has been theoretically and numerically shown in @xcite that hamiltonian semiclassical chaos may arise with cold atoms with recoil in a standing - wave microcavity .    in trying to describe adequately dynamics of single atoms and photons in a high - quality cavity
, one should go beyond the semiclassical approximation and treats the atom - photon interaction on a quantum ground .
the fully quantum model of the interaction between a single two - level atom ( without a recoil ) with a single - mode quantized field in an ideal cavity is known as the jaynes - cummings model @xcite .
it describes the atom - field system as a quantum - electrodynamical object whose evolution in time is ( quasi)periodic . in schrdinger picture , it may be described by an infinite set of linear ordinary differential equations for the probability amplitudes to find the atom in the ground / excited state and the field in the state with @xmath0 photons , where @xmath0 runs from zero to infinity . when adopting the semiclassical approximation , we reduce , by hook or by crook , this infinite set to a small number of equations for atomic and field variables .
in fact , we decouple the quantum atom - field system into an atom and field parts which may exchange excitations with each other . in doing so
, we get automatically products of the atom and field variables in the equations of motion .
the respective maxwell
bloch equations may , under appropriate conditions , produce chaos in the classical sense of sensitive dependence on initial conditions in the reduced classical phase space spanned by the atomic and field expectation values .
it resembles the procedure of deriving the famous lorenz equations from an infinite hierarchy of mode equations in fluid convection by reducing it to three main modes only .
the reduction of an infinite set of linear equations , that are not chaotic in the classical sense , to a finite set of nonlinear equations , that may be chaotic , makes an impression that it is only a mathematical trick .
it seems that it is not only a useful trick enabling us to handle with the equations of motion but the reduction is a model of processes that may occur in nature . coherence loss caused by inevitable interaction with environment ( decoherence ) breaks quantum unitarity suppressing some quantum properties of motion and revealing its classical properties .
such a situation can be modeled with the use of a quantum - classical hybrid , a system with quantum and classical degrees of freedom dynamically coupled to each other . which part of the whole system under consideration is quantized and which one is treated classically depends on the physical situation .
if we have a large number of atoms ( or photons ) we may adopt with a good accuracy the semiclassical approximation treating the cavity field as a classical one .
when considering a single atom with a recoil in a high - quality microcavity it is necessary to quantize the internal electronic atomic and field degrees of freedom but the translational atomic motion may be treated classically if the values of the atomic momentum are greater than the photon momentum .
the interaction between matter and radiation , commonly evidenced by spontaneous emission and induced emission and absorption of photons by atoms , is one of the most fundamental dynamical interactions in nature . in free space
, an excited atomic state decays irreversibly because an infinity of vacuum electromagnetic states is available to the radiated photons which disappears in the free - space mode continuum .
the situation may be changes if the mode structure and the density of modes are modified by placing the radiating atom into a high - quality cavity .
a cavity - induced enhancement and inhibition of spontaneous emission has been observed in a number of experiments in different ranges of the electromagnetic field , from microwaves to visible light ( for a review see @xcite ) .    in a near - resonant ideal cavity with a single - mode frequency closed to the atomic transitional frequency , the radiation emitted by the atom is reflected at the cavity walls and is reabsorbed by the atom many times before it dissipates . in
particularly , spontaneous emission becomes reversible . in cavity quantum electrodynamics
near - resonant interaction of a single atom with a single - mode cavity field is commonly modeled by the hamiltonian @xmath1 where @xmath2 represents the free motionless two - level atom with @xmath3 being the energy separation between the electronic states , the ground @xmath4 and excited @xmath5 ones .
the pauli spin operator acts on the states as @xmath6 the hamiltonian @xmath7 represents the mode of the isolated field .
the photon annihilation @xmath8 and creation @xmath9 operators with the commutation rule @xmath10 = 1 $ ] act on the photon number states as follows : @xmath11 the electric dipole interaction between the atom and the field mode is described by the operator @xmath12 where @xmath13 is the atom - field coupling constant at the atomic position @xmath14 in the one - dimensional cavity , which is known under the names `` vacuum or single - photon rabi frequency '' .
the pauli spin operators @xmath15 with the commutation rules @xmath16=\mp 2\hat \sigma_\pm$ ] , @xmath17=\hat \sigma_z$ ] describe transitions between the ground and excited electronic states @xmath18 the first form of the hamiltonian ( [ 6 ] ) follows from the quantization of the classical electric dipole atom - field interaction by replacing the atomic dipole moment by the operator @xmath19 and the classical electric field by the operator @xmath20 .
the products @xmath21 and @xmath22 describe the transition of the atom from the ground ( excited ) to the excited ( ground ) state and simultaneous annihilating ( creating ) a photon in the field mode , respectively .
they are usual energy - conserving processes .
the product @xmath23 describes the transition of the atom from the ground to the excited state and simultaneous creating a photon in the mode .
the reverse process of simultaneous atom and field de - excitation is described by the term @xmath24 .
the last two processes are not energy - conserving and may be neglected while we deal with a single atom .
omitting the energy - nonconserving terms defines the jaynes - cummings model @xcite which was proposed originally in 1963 as a purely theoretical tool for studying fundamentals of the interaction between matter and radiation .
it seemed to be far from reality because neither the atom nor the cavity mode can hardly be perfectly isolated from environment with those - days technique .
exciting progress in experimental techniques has drastically changed the situation .
experiments , using very high - quality microwave cavities with @xmath25 and optical microcavities with @xmath26 , have now achieved the exceptional circumstance of strong coupling between atoms and a cavity field with the strength of the coupling exceeding the atomic and cavity decays that provides manipulations with single atoms and photons and experimental proving of some predictions of the jaynes - cummings model @xcite .
the art of experimentalists and theoretical efforts have opened exciting perspectives for realizing gedanken experiments on fundamentals of quantum mechanics and implementing quantum communications and computing . as we will show in this lecture , an atom , interacting with a radiation field in a high-@xmath27 cavity ,
provides an example of strongly - coupled microscopic nonlinear system whose dynamics may be very complicated and even chaotic .
the jaynes - cummings hamiltonian has the form @xmath28 where @xmath29 is an amplitude value of the vacuum rabi frequency and @xmath30 is a shape function of the cavity mode .
let us expand the atom - field state vector over the electronic atomic states and fock ( or photon - number ) field states @xmath31 @xmath32 where @xmath33 are probability amplitudes to find the atom in the excited ( ground ) state and @xmath0 photons in the field mode , respectively .
substitution of eqs.([8 ] ) and ( [ 9 ] ) into the schrdinger equation @xmath34 results in an infinite - dimensional set of coupled linear equations for the complex - valued probability amplitudes @xmath35,&\\ \dot b_{n+1}&=\phantom{- } i\left[\delta_b b_{n+1}^ *   + \sqrt{n+1}\,f(x ) a_n^ * \right],&\\ \end{aligned } \label{11}\ ] ] where dot denotes differentiation with respect to dimensionless time @xmath36 , @xmath37 , and @xmath38/\omega_0 $ ] .
it is convenient to define new real - valued combinations of the probability amplitudes @xmath39 for which we get from ( [ 11 ] ) the quantum bloch - like equations of motion @xmath40 where @xmath41 is the dimensionless detuning between the atomic transition , @xmath42 , and the field - mode , @xmath43 , frequencies . for each specified photon number @xmath0 ,
the quantity @xmath44 is conserved .
two global integrals of motion @xmath45 reflect conservation of the total energy and of the total probability , respectively .
eqs.([13 ] ) can be easily solved for each @xmath0 in terms of trigonometric functions ( we assume that @xmath30 does not depend on time ) .
the general solution of the schrdinger equation is the sum of the solutions of an infinite set of the independent bloch - like equations .
the exact general solution for one of the measured quantities , the atomic population inversion @xmath46 is the following : @xmath47 where @xmath48 is the @xmath0-photon raby frequency , @xmath49 , @xmath50 and @xmath51 are the respective initial values which are determined by the initial atomic and field states . as it is seen from ( [ 16 ] ) , the solution with initially sharply defined atomic and field energy states ( @xmath4 or @xmath5 and @xmath31 ) describes a periodic exchange of one quantum of energy between the atom and the field .
the fock state @xmath31 with a specified number of photons in the mode is an exotic field state . in general
, a pure quantized field state is an infinite superposition of the photon - number states @xmath52 where @xmath53 is the probability for observing @xmath0 photons .
the most classical of single - mode quantum states is a coherent state @xmath54 whose photon probabilities follow a poisson distribution with the mean number of photons @xmath55 and the root mean square spread @xmath56
. evolution of the population inversion of an atom , @xmath57 , in the field , initially prepared to be in a coherent state , contains all the rabi frequencies @xmath58 .
the incommensurability of rabi frequencies for different @xmath0 inevitably washes out the periodicity of population transfer resulting in a collapse of the excited - state occupation probability .
since the frequencies @xmath58 form a    * ( a ) * + * ( b ) * +    discrete set , the signal reappears after a time ( see fig .  [ fig1]a computed with @xmath59 , @xmath60 and with the initial conditions ( [ 33 ] ) ) .
the collapse - revival phenomenon , a direct demonstration of quantum nature of the radiation field , has been demonstrated experimentally @xcite .
let us consider as well another pure field state @xmath61 whose photon probabilities follow a bose - einstein distribution @xmath62 collapses and revivals also occur with this field state ( see fig .
[ fig1]b computed with @xmath59 , @xmath60 and with the initial conditions ( [ 34 ] ) ) , but the collapse time is much shorter as compared with a coherent state because the bose - einstein spread in the photon number is far larger than for a coherent field with the same mean number of photons .
[ fig1 ] demonstrates that the simple jaynes - cummings model with @xmath59 may produce very complicated time evolution
. however , the respective oscillations , of course , are not chaotic in the classical sense of exponential sensitivity to small changes in initial conditions .
it is clear from the very structure of the equations of motion ( [ 13 ] ) which are infinite - dimensional but linear . in classical mechanics
the phase space is continuous and states of a classical system can be arbitrarily close to each other providing a possibility of exponential diverging of initially close trajectories in a bounded region of the phase space . in quantum mechanics
there is no notion of a trajectory , and the quantum phase space is not continuous due to the heisenberg uncertainty principle .
evolution of _ an isolated quantum system _ is unitary and no dynamical chaos in the classical sense is possible with such systems .
the notion `` quantum chaos '' just refers to the behavior of quantum systems whose classical counterparts behave chaotically @xcite .
real quantum systems are not isolated , they interact with their environment and , under attempts to measure their states , with classical devices which , by virtue of their purpose , should be in unstable states .
loosing of coherence , due to inevitable interaction with environment ( decoherence ) , breaks down quantum unitarity suppressing quantum properties of motion and manifesting classical ones . a quantum rotor with periodic kicks provides an illustrative example of suppressing quantum dynamical localization when one takes into account a thermal `` bath '' in the respective equations of motion @xcite .
in the next section we will elaborate these ideas with an extended version of the strongly coupled atom - field system considered above .
in the process of emitting and absorbing photons of the cavity - field mode , the atom not only changes its internal electronic states but its external translational state is changed as well due to the photon recoil effect . in this section
we consider a single two - level atom with mass @xmath63 moving in an ideal cavity which sustains a single standing - wave mode along the axis @xmath14 with the wave vector @xmath64 and the shape function @xmath65 .
the respective hamiltonian is the following extension of the jaynes - cummings hamiltonian ( [ 8 ] ) @xmath66 where the momentum @xmath67 and position @xmath68 operators satisfy the standard commutation relation @xmath69=i\hbar$ ] .
operators , belonging to different degrees of freedom , commute with each other at the same time moment .
we have now three degrees of freedom , the internal atomic and the field ones , and the external ( or translational ) atomic degree of freedom .
the first two degrees of freedom are treated as fully quantum ones in the schrdinger picture .
in fact , there are an infinite number of quantum degrees of freedom ( see eq .
( [ 13 ] ) ) which are entangled
. the external degree of freedom will be treated as the classical one that may be justified by large values of the atomic momentum as compared with the photon momentum @xmath70 .
the hamilton equations of motion for the classical external degree of freedom is easily found from the hamiltonian ( [ 22 ] ) @xmath71 where @xmath72 denotes an expectation value of the corresponding operator over a quantum state @xmath73 of the electronic - field hamiltonian . using the normalizations @xmath74 , @xmath75 , and @xmath76
, we get from ( [ 23 ] ) @xmath77 where @xmath78 is the normalized recoil frequency which characterizes the average change in kinetic energy of the atom , @xmath79 , in the process of emission and absorption of a photon . after computing the expectation value in ( [ 24 ] ) with the state vector ( [ 9 ] ) , we obtain our basic hamilton - schrdinger equations @xmath80 which describe the atom - field quantum - classical hybrid with the quantum degrees of freedom dynamically coupled to the classical degree of freedom
. this infinite set of _ nonlinear ordinary differential equations
_ possesses an infinite number of the integrals of motion , the total energy integral @xmath81 the bloch - like integral for each @xmath0 @xmath82 and the global integral @xmath83 the infinite - dimensional nonlinear dynamical system ( [ 25 ] ) is a quantum generalization of the five - dimensional semiclassical set of the equations of motion for the same problem that was derived in refs @xcite with a classical field .
the later one was shown @xcite to be chaotic with positive values of the maximal lyapunov exponent in some ranges of the system s control parameters .
it is rather difficult to compute this quantitative indicator of dynamical chaos with an infinite dimensional set of ode s .
however , we have found another signatures of hamiltonian chaos with eqs .
( [ 25 ] ) , the main of which are fractals which will be demonstrated in the next section .    *
( a ) * + * ( b ) * +    let us compare the rabi - oscillation signals with the atom - field quantum - classical hybrid ( [ 25 ] ) and with the fully quantum jaynes - cummings model without recoil ( [ 13 ] ) . in fig .
[ fig2]a and b the population inversion @xmath57 ( see eq .
( [ 15 ] ) ) for an initially excited atom interacting in resonance with the field that is initially in a coherent state ( [ 19 ] ) and in a bose - einstein state ( [ 20 ] ) , respectively , is shown .
this figure should be compared with fig .
[ fig1 ] where the rabi oscillations with the jaynes - cummings model have been computed for the same initial conditions and the field states with @xmath84 but with a motionless atom .
more pronounced collapses and revivals occur with the quantum - classical hybrid as compared with the fully quantum model without recoil . in order to understand peculiarities of the rabi oscillations in fig .
[ fig2 ] look at fig .
[ fig7]a demonstrating the scheme of gedanken experiments .
an atom starts at @xmath85 and moves to the right with the initial momentum @xmath86 .
it reaches the first node of the standing wave , where the coupling coefficient with the field mode is zero at time moment @xmath87 in dimensionless units ( we choose @xmath88 in our computer simulations ) .
for both the field states , one can see the respective slowing down of the rabi oscillations at the moments of time @xmath89 , when the atom transverses the @xmath90-th node .
the pronounced peaks of revivals of the rabi oscillations occur at the moments @xmath91 , when the atom transverses the @xmath92-th antinode of the standing wave where its coupling with the field is maximal .    in spite of a complicated character of the rabi oscillations , shown in fig .
[ fig2 ] , they are regular because at exact resonance , @xmath93 , all the values of the real - valued amplitudes @xmath94 are conserved during the evolution . with initially excited atom we have @xmath95 for each @xmath0
. it immediately follows from eqs.([25 ] ) that the atom moves with a constant velocity through the cavity .
it means that @xmath14 is a linear function of time , and the hamilton - schrdinger equations([25 ] ) reduce to the periodically modulated linear bloch - like equations ( [ 13 ] ) with periodic solutions .
the exact solution for the atomic population inversion can be easily found @xmath96 it is a periodic function ( shown in fig .
[ fig2 ] ) with the period to be equal to @xmath97 and the maxima at @xmath98 @xmath99 .    out off resonance , @xmath100
, the quantum - classical hybrid may demonstrate chaos . to diagnose chaos
it is instructive to compute the maximal lyapunov exponent @xmath101 whose values depend on initial conditions , on the detuning @xmath102 , the mean number of photons in the mode , and on the recoil frequency @xmath103 .
if the quantized field is initially prepared in the fock state @xmath31 with exactly @xmath0 quanta in the mode , the infinite - dimensional set ( [ 25 ] ) reduces to 8 equations @xmath104 with initial conditions @xmath105 describing the atom - field system initially prepared in a superposition state @xmath106    the so - called topographic @xmath101-maps , showing by color modulation values of @xmath101 in the ranges of values of two control parameters with the third one to be fixed , have been computed with the set ( [ 29 ] ) . in fig .
[ fig3 ] we show the @xmath101-map with initial zero atomic population , @xmath107 , @xmath108 , in the ranges of the detuning @xmath102 and the recoil frequency @xmath103 with the fixed value of the initial number of photons @xmath109 . at exact resonance ,
the set ( [ 29 ] ) is integrable and @xmath110 at @xmath93 . in the range @xmath111
, that corresponds to realistic values of the recoil frequency , the maximal lyapunov exponent may be positive and one may expect chaotic atomic motion in the respective ranges of @xmath103 and @xmath102 .
another @xmath101-map with initially excited atom , @xmath112 , @xmath113 , demonstrates in fig .
[ fig4 ] the values of @xmath101 in dependence on @xmath103 and @xmath0 at @xmath114 in double logarithmic scale .
the magnitude of the maximal lyapunov exponent grows , in average , with increasing the initial number of photons in the cavity mode .    *
( a ) * + * ( b ) *    computing the values of @xmath101 is difficult with the cavity field initially prepared in a superposition state ( [ 18 ] ) ( including coherent and bose - einstein states ) that generates an infinite number of the equations of motion ( [ 25 ] ) .
we have found another signatures of chaos . in fig .
[ fig5 ] an atomic trajectory , @xmath115 , and oscillations of the atomic momentum , @xmath116 , are shown to illustrate manifestations of chaos in the set ( [ 25 ] ) with the initially coherent field and the excited atom with @xmath117 , @xmath118 , @xmath88 , and @xmath119 .
a typical weakly chaotic atomic trajectory represents a kind of a random walking with small oscillations of the atom in potential wells , long ballistic flights of the lvy type with almost constant velocity and erratic turnbacks .
just the presence of lvy flights fluenses strongly the statistical properties of the system generating power - like statistical laws that have been extensively studied in the semiclassical approximation in refs.@xcite .    *
( a ) * + * ( b ) *    fig .  [ fig6 ] presents poincar sections of motion in the set ( [ 25 ] ) with five different values of the initial atomic momentum to be projected on the plane of the atomic external variables @xmath120 .
[ fig6]a demonstrates an example of chaotic poincar section @xmath121 whereas fig .
[ fig7]b presents a regular motion under the same conditions except for the detuning @xmath114 .
a feasible scheme for detecting manifestation of chaos with hot two - level rydberg atoms moving in a high - q microwave cavity has been proposed in @xcite .
the same idea could be realized with cold usual atoms in a high - q microcavity .
consider a 2d - geometry of a gedanken experiment with a monokinetic atomic beam propagating almost perpendicularly to the cavity axis @xmath14 . in a reference frame
moving with a constant velocity in the perpendicular direction , there remains only the transverse atomic motion along the axis @xmath14 .
one measures atomic population inversion after passing the interaction zone . before injecting atoms in the cavity ,
it is necessary to prepare all the atoms in the same electronic state , say , in the excited state , with the help of a @xmath122-pulse of the laser radiation .
it may be done with _ only a finite accuracy _
, say , equal to @xmath123 for the initial population inversion @xmath124 .
the values of the population inversion @xmath125 are measured with detectors at a fixed time moment .
if we would work with the values of the control parameters corresponding to the regular atom - field dynamics , we would expect to have a regular curve @xmath125@xmath124 . in the chaotic regime ,
the atomic inversion at the output can be predicted ( within a certain confidence interval @xmath126 ) for a time not exceeding the so - called predictability horizon @xmath127 which depends weakly on @xmath128 and @xmath126 .
since the maximal confidence interval lies in the range @xmath129 and @xmath101 may reach the values of the order of @xmath130 ( see @xmath101-maps ) , the predictability horizon in accordance with the formula ( [ predtime ] ) can be very short : with @xmath101 = 0.5 the predictability horizon @xmath131 may be of the order of 10 in units reciprocal of the vacuum rabi frequency @xmath29 that corresponds to @xmath132 s with the realistic value of @xmath133 rad@xmath134 s@xmath135 .     on its initial values @xmath124 with @xmath118 ( a ) at @xmath136 and ( b ) at @xmath137 with the inset showing this dependence at the exact atom - field resonance , @xmath93.,scaledwidth=45.0% ]    in the regular regime , the inevitable errors in preparing @xmath128 produce the output errors @xmath138 of the same order .
in the chaotic regime , the initial uncertainty increases exponentially resulting in a complete uncertainty of the detected population inversion in a reasonable time . it is demonstrated in fig .
[ fig7 ] , where we plot the dependence of the values of @xmath139 at @xmath136 ( fig .
[ fig7]a ) and @xmath137 ( fig .
[ fig7]b ) on the values of @xmath140 in the chaotic regime with an initially fock field at @xmath118 and @xmath141 .
simulation shows that an initial error @xmath142 in preparing the atomic electronic state leads to complete uncertainty @xmath143 in a rather short time . to see the difference
, it is desirable to carry out a control experiment at the exact resonance ( @xmath93 ) when the atomic motion is fully regular with any initial values .
the dependence @xmath125@xmath124 with @xmath93 is demonstrated in the inset in fig .
[ fig7]b with all the other control parameters and initial values being the same .
in this section , we treat the atom - photon interaction in a high-@xmath27 cavity as a chaotic scattering problem .    * ( a ) * + * ( b ) *    let us consider the scheme of scattering of atoms by the standing wave shown in fig .
atoms , one by one , are placed at the point @xmath85 with different initial values of the momentum @xmath144 along the cavity axis . for simplicity
, we suppose that they have no momentum in the other directions ( 1d - geometry ) .
we compute the time the atoms need to reach one of the detectors placed at the cavity mirrors .
the dependence of this exit time @xmath145 on the initial atomic momentum @xmath144 is studied under the other initial conditions and parameters being the same . to avoid complications that are not essential to the main theme of this section , we consider the cavity with only two standing - wave lengths .      in this subsection
the field is supposed to be initially prepared in a fock state . before injecting into a cavity , atoms are suppose to be prepared in the superposition state with @xmath146 , @xmath147 , i.  e.  in the state with zero population inversion @xmath148 .    at exact resonance ( @xmath93 ) with @xmath149 ,
the optical potential @xmath150 is equal to zero , and the analytical expression for the dependence in question can be easily found to be the following : @xmath151 atoms simply fly through the cavity in one direction with their initial constant velocity and are registered by one of the detectors . out off resonance ( @xmath152 ) , the atomic motion has been numerically found in preceding section to be chaotic with positive values of the maximal lyapunov exponent in the following ranges of the values of the control parameters : the detuning @xmath153 and the recoil frequency @xmath154 .    * ( a ) * + * ( b ) * + * ( c ) *    fig .
[ fig9 ] shows the function @xmath155 with the normalized detuning @xmath118 , the recoil frequency @xmath156 , and the initial number of cavity photons @xmath157 .
the exit - time function demonstrates an intermittency of smooth curves and complicated structures that can not be resolved in principle , no matter how large the magnification factor .
[ fig9]b shows magnification of the function for the small interval @xmath158 .
further magnification in the range @xmath159 shown in fig .
[ fig9]c reveals a beautiful self - similar structure . some structures in fig .
[ fig9]a that look like fractal are not , in fact , unresolvable and self - similar .
magnification of the structure in the range @xmath160 demonstrates quite a smooth function without unresolvable substructures and with only two singular points on the borders of the respective momentum interval . beating in all the structures of the atomic fractal in fig .
[ fig9 ] should be attributed to the structure of the hamilton - schrdinger equations ( [ 29 ] ) which describe two atom - field oscillators with slightly different frequencies .
the exit time @xmath145 , corresponding to both smooth and unresolved @xmath144 intervals , increases in average with increasing the magnification factor .
it follows that there exist atoms never reaching the detectors in spite of the fact that they have no obvious energy restrictions to leave the cavity .
tiny interplay between chaotic external and internal dynamics prevents these atoms from leaving the cavity . the similar phenomenon in hamiltonian systems
is known as _ dynamical trapping _ @xcite .
different kinds of atomic trajectories , which are computed with the system ( [ 29 ] ) , are shown in fig .
a trajectory with the number @xmath161 transverses the central node of the standing - wave , before being detected , @xmath161 times and is called @xmath161-th trajectory .
there are also special separatrix - like @xmath162-trajectories following which atoms in infinite time reach the stationary points @xmath163 ( @xmath164 ) , @xmath165 , transversing @xmath161 times the central node .
these points are the anti - nodes of the standing wave where the force acting on atoms is zero .
a detuned atom can asymptotically reach one of the stationary points after transversing the central node @xmath161 times .
the trajectory with number 1 , showing in fig .
[ fig8]b , is close to a separatrix - like @xmath166-trajectory .
the smooth @xmath144 intervals in the first - order structure in fig .
[ fig9]a correspond to atoms transversing once the central node and reaching the right detector .
the unresolved singular points in the first - order structure with @xmath167 at the border between the smooth and unresolved @xmath144 intervals are generated by the @xmath166-trajectories .
analogously , the smooth @xmath144 intervals in the second - order structure in fig .
[ fig9]b correspond to the 2-nd order trajectories with singular points between them corresponding to the @xmath168-trajectories and so on .
there are two different mechanisms of generation of infinite exit times , namely , dynamical trapping with infinite oscillations ( @xmath169 ) in a cavity and the separatrix - like motion ( @xmath170 ) .
the set of all initial momenta generating the separatrix - like trajectories is a countable fractal .
each point in the set can be specified as a vector in a hilbert space with @xmath161 integer nonzero components .
one is able to prescribe to any unresolved interval of a @xmath161-th order structure a set with @xmath161 integers , where the first integer is a number of a second - order structure to which trajectory under consideration belongs in the first - order structure , the second integer is a number of a third - order structure in the second - order structure mentioned above , and so on .
such a number set is analogous to a directory tree address : `` @xmath171a subdirectory of the root directory@xmath172/@xmath171a subdirectory of the 2-nd level@xmath172/@xmath171a subdirectory of the 3-rd level@xmath172/ ... '' .
unlike the separatrix fractal , the set of all initial atomic momenta leading to dynamically trapped atoms with @xmath169 seems to be uncountable .    *
( a ) * + * ( b ) *    nonlinear fractal dynamics implies specific statistical properties of chaotic motion in hamiltonian systems ( for a review see @xcite ) .
even a finite - dimensional atom - field system s phase space has a complicated topology whose simplified image is given by projections of poincar sections and fractal scattering functions .
tiny interplay between all the degrees of freedom is responsible for dynamical trapping of atoms even in a very short microcavity and for anomalous statistical properties of the fundamental atom - photon interaction .
the probability distribution of exit times @xmath173 for @xmath174 events with initial atomic momenta in the range @xmath175 is shown in fig .
[ fig10]a in double logarithmic scale .
it is close to a poissonian distribution with comparatively short exit times up to @xmath176 and demonstrates an algebraic decay with the characteristic exponent @xmath177 at its tail ( @xmath178 ) . in the range of initial atomic momenta @xmath179 , the respective pdf for @xmath180 events
is computed under the same other conditions to be a poissonian - like .      *
( a ) * + * ( b ) * + * ( c ) *    following to the scheme in fig .
[ fig8]a , let us consider in this section scattering of atoms by the standing wave field initially prepared in the coherent state ( [ 19 ] ) .
before injecting into a cavity , atoms supposed to be prepared in the excited state @xmath5 with @xmath113 . in this case the initial conditions are @xmath181 fig .
[ fig11 ] demonstrates the respective exit - time function @xmath155 with @xmath182 , @xmath88 , and @xmath84 .
a coherent state is not a sharply defined state , as a fock one , but a superposition of an infinite number of fock states .
we used a truncated basis of @xmath183 fock states for the cavity mode in our simulation .
resolution of one of the unresolved structures in fig .
[ fig11]a is shown in fig .
[ fig11]b .
it resembles the respective fragment of the fock fractal in fig .
further magnification of the function in fig .
[ fig11]b , which is shown in fig .
[ fig11]c , reveals again a self - similar structure with smooth and singular zones .
we collect an exit time statistics with @xmath184 events by counting atoms with initial momenta in the range @xmath185 reaching the detectors . the plot of the respective histogram of exit times , shown in fig .
[ fig12 ] , demonstrates a few local maxima .
the corresponding times of exit of atoms can be approximately estimated with the help of the formula ( [ 33 ] ) for resonant atoms .
the first maximum around @xmath186 corresponds , mainly , to atoms which , being initially directed to the right with comparatively low momenta @xmath187 , turn back before reaching the central node and are registered by the left detector .
it follows from ( [ 33 ] ) that for such atoms @xmath188 at @xmath189 .
however , the atoms , to be injected initially with the momentum @xmath190 and registered by the right detector , may contribute to the first maximum as well because , after transversing the central node , they can be accelerated and gain the values of the momentum @xmath191 up to @xmath192 .
the second local maximum around @xmath193 corresponds , mainly , to atoms with @xmath194 which transverse the central node and are registered by the right detector for the time @xmath195 .
the other local maxima of the pdf are not so pronounced as the first ones , they are formed by the atoms transversing the central node a few times .
the pdf in fig .
[ fig12 ] demonstrates an exponential decay at the tail up to the exit times @xmath196 .      in conclusion of this section we present in fig .
13 the fractal computed with the field initially prepared in a bose - einstein state ( [ 20 ] ) and atoms supposed to be prepared in the excited state @xmath5 with @xmath113 .
the initial conditions are @xmath197 the values of the control parameters are the same as for the coherent fractal . fig .
13 demonstrates again a self - similar structure of the exit - time distribution .
we want to stress that the hamilton - schrdinger sets of equations generating the fock and the other fractals differ strongly in the number of equations . nevertheless , the respective @xmath155 functions are rather similar in their main features .
the atomic fractals we have found are generated by dynamical chaos in the atom - photon interaction but not by noise .
* ( a ) * + * ( b ) * + * ( c ) * +    [ fig13 ]
we have studied in the strong - coupling regime the fundamental interaction between a two - level atom with recoil and a quantized radiation field in a single - mode cavity in the mixed quantum - classical formalism modelling quantum evolution of the electronic - field purely quantum system to be disturbed by translational atomic motion .
we have managed to derive an infinite set of the hamilton - schrdinger nonlinear ode s for real valued electronic - field probability amplitudes and expectation values of the atomic center - of - mass position and momentum .
these equations of motion are able to demonstrate known features of the atom - field quantum evolution in the strong - coupling limit such as collapses and revivals of the atomic population inversion modulated by atomic motion through the nodes and antinodes of the standing wave .
we have shown that even in the absence of any other interaction with environment the hamiltion - schrdinger dynamical system provides the emergence of classical hamiltonian dynamical chaos from cavity quantum electrodynamics .
we have investigated the atom - photon nonlinear dynamics with fock , coherent and bose - einstein quantum states of the initial cavity field .
positive values of the maximal lyapunov exponent @xmath101 have been found with reasonable values of the control parameters , the detuning of the atom - field resonance , @xmath102 , the atomic recoil frequency , @xmath103 , and the initial mean number of photons .
exponential sensitivity to initial conditions may manifest itself in the chaotic dependence of the output values of the atomic population inversion on its input values that may , in principle , be measured in real experiments .
new manifestations of hamiltonian chaos in the atom - photon interaction , we have found numerically in the dependence of the atomic exit times on the initial momentum , are atomic fractals .
they were classifield by the names of the respective quantum field states , i.e. , the fock , coherent and bose - einstein fractals .
we have found anomalous statistical properties of the chaotic atom - photon interaction : lvy atomic flights through a cavity and a power law at the tail of the exit - time distribution function in the fock field .
power - law tails should appear ( at least , in some ranges of the control parameters ) in the pdf s with coherent and bose - einstein fields , as well , but their detection would require a longer computation time than we have used .
i am grateful to leonid konkov and michael uleysky for preparing the figures .
this work was supported by the russian foundation for basic research under grant no .
020217796 , by the program  mathematical methods in nonlinear dynamics " of the russian academy of sciences , and by the program for basic research of the far eastern division of the russian academy of sciences .
chirikov , phys .
52 ( 1979 ) 263 .
zaslavsky , phys .
80 ( 1981 ) 157 .
f. haake , _ quantum signatures of chaos _ , springer , berlin ( 1991 ) .
gutzwiller , _ chaos in classical and quantum mechanics _ , springer , new york ( 1990 ) . , ed . by p.r .
berman , academic press , new york ( 1994 ) .
h. walther , phys .
rep . 219 ( 1992 ) 263 .
j. m. raimond , m. brune , s. haroche , rev .
73 ( 2001 ) 565 . w. vogel , d .-
welsch , s. wallentowitz , _ quantum optics : an introduction _ , wiley - vch , berlin ( 2001 ) .
ackerhalt , p.w .
milonni , m .-
shih , phys .
128 ( 1985 ) 205 .
harrison , contemp .
29 ( 1988 ) 341 .
khanin , _ foundations of laser dynamics _ ,
nauka , moscow ( 1999 ) [ in russian ] .
h. haken , phys .
lett . a 53 ( 1975 ) 77 .
belobrov , g.m .
zaslavskii , g.kh .
tartakovskii , zh .
71 ( 1976 ) 1799 [ sov . phys .
jetp 44 ( 1976 ) 945 ] .
dicke , phys .
93 ( 1954 ) 493 .
prants , l.e .
konkov , pisma zh .
65 ( 1997 ) 801 [ jetp lett .
65 ( 1997 ) 833 ] .
prants , l.e .
konkov , phys . lett .
a 225 ( 1997 ) 33 .
prants , l.e .
konkov , zh . eksp . teor .
( 1999 ) 740 [ jetp 88 ( 1999 ) 406 ] .
prants , l.e .
konkov , i.l .
kirilyuk , phys .
e 60 ( 1999 ) 335 .
prants , l.e .
konkov , phys .
e 61 ( 2000 ) 3632 .
ioussoupov , l.e .
konkov , s.v .
prants , physica  d 155 ( 2001 ) 311 . s.v .
prants , l.e .
konkov , pisma zh .
73 ( 2001 ) 200 [ jetp lett .
73 ( 2001 ) 180 ] .
prants , v.yu .
sirotkin , phys .
a 64 ( 2001 ) 033412 .
jaynes , f.w .
cummings , proc .
ieee 51 ( 1963 ) 89 .
raimond , s. haroche , in : _ confined electrons and photons _ , ed .
by e.  burstein and c.  weisbuch , plenum press , new york , 1995 .
g. rempe , h. walther , n. klein , phys .
58 ( 1987 ) 353 .
prants , pisma zh . eksp .
( 2002 ) 777 [ jetp letters 75 ( 2002 ) 651 ] .
prants , m. edelman , g.m .
zaslavsky , phys .
e 66 ( 2002 ) 046222 .
argonov , s.v .
prants , zh .
( 2003 ) 946 .
[ jetp 96 ( 2003 ) 832 ] .
zaslavsky , _ physics of chaos in hamiltonian systems _ , academic press , oxford , 1998 .


during the past two decades much attention has been devoted to the study of disordered bulk metals with electron - electron interactions  @xcite from the one hand and arrays of normal metallic grains connected by tunnel junctions  @xcite from the other .
the coulomb interaction plays a very important role in both types of system and many non - trivial physical effects occur due to it . in disordered metals , _ e.g. _ , the coulomb interaction results in quantum corrections to conductivity and reduces the density of states , a phenomenon known as zero bias anomaly .
qualitatively similar effects were predicted  @xcite and studied experimentally  @xcite for arrays of tunnel junctions in the normal state , where a sufficiently strong coulomb interaction can suppress the conductivity ( coulomb blockade of tunneling ) .    in spite of the similarities in the physical behavior of disordered metals and normal arrays of tunnel junctions ,
completely different approaches are used for the theoretical description of these respective systems . as concerns disordered metals with interaction , either diagrammatic methods  @xcite or non - linear @xmath0-models @xcite are employed .
the starting point of these theories is a model of a weakly disordered bulk metal . on the other hand , for the description of coulomb blockade effects in a junction array or a granular material , a completely different approach based on the theory of ambegaokar , eckern and schn ( aes )  @xcite
is commonly used .
the free energy functional proposed by aes does not contain any disorder : the finite conductance obtained within this model is due to tunneling of electrons between grains and can be incorporated into the theory through a term describing the so - called quantum dissipation , first introduced phenomenologically in ref .  .
although both approaches have been introduced almost 20 years ago , only a few indirect attempts have been made to reconcile them .
we mention the work by nazarov  @xcite and levitov and shytov  @xcite that treat the diffusive zero bias anomaly and the coulomb blockade of tunneling on the same footing .
more recently , nazarov  @xcite showed how coulomb blockade phenomena survive in diffusive systems .
indeed , the properties of a granular metal without electron - electron interactions are qualitatively similar to those of a disordered bulk metal ( see e.g. ref .  ) and it would be natural to expect that this qualitative similarity persists even in the presence of interactions .
in fact , one may conjecture that bulk disordered metals from the one hand , and granular metals or arrays of tunnel junctions from the other , can be described within one unifying scheme .    in order to develop such a scheme , we consider in this paper a granular normal metal with coulomb interaction between electrons .
the consideration is simplified by taking into account only the long range part of the interaction that leads to coulomb blockade effects .
we assume that the macroscopic dimensionless conductance @xmath1 of the granular metal is large , which enables us to develop a perturbative theory with @xmath2 as the small parameter . within this theory ,
physical quantities of interest can be calculated at arbitrary temperatures .
the diagrammatic analysis is supplemented by the derivation of a @xmath3-model that can be considerably simplified provided one does not take into account weak localization effects .
we demonstrate below that granular metals and networks of tunnel junctions can always be described using the standard techniques developed in the theory of disordered metals ( expansions in cooperons and diffusons , @xmath3- model calculations ) . with the help of these techniques
one can in principle calculate any physical quantity without using the notion of quantum dissipation .
disorder is inevitably present in the system even if electrons move ballistically within the grains . in the latter case
any small irregularity in the shape of the grains causes the electron motion within the grains to be chaotic and this assumption is sufficient for our theory to be applicable .
an integrable shape of the grains does not seem to be realistic , and even if it happened for isolated grains , tunneling from grain to grain would add an additional chaoticity .
these ideas about intrinsic disorder in granular metals have been used in previous studies  @xcite .
we find that , at not too low a temperature , the aes free energy functional gives correct results .
however in the limit @xmath4 aes theory is no longer applicable .
the region of validity of the aes free energy functional is determined by the inequality @xmath5 , where @xmath6 is the mean level spacing , @xmath7 is the density of states at the fermi surface and @xmath8 is the volume of a single grain .
we conclude in particular that the aes notion of quantum dissipation at @xmath9 can only be applied to tunnel junctions connecting infinitely large conductors .
if the volume of the conductors is finite the aes picture can be applied at finite temperatures only .
apparently , the notion of quantum dissipation should be treated with care .
a qualitative discussion of the problem `` dissipation versus dynamical screening '' in the case of superconducting grains can be found in a recent lecture by altshuler@xcite .
the remainder of the paper is organized as follows . in sec .
[ choice ] we formulate the model and write physical quantities in terms of functional integrals over anticommuting fields .
we decouple the coulomb interaction by an integration over auxiliary fields . in sec .
[ aes ] we consider the aes action and show that it does not correspond to the physics of normal metals with a screened coulomb interaction . in sec .
[ metal ] we develop a diagrammatic technique .
we show what kind of diagrams correspond to the aes action and demonstrate that , at low temperatures , additional contributions arise . summing up these additional diagrams we arrive at expressions corresponding to the dynamically screened coulomb interaction of a normal metal .
[ sigma ] contains a derivation of a new action for a granular metal with electron - electron interaction which is applicable at arbitrary temperatures .
the results are compared with the results obtained diagrammatically and with those obtained on the basis of the aes action .
our results are summarized in the conclusion .
we consider an array of normal metallic grains coupled to each other .
the dimensionality of the array may be arbitrary .
the grains are assumed to contain imperfections : there can be impurities inside them as well as on their surface .
this implies that the electron motion in the grains is chaotic .
the mean level spacing in a single grain is equal to @xmath6 .
however , the shapes of the grains need not considerably differ from each other and we assume for simplicity that the grains are arranged in a cubic lattice .
the electrons can hop from grain to grain ; moreover , they interact with each other .
we are interested only in the long range part of the coulomb interaction and write it in a simplified form describing charging of the grains .
in such a formulation the spin of the electrons is not important and can be taken into account at the end of the calculations when writing proper densities .    under these assumptions
the electron motion can be conveniently described by a functional integral of the type @xmath10\right ) { \cal d}\psi , \ ] ] with an action @xmath11 $ ] that can be written in the form @xmath12=s_{{\rm g}}[\psi ] + s_{{\rm t}}[\psi ] + s_{{\rm c}}[\psi ] .
\label{a1}\ ] ] the action @xmath13 $ ] in eq .
( [ a1 ] ) describes the electron motion within the grains and we write it in the form @xmath14=\sum_{i}\int_{0}^{\beta } drd\tau \psi _ { i}^{\ast } ( \tau , r)\left ( \partial _ { \tau } -\mu + \hat{h}_{i}\right ) \psi _ { i}(\tau , r ) , \label{a2}\ ] ] where @xmath15 is a fermion field in grain @xmath16 at imaginary time @xmath17 and coordinate @xmath18 ( @xmath19 is its complex conjugate ) , @xmath20 is the chemical potential , and @xmath21 is the inverse temperature . the fermion fields @xmath22 satisfy the condition @xmath23    the operator @xmath24 in eq .
( [ a2 ] ) includes scattering by impurities and can be written as @xmath25 where @xmath26 is the impurity potential .
we assume that it is gaussian distributed with the correlation @xmath27    the action @xmath28 $ ] stands for tunneling between the grains , @xmath29=\sum_{ij}\int_{0}^{\beta } d\tau \psi _ { i}^{\ast } ( \tau ) t_{ij}\psi _ { j}(\tau ) ,   \label{a4}\ ] ] and @xmath30 $ ] describes the charging of the grains @xmath31=\frac{e^{2}}{2}\sum_{ij}\int_{0}^{\beta } d\tau n_{i}(\tau ) c_{ij}^{-1}n_{j}(\tau ) , \label{a5}\ ] ] where @xmath32 is the density field in grain @xmath16 , @xmath33 is the electron charge , and @xmath34 is the capacitance matrix .
( [ a1 ]  -  [ a5 ] ) describe the model completely and one can start explicit calculations .    due to the coulomb interaction ,
the action @xmath35 $ ] is not quadratic in @xmath36 and the integration over @xmath36 can not be performed immediately . a convenient way to proceed in such a case is to decouple the term @xmath37 $ ] by a gaussian integration over auxiliary fields .
this transformation can be written as follows : @xmath38}={\cal n}\int { \cal d } ve^{-s_{2}[v]}\int { \cal d}(\psi ^{\ast } , \psi ) e^{-s_{1}[\psi , v ] } , \label{a6}\ ] ] where @xmath39 is a normalization factor and @xmath40 is the decoupling field .
the action @xmath41 $ ] in eq .
( [ a6 ] ) has the form @xmath42=\frac{1}{2}\int_{0}^{\beta } d\tau \sum_{ij}v_{i}(\tau ) c_{ij}v_{j}(\tau ) ,   \label{a7}\ ] ] which shows that the variable @xmath43 has the meaning of a voltage on grain @xmath16 at the time @xmath17 .
the new effective action @xmath44 $ ] reads @xmath45=s_{{\rm g}}[\psi , v]+s_{{\rm t}}[\psi ] , \ ] ] with @xmath46=\sum_{i}\int_{0}^{\beta } d\tau \psi _ { i}^{\ast } ( \tau ) \left ( \partial _ { \tau } -\mu + \hat{h}_{i}+iev_{i}(\tau ) \right ) \psi _ { i}(\tau ) .
\label{a8}\ ] ] one more transformation can be performed exactly . following the procedure of ref .
where a single grain was considered , we represent @xmath43 in the form @xmath47 where @xmath48 is the static part and @xmath49 .
the replacement @xmath50 with @xmath51 does not violate the condition given by eq .
( [ a51 ] ) and we can remove the variable @xmath52 from the action @xmath53 $ ] , eq .
( [ a8 ] ) . however , this variable appears in the action @xmath54 $ ] as an additional phase of the field @xmath55 . as concerns the static part @xmath48 ,
its fluctuations can be neglected even in a single isolated grain , provided the temperature is high @xmath56  @xcite .
this restriction is not necessary for the system of the coupled grains in the limit of large conductance @xmath57 considered here .
using eqs .
( [ a1]- [ a10 ] ) we finally reduce the calculation of physical quantities to the computation of a functional integral of the form @xmath58-s_{1}[\psi , \phi ] -s_{2}[\phi ] \right ) { \cal d } \psi { \cal d}\phi ,   \label{a11}\ ] ] where @xmath59=\sum_{i}\int_{0}^{\beta } d\tau \psi _ { i}^{\ast } ( \tau ) \left ( \partial _ { \tau } -\mu + \hat{h}_{i}\right ) \psi _ { i}(\tau ) ,   \label{a12}\ ] ] @xmath60=\sum_{ij}\int_{0}^{\beta } d\tau t_{ij}\psi _ { i}^{\ast } ( \tau ) \psi _ { j}(\tau ) \exp \left ( i\phi _ { ij}\left ( \tau \right ) \right ) , \label{a13}\ ] ] @xmath61=\int_{0}^{\beta } d\tau \sum_{ij}\frac{c_{ij}}{2e^{2}}\frac { d\phi _ { i}(\tau ) } { d\tau } \frac{d\phi _ { j}(\tau ) } { d\tau } ,   \label{a14}\ ] ] and @xmath62 .    eqs . ( [ a11]-[a14 ] ) completely specify the model that will be studied in the subsequent sections .
disorder is still present in the operators @xmath63 , eq .
( [ a3 ] ) , but the mean free path within the grains is assumed to be large , @xmath64 , where @xmath65 is the fermi momentum . in this limit
, the macroscopic conductivity is determined mainly by the tunneling conductance @xmath1 , @xmath66 where @xmath67 is the tunneling resistance .
a functional integral formulation was used in ref .   to treat the quantum dynamics of a josephson junction .
the action derived in that work ( aes action ) was used later in the context of a normal tunnel junction  @xcite , as well as in connection with arrays of tunnel junctions in a number of papers  @xcite .
the aes action can be obtained from eqs .
( [ a11]- [ a14 ] ) .
first , one should average the green function @xmath68 corresponding to the action @xmath69 $ ] , eq .
( [ a12 ] ) , over impurities in the first born approximation , which gives ( in fourier representation ) @xmath70 where @xmath71 is defined in eq .
( [ a30 ] ) and @xmath72 . after
that one should make a cumulant expansion in the term @xmath73 $ ] , eq .
( [ a13 ] ) . keeping only the second order of the expansion and assuming for simplicity that the capacitance matrix @xmath74 is diagonal one arrives at the aes action @xmath75 $ ]  @xcite , @xmath76=\frac{1}{4e_{c}}\sum\limits_{i}\int\limits_{0}^{\beta } d\tau \left ( \frac{d\phi _ { i}}{d\tau } \right ) ^{2}+\frac{g_{t}}{2 } \sum\limits_{i , j}\int\limits_{0}^{\beta } d\tau d\tau ^{\prime } \alpha ( \tau -\tau ^{\prime } ) \left ( 1-\cos \left ( \phi _ { ij}(\tau ) -\phi _ { ij}(\tau ^{\prime } ) \right ) \right ) , \label{s}\ ] ] where @xmath77 in eq .
( [ s ] ) @xmath78 is the charging energy @xmath79 . the second term in eq .
( [ s ] ) contains the sum over neighboring grains and @xmath80 is the tunneling energy between them . in the limit of small phase fluctuations
we can expand the second term in eq .
( [ s ] ) with respect to @xmath81 and obtain @xmath82=\frac{1}{4e_{c}}\sum\limits_{i}\int\limits_{0}^{\beta } d\tau \left (   \frac{d\phi _ { i}\left ( \tau \right ) } { d\tau } \right ) ^{2}+ \frac{g_{t}}{2 } \sum\limits_{i , j}\int\limits_{0}^{\beta } d\tau   d\tau ^{\prime } \alpha ( \tau -\tau ^{\prime } ) \left ( \phi _ { ij}(\tau ) -\phi _ { ij}(\tau ^{\prime } ) \right ) ^{2}.   \label{s2}\ ] ] using the periodicity of all functions in @xmath83 we write the fourier expansion as @xmath84 where @xmath85 are matsubara frequencies .
hence we obtain for the action @xmath82=\frac{1}{4e_{c}}\sum\limits_{i}t\sum\limits_{\omega _
{ n}}\omega _ { n}^{2}\mid \phi _ { i}(\omega _ { n})\mid ^{2}+\frac{g_{t}}{\pi } \sum\limits_{i , j}t\sum\limits_{i\omega _ { n}}\mid \omega _ { n}\mid \mid \phi _ { ij}(i\omega _ { n})\mid ^{2}. \label{s3}\ ] ] the second term in eq .
( [ s3 ] ) describing the tunneling is linear in frequency and keeps its form down to @xmath9 .
therefore , it was attributed to quantum dissipation  @xcite . performing a fourier transformation in space @xmath86 we rewrite the action @xmath87 $ ] as @xmath82=\frac{1}{2}\sum\limits_{{\bf k}}t\sum\limits_{\omega _ { n}}\left ( \frac{\omega _ { n}^{2}}{2e_{c}}+\frac{4g_{t}}{\pi } |\omega _ { n}|\sum_{a}(1-\cos { \bf k}{\bf d}_{a})\right ) |\phi _ { \omega _ { n } } ( { \bf k})|^{2 } \label{smallfluct } , \ ] ] where @xmath88 , @xmath89 , @xmath90 , @xmath91 @xmath92  are vectors connecting the centers of the neighboring grains .
the phase - phase correlation function @xmath93 can be immediately calculated using eq .
( [ smallfluct ] ) and we obtain @xmath94    suppose that one assumes , following refs .
, that eq .
( [ s ] ) is applicable to a system of tunnel junctions down to @xmath9 .
this would imply in particular that eq .
( [ cor ] ) is correct for @xmath95 , too .
however , the derivative @xmath96 is proportional to the voltage on the grains which in turn is linearly related to the charge density . multiplying eq .
( [ cor ] ) by @xmath97 , one thus obtains the density - density correlation function @xmath98 .
let us consider limit of small frequencies and wave vectors @xmath99 . in this limit
, the network should correspond to a bulk normal metal , where one should have the propagator corresponding to a screened coulomb interaction .
however , the limit @xmath100 , @xmath101 in eq .
( [ cor ] ) does not reproduce the screened coulomb interaction .    how to resolve this discrepancy ? in the next section
we carry out an expansion in @xmath102 and demonstrate that keeping only the second order , which led us to eq .
( [ s ] ) is not sufficient and that , in the limit of small frequencies or temperatures , one should sum up an infinite class of diagrams . in section [ sigma ]
we will then show how these processes can be included into an effective action formulation structurally similar to the aes approach .
in this section we present a detailed derivation of the density - density correlation function @xmath103 summing an infinite series of diagrams . for the computation we use the standard diagram technique  @xcite with a modification suggested in refs .   to include tunneling between the grains . within this technique
one should make an expansion in the tunneling @xmath102 , the phases @xmath104 and in the impurity potential @xmath26 .
below we denote the electron green functions by solid lines , the phases @xmath105 by wavy lines and the tunneling elements @xmath102 by crossed circles .
impurity propagators corresponding to eq .
( [ a3 ] ) are denoted by dashed lines .    before turning to the quantitative discussion ,
let us briefly outline the skeleton of the analysis .
to lowest order in the tunneling matrix elements , the polarization operator of the system is represented by the process depicted in fig .
[ grainver ] .
this diagram describes the ( absolute square of ) amplitude for tunneling from one grain into the next and the subsequent relaxation inside the target grain ( represented through the green function ) .
the quantitative evaluation of this process leads to the coupling constant @xmath106 of the aes approach .    on this level ,
coherence effects associated to mulitiple tunneling and/or impurity scattering are completely neglected ; the tunneled electron forgets about its phase memory implying that the dissipative picture of the aes approach obtains .
indeed , the two green functions depicted in fig . [ grainver ] live in different grains which means that no phase correlation is possible .
how can this picture change in principle ?
including one more order in the tunneling , leads to corrections of the structure shown in fig .
[ hightun ] .
what makes these processes qualitatively different from those discussed above is that they include the possibility of _ phase coherent _ multiple impurity scattering : the two green functions labeled @xmath16 or @xmath107 scatter off the same impurities and/or chaotic potential where the associated scattering phases may cancel due to the fact that one of the green functions is advanced , the other retarded .
the net two - particle modes emanating from such processes , diffusons in the context of disordered systems , are represented through the hatched regions of fig .
[ diffuson ] , right .
the key question now is , are such processes relevant or not ?
a crude estimate can be given as follows : it is known that for low temperatures ( for a precise definition of low , see below ) , each diffuson mode diverges as @xmath108 .
the fact that we are considering processes of higher order in the tunneling amplitude introduces one more power of @xmath106 .
thus , for @xmath109 , higher order processes become as important as the first order contribution . for lower temperatures , these processes
_ must _ be taken into account to obtain a physically correct picture .
in the following we will put this discussion onto a firm basis and discuss how the coherent tunneling series can be summed in a controlled way .    before starting the computation for a granular system
we want to mention an important difference between the diagrams for the current - current and density - density correlation functions . calculating the classical conductivity for a granular system in the lowest order in the tunneling we need not renormalize the current vertices by impurities .
to understand this fact let us recall a well known result for bulk metals .
the diagram for the drude conductivity is shown in fig .
[ vertex ] .    the second diagram in fig .
[ vertex ] , renormalizing the current vertex , gives an additional integral over momentum that yields zero , @xmath110 in eq .
( [ b1 ] ) , the green function @xmath111 is determined by eq .
( [ a15 ] ) .
diagrams with many impurity lines yield zero as well and the current vertices are not renormalized .
the same is true for granular metals , which can be seen after a proper replacement of the current @xmath112 as concerns the polarization operator or the density - density correlation function , the situation is different . in this case , instead of vector vertices we have scalar ones and the renormalization due to impurities can be important .    expanding the functional integral , eq .
( [ a11 ] ) in the tunneling term @xmath113 , eq .
( [ a13 ] ) , we obtain to second order in @xmath102 and in @xmath105 the diagrams represented in fig .
[ grainver ] .
the two electron lines relate to different grains . as the impurities in the different grains
are not correlated , impurity lines connecting two green functions are do not exist ( although the diagram represents a polarization loop . ) the second and third diagrams in fig .
[ grainver ] are equal to each other and have an opposite sign with respect to the first diagram .
the analytical result for the sum of the three diagrams reads @xmath114 where the dimensionless tunneling conductance @xmath1 is given by eq .
( [ sg ] ) .
we see from eq .
( [ b3 ] ) that , although the green functions @xmath115 determined by eq .
( [ a15 ] ) contain @xmath71 , all information about the disorder drops out of the polarization loop . to calculate the phase - phase correlation function
one should notice that the bare propagator corresponding to @xmath116 , eq . ( [ a14 ] ) , has the form @xmath117 and diverges at @xmath118 .
therefore , one should sum up the infinite geometrical series represented in fig .
[ series ] .
= 1.2 cm    the summation can be easily performed and we arrive at the propagator @xmath119 , eq .
( [ cor ] ) . summarizing ,
we conclude that , in order to reproduce the results obtained in the quadratic approximation of the aes action , eqs .
( [ s2 ] , [ s3 ] , [ smallfluct ] ) , one should sum up the diagrams of fig .
[ series ] . on the other hand
, the form of the propagator @xmath120 does not correspond to the propagator of a normal metal with a screened coulomb interaction ( see _ e.g. _  ) , which signals that something is missing in this analysis .    to resolve this paradox
one should consider diagrams of higher order in tunneling , which complicates the calculations .
first of all , let us note that in order to consider higher order diagrams one should specify tunneling more explicitly .
we assume that the radius of the area at which the grains contact each other much exceeds atomic distances ( large number of conducting channels ) . at the same time
, the potential barrier between the grains can be arbitrarily large , which can lead to an arbitrary tunneling conductance .
such a situation can be conveniently modeled by random tunneling matrix elements @xmath121 correlated as @xmath122 where @xmath121 are written for the same contact .
correlations of the matrix elements of different contacts are assumed to be zero . to simplify notations
we do not write subscripts numerating the contacts . taking into account eq .
( [ average ] ) is indeed very important . in order to illustrate this fact we consider two diagrams represented in fig .
[ hightun]a .
we see immediately the difference between them : the first diagram is not averaged over @xmath121 and can contains four different momenta . after averaging over
@xmath121 one obtains the second diagram containing only three different momenta .
as it has been mentioned , we will expand not only in the tunneling amplitude @xmath123 but also in the phases @xmath124 . performing the expansion in the phases and denoting them by wavy lines we have three different possibilities as depicted in fig .
[ hightun]b .
now we can start averaging over impurities . as usual  @xcite , we average first the green functions and obtain eq .
( [ a15 ] ) for the average green function @xmath115 .
then , we can consider averages with impurity lines connecting the green functions on opposite sides of the diagrams in fig .
[ hightun]b .
diagrams with non - intersecting impurity lines are most important and we recognize the well - known diffusons .    in the present paper
we consider the limit of sufficiently small grains , such that all relevant energies like temperature @xmath125 , frequency @xmath100 , tunneling energy @xmath80 , etc .
are much smaller than the thouless energy @xmath126,where @xmath127 is the diffusion coefficient of a single grain and @xmath128 is the radius of the grain . in this limit
the spectrum of the diffuson in a single grain is zero - dimensional ( no dependence on momenta inside the grain ) .
coupling between the grains leads to a dependence on quasi - momentum @xmath129 .
diagrams that should be summed in order to give the complete form of the diffuson are represented in fig .
[ diffuson]a .
now let us present some details of the calculation of the diagrams in fig .
[ hightun ] .
the different classes of diagrams contributing to the polarization operator of the granular metals before impurity averaging are shown in fig .
[ hightun]b .
before impurity averaging the sum over matsubara frequencies for the first diagram in fig .
[ hightun]b is @xmath130 writing this expression we used eq .
( [ average ] ) .
after impurity averaging of the first diagram in fig .
[ hightun]b , we obtain the diagram represented in fig .
[ diffuson]a .
now the disorder averaging can be done and we obtain , in a standard way , the average of the product of two green functions for a granular metal , @xmath131 . here
@xmath133 being the diffuson for a single isolated grain .
the equation for the complete diffuson  @xcite is shown in fig .
[ diffuson]b .
the self - energy @xmath134 in eq .
( [ dyson ] ) is the sum of the three diagrams in fig .
[ diffuson]c .
the second and third diagrams in fig .
[ diffuson]c are equal to each other and have the opposite sign with respect to the first diagram . solving eq .
( [ dyson ] ) we obtain @xmath135    the calculation of the second diagram in fig .
[ hightun]b is analogous to the calculation of the first one .
the third diagram in fig .
[ hightun]b is independent ofthe frequency @xmath136 and yields a constant .
the diagram in fig .
[ diffuson]a contains only one diffuson .
more complicated diagrams can be drawn and one of them is depicted in fig .
[ dif ] .
= 2.5 cm    however , all such many - diffuson contributions can be neglected provided the tunneling conductance @xmath1 , eq .
( [ sg ] ) , is large .
a direct calculation of this diagram shows that it involves the small parameter @xmath137 .
this is similar to expansions in the diffusion modes for a bulk metal where the inverse conductivity is the expansion parameter .
the condition @xmath57 means that we are far from the anderson metal - insulator transition .
the result of averaging and summation of the diagrams represented in fig .
[ hightun ] can be written as @xmath138 we note here that the functions @xmath139 and @xmath140 have a different sign . comparing the function @xmath141 , eq .
( [ newdif ] ) , with the quantum dissipation part @xmath142 , eq .
( [ b3 ] ) , represented in fig .
[ grainver ] we see immediately that the contribution @xmath141 , eq . ( [ newdif ] ) , can be neglected only if the temperature @xmath125 is sufficiently high @xmath143 ( we recall that the static component of the phase @xmath124 has been neglected and therefore @xmath144 in eq .
( [ newdif ] ) .
so , the notion of quantum dissipation is applicable essentially at non - zero temperatures only for a granular material .
however , it can be used for a contact connecting two bulk metals with an infinite volume .    at low temperatures @xmath145
, there is no reason to neglect the contribution @xmath146 , eq .
( [ newdif ] ) , and its presence changes completely the form of the propagator . adding the functions @xmath147 , eq .
( [ b3 ] ) , and @xmath148 , eq . ( [ newdif ] ) , we obtain the total polarization @xmath149 @xmath150 proper diagrams contributing to this function are represented in fig .
[ complete ] .
= 2.5 cm    now , we can write the final expression for the phase correlation function @xmath151 that should be written instead of @xmath120 , eq . ( [ cor ] ) . using the bare propagator @xmath152 , eq .
( [ b4 ] ) , and the self - energy part @xmath153 , eq .  ( [ total ] ) , we write the phase propagator @xmath154 in the form @xmath155    using the relation between the phase @xmath124 and the voltage @xmath8 , eq .
( [ a10 ] ) , we derive an equation for the effective coulomb interaction @xmath156 on the other hand , the dynamically screened coulomb interaction in a disordered bulk metal can be written in the form @xmath157 where @xmath158 is the bare coulomb potential .
we see that eqs .
( [ f1 ] , [ f2 ] ) written in the continuum limit correspond to eq .
( [ f3 ] ) .
thus , we conclude that , when calculating physical quantities , the main contribution is due to small @xmath99 only in the limit @xmath159 . provided
this inequality is satisfied , the properties of the tunnel junction array are equivalent to the properties of a weakly disordered bulk metal with coulomb interaction . in the opposite limit @xmath160 , one can describe the system in the `` quantum dissipation '' approximation .
calculations were performed in this section using a diagrammatic expansion in the lowest order in the field @xmath124 . on the other hand , the aes free energy functional , eq .
( [ s ] ) , is written for arbitrary values of @xmath124 .
can one generalize eq .
( [ f2 ] ) to arbitrary values of @xmath124 ? we will try to do this in the next section , deriving a simplified version of a replica @xmath3-model and using a non - trivial saddle point .
in the preceding sections we analyzed effects of the coulomb interaction on the electron motion in granular materials using a diagrammatic expansion .
this approach works well when fluctuations of the phase @xmath124 or , in other words , fluctuations of the voltages on the grains are small .
however , one can imagine situations when the fluctuations are large .
for example , one can try to consider non - perturbative excitations like instantons  @xcite corresponding to an integer charge transfer between grains . in principle , the aes action , eq .
( [ s ] ) , can describe such excitations very well but , as we have seen previously , it can only be used at sufficiently high temperatures @xmath160 .
so , our task now is to derive an action that will be applicable at lower temperatures .
the derivation at temperatures @xmath145 is more complicated because now we should explicitly take into account disorder , which was not necessary for the derivation of the aes action ( strictly speaking , disorder determines the mean free time @xmath161 in the green functions used for derivation of the aes functional but this is a trivial contribution ) . on the other hand ,
the problem is not as complicated as the problem of anderson localization in the presence of interactions considered in refs .  .
in our diagrammatic expansions we considered the limit of large conductances @xmath57 and neglected all weak localization corrections .    nevertheless , we start the derivation of the action using the replica @xmath3-model approach for interacting systems suggested by finkelstein @xcite .
necessary simplifications will be made in the @xmath3-model .
although there are difficulties in using the replica @xmath3-models for non - perturbative calculations , our goal is more modest and the final results obtained in the limit @xmath57 neglecting all weak localization corrections will not depend on replica indices at all .
the model under consideration has been formulated in section  [ choice ] .
the derivation of the @xmath3-model can be carried out starting from eqs .
( [ a1]-[a5 ] ) . instead of the fields @xmath36 and @xmath8 one should write fields @xmath162 and @xmath163 carrying the replica index @xmath164 upon calculating physical quantities for an arbitrary @xmath18 , one should put @xmath165 .
a proper @xmath3-model for bulk systems with interaction has been derived in refs .  .
a generalization to granular metals can been made without difficulties following ref .  , where a @xmath3-model was written for a granular metal without interaction .
although one can proceed along the lines of refs .  , _
i.e. _ expanding the action obtained after integration over @xmath36 in the field @xmath8 , a more economic way is to use the substitution given by eqs .
( [ a90 ] , [ a10 ] ) , as well as the subsequent eqs .
( [ a11]-[a14 ] ) . of course
, this is already an approximation because we neglect the static component of the phase @xmath124 .
however , for our purposes this is not an essential restriction because we do not want to consider effects of localization .
moreover , working in the limit @xmath166 or @xmath57 allows us to ignore the static component of @xmath124 .
the derivation of the @xmath3-model starting from eqs .
( [ a11]-[a14 ] ) is practically the same as the one presented in ref .  .
as usual , one decouples the `` effective interaction '' of the type @xmath167 that appears after averaging over impurities by a gaussian integration over matrices @xmath168 .
these matrices contain as elements both the replica indices and matsubara frequencies ( or two imaginary times @xmath17 and @xmath169 ) . neglecting the tunneling term , eq .
( [ a13 ] ) , one would obtain a zero - dimensional @xmath3-model ( we assume that all relevant energies are smaller than the thouless energy @xmath170 of a single grain ) .
the relevant tunneling term in the @xmath3-model is obtained by a cumulant expansion in the tunneling . the only difference with respect to the model without interaction of ref .
is the presence of the phases @xmath124 in eq .
( [ a13 ] ) , which leads to additional phase factors in the term describing the coupling between the grains .
the final result for the free energy @xmath171 $ ] can be written as @xmath172=f_{2}[\phi ] + f_{{\rm \omega } } + f_{{\rm t}}[q,\phi ] , \label{full5 } \\ & & \qquad f_{2}[\phi ] = \frac{1}{2e^2}\sum_{i , j}\int { \rm tr}c_{ij}\frac{\partial \phi _ { i}\left ( \tau \right ) } { \partial \tau } \frac{\partial \phi _ { j}\left ( \tau \right ) } { \partial \tau } d\tau , \label{full6 } \\ & & \qquad f_{\omega } [ q]=\frac{\pi } { \delta } \sum_{i}\int { \rm tr } ( \hat{\omega}_{\tau } q_{i}\left ( \tau , \tau \right ) ) d\tau , \label{full7 } \\ & & \qquad f_{{\rm t}}[q,\phi ] = -\frac{g_{t}}{4}\sum_{i , j } \int { \rm tr}\left ( e^{i\phi _ { ij}(\tau ) } q_{i}\left ( \tau , \tau ^{\prime } \right ) e^{-i\phi _ { ij}(\tau ^{\prime } ) } q_{j}\left ( \tau ^{\prime } , \tau \right ) \right ) d\tau d\tau ^{\prime } , \label{full8}\end{aligned}\ ] ] where @xmath62 and the symbol tr implies a trace over both replica indices .
the field @xmath124 is diagonal in the replica indices , while the matrix @xmath173 is a @xmath174 matrix with the constraints @xmath175    the action of the operator @xmath176 in eq .
( [ full7 ] ) on an arbitrary function @xmath177 is given by the relation @xmath178 in frequency representation this operator is equal to a vector of fermionic matsubara frequencies .
the matrix @xmath173 , satisfying the constraints , eq .
( [ e1 ] ) , can be conveniently parametrized as @xmath179 where @xmath180 and @xmath181 are unitary matrices @xmath182 the matrix @xmath183 in frequency representation has the form @xmath184 the limiting cases of the free energy @xmath171 $ ] , eqs .
( [ full5]-[e2 ] ) are rather simple . without electron - electron interaction
the phase difference between neighboring grains equals to zero @xmath185 .
such an action for granular metals corresponds to that derived within the supersymmetry scheme by one of the authors  @xcite .
the second limiting case is achieved at sufficiently high temperatures . in this limit disorder
is not important and the matrix @xmath168 does not fluctuate , being equal to the value @xmath183 , eq .
( [ e4 ] ) . inserting @xmath186 into eqs .
( [ full5]-[full8 ] ) we reproduce immediately the aes action , eq .
( [ s ] ) ( in fact we obtain the @xmath18 times replicated aes functional but only one replica field is important for us ) .    at
not too high a temperature , eqs .
( [ e3]-[e5 ] ) lead to a non - trivial behavior even if the disorder is weak and all effects related to weak localization can be neglected .
the limit of weak disorder enables us to neglect fluctuations of @xmath168 at given @xmath124 and take its value from a saddle - point equation .
we thus reduce the computation of the functional integral with the free energy functional to the solution of a saddle - point equation , substituting the solution to the free energy @xmath171 $ ] , eqs .
( [ full5]-[full8 ] ) , calculating a functional integral with the reduced free energy @xmath187 $ ] @xmath188=f[\bar{q},\phi ]   \label{e6}\ ] ] where @xmath189 is the solution of the saddle - point equation for a given @xmath124 .    minimizing the functional @xmath171 $ ] , eqs .
( [ full5]-[full8 ] ) , with the constraints , eqs .
( [ e1 ] ) , we obtain the following equation @xmath190d\tau ^{\prime \prime } + { \frac{\pi } { \delta } } [ \bar{q}_{i}\left ( \tau , \tau \right ) , \hat{\omega})]=0 , \label{e7}\ ] ] where the summation over @xmath107 is performed over the nearest neighbors of @xmath16 . in the commutator
@xmath191 $ ] one should exchange properly the times when changing the order of the functions .
( [ e7 ] ) represents an integral equation that enables us to find , in principle , the solution @xmath192 for arbitrary values of the parameter @xmath193 .
in general , this equation can be solved only numerically . however , for @xmath194 drastic simplifications arise . for large intergranular coupling
, the phase @xmath195 fluctuates only weakly from grain to grain which means that @xmath196 is small . under these conditions ,
a solution of the diagonal form @xmath197 can be sought for . as the solution
is assumed to have a diagonal form , all replica indices decouple and the equations can be solved for each replica separately .
therefore we can drop the replica indices and consider eq .
( [ e7 ] ) as an equation for a single function @xmath192 for a given function @xmath198 .    in the high temperature limit @xmath160 the second term in eq .
( [ e7 ] ) is much larger than the first one and we arrive at the solution @xmath199 which leads to the aes free energy functional , eq .
( [ s ] ) .    in the opposite limit of low temperatures @xmath159 eq .  ( [ e7 ] ) can be simplified using the assumption that @xmath192 varies slowly in space .
then , one may expand the exponentials in @xmath200 , which leads to the equation @xmath201 + \frac{i\pi } { \delta } [ \frac{\partial } { \partial \tau } , k_{i}]=0 .
\label{e10}\ ] ]    eq .
( [ e10 ] ) is presented without writing the integration over imaginary times @xmath17 explicitly , but of course this integration is implied . transforming this equation to frequency and momentum space
we obtain @xmath202 which gives the solution @xmath203 where @xmath204 and @xmath205 are the fourier transforms of the functions @xmath192 and @xmath43 respectively .
substituting this ansatz for @xmath206 back into the action we obtain @xmath207 = \frac{1}{4e_c } \sum_n \sum_k v_{k , n}\left(1 + 4e_c \frac{(2/\pi)g_t\sum_a(1-\cos k_ad)}{|\omega_n|+(2/\pi)g_t\delta \sum_a(1-\cos k_ad)}\right)v_{-k ,- n}\ ] ] where @xmath208 is the charging energy . in the limit of small characteristic momentum , @xmath209
, we recover the diffusively screened coulomb interaction in a disordered evironment .    from here
one can proceed to the calculation of physical observables such as the conductance , the density of states or others . to this end
one should introduce configurations @xmath210 fluctuating around the mean field configuration discussed above .
next one would follow the standard algorithm of doing calculations within the @xmath211-model approach : ( i ) express the quantity of interest in terms of @xmath168 , ( ii ) compute the functional average over the effective action @xmath212=f[v , t]$ ] as good as is possible .
let us outline how the connection between the enlarged aes - type formulation , ( [ full5 ] ) - ( [ full8 ] ) and the finkelstein approach for interacting disordered media can be made explicit . to do so , we subject our @xmath168 s to a gauge transformation , @xmath213 as a result ( i ) the hopping part of the action
, @xmath214 becomes @xmath195-independent .
however ( ii ) in the frequency part , @xmath215 .
gaussian integration over @xmath216 then produces an effective action @xmath212 $ ] which , after taking a continuum limit , is identified as the action of the finkelstein approach .
in this paper we constructed a theoretical framework to describe electronic transport in arrays of tunnel junctions or granular metals , in the limit of large inter - granule conductance , @xmath106 .
both disorder and the electron - electron interaction were taken into account .
the prime objective of this enterprise was to unify two large , and seemingly non - overlapping theories of interacting metallic compounds : the aes approach , focusing on charging phenomena in individual grains and their impact on large scale transport behaviour , and the finkelstein theory of disordered interacting metals with its emphasis on the interplay interaction / disorder .
the key to the reconciliation of these two approaches was to observe that for low enough temperatures , @xmath217 , the effective action underlying the aes theory becomes incomplete .
the physical reason is that for low temperatures the electrons and holes participating in the tunneling processes between individual grains maintain their quantum phase memory for a long time , largely in excess of the average tunneling time .
this means that the particle tunneling is not only accompanied by charging and dissipation ( as in the aes approach ) , but also by other , more long ranged physical processes .
specifically , we identified the quantum dissipation contained in the aes action as the high temperature limit of the long ranged screened electron - electron interaction , an observation first made in@xcite .
we re - emphasize that these phenomena were not bound to the presence of a significant disorder concentration ; indeed , none of the results discussed above , displayed dependence on some disorder concentration. the only thing that mattered was chaoticity of the electron motion on scales set by the phase coherence length , a condition that is practically always met in real life systems .    technically , two different routes for including these processes into the theory were proposed .
one consisted of the perturbative summation of diagram classes associated to multiple chaotic scattering . the other , based on an effective action formulation , indentified what is missing in the aes action . remarkably , this second approach readily led to a unification of the aes approach and the finkelstein nonlinear @xmath211-model for weakly disorderd interacting metals . with the benefit of hindsight ,
this fact is easy to understand : it is a well known fact in mesoscopic physics , that the long processes resulting from multiple scattering can be interpreted as a certain type of goldstone modes . within an effective action approach
, these modes must be described by an own degree of freedom , the @xmath168-matrices of the nonlinear @xmath211-model .
thus it is no surprise that we obtained an action of the type ( [ full5 ] - [ full8 ] ) , involving @xmath168-fields and the coulomb phase fields of the aes approach , as an effective description of the low energy phase of the system . for large temperatures ,
the fluctuations of the @xmath168-matrices became inessential and the aes action was retrieved . in contrast , the continuum limit of this action was identified as the finkelstein model and the @xmath168-matrices smoothly integrated into the aes approach , without leading to conceptual complications .
let us make some remarks on potential experimental ramifications of our findings .
it has been predicted in refs .   that two - dimensional normal arrays of tunnel junctions should undergo a kosterlitz - thouless - berezinskii ( ktb ) phase transition at a temperature @xmath218 of the order of the charging energy @xmath208 . at low temperatures , @xmath219
, the array should be in an insulating state ( coulomb blockade ) , whereas for @xmath220 the array is conducting .
in particular , at @xmath221 , the conductivity @xmath211 of the array is predicted to increase with temperature according to a square - root cusp dependence , @xmath222 , where @xmath223 is a constant of order unity .
the experiments  @xcite mentioned in the introduction searched for this ktb transition . in view of our results ,
the applicability of the ktb scenario should be governed by the parameter @xmath224 .
the experiments of refs . were done on arrays with relatively large grains , such that @xmath225 .
indeed , the results obtained in refs .   were in agreement with the ktb scenario . however , ref .
found thermally activated behavior of @xmath226 , rather than the predicted square root cusp dependence on temperature . on the other hand , the experiment by yamada _
et al . _
@xcite was done on arrays consisting of relatively small cu grains with a size of about @xmath227 , such that @xmath228 .
clearly this is beyond the range of applicability of refs . and the ktb scenario should be treated with care in this case .
remarkably , ref .   concludes good agreement of the results with the theory  @xcite .
given these discrepancies , we conclude that a careful analysis of the temperature - dependent conductivity of junction arrays in the framework of the theory presented in this paper would be of interest .
the authors thank a. andreev , r. fazio and a. tschersich for helpful discussions . a support of the graduiertenkolleg 384 and the sonderforschungsbereich 237
is greatly appreciated .