Next Article in Journal
Three-Dimension Calculation for the Scattering Problem for Non-Spherical Potential
Next Article in Special Issue
Electric Octupole-Dependent Contributions to Optical Binding Energy
Previous Article in Journal
Nanocomposite-Supported Polymeric Composites Prepared with Different Deposition Bases: Characterization and Application in X-ray Shielding
Previous Article in Special Issue
Large Angular Momentum States in a Graphene Film
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Scalar QED Model for Polarizable Particles in Thermal Equilibrium or in Hyperbolic Motion in Vacuum

1
College of Optical Sciences and Department of Physics, University of Arizona, Tucson, AZ 85721, USA
2
Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA
*
Author to whom correspondence should be addressed.
Physics 2024, 6(1), 356-367; https://doi.org/10.3390/physics6010023
Submission received: 7 November 2023 / Revised: 5 January 2024 / Accepted: 15 January 2024 / Published: 5 March 2024

Abstract

:
We consider a scalar QED (quantum electrodynamics) model for the frictional force and the momentum fluctuations of a polarizable particle in thermal equilibrium with radiation or in hyperbolic motion in a vacuum. In the former case the loss of particle kinetic energy due to the frictional force is compensated by the increase in kinetic energy associated with the momentum diffusion, resulting in the Planck distribution when it is assumed that the average kinetic energy satisfies the equipartition theorem. For hyperbolic motion in vacuum the frictional force and the momentum diffusion are similarly consistent with an equilibrium with a Planckian distribution at the temperature T = a / 2 π k B c . The quantum fluctuations of the momentum imply that it is only the average acceleration a that is constant when the particle is subject to a constant applied force.

1. Introduction

The question of whether there is radiation from a uniformly accelerated oscillator in vacuum continues to be of interest [1] (and refernces therein). In this paper, we revisit the question within the context of a 2D (2-dimensional) scalar QED (quantum electrodynamics) model [2,3]. Our approach differs from previous study with this model in that we focus on the force on the particle and its momentum fluctuations as it exchanges energy with the field. Since it is well established that the particle will perceive itself to be in a field with a Planck distribution at a temperature T = a / 2 π k B c , where a is the acceleration, is the reduced Planck constant, c is the speed of light, and k B is the Boltzmann constant [4,5], we first consider, within the 2D scalar QED model, a polarizable particle in equilibrium with thermal radiation in the absence of any external force that would accelerate the particle. We follow the original approach of Einstein and Hopf [6] and later Einstein [7] in which it is shown that the frictional force and the momentum diffusion experienced by the particle are consistent with the Planck distribution for the field. The 2D scalar QED model exhibits the basic physics of the problem while permitting some simplification compared to the complete QED theory [8,9]. Essentially the same approach in the case of a uniformly accelerated particle in vacuum yields the expected result: the particle undergoes a frictional force and momentum fluctuations as if it were in thermal equilibrium at the temperature T = a / 2 π k B c . The result for the momentum fluctuations seems particularly interesting in the way that the Bose–Einstein form of the fluctuations emerges. The coupling of the particle to the quantum field results in quantum fluctuations about its classical hyperbolic trajectory. These fluctuations are superposed on an otherwise uniformly accelerated motion. In the following Section we present the basic equations of motion needed for our analysis. In Section 3 we derive, based on the Lorentz transformation of the spectral energy density of the field, the frictional force on a polarizable particle in a field of blackbody radiation, and in Section 4 we derive the momentum fluctuations experienced by such a particle. The expressions obtained for the frictional force and the momentum fluctuations are shown to be consistent with a Planck distribution for the spectral energy density of the field. In Section 5 we extend these considerations to the case of hyperbolic motion in a vacuum. Our results are summarized in Section 6.

2. Hamiltonian and Equations of Motion

