Quantum radiation by an Unruh-DeWitt detector in oscillatory motion

Quantum radiated power emitted by an Unruh-DeWitt (UD) detector in linear oscillatory motion in (3+1)D Minkowski space, with the internal harmonic oscillator minimally coupled to a massless scalar field, is obtained non-perturbatively by numerical method. The signal of the Unruh-like effect experienced by the detector is found to be pronounced in quantum radiation in the highly non-equilibrium regime with high averaged acceleration and short oscillatory cycle, and the signal would be greatly suppressed by quantum interference when the averaged proper acceleration is sufficiently low. An observer at a fixed angle would see periods of negative radiated power in each cycle of motion, while the averaged radiated power over a cycle is always positive as guaranteed by the quantum inequalities. Coherent high harmonic generation and down conversion are identified in the detector's quantum radiation. Due to the overwhelming largeness of the vacuum correlators of the free field, the asymptotic reduced state of the harmonics of the radiation field is approximately a direct product of the squeezed thermal states.


Introduction
A uniformly, linearly accelerated point-like detector moving in the Minkowski vacuum will experience thermal fluctuations at a temperature proportional to the proper acceleration of the detector [1]. This is called the Unruh effect and the temperature is called the Unruh temperature. While the derivation in time-dependent perturbation theory is well established, direct experimental evidence for the Unruh effect is still lacking. One closely related observation in laboratories is the electron depolarization in storage rings, namely, the Sokorov-Ternov effect [2], which can be connected to the "circular Unruh effect" [3][4][5][6]. Nevertheless, the centripetal acceleration in the circular Unruh effect is quite different in nature from the original linear, uniform acceleration in the Unruh effect [7]. To get closer to the original conditions in Unruh's derivation, there have been existing proposals to look at the quantum correction by the Unruh effect to the radiation emitted by a linearly accelerated charge or atom [8][9][10], which is called the "Unruh radiation".
Seeking the evidence of the Unruh temperature in quantum radiation is, however, not as straightforward as it appears. To well define a finite temperature in an atom-field state, the atom should be in equilibrium with the field. Unfortunately a uniformly, linearly accelerated Unruh-DeWitt (UD) detector (analogous to an atom) [1,11] derivatively coupled to a massless scalar field in the Minkowski vacuum emits no radiation in equilibrium conditions in (1+1)D Minkowski space [12][13][14][15][16]. In (3+1)D Minkowski space there will be radiation by a uniformly accelerated UD detector in steady state at late times (at a constant radiation rate with respect to the proper time of the detector), but the radiated energy is not converted from the one experienced by the detector in the Unruh effect [17]. The physical reason for these results is that quantum interference between the vacuum fluctuations driving the detector and the radiation emitted by the driven detector is perfectly destructive in equilibrium conditions.
In laboratories, producing an eternal, constant linear acceleration for a charge or an atom is impossible, anyway. In Ref. [8] and similar proposals the charge motion would be driven by an intense laser field, which can make the acceleration linear, but not uniform, thus the radiation corresponding to the Unruh effect may not be totally suppressed by interference. The only concern is that the Unruh temperature is not well defined in these non-equilibrium setups. Fortunately, a time-varying effective temperature whose value is close to the Unruh temperature of the averaged acceleration can be defined for the detectors in oscillatory motion [18]. As we will demonstrate later, in the regime of high acceleration and short oscillating cycle of the motion, the signal of the effective Unruh temperature can be pronounced in the Unruh radiation 1 .
This paper is organized as follows. To get non-perturbative, time-dependent results of the radiation field emitted by a point-like detector in oscillatory motion, in Section 2 we introduce the UD harmonic-oscillator (HO) detector model considered in Refs. [17,18] and then address some technical issues. We determine the radiation in the radiation zone defined in the Minkowski coordinates for a laboratory observer [17,19,20]. Then we present our numerical results of the radiated power in Section 3. We show that when observed at a fixed angle, there will be negative radiated power in some periods during each cycle of the oscillatory motion [21,22]. This indicates that the Unruh radiation observed at the null infinity may correspond to a multi-mode squeezed state of the field [9,10,23,24] or something similar. So we further study the correlations in the radiation field in Section 4, where we identify the nonlinear optical effects such as the coherent highharmonic generation and down-conversion [10]. Then in Section 5 the asymptotic reduced state of the field harmonics is constructed using the correlators of the field in the radiation zone. Our analysis shows that this asymptotic reduced state of the radiation field looks like a direct product of the squeezed thermal states. Finally, a summary of our results is given in Section 6. A few analytic results for the two-point correlators are given in Appendix A, which can help to control the singularities in our numerical calculation.

Renormalized expectation values of stress-energy tensor
Consider an Unruh-DeWitt detector with its internal degree of freedom acting as a harmonic oscillator and minimally coupled to a massless scalar field Φ in (3+1)D Minkowski space, described by the action where g µν = diag(−1, 1, 1, 1), z µ is the worldline of the detector parametrized by its proper time τ , Ω 0 is the bare natural frequency of the internal HO, and λ 0 is the coupling constant of the detector and the field. Here we take c = = G = 1. From (2.1) one can derive the conjugate momentaP = m 0 dQ/dτ andΠ = ∂ tΦ of the detector and the field, respectively. Below we set m 0 = 1 for simplicity. As we discussed in Ref. [17], by virtue of the linearity of this UD detector theory, the field operatorΦ x after the detector-field coupling is switched on will become a linear combination of the mode functions each associated with a creation (b † k ) or annihilation operator (b k ) of the free field mode with wave vector k, or a raising (â † ) or lowering operator (â † ) of the free internal HO of the detector Q. Also due to the linearity, a mode function of the field has the form φ κ φ a x and φ k x are associated withâ andb k , respectively), which is the superposition of the homogeneous solution φ [0] κ x corresponding to vacuum fluctuations of the free field and the inhomogeneous solution φ [1] κ x sourced from the point-like detector. One can group the homogeneous solutions of the mode functions with the associated operators intoΦ [0] x and the inhomogeneous solutions intoΦ [1] x such that the field operator is in the formΦ x =Φ [0] x +Φ [1] x . Suppose the detector-field coupling is switched on at t = t I , when the combined system is initially in the factorized state which is a product of the ground state of the free UD detector |g A and the Minkowski vacuum of the field |0 M . The wavefunction (or density matrix) is thus in a Gaussian form.
In the Heisenberg picture, the correlators of the field amplitude at different spacetime points x and x after t = t I are given by x with respect to the initial state |ψ(0) in (2.2). Then the expectation value of the stress-energy tensor (minimal, ξ = 0) can be written as where T [ij] µν is contributed by G [ij] .
is the Green's function of the free field, it diverges as x → x, so does T [00] µν . Nevertheless, there is no physical effect from this part of the stress-energy tensor in Minkowski space, and so it can be subtracted in the spirit of the normal ordering in obtaining the vacuum energy in the conventional quantum field theory. We thus define the renormalized stress-energy tensor as µν (x) . Doing this is nothing but setting the zero point of vacuum stress-energy.
Suppose a UD detector is oscillating about the spatial origin of the Minkowski coordinates, which is chosen as the laboratory frame. Suppose a set of the radiation-detecting apparatus are located at a large constant radius r at different angles from the spatial origin, namely, at x µ = (x 0 , x 1 , x 2 , x 3 ) = (t, r sin θ cos ϕ, r sin θ sin ϕ, r cos θ) in the Minkowski coordinates. Then the differential radiated power per unit solid angle measured in laboratories can be written as dP dΩ II (t, θ, ϕ) = − lim r→∞ r 2 T tr (x) ren = dP [11] dΩ II + dP [01] dΩ II + dP [10] dΩ II , (2.5) where dΩ II is the element of the solid angle and the [ij] component is defined by In this paper we are considering the cases with the UD detector in linear oscillatory motion in x 3 -direction, namely, z µ (τ ) = (z 0 (τ ), 0, 0, z 3 (τ )), and the radiation will be independent of the azimuth angle ϕ by symmetry.
The whole calculation will be started with the subtracted two-point correlators of the field. However, as many quantities for quantum fields with infinite degrees of freedom, each of G [11] , G [10] and G [01] is still singular in the coincidence limit. One has to control the singularities with the hope that some of them could cancel in the measurable quantities while others could be tamed by introducing physical cutoffs. To identify the problem, let us look into more details of the correlators.

