# Charge and Current Sum Rules in Quantum Media Coupled to Radiation

###### Abstract

This paper concerns the equilibrium bulk charge and current density correlation functions in quantum media, conductors and dielectrics, fully coupled to the radiation (the retarded regime). A sequence of static and time-dependent sum rules, which fix the values of certain moments of the charge and current density correlation functions, is obtained by using Rytov’s fluctuational electrodynamics. A technique is developed to extract the classical and purely quantum-mechanical parts of these sum rules. The sum rules are critically tested in the classical limit and on the jellium model. A comparison is made with microscopic approaches to systems of particles interacting through Coulomb forces only (the non-retarded regime). In contrast with microscopic results, the current-current density correlation function is found to be integrable in space, in both classical and quantum regimes.

###### Keywords:

Sum rules fluctuations radiation classical limit jellium^{†}

^{†}journal: J Stat Phys

∎

## 1 Introduction

This paper is about the equilibrium statistical mechanics of infinite (bulk) media, conductors and dielectrics, fully coupled to the radiated electromagnetic (EM) field. The crucial problem of statistical analysis of such models is to determine how fluctuations of various statistical quantities, like charge and current densities, scalar and vector potentials, electric and magnetic fields, etc., around their mean values are correlated in space and time. A special attention is devoted to the behavior of correlation functions at asymptotically large distances and to the so-called sum rules, which fix the values of certain moments of the correlation functions.

The studied models are composed of spinless charged particles which are quantum, but non-relativistic (i.e., they obey the Schrödinger equation and not the Dirac particle-antiparticle formalism); according to the correspondence principle, a quantum system admits the classical description in the high-temperature limit. On the other hand, the interaction of charged particles via the radiated EM field can be considered either non-relativistic (non-retarded) or relativistic (retarded). In the non-retarded regime, the speed of light is taken to be infinitely large, , ignoring in this way magnetic forces acting on the particles; in other words, charges interact pairwisely only via instantaneous Coulomb potentials. In the retarded regime, is assumed finite and the particles are fully coupled to both electric and magnetic (transverse) parts of the radiated field. The ultimate aim is to describe quantum particles in the retarded regime, but due to serious technical difficulties one usually starts with the simpler non-retarded regime.

Two complementary approaches exist in the theoretical study of charged systems. The microscopic one is based on the explicit solution of models defined by their microscopic Hamiltonians. The macroscopic one is based on the assumption of validity of macroscopic electrodynamics. Microscopic description is complicated and laborious, but if some explicit results are available their value is considerable and, e.g., they can reveal a restricted applicability of macroscopic methods.

Microscopic approaches are usually restricted to the non-retarded regime. A particular interest was devoted in the past to classical and quantum models of conductors. Classical particles interacting by the instantaneous Coulomb potential have been treated in the pioneering work by Debye and Hückel Debye23 (). An exponential shielding of the Coulomb potential, typical for the classical regime, was observed. One of the first classical analyses of a system of charged particles fully coupled to the radiation, based on the linearized Vlasov equation and Maxwell’s equations, was performed by Felderhof Felderhof64 (); Felderhof65 (), see also the monograph VanKampen67 (). The existence of general sum rules for certain moments of the classical charge density correlation functions was established by Stillinger and Lovett Stillinger68 (). The quantum sum rules, static and time-dependent, were the subject of numerous studies Pines66 (); Martin85 (); Martin86 (); Martin88 (); John93 (). A treatment of quantum Coulomb fluids, based on the use of the Feynman-Kac path integral representation of the thermal Gibbs weight Alastuey89 (), indicates that the charge-charge density correlation functions exhibit a long-range decay (see review Brydges99 ()). A microscopic analysis of conductor systems when transverse EM interactions are added to the Coulomb ones, was done recently Boustani06 (); Buenzli07 (); a conceptual problem related to the use of classical radiation was pointed out in Jancovici06 (). Microscopic models of dielectrics are less developed; we mention a fluid model of quantized polarizable particles Hoye81 (); Thompson82 () with the static dipole-dipole interaction.