We consider a point particle interacting with a (1 + 1)D scalar field, ϕ ( y , t ) , with y the coordinate and t the time, via a harmonically bound charged oscillator, as shown in Figure 1. The Lagrangian of the system is [10]
L = 1 2 M Y ˙ 2 + 1 2 m x ˙ 2 1 2 m ω 0 2 x 2 + d y t ϕ ( y , t ) 2 y ϕ ( y , t ) 2 + e ϕ ( y , t ) x ˙ δ ( y Y ) .
Here, M denotes the mass of the particle, m and e are the mass and charge of the oscillator, respectively, ω 0 is the resonant frequency of the oscillator, a = / a , δ ( · ) is the Dirac delta function, and the dot over a letter denotes the time derivative. We denote the center-of-mass position of the particle by Y ( t ) and the charged oscillator coordinate by x ( t ) . The charged oscillator interacts with the field and its y coordinate is constrained to be at Y ( t ) . The bilinear coupling term is motivated by the minimal coupling form of interaction, considering that one can draw a correspondence between ϕ and the vector potential, as discussed in similar models used in microscopic descriptions of mirrors and optomechanics [10,11,12]. Defining the canonical conjugate momenta associated with the various degrees of freedom as P y = M Y ˙ , p x = m x ˙ + e ϕ ( Y , t ) and Π ( y , t ) = ϕ ˙ ( y , t ) , we obtain the corresponding Hamiltonian in the minimal coupling form as
H = P y 2 / 2 M + 1 2 m ω 0 2 x 2 + 1 2 m p x e ϕ ( Y , t ) 2 + k ω k a k a k ,
with
ϕ ( y , t ) = k C k a k ( t ) e i k y + a k ( t ) e i k y , C k = 2 π ω k L 1 / 2 .
a k and a k are the bosonic annihilation and creation operators, respectively, following the canonical commutation relations [ a k , a k ] = δ k k , with δ i j the Kronecker delta, and the quantization volume is of length L. (We omit the speed of light c in Equations (1)–(3)) From the Heisenberg equations of motion that follow from the Hamiltonian (2) we obtain the following equations for the operators x , ϕ 0 , and (see Appendix A)
x ¨ + 2 β x ˙ + ω 0 2 x = e m ϕ 0 ( Y , t ) t ,
ϕ 0 ( y , t ) = k C k a k e i ( ω k t k y ) + a k e i ( ω k t k y ) ,
P ˙ y = e 2 x ˙ ϕ 0 ( Y , t ) y + ϕ 0 ( Y , t ) y x ˙ .
ϕ 0 ( y , t ) is the source-free part of ϕ ( y , t ) ; the source (“radiation reaction”) part of ϕ ( y , t ) results in the damping force 2 β m x ˙ in Equation (4). β is found without approximation to be π e 2 / m , consistent with the commutation relation [ x ( t ) , m x ˙ ( t ) ] = i maintained by the formal solution of Equation (4); see, for instance, Ref. [13].
Equations (4) and (5), together with the equation for the radiation field ϕ ( y , t ) , are the basis for the analyses by Raine, Sciama and Grove [2] and Ford and O’Connell [3] and their conclusion that an oscillator in hyperbolic motion does not radiate. We focus instead on Equation (6). In other words, we focus not on the radiation field but on the force acting on the particle. In the case of hyperbolic motion, this force results in fluctuations superposed on the otherwise uniform acceleration.
For the linear response of x ( t ) to ϕ 0 ( Y , t ) / t we define the polarizability
α ( ω ) = e 2 / m ω 0 2 ω 2 2 i β ω
such that the steady-state expression for e x ( t ) for the particle at rest is
e x ( t ) = i k C k ω k [ α ( ω k ) a k e i ω k t α ( ω k ) a k e i ω k t ] ,
where we have assumed without loss of generality that Y = 0 specifies the particle coordinate.

3. Frictional Force