Two-point correlators at late times
At late times, the retarded field has carried the initial information in the detector away to the null infinity, so the behavior of the combined system around the detector is dominated by vacuum fluctuations of the field as well as the detector's response to them. From Ref. [17] and [25], one has where the v-part of the correlator Q(τ )Q(τ ) v is defined by the mode functions associated with the field operators and the initial data in the field state only. Here γ ≡ λ 2 0 /(8π) is the coupling strength, Ω ≡ Ω 2 r − γ 2 is the natural frequency with Ω r renormalized from the bare natural frequency of the detector Ω 0 , K(X) ≡ e −γX sin ΩX is the propagator, τ I is the proper time of the detector at the initial moment, and is the positive-frequency Wightman function of the free massless scalar field in the Minkowski vacuum state, with a proper choice of the integration contour understood [26]. While one could get rid of the divergence of the integrand asτ →τ in (2.7) by choosing the integration contour, the divergence of Q(τ )Q(τ ) in the coincidence limit τ → τ is unavoidable. From (A1) in Ref. [17], denoting x µ = (t, x) and Φ x = Φ x (t), one has G [11] (x, x ) = Φ [1] x (t)Φ [1] x (t ) = with the singular behaviors of (2.7) and other correlators of the detector. Also at late times, (2.12) whose integrands and integrals diverge asτ → τ − (x ) and x → x, respectively. Here we is the retarded distance determined by the local frame of the detector as and η − ≡ τ − (x) − τ I , while R ≡ R(x ) and η − ≡ η − (x ). Note that the retarded distance aX/2 and the retarded proper time for a uniformly accelerated detector in Ref. [17] has been generalized to R and η − for a detector in oscillatory motion here 2 . For an observer at the null infinity the more the 4-velocity of the detector is pointing towards the observer, the smaller R/r is.

Controlling the singularities
From the correlators of the detector QQ , QP , P Q , and P P in (2.10) and thus (2.6), one could extract the Unruh or the effective temperature experienced by the detector [27]. So we call the all-retarded-field part of the differential radiated power dP [11] /dΩ II as the naive Unruh radiation. It diverges when one takes the coincidence limit on the two-point correlators of the detector, as one can see from (2.6), (2.7) and (2.10). When the trajectory of the detector is not as simple as those in uniform motion or uniform acceleration, setting consistent cutoffs in the double integral for the correlators such as (2.7) is not easy. In Refs. [18] and [25] we have dealt with these singularities carefully. We subtract the integral for the two-point correlators of the detector in oscillatory motion by those for a uniformly accelerated detector. The subtracted integral gives a finite result. Then we add the analytic results for the uniformly accelerated detector back, whose singular behavior are well understood and under control once the UV cutoff is introduced.
For the interference terms of the differential radiated power, dP [01] /dΩ II + dP [10] /dΩ II , the situation is similar. As we mentioned, in the integrand of (2.12), (2.14) diverges asτ → τ − (x) and → 0+. When is positive and non-zero, expanding ). To subtract out the divergent −2 and −1 terms, one needs to introduce a reference worldline,z µ (τ ) withz µ (τ − ) = z µ (τ − ),ż µ (τ − ) =ż µ (τ − ), and z µ (τ − ) =z µ (τ − ). For a general worldline z µ at a specific moment τ , the simplest reference worldline for subtraction is again the one for a uniformly accelerated detector, and luckily, we have also obtained the analytic results of the interference terms for the uniformly accelerated detector in closed form in Ref. [17]. Similar to what we did for dP [11] /dΩ II , after we get the finite result for the subtracted interference terms, we add the analytic result back in the final step to get the complete result with the divergences well controlled. 2 In Eq. (A1) in Ref. [17], ∂t∂ r G [11] v (x, x ) only counts the contribution by the v-parts of the correlators of the detector .. v. The expression for the a-part, ∂t∂ r G [11] a (x, x ), has the same form as Eq. (A1) in Ref. [17] except the v-parts of the detector-detector correlators .. v are replaced by the a-parts .. a. Since Φx = 0 in the cases we are considering, the expression for ∂t∂ r G [11] (x, x ) = ∂t∂ r [G [11] a (x, x )+G [11] v (x, x )] is simply the same expression as Eq. (A1) in Ref. [17] except all the v-part of the detector-detector correlators there are replaced by the complete one, namely, .. v → .. = .. a + .. v.
In the cases withz µ (τ − ) = 0, the reference worldline and the analytic result to be added reduce to those for the detector in uniform motion withz µ (τ − ) = z µ (τ − ) anḋ z µ (τ − ) =ż µ (τ − ). Some analytic expressions of the correlators for an UD detector in uniform acceleration and uniform motion are given in Appendix A, for adding back to the subtracted numerical results.
For the reference worldlines either in uniform acceleration or in uniform motion, the UV divergence (Λ 1 in Section A) in T [11] µν will be exactly canceled by the ones in the interference terms T [10] µν + T µν [17]. Thus, after combining the numerical result of the subtracted power and the exact analytic result from the reference worldlines, the final result will be regular and independent of the UV cutoff for the detector.