Macroscopic phenomenological approaches usually allow us to predict, with less effort, basic features of relatively complicated complex physical systems in the retarded regime. A macroscopic theory of equilibrium thermal fluctuations of the EM field in conductors and dielectrics was proposed by Rytov Rytov53 (); Rytov58 (); Levin67 (), see also the Landau and Lifshitz Course of Theoretical Physics LP () and the book Sitenko95 (). For bulk systems, an interesting phenomenon occurs in two-point correlation functions of microscopic quantities like the induced electric potential or field: the presence of transverse interactions beyond Coulomb cancel out some (spatially nonintegrable) long-ranged terms observed in the non-retarded regime Jancovici06 (); LP (). Main applications of Rytov’s theory were related to inhomogeneous situations. In the famous Lifshitz paper Lifshitz56 (), Rytov’s theory was used in the context of the Casimir attraction between parallel dielectric plates. Another application appeared in the field of thermally excited surface EM waves Joulain05 (). Recently, the theory was used in the study of retardation effects on the long-range surface charge density correlation function between two distinct media Samaj08 (); Jancovici09 ().

Macroscopic approaches, being essentially of mean-field type, are expected to provide reliable results for the leading terms in the asymptotic long-wavelength behavior of correlations. Those results can be recast as sum rules for moments of correlations integrated in real space, or as asymptotic large-distance behaviors. As concerns these asymptotic behaviors, macroscopic mean-field teories fail when fluctuations become crucial. The failure of mean-field predictions is well illustrated by the existence of tails in static charge-charge correlations of quantum Coulomb systems, as explained in Alastuey99 (); see also the exact calculation of the tail for hydrogen plasma at low densities Cornu97 (). Another implicit assumption for the validity of macroscopic approaches is that correlations decay sufficiently fast. It is then legitimate to perform local averages for EM fields slowly varying in space, which reflects itself in the introduction of the frequency dependent dielectric function in Maxwell’s equations.

The present paper concerns sum rules for the bulk charge and current density correlation functions in quantum media, conductors and dielectrics, fully coupled to the radiation. Although this subject (in the non-retarded version) was in the center of interest of microscopic theories in the past, it remained practically unheeded by macroscopic approaches. Here, a sequence of static or time-dependent sum rules, known or new, is obtained by Rytov’s fluctuational electrodynamics. A technique is developed to extract the classical and purely quantum-mechanical parts of these sum rules. The sum rules are critically tested on a jellium model of conductors. A comparison is made with microscopic approaches to systems of particles interacting through Coulomb forces only; in contrast to these microscopic approaches, the current-current density correlation function is found to be integrable in space, in both classical and quantum regimes.

The paper is organized as follows. In Sect. 2, we review briefly Rytov’s theory of EM field fluctuations and present basic expressions for the charge and current densities. The sum rules for the charge-charge, charge-current and current-current density correlation functions are derived and discussed in Sects. 3, 4 and 5, respectively. Sect. 6 deals with problematic points of both the Rytov fluctuational theory and microscopic approaches. Sect. 7 is devoted to an interesting comparison of the present results to those of Felderhof Felderhof64 (); Felderhof65 (), in the classical limit. Sect. 8 is the Conclusion.

## 2 Macroscopic fluctuational formalism

We consider the -dimensional space, defined by Euclidean vectors and time . The two-point functions considered in this paper will be translationally invariant in both Euclidean space and time, so we can use the spectral (Fourier) representation

(2.1) |

where denotes the frequency and the wave vector. The physical system of interest is a medium (conductor or dielectric) and an EM field present in it, which are in thermal equilibrium.

The medium is composed of moving charged particles which are assumed to be non-relativistic. In the long-wavelength scale much larger than the mean interparticle distance, the macroscopic characteristic of the medium is the frequency dependent dielectric function . We shall assume that the matter has no magnetic structure, i.e. it is not magnetoactive, and permeability . The medium is coupled to the EM field generated by the particles. The classical EM field is described by the scalar potential and the vector potential with component . In the considered Weyl gauge , the microscopic electric and magnetic fields are given by

(2.2) | |||||

(2.3) |