A polarizable particle moving in a blackbody field experiences a drag force, as shown nonrelativistically by Einstein and Hopf [6] and Einstein [7]. To obtain a relativistic expression for the drag force we can proceed as in Ref. [8] or Ref. [9], for instance. Here we take a different, more heuristic approach based on the Lorentz transformation of the Planck spectrum in the 2D scalar QED model. Our starting point is the expression for the free-space energy density of the 2D massless scalar field:
E = 1 2 ϕ t 2 + ϕ y 2 ,
which is proportional to ω 2 for a mode of frequency ω . In a frame S′ moving with velocity v with respect to the frame in which the energy density is E, the energy density E ω 2 , with
ω = γ ω ( 1 + v ) E = γ 2 ( 1 + v ) 2 E , γ = ( 1 v 2 ) 1 / 2 .
We define the spectral energy density ρ ( ω ) such that ρ ( ω ) d ω is the energy per unit length (or “volume”) in the frequency interval [ ω , ω + d ω ] . In the moving frame ρ ( ω , k ^ ) ρ ( ω , k ^ ) ,
ρ ( ω , k ^ ) d ω = E E ρ ( ω , k ^ ) d ω = 1 γ 2 ( 1 + v ) 2 ρ ( ω , k ^ ) d ω
for a plane wave of frequency ω propagating in the direction of the unit vector k ^ , and therefore
ρ ( ω , k ^ ) = 1 γ 2 ( 1 + v ) 2 ρ ( ω , k ^ ) d ω d ω = 1 γ ( 1 + v ) ρ ( ω , k ^ ) .
For plane waves propagating in the direction parallel to the velocity v ,
ρ ( ω ) = ρ [ γ ω ( 1 + v ) ] γ ( 1 + v ) .
The energy eigenvalues for the free field are ( ω k = | k | )
k ω k n ( ω k ) + 1 2 = L 2 π d k ω n ( ω ) + 1 2 L 0 d ω ρ ( ω ) ,
where the spectral energy density
ρ ( ω ) = ω π n ( ω ) + 1 2 ,
with n ( ω ) the average photon number at frequency ω . From Equations (13) and (15),
ω n ( ω ) + 1 2 = ω n γ ω ( 1 + v ) + 1 2 ,
and so the Lorentz-transformed spectral energy density is
ρ ( ω ) = ω π n γ ω ( 1 + v ) + 1 2 .
In the case of thermal equilibrium at temperature T,
n γ ω ( 1 + v ) + 1 2 = 1 2 coth γ ω 2 k B T ( 1 + v )
and
ρ ( ω ) = ω 2 π coth γ ω 2 k B T ( 1 + v ) .
This is the 2D version of the Lorentz-transformed Planck spectrum obtained in earlier work [14,15].
The force on the particle may be expressed as
F = 1 2 0 d ω R ( ω ) ρ ( ω ) = 2 π 0 d ω R ( ω ) ω n γ ω ( 1 + v ) + 1 2 ,
where R ( ω ) is the rate at which the field loses energy (and momentum p = E ) in a given direction ( k > 0 or k < 0 ), this momentum being taken up by the particle. Perhaps the simplest way to obtain an expression for R ( ω ) is to start from the wave equation for ϕ [10]:
2 ϕ y 2 2 ϕ t 2 = 4 π e x ˙ δ ( y Y ) ,
which has the same form as the wave equation for an electric field propagating in one direction in a medium with number density N ( = δ ( y Y ) here) of particles with polarizability α ( ω ) . In that case the energy dissipation rate R ( ω ) = 2 k I ( ω ) 4 π N ω α I ( ω ) in the dilute-medium approximation, where k I and α I ( ω ) are respectively the imaginary parts of the wave number and the polarizability [16].
In the present case we have one particle in our volume L, implying R ( ω ) = 4 π ω α I ( ω ) and
F = 2 0 d ω ω 2 α I ( ω ) n γ ω ( 1 + v ) + 1 2 .
This expresses the average force in terms of the linear response of the particle to a broadband field with Doppler-shifted frequencies ω ( 1 + v ) / ( 1 v ) . Including field modes with k either parallel or anti-parallel to the particle velocity v , we obtain the net force
F = 2 0 d ω ω 2 α I ( ω ) n γ ω ( 1 + v ) n γ ω ( 1 v ) 4 v 0 d ω ω 3 α I ( ω ) n ω = 4 π v 0 d ω ω α I ( ω ) ρ ω ρ ω .
In the full QED theory, in contrast, the nonrelativistic approximation to the friction force is [14,15]
F = 4 π v 0 d ω ω α I ( ω ) ρ ω 3 ρ ω ,
with ρ ( ω ) = ω 3 n ( ω ) / π 2 c 3 . It is noteworthy that Einstein and Hopf [6] and Einstein [7] obtain ρ ( ω ) ω 3 ρ ω in their formula for the friction force from essentially just kinematic effects involving Doppler shifts and aberration. The factor 1/3 that appears in (24) is due to the vector character of the three-dimensional electromagnetic field, whereas our expression (23) assumes a one-dimensional scalar field.

4. Momentum Fluctuations and Thermal Equilibrium