On-resonance case
When the period of a cycle of oscillatory motion in the proper time of the UD detector is integer times of the natural period of the internal HO (τ p = n × (2π/Ω), n integer), it is possible to get the late-time result with a finite domain of integration to make the numerical calculation more economic. Our experience in calculating the effective temperature in a UD detector in oscillatory motion shows that such kind of the resonance condition is not catastrophic [18]. In these cases, since z 0 (τ − τ p ) = z 0 (τ ) − t p and z 3 (τ − τ p ) = z 3 (τ ) for a detector in oscillatory motion of period τ p in proper time and t p in the coordinate time, (2.16) The integrand can be written in closed form by noting that ∞ n=0 z n (a + n) −s = Φ(z, s, a), which is the Hurwitz-Lerch transcendent [28]. This reduces the domain of the integral from R 2 to a finite square, though the integrand diverges atτ =τ + nτ p for some integer n and have to be treated in the way given in Section 2.2.
Similarly, the late-time interference terms (2.11) and (2.12) in the on-resonance cases can be written as , (2.18) whose domains are reduced to finite intervals, though the integrands also diverge asτ − nτ p → τ − (x ) for some positive integer n and should be properly treated.
When nτ p = 2π/Ω with integer n, the situation is similar to the above case with τ p = n × (2π/Ω). Only minor modifications on the above integrals are needed.

Radiated power
As an example, let us consider a detector moving along the worldline given by Chen and Tajima in Ref. [8], which is the trajectory of a charge at a nodal point of magnetic field in a cavity. The effective temperature of the detector in this worldline has been studied in [18]. Let the coupling is switched on at τ I = 0 (when t = 0). Then η − = τ − . Since τ I does not go to −∞, we are not really at late times here, and in the results of this section we have actually included the a-parts of the correlators [17] in addition to the v-parts discussed in Section 2.1, though the contributions by the a-parts are small in the figures we are going to present. Suppose the observer is located at x µ = (t 0 + r, rr). For r /ω 0 and in the radiation zone, where v µ (τ ) ≡ż µ (τ ) = ( 1 + 4a 2 0 sin 2 ω 0 t, 0, 0, 2a 0 sin ω 0 t) and a µ (τ ) ≡z µ (τ ) are the fourvelocity and four-acceleration of the detector moving along the worldline z µ , respectively, while F (φ, m) is the elliptic integral of the first kind [18]. (3.4) can be quickly derived by partially differentiate the equation σ(x, z(τ − (x))) = 0. Let the period of the oscillatory motion be t p = 2π/ω 0 in the coordinate time and τ p = ω −1 0 F (2π, −4a 2 0 ) in the proper time of the detector. We define the directional proper acceleration √ a µ a µ is the proper acceleration [18]. Then one has a µ = (αv 3 , 0, 0, αv 0 ) for z µ in (3.1), so that |∂ µ R/R| : |∂ µ τ − | ≈ |α| = |a|, µ = 0, 3 around θ = 0 and π in the radiation zone. In the cases with large a 0 ω 0 , |∂ µ R/R| can dominate over |∂ µ τ − | at most of the observing angles. Also the reference worldline to control the singularities for an observer at x would be, The advanced time for this reference worldline or its image reads [17,29] τ We show our numerical results for the renormalized differential radiated power (2.5) emitted by a UD detector moving along (3.1) in Figures 1, 2 and 3.

Negative radiated power
In Figure 1, we demonstrate the time evolution of the full differential radiated power in one cycle at large times at fixed angles. One can see that the time evolution of the full differential radiated power at each fixed angle with the interference terms included has a double-peak structure. Each pulse consists of two main peaks of positive flux, and valleys of negative flux in-between.
In the moving-mirror model in (1+1)D, Fulling and Davies found that a negative radiated power will be observed when the acceleration of the mirror is time varying while the mirror is moving towards the observer [30]. Our results at θ = 0 and π (on the plane of the oscillatory motion of the detector) are consistent with this observation. However, when observed far off the oscillation plane, the negative radiated power in our examples does not occur around the moment when the observed proper acceleration has the most significant change. This is clear in the plots with θ = 0 or π in Figure 1, where the double-peak pulse occurs around the maximum of the observed proper acceleration |α(t − (x))|, while the value of |∂ t α| can be small in the period of negative flux between the peaks. Actually, the doublepeak structure in Figure 1 has been obvious in the naive term dP [11] /dΩ II (green dashed curve). We find this behavior in our example is dominated by the R ,t R ,r /R 4 factor in (2.10) (here a 0 ω 0 = 6.554 in α, see the statement below (3.5).) When the scaled retarded distance R/r for the observer at some fixed angle become very small, the observed radiated power will be amplified and a pulse emerges. However, around the moment that the retarded distance reaches the minimum, one has R ,t = 0 (i.e. a 0 (τ − (x)) = a 3 (τ − (x)) cos θ by (3.5)) Figure 1. The differential radiated power per solid angle dP/dΩ II against t 0 /t p emitted by a detector in the worldline (3.1) and observed at angle θ and time t = t 0 + r in the radiation zone r 2π/ω 0 . Here γ = 0.01, Ω = 2.3, Λ 1 = Λ 0 = 20, r = 10 40 , ω 0 = 3.277, a 0 = 2, and so the averaged proper accelerationā = 10, the period of each cycle t p ≡ 2π/ω 0 = 1.917 in the rest frame, and τ p = 0.838 in the comoving frame. The solid and green dashed curves represent the differential radiated power with and without the interference terms, respectively, while the blue dashed curves represent the naive differential radiated power dP [11] /dΩ II without the Unruh effect, namely, the two-point correlators for the accelerated detector have been replaced by the ones for an inertial detector. The black and red sections of the solid curve represent the positive and negative radiated power. The gray dashed curves represent the scaled directional acceleration of the detector α(t − (x)) at the moment it emitted the observed radiation. and so a valley between two positive main peaks is formed in a pulse. On the other hand, the negative correction from the interference terms become the most negative at some moment a little bit ahead of the valley of the naive term, so the total differential radiated power around the valley become negative (also see the inset of Figure 4 (left)).
At each fixed time, there will always exist negative radiated power around some observing angle, as shown in Figure 2. Nevertheless, the averaged radiated power over a cycle of oscillation at each fixed angle must be positive (Figure 2 (lower-right)), as guaranteed by the quantum inequalities [31].
Looking more closely (also see Figure 3) one can see that, when the averaged proper acceleration [18]ā is sufficiently small (ā = 10 in Figure 1), and the observing angle is in the vicinity of θ = 0 or π, there may be a longer period of negative radiated power between the second positive peak of one pulse and the first positive peak of the successive pulse, while the value of the negative radiated power is very close to zero. In some cases the period of this negative radiated power can be longer than a half of the period of a cycle (e.g. Figure 4). The differential radiated power dP/dΩ II against θ at fixed t 0 for the same detector in Figure 1. The bottom-right plot is the differential radiated power averaged over a cycle of oscillatory motion. The radiated energy is positive for all θ while concentrated around θ = 0 and π, namely, the directions of the linear oscillatory motion.
In our examples the period of the oscillatory motion t p = 2π/ω 0 are small compared to the time scale of the relaxation of the detector 1/γ. Under this condition, if we increase the coupling strength γ with other parameters fixed, the time evolution of the angular distribution of the differential radiated power will be similar to Figures 1 and 2, except that the peaks will be roughly amplified as γ 1 . The ratio of the maximum amplitudes of the negative radiated flux to the positive one does not change significantly as we increase a 0 or 1/ω 0 but keepā fixed.
The negative radiated power in the above result does not imply absorption, or radiation in the opposite direction. It can excite an UD detector at a rate lower than the case in zero energy density [32], and produce no decrease of entropy [33]. Our result reveals another resemblance between the detector theory and the moving-mirror models in quantum field theory in curved spacetime. Actually, some Unruh-DeWitt detector theories in (1+1)D have been used to describe mirrors in a more realistic way than those simply introducing boundary conditions for the fields at the mirror's position [34][35][36][37].