with being the velocity of light and being the unit antisymmetric pseudotensor. The elementary excitations of the quantized EM field are described by the photon operators of the vector potential components which are self-conjugate Bose operators. The retarded photon Green function tensor is defined by

(2.4) |

Here, denotes the vector-potential operator in the Heisenberg picture and the angular brackets represent the equilibrium averaging at temperature . For non-magnetoactive media, the Green function tensor possesses the symmetry

(2.5) |

The EM fields are random variables which fluctuate around their mean values governed by Maxwell’s equations. To describe fluctuations of the EM fields, Rytov Rytov53 (); Rytov58 (); Levin67 (); LP () studied the effect of a weak classical current due to the random particle motion in the medium. This current acts as an “external force” on the vector-potential operator in the interaction Hamiltonian

(2.6) |

The mean values of the components of the vector-potential operator can be thus expressed in terms of the Green function (2.4) by using Kubo’s linear response in current,

(2.7) |

The mean value satisfies macroscopic Maxwell’s equations due to the classical current . This fact automatically implies a set of differential equations of dyadic type fulfilled by the Green function tensor. In particular, in the frequency Fourier space, we have

(2.8) |

This set of equations has to be supplemented by the obvious boundary condition of regularity for the distance going to infinity. In the case of a spatially homogeneous infinite medium, , under the Fourier transform with respect to the difference Eq. (2.8) reduces to a set of algebraic equations, whose solution reads

(2.9) |

Here, . The first/second term on the right-hand-side (r.h.s.) of this expression represents the longitudinal/transverse component of the Green function. Formula (2.9) does reduce to its well-known form in the vacuum when one sets .

Eq. (2.7) defines as the tensor of susceptibilities corresponding to the random variable . The fluctuation-dissipation theorem tells us that the fluctuations of random variables are expressible in terms of the corresponding generalized susceptibilities. For the present symmetry (2.5), the theorem implies

(2.10) |

where is the inverse temperature and is the Fourier transform of the symmetrized correlation function

(2.11) |

We note that the symmetrization (2.11) is needless in the classical case.

The fluctuation formula (2.10) enables us to calculate the correlations of arbitrary statistical quantities provided that they are expressible in terms of the components of the vector-field potential, within the classical EM field and in the gauge defined by Eqs. (2.2) and (2.3). Let a scalar quantity be expressible as , where are operators acting on and variables. Within the spectral representation with a single frequency and wave vector, , this relation takes an algebraic form . It follows from the definition (2.1) that the Fourier transform of the symmetrized correlation function of statistical quantities and , , is then determined by

(2.12) | |||||

where given by (2.10).

In this paper, we shall concentrate on the volume charge density and the electric current density . The charge density is expressible as

(2.13) |

or, in the spectral representation, as

(2.14) |

The electric current density is expressible in terms of the electric and magnetic fields as follows

(2.15) |

Substituting here the gauge relations (2.2) and (2.3), and using the summation formula for the elements of the unit antisymmetric pseudotensor

(2.16) |

we obtain

(2.17) |

or, in the spectral representation,

(2.18) |

It is a simple task to verify that both versions of the continuity equation

(2.19) |

are satisfied.

Henceforth, in order to simplify the notation we set and , i.e. the time and position differences between the two points in the correlation functions will be and , respectively.

## 3 Charge-charge density correlations

We start with the charge-charge density correlation function. Using the spectral representation (2.14) it is easy to show that

(3.1) |

where

(3.2) |

Note that ; the equality takes place in the limit . The fact that there is only one term of the order in the Fourier representation (3.1) is related to the applicability of the fluctuational theory in the long-wavelength limit (characteristic length is usually the mean interparticle distance); an exact formula for , valid for any value of , would contain also higher powers of .

The inverse Fourier transform of (3.1) leads to the zeroth-moment (neutrality) condition

(3.3) |

and the second-moment condition

(3.4) |

These sum rules do not depend on and so they are the same in both nonretarded and retarded regimes. Note that the sum rule (3.4) is not universal in the sense that its r.h.s. is a complicated function of the temperature and the medium characteristics.

### 3.1 Static case