From Equation (6),
P y ( t ) = i K k C K C k ω K 2 k [ α ( ω K ) a K a k g K k ( t ) α ( ω K ) a K a k f K k ( t ) + α ( ω K ) a K a k f K k ( t ) α ( ω K ) a K a k g K k ( t ) ] ,
where
f K k ( t ) = 0 t d τ e i ( ω K ω k ) τ and g K k ( t ) = 0 t d τ e i ( ω K + ω k ) τ .
The calculation of momentum fluctuations Δ P y 2 ( t ) = P y 2 ( t ) , averaged over the state of the field, is simplified if we assume that ϕ 0 and its derivatives may be regarded as Gaussian random processes. This allows us to treat the K and k modes ( K k ) in Equation (25) as uncorrelated. This property was originally employed in a classical context by Einstein and Hopf [6,17]. A simplified proof of this property for Gaussian random processes X and Y with X = Y = X Y = 0 is given in Ref. [9]. We further simplify by dropping terms involving g K k ( t ) , since their nonresonant character results in no contribution to P y 2 ( t ) . Then
P y 2 ( t ) = K k C K 2 C k 2 ω K 4 k 2 | α ( ω K ) | 2 a K a K a k a k + a K a K a k a k | f K k ( t ) | 2 L 2 π 2 d K 2 π ω K L ω K 4 | α ( ω k ) | 2 d k 2 π ω k L k 2 × [ n ( ω K ) + 1 ] n ( ω k ) + n ( ω K ) [ n ( ω k ) + 1 ] sin 2 1 2 ( ω K ω k ) t ( ω K ω k ) 2 / 4 4 2 0 d ω ω 4 | α ( ω ) | 2 0 d ω [ 2 n ( ω ) n ( ω ) + n ( ω ) + n ( ω ) ] 2 π t δ ( ω ω ) = 16 π 2 t 0 d ω ω 4 | α ( ω ) | 2 [ n 2 ( ω ) + n ( ω ) ] .
In the penultimate step we have used the long-time approximation d K sin 2 ( ω K ω k ) t / 2 ( ω K ω k ) / 2 2 d K 2 π t δ ( ω K ω k ) .
In thermal equilibrium the average rate of increase of the kinetic energy must balance the average rate F v at which the particle loses energy because of the frictional force:
d d t P y 2 2 m + F v = 0 ,
or, from (23) and (27), and the energy equipartition assumption, 1 2 m v 2 = 1 2 k B T , for thermal equilibrium at temperature T,
8 π 2 m 0 d ω ω 4 | α ( ω ) | 2 [ n 2 ( ω ) + n ( ω ) ] + 4 k B T m 0 d ω ω 3 α I ( ω ) n ω = 0 .
Using the relation | α ( ω ) | 2 = α I ( ω ) / 2 π ω (Equation (7)), and the fact that α I ( ω ) > 0 for all ω > 0 , we obtain
n ω = k B T ( n 2 + n ) ,
which is satisfied as expected by
n ( ω ) = 1 e ω / k B T 1 .

5. Hyperbolic Motion in Vacuum

We expect the results above to apply to the case of hyperbolic motion in the zero-temperature vacuum if we simply replace k B T by a / 2 π . In particular, we expect the coupling of the particle to the quantum vacuum field to cause quantum fluctuations about its classical trajectory. A calculation showing this seems nevertheless worthwhile. We now present such a calculation for the momentum fluctuations experienced by a polarizable particle in hyperbolic motion in vacuum.
We start with Equations (5), (6) and (8), but now with
t ( τ ) = 1 a sinh a τ and y ( τ ) = 1 a cosh a τ
for hyperbolic motion, and therefore
ω t ( τ ) + k y ( τ ) = ω ( t ( τ ) + y ( τ ) ) = ω a e a τ
for a particle accelerated in the y direction and for modes k = ω > 0 , which we assume for now. We will carry out the calculation with (proper) time durations τ registered by a clock moving with the particle characterized by a polarizability α ( ω ) . Since ϕ 0 ( y , t ) may be expressed as a function of τ , we write it in terms of operators g ( Ω ) and g ( Ω ) as follows:
ϕ 0 ( τ ) = 0 d Ω g ( Ω ) e i Ω τ + g ( Ω ) e i Ω τ .
Then
e x ˙ = 0 d Ω Ω 2 α ( Ω ) g ( Ω ) e i Ω τ + α ( Ω ) g ( Ω ) e i Ω τ ,
ϕ 0 y = i 0 d Ω Ω g ( Ω ) e i Ω τ g ( Ω ) e i Ω τ ,
and
P y ( τ ) = i 2 0 d Ω 1 0 d Ω 2 Ω 1 2 Ω 2 [ α ( Ω 1 ) g ( Ω 1 ) g ( Ω 2 ) G 12 ( τ ) α ( Ω 1 ) g ( Ω 1 ) g ( Ω 2 ) F 12 ( τ ) + α ( Ω 1 ) g ( Ω 1 ) g ( Ω 2 ) F 12 ( τ ) α ( Ω 1 ) g ( Ω 1 ) g ( Ω 2 ) G 12 ( τ ) ] + h . c . ,
where “h.c.” stands for Hermitian conjugate and we have defined
F 12 ( τ ) = 2 e i ( Ω 1 Ω 2 ) τ / 2 sin 1 2 ( Ω 1 Ω 2 ) τ Ω 1 Ω 2 , G 12 ( τ ) = 2 e i ( Ω 1 + Ω 2 ) τ / 2 sin 1 2 ( Ω 1 + Ω 2 ) τ Ω 1 + Ω 2 .
Now we proceed as in the calculation for thermal equilibrium, regarding ϕ 0 and its derivatives as Gaussian random processes. Then, in taking expectation values, we can treat the operators depending on Ω 1 and Ω 2 as uncorrelated.

5.1. Momentum Fluctuations