Evidence of Unruh effect in radiation
In Ref. [18] we observed that, at a lower (higher) value of the averaged proper acceleration a, the effective temperature of an UD detector in oscillatory motion tends to be higher (lower) than the naive Unruh temperatureā/(2π) experienced by a uniformly accelerated detector at the proper accelerationā. Indeed, the effective temperature in the example shown in Figures 1 and 2 is about T eff ≈ 1.6754 to 1.6767, which is higher than the naive Unruh temperatureā/(2π) ≈ 1.5916 (ā = 10), while in Figure 3 (right), the effective temperature is about 3.156, which is lower thanā/(2π) ≈ 3.183 (ā = 20). Anyway, in Figure 3. From left to right are the time evolutions of the differential radiated power at θ = π/3 with the averaged accelerations from low to high values (ā = 2.5, 5, and 20, respectively; The result withā = 10 has been shown in the lower-left plot, Figure 1). We compare the full result dP/dΩ II (solid) with the (modified) naive terms dP [11] /dΩ II : the green dashed, red dotted, and blue dashed curves are contributed by the correlators of the detector in oscillatory motion, in uniform acceleration at the proper accelerationā, and at rest, respectively. One can see that dP/dΩ II gets closer to the naive result dP [11] /dΩ II at the effective or Unruh temperature as the averaged accelerationā increases, while the red "tail" of the small negative radiated power after the second pulse of the full result gradually disappears. Also the green dashed curve is much closer to the red dotted curve than to the solid curve, indicating that the deviation of the detector's effective temperature from the Unruh temperature is not important for the radiated power, compared with the correction from the interference terms. Figure 3 one can see that the deviation of the effective temperature from the naive Unruh temperatureā/(2π) due to non-uniform acceleration is negligible in the radiated power, compared with the correction from the interference terms.
A detector at rest still has non-zero correlators Q 2 and P 2 contributed by vacuum fluctuations at zero temperature, such that the naive Unruh radiation dP [11] /dΩ II from (2.10) is positive even at zero averaged acceleration. Since we expect that the radiation by the detector should cease as its averaged proper accelerationā → 0, the negative interference terms should be able to cancel the naive differential radiated power in this case. Indeed, we find the radiated power tends to be suppressed larger by the interference terms when ω 0 orā gets smaller (Figure 3 (left)). Here, for a fixed a 0 , a smaller ω 0 on the one hand gives a smaller averaged proper accelerationā, on the other hand it implies a longer period of oscillatory motion, so that the detector has more time to approach to the equilibrium conditions studied in Ref. [17]. Both suppress the radiated power.
In contrast, as ω 0 orā increases, the importance of the interference terms decreases, and the full result of the differential radiated power get much closer to the naive result dP [11] /dΩ II at the effective or Unruh temperature than the naive result at zero temperature (Figure 3 (middle) to (right)). This suggests that the Unruh-like effect experienced by a UD detector could be observed in the Unruh radiation in highly non-equilibrium conditions, with a very short period of oscillatory motion and a very high averaged proper acceleration.

Correlations in radiation
To obtain the late-time two-point correlators more efficiently, we consider the on-resonance case below. late-time differential radiated power dP/dΩ II observed at θ = 0 (black/red). Here γ = 0.01, Ω = 4.3, Λ 1 = Λ 0 = 20, a 0 = 4,ā ≈ 7.6, and the period of each cycle τ p = 2π/Ω (on-resonance case) in the comoving frame. The black and red parts of the solid curve represent the differential radiated power with positive and negative values, respectively. In the inset, the green-dashed and orange-dotted curves represent the naive term dP [11] /dΩ II and the interference terms d(P [01] + P [10] )/dΩ II , respectively. (Right) |dP/dΩ II | in the frequency domain for the harmonics with ω = nω 0 , n = 1, 2, 3, · · · , after a Fourier transform (black). The maximum of the right hump is located at ω = 313ω 0 , at which frequency the sine wave is plotted as the gray solid curve in the inset of the left plot. Between the two humps the local minimum is located at ω = 115ω 0 . The gray dashed curve represents |d(P [01] + P [10] )/dΩ II | (orange dotted) subtracted by |dP [11] /dΩ II | (green dashed), showing that the interference terms contribute the major part for 131 ≤ n ≤ 597. Note that only the discrete data for n ∈ Z are shown in this plot.

