Thermal Fluctuations and Electromagnetic Noise Spectra in Quantum Statistical Mechanics

We derive the thermal noise spectrum of the longitudinal and transverse electric field operator of a given wave vector starting from the quantum-statistical definitions and relate it to the frequency and wave vector dependent complex conductivity in a homogeneous, isotropic system of electromagnetic interacting charged particles in the frame of the non-relativistic QED. No additional assumptions except the validity of linear response are used in the proof. The Nyquist formula for vanishing frequency, as well as the noise spectral density of Callen-Welton follow as byproduct. Furthermore we discuss also the noise of the photon occupation numbers.


Introduction
There are in the literature several theoretical approaches to the thermal noise in solids. The basic result is Nyquist's law from 1928 [1] for voltage noise at zero frequency, confirmed simultaneously by the experiment of Johnson [2]. Later in 1951 Callen and Welton [3] established a formula relating the integrated voltage noise spectrum to the frequency Ladislaus Alexander Bányai banyai@itp.uni-frankfurt.de Paul Gartner gartner2@gmail.com Mircea Bundaru bundaru@gmail.com 1 dependent resistance based on a fluctuation-dissipation theorem. In 1985 Bányai, Aldea and Gartner Ref. [4] gave a derivation of the Callen-Welton spectral density based directly on the spectrum of the voltage fluctuations within a mean-field approximation of the linear response (as it was actually done also in Ref. [3]), as well as a more exact one taking into account explicitly the Coulomb interactions. In this improved latter case, it turned out that the frequency dependent resistance is replaced by a certain frequency-dependent impedance. Thereafter, derivations of the noise spectral density under different stochastic, thermodynamic electric circuitry approaches with sometimes differing results were published and finally Ginzburg and Pitaevskii in 1988 [5] also confirmed the spectral density of Callen-Welton within a Langevin approach to the electrical circuit, however only with constant resistance R, inductance L and capacity C.
Later review articles on the fluctuation-dissipation theorem [6] and noise [7] are still mean-field theories and Ref. [7] concentrates more on mezoscopic systems.
The results of Ref. [4] went unmentioned by both Ref. [5] and [6,7]. The fluctuation-dissipation theorem was extended byŜpiĉka, Velický and Kalvová [8] to non-equilibrium within the non-equilibrium Green function formalism via the generalized Master equation and derived certain noise spectra for given physical systems.
In a recent book of one of the authors [9] an alternative derivation of the old result of [4] was given, starting from the longitudinal electric field fluctuations. The extension of this approach is the aim of this paper. Within the frame of the non-relativistic QED we correlate the thermal noise spectral densities of longitudinal and transverse electric fields of a given wave vector k with the frequency and wave-vector dependent complex conductivity in a homogeneous, isotropic macroscopic system. We derive also an expression for the thermal noise spectral density of the photon occupation numbers of wave vector k that may be useful for mesoscopic systems.
One starts from the quantum-mechanical definition [10] of the fluctuation of an observable in a system in macro-canonical equilibrium defined as the average quadratic deviation in time. The noise spectrum is just its Fourier transform. We use the definition of the complex frequency and wave-vector dependent dielectric function (or conductivity) by the well-known modification of Kubo's linear response theory [11] in an electromagnetic interacting system due to [12] and [13]. As in [4] we stress the role of the internal electromagnetic fields (i.e. including also the magnetic component too) produced by the particles themselves. We treat ultimately the fluctuations of the local charge density and local current density as true sources of the electromagnetic equilibrium noise. Ref. [4] has already shown for the voltage noise in an isolated sample that this leads to differences with respect to results obtained within the mean-field or the circuitry approaches. Even though quantitatively they may be not so important, differences in the analytic expressions do matter.
Along this line we derive first the thermal noise spectrum of a longitudinal electric field of wave-vector k in a homogeneous isotropic system. We obtain a generalized Nyquist formula relating it for non-vanishing frequency (of the noise) and wave-vector k (of the field) to the frequency and wave-vector dependent longitudinal dielectric function (or conductivity).
Thereafter, we derive the never touched noise spectrum of the transverse electric field (radiation field) starting from the non-relativistic QED by neglecting the retardation (i.e. in the 1 2 approximation). A Nyquist -like relation occurs again. Though, it has certain peculiarities due to the fact that unlike in the longitudinal case, the noise may be related to the conductivity only when the frequency of the noise is related to the wave vector of the field as . Within the same 1 2 approximation of the QED we discuss also the noise spectrum of the photon occupation numbers of wave vector k. This also may be expressed through the electronic transverse current-current correlator. However due to the photonic factors under the integrals it cannot be related to the transverse conductivity. Nevertheless, the formula may be useful for mesoscopic systems.
The only implicit thermodynamic element in our proofs is the expectation, that linear response theory is correct i.e. the thermodynamic limit introduces irreversibility. No supplementary assumptions were needed in deriving these results.