The momentum fluctuations, Δ P y 2 ( τ ) = P y 2 ( τ ) , of the particle can be obtained as:
P y 2 ( τ ) = 4 Re 0 d Ω 1 0 d Ω 2 0 d Ω 3 0 d Ω 4 Ω 1 2 Ω 3 2 Ω 2 Ω 4 × α ( Ω 1 ) α ( Ω 3 ) g ( Ω 1 ) g ( Ω 3 ) g ( Ω 2 ) g ( Ω 4 ) F 12 ( τ ) F 34 ( τ ) .
We have dropped terms involving G i j and G i j . Because of their nonresonant character they make no contribution to P y 2 .
From Equation (34), again for a particle accelerated along the y direction,
g ( Ω ) = 1 2 π d τ e i Ω τ ϕ 0 ( τ ) = 1 2 π d τ e i Ω τ k C k a k e i ω k a e a τ + a k e i ω k a e a τ = 1 2 π k C k a k ξ ( ω k , Ω ) + a k η ( ω k , Ω ) ,
where we have defined the functions ξ ( ω k , Ω ) and η ( ω k , Ω ) , corresponding to the time-dependent Doppler shift observed by the accelerated observer [18], as
ξ ( ω , Ω ) = d τ e i Ω τ e i ω a e a τ = 1 a Γ i Ω a ω a i Ω / a e π Ω / 2 a ,
η ( ω , Ω ) = d τ e i Ω τ e i ω a e a τ = 1 a Γ i Ω a ω a i Ω / a e π Ω / 2 a ,
where Γ ( x ) is the Gamma function (see Appendix B for details).
We must now evaluate the correlation functions appearing in (39) that relate to the power spectrum of the vacuum fluctuations as seen by the accelerated particle. Consider first
g ( Ω 1 ) g ( Ω 3 ) = 1 2 π 2 k C k 2 ξ ( ω k , Ω 1 ) ξ ( ω k , Ω 3 ) .
We have used the vacuum expectation values a k a k = a k a k = a k a k = 0 and a k a k = δ k k . Now, in the usual mode continuum limit,
k C k 2 ξ ( ω k , Ω 1 ) ξ ( ω k , Ω 3 ) = a 2 Γ i Ω 1 a Γ i Ω 3 a e π ( Ω 1 + Ω 3 ) / 2 a e i ( Ω 1 Ω 3 ) / a × 0 d ω 1 ω ω a i ( Ω 1 Ω 3 ) / a .
The integral may be evaluated with the change of variable x = ln ( ω / a ) to obtain 2 π a δ ( Ω 1 Ω 3 ) , and so
g ( Ω 1 ) g ( Ω 3 ) = 2 π a Γ i Ω 1 a 2 e π Ω 1 / a δ ( Ω 1 Ω 3 ) .
We find similarly that
g ( Ω 2 ) g ( Ω 4 ) = 2 π a Γ i Ω 2 a 2 e π Ω 2 / a δ ( Ω 2 Ω 4 )
and g ( Ω 1 ) g ( Ω 3 ) = g ( Ω 1 ) g ( Ω 3 ) = 0 . After some algebra we obtain
P y 2 ( τ ) = 4 0 d Ω 1 0 d Ω 2 0 d Ω 3 0 d Ω 4 Ω 1 2 Ω 3 2 Ω 2 Ω 4 α ( Ω 1 ) 2   × 2 π a 2 Γ i Ω 1 a 2 Γ i Ω 3 a 2 δ ( Ω 1 Ω 3 ) δ ( Ω 2 Ω 4 )   × e π ( Ω 1 Ω 2 ) / a sin 2 1 2 ( Ω 1 Ω 2 ) τ 1 2 ( Ω 1 Ω 2 ) ] 2 4 2 π a 2 ( 2 π τ ) 0 d Ω 1 Ω 1 6 | α ( Ω 1 ) | 2 Γ i Ω 1 a 4 .
Finally we use the identity
Γ i Ω a 4 = π 2 a 2 / Ω 2 sinh 2 ( π Ω / a ) = π 2 a 2 Ω 2 4 e 2 π Ω / a ( e 2 π Ω / a 1 ) 2
to obtain
P y 2 ( τ ) = 16 π 2 τ 0 d ω ω 4 | α ( ω ) | 2 [ n 2 ( ω ) + n ( ω ) ] ,
where a factor 2 is included to account for modes with k < 0 as well as k > 0 , and have defined
n ( ω ) = ( e ω / k B T 1 ) 1 , T a / 2 π k B .
The result (49) has the same form as Equation (27), as expected. Interestingly, the Bose–Einstein factor n 2 + n in Equation (27) was obtained from the algebra of annihilation and creation operators, whereas the derivation of Equation (49) depended mainly on the analytical properties of ξ ( ω , Ω ) and η ( ω , Ω ) , as well as the vacuum expectation values a k a k = δ k k .