High harmonic generation
The frequency spectrum of the differential radiated power with respect to the observer's clock shows that there are quanta at high harmonics generated in a broad range of frequencies. Indeed, in Figure 4 (right) one can see two main humps in the frequency spectrum of the differential radiated power at θ = 0, the left one has the maximum around ω 0 , and the right one around 313ω 0 . The former is the frequency of the oscillatory motion of the detector, while the latter corresponds to the time-scale of the double-peak pulse around t − ≈ t p /4, when the detector reaches its maximum speed v max = 2a 0 / 1 + 4a 2 0 toward the observer at θ = 0. The major part of the right hump around the maximum is contributed by the interference terms d(P [01] + P [10] )/dΩ II . The two humps are quite broad, which is similar to the high harmonic generation (HHG) in atomic and plasma physics. The HHG here, however, is caused by the relativistic motion [38][39][40], which is not simple-harmonic in spacetime, rather than the bremsstrahlung during the recombination in the three-stage model [41][42][43]. The pulses are compressed in time and concentrated in a narrow angular distribution by the Lorentz boost of the radiated field, similar to the Larmor radiation of a classical charge in relativistic motion in electrodynamics [20].
While the origin of the HHG sounds classical in our linear system, the harmonics generated here are quantum coherent. Later in the scaled two-point correlators of the field r 2 Φ x (t),Φ x (t ) | r→∞ in the radiation zone, we will see that the harmonics are correlated in the frequency domain. Thus the two short pulses of the radiated energy around θ = 0 and π in the upper-left and lower-right plots of Figure 1 in each cycle of the motion is associated with the coherent HHG. In this aspect the highly compressed pulses we found are similar to the attosecond laser.
Interesting enough, a wide range of the right slope of the right hump in Figure 4 (right) behaves like ω 3 /(e ω/(k B T B ) − 1), which is the spectrum of the black-body radiation. The parameter T B ≡ T B (c + v max )/(c − v max ) is blue-shifted from an effective temperature T B . We find T B ≈ 1.08ā ( = k B = c = 1) in the interval 500 ≤ n ≤ 1500 in Figure  4 (right). As we increase the value of a 0 from 4 to 40 and so v max goes deeper into the relativistic regime, we find that T B goes from 1.08ā to 1.005ā, which suggests T B →ā as a 0 → ∞. Further analysis on this observation is ongoing. Note that the value of T B here is not equal to the effective temperature T eff experienced by the UD detector, which is a little less than the naive Unruh temperatureā/2π at large averaged proper accelerationā. Note also that in our fitting for T B we only count the harmonics (ω = nω 0 , n ∈ Z) in Figure 4 (right). The absolute value of the full, continuous spectrum of the radiated power in our case is not as smooth or isotropic as the spectrum of the black-body radiation: it depends on θ, and at each θ it looks like a comb with the peak values located at the harmonics in the frequency domain (cf. Figure 5 (upper left)).

Correlations of pulses
In Figure 5 (upper row) we show the scaled, symmetrized, late-time renormalized correlators r 2 Φ x (t),Φ x (t ) ren | r→∞ , where A, B ≡ (AB + BA) /2, of the field amplitudes observed at different times t and t on the same side (θ = θ = 0, upper-left plot) and opposite sides (θ = 0 and θ = π, upper-right plot) of the axis of the detector motion. Denoting n t ≡ (t 0 /t p ) − (1/4) − N 0 and n t ≡ (t 0 /t p ) − (1/4) − N 0 for t 0 ≡ t − r and t 0 ≡ t − r with some very large r 2π/ω 0 (in the radiation zone) and some large integer N 0 (γt p ) −1 (at late times) 3 , one can see that the correlations of the field amplitudes are amplified and peaked around n t ∈ Z and n t ∈ Z for (θ, θ ) = (0, 0), and around n t ∈ Z and n t − 1/2 ∈ Z for (θ, θ ) = (0, π), respectively. The corresponding t and t are the moments that the pulses reach the observers at θ = 0 (t − r = [n + N 0 + (1/4)]t p , n = 1, 2, 3, · · · ) and θ = 0 (t − r = [n + N 0 + (1/4)]t p ) or π (t − r = [n + N 0 + (3/4)]t p ). Thus the pulses of the radiated power are correlated. The same scaled correlators as functions of the retarded proper times τ − are shown in Figure 6. It is obvious that the peaks in the upper plots of Figure 5 are compressed from the ones in Figure 6 through the nonlinear relation between the observer's time t and the retarded proper time τ − of the detector at fixed r and θ.
Looking more closely to the peaks, one can see that the shape of the peaks on the diagonal axes (t = t ) in Figure 5 (upper-left) is different from the shape of the off-diagonal peaks. Around the diagonal peaks, the value of r 2 Φ x (t),Φ x (t ) ren at θ = θ = 0 is positive for both (t mod t p ) and (t mod t p ) are a little less than τ p /4, and is negative for both (t mod n) and (t mod n) are a little greater than τ p /4 ( Figure 5 (lower-middle)). Around the off-diagonal peaks, if the peak value of the renormalized field correlator is positive (negative), the value will be always positive (negative) in t − t direction, and will become Figure 5. (Upper row) The scaled, symmetrized, late-time renormalized field correlators r 2 Φ x (t),Φ x (t ) ren evaluated at x with θ = 0 and x with θ = 0 (left) and θ = π (right), both in the radiation zone. Here γ = 0.01, Ω = 2.3, a 0 = 4, the period of a cycle of oscillatory motion in the proper time of the detector is exactly two times of the natural period of the internal HO (τ p = 2 × (2π/Ω)), n t ≡ (t 0 /t p ) − (1/4) − N 0 , n t ≡ (t 0 /t p ) − (1/4) − N 0 with some large integer N 0 . (Lower row) Close-ups of the peaks in the upper-left plot around (n t , n t ) = (1, 1) (middle) and (n t , n t ) = (1, 3) (right) for θ = θ = 0. All the peaks in the upper-right plot for (θ, θ ) = (0, π) have similar shape to the one in the lower-right plot. The lower-left plot is the symmetrized two-point correlator Q (τ ),Q(τ ) of the detector, which dominates the signs of the above field correlators. negative (positive) in t+t direction when getting close to but not immediately neighboring to the peak ( Figure 5 (lower-right)). These features are even clearer in Figure 6.
In the example we consider here, the sign of the off-diagonal peaks of the field correlation is dominated by the sign of Q (τ − (t)),Q(τ − (t )) in G [11] in Eq. (2.9). This is evident by comparing Figure 5 (upper-left) and Figure 6 (left) with Figure 5 (lower-left). As |n−n | increases, the absolute peak value of the correlation roughly decreases like e −γ|n−n |τp . Figure 7 shows the correlators of the radiation field transformed from those in Figure 5 to the frequency domain. The real parts of the Fourier transformed correlators of the field at (θ, θ ) = (0, 0) and (θ, θ ) = (0, π) are presented in the upper-left and the upper-right plots, respectively. One can immediately see that in both cases there are strong correlations at the lattice points (ω, ω ) = (mω 0 , m ω 0 ) with m, m = 0, 1, 2, 3, · · · , indicating that different harmonics are coherent, as we claimed earlier.