In the static case, using the notation , the general result (3.4) can be rewritten in the form

(3.5) |

note that the subtraction of 1 from has no effect due to the operation .

The integral on the r.h.s. of Eq. (3.5) can be formally evaluated by using integration techniques and the general analytic properties of dielectric functions in the complex frequency upper half-plane, in close analogy with the treatment of the surface charge density correlation function in Sect. IV of Ref. Jancovici09 (). We shall not repeat the derivation procedure and only write down the final result. Although our -function defined in (3.5) differs from that considered in Jancovici09 (), it possesses the necessary symmetry for real and goes to zero in the asymptotic limit due to the asymptotic relation Jackson (); LL ()

(3.6) |

Here, is the plasma frequency of the medium composed of species (electrons and ions) with the number density , charge and mass ; since the masses of ions is much bigger than the mass of electron, it is usually sufficient to approximate the plasma frequency of the medium by with being the number of all electrons in all atoms in unit volume, and are the electron charge and mass, respectively. Let be a closed contour in the upper half-plane defined as , where is the path following the real axis, except it goes around the origin in a small semicircle whose radius tends to zero and is a semicircle at infinity (see Fig. 1). If we denote by all poles of the function in the region bounded by the contour , the procedure described in Ref. Jancovici09 () implies

(3.7) |

where means the residue. The dielectric function has no zeros in the upper half-plane, so the only source of the poles of inside is the function given by (3.2). Since can be expanded in as follows Gradshteyn ()

(3.8) |

it exhibits in the upper half-plane an infinite sequence of simple poles at the (imaginary) Matsubara frequencies

(3.9) |

We thus conclude that the second moment of the static charge-charge density correlation function (3.5) is expressible as

(3.10) |

In the upper half-plane, the dielectric function takes real values only just on the imaginary axis and this fact ensures that the second moment (3.10) is real. Moreover, decreases on the imaginary axis monotonically from (for dielectrics) or from (for conductors) at to 1 at ; the convergence of the sum on the r.h.s. of (3.10) is ensured by the asymptotic behavior (3.6).

The representation (3.10) is very useful in the high-temperature (classical) limit when each of the frequencies is much larger than , the corresponding terms in the sum on the r.h.s. of Eq. (3.10) vanish and therefore

(3.11) |

In the case of conductors, setting implies the universal Stillinger-Lovett sum rule Stillinger68 (). In the case of a dielectric of static dielectric constant , the above nonuniversal result has already been obtained by Chandler Chandler77 (); see also the article Alastuey00 (). It should be noted that, at finite temperature, a quantum system always contains a finite fraction of ionized charges and thus, in a strict sense, it is conductor. In the case of the so-called dielectrics, the corresponding screening length is huge, much larger than the size of considered macroscopic samples, and therefore one does observe dielectric properties at experimental length scales. This point is well illustrated by exact results for atomic hydrogen Ballenegger03 ().

The formula (3.11) can be rederived for dielectrics in an alternative way by extending screening and linear response ideas in classical Coulomb fluids Jancovici95 (); Pines66 () to dielectrics as follows Janco (). If we put an infinitesimal charge at position , the dielectric - test charge interaction Hamiltonian reads , where is the potential at point created purely by the dielectric. By linear response theory, the average of this potential at some point is changed by the presence of by

(3.12) |

where we suppose that . The basic physical assumption about screening in dielectric media is that only a fraction of the external charge is screened. Consequently, surrounds itself with a polarization cloud of microscopic size carrying charge . For much larger than any microscopic scale, the l.h.s. of (3.12) is thus the potential at due to a pointlike charge at , i.e.

(3.13) |

The charge density is related to the potential via the Poisson equation . Taking first the Laplacian on and then on in both sides of Eq. (3.13), we obtain

(3.14) |

The second-moment sum rule (3.11) follows from this expression after an integration by parts,

(3.15) | |||||