5.2. Frictional Force

Finally we consider briefly the frictional force on the particle’s center-of-mass, starting with
F = d P y ( τ ) d τ .
Using Equation (37) and ignoring the non-resonant terms, we obtain the simplified expression
F = i 2 0 d Ω 1 0 d Ω 2 Ω 1 2 Ω 2 α Ω 1 g Ω 1 g Ω 2 d F 12 ( τ ) d τ + α Ω 1 g Ω 1 g Ω 2 d F 12 ( τ ) d τ + h . c . .
Noting that d F 12 ( τ ) d τ = e i ( Ω 1 Ω 2 ) τ , and using the expectation values of the correlation functions from Equations (45) and (46), yields
F = i 4 π a 0 d Ω Ω 3 Γ i Ω a 2 α Ω e π Ω / a + α Ω e π Ω / a + i 4 π a 0 d Ω Ω 3 Γ i Ω a 2 α Ω e π Ω / a + α Ω e π Ω / a = π a 0 d Ω Ω 3 cosh π Ω / a Γ i Ω a 2 α I Ω
= 0 d Ω Ω 2 α I Ω coth π Ω / a = 2 0 d ω ω 2 α I ( ω ) [ n ( ω ) + 1 2 ] ,
where in the last step we have again defined n ( ω ) = ( e 2 π ω / a 1 ) 1 and have expressed Γ i Ω a 2 as π a / Ω sinh ( π Ω / a ) . This expression has the form F = 1 2 0 R ( ω ) ρ ( ω ) , as in Equation (20) but with ρ ( ω ) appearing instead of the Lorentz-transformed ρ ( ω ) . In our heuristic approach we now argue, as in Section 3, that if the particle has a velocity v the spectral density ρ ( ω ) experienced by it should be ρ ( ω ) , just as in Equation (17), such that
F = 2 0 d ω ω 2 α I ( ω ) n ( γ ω ( 1 + v ) ) + 1 2 .
The net frictional force is then, exactly as in Section 3,
F = 4 v 0 d ω ω 3 α I ( ω ) n ω .
Assuming an equilibrium power balance according to Equation (28) we obtain, from Equations (56) and (49), Equation (30), but now with k B T = a / 2 π if we assume the energy equipartition relation 1 2 m v 2 = a / 2 π for thermal equilibrium at the Davies-Unruh temperature.

6. Summary

The 2D scalar QED model was used previously [2,3] to address the question of whether there is radiation from a uniformly accelerated oscillator in vacuum. We have focused here not on the radiation, or lack thereof, but on the force on a polarizable particle in thermal equilibrium (at some temperature T > 0 ) with the field or in uniform acceleration in a vacuum. In both cases the momentum fluctuations of the particle depend on the photon number variance having the Bose–Einstein form n 2 ( ω ) + n ( ω ) , where n ( ω ) is the mean photon number. In the former case this variance reflects the blackbody photon statistics of real photons undergoing absorption and emission, whereas in the latter case it emerges, along with the Davies-Unruh temperature, as a consequence of the different vacuum field experienced by the particle in hyperbolic motion [19]. The quantum-mechanical momentum fluctuations of the particle imply that only the average acceleration is constant when a constant external force acts on the particle. We also demonstrate that the frictional force and momentum fluctuations obtained for hyperbolic motion in vacuum are related to each other via the power balance (fluctuation-dissipation) relation (28).

Author Contributions

Conceptualization, methodology, writing—original draft and preparation, and writing—review and editing, K.S. and P.W.M.; appendices, K.S. All authors have read and agreed to the published version of the manuscript.

Funding

Kanu Sinha’s research was supported in part by the National Science Foundation under Grant No. NSF PHY-2309341 and by the John Templeton Foundation under Award No. 62422.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

Peter Milonni has fond memories of discussions with Franco Persico many years ago in Trieste and Erice and thanks the editors for the opportunity to contribute to this Special Issue in his memory.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Equations of Motion