Correlation of harmonics, and down conversion
The plots in the lower row of Figure 7 show the imaginary parts of the Fourier transformed correlators. In the lower-left plot for θ = θ = 0, one can clearly see the antidiagonal, negative valued "thin walls" in 45 degrees in the unit cells of the lattice of the harmonics. They satisfy the conditions ω + ω = nω 0 , n = 1, 2, 3, · · · . The one with ω + ω = ω 0 indicates that in quantum radiation there exist correlated pairs of field quanta with both ω and ω even lower then the frequency ω 0 of the oscillatory motion driven by some external agent. This may be interpreted as a down conversion process. The similar anti-diagonal thin walls with ω + ω = nω 0 , n = 2, 3, 4, · · · , correspond to the down conversions from the higher harmonics.
As pointed out in Ref. [23], in the leading-order perturbative result for the expectation value of the stress energy tensor, there exist interference terms corresponding to the creation of a pair of the Minkowski particles. In Refs. [9,10], it is found that electrons driven by a strong, periodic electromagnetic field can also convert vacuum fluctuations into entangled photon pairs in the first order perturbation, which can be understood in terms of the Unruh effect. The whole process can be viewed as a down conversion from the photons at the frequency ω 0 of the intense laser to the lower frequencies ω 1 , ω 2 of the entangled photon pair with ω 1 + ω 2 = ω 0 . Since the perturbative results in the time-dependent perturbation theory is valid in transient rather than in equilibrium conditions [45], it is natural to see a similar down conversion in our full result for a detector in linear oscillatory motion, which is non-equilibrium.
In (2.9) and (2.10), the factors such as R ,t R ,r /(RR ) 2 and η −,µ η −,ν /(RR ) are prod- ucts of periodic functions ∼ f (t)g(t ). If a correlation function of t and t can be arranged in the form C(t, t ) = p=1,2,··· f p (t) × g p (t ), where all the functions f p (t) and g p (t ) are periodic at the frequency ω 0 = 2π/t p , then C(t, t ) = p n∈Zfpn e inω 0 t × m∈Zgpm e imω 0 t = p,n,mfpngpm e iω 0 (nt+mt ) will not contribute any nonzero frequency lower than ω 0 with respect to t or t . Since the field correlators have the property G(t, t ) = G(t + t p , t ) or G(t, t + t p ) in general, they are not in the form of p=1,2,··· f p (t)g p (t ), and in (2.9) and (2.10) the correlators of the detector such as Q (η − (t)),Q(η − (t )) could produce down conversion. Similarly in (2.12), the Wightman function D + could do. Thus the cause of down conversion in this model must be quantum.
Performing inverse Fourier transforms on the imaginary parts of the renormalized correlator of the field with ω + ω = ±ω 0 only, we obtain the results in the time domain in the upper row of Figure 8. One can clearly see that the anti-diagonal structures in the lowest unit cells in the lower plots of Figure 7, namely, the down conversions of the fundamental frequency ω 0 , correspond to the periodicity in t + t direction (+45 • on the n t n t -plane) with period 2t p in t + t for the renormalized correlators of the field in Figure  5. The symmetric and antisymmetric behavior of the anti-diagonal structures in the lowest unit cells about the line ω = ω in the frequency domain also propagate to the symmetric and antisymmetric structures about the line t = t in the time domain in the upper plots of Figure 8.
The thin walls of the field correlators in the frequency domain in the lower-left plot of Figure 7 for θ = θ = 0 spread widely in ω − ω direction, thus when the inverse Fourier transform includes all the imaginary part of the renormalized field correlators in the frequency domain, the result will be highly concentrated around the diagonal lattice points in the time domain as shown in Figure 8 (lower-left). The negative correlation around (n t , n t ) = (n + , n + ) with n ∈ Z and > 0 in the vicinity of the diagonal peak centered around (n t , n t ) = (n, n) (see the lower-middle plot of Figure 5) are mainly contributed by these imaginary parts.
The radiation field amplitudes at θ = 0 and θ = 0 can satisfy the phase-matching con- Figure 9. The sampling regions A and A for the asymptotic state, with boundaries colored in purple, in the radiation zone (r 2π/ω 0 ). The blue curve represents the worldline of the detector and the red arrows represent the emitted pulses in the oscillatory plane.
ditions ω + ω = nω 0 and k + k = nω 0x 3 simultaneously. This may explain the significance of the correlation along the anti-diagonal thin walls. In contrast, for the correlators of the field amplitudes at θ = 0 and θ = π, the phase matching conditions cannot be fulfilled, and the anti-diagonal thin walls are significant only in the unit cells of the lowest few harmonics, as shown in Figure 7 (lower right). The contribution by the whole imaginary parts of the Fourier-transformed field correlator at (θ, θ ) = (0, π) (Figure 8 (lower-right)) to the original correlator in the time domain is small and widespread compared with the case of (θ, θ ) = (0, 0) ( Figure 8 (lower-left)).
In the off-resonance cases, the pattern of the late-time field correlators has no exact periodicity in t + t direction. While the details of different pulses would not be exactly the same at late times, each pulse would still be concentrated around t 0 mod t p = t p /4 at θ = 0 or 3t p /4 at θ = π. As down conversions involve low frequency quanta which are not sensitive to the detail of each pulse, we expect down conversions would still appear in the off-resonance cases.

Squeezing in the asymptotic reduced state
It has been argued that the presence of the negative radiated power in Section 3.1 indicates that the Unruh radiation corresponds to a multi-mode squeezed state of the field [9,10]. However, negative energy density can also arise in many field states with even particle numbers, the two-particle Fock state of some field mode could be the simplest [24]. It is thus interesting to examine more details of the field state for the Unruh radiation observed in the radiation zone.