The formula (3.10) provides a split of the integral value onto its classical and purely quantum-mechanical parts. In the quantum regime, it can be used to generate a systematic semiclassical expansion of the second charge-charge moment in powers of . In general, the r.h.s. of (3.10) is not expressible in terms of elementary functions. Perhaps the only exception is the jellium model of conductors whose dielectric function is adequately described, in the long-wavelength limit , by the Drude formula with the dissipation constant taken as positive infinitesimal Jackson (),

(3.16) |

Inserting (3.16) with into (3.10) and using the analog of the summation formula (3.8)

(3.17) |

for , we obtain

(3.18) |

This result coincides with the microscopic finding Martin85 (); Martin86 (). Of course, the result (3.18) could be derived much simpler by inserting the Drude formula for the jellium (3.16) directly into (3.5) and then calculating the ensuing integral by using the Weierstrass theorem, see Eq. (3.19) which follows. The importance of the formula (3.10) consists in the fact that it provides the division of the integral onto its classical and purely quantum-mechanical parts for all media.

### 3.2 Time-dependent case

In the case of the jellium model with the dielectric function (3.16), the Weierstrass theorem

(3.19) |

( denotes the Cauchy principal value) implies

(3.20) |

Substituting this expression into Eq. (3.4), we arrive at

(3.21) |

The time-dependent second-moment sum rule for the unsymmetrized charge-charge density correlation function has been obtained microscopically in Refs. Martin88 (); John93 (),

(3.22) |

It is easy to verify that the symmetrization of this expression leads directly to the result (3.21).

The general formula (3.4) is applicable to any conducting or dielectric medium, but since no other than jellium microscopic results are available we shall omit its analysis.

## 4 Charge-current density correlations

As concerns the charge-current density correlation function, using the spectral representations (2.14) and (2.18) we obtain

(4.1) | |||||

As a check of the validity of this relation, we first multiply its both sides by , then sum over and finally use the spectral version of the continuity equation (2.19) to arrive at (3.1).

The inverse Fourier transform of (4.1) leads to two conditions:

(4.2) |

and

(4.3) |

Note that the sum rule (4.3) fixes the correlation function of the th current component at time and the corresponding th component of the total dipole moment at time . The conditions (4.2) and (4.3) have the same form in both nonretarded and retarded regimes.

The static version of the sum rule (4.3) is trivial. Since is an odd function of , we have

(4.4) |

Let us now consider the time-dependent version of the condition (4.3) for the jellium. Substituting the formula (3.20) into (4.3), we get

(4.5) |

The same result has been obtained, after the symmetrization, by microscopic approaches John93 (); Alastuey89 ().

## 5 Current-current density correlations

For the current-current density correlation function, using the spectral representation (2.18) we obtain after some algebra that

(5.1) |

As a check, multiplying both sides of this relation by , summing over and using the spectral version of the continuity equation (2.19) we arrive at (4.1).

The Fourier component (5.1) has a well defined value at since the two “dangerous” terms of type just cancel with one another at this point, in both quantum and classical regimes. As a result, the current-current density correlation function is integrable in space and its zeroth moment reads

(5.2) |

The spatial integrability of the current-current density correlation function observed here is, specifically in the quantum regime, in contradiction with findings of microscopic theories for particles interacting through Coulomb forces only Martin85 (); Martin88 (); John93 (); a discussion about this subject will be given below in Sect. 6.

### 5.1 Static case

In the static case, the formula (5.2) can be rewritten in the form

(5.3) |

where

(5.4) |

the subtraction of the real sum from has no effect due to the operation Im. The -function possesses the symmetry and vanishes at due to the asymptotic relation (3.6). We can therefore apply the integration techniques together with the general analytic properties of dielectric functions in the complex frequency upper half-plane presented in the second paragraph of Sect. 3.1, see Eqs. (3.7)-(3.10). The final result reads

(5.5) |

where are the (real) Matsubara frequencies.

In the high-temperature (classical) limit , each of the frequencies is much larger than , the corresponding terms in the sum on the r.h.s. of Eq. (5.5) vanish and thus, for all conductor and dielectric media, it holds

(5.6) |

This classical result can be obtained in a simpler alternative way. In the classical regime, the position-dependent vector potential of the EM field can be eliminated from the momentum part of the Gibbs measure as well as from the definition of the current density by the Bohr-van Leeuwen theorem Bohr (); Leeuwen (). Consequently,