The equations of motion for the charge at Y = 0 are given by:
d x d t = 1 2 m p x e ϕ ( 0 , t ) ,
d p x d t = m ω 0 2 x ,
such that
d 2 x d t 2 + ω 0 2 x = e m ϕ ( 0 , t ) t .
Dividing the field into source-free ( ϕ 0 ) and radiation reaction parts ( ϕ R R ) , we get:
d 2 x d t 2 + ω 0 2 x + e m ϕ R R ( 0 , t ) t = e m ϕ 0 ( 0 , t ) t ,
where a friction term arises explicitly from the radiation reaction part.
The Heisenberg equation of motion for the field operator a k is
i d a k d t = a k , k ω k a k a k e m ϕ ( 0 , t ) p x + e 2 m ϕ 2 ( 0 , t )
d a k d t = i ω k a k + i e C k p x e ϕ ( 0 , t ) m
Using Equation (A1) in the above and integrating over time, we obtain:
a k ( t ) = a k ( 0 ) e i ω k t + i e C k 0 t d τ x ˙ ( τ ) e i ω k t τ
Substituting the above in Equation (3), one gets the free and source parts of the field as follows
ϕ ( y , t ) = k C k a k ( 0 ) e i ω k t e i k y + a k ( 0 ) e i ω k t e i k y ϕ 0 ( y , t ) ( Free field ) + k 2 C k 2 e 0 t d τ x ˙ ( τ ) sin t τ Radiation reaction
The radiation reaction part of the field as it appears in the equation of motion for the charge is
ϕ R R ( 0 , t ) t = 2 e k C k 2 ω k 0 t d τ x ˙ τ cos ω k ( t τ ) 2 π e x ˙ t ,
where in the last step we have taken the continuum limit of the summation ( k L 2 π d k ) . Substitution of the above in Equation (A4) yields Equation (4).

Appendix B. Derivation of ξ(ωk, Ω) and η(ωk, Ω)

In this Appendix, we present a derivation of the functions ξ ( ω , Ω ) and η ( ω , Ω ) that are related to the frequency spectrum seen by the accelerated observer [18]. We define the variable z = e a τ , such that, d z = a e a τ d τ , yielding:
ξ ( ω , Ω ) = d τ e i Ω τ e i ω a e a τ = 1 a 0 d z z i Ω / a 1 e z / ( i a / ω )
Substituting z ˜ = z / ( i a / ω ) , we obtain
ξ ( ω , Ω ) = 1 a i a ω i Ω / a 0 d z ˜ z ˜ i Ω / a 1 e z ˜ = 1 a Γ i Ω a ω a i Ω / a e i π / 2 i Ω / a ,
which readily gives Equation (41).
Similarly,
η ( ω , Ω ) = d τ e i Ω τ e i ω a e a τ = 1 a Γ i Ω a ω a i Ω / a e π Ω / 2 a = 1 a Γ i Ω a e i π ω a i Ω / a e π Ω / 2 a ,
which yields Equation (42).