Asymptotic reduced states of the field
A quantum state of the combined detector-field system (2.1) can be described by the density matrixρ[(Q, Φ x ), (Q , Φ x ); x 0 ], or equivalently in the (K, ∆)-representation [44] ρ[K, ∆; where we write (Q, Φ x ) = Σ − (∆/2) and (Q , Φ x ) = Σ + (∆/2) with the boldface letters Σ and ∆ denoting the combined detector and field vectors in the configuration space. Our combined system is linear and started with a Gaussian state (2.2) with Φ µ = Π µ = 0, so the reduced state of the field (with the detector's degree of freedom traced out) simply reads where the indices i, j run over a time-slice in the position or wave-vector space depending on the representation of the field we chose, K j = (K j , ∆ j ) and K j = (K j ) T , and V is the covariance matrix with each element a 2 × 2 matrix of the symmetrized two-point correlators of the field [45]. Then one can extract quantum information such as squeezing and entanglement from the covariance matrix. Nevertheless, no apparatus localized in the radiation zone is possible to get the full information of the field in the whole space, nor do our numerical calculations for the Unruh radiation. What can be measured or numerically calculated in our setting are the asymptotic correlators of the field defined in a finite spatial or temporal region in the radiation zone (Figure 9), with the bound field (counterpart of the velocity field in electrodynamics [19,20]) totally ignored. Inserting these asymptotic correlators of the field into (5.3) and (5.2), what we get is the asymptotic reduced state of the field, which is certainly not the true reduced state of the field.

Mode decomposition
The Minkowski vacuum state of a massless scalar field is the simplest in terms of the field modes: it is a direct product of the ground states for the (complex) HOs Φ k = d 3 xe −ik·x Φ x with the natural frequencies ω = |k| [46,47], thus separable in terms of field modes Φ k . In quantum optics, conventional discussions on the squeezed state are also based on the field modes. To compare with the vacuum state and look into the squeezing in the radiation, therefore, we would represent the asymptotic reduced state of the Unruh radiation in terms of the field modes.
Since the motion of our source (the UD detector) is quite localized in space, if we are looking at the complete field state, we may have to expand the retarded field in terms of the solutions of the massless Klein-Gordon equation in spherical coordinates, which is not easy to deal with in calculations. Fortunately, we are only interested in the asymptotic states here, and so the partial wave analysis applies: The spherical waves with the spherical Hankel function h (1) l (kr) ≈ e ikr /kr locally look like the plane waves in the radiation zone (kr → ∞), namely, e ik(r+∆r) /[k(r + ∆r)] ≈ f (r)e ik∆r for a large constantr. Below we are considering the asymptotic reduced state of the field in the conventional plane-wave field modes, which has the translational symmetry that a true state of the field scattered by a localized source is lacking. As long as the decomposition is linear, the asymptotic reduced state of the field (5.2) will always be Gaussian.

Squeezed thermal harmonic modes
To get a clearer picture and make the numerical calculation more economic, we will further restrict ourselves to the most significant subset of the field degrees of freedom. We have learned in the previous sections that the quantum radiation emitted by the detector is concentrated around θ = 0 and π, or +x 3 and −x 3 directions, respectively. In both directions the field modes of the harmonics with frequency nω 0 , n = 1, 2, 3, · · · are the most interesting. We can construct the covariant matrix of these field modes using the late-time correlators R ±nk 0 ,R ±n k 0 , R, R = Φ, Π, and k 0 = (0, 0, ω 0 ) in the radiation zone. Note that one should not subtract out the correlators of the vacuum state of the free field with ω n ≡ (nω 0 ) 2 + k 2 ⊥ from the elements of the covariant matrix [47], though they are singular as the volume V of the sampling region goes to infinity (V δ nn → (2π) 3 δ 3 (k − k ) as V → ∞).
In our example, we first calculate the sum of the scaled renormalized field correlators at a late time x 0 =t ≡ [N 0 + (3/4)]t p +r with some integer N 0 (γt p ) −1 and somē r 2π/ω 0 in the two finite sampling regions A around θ = 0 and A around θ = π as rr R r ,R r ren ≡ rr where R, R = Φ or Π,r = r − (t p /2),r = r − (t p /2), r, r 2π/ω 0 , rx 3 , r x 3 ∈ A and −rx 3 , −rx 3 ∈ A , to suppress the imaginary part produced by the Fourier transform to be done below. Here the rr factor is to balance the 1/r and 1/r dependence ofR x andR x in the radiation zone. We write r =r +∆r and r =r +∆r with ∆r, ∆r r. Since A and A are finite regions, the wave-vector space is discrete with the minimal line element dk = 2π/L, where L is the radial length scale of the sampling regions. We further assume that dΩ is so small andr is so large thatr x is constant in θ and ϕ directions. Although the Fourier transformed correlator should be integrated over the whole space, we replace the integral by an approximated one to obtain r+∆r could be of the same order. We then arrange the covariance matrix V as Inserting R k ,R k ren in (5.8) and the vacuum correlators (5.4) and (5.5) into (5.9) and then diagonalizing it, we find that each eigen-vector is dominated by a pair of Φ ±k or Π ±k with the same |k| = |n|ω 0 because the correlators (5.4) and (5.5) for the vacuum state of the free field are much greater than the corrections by the radiation field (5.8). Denoting the approximated eigenvectors Φ −nk 0 +r −2 ± n (Π), (5.12) for n = 1, 2, 3, · · · , respectively, where with s R ≡ sign{Re R nk 0 ,R nk 0 ren } taking the value +1 or −1. In fact, the amplitude of a real scalar field in the k-space is complex subject to the relation Φ −k = Φ * k . To fit the standard form of a Gaussian state (5.2) for a collection of quantum mechanical HOs, with different indexes i and j referring to different degrees of freedom, it would be convenient to split Φ k into the real and imaginary parts, namely, , which turn out to be the approximated eigenvectors Φ + k and Φ − k for the above covariance matrix V up to a factor i. They can be thought of as two HOs with index k in a half of the k-space while Φ ± −k (= ±Φ ± k ) are not independent degrees of freedom. According to (5.11) and (5.12), the asymptotic reduced state of the field harmonics in θ = 0 and π directions is approximately a product state of each (nk 0 , −nk 0 ) mode-pair in terms of Φ + nk 0 and Φ − nk 0 . We will see that each asymptotic reduced state of a mode-pair Φ ±nk 0 looks like a two-mode squeezed thermal state (so Φ nk 0 and Φ −nk 0 are entangled while Φ + nk 0 and Φ − nk 0 are approximately separable). For a HO with the natural frequency Ω in a steady squeezed thermal state with the squeeze parameter r and temperature T , one has Q 2 = e −2r /(2Ω) coth(Ω/(2T )) and P 2 = e 2r ( Ω/2) coth(Ω/(2T )) provided that Q, P = 0 [48]. Thus we extract the squeeze parameter and the effective temperature by computing , (5.14) is the uncertainty function satisfying the uncertainty relation U ± ≥ /2. Note that the above definition of T ± eff (n) is identical to the effective temperature defined in Eq.(33) of Ref. [27], though the latter is for the UD detector rather than the field modes. Since V andr are large parameters, one has the approximated values (5.18) in the radiation zone. To eliminate the dependence on the observer's parameters V and r, we define the normalized squeeze parameter and the normalized uncertainty relation in the radiation zone as Their values for the same case in Figure 4 are shown in Figure 10.
In the left plot of Figure 10, we find that as n increases, [r + ] drops from the maximum value at n = 1, which is positive, then becomes negative for n ≥ 99, and reaches the minimum at n = 141. [r − ] increases from a very small value at n = 1, reaches the maximum value at n = 35. It is positive for all n. [U + − ( /2)] starts from the maximum value at The values of the effective temperature T ± eff (n) depend on the parameters of the combined system as well as those of the observer (such as V andr). When these parameters are fixed, in the inset of Figure 10 (left), one can see that T ± eff (n) behaves differently from [r ± ] and [U ± − ( /2)] (monotonically increasing in the same domain of n). From (5.15) the effective temperature T ± eff (n) for U + is higher than the one for U − if and only if U + > U − . Compare Figure 10 with Figure 4, we find that the most (or least) squeezed or uncertain harmonics (n > 1) are quite different from those harmonics contributing the maximal or minimal radiated power. So far we do not have any simple equality connecting the effective temperature of the harmonics T ± eff (n) and the effective temperature read off from the frequency spectrum of the radiated power (T B ≈ā) in Section 4.1, or the effective temperature experienced by the UD HO detector (T eff ≈ā/(2π)) in [18].

Summary
We calculated the radiated power of an Unruh-DeWitt harmonic-oscillator detector in a linear oscillatory motion in the Minkowski vacuum of a massless scalar field in (3+1)D. The radiation field is determined in the radiation zone for a laboratory observer and is independent of the cutoff corresponding to the time resolution of the detector. We find that there are two pulses in a cycle of the oscillatory motion, one is concentrated around θ = 0, the other around θ = π, and is highly compressed in the time domain when the detector motion is in the relativistic regime. There are periods of negative radiated power in each cycle at each fixed observing angle in the laboratory frame. The negative radiated power is due to the interference between vacuum fluctuations driving the HO and the retarded field emitted by the driven HO, which tends to suppress the signal of the Unruh effect in the radiation. Anyway, the averaged radiated power over a cycle is always positive as guaranteed by the quantum inequalities.
While the detector in the oscillatory motion is emitting the radiation, the effective temperature inside the detector is close to the Unruh temperature with the averaged proper acceleration. If we replace the correlators at the effective temperature by those at zero temperature, the result of the radiated power will be lower. The difference becomes significant when the detector is moving at a higher averaged acceleration and a shorter oscillating cycle, even with the interference terms included. This suggests that one may be able to distinguish the Unruh effect in the radiation by an Unruh-DeWitt detector in highly non-equilibrium conditions.
The compression of the pulse shape in the time domain corresponds to the coherent high harmonic generation in the frequency domain. Strong correlations can be found between the pulses in the time domain, and between the harmonics in the frequency domain. The correlation between two radiation fields emitted at two different retarded times τ and τ decreases as the interval of the retarded times |τ − τ | increases, while the correlation between two radiation field harmonics at two different frequencies nω 0 and n ω 0 decreases as |n−n | increases. At late times the correlations shows quasi-periodicity in t+t direction in the time domain, which can be related to the down conversion of the field quanta in the frequency domain.
We further construct an asymptotic reduced state of the harmonics of the field around θ = 0 and π in the radiation zone. We find that the harmonic modes of different wavelengths are approximately separable, due to the overwhelming largeness of the vacuum correlators of the free field. For each pair of field modes of the same harmonic frequency but in opposite directions, the reduced state looks like a two-mode squeezed thermal state.
distances from the detector are dominated by the v-parts [17], with the initial moment τ I → −∞ and the interaction duration η = τ − τ I → ∞.
For a UD detector in the model (2.1) uniformly accelerated along the worldline z µ U A (τ ) = (a −1 sinh aτ, 0, 0, a −1 cosh aτ ) with the proper acceleration a, the v-parts of the symmetrized correlators of the detector (with "mass" m 0 = 1 in [17]) in the Minkowski vacuum are given by [25] Q (η),Q(η ) U A Re ae −aτ −γτ − ie −iΩτ a + γ + iΩ F γ+iΩ (e −aτ ) and we have the initial moment τ I → τ I [17,27]. The same logarithmic divergence arises in the last line of (A.1) in the coincidence limit τ → τ , and we introduce another cutoff Λ 1 = −γ E − ln Ω|τ − τ | to control it. In the coincidence limit τ → τ for Q (η),Q(η ) U A v in (A.1), however, the Λ 1 term is purely imaginary and so should be dropped. For P (η),Q(η ) U A v and Q (η),P (η ) U A v , the situation is similar, and Λ 1 is present only in lim τ →τ P (η),P (η ) U A v . Then the above correlators in the coincidence limit reduce to the results given in [17,27].
The two-point correlators of a detector in uniform motion can be obtained by letting a → 0 and applying with finite η > 0 to (A.1). The two-point correlators of the field in the presence of the above uniformly accelerated detector going along z µ U A (τ ) can also be written in closed form. For the naive terms R [1] (x),R [1] (x ) U A v and their derivatives (R = Φ, Π), one can insert the above two-point correlators of the detector to Eqs. (2.9) and (2.10). For the interference terms, the simplest one is where η ± (x) ≡ τ ± (x) − τ I with the advanced time τ + (x) and the retarded time τ − (x) of the worldline z µ U A (τ ) for the observer at x [17], F = 1 if x is in the Rindler wedge F for the worldline z µ U A (τ ), and F = 0 if x is in the L-or R-wedge. The function h f (X) is given by for X < 0. (A.7) According to (A.3), h 0 (X) has a logarithmic divergence as X → 0, while h 1 (0) is regular. Other two-point correlators of the field can be obtained straightforwardly from (A.6) by Π [1]  with y = ±e −X . In the coincidence limit X → 0 the first term in the right hand side of (A.8) diverges like 1/X and so the derivative of (A.8) diverges like 1/X 2 . Fortunately these divergences are either purely imaginary thus absent in the symmetrized two-point correlators of the field, or canceled in the sum Φ [1]