(5.7) | |||||

The momentum part of the one-particle Boltzmann factor is proportional to , so that . With regard to the definition of the plasma frequency (3.6) we conclude that

(5.8) |

This result has been obtained for the classical jellium longtime ago by Felderhof Felderhof64 (); Felderhof65 (). The classical sum rule (5.6) is a trivial consequence of the exact relation (5.8). Moreover, higher moments of the current-current density correlation function vanish in the classical limit,

(5.9) |

In the quantum regime, the formula (5.5) can be used to generate a systematic semiclassical expansion of the zeroth moment of the current-current density correlation in powers of . For the jellium conductor with the Drude dielectric function (3.16), using the summation formula (3.17) we obtain explicitly

(5.10) |

This result is not confirmed by microscopic theories for particles interacting only through Coulomb forces Martin85 (); Martin88 (); John93 () which predict the spatial nonintegrability of the current-current density correlation function in the quantum regime (see Sect. 6).

### 5.2 Time-dependent case

## 6 Problematic points in current-current fluctuations

There are certain problematic points in the results for the current-current density correlation function, in the fluctuational theory presented here as well as in microscopic approaches Martin85 (); Martin88 (); John93 ().

We start with an obvious inconsistency having the origin in the formula (5.1). We shall consider the jellium with the dielectric constant (3.16). Using the relation (3.20) and the similar one

(6.1) |

with a modified plasma frequency , the partial inverse Fourier transformation of (5.1) in time leads to

(6.2) |

In the static case and in the classical limit, which corresponds to setting , this formula is equivalent to

(6.3) | |||||

which is an infinite series in the wave-vector components. On the other hand, as follows from the exact analysis between Eqs. (5.6)-(5.9), only the zeroth moment of the classical current-current density correlation function is nonzero, so that it must hold

(6.4) |

Our expression (6.3) coincides with this exact formula only in the leading absolute term which was in the center of interest in Sect. 5. The erroneous presence of higher-order terms is intuitively related to the fact that Rytov’s fluctuational theory works with the local dielectric function which is the long-wavelength limit of the -dependent one , . A complete theory should involve the whole function . We can anticipate the validity of Rytov’s fluctuational theory only in the leading order, in which the -dependence of can be ignored without any approximation. In fact, without saying it, this strategy has already been adopted.

In microscopic approaches for purely Coulomb systems Martin85 (); Martin88 (); John93 (), the obtained result for the (partial) Fourier transform of the current-current density correlation function

(6.5) |

differs from the present one (6.2) in the prefactor to the transverse component which does not depend neither on time nor on . The expression (6.5) works well in the static classical regime, where it reproduces the exact formula (6.4). As soon as the time difference or the regime is quantum , the two “dangerous” terms of type do not cancel with one another in (6.5). Since

(6.6) |

the current-current density correlation function has a nonintegrable algebraic decay at large distances.

The discrepancy between our (retarded) and microscopic (non-retarded) results for the current-current density correlations is probably analogous to that for the induced electric potential and field correlation functions Jancovici06 (); LP (). Our analysis does not exclude an algebraic long-range behavior of the current-current density correlation function, but this decay must be fast enough to ensure the spatial integrability of the correlation function.

## 7 Classical limit of current-current fluctuations

The classical limit of the current-current density correlation function for the jellium was investigated longtime ago by Felderhof, in the nonretarded regime Felderhof64 () as well as for the retarded plasma fully coupled to the classical radiation Felderhof65 (). Changing the notation and using the substitution , his final formula (13.5) in Ref. Felderhof65 () can be rewritten after some straightforward algebra as follows

(7.1) |

where

(7.2) |

and

(7.3) |

with ; the function and the parameter are defined by

(7.4) |

In the limit (), becomes one of the standard definitions of the -function. Consequently,

(7.5) |