Fluctuations and their Spectral Density
The time fluctuation of a given observable (hermitian operator) in thermal equilibrium is defined [10] as the average square deviation where the average is taken over the macro-canonical equilibrium density matrix 0 (2) and the time evolution of the operator is given by . ( This is analogous to the general definition of a fluctuation both in classical-or quantum statistics. Expanding (1) one has 2 0 2 0 0 Leaving apart the constant 2 0 2 that has a trivial contribution, the time dependent entity of interest is 0 0 which is a real and even function of . As a consequence its Fourier transform 2 0 cos (6) is also real and even. Moreover, according to the Wiener-Khinchin Theorem [14][15][16] it is positive. This is defined [10] as the "noise" spectrum of in quantum statistics 1 . It is easy to show, by expansion in the basis of the eigenfunctions of for any observable (here 0 ) the following two identities Correlators of both types appear in the linear response theory. Our purpose is to relate the electromagnetic thermal noise to the coefficients of that formalism.

Linear Response to an Applied Electromagnetic Field
In a system of electromagnetic interacting charged particles the perturbing Hamiltonian contains only the applied (external) fields. However, the charged particles themselves are sources of electromagnetic fields. This makes the formulation of the problem more delicate. The linear part of the perturbation due to applied electromagnetic scalar x and vector potentials A k is The average internal e.m. (electromagnetic) field has yet to be calculated. Together with the external field it gives rise to the total e.m. field in the system (matter). The conductivity is defined through the relationship between the average current and this total field. This relationship may be given explicitly only in homogeneous, isotropic systems after Fourier transformation in time and space variables.
The average of the current density operator given by the generalization 2 of Kubo's adiabatic linear response theory [11] to the above general electromagnetic perturbation is where 1 2 3 are vector-indices and summation over repeated indices is understood, while is the adiabatic parameter we shall omit for simplicity in the following. (The current density operators in the correlator being already those in the absence of the external field!) This relationship is gauge invariant and valid for any external electric field stemming from the e.m. potentials It may be obtained by linear response to the perturbation (9), which involves its timederivative, followed by the use of the continuity equation and partial integration assuming the vanishing of the correlator at infinite space-and time-distances.
If the e.m. interaction between the charged particles is ignored, one may reinterpret the e.m. field E as the mean-field E . Only within this approximation the kernel of (11) may be related directly to the conductivity of that system.
As it stays (11) actually needs the formulation of the whole problem within the nonrelativistic QED (see for example Ref. [9] for the definition of its Hamiltonian ), not just the usual Coulomb Hamiltonian used in condensed matter theory. This emerges from the necessity to consider the averages of both the longitudinal and transverse electric fields.
However, after the first few steps we shall resort to a simplified approach by neglecting the retardation i.e. in an 1 2 approximation of the QED in which the motion of the electrons and photons are separated. The photons are considered to be free, while the charged particles interact directly by charge-charge and current-current forces. Their e.m. Hamiltonian [17], [18] (here for sake of simplicity just for electrons) is where j x is the transverse part of the current density (10) and the symbol ... means taking the normal product of the operators. Of course, more generally all charged particles (electrons and ions) should be included in the Hamiltonian, as well as their charge and current densities, but we limit here the discussion to the electron dynamics.

Frequency and Wave Vector Dependent Conductivity in a Homogeneous, Isotropic System
In a homogeneous (translation invariant) system one obtains by Fourier transformation in time and coordinate of (11) i.e. one has a local relationship in the Fourier space If the system under consideration is isotropic, then one may separate the longitudinal and transverse parts and It is important to remark here that the frequency and the wave vector in the transverse case are not independent ( ).

The Longitudinal Dielectric Function and Conductivity
In the longitudinal case it is convenient to consider the dielectric function k using the relationship between the charge density and the scalar potential in the Coulomb gauge It is important to have in mind that in the mean-field approximation that considers Coulomb forces only in a self-consistent manner, the longitudinal conductivity is just proportional to the longitudinal current correlator of the non-interacting charges.

The Transverse Conductivity
and therefore with the definition

The Noise Spectrum of the Longitudinal Electric Field
From the hermiticity of the operator E x it follows that after Fourier transforming in the space variable E k E k and we may define two hermitian scalar operators (observables) as the "real" and "imaginary" parts of the longitudinal field operator k E k The noise of any of these observables (omitting the upper indices R,I) according (7) is Here the inverted order of arguments in the notation underlines that k is the wavevector of the electric field, while is the frequency of the noise spectrum.
The Fourier transform of the longitudinal electric field is related by the Poisson equation ( E x  4 x ) to the Fourier transform of the charge density by It is obvious that the quantities of interest are k q with q k. One has, using translation invariance in homogeneous systems with being the volume (under cyclic boundary conditions).
In the case q k the -condition implies k 0, and one is left with the commutator of the total charge with the charge density at the origin. The former is a conserved quantity and therefore it can be taken at any time, in particular at the same time as the local charge, whith the result that their commutator is zero. Therefore only the case q k contributes, and one gets This last result coincides (up to slight differences in the definitions and unit systems) with those obtained earlier in Ref. [4] 3 and [9].
In the 0 limit, both the resistance and the capacity remain finite, thus the contribution of the capacity to the impedance vanishes. As a result, in this limit the old Nyquist theorem [1] 0 0 4 (51) emerges.
In the mean-field version the capacity that stems from the explicit Coulomb interaction is absent and one gets the standard Callen-Welton spectral density. This occurs implicitly if 2 1. 4 It is worth to mention that Ref. [4] starts directly from the definition of the macroscopic voltage where the isolation of the sample (no current flow perpendicular to the boundary cross-section) is explicitly included. Nevertheless, in both derivations the origin of the voltage fluctuation is the fluctuation of the charge density in the system.

The Noise Spectrum of the Transverse Electric Field
Let us consider now the noise of the transverse (radiation) electric field E x in a homogeneous, isotropic sample. There are two possible polarizations of this field defined by the two orthogonal to each other and to the wave vector k unit vectors e k e k ( 1 2). The transverse electric field operator in the non-relativistic QED is given by with the transverse vector potential operator (radiation field) in the Heisenberg picture of the free Hamiltonian being We consider analogously to the treatment of the noise of the longitudinal electric field the two hermitian (scalar) operators constructed from the Fourier transform of the electrical field in the space variable for an arbitrary polarization direction e k . The noise of the observable according to (8)  and analogously for . From this step on we shall neglect the retardation thus are in the frame of the Hamiltonian (13). Then since in an electrically isolated system only the internal electric field fluctuates one may express the transverse electric field directly by the transverse current density where of course k is the projection of the current density on the same chosen polarization direction. Consequently we consider the time evolution also in the 1 2 approximation i.e. in the frame of the Hamiltonian (13).
Thereafter we get, as in the longitudinal case the only surviving contribution The integrated parts do not contribute, and the simplest way to see this is by recalling that fluctuations were originally expressed as a time integral over the whole real axis, see Eq. 6, and the correlations are supposed to vanish at . Alternatively, a direct proof is obtained by using the invariance against space inversion to show that the real part of the integral over is even in time, while the imaginary part is odd. As a consequence, the time integration may be extended to .
At first glance it seems that we found a simple relationship between the noise of the transverse field and the transverse conductivity. However the correct interpretation is more subtle. In the case of the transverse conductivity the variables and are not independent: . Therefore, according to (36) one has just a frequency dependent conductivity and In the transverse noise spectral density is the wave vector of the electric field while is the frequency of the noise. They are independent variables. A direct relationship to the transverse conductivity however exists only along the line . Thus This relation is the analogue of (44) for the noise spectrum of longitudinal fields. Actually our definition and discussion of the noise spectrum of the electric field follows that of Ref. [9] including now the dependence on the wave vector k and brings it farther through the 1 2 approximation for the transverse conductivity.
From (44), (66) and the positivity of the noise spectral density follows also the positivity of the real part of the longitudinal/transversal conductivity ( 0) within the frame of the linear response theory.

The Noise Spectrum of the Photon Occupation Numbers
Beside the field noises one might consider also the noise of the photon occupation numbers. The noise spectral density of the photon occupation numbers k (with an arbitrary chosen polarization not mentioned) is given by In what follows we neglect as before the retardation and ignore consequently terms higher order than 1 2 . Within this approximation we remain only with an approximate Hamiltonian being the sum of the 1 2 e.m Hamiltonian (13) (as well as some other interactions) for electrons and of the free Hamiltonian of the photons both in the averaging over equilibrium as in the time evolution.
Performing the commutations and taking again into account the translation, rotation and reflection invariance in the coordinate space we get after some algebra is the Bose distribution of photons. The above expression resembles (62) for the noise of the transverse electric field, however they differ essentially due to the presence of the photonic factors alongside the electronic current-current correlator under the integrals. Therefore, it cannot be related to the transverse conductivity. Nevertheless, (72) may be the starting point for the direct computation of the photon number noise spectrum within some approximations for homogeneous, isotropic solid state models.
On the other hand, photon number fluctuations may be very important in mesoscopical systems. The above equation, whose derivation does not depend on the thermodynamic limit may serve for their discussion.

Conclusions
We started from the generally accepted quantum statistical definition [10] of the thermal fluctuation of an observable as its average square deviation in time and the noise as its Fourier transform. We have chosen as relevant electromagnetic observables the real and imaginary (hermitian) parts of the Fourier transforms of the longitudinal k and transverse k electric fields in a homogeneous isotropic system. We use throughout the improved linear response theory according to the proper definition of the total field in a system of electromagnetically interacting charged particles according to [12] and [13] (see also Ref. [9]) to relate the thermal noise spectrum of these operators to the frequency and wave-vector dependent complex dielectric function (or conductivity).
Our approach may be considered in the longitudinal case as the quantum-mechanical version of the old analysis of electric fluctuations [19] in the classical plasma theory.
In the case of a longitudinal field in the zero wave vector and zero frequency limit the classical Nyquist theorem [1] results, while for non-vanishing frequencies a modification of the result of Ref. [3] as it was predicted earlier in Ref. [4] emerges. Actually we got as well in the longitudinal as in the transverse case more general results (44) and (66), including 0. However, in the transverse case the relationship to the transverse conductivity occurs obviously only for and it was obtained from the QED within the 1 2 approximation i.e. after neglecting the retardation. These relations may be looked at as generalized Nyquist theorems for the noise spectra of intensive entities.
Finally, we discussed also the noise spectrum of the photon occupation numbers within the same 1 2 approximation. This noise spectrum (72) is again related to the transverse electronic current-current correlator. However, due to the presence of the Bose function under the integrals it cannot be related to the transverse conductivity. Nevertheless, this result may be useful for mesoscopic systems.
These three kind of noise spectra we discussed exhaust the experimental configurations of electromagnetic thermal fluctuations in (macroscopic) condensed matter.
It is important to mention, that in our derivations no other ingredients were introduced except the implicit assumption of linear response theory about the time-decay of correlations (i.e. irreversibility in the thermodynamic limit).
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.