References

  1. O’Connell, R.F. The demise of Unruh radiation. Mod. Phys. Lett. A 2020, 35, 2050329. [Google Scholar] [CrossRef]
  2. Raine, D.J.; Sciama, D.W.; Grove, P.G. Does a uniformly accelerated charge radiate? Proc. R. Soc. Lond. A Math. Phys. Eng. Sci. 1991, 435, 205–215. Available online: https://www.jstor.org/stable/52011 (accessed on 10 January 2024).
  3. Ford, G.W.; O’Connell, R.F. Is there Unruh radiation? Phys. Lett. A 2006, 350, 17–26. [Google Scholar] [CrossRef]
  4. Davies, P.C.W. Scalar production in Schwarzschild and Rindler metrics. J. Phys. A Math. Gen. 1975, 8, 609–616. [Google Scholar] [CrossRef]
  5. Unruh, W.G. Notes on black hole evaporation. Phys. Rev. D 1976, 14, 870–892. [Google Scholar] [CrossRef]
  6. Einstein, A.; Hopf, L. Statistische Untersuchung der Bewegung eines Resonators in einem Strahlungsfeld. Ann. d. Phys. 1910, 338, 1105–1115, English translation: Statistical investigation of a resonator’s motion in a radiation field. In The Collected Papers of Albert Einstein. Volume 3: The Swiss Years: Writings 1909–1911 (English Translation Supplement); Schulmann, R., Kox, A.J., Renn, J., Eds.; Princeton University Press: Princeton, NJ, USA, 1994; pp. 220–230. Available online: https://einsteinpapers.press.princeton.edu/vol3-trans/234 (accessed on 10 January 2024). [CrossRef]
  7. Einstein, A. Zur Quantentheorie der Strahlung. Phys. Z. 1917, 18, 121–128, English translation: On the quantum theory of radiation. In The Collected Papers of Albert Einstein. Volume 6: The Berlin Years: Writings, 1914–1917 (English Translation Supplement); Kox, A.J., Schulmann, R., Eds.; Princeton University Press: Princeton, NJ, USA, 1997; pp. 220–233. Available online: https://einsteinpapers.press.princeton.edu/vol6-trans/232 (accessed on 10 January 2024).
  8. Milton, K.A.; Day, H.; Li, Y.; Guo, X.; Kennedy, G. Self-force on moving electric and magnetic dipoles: Dipole radiation, Vavilov-Cerenkov radiation, friction with a conducting surface, and the Einstein-Hopf effect. Phys. Rev. Res. 2020, 2, 043347. [Google Scholar] [CrossRef]
  9. Sinha, K.; Milonni, P.W. Dipoles in blackbody radiation: Momentum fluctuations, decoherence, and drag force. J. Phys. B At. Mol. Opt. Phys. 2022, 55, 204002. [Google Scholar] [CrossRef]
  10. Wang, Q.; Unruh, W.G. Motion of a mirror under infinitely fluctuating quantum vacuum stress. Phys. Rev. D 2014, 89, 085009. [Google Scholar] [CrossRef]
  11. Sinha, K.; Lin, S.-Y.; Hu, B.L. Mirror-field entanglement in a microscopic model for quantum optomechanics. Phys. Rev. A 2015, 92, 023852. [Google Scholar] [CrossRef]
  12. Sinha, K. A Microscopic Model for Quantum Optomechanics. Ph.D Thesis, University of Maryland, College Park, MD, USA, 2015. [Google Scholar] [CrossRef]
  13. Milonni, P.W. Radiation reaction and the nonrelativistic theory of the electron. Phys. Lett. A 1981, 82, 225–226. [Google Scholar] [CrossRef]
  14. Pauli, W. Theory of Relativity; Pergamon Press Ltd.: New York, NY, USA, 1958; Section 49; Available online: https://archive.org/details/theoryofrelativi00paul/ (accessed on 10 January 2024).
  15. Ford, G.W.; O’Connell, R.F. Lorentz transformation of blackbody radiation. Phys. Rev. E 2013, 88, 044101. [Google Scholar] [CrossRef] [PubMed]
  16. Milonni, P.W.; Eberly, J.H. Laser Physics; John Wiley & Sons, Inc.: New York, NY, USA, 2010; Section 8.2. [Google Scholar] [CrossRef]
  17. Einstein, A.; Hopf, L. Über einen Satz der Wahrscheinlichkeitsrechnung und seine Anwendung in der Strahlungstheorie. Ann. d. Phys. 1910, 338, 1096–1104, English translation: On a theorem of the probability calculus and its application in the theory of radiation. In The Collected Papers of Albert Einstein. Volume 3: The Swiss Years: Writings 1909–1911 (English Translation Supplement); Schulmann, R., Kox, A.J., Renn, J., Eds.; Princeton University Press: Princeton, NJ, USA, 1994; pp. 211–219. Available online: https://einsteinpapers.press.princeton.edu/vol3-trans/225 (accessed on 10 January 2024). [CrossRef]
  18. Alsing, P.M.; Milonni, P.W. Simplified derivation of the Hawking–Unruh temperature for an accelerated observer in vacuum. Am. J. Phys. 2004, 72, 1524–1529. [Google Scholar] [CrossRef]
  19. Cetto, A.M.; de la Peña, L. Real vacuum fluctuations and virtual Unruh radiation. Fortschr. Phys. 2017, 65, 1600039. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of a point particle of mass M and harmonically bound oscillator with a charge e and mass m. The charged oscillator moves along the x-axis with a frequency ω 0 and interacts with a scalar field ϕ ( y , t ) at y = Y , where Y is the center-of-mass position of the point particle. See text for more details.
Figure 1. Schematic representation of a point particle of mass M and harmonically bound oscillator with a charge e and mass m. The charged oscillator moves along the x-axis with a frequency ω 0 and interacts with a scalar field ϕ ( y , t ) at y = Y , where Y is the center-of-mass position of the point particle. See text for more details.
Physics 06 00023 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sinha, K.; Milonni, P.W. Scalar QED Model for Polarizable Particles in Thermal Equilibrium or in Hyperbolic Motion in Vacuum. Physics 2024, 6, 356-367. https://doi.org/10.3390/physics6010023

AMA Style

Sinha K, Milonni PW. Scalar QED Model for Polarizable Particles in Thermal Equilibrium or in Hyperbolic Motion in Vacuum. Physics. 2024; 6(1):356-367. https://doi.org/10.3390/physics6010023

Chicago/Turabian Style

Sinha, Kanu, and Peter W. Milonni. 2024. "Scalar QED Model for Polarizable Particles in Thermal Equilibrium or in Hyperbolic Motion in Vacuum" Physics 6, no. 1: 356-367. https://doi.org/10.3390/physics6010023

APA Style

Sinha, K., & Milonni, P. W. (2024). Scalar QED Model for Polarizable Particles in Thermal Equilibrium or in Hyperbolic Motion in Vacuum. Physics, 6(1), 356-367. https://doi.org/10.3390/physics6010023

Article Metrics

Back to TopTop