where is the dielectric function of the jellium (3.16). Inserting these expressions into (7.1) and comparing the resulting formula with the classical version () of the current-current density correlation function (5.1), we see the equivalence of the two approaches in the considered limit . However, when is finite, the approaches give different results. Felderhof’s classical formula turns out to be correct for an arbitrary value of since it leads to the exact relation (5.8) or its equivalent (6.4) Felderhof65 (). From the discussion in Sect. 6 we know that this is no longer true for the classical version of our formula (5.1). This fact might be a motivation for searching new improvements of Rytov’s fluctuational theory which reproduce classical results for any value of and simultaneously describe adequately correlation functions in the quantum regime.

## 8 Conclusion

In this paper, by using Rytov’s fluctuational theory we have derived a series of known or new sum rules for the bulk correlation functions of charge and current densities in conductor and dielectric media, fully coupled to the radiation field (retarded regime). Few important facts are worth to mention.

In the Fourier space, we can anticipate the validity of Rytov’s macroscopic approach only in the leading order. In this order, the -dependence of the dielectric function can be ignored without any approximation.

In the case of the static sum rules, an analysis in the complex frequency upper half-plane was used to split these sum rules onto their classical and purely quantum-mechanical parts. The classical part coincides, for both conductors and dielectrics, with the results obtained previously by other methods. The quantum-mechanical part was checked on the jellium model of conductor.

As concerns the time-dependent sum rules, their validity was tested on the jellium, too. The retarded results for the charge-charge and charge-current density correlation functions coincide with the non-retarded ones obtained previously by microscopic approaches to the jellium model. The general formulas are applicable to any media, but no other than jellium microscopic results are available to make a comparison.

The most interesting quantity is the current-current density correlation function. Within microscopic theories in the non-retarded regime, this correlation function turns out to be spatially integrable in the classical limit, but it becomes nonintegrable in the quantum case due to an algebraic decay at large distances. This contradicts our retarded prediction (5.2) about a finite value of the integral over the space, in both classical and quantum cases.

A generalization of the present results to inhomogeneous situations, like a surface contact between two distinct media, is left for the future.

###### Acknowledgements.

It is a pleasure to thank Bernard Jancovici for fruitful discussions on the subject of this work and careful reading of the manuscript. I am grateful to LPT for very kind invitation and hospitality. The support received from the European Science Foundation (“Methods of Integrable Systems, Geometry, Applied Mathematics”), Grant VEGA No. 2/0113/2009 and CE-SAS QUTE is acknowledged.## References

- (1) Debye P., Hückel E.: The theory of electrolytes I. Lowering of freezing point and related phenomena. Phys. Z. 24, 185-206 (1923)
- (2) Felderhof B.U.: A theory of fluctuations in plasmas. Physica 30, 1771-1792 (1964)
- (3) Felderhof B.U.: On fluctuations and coherence of radiation in classical and semi-classical plasmas. Physica 31, 295-316 (1964)
- (4) Van Kampen N.G., Felderhof B.U.: Theoretical Methods in Plasma Physics. North Holland, Amsterdam (1967)
- (5) Stillinger F.H., Lovett R.: General restriction on the distribution of ions in electrolytes. J. Chem. Phys. 49, 1991-1994 (1968)
- (6) Pines D., Nozières Ph.: The theory of quantum liquids. Benjamin, New York (1966)
- (7) Martin Ph.A., Oguey Ch.: Dipole and current fluctuations in the quantum one-component plasma at equilibrium. J. Phys. A: Math. Gen. 18, 1995-2016 (1985)
- (8) Martin Ph.A., Oguey Ch.: Screening of classical charges in quantum Coulomb systems. Phys. Rev. A 33, 4191-4198 (1986)
- (9) Martin Ph.A.: Sum rules in charged fluids. Rev. Mod. Phys. 60, 1075-1127 (1988)
- (10) John P., Suttorp L.G.: Equilibrium fluctuation formulae for the quantum one-component plasma in a magnetic field. Physica A 192, 280-304 (1993)
- (11) Alastuey A., Martin Ph.A.: Absence of exponential clustering in quantum Coulomb fluids. Phys. Rev. A 40, 6485-6520 (1989)
- (12) Brydges D.C., Martin Ph.A.: Coulomb systems at low density: a review. J. Stat. Phys. 96, 1163-1330 (1999)
- (13) El Boustani S., Buenzli P.R., Martin Ph.A.: Equilibrium correlations in charged fluids coupled to the radiation field. Phys. Rev. E 73, 036113 (2006)
- (14) Buenzli P.R., Martin Ph.A., Ryser M.D.: Thermal quantum electrodynamics of non relativistic charged fluids. Phys. Rev. E 75, 041125 (2007)
- (15) Jancovici B.: Electric-field correlation in quantum charged fluids coupled to the radiation field. Phys. Rev. E 74, 052103 (2006)
- (16) Høye J.S., Stell G.: Quantum statistical mechanics for polarizable fluids. J. Chem. Phys. 75, 5133-5142 (1981)
- (17) Thompson M.J., Schweizer K.S., Chandler D.: Quantum theory of polarization in liquids: Exact solution of the mean spherical and related approximations. J. Chem. Phys. 76, 1128-1135 (1982)
- (18) Rytov S.M.: Theory of Electrical Fluctuations and Thermal Radiation. Publishing House of AS USSR, Moscow (1953)
- (19) Rytov S.M.: Correlation theory of thermal fluctuations in an isotropic medium. Sov. Phys. JETP 6, 130-140 (1958)
- (20) Levin M.L., Rytov S.M.: Theory of Equilibrium Thermal Fluctuations in Electrodynamics. Science, Moscow (1967)
- (21) Lifshitz E.M., Pitaevskii L.P.: Statistical Physics. Pergamon, Oxford (1981)
- (22) Sitenko A., Malnev V.: Plasma Physics Theory. Kluwer, London (1995)
- (23) Lifshitz E.M.: The theory of molecular attractive forces between solids. Sov. Phys. JETP 2, 73-83 (1956)
- (24) Joulain K., Mulet J.P., Marquier F., Carminati R., Greffet J.J.: Surface electromagnetic waves thermally excited: Radiative heat transfer, coherence properties and Casimir forces revisited in the near field. Surf. Sci. Rep. 57, 59-112 (2005)
- (25) Šamaj L., Jancovici B.: Equilibrium long-ranged charge correlations at the surface of a conductor coupled to electromagnetic radiation. Phys. Rev. E 78, 051119 (2008)
- (26) Jancovici B., Šamaj L.: Equilibrium long-ranged charge correlations at the surface between media coupled to electromagnetic radiation. Phys. Rev. E 80, 031139 (2009)
- (27) Alastuey A.: Breakdown of Debye screening in quantum Coulomb systems and van der Waals forces. Physica A 263, 271-283 (1999)
- (28) Cornu F.: Exact algebraic tails of static correlations in quantum plasmas at low density. Phys. Rev. Lett. 78, 1464-1467 (1997)
- (29) Jackson J.D.: Classical Electrodynamics. Wiley, New York (1975)
- (30) Landau L., Lifshitz E.M.: Electrodynamics of Continuous Media. Pergamon, Oxford (1960)
- (31) Gradshteyn I.S., Ryzhik I.M.: Table of Integrals, Series, and Products. Academic Press, London (2000)
- (32) Jancovici B.: Classical Coulomb systems: screening and correlations revisited. J. Stat. Phys. 80, 445-459 (1995)
- (33) Chandler D.: The dielectric constant and related equilibrium properties of molecular fluids: interaction site cluster theory analysis. J. Chem. Phys. 67, 1113-1124 (1977)
- (34) Alastuey A., Ballenegger V.: Statistical mechanics of dipolar fluids: dielectric constant and sample shape. Physica A 279, 268-286 (2000)
- (35) Ballenegger V., Martin Ph. A.: Dielectric versus conductive behaviour in quantum gases: exact results for the hydrogen plasma. Physica A 328, 97-144 (2000)
- (36) Jancovici B.: private communication
- (37) Bohr N.: Dissertation. Copenhagen (1911)
- (38) Van Leeuwen, H.J.: Problèmes de la théorie électronique du magnétisme. J. Phys. Radium 2, 361-377 (1921).