Oct 17, 2017 - to photo-detection amplitudes in the vacuum state, they do effect the ... Together with the laws of space and time the quantum ... cond...

0 downloads 0 Views 1010KB Size

arXiv:1710.06207v3 [quant-ph] 16 Jan 2018

January 17, 2018

The solution to the wave equation as a Cauchy problem with prescribed fields at an initial time t = 0 is purely retarded. Similarly, in the quantum theory of radiation the specification of Heisenberg picture photon annihilation and creation operators at time t > 0 in terms of operators at t = 0 automatically yields purely retarded source-fields. However, we show that twotime quantum correlations between the retarded source-fields of a stationary dipole and the quantum vacuum-field possess advanced wave-like contributions. Despite their advanced nature, these correlations are perfectly consistent with Einstein causality. It is shown that while they do not significantly contribute to photo-detection amplitudes in the vacuum state, they do effect the statistics of measurements involving the radiative force experienced by a point charge in the field of the dipole. Specifically, the dispersion in the charge’s momentum is found to increase with time. This entails the possibility of obtaining direct experimental evidence for the existence of advanced waves in physical reality, and provides yet another signature of the quantum nature of the vacuum.

1 Introduction Electrodynamics describes the production, propagation and absorption of electromagnetic waves. The problem most commonly encountered in classical electrodynamics is that in which the source is specified at some initial time, and then the electromagnetic fields are sought at later times. The unique solution to this problem is fully retarded, which means that the electromagnetic source-fields depend in a causal way on the state of the source in the past [1, 12, 36]. Advanced source-fields, which depend on the source in the future are obtained when Maxwell’s equations are solved running into the past [1, 31]. These solutions are not viewed as physical within the context of an initial value problem. Accepted in

Quantum 2017-12-15, click title to verify

It is well known that in quantum optics the commutators of the free electromagnetic fields involve both retarded and advanced green’s functions for the waveoperator [6]. On the other hand the subject of advanced waves in the presence of sources has received more limited attention in quantum theory. The reason for this may be that the problem of finding the Heisenberg picture source-fields at time t > 0 in terms of Schr¨ odinger picture operators at time t = 0 is an initial value problem whose solution entails purely retarded source-fields. However, due to differences between quantum and classical theories the complete absence of advanced contributions in quantum optics is not immediate. The fluctuating quantum vacuum offers one of the most striking examples of such a difference. The theories also differ in their treatment of electromagnetic correlations. The quantum vacuum is used to interpret various phenomena in quantum optics including spontaneous emission [18], atomic energy-shifts [24, 35], and Casimir forces [5]. Together with the laws of space and time the quantum theory of the vacuum predicts intriguing physical phenomena such as Hawking radiation associated with black holes [11] and the Unruh-Davies effect whereby an accelerated detector in the vacuum of Minkowski spacetime registers thermal excitations [8, 34]. The Unruh-Davies effect and related effects [20] can be understood in terms of correlations between the quantum vacuum at timelike separated points in spacetime. Vacuum induced correlations across spacelike separations have also received attention [25]. In atomic and optical physics many effects attributed to the quantum vacuum can often also be understood, at least partially, in terms of radiation reaction [7, 13, 17]. Thus, identifying specific signatures associated with the quantum vacuum remains of broad interest in quantum optics. Recently, direct observation of vacuum fluctuations has been achieved [19, 26, 27], which opens up new possibilities for investigating the interplay between vacuum fluctuations and the electromagnetic fields produced by charged sources. In this paper we consider the fields produced by a quantum dipole. We show that the presence of the 1

quantum vacuum and the distinct nature of quantum correlations imply that certain physical predictions involve advanced waves. More specifically, we show that the combination of purely retarded source-fields and quantum vacuum fields result in advanced-wave contributions to two-time quantum correlations. There are seven sections to this paper. We begin in section 2 by briefly reviewing the interpretation of waves in classical physics. In section 3 we introduce a model of the dipole-field system and derive the electromagnetic source-fields of the dipole. In section 4 we discern the role of interference between vacuum and source-fields within the process of spontaneous emission. In section 5 we consider general two-time quantum correlations functions of the electromagnetic field, and show that vacuum-source correlations receive advanced-wave contributions. Our subsequent aim is to determine whether such correlations might be physically measurable. In section 6 we briefly verify that vacuum-source correlations do not significantly contribute to photo-detection amplitudes. In section 7 we consider the radiative force experienced by a test charge in the field of a quantum dipole. By considering a charge with relatively slow dynamics, we are able to develop a quantum description of the free charge, which is analogous to the multipolar description of bound charge systems. This description enables us to properly identify the relevant vacuum and source-fields. We then show that vacuum-source correlations make signficant contributions to the statistical predictions involving the force experienced by the test charge. Finally in section 8 we summarise our findings.

2 Classical solutions of the wave equation We begin by briefly reviewing the interpretation of solutions to the wave equation in classical physics. This interpretation may be contrasted with the interpretation given to the advanced-wave like contributions to quantum correlations that are reported in subsequent sections. The inhomogeneous wave equation for a scalar field ψ over Minkowski spacetime is ∂t2 − ∇2 ψ(t, x) = f (t, x), where we have assumed units such that the speed of the wave is one and where f (t, x) describes the sources of the waves. A well-posed problem is one in which boundary conditions and the source data f are specified, along with the fields ψ(t0 , x) and ψt0 (t0 , x) = ∂t0 ψ(t0 , x) at some time t0 [1]. The fields at some other time t are then sought. Methods of solving the wave equation in classical electrodynamics invariably involve green’s functions for the wave-operator ∂t2 − ∇2 [12, 31, 36]. The retarded and Accepted in

Quantum 2017-12-15, click title to verify

advanced green’s functions G+ and G− admit the representations G± (t, x|t0 , x0 ) =

1 δ (t0 − t ± |x − x0 |) . 4π|x − x0 |

(1)

In conjunction with either initial or final conditions the wave equation is easily solved in reciprocal space [31]. This method of solution is closest to that used in quantum theory. The quantum electric field associated with a point source within the Coulomb gauge is found using essentially this method in appendix 9.4. If initial conditions are specified such that before any sources are activated at t = 0 there is a known input configuration ψ in , which satisfies the homogenous equation ∂t2 − ∇2 ψin (t, x) = 0, the field at time t > 0 is given by [31] ψ(t, x) =ψin (t, x) Z Z t 3 0 + d x dt0 f (t0 , x)G+ (t, x|t0 , x0 ).

(2)

0

Significantly, the advanced green’s function G− does not contribute to the solution for t > 0 of the initial value problem. Moreover, when initial conditions are specified for the source function f , the homogeneous component ψin can be taken to vanish [36]. Therefore, in classical electrodynamics only the retarded source-field plays a role in the physical predictions associated with a source that is initially fully specified. The retarded solution is replaced by the advanced solution when final conditions are specified instead of initial conditions, and the Maxwell equations are solved running into the past [1, 12, 31]. In contrast to the classical situation, in quantum theory, even when all sources are fully specified the vacuum plays an important role in physical phenomena [17]. We will show further that although all sourcefields remain fully retarded, certain predictions feature advanced waves associated with the quantum vacuum.

3 Dipole-field Hamiltonian and electric and magnetic source-fields 3.1 Dipole-field Hamiltonian Consider a single atom coupled to the electromagnetic field. Within the electric dipole approximation the Hamiltonian describing this system is H = H0 + V with Z X X † 1 3 H0 = ωn |ni hn| + d k ω aλ (k)aλ (k) + 2 n λ

= Hd + Hf ˆ · DT (0) V = −d

(3) 2

where |ni denotes an energy eigenstate of the dipole with associated energy ωn , aλ (k) denotes the annihilation operator for a photon with momentum k and polarˆ denotes the dipole moment operator, and isation λ, d DT denotes the transverse electric displacement field. The dipole is assumed to be located at the origin 0. In terms of photonic operators aλ (k) the transverse displacement field admits the mode expansion

3.2 Electric and magnetic source-fields The equation of motion for the photon annihilation operator is found using the Hamiltonian (3) and once formally integrated reads aλ (t, k) = aλ (k)e−iωt Z t r 0 −iω(t−t0 ) + dt e 0

DT (t, x) =

(+) DT (t, x)

Z

(+)

DT (t, x) = i =

d3 k

(−) DT (t, x),

+ Xr

λ (−) DT (t, x)†

ω eλ (k)aλ (t, k)eik·x 2(2π)3 (4)

where the e1,2 (k) are mutually orthogonal unit polarisation vectors orthogonal to k, and ω = |k|. The vacuum (free) components of the electric field and transverse displacement field are identical for all x and are given by (+)

(+)

DT,0 (t, x) = E0 (t, x), Z Xr ω 3 eλ (k)aλ (k)e−iωt+ik·x . =i d k 2(2π)3 λ

(5) For x 6= 0 the source components also coincide DT,s = Es . In terms of photonic operators the electric field for x 6= 0 is therefore given by DT (t, x) = E(t, x) = E(+) (t, x) + E(−) (t, x) (+)

E(+) (t, x) = DT (t, x).

(6)

The situation is different at the position of the dipole x = 0. This is most easily seen by noting that D(x) − E(x) = P(x) where P(x) = dδ(x) in the electric dipole approximation. Thus, at the position of the dipole E differs from D by an infinite field. Since radiation is detected by an absorbing atom at the position of the atom, it is important to recognise the distinction between the electric field and the displacement field at this position. The magnetic field admits the modeexpansion B(t, x) = B(+) (t, x) + B(−) (t, x), B(+) (t, x) Z Xr = i d3 k λ

ω ˆ k × eλ (k) aλ (t, k)eik·x 2(2π)3

(7)

and similarly to the transverse displacement field it can be split into vacuum and source components via the corresponding partitioning of aλ (t, k). Accepted in

Quantum 2017-12-15, click title to verify

ω ˆ 0 ). eλ (k) · d(t 2(2π)3 (8)

Away from the dipole, the electric and magnetic dipole source-fields are found using Eqs. (6) and (80) and are exactly analogous to classical fields of an oscillating dipole; i o 1 n h ˆ ˆ r) ˆ · d(tr ) − d(t 3ˆ x x Es (t, x) = 3 4πx i o 1 n h ˆ˙ ˆ˙ r ) ˆ · d(tr ) − d(t + 3ˆ x x 2 4πx o i 1 n h ¨ˆ ¨ˆ ˆ x ˆ · d(tr ) − d(t x + r) , 4πx 1 1 ¨ˆ ˆ˙ r ) ˆ × d(t ˆ × d(t Bs (t, x) = − x x (9) r) − 4πx 4πx2 where tr = t − x > 0. For t < x the source-fields vanish. The source-fields can be split into radiative and nonradiative components. The radiation fields are responsible for irreversible energy loss from the dipole. They vary as x−1 , are orthogonal to x and depend on the dipole moment acceleration. They are related via Erad,s (t, x) = −ˆ x × Brad,s (t, x). The remaining parts of the source-fields are non-radiative. The non-radiative electric field includes both velocity-dependent and static Coulomb-like parts whereas the non-radiative magnetic field is entirely velocity-dependent. (±) The creation and annihilation source-fields Es are not causal, that is, they do not generally vanish for t < x. However, as in references [4, 16] we can perform both a rotating-wave approximation and an associated Markov approximation, which yield for t > x E(+) s (t, x) = X iωnm 1 {3ˆ x (ˆ x · dnm ) − dnm } + 4πx2 4πx3 n,m n

2 ωnm ˆ (ˆ {dnm − x x · dnm )} σnm (tr ), + 4πx 2 X ωnm iωnm ˆ B(+) (t, x) = x × d − σnm (tr ) nm s 4πx 4πx2 n,m n

(10) where σnm (t) = eiHt |ni hm| e−iHt . For t < x the creation and annihilation fields vanish within the approximations made. These approximations also ensure that 3

the positive and negatve frequency fields are associated with dipole lowering and raising operators respectively. Eq. (10) can be further simplified through the approximation σnm (t) ≈ σnm (0)eiωnm t . The latter is appropriate for perturbative evaluations of quadratic field functionals correct to second order. The steps leading to the final result in Eq. (10) are given in appendix 9.1. In later sections we we will focus on the case of a twolevel dipole with ground and excited states denoted |gi and |ei respectively. The two-level transition frequency is denoted ω0 , and the transition dipole moment, which is assumed to be real, is denoted d. We also introduce the two-level operators σ + = |ei hg| ,

σ − = (σ + )† ,

σ z = [σ + , σ − ] (11)

in terms of which the source-fields can be written ∗ − E(+) s (t, x) = E (x)σ (tr ), ∗ − B(+) s (t, x) = B (x)σ (tr )

wherein for convenience we have defined 1 iω0 + {3ˆ x (ˆ x · d) − d} E(x) = 4πx2 4πx3 ω2 ˆ (ˆ + 0 {d − x x · d)} , 4πx ω02 iω0 ˆ × d. − B(x) = x 4πx 4πx2

(12)

the radiated power due to the spontaneous photon emission of an initially excited dipole. The S-matrix element representing the probability amplitude for the spontaneous emission process |e, 0i → |m, kλi with m < e and |kλi the state of a single photon, can be calculated easily using second order perturbation theory. The associated rate at which this process occurs is found using Fermi’s golden rule to be Γem =

3 ωem |dem |2 3π

(14)

where ωem = ωe − ωm and dem = he| d |mi. The Smatrix consists of probability amplitudes for processes that conserve the bare energy H0 . Thus, in spontaneous emission the energy of the resulting photon is precisely the energy of the atomic transition through which it was emitted. The energy flux in the process |e, 0i → |m, kλi , m < e is therefore ωem Γem , which implies that the total radiated power is to second order given by X P = ωem Γem . (15) m

(13)

We remark finally that it is not obvious that the equal-time commutation relations between the fields whose source components are given in Eq. (12) are the same as those for the fields in Eq. (6), because approximations have been made. That the commutation relations do remain the same is however necessary in order that the resulting theory is formally consistent. We prove that the equal-time commutators are in fact preserved in appendix 9.3.

The above derivation of P relies entirely on the calculation of matrix elements of the interaction Hamiltonian, so it does not involve using any electromagnetic fields away from the source at 0. In fact, since DT (0) 6= E(0), the electric and magnetic fields do not directly feature anywhere in the above calculation.

4.2 Derivation using electromagnetic fields An alternative expression for the radiated power that does involve the electromagnetic fields, can be obtained from Poynting’s theorem. The energy flux radiated by the dipole located at the centre of a sphere with radius x is given by Z Z 2 P = lim dΩ x I(t, x) = dΩ x2 Irad (t, x) (16) x→∞

4 Vacuum-source correlations and radiation intensity Before considering general quadratic functionals (correlation functions) of the electromagnetic fields we consider in detail a particular example - the radiation intensity. Our aim is to identify the contribution made by the quantum vacuum, and in particular, the role played by vacuum source-field interference. The treatment in this section extends previous perturbative treatments found in [21, 28, 29, 32].

4.1 Radiated power In order to identify the radiation intensity, which is the radiated power per unit area, we begin by calculating Accepted in

Quantum 2017-12-15, click title to verify

where the integration is performed over all directions. Eq. (16) is definitive of the radiation intensity Irad , which is the radiated power per unit area. The intensity component Irad is a quadratic functional of the radiation fields, which vary as x−1 . We now proceed to determine the radiated power using the fields in section 3.1. In order to compare the result of Eq. (16) with Eq. (15), which was obtained using second order perturbation theory we make use of Eq. (10) coupled with the perturbative approximation σnm (t) ≈ σnm (0)eiωnm t . First let us suppose that the intensity Irad in Eq. (16) is of the Glauber type (−)

(+)

IG (t, x) = hErad (t, x) · Erad (t, x)i0,e (−)

(+)

= hErad,s (t, x) · Erad,s (t, x)ie;0

(17)

4

where the average is taken in the state |e, 0i, and we have retained only the radiative contribution dependent on Erad = E0 + Erad,s . Letting σnm (t) ≈ σnm (0)eiωnm t in Eq. (10) we find that 2 X 1 ˆ (ˆ IG (t, x) = ω 4 [dem − x x · dem )]2 . (18) 4πx m

4.3 General definition of the radiation Intensity In order to motivate a general definition of the radiation intensity let us first consider the case of a free field confined to a volume V . The vacuum energy Hvac and the associated density Hvac /V are not measurable. Only differences in energy-density from the vacuum value are measurable. The integral of the measurable local energy-density over all space should therefore yield the energy difference Hf − Hvac [3, 30]. Quite generally, the vacuum energy density can be expressed as Hvac (+) (−) = hE0 (t, x) · E0 (t, x)i0 V (+) (−) = [E0,i (t, x), E0,i (t, x)] (+)

(−)

(20)

Thus, in the general case in which sources are present the average electromagnetic energy-density can be written Uem (t, x) − Uvac ≡ 1 (+) (−) E(t, x)2 + B(t, x)2 − [E0,i (t, x), E0,i (t, x)]. (21) 2 (±)

Since the total fields E(±) and the vacuum fields E0 satisfy the same equal-time commutation relations (cf appendix 9.3), (+)

(−)

(+)

[E0,i (t, x), E0,j (t, x0 )] = [Ei (+)

= [Bi

(−)

(t, x), Bj

(−)

(t, x), Ej

(+)

(t, x0 )]

(−)

(t, x0 )] = [B0,i (t, x), B0,j (t, x0 )], (22)

for the total energy density U = Uem + Umatter , where Umatter is the energy density of sources, we have U (t, x) − Uvac = : U (t, x) :

Accepted in

Quantum 2017-12-15, click title to verify

S(t, x) =

1 [E(t, x) × B(t, x) + H.c.] . 2

(23)

(24)

Analogously to the vacuum energy-density we define the vacuum Poynting vector by Svac,i (t, x) (+)

(−)

= hS0,i (t, x)i0 = ijk [E0,j (t, x), B0,k (t, x)] (+)

= ijk [Ej

(−)

(t, x), Bk (t, x)]

(25)

where ijk denotes the Levi-Civita symbol. Since Uvac and Svac are independent of time and space Poynting’s theorem can be written in terms of : U (t, x) : = U (t, x)− Uvac and : S(t, x) : = S(t, x) − Svac as ∂t h: U (t, x) :i = −∇ · h: S(t, x) :i.

(26)

The practice of normally-ordering operator products in this way constitutes a basic application of Wick’s theorem. For a dipole at 0 the radiated power can be obtained by integrating the right-hand-side of Eq. (26) over a sphere centered at 0 with radius x. Using the divergence theorem and subsequently taking the limit x → ∞ the associated intensity is then seen to be ˆ · hSrad (t, x) − Svac i, Irad (t, x) = x 1 Srad = [Erad × Brad + H.c.] 2

Uvac =

= [B0,i (t, x), B0,i (t, x)].

where : · : denotes normal-ordering of the creation and annihilation operators within Uem . By Poynting’s theorem the local energy-flux can be written in terms of the Hermitian Poynting vector

(27)

where Erad (t, x) = E0 (t, x) + Erad,s (t, x) and Brad (t, x) = B0 (t, x)+Brad,s (t, x). Assuming the initial state |e, 0i and making use of the identities ˆ · Erad,s (t, x) = 0 = x ˆ · Brad,s (t, x), x ˆ × Erad,s (t, x) = Brad,s (t, x), x ˆ × Brad,s (t, x) = −Erad,s (t, x), x B2rad = E 2rad

(28)

it is straightforward to show that the radiation intensity can also be written Irad (t, x) = hErad (t, x)2 − E0 (t, x)2 ie;0 = 2IG (t, x) (29) where the second equality follows within the rotatingwave approximation. This explains why 2IG rather than IG yields the correct radiated power. The anti-normally (+) (−) ordered combination Erad (t, x)·Erad (t, x), which is such (+) (−) (−) (+) that h: Erad (t, x)·Erad (t, x) :i = hErad (t, x)·Erad (t, x)i, is slowly-varying (resonant), so it is not eliminated by the rotating-wave and Markov approximations. 5

4.4 Vacuum source-field correlations

found to second order to be (appendix 9.2)

We now turn our attention to the contribution of vacuum source-correlations to the radiation intensity. In references [21, 28, 29, 32] Eq. (27) is used to perturbatively calculate the radiated power. We provide a simpler derivation of the same result, and in section 4.5 we provide a non-perturbative extension. The physical intensity is determined relative to the vacuum value, so assuming the initial state |e, 0i the intensity can be written according to Eq. (29) as Irad (t, x) =hErad,s (t, x)2 ie,0 + hE0 (t, x) · Erad,s (t, x) + Erad,s (t, x) · E0 (t, x)ie,0 .

(+)

(−)

where all dipole levels m are included within the summation. It remains to calculate the contribution of the vacuum-source correlations on the first line in Eq. (30). This can be done a number of ways, with varying complexity. The authors [21, 28, 29] use the full sourcefield soultions of the electric and magnetic fields and separate out time-independent and oscillatory contributions to the power. The time-independent contributions give our final result Eq. (34) below, while the time-dependent contributions vanish on average, and can be interpreted as virtual transients. Reference [32] gives a simplified derivation of Eq. (34), which involves using only the radiation source-fields rather than the full source-fields. Here we give a yet simpler derivation which makes use of the Markov and rotating-wave approximations. The latter eliminate virtual photon contributions, and therefore yield Eq. (34) directly, without requiring a time-average. The vacuum-source correlation terms within (30), are Quantum 2017-12-15, click title to verify

in agreement with references [21, 28, 29, 32]. The corresponding contribution to the radiated power is Z Pvac−s = dΩ x2 hE0 (t, x) · Erad,s (t, x)i + c.c. =

1X sgn(ωem )ωem Γem . 2 m

(34)

By adding the two contributions Eqs. (32) and (34) one obtains the correct result Eq. (15). Vacuum-source correlations are seen to be responsible for supplying half of the total radiated power, as well as for eliminating, in the final result, the contributions of energy nonconserving virtual transitions for which m > e.

(31)

to second order in the transition dipole moments. By substituting Eqs. (18) and (31) into Eq. (16) we obtain, within the rotating-wave approximation, the contribution of E2rad,s in Eq. (30) to the radiated power to second order as Z 1X ωem Γem (32) Ps = dΩ x2 hErad,s (t, x)2 i = 2 m

Accepted in

(−)

(33)

(30)

First we calculate the anti-normally ordered combina(+) (−) tion hErad,s (t, x) · Erad,s (t, x)ie;0 , which is easily found using Eq. (10) to be hErad,s (t, x) · Erad,s (t, x)ie;0 = 2 X 1 ˆ (ˆ ω 4 [dem − x x · dem )]2 4πx m>e em

(+)

hE0 (t, x) · Erad,s (t, x)i0;e + c.c. = 1 X 4 ˆ (ˆ sgn(ωem )ωem [dem − x x · dem )]2 (4πx)2 m

4.5 Non-perturbative extension The results of section 4.4 can be extended beyond the use of second order perturbation theory. For this purpose we restrict our attention to a two-level emitter. Judicious use of the Markov and rotating-wave approximations within the Heisenberg equations, leads to a tractable system of coupled Heisenberg-Langevin equations for the dipole operators, which can then be formally integrated. We begin by using the solution (80), and the rotatingwave and Markov approximations to calculate the dipole’s own reaction field as (+)

Z

DT,s (t, 0) = i

d3 k

r

ω X eλ (k)aλ,s (t, k) 2(2π)3 λ

ˆ Γd =i σ − (t). 2d

(35)

Substituting Eq. (100) into the Heisenberg equation for each of σ + and σ z , making the rotating-wave approximation, and formally integrating the result yields the solutions Z t 0 + σ + (t) = e(iω0 −Γ/2)t σ + + dt0 e(iω0 −Γ/2)(t−t ) σvac (t0 ) 0

(36) σ z (t) + 1 = e−Γt [σ z + 1] +

Z

t

0

z dt0 e−Γ(t−t ) σvac (t0 ).

0

(37) 6

where + σvac (t)

= i[d ·

(−) DT,0 (t, 0)]σ z (t) (−)

z σvac (t) = −2i[d · DT,0 (t, 0)]σ − (t) + H.c..

(38) (39)

In addition to the solutions to the Bloch equations (without driving), the solutions (36) and (37) also contain contributions from the quantum vacuum field. Using Eqs. (36), (37) and (12) we find the Glauber intensity function in the state |0, ei to be 2 2 ω0 ˆ (ˆ [d − x x · d)] e−Γtr (40) IG (t) = 4πx The corresponding power is 1 (41) PG (t) = ω0 Γe−Γtr 2 whose short-time limit coincides with the power in Eq. (19) restricted to two levels. Next, following the same method as in section 4.4 we find the anti-normally ordered combination (+) (−) hErad,s (t, x) · Erad,s (t, x)ie;0 to be (+)

(+)

(42)

Quantum 2017-12-15, click title to verify

≈ hXi (t, x)Yj (t0 , x0 ) − [Xi

(t0 , x0 )i (−)

(t, x), Yj

(t0 , x0 )]i. (47)

In the vacuum state the commutator (+) (−) [Xi (t, x), Yj (t0 , x0 )], which is subtracted from (

hXi (t, x)Yj t0 , x0 )i in Eq. (47) includes the pure vacuum value hX0,i (t, x)Y0,j (t0 , x0 )i0 as well as sourcedependent components. If instead we subtract only the vacuum value from hXi (t, x)Yj (t0 , x0 )i we obtain the alternative correlation function

(+)

(44)

(45)

(46)

whose short time limit coincides with the power in Eq. (15) for a two-level system. These results mirror the previous perturbative results in that neither Eq. (41) nor Eq. (43) yield the correct radiated power. Furthermore, since the pure source power in Eq. (43) does not decay, it is unphysical for long times when considered on its own. As with the previous perturbative results only the combination of the individually non-physical contributions (43) and (45), yields the correct physical power. Accepted in

(+)

(t, x)Yj

CXi Yj (t, x|t0 , x0 )

By adding Eqs. (43) and (45) we obtain the total radiated power P (t) = ω0 Γe−Γtr = 2PG (t)

(−)

(+)

(−)

The corresponding contribution to the power is 1 −Γtr Pvac−s (t) = ω0 Γ e − . 2

(−)

hXi (t, x)Yj (t0 , x0 )i are taken as rapidly oscillating and are negligible. The corresponding Glauber-type (−) (+) correlation function hXi (t, x)Yj (t0 , x0 )i is therefore related to hXi (t, x)Yj (t0 , x0 )i by (−)

which when combined with Eq. (40) gives the total contribution of the source-field to the radiated power as 1 Ps (t) = ω0 Γ. (43) 2 Finally, the vacuum-source correlations are found using Eqs. (36), (37), (12) and (5) to be (appendix 9.2) hE0 (t, x) · Erad,s (t, x)i0,e + c.c. 2 2 ω0 ˆ (ˆ = 2e−Γtr − 1 . [d − x x · d)] 4πx

In section 4 we showed that a specific quadratic field functional, namely, the intensity functional is the difference between the quantum version of the classical definition and the quantum vacuum value. Accounting for vacuum-source correlations was essential in obtaining the correct physical results. The most general quadratic field functionals are correlation functions defined in terms of the components of the electromagnetic field tensor at different points in spacetime. In a designated (laboratory) inertial frame all such correlations have the form hXi (t, x)Yj (t0 , x0 )i where X, Y = E, B. In the rotating-wave (+) (+) approximation the terms hXi (t, x)Yj (t0 , x0 )i and

GXi Yj (t, x|t0 , x0 ) = 2hXi

(−)

hErad,s (t, x) · Erad,s (t, x)ie;0 = 2 2 ω0 ˆ (ˆ [d − x x · d)] (1 − e−Γtr ), 4πx

5 Vacuum-source correlations and advanced waves

(−)

= hXi (t, x)Yj (t0 , x0 ) − [X0,i (t, x), Y0,j (t0 , x0 )]i. (48) For t 6= t0 the commutators between the total fields E(±) and B(±) are not equal to their counterparts defined in (±) (±) terms of the vacuum fields E0 and B0 . As a result even within the rotating-wave approximation GXi Yj and CXi Yj do not generally coincide. Unlike GXi Yj the vacuum-fields contribute to CXi Yj within the vacuum state. We define the difference operator (+)

∆Xi Yj (t, x|t0 , x0 ) =[Xi

(−)

(t, x), Yj

(t0 , x0 )]

(+)

(−)

− [X0,i (t, x), Y0,j (t0 , x0 )]

(49)

in terms of which, in the rotating-wave approximation, we have CXi Yj = GXi Yj + h∆Xi Yj i.

(50)

A perturbative calculation of certain correlation functions of the form given in Eq. (48) in the particular 7

case that t = t0 has been carried out without using the rotating-wave or Markov approximations within the non-relativistic QED framework by Power and Thirunamachandran [21, 22] (see also [28, 29, 32]). Their results are consistent with the results presented in section 4.4, and differ only by the inclusion of oscillatory transient contributions. The results in references [21, 22] are also consistent with the short-time limiting behaviour of the non-perturbative results given in section 4.5. Here we consider the the more general case in which t and t0 are arbitrary, which is fundamentally different, because ∆Xi Yj (t, x|t0 , x0 ) is generally non-vanishing. While we utilise the rotating-wave and Markov approximations our approach is not perturbative and is similar to that used in section 4.5. More precisely, we make use of the Heisenberg-Lengavin equations for a two-level dipole. Details of the calculations leading to the results below are given in appendix 9.3. (−) Using Eq. (12) we write Xs (t, x) = X (x)σ + (t − x) where X = E, B and correspondingly X = E, B. The GXi Yj are easily found to be GXi Yj (t, x|t0 , x0 ) = 2Xi (x)Yj∗ (x0 )θ(tr )θ(t0r )hσ + (tr )σ − (t0r )i.

0

0

(52)

Next we consider ∆Xi Yj , which can be partitioned into three terms as (+)

(−)

∆Xi Yj (t, x|t0 , x0 ) = [Xs,i (t, x), Ys,j (t0 , x0 )] (−)

(+)

(−)

+ [X0,i (t, x), Ys,j (t0 , x0 )] + [Xs,i (t, x), Y0,j (t0 , x0 )]. (53) The pure source component is easily found using Eq. (12) to be

(0, 0) Figure 1: The world line of the dipole, represented by the central vertical arrow, passes through the intersection of the backward lightcone of the event (t0 , x0 ) and the forward lightcone (+) (−) of (t, x). The commutator [X0,i (t, x), Ys,j (t0 , x0 )] contains a component that depends on the source operator σ − (ta ) where ta is the advanced time associated with the event (t, x) at (+) which the vacuum field X0,i (t, x) is evaluated.

where ta = t + x and t0a = t0 + x0 , and the notation ± means that + holds for the case X = E while − holds when X = B. Similarly we find that (−)

= Xi∗ (x)Yj (x0 )θ(tr )θ(t0r )[σ − (tr ), σ + (t0r )].

(54)

The remaining components involve both the vacuum and source-fields, and they each consist of both retarded and advanced parts. After fairly lengthy calculations we obtain the compact expression (appendix 9.3) (−)

[X0,i (t, x), Ys,j (t0 , x0 )] = − Xi∗ (x)Yj (x0 )θ(t0r − tr )θ(tr )θ(t0r )[σ − (tr ), σ + (t0r )] ± Xi (x)Yj (x0 )θ(t0r − ta )θ(t0r )[σ − (ta ), σ + (t0r )], Quantum 2017-12-15, click title to verify

− Xi∗ (x)Yj (x0 )θ(tr − t0r )θ(tr )θ(t0r )[σ − (tr ), σ + (t0r )] ± Xi∗ (x)Yj∗ (x0 )θ(tr − t0a )θ(tr )[σ − (tr ), σ + (t0a )]

(56)

where + holds for the case Y = E while − holds when Y = B. The sum of the purely retarded parts from Eqs. (55) and (56) exactly cancels the pure source component given in Eq. (54). We therefore obtain the final result ∆Xi Yj (t, x|t0 , x0 ) = (−1)α Xi (x)Yj (x0 )θ(t0r − ta )θ(t0r )[σ − (ta ), σ + (t0r )]

(−)

[Xs,i (t, x), Ys,j (t0 , x0 )]

Accepted in

(t, x)

(+)

2Xi (x)Yj∗ (x0 )θ(tr )θ(t0r )e−[Γ/2](tr +tr ) eiω0 (tr −tr ) .

(+)

ta

[Xs,i (t, x), Y0,j (t0 , x0 )] =

GXi Yj (t, x|t0 , x0 ) =

(+)

(t0 , x0 ) t0r

(51)

where tr = t − x and t0r = t0 − x0 and where θ denotes the Heaviside step-function. Assuming the initial state |e, 0i we obtain

(+)

time

(55)

+(−1)β Xi∗ (x)Yj∗ (x0 )θ(tr − t0a )θ(tr )[σ − (tr ), σ + (t0a )]. (57) where α = 0 for X = E and α = 1 for X = B while β = 0 for Y = E and β = 1 for Y = B. The Heaviside functions in this expression imply that ∆Xi Yj (t, x|t, x0 ) = 0 as required. However, for t 6= t0 , h∆Xi Yj (t, x|t0 , x0 )i generally consists of non-vanishing two-time correlations, which are due to the advanced components of the commutators of the free and sourcefields in Eqs. (55) and (56). Assuming the initial state 8

perturbative approximation σdz (t0 ) ≈ σdz (0) in Eq. (60), which for a detector in the ground state yields Z t 0 p˙d (t) =di dj dt0 hDT,i (t, x)DT,j (t0 , x)iψ e−iω0 (t−t )

|e, 0i, with some work we obtain (appendix 9.3) h∆Xi Yj (t, x|t0 , x0 )ie,0 = (−1)α Xi (x)Yj (x0 )θ(t0r − ta )θ(t0r )θ(ta )

0

0

× e(iω0 −Γ/2)tr e(−iω0 −Γ/2)ta (eΓta − 2)

+ c.c. .

+ (−1)β Xi∗ (x)Yj∗ (x0 )θ(tr − t0a )θ(tr )θ(t0a ) 0

0

× e(−iω0 −Γ/2)tr e(iω0 −Γ/2)ta (eΓta − 2).

(58)

Since ∆Xi ,Yj (t, x|t0 , x0 ) is non-local in both space and time the presence of advanced contributions is not acausal. Supposing for definiteness that t0 > t then all fields at events within the backward lightcone of (t0 , x0 ) may contribute to ∆Xi ,Yj (t, x|t0 , x0 ). As figure 1 illustrates, provided a portion of the dipole’s world line belongs to the intersection of the backward lightcone of (t0 , x0 ) and the forward lightcone of (t, x) we have t0r ≥ ta . In this case ∆Xi ,Yj (t, x|t0 , x0 ) receives a non-zero advanced contribution coming from the first line in Eq. (57).

6 Photo-detection theory The theory of photo-detection is fundamental within quantum optics and was initiated some time ago [10, 14, 15]. In this section we briefly verify that when source and detector are taken as identical two-level dipoles, advanced wave correlations do not contribute to photo-detection amplitudes within the rotating-wave and Markov approximations. The source is located at the origin 0 and the detector is located at a position x with x ω0−1 . The interaction Hamiltonian for the model is given by ˆ s · DT (0) − d ˆ d · DT (x) V = −d

(59)

where the subscripts s and d stand for source and detector respectively. The state of the composite system in which the source is in the excited state, the detector is in the ground state, and the field contains no photons is denoted |ψi. The rate of excitation of the detector is found using the equation of motion for hσdz (t)iψ to be Z p˙d (t) = − di dj

t

0

dt0 e−iω0 (t−t )

0

× hDT,i (t, x)DT,j (t0 , x)σdz (t0 )iψ + c.c. .

(60)

Since actual photo-detection events depend on photoionisation of an electron within the continuum of final detector states, it is argued in [16] that the rate of stimulated emission from an excited state should be assumed to be much smaller than the rate of absorption from the ground state. This assumption allows us to make the Accepted in

Quantum 2017-12-15, click title to verify

(61)

The total transverse displacement field at the position of the detector includes the detectors own reaction field, and is given by DT (t, x) = E0 (t, x) + Es (t, x) + DT,d (t, x).

(62)

Using again the perturbative approximation σd± (t) ≈ for the detector operators, and the rotatingwave approximation, all terms in Eq. (61) that involve the detector’s own reaction field are seen to vanish under the assumption that the detector is in the ground state. Moreover, while the pure vacuum correlation (+) (−) function hE0,i (t, x)E0,j (t0 , x)iψ can be shown to contribute to the detector level-shifts, it does not contribute within the rotating-wave and Markov approximations to the detectors excitation [16]. We therefore obtain Z t 0 p˙d (t) = di dj dt0 CEi Ej (t, x|t0 , x)e−iω0 (t−t ) + c.c. σd± e±iω0 t

0

(63) The advanced contributions to this expression are found using Eq. (58) with x = x0 . Since t0 takes values from 0 upto t > 0 the first term on the right-hand-side of Eq. (58) does not contribute. It now suffices to note that the second term on the right-hand-side of Eq. (58) pos0 sesses an oscillatory dependence of the form e−iω0 (t−t ) , so that when substituted into Eq. (63) this term gives only rapidly oscillating contributions to p˙d (t), which are negligible using the rotating-wave approximation. Within the rotating wave approximation we therefore have that Z t 0 p˙d (t) = di dj dt0 GEi Ej (t, x|t0 , x)e−iω0 (t−t ) + c.c., 0

(64) which is essentially the same result originally derived by Glauber [10]. A few concluding remarks are in order regarding the final result Eq. (64). We note that the replacement of CEi Ej by GEi Ej in the last step is a result of the rotating-wave approximation. Therefore, advanced wave components of CEi Ej may in principle contribute to the detector excitation, but such contributions are no larger than conventional counter-rotating contributions involving, for example, the detector’s own reaction field. The rotating-wave approximation is generally considered to be valid for optical frequencies and realistic detector response times [16]. On the other hand, 9

in the more realistic case of multi-level non-identical dipoles, the validity of the rotating-wave approximation is less clear. One could also consider detection on much shorter time scales. However, it is not entirely clear how regimes which require moving beyond the rotating-wave approximation should be modelled. In such regimes a retention of the conventional ontology consisting of atomic excitations and photons defined in terms of a bare energy H0 , which is not even approximately conserved, seems physically dubious. Little more can be said presently about the role of advanced waves within such regimes. We remark that the techniques found in [2, 9, 33], which allow one to eliminate counter-rotating terms from the dipole-field interaction Hamiltonian may offer a possible recourse.

7 Radiative forces In this section we show that certain measurement statistics involving the radiative force acting on an elementary quantum point charge q due to the field of a quantum dipole, depend on the two-time correlation functions found in section 5. The point charge has position and momentum operators r and p. The Lorentz force F(t) = m¨r(t) experienced by the charge due to an electromagnetic field is m¨r(t) =qE(t, r(t)) q + [˙r(t) × B(t, r(t)) − B(t, r(t)) × r˙ (t)]. 2 (65) The momentum π = m˙r is found by integrating Eq. (65). As an example of a statistical quantity that depends on two-time correlations in the external fields we consider the rate of change of the dispersion of the momentum, i.e., the momentum diffusion, which is d ∆π(t) = (hπ(0) · F(t)i + c.c.) − 2hπ(0)i · hF(t)i dt Z t + dt0 [hF(t) · F(t0 )i + c.c.] − 2hF(t)i · hF(t0 )i . 0

(66) Substitution of the Lorentz force given by Eq. (65) into Eq. (66), leads to a dependence of the momentum diffusion on integrated two-time correlation functions of the electromagnetic field. Unlike in photo-detection theory, these correlation functions are not accompanied by oscillatory factors due to the detector itself. As a result the advanced components of the relevant correlation functions are not generally negligible. We now derive an explicit expression for the momentum diffusion based on a series of simplifying assumptions and approximations. We begin by noting that a Accepted in

Quantum 2017-12-15, click title to verify

free charge q can be described using the Coulomb-gauge Hamiltonian H0 =

1 [p − qAT (r)]2 + Vself (r) + Hf 2m

(67)

where AT denotes the transverse vector potential, Vself denotes the Coulomb self-energy of the charge, and the Coulomb gauge free field Hamiltonian is defined in ˙ T as terms of the transverse electric field ET = −A Z 1 d3 x [ET (x)2 + B(x)2 ]. (68) Hf = 2 The total electric field is found by adding to ET the longitudinal field EL = −∇φcoul where φcoul is the Coulomb potential of the point charge. The vacuum and source components of the transverse electric field ET,0 and ET,s , which are defined by ET = ET,0 + ET,s are such that changes in ET,0 are generated by the Hamilto˙ T,0 (t, x) = −i[ET,0 (t, x), Hf (0)]. nian Hf alone, i.e., E The source component ET,s (t, x) includes two separate acausal components (appendix 9.4). The first is equal to −EL (t, x), which is independent of the definition of Hf and which is cancelled out within the source-field Es (t, x). The second instantaneous component depends on the source at the initial time t = 0, and arises because in the Coulomb gauge Hf is defined in terms of ET rather than E itself. Due to this component the electric source-field Es is not causal within the Coulomb gauge. A complete derivation of the Coulomb gauge electric source-field is given in appendix 9.4. The occurrence of a static precursor within the electric source-field also arises in the Coulomb gauge treatment of a dipole system [23]. On the other hand, in the multipolar gauge the free field Hamiltonian Hf is defined in terms of the transverse displacement field DT , which possesses a fully retarded source component and is equal to the total electric field away from the dipole [6, 23]. Using superscripts C and m to refer to the Coulomb and multipolar gauges respectively, we have m C m that although EC 0 6= E0 and Es 6= Es it can be shown C m m C that E0 − E0 = Es − Es , so that the total field E is indeed gauge-invariant [23]. Only the multipolar gauge however, correctly identifies the physical vacuum and source components E0 and Es , such that Es is fully causal and is the quantised version of the classical field of an oscillating dipole as given in Eq. (9). To correctly identify the vacuum and source components of the electromagnetic fields when considering the free charge q, we seek a description analogous to the multipolar description of the dipole. To this end we consider a heavy charge with relatively slow dynamics characterised by the time scale tq Γ−1 . We assume a state of the charge such that initially the spread in position is small compared with an optical 10

√ wavelength ∆r ω0−1 , and such that the spread in Dopplerp shifts is small compared to the dipole’s decay rate; ω0 ∆p/m Γ. The initial state of the charge is therefore well-localised in position and momentum. Since the coupling of the charge to the field is weak and t tq we have that AT (t, r(t)) ≈ AT (t, r(0) + p(0)t/m). Furthermore since the atom is initially well-localised in position and momentum we have AT (t, r(0) + p(0)t/m) ≈ AT (t, r0 + p0 t/m) where r0 = hr(0)i and p0 = hp(0)i. Finally we restrict our attention to a charge for which p0 = 0, which implies that AT (t, r(t)) ≈ AT (t, r0 ). In summary, for times t tq the initial spatial localisation of the charge is maintained, which means that hr(t)i − r0 remains small compared with ω0−1 . For times t tq the kinetic energy of the charge is therefore Hq (t) ≈

1 [p(t) − qAT (t, r0 )]2 . 2m

(69)

In this Hamiltonian the charge’s initial wave-packet centre r0 acts like the reference centre of a dipole, which allows us to construct a description of the charge analogous to the multipolar description of a dipole . Specifically, we are now able to transform out the vector potential dependence of Hq (t) via a unitary transformation e−iqr(t)·AT (t,r0 ) . We transform the Coulomb gauge Hamiltonian for the whole system consisting of the charge, dipole and field using the unitary e−iqr(t)·AT (t,r0 ) e−id(t)·AT (t,0) . We then make the two-level and rotating-wave approximations for the dipole and in the usual way, we absorb all dipole and charge self-energy terms into the definitions of their respective bare energies. This yields the final Hamiltonian H = H0 + V where H0 p2 + = ω0 σ σ + 2m + −

Z

(−)

X † 1 d k ω aλ (k)aλ (k) + , 2 3

λ

(+)

V = d · ΠT (0)σ − + σ + d · ΠT (0) + qr · ΠT (r0 ) (70) with (+)

ΠT (t, x) Z Xr = −i d3 k λ

ω eλ (k)aλ (t, k)eik·x . 2(2π)3

(71)

The field canonical momentum ΠT coincides with minus the total electric field away from the dipole and the charge; ΠT (x) = −E(x), x 6= 0, r0 . In the description that uses the Hamiltonian in Eq. (70) the electric source-field is a properly causal field, and moreover, the charge’s canonical momentum p coincides with the mechanical momentum m˙r. Accepted in

Quantum 2017-12-15, click title to verify

Since we can replace π(0) by p0 = 0 in the first term in Eq. (66) this term is seen to vanish. To obtain the momentum diffusion correct to second order in q we ˙ require the force p(t) correct to first order in q, which is found using the Hamiltonian (70) to be ˙ p(t) = qEe (t, r0 ),

Ee = E0 + Es ,

(72)

where E0 and Es are given in Eqs. (5) and (12) respectively. By only considering terms upto order q we find that the charge does not influence the dipole source operators on which the external field in Eq. (72) depends. We may therefore also use Eqs. (36) and (37). Assuming the initial state |e, 0i for the dipole and the radiation field, the momentum diffusion is to O(q 2 ) found using Eq. (72) to be Z t d dt0 hEe (t, r0 ) · Ee (t0 , r0 )i0,e + c.c. ∆p(t) = q 2 dt 0 (73) The correlation function in this expression includes an infinite pure-vacuum contribution which is independent of the dipole. A similar contribution was found in the context of photo-detection theory in section 6 where it was found to be equivalent to the detectors own radiation-reaction. This can be understood as a fluctuation-dissipation relation [16]. Neglecting the pure-vacuum contribution in Eq. (73) while retaining all components that depend on the source dipole we obtain Z t d ∆p(t) = q 2 dt0 CEi Ei (t, r0 |t0 , r0 ) + c.c. (74) dt 0 where CEi Ei is defined in Eq. (48). Using Eq. (50) we write Eq. (74) as the sum of two contributions; d∆p/dt = d∆ps /dt + d∆pvac−s /dt. Assuming the initial state |e, 0i and restricting our attention to the radiation zone x ω0−1 , the first contribution is found using Eq. (52) to be Z t d ∆ps (t) = q 2 dt0 GEi Ei (t, r0 |t0 , r0 ) + c.c. dt 0 2θ(tr ) = N Γ × e− 2 tr (Γ cos[ω0 tr ] + 2ω0 sin[ω0 tr ]) − Γe−Γtr (75) where 1 q 2 E rad (r0 )2 = N (Γ/2)2 + ω02

(76)

with E rad the radiative component of E given in Eq. (13). What remains is the advanced-wave contribution arising due to vacuum source-field correlations. 11

10

10

N [∆pvac−s (t) − ∆pvac−s ] N [∆ps (t) − ∆ps ]

0

2r0 c

N [∆pvac−s (t) − ∆pvac−s ] N [∆ps (t) − ∆ps ]

6

0

Γt

Figure 2: The contributions N [∆ps (t) − ∆ps ] and N [∆pvac−s (t)−∆pvac−s ] are plotted. The decay rate Γ = 108 is chosen in the optical regime. To clearly illustrate the oscillatory character the transition frequency has been chosen as ω0 = 10Γ. The time r0 /c where c denotes the speed of light in the vacuum has been chosen as 1/(3Γ).

This contribution is found using Eq. (58) to be Z t d 2 ∆pvac−s (t) = q dt0 h∆Ei Ei (t, r0 |t0 , r0 )i + c.c. dt 0 1 = θ(tr − r0 ) Γ(2e−Γtr + 1) N

2r0 c

4 Γt

Figure 3: The contributions N [∆ps (t) − ∆ps ] and N [∆pvac−s (t) − ∆pvac−s ] are plotted as in figure 2, but with the larger transition frequency ω0 = 100Γ. These parameters are well within the Markovian regime ω0 Γ. Comparing figures 2 and 3 we see that for fixed Γ increasing ω0 /Γ results in more rapid oscillations, but the same overall behavior with increasing time.

Γ

− e− 2 t Γ(eΓr0 + 2) cos[ω0 (t − 2r0 )] Γr0 − 2ω0 (e − 2) sin[ω0 (t − 2r0 )] .

(77)

Comparing Eqs. (75) and (77) we see that the contribution from GEi Ei is non-zero for times t ≥ r0 while the contribution from h∆Ei Ei i is only non-zero for times t ≥ 2r0 . This allows the two contributions to be clearly distinguished. For t < 2r0 the world line of the dipole does not pass through the intersection of the backward lightcone of (t, r0 ) and the forward lightcone of any point on the charge’s world line with time coordinate in the interval [0, t]. In this case there are no advanced contributions to CEi Ei (t, r0 |t0 , r0 ). The change in momentum dispersion is Z ∆p(t) − ∆p = 0

t

dt0

d ∆p(t0 ), dt0

(78)

which can be partitioned in the same way as the momentum diffusion as the sum of two terms ∆ps (t)−∆ps and ∆pvac−s (t)−∆pvac−s . These two contributions are shown separately in figures 2 and 3. Both contributions have an oscillatory character, which is increasingly suppressed as time increases. Accepted in

Quantum 2017-12-15, click title to verify

It is instructive to compare this scenario with that of a free particle whose initial Gaussian wave-packet spreads in time, but whose momentum dispersion is conserved. We see that since the particle is not free the momentum dispersion is oscillatory and is not conserved. However, if we retain only the contribution ∆ps (t) − ∆ps that comes from GEi Ei then in the long time limit when the dipole has completely decayed the charge’s momentum dispersion reaches a steady value, effectively behaving like that of a free particle possessing non-zero average momentum. Of course, even if initially saturated the uncertainty principle for a free charge is highly unsaturated for times on the order of tq due to the spreading of the initial wave-packet. If the contribution ∆pvac−s (t)−∆pvac−s coming from ∆Ei Ei is also included we see that the behaviour of the momentum dispersion is quite different. In contrast to the case of a free particle the momentum dispersion increases linearly in time in the long-time limit, despite the decay of the initially excited dipole. The change in the dispersion of the position, which quantifies the spreading of the wave-packet is found via integration of Eq. (72). Assuming the initial dipole-field state |e, 0i, 12

case of a free charge the product ∆r(t)∆p(t) increases with time such that the uncertainty principle becomes increasingly unsaturated.

13

N [∆p(t) − ∆p]

8 Conclusions

2r0 c

0

4 Γt

Figure 4: The complete change in momentum dispersion N [∆p(t) − ∆p] is plotted with all parameters chosen as in figure 3. Due to the advanced-wave contribution coming from vacuum-source correlations the momentum dispersion increases with time. For large t a good fit for the behaviour of N [∆p(t) − ∆p] is N [∆p(t) − ∆p] = k + Γt where k is a constant.

it is given by q2 t2 ∆p + ∆r(t) − ∆r = 2m2 m2 Z t Z t Z t1 Z t2 × dt1 dt2 dt3 dt4 CEi Ei (t3 , r0 |t4 , r0 ) 0

0

0

0

+ c.c.

(79)

where we have assumed that there are no initial correlations between the position and momentum of the charge, and where we have again neglected the pure vacuum contribution. As well as the free component t2 ∆p/m2 that results in the usual spreading of the wave packet, there is also a dipole-dependent contribution in Eq. (79), which oscillates in time. In the radiation zone r0 ω0−1 the free component dominates, because the remaining contribution decreases as r0−2 . As in the

Accepted in

Quantum 2017-12-15, click title to verify

In this paper we have focused on quadratic functionals of the electromagnetic field associated with a single stationary dipole. We have extended previous perturbative results, which demonstrate that vacuum-sourcefield correlations provide significant contributions to the radiation intensity of the dipole. We have derived general non-perturbative expressions for arbitrary quadratic field functionals using standard optical approximations. Due to correlations between the vacuumfield and the source-fields of the dipole, contributions coming from the advanced green’s function for the waveoperator are generally non-vanishing within unequaltime field correlation functions. This lies in marked contrast to the derivation of the source-fields themselves. The contribution of vacuum-source correlations to photo-detection amplitudes was shown to be insignificant. However, by developing a description of a free charge q analogous to the multipolar description of a dipole, it was shown that vacuum-source correlations yield significant contributions to statistical predictions involving the force experienced by a free charge in the field of a dipole. In classical electrodynamics, while advanced solutions to the wave equation can be used as a theoretical tool, advanced waves are not usually thought to posses any basis in physical reality. Our results indicate that advanced waves associated with the quantum vacuum do exist, and that it should in principle be possible to verify this using an experiment. These results offer yet another signature of the quantum nature of the vacuum and open up interesting prospects for further investigation of advanced-wave like correlations in various branches of quantum optics. Acknowledgement. This work was supported by the UK engineering and physical sciences research council grant number EP/N008154/1. I thank Dr. A. Nazir for useful discussions relating to this work.

13

9 Appendix 9.1 Derivation of source-fields The equation of motion for the photon annihilation operator is found using the Hamiltonian (3) and once formally integrated reads Z t r 0 ω X eλ (k) · dnm σnm (t0 ) ≡ aλ,0 (t, k) + aλ,s (t, k). (80) dt0 e−iω(t−t ) aλ (t, k) = aλ (k)e−iωt + 3 2(2π) 0 nm where σnm (t) = eiHt |ni hm| e−iHt . To go further we express Eq. (80) in terms of interaction picture dipole operators σ ˜nm (t) = σnm (t)e−iωnm t and perform a rotating-wave approximation, which neglects terms oscillating rapidly at frequencies ω + ωnm , n > m. Substituting the resulting expression into Eq. (6), and evaluating the angular integral and polarisation summation yields Z t Z 0 i X ∞ (+) 3 j dt0 e−i(ω+ωnm )(t−t ) σ ˜nm (t0 )eiωnm t (81) dω ω τ (ωx)d Erad,s,i (t, x) = ij nm 2 4π n

sin(ωx) cos(ωx) sin(ωx) + (δij − 3ˆ xi x ˆj ) − ωx (ωx)2 (ωx)3

(82)

The near, intermediate and far zone components of expression (81) are evaluated separately. Using the identity sin x = (eix − e−ix )/2i the far zone component can be written Z t Z ∞ 0 1 X (+) 0 0 ˆ Erad,s (t, x) = [d − x (ˆ x · d )] dt σ ˜ (t ) dω ω 2 ei(ω+ωnm )t e−iωtr − e−iωta (83) nm nm nm 2 8π x n

0

= 2πf (ωmn )[δ(t0 − tr ) ± δ(t0 − ta )], where f is a suitably behaved function, to obtain the radiation source-field 1 X 2 (+) (−) (+) ˆ (ˆ ω [dnm − x x · dnm )]σnm (tr ). Erad,s (t, x) = Erad,s (t, x) + Erad,s (t, x), Erad,s (t, x) = 4πx n

(84)

(85)

The derivations of the near and intermediate-zone components of the electric source-field, and the magnetic sourcefield do not involve any essentially different steps to those above. The final result including both full source-fields is given by Eq. (10).

9.2 Contribution of vacuum source-field correlations to radiated power 9.2.1

Perturbative calculation of vacuum-source-field correlations and associated radiated power

The equation of motion for the operator σnm (t) can be integrated to give Z t o Xn 0 0 σ ˜nm (t) = σnm + i dt0 [dmp · DT (t0 , 0)]˜ σnp (t0 )eiωmp t − [dpn · DT (t0 , 0)]˜ σpm (t0 )eiωpn t . 0

(86)

p

where Z DT (t, 0) = i

d3 k

Xr λ

Accepted in

Quantum 2017-12-15, click title to verify

h i ω † e (k) a (t, k) − a (t, k) . λ λ λ 2(2π)3

(87)

14

The first order component of Eq. (86) is Z t o Xn 0 0 (1) σ ˜nm (t) = i dt0 [dmp · DT,0 (t0 , 0)]σnp eiωmp t − [dpn · DT,0 (t0 , 0)]σpm eiωpn t 0

(88)

p

where Z DT,0 (t, 0) = i

d3 k

Xr λ

i h ω † −iωt iωt . e (k) a (k)e − a (k)e λ λ λ 2(2π)3

(89)

Using Eqs. (85) and (88) the radiation source-field correct to second order can now be completely expressed in terms of operators at t = 0; 1 X 2 (2) (1) ˆ (ˆ Erad,s (t, x) = ω [dnm − x x · dnm )]˜ σnm (tr )eiωnm tr . (90) 4πx n,m nm Meanwhile the vacuum electric field is Z E0 (t, x) = i

3

d k

Xr λ

ω eλ (k)aλ (k)e−iωt+ik·x + H.c.. 2(2π)3

(91)

Using Eqs. (90) and (91) we obtain hE0 (t, x) ·

(2) Erad,s (t, x)i0;e

=

Z

i (8π 2 )2 x

d3 k ω eik·x

X

2 ˆ (ˆ ωen [den · eλ (k)] eλ (k) · [dne − x x · dne )]

n

Z ×

tr

h i 0 0 0 dt0 eiω(t −t) eiωen (t −tr ) − e−iωen (t −tr ) .

(92)

0

We now evaluate the angular integral and sum over polarisations, and we retain only terms which vary as x−2 , to give Z tr X 1 (2) 0 2 2 iωen (t0 −tr ) −iωen (t0 −tr ) ˆ dt ω [d − x (ˆ x · d )] e − e hE0 (t, x) · Erad,s (t, x)i0;e = − en en en 4(2π)3 x2 0 n Z ∞ 0 0 dω ω 2 eiω(t −tr ) − eiω(t −ta ) . (93) 0

Next we perform a rotating wave-approximation as in the derivation of Eq. (83), which yields Z tr Z ∞ 1 (2) 0 2 iω(t0 −tr ) iω(t0 −ta ) dt dω ω e − e hE0 (t, x) · Erad,s (t, x)i0;e = − 4(2π)3 x2 0 0 " # X X 2 2 iωen (t0 −tr ) 2 2 −iωen (t0 −tr ) ˆ (ˆ ˆ (ˆ × ωen [den − x x · den )] e − ωen [den − x x · den )] e . n>e

n

(94) Since the integrand is now dominated by positive resonant frequencies ω = ωen , n < e and ω = −ωen , n > e, we can perform the Markov approximation, Eq. (84), to obtain 1 1 X (2) 4 ˆ (ˆ sgn(ωen )ωen [den − x x · den )]2 (95) hE0 (t, x) · Erad,s (t, x)i0;e = 2 (4πx)2 n where we have used Z 0

t

dt0 δ(t0 − t)f (t0 ) =

1 f (t). 2

From Eq. (95) we obtain the corresponding contribution to the radiated power; Z 1X (2) Pvac−s = dΩ x2 hE0 (t, x) · Erad,s (t, x)i0;e + c.c. = sgn(ωen )ωen Γen . 2 n

(96)

(97)

which is Eq. (34). Accepted in

Quantum 2017-12-15, click title to verify

15

9.2.2

Radiation reaction: derivation of the Heisenberg-Langevin equations

Restricting our attention to a two-level dipole with ground state |gi and excited state |ei, we define the dipole operators σ − |gi he| ,

σ + = |ei hg| ,

σ z = [σ + , σ − ].

(98)

The two-level transition frequency, dipole moment and decay rate are denoted ω0 , d and Γ respectively. From the solution (80) we calculate the dipole’s own reaction field as (+) DT,s (t, 0)

Z

r

ω X eλ (k)aλ,s (t, k) 2(2π)3 λ 3 Z t h i ˆZ ∞ 0 0 0 ω Γ d ˜ + (t0 )eiω0 t + σ ˜ − (t0 )e−iω0 t . dω dt0 eiω(t −t) σ =i 2π d 0 ω0 0 =i

3

d k

(99)

where we have evaluated the angular integral and polarisation summation. We now perform the rotating-wave approximation and then the Markov approximation, Eq. (84), which give (+) DT,s (t, 0)

Z =i

d3 k

r

ˆ Γd ω X eλ (k)aλ,s (t, k) = i σ − (t) 3 2(2π) 2d

(100)

λ

where we have used used Eq. (96). Substituting Eq. (100) into the Heisenberg equation for each of σ + and σ z yields the following equations of motion Γ Γ 0+ σ˙ + (t) = iω0 − σ + (t) + σ − (t) + σvac (t) (101) 2 2 0z σ˙ z (t) = −Γ[σ z (t) + 1] + σvac (t)

(102)

where 0+ σvac (t) =

Z

d3 k

Xr λ

0z σvac (t) =

Z

s d3 k

X λ

h i ω † z iωt z −iωt d · e (k) a (k)σ (t)e − σ (t)a (k)e λ λ λ 2(2π)2

(103)

h i 2ω † + − iωt + − −iωt d · e (k) a (k)[σ (t) − σ (t)]e − [σ (t) − σ (t)]a (k)e . λ λ λ (2π)2

(104)

Making the rotating-wave approximation in Eqs. (101) and (102) yields optical Heisenberg-Langevin type equations (with no external driving). Formally integrating these equations then yields σ + (t) = e(iω0 −Γ/2)t σ + +

Z

t

0

+ dt0 e(iω0 −Γ/2)(t−t ) σvac (t0 ) 0 Z t 0 z −Γt z z σ (t) + 1 = e [σ + 1] + dt0 e−Γ(t−t ) σvac (t0 ).

(105) (106)

0

where Z

Xr

ω (−) d · eλ (k)a†λ (k)σ z (t)eiωt = i[d · DT,0 (t, 0)]σ z (t) (107) 2(2π)2 λ s Z h i X 2ω (−) † z 3 − iωt + −iωt σvac (t) = − d k d · e (k) a (k)σ (t)e + σ (t)a (k)e = −2i[d · DT,0 (t, 0)]σ − (t) + H.c. λ λ λ (2π)2 + σvac (t) =

d3 k

λ

(108) Making the rotating-wave approximation yields the solutions in Eqs (36) and (37). Accepted in

Quantum 2017-12-15, click title to verify

16

9.2.3

Non-perturbative calculation of vacuum-source-field correlations and associated radiated power

We now use the solutions (36) and (37) to calculate the radiated power. We note that although Markov and rotating wave approximations have been employed we have not used perturbation theory. The pure source-field contribution to the intensity is immediately found using the solution (85) restricted to two dipole levels; 2 ω02 ˆ (d · x ˆ )] hσ + (tr )σ − (tr ) + σ − (tr )σ + (tr )i0;e = [d − x 4πx 2 2 ω0 ˆ (d · x ˆ )] . = [d − x 4πx

hErad,s (t, x) · Erad,s (t, x)i0,e

(109)

The corresponding contribution to the power is then easily obtained as Ps = ω0 Γ/2. The vacuum source-field correlation function is found using Eqs. (91) and (85) to be (+)

hE0 (t, x) · Erad,s (t, x)i0,e = hE0 (t, x) · Erad,s (t, x)i0,e =

ω02 (+) ˆ (d · x ˆ )] · hE0 (t, x)[σ + (tr ) + σ − (tr )]i0;e . [d − x 4πx (110)

From Eqs. (36) and (38) it follows that (+)

hE0 (t, x)[σ + (tr ) + σ − (tr )]i0;e Z

r

ω X eλ (k)haλ (k)σ + (tr )i0;e e−iωt+ik·x 2(2π)3 λ Z Z tr X ω 0 0 = i d3 k eλ (k)[eλ (k) · d]eik·x dt0 hσ z (t0 )i0;e e−iω(t−t ) e(iω0 −Γ/2)(tr −t ) 3 2(2π) 0 (+)

= hE0 (t, x)σ + (tr )i0;e = i

d3 k

(111)

λ

Substituting this result into Eq. (110) we obtain (+)

(+)

(−)

hE0 (t, x) · Erad,s (t, x)i0,e = hE0 (t, x) · Erad,s (t, x)i0,e Z Z tr X ω 0 0 0 iω 2 ik·x ˆ ˆ d3 k (e (k) · d)(e (k) · [d − x (d · x )])e dt0 (2e−Γt − 1)e(iω0 −Γ/2)(tr −t ) e−iω(t−t ) (112) = 0 λ λ 4πx 2(2π)3 0 λ

where we have used hσ z (t)i0;e = 2e−Γt −1, which follows from Eq. (37). We see from Eq. (112) that the rotating-wave (+) (−) and Markov approximations do not eliminate the anti-normally ordered contribution hE0 (t, x) · Erad,s (t, x)i0,e , which is a slowly varying resonant contribution. It is nonetheless not included in the correlation function hE(−) (t, x)· E(+) (t, x)i0,e . Evaluating the angular integral and polarisation summation now yields (+)

(−)

hE0 (t, x) · Erad,s (t, x)i0,e ω 2 0 ˆ (d · x ˆ )] = [d − x 4πx

Z

tr

0

dt

2e

−Γt0

0

Z ∞ h i 0 0 (iω0 −Γ/2)(tr −t0 ) 1 −1 e dω ω 2 e−iω(tr −t ) − e−iω(ta −t ) . (113) 2π 0

We now make use of the Markov approximation, Eq. (84), which gives (+) hE0 (t, x)

·

(−) Erad,s (t, x)i0,e

=

2 ω02 1 −Γtr ˆ (d · x ˆ )] [d − x e − 4πx 2

(114)

where we have used Eq. (96). From Eq. (114) we obtain the corresponding contribution to the radiated power; Z 1 (+) (−) Pvac−s = dΩ x2 hE0 (t, x) · Erad,s (t, x)i0;e + c.c. = ω0 Γ e−Γtr − , (115) 2 which is Eq. (45). Accepted in

Quantum 2017-12-15, click title to verify

17

9.3 Unequal-time commutators of electromagnetic fields and advanced-wave correlations 9.3.1

Unequal-time commutators

We assume that dipole operators σ ± (t) and σ z (t) commute with the field operators aλ (t, k) and a†λ (t, k). We demonstrate here how various commutation relations between the fields E± , B± at arbitrary space-time points can (+) (−) be calculated. As our example we take the electric field commutator [Ei (t, x), Ej (t0 , x0 )], which can be written as the sum of four terms; (+)

[Ei

(−)

(t, x), Ej

(t0 , x0 )]

(+)

(−)

(+)

(−)

(+)

(−)

(+)

(−)

= [E0,i (t, x), E0,j (t0 , x0 )] + [E0,i (t, x), Es,j (t0 , x0 )] + [Es,i (t, x), E0,j (t0 , x0 )] + [Es,i (t, x), Es,j (t0 , x0 )]. (116) The pure vacuum term is easily found using Eq. (5) to be (+) (−) [E0,i (t, x), E0,j (t0 , x0 )]

Z

d3 k

=

X λ

0 0 ω e−iω(t−t )+ik·(x−x ) . 2(2π)3

(117)

The pure source term is also easily found with the help of Eq. (12) to be (+)

(−)

[Es,i (t, x), Es,j (t0 , x0 )] = Ei∗ (x)Ej (x0 )θ(tr )θ(t0r )[σ − (tr ), σ + (t0r )]

(118)

where tr = t − x and t0r = t0 − x0 . Simplification of the remaining two terms on the right-hand-side of Eq. (116), which involve both the vacuum and source-fields requires more work. We begin by noting that by rearranging the rotating-wave approximated integrated Heisenberg equation for the operator aλ (t, k), we can express aλ (k) as aλ (k) = eiωt aλ (t, k) −

Z

t

0

dt0 eiωt

r

0

ω [eλ (k) · d]σ − (t0 ). 2(2π)3

(119)

(−)

Noting that Es,j (t0 , x0 ) is a function of σ + (t0r ) we use Eq. (119) to express the vacuum field as a function of aλ (t0r , k), (+)

(−)

which allows us to write the term [E0,i (t, x), Es,j (t0 , x0 )] as (+) (−) [E0,i (t, x), Es,j (t0 , x0 )]

= −iEj (x

0

)θ(t0r )

Z

3

d k

X λ

ω eλ,i (k)[eλ (k) · d]eik·x 2(2π)3

Z

t0r

ds e−iω(t−s) [σ − (s), σ + (t0r )]

0

(120) where we have used the fact that equal-time dipole and field operators commute. We now evaluate the polarisation summation and integration over solid angle to obtain (+) (−) [E0,i (t, x), Es,j (t0 , x0 )]

i = − 2 Ej (x0 )dk θ(t0r ) 4π

Z

∞ 3

Z

dω ω τik (ωx) 0

t0r

ds e−iω(t−s) [σ − (s), σ + (t0r )]

(121)

0

where τik is defined in Eq. (82). Using Eq. (82) yields (+)

(−)

[E0,i (t, x), Es,j (t0 , x0 )] 2 Z ∞ Z t0r 1 ω 0 0 − + 0 ds [σ (s), σ (tr )] = − 2 Ej (x )dk θ(tr ) dω (δik − x ˆi x ˆk ) eiω(s−tr ) − eiω(s−ta ) 8π x 0 0 iω 1 iω(s−tr ) iω(s−ta ) iω(s−tr ) iω(s−ta ) + 2 (δik − 3ˆ xi x ˆk ) e +e − 3 (δik − 3ˆ xi x ˆk ) e −e , x x

Accepted in

Quantum 2017-12-15, click title to verify

(122)

18

and using the Markov approximation we obtain (+)

(−)

[E0,i (t, x), Es,j (t0 , x0 )] ≈− + =− −

2 Z t0r 1 ω 0 0 − + 0 Ej (x )dk θ(tr ) ˆi x ˆk ) [δ(s − tr ) − δ(s − ta )] ds [σ (s), σ (tr )] 0 (δik − x 4π x 0 1 iω0 (δik − 3ˆ xi x ˆk ) [δ(s − tr ) + δ(s − ta )] − 3 (δik − 3ˆ xi x ˆk ) [δ(s − tr ) − δ(s − ta )] x2 x 2 1 1 ω0 iω0 0 0 0 − + 0 − 3 (δik − 3ˆ xi x ˆk ) Ej (x )dk θ(tr )θ(tr )θ(tr − tr )[σ (tr ), σ (tr )] (δik − x ˆi x ˆk ) + 4π x x2 x 1 1 ω02 iω0 0 0 − + 0 + 3 (δik − 3ˆ xi x ˆk ) Ej (x )dk θ(ta )θ(tr − ta )[σ (ta ), σ (tr )] − (δik − x ˆi x ˆk ) + 4π x x2 x

= − Ei∗ (x)Ej (x0 )θ(t0r )θ(tr )θ(t0r − tr )[σ − (tr ), σ + (t0r )] + Ei (x)Ej (x0 )θ(t0r )θ(ta )θ(t0r − ta )[σ − (ta ), σ + (t0r )]. (123) (+)

(−)

This result coincides with the result given in Eq. (55) with X = E = Y. The remaining term [Es,i (t, x), E0,j (t0 , x0 )] on the right-hand-side of Eq. (116) can be calculated in the same way as the term calculated above and is given by the right-hand-side of Eq. (56) with X = E = Y. Thus, having calculated all source-dependent terms on the right-hand-side of Eq. (116) we have using Eqs. (118) and (123) that (+)

∆Ei Ej (t, x|t0 , x0 ) = [Ei

(−)

(t, x), Ej

(+)

(−)

(t0 , x0 )] − [E0,i (t, x), E0,j (t0 , x0 )]

= Ei (x)Ej (x0 )θ(t0r )θ(ta )θ(t0r − ta )[σ − (ta ), σ + (t0r )] + Ei∗ (x)Ej∗ (x0 )θ(tr )θ(t0a )θ(tr − t0a )[σ − (tr ), σ + (t0a )],

(124)

which coincides with Eq. (57) in the case X = E = Y. Calculations of the same type as above can be used to find the (+) (−) remaining unequal time commutators [Xi (t, x), Yj (t0 , x0 )], X, Y = E, B, which eventually yields Eqs. (55), (56) (+)

and (57). By setting t = t0 in Eq. (57) the right-hand-side vanishes, which proves that [Xi (+)

(−)

(t, x), Yj

(t, x0 )] =

(−)

[X0,i (t, x), Y0,j (t, x0 )] within the Markov and rotating-wave approximations. Thus, the approximate theory is formally consistent as claimed at the end of section 3.2. 9.3.2

Advanced-wave correlations

Supposing that the initial state of the dipole and field is |0, ei the average h∆Xi Yj (t, x|t0 , x0 )ie,0 depends, according to Eq. (57), on the averages h[σ − (ta ), σ + (t0r )]ie,0 and h[σ − (tr ), σ + (t0a )]ie,0 . These averages can be found using Eqs. (36) and (37). The averages hσ + (t0r )σ − (ta )ie,0 and hσ + (t0a )σ − (tr )ie,0 are immediately found to be 0

hσ + (t0r )σ − (ta )ie,0 = e(iω0 −Γ/2)tr e(−iω0 −Γ/2)ta ,

0

hσ + (t0a )σ − (tr )ie,0 = e(iω0 −Γ/2)ta e(−iω0 −Γ/2)tr .

(125)

The average hσ − (ta )σ + (t0r )ie,0 is hσ − (ta )σ + (t0r )ie,0 = di dj

ta

Z

t0r

Z ds

0

0

0

0

(+)

(−)

ds0 e(−iω0 −Γ/2)(ta −s) e(iω0 −Γ/2)(tr −s ) hσ z (s)DT,0,i (s, 0)DT,0,j (s0 , 0)σ z (s0 )ie,0 . (126)

(+)

(−)

Using the commutator [DT,0,i (s, 0), DT,0,j (s0 , 0)] Eq. (126) can be written as the sum of two terms; hσ − (ta )σ + (t0r )ie,0 Z ta Z t0r Z X 0 = ds ds d3 k 0

0

Z + di dj

Z ds

0

Accepted in

λ

ta

0

t0r

0 0 0 ω [eλ (k) · d]2 e(−iω0 −Γ/2)(ta −s) e(iω0 −Γ/2)(tr −s ) e−iω(s−s ) hσ z (s)σ z (s0 )ie,0 2(2π)3

0

0

(−)

(+)

ds0 e(−iω0 −Γ/2)(ta −s) e(iω0 −Γ/2)(tr −s ) hσ z (s)DT,0,j (s0 , 0)DT,0,i (s, 0)σ z (s0 )ie,0 .

Quantum 2017-12-15, click title to verify

(127)

19

Performing the angular integration and polarisation summation and employing the Markov approximation in the first term on the right-hand-side yields Z

ta

Z

t0r

ds

Γ 0

0

0

ds0 δ(s − s0 )e(−iω0 −Γ/2)(ta −s) e(iω0 −Γ/2)(tr −s ) hσ z (s)σ z (s0 )ie,0 .

(128)

0

Since within ∆Xi Yj (t, x|t0 , x0 ) this term is multiplied by θ(t0r − ta ) we can assume that t0r ≥ ta in the above. Performing the s integral then yields Z ta 0 0 0 ds0 eΓs = e(−iω0 −Γ/2)ta e(iω0 −Γ/2)tr (eΓta − 1). (129) Γe(−iω0 −Γ/2)ta e(iω0 −Γ/2)tr 0

To show that the second term on the right-hand-side of Eq. (127) vanishes one can show that for s 6= s0 it is (−) (+) always possible to change the order of σ z (s)DT,0,j (s0 , 0) or of DT,0,i (s, 0)σ z (s0 ) within the expectation value. The integral in Eq. (127), which involves this expectation value must therefore vanish. First consider the commutator (−) [σ z (s), DT,0,j (s0 , 0)]. Using Eq. (119) we can express this as (−)

[σ z (s), DT,0,j (s0 , 0)] = i

Z

d3 k

X λ

ω eλ,j (k)[eλ (k) · d] 2(2π)3

Z

s

0

00

ds00 [σ z (s), σ + (s00 )]eiω(s −s

)

(130)

0

Writing this expression in terms of the interaction picture operator σ ˜ + (s00 ) allows us to use the Markov approximation, which gives (−)

[σ z (s), DT,0,j (s0 , 0)] = i

Z s dˆj dˆj Γ ds00 δ(s00 − s0 )[σ z (s), σ + (s00 )] = i Γ[σ z (s), σ + (s0 )]θ(s − s0 )θ(s0 ). d d 0

(131)

(+)

Similarly the commutator [DT,0,i (s, 0), σ z (s0 )] is found to be i (−)

dˆi Γ[σ + (s), σ z (s0 )]θ(s0 − s)θ(s). d

(+)

(−)

(132)

(+)

Since hσ z (s)DT,0,j (s0 , 0)DT,0,i (s, 0)σ z (s0 )ie,0 = h[σ z (s), DT,0,j (s0 , 0)][DT,0,i (s, 0), σ z (s0 )]ie,0 using Eqs. (131) and (132) the second term on the right-hand-side of Eq. (127) can be expressed in the form Z

ta

Z ds

0

t0r

ds0 f (s, s0 )θ(s − s0 )θ(s0 − s) = 0.

(133)

0

Thus, the second term on the right-hand-side of Eq. (127) vanishes within the Markov approximation. As a result hσ − (ta )σ + (t0r )ie,0 is given by the right-hand-side of Eq. (129). Using Eq. (125) and (129) we therefore have 0

h[σ − (ta ), σ + (t0r )]ie,0 = e(−iω0 −Γ/2)ta e(iω0 −Γ/2)tr (eΓta − 2)

(134)

in agreement with Eq. (58). The remaining commutator h[σ − (tr ), σ + (t0a )]ie,0 can be found in a similar fashion and this yields the remaining part of Eq. (58).

9.4 Static precursors within Coulomb gauge electric source-field of a charge The Hamiltonian (67) yields the following solution for the Coulomb gauge operator aλ (t, k) Z t 0 1 −iωt ˜ 0 , k) aλ (t, k) = e aλ (k) + i dt0 e−iω(t−t ) √ eλ (k) · J(t 2ω 0 ˜ denotes the Fourier transform of the current J; where J Z q ˜ k) = p 1 J(t, d3 x J(t, x)e−ik·x , J(t, x) = [r˙ (t)δ(x − r(t)) + δ(x − r(t))˙r(t)] . 2 (2π)3

Accepted in

Quantum 2017-12-15, click title to verify

(135)

(136)

20

When substituted into the mode-expansion for the transverse electric field Eq. (135) yields the source-field Z ET,s (t, x) = −

3

Z

d k 0

t

1 ˜ 0 , k) − k[ ˆk ˆ · J(t ˜ 0 , k)] eik·x e−iω(t−t0 ) + H.c. J(t dt0 p 2 (2π)3

(137)

Using integration by parts and the continuity equation ρ˙ = −∇ · J where ρ(t, x) = qδ(x − r(t)) is the charge density, this can be written "

Z

i

ˆ ρ(t0 , k) k˜ d k cos[ω(t − t0 )]eik·x ω

#t0 =t

3

ET,s (t, x) = p (2π)3 t0 =0 Z t Z Z t Z 0 d sin[ω(t − t0 )] ik·x 1 ˜ 0 , k) sin[ω(t − t )] eik·x − p ∇ dt0 dt0 ρ˜(t0 , k) d3 k −p d3 k J(t e . dt ω ω (2π)3 0 (2π)3 0

(138)

p The function sin[ω(t−t0 )]/ (2π)3 ω is essentially the Fourier transform of the green’s function for the wave-operator. ˜ L (t, k) = −ik˜ ˆ ρ(t, k)/ω Eq. (138) can be written Using the convolution theorem and Gauss’ law E Z d d3 x0 EL (0, x0 )G(t, x|0, x0 ) ET,s (t, x) = − EL (t, x) + dt Z Z t Z Z t d 3 0 0 0 0 + 0 0 3 0 d x dt J(t , x )G (t, x|t , x ) − ∇ d x dt0 ρ(t0 , x0 )G+ (t, x|t0 , x0 ) (139) − dt 0 0 where G = G+ + G− , and where we have used −p

Z ˆ ρ(0, k) k˜ 1 ˜ L (0, k) cos[ωt]eik·x d3 k E cos[ωt]eik·x = p 3 3 ω (2π) (2π) Z Z d 1 sin[ωt] d 3 ˜ ik·x p = d k EL (0, k) e = d3 x0 EL (0, x0 )G(t, x|0, x0 ). dt (2π)3 ω dt i

Z

d3 k

(140)

By adding the longitudinal field to ET,s (t, x) we obtain the total electric source-field Es (t, x). Along with the well-known retarded source-field given by the second line in Eq. (139), the electric source-field also possesses a term given in Eq. (140), which is dependent on the charge density at the initial time t = 0.

References [1] G. Barton. Elements of Green’s Functions and Propagation: Potentials, Diffusion, And Waves. Oxford University Press, Oxford : New York, new ed edition edition, July 1989. ISBN 978-0-19-851998-0. [2] C. Baxter, M. Babiker, and R. Loudon. Gauge Invariant QED with Arbitrary Mixing of p.a and q.e Interactions. Journal of Modern Optics, 37(4):685–699, 1990. ISSN 0950-0340. DOI: 10.1080/09500349014550761. URL http://www.tandfonline.com/doi/abs/10.1080/09500349014550761. [3] Iwo Bialynicki-Birula. V Photon Wave Function. In E. Wolf, editor, Progress in Optics, volume 36, pages 245–294. Elsevier, January 1996. URL http://www.sciencedirect.com/science/article/pii/ S0079663808703160. DOI: 10.1016/S0079-6638(08)70316-0. [4] Howard J. Carmichael. Statistical Methods in Quantum Optics 1. Springer Berlin Heidelberg, Berlin, Heidelberg, 1999. ISBN 978-3-642-08133-0 978-3-662-03875-8. URL http://link.springer.com/10.1007/ 978-3-662-03875-8. DOI: 10.1007/978-3-662-03875-8. [5] H. B. G. Casimir and D. Polder. The Influence of Retardation on the London-van der Waals Forces. Physical Review, 73(4):360–372, February 1948. DOI: 10.1103/PhysRev.73.360. URL https://link.aps.org/doi/10. 1103/PhysRev.73.360. [6] Claude Cohen-Tannoudji, Jacques Dupont-Roc, and Gilbert Grynberg. Photons and atoms: introduction to quantum electrodynamics. Wiley VCH, March 1997. ISBN 0-471-18433-0. Accepted in

Quantum 2017-12-15, click title to verify

21

[7] J. Dalibard, J. Dupont-Roc, and C. Cohen-Tannoudji. Vacuum fluctuations and radiation reaction : identification of their respective contributions. Journal de Physique, 43(11):1617–1638, November 1982. ISSN 0302-0738. DOI: 10.1051/jphys:0198200430110161700. URL http://dx.doi.org/10.1051/jphys: 0198200430110161700. [8] P. C. W. Davies. Scalar production in Schwarzschild and Rindler metrics. Journal of Physics A: Mathematical and General, 8(4):609, 1975. ISSN 0305-4470. DOI: 10.1088/0305-4470/8/4/022. URL http://stacks.iop. org/0305-4470/8/i=4/a=022. [9] P. D. Drummond. Unifying the p.a and q.e interactions in photodetector theory. Physical Review A, 35 (10):4253–4262, May 1987. DOI: 10.1103/PhysRevA.35.4253. URL http://link.aps.org/doi/10.1103/ PhysRevA.35.4253. [10] Roy J. Glauber. The Quantum Theory of Optical Coherence. Physical Review, 130(6):2529–2539, June 1963. DOI: 10.1103/PhysRev.130.2529. URL http://link.aps.org/doi/10.1103/PhysRev.130.2529. [11] S. W. Hawking. Particle creation by black holes. Communications in Mathematical Physics, 43(3):199–220, August 1975. ISSN 0010-3616, 1432-0916. DOI: 10.1007/BF02345020. URL https://link.springer.com/ article/10.1007/BF02345020. [12] John David Jackson. Classical electrodynamics. Wiley, 3 edition, August 1998. ISBN 0-471-30932-X. [13] R. L. Jaffe. Casimir effect and the quantum vacuum. Physical Review D, 72(2):021301, July 2005. DOI: 10.1103/PhysRevD.72.021301. URL https://link.aps.org/doi/10.1103/PhysRevD.72.021301. [14] L. Mandel. Fluctuations of Photon Beams and their Correlations. Proceedings of the Physical Society, 72(6): 1037, 1958. ISSN 0370-1328. DOI: 10.1088/0370-1328/72/6/312. URL http://stacks.iop.org/0370-1328/ 72/i=6/a=312. [15] L. Mandel, E. C. G. Sudarshan, and E. Wolf. Theory of photoelectric detection of light fluctuations. Proceedings of the Physical Society, 84(3):435, 1964. ISSN 0370-1328. DOI: 10.1088/0370-1328/84/3/313. URL http: //stacks.iop.org/0370-1328/84/i=3/a=313. [16] P. W. Milonni, D. F. V. James, and H. Fearn. Photodetection and causality in quantum optics. Physical Review A, 52(2):1525–1537, August 1995. DOI: 10.1103/PhysRevA.52.1525. URL http://link.aps.org/ doi/10.1103/PhysRevA.52.1525. [17] Peter W Milonni. The quantum vacuum: an introduction to quantum electrodynamics. Academic Press, Boston, 1994. ISBN 0-12-498080-5 978-0-12-498080-8. [18] Peter W. Milonni, Jay R. Ackerhalt, and Wallace Arden Smith. Interpretation of Radiative Corrections in Spontaneous Emission. Physical Review Letters, 31(15):958–960, October 1973. DOI: 10.1103/PhysRevLett.31.958. URL http://link.aps.org/doi/10.1103/PhysRevLett.31.958. [19] A.S. Moskalenko, C. Riek, D.V. Seletskiy, G. Burkard, and A. Leitenstorfer. Paraxial Theory of Direct Electrooptic Sampling of the Quantum Vacuum. Physical Review Letters, 115(26):263601, December 2015. DOI: 10.1103/PhysRevLett.115.263601. URL https://link.aps.org/doi/10.1103/PhysRevLett.115.263601. [20] S. Jay Olson and Timothy C. Ralph. Extraction of timelike entanglement from the quantum vacuum. Physical Review A, 85(1):012306, January 2012. DOI: 10.1103/PhysRevA.85.012306. URL https://link.aps.org/ doi/10.1103/PhysRevA.85.012306. [21] E. A. Power and T. Thirunamachandran. Quantum electrodynamics with nonrelativistic sources. IV. Poynting vector, energy densities, and other quadratic operators of the electromagnetic field. Physical Review A, 45(1): 54–63, January 1992. DOI: 10.1103/PhysRevA.45.54. URL http://link.aps.org/doi/10.1103/PhysRevA. 45.54. [22] E. A. Power and T. Thirunamachandran. Quantum electrodynamics with nonrelativistic sources. V. Electromagnetic field correlations and intermolecular interactions between molecules in either ground or excited states. Physical Review A, 47(4):2539–2551, April 1993. DOI: 10.1103/PhysRevA.47.2539. URL http://link.aps.org/doi/10.1103/PhysRevA.47.2539. [23] E. A. Power and T. Thirunamachandran. Time dependence of operators in minimal and multipolar nonrelativistic quantum electrodynamics. I. Electromagnetic fields in the neighborhood of an atom. Physical Review A, 60(6):4927–4935, December 1999. DOI: 10.1103/PhysRevA.60.4927. URL http://link.aps.org/doi/10. 1103/PhysRevA.60.4927. [24] Edwin A. Power. Zero-Point Energy and the Lamb Shift. American Journal of Physics, 34(6):516–518, 1966. DOI: 10.1119/1.1973082. URL http://link.aip.org/link/?AJP/34/516/1. [25] Benni Reznik, Alex Retzker, and Jonathan Silman. Violating Bell’s inequalities in vacuum. Physical Review A, Accepted in

Quantum 2017-12-15, click title to verify

22

[26]

[27]

[28]

[29] [30] [31] [32]

[33]

[34] [35]

[36]

71(4):042104, April 2005. DOI: 10.1103/PhysRevA.71.042104. URL https://link.aps.org/doi/10.1103/ PhysRevA.71.042104. C. Riek, D. V. Seletskiy, A. S. Moskalenko, J. F. Schmidt, P. Krauspe, S. Eckart, S. Eggert, G. Burkard, and A. Leitenstorfer. Direct sampling of electric-field vacuum fluctuations. Science, 350(6259):420–423, October 2015. ISSN 0036-8075, 1095-9203. DOI: 10.1126/science.aac9788. URL http://science.sciencemag.org/ content/350/6259/420. C. Riek, P. Sulzer, M. Seeger, A. S. Moskalenko, G. Burkard, D. V. Seletskiy, and A. Leitenstorfer. Subcycle quantum electrodynamics. Nature, 541(7637):376–379, January 2017. ISSN 0028-0836. DOI: 10.1038/nature21024. URL http://www.nature.com/nature/journal/v541/n7637/full/nature21024.html. A. Salam. Molecular quantum electrodynamics in the Heisenberg picture: a field theoretic viewpoint. International Reviews in Physical Chemistry, 27(3):405–448, 2008. ISSN 0144-235X. DOI: 10.1080/01442350802045206. URL http://www.tandfonline.com/doi/abs/10.1080/01442350802045206. Akbar Salam. Molecular Quantum Electrodynamics: Long-Range Intermolecular Interactions. Wiley, 1 edition, November 2009. ISBN 0-470-25930-2. J. E. Sipe. Photon wave functions. Physical Review A, 52(3):1875–1883, September 1995. DOI: 10.1103/PhysRevA.52.1875. URL https://link.aps.org/doi/10.1103/PhysRevA.52.1875. Herbert Spohn. Dynamics of charged particles and their radiation field. Cambridge University Press, June 2007. ISBN 0-521-03707-7. Adam Stokes. Quantum optical dipole radiation fields. European Journal of Physics, 37(3):034001, 2016. ISSN 0143-0807. DOI: 10.1088/0143-0807/37/3/034001. URL http://stacks.iop.org/0143-0807/37/i=3/ a=034001. Adam Stokes, Andreas Kurcz, Tim P. Spiller, and Almut Beige. Extending the validity range of quantum optical master equations. Physical Review A, 85(5):053805, May 2012. DOI: 10.1103/PhysRevA.85.053805. URL http://link.aps.org/doi/10.1103/PhysRevA.85.053805. W. G. Unruh. Notes on black-hole evaporation. Physical Review D, 14(4):870–892, August 1976. DOI: 10.1103/PhysRevD.14.870. URL https://link.aps.org/doi/10.1103/PhysRevD.14.870. Theodore A. Welton. Some Observable Effects of the Quantum-Mechanical Fluctuations of the Electromagnetic Field. Physical Review, 74(9):1157–1167, November 1948. DOI: 10.1103/PhysRev.74.1157. URL http://link. aps.org/doi/10.1103/PhysRev.74.1157. Andrew Zangwill. Modern Electrodynamics. Cambridge University Press, Cambridge, December 2012. ISBN 978-0-521-89697-9.

Accepted in

Quantum 2017-12-15, click title to verify

23