Quantum Larmor radiation in de Sitter spacetime

We study the radiation emitted by inertial charge evolving on the expanding de Sitter spacetime. Performing a perturbative calculation, within scalar quantum electrodynamics (sQED), we obtain the transition amplitude for the process and using this we define the energy radiated by the source. In the non-relativistic limit we find that the leading term is compatible with the classical result (Larmor formula). The first quantum correction is found to be negative, a result which is in line with a number of similar quantum field theory results. For the ultra-relativistic case we find a logarithmic divergence of the emitted energy for large frequencies, which we link to the nature of the spacetime. We compare our results with that of Nomura et al. (2006), where the authors make a similar calculation for a general conformally flat spacetime.


Introduction
It is a well known result in classical electrodynamics that accelerated charges radiate. The emitted power is given by the famous Larmor formula [1]. The radiated energy in the case of non-relativistic motion of the source and with acceleration parallel to the velocity, adjusted for units, can be written as It is expected that the same result can be recovered from quantum theory in the limith → 0. Indeed, in Ref. [2], the authors obtained from sQED the lowest order contribution as being in agreement with the Larmor formula. The authors considered two distinct cases a e-mail: robert.blaga@e-uvt.ro b e-mail: sergiu.busuioc@e-uvt.ro of external electromagnetic fields that give rise to the same classical acceleration. Interestingly, although the leading term in both cases agrees with the classical result, the main quantum corrections differ. Similar results were obtained in Refs. [3,4] in a spatially homogeneous time-dependent electric field and electromagnetic plane-wave background. A distinct problem is the radiation of a charge in a time-dependent spacetime. In this case, in the GR picture, the source is inertial and the dynamic background plays the role of the external field. The problem was tackled in Refs. [5,6], in the general case of a conformally flat spacetime, by using the WKB approximation for the mode functions. The authors found that the leading term reproduces exactly the relativistic version of (1), when the trajectory is expressed in terms of conformal time.
Because of its privileged position in cosmological physics, the case of the de Sitter spacetime deserves a separate, more detailed treatment. Our goal in this paper is to obtain the term corresponding to the classical radiation and to calculate the leading quantum corrections for the energy radiated by a charge evolving on the expanding de Sitter spacetime (dS). We approach the problem with a perturbative calculation within sQED. We derive the radiated energy from the 1st order transition amplitude of the process which is analogous to the classical one.
One might wonder why there is radiation at all, given that the source is inertial (i.e. follows a geodesic trajectory). The motion of charges in gravitational fields has produced some controversy over the past decades, resulting in a considerable amount of literature on the subject [7,8,9]. The peculiarities of the problem are nicely illustrated by Chiao's paradox [10,11]. The question asked by Chiao is the following: will a charge on a circular orbit around a planet radiate and thus spiral inwards, as Newtonian intuition predicts, or continue moving along the geodesic, in accordance with the equivalence principle ? The paradox can be solved by noting that the equivalence principle has only local validity, while an electromagnetic charge along with its field is an extended object. "The Coulomb field of the particle, as it sweeps over the 'bumps' in spacetime, receives 'jolts' that are propagated back to the particle.
[...] The radiated effect comes from the work performed by this force." [12] An important feature of this radiation is that it is observer dependent. The classical example is that of the uniformly accelerated charge in flat space. While an inertial observer sees the charge radiating according to the Larmor formula, a co-accelerated observer will detect no radiation. [13] A similar situation arises in de Sitter space for comoving versus noncomoving observers. On physical grounds we expect that similarly to the uniformly accelerated case [12], the radiation reaction on a charge in dS cancels out, leaving the particle on the initial (geodesic) trajectory. The rule of thumb is: if there is variation in the local (physical) momentum of the charge in the relative motion with respect to the observer, there will be radiation.
The paper is structured as follows: In sec.II we gather the basics of sQED on de Sitter spacetime. In sec.III, starting from the transition amplitude, we define the energy radiated through the process. We obtain an asymptotic form for the energy in a weak gravitational field and proceed by expanding the result for different regimes of motion of the source. In the non-relativistic limit, we obtain the leading term and the first quantum contributions to the radiated energy. We obtain also a closed form for the energy in the ultra-relativistic limit. We find that the total emitted energy in this case is plagued by divergences, which is a typical feature of dS. In sec.IV we summarize and discuss our main results.
We work in natural units whereh = c = 1.

Basics of sQED on dS
The expanding patch of the de Sitter spacetime is described by the line element where ω is the Hubble constant, and we have introduced for convenience the time parameter ωη = e −ωt , η ∈ (0, ∞), with opposite sign as compared to what one usually calls conformal time.
The scalar modes that define the Bunch-Davies vacuum are [14,21] where ν = i µ 2 − 9/4, and m = ωµ is the mass of the scalar field. The momenta appearing in (3)  Given that the electromagnetic field is conformally invariant, the covariant components of the field are identical to their Minkowskian counterparts, while the contravariant components can be obtained by raising the indices with the metric tensor. The mode functions for the Maxwell field are 1 [23]: Following Ref. [17] we work with the electromagnetic field in the Coulomb gauge, given by A 0 = 0, √ −gA i ;i = 0. Given the mode functions (3) and (4), we can write the usual mode expansion for the scalar and Maxwell fields as: where a † (p), b † (p), c † (k) are the respective creation operators for the scalar particles, antiparticles and photons.
We are interested here only in tree-level QED processes, that are generated by the 1st order term in the expansion of the S-matrix: We have dropped the four-point interaction term because it does not contribute to the process studied here. We use a general prescription for interacting fields, following the classical textbook [14]. A detailed treatment of sQED on dS, including the reduction mechanism, can be found in Ref. [22]. The transition amplitude can be written as:

Quantum radiation of scalar charges
We are interested in the quantum theoretical counterpart of a charge emitting electromagnetic radiation given an external influence. In our case the charge is inertial, the expanding background playing the role of the external influence. The setup is as follows: in the initial state there is one scalar particle with momentum p , and in the final state we have a photon with momentum k and an arbitrary state of the scalar field. We average over all configurations that are indistinguishable from the point of view of a detector measuring the emitted radiation. This means basically summing over all possible final states of the scalar field.
where a and b * represent the number of scalar particles and antiparticles. Notice that the quantity (8) is independent of the definition of particles in the out state. Indeed this is the case because we can factor out an identity in (8), which can in turn be replaced by any complete orthonormal basis. The most natural way to proceed is in fact to insert an in basis (built from the Bunch-Davies modes (3)), which then truncates the sum at a finite number of terms. In our case we are left with the following terms: The first term represents a particle emitting a photon, while the second term represents the particle passing through without interacting, accompanied by the production of a pair and a photon from the vacuum 2 . The 2 Note that momentum conservation constrains the momenta p + k = p in the first, and k + q + q = 0 in the second process.
two configurations are illustrated by the Feynman diagrams Fig.1 . A very important observation is that the second process yields homogeneous and isotropic radiation. In an experimental context we can imagine that the detector can be adjusted to account for this background radiation. We can then drop this contribution and focus on the first term only.

Transition amplitude and radiated energy
The amplitude corresponding to the process depicted in Fig.1a was obtained in ref. [19] by one of the authors, and equals: In Ref. [19] we numerically analysed the result of the temporal integration and found a closed analytical form for the amplitude. Here instead we are interested only in the weak gravitational field limit (m/ω → ∞). With this in mind, we search for an asymptotic expression of the amplitude, in order to obtain the emitted energy as a power series in the Hubble constant ω. The energy emitted through the process can be computed as the energy of a photonhk, weighted with the probability of emitting a photon with the corresponding momentum. The expression for the energy can be written as: where V is the conformal volume, that will cancel the δ(0) term from the amplitude via the usual trick. Making use of momentum conservation p = p + k, the polarization term in (10) gives: Integrating over the final momentum of the source and with the use of (12), the radiated energy becomes: In what follows we attempt to find an asymptotic form for the temporal integral in (13).

Asymptotic expression in weak gravitational field
We seek an asymptotic expression for the radiated energy in the weak gravitational field regime. The idea is to obtain the energy as a series in powers of the Hubble constant ω. The leading term should be independent of h, so that we can consider it the "classical" radiation, i.e. it should reproduce the result obtained from classical electrodynamics. Our expectation is enforced by the results obtained in ref. [5] for a conformally flat universe (of which dS is a particular instance of). The calculation in [5] was performed in the WKB approximation, and the condition for weak gravitational field (µ → ∞) indeed assures that the WKB condition is fulfilled in our case also.
To obtain an asymptotic expansion of (10) we start by writing the Hankel functions H (1) ν , H (2) ν , in terms of modified Bessel functions K ν [20]: Using the property of the modified Bessel functions we write the product of Hankel functions as follows: Next, we use a large argument expansion [20]: where: which holds uniformly for | arg z| < 1 2 π when ν → ∞ 3 . The sign of the indices in (16) has been taken so that the condition on arg z is always fulfilled.
Substituting the expansion (17) into (16) and keeping only terms up to order O( 1 µ 2 ), we obtain: where z = p η µ , z = pη µ and we have considered ν iµ. The temporal integral with the expansion (18) can not be solved analytically. To continue, we need to further expand the asymptotic formula for small and large values of z. By observing that we can properly consider p phys m to be a nonrelativistic approximation, while p phys m represents an ultra-relativistic limit.

Radiation in the non-relativistic limit
First we discuss the radiation in the non-relativistic limit (z 1). Expanding all functions around small z and again keeping terms only up to order O 1 µ 2 , the asymptotic expression (18) reduces to: With the help of (20), we can now compute the squared absolute value of the temporal integral from the expression of the energy (13): Gathering all terms we can now obtain via eq.(11) the energy emitted under a unit solid angle and frequency: where the integration over the final momentum p was rendered trivial due to the Dirac delta function (p 2 = p 2 + k 2 − 2kp cos θ).
We note that all integrals are of the following form: This leads us to the resulted angular distribution of emitted energy: Plotting eq.(24) we observe: a) the characteristic sin 2 distribution for the radiation in the case of vanishing momentum of the source, b) as we increase the source momentum the radiation is emitted in a narrowing cone around the direction of motion, c) increasing amount of radiation in the backward direction. In order for formula (24) to remain valid, the 1/µ corrections must remain small as compared to the leading term. We require thus that p / µ. A further integration over dΩ gives us the total energy emitted in the process: If we write E = E cl + E (1) + E (2) we can identify from eq.(25) the lowest order term as Guided by the results of Ref. [5], we consider the acceleration to bë For our non-relativistic approximation this gives The remaining factor of 1 4 is due to the presence of the adiabatic cut-off. When we take the limit → 0 the energy diverges. This can be understood as follows: the role of the cut-off is to decouple the fields and thus halt the interaction on time scales larger than 1/ . When we take the vanishing limit this is equivalent to considering an infinite interaction time. Then, the energy radiated with a constant rate, under an infinite time, will be infinite. This also holds for a constantly accelerated charge in flat space. The results are consistent with that of Ref. [2].
Interestingly, if we naively take the non-relativistic limit in the results of Ref. [5] and also consider the adiabatic cut-off, we would obtain a result that is twice larger than (26). This is due to the fact that their calculation was tailored for a conformally flat spacetime with the conformal time ranging over the complete real axis. For the particular case of dS, this would mean the global de Sitter space. A similar situation was reported in Refs. [21,22,23] for Coulomb scattering in the expanding de Sitter space. For the expanding patch of dS, described by the line element (2), the calculation in Ref. [5] breaks down in eq. (30) where the boundary terms were neglected and in the subsequent integration over frequencies. If we were instead to consider the nonrelativistic approximation (z 1), by neglecting from the beginning quantities of order (p/p 0 ) 2 z 2 and with the adiabatic cut-off, the results would be identical to ours.
The leading quantum correction to the emitted energy is A negative quantum correction was also reported in all similar studies [2,3,4,5,6], for charges evolving in external electromagnetic and gravitational fields. The fact that the quantum effect suppresses the classical result thus seems to be a generic feature in such contexts. In Refs. [3,5,6] it is noted that the quantum corrections arise due to a non-local integration in time over the classical trajectory. In our case the trajectory is fixed, with constant accelerationẍ(η) = pω m and the non-locality is implicit in the result. On the other hand, in Ref. [4], the authors do not find the aforementioned nonlocality for the case of a charge moving in an electromagnetic plane-wave background. The difference is that this calculation is performed using the Schwinger-Keldysh (inin) formalism. It remains an open question why this difference arises. It will be an interesting subject for future work to calculate the radiation of a charge in de Sitter space using the in-in formalism and to compare with the results obtained in this paper .
In Ref. [6] it is found that the first correction to the radiation of a charge moving in a conformally flat background contains third derivative terms. Up to the orders that we have considered in our case we have ...
x (η) 0. The fact that we have a non-zero first order contribution thus suggests that our method captures terms that the WKB approximation misses.

Radiation in the ultra-relativistic limit
In this section we examine the behavior of the probability and the emitted energy through the process in the ultra-relativistic limit. Starting from (18) and imposing the condition for ultra-relativistic motion of the source (z 1), we obtain: The energy radiated under a unit solid angle and in unit frequency thus becomes: Integrating over the momenta of the photon we obtain the total energy emitted in the process: By changing the integration variable to p = p 2 + k 2 − 2p k cos θ, the angular integral becomes: A further change of variable to z = p − p + k results in: The indefinite integral over z has the following result: Using the notation introduced above we can write the energy emitted in the process as: In fig.(3) we have plotted the integrand of (36), which is the frequency distribution of the energy. The bulk of the radiation is emitted under frequencies k ≤ p as one would expect on physical grounds. For a small cut-off we see that most of the radiation is emitted for small frequencies. Increasing the parameter reveals that there is actually another competing channel around k p . This is not present in the non-relativistic case. We can understand this as follows: because we are investigating the process under weak gravitational field conditions, there is a loose energy conservation principle at action, which is reminiscent from flat space. For the non-relativistic case, where the energies go as ∼ p 2 , the photon momentum cannot compete with the source, and thus the only route towards energy conservation is p p , k → 0. On the other hand in the ultrarelativistic limit, because the energies go as ∼ p, the energy of the photon is on the same footing as the energy of the source, and the channel with k p , p → 0 becomes relevant. Thus we understand the peak at k p as arising from an interplay between the gravitational field, which gently lifts the energy conservation constraint, and the relativistic regime, which puts the energy of the radiation on a par with that of the source.
For large frequencies k > p we have a tail that fallsoff as 1/k, which leads to a logarithmic divergence when integrated over. The presence of this divergence is intimately linked to the famous divergence problem of de Sitter space [24]. We can understand it as a symptom of the finite integration over conformal time in (13). Because the Maxwell field is conformal, the photon effectively "lives" in conformal time and "feels" the limit η → 0 as being abrupt, although in the physical picture everything seems to be diluted away smoothly by the expansion of space. The finite limit for the temporal integration manifests like a finite-time sudden cutoff which leads to transitory effects, undesirable divergences and other artifacts [25,26,27,28,29].

Discussion
In this paper we have analysed the quantum radiation of a charge evolving on the expanding de Sitter spacetime. The emitted energy was derived from the transition amplitude of the corresponding sQED process. We compared the results of our perturbative calculation with that of Ref. [5], which was done in the WKB approximation. We have obtained the radiated energy as a power series in the Hubble constant, in the asymptotic case of a weak gravitational field. For a non-relativistic motion of the source, we have found the leading term to be compatible with the expected classical results. This is also identical to the results of Ref. [5], within the same approximation. Furthermore, the leading quantum correction was found to be negative, a result also reported in all similar studies. In the ultra-relativistic limit we expected to obtain a result which takes the form of the relativistic generalization of the Larmor formula. Instead we found that the energy has a logarithmic divergence for large frequencies. We interpret this as follows: the finite integration limit for the conformal time mimics a sudden decoupling of the interaction at time η → 0. Because this "event" happens under an arbitrarily small time interval, arbitrarily high frequency modes can get excited. Thus we also understand why this effect does not show up in the non-relativistic case, where only small frequency photons are emitted.
It would be interesting to see whether the above mentioned pathological fingerprint also shows up in a classical calculation, for the same setup. There are a number of papers that deal with the radiation of classical charges evolving on the global dS [30,31,32,33]. For the expanding patch of dS, the only study that we are aware of is done in Ref. [34]. We note that our results are compatible with that of Ref. [34], in that we find that comoving observers see no radiation. Indeed if we set p = 0 in eq.(25) and eq.(36) we find vanishing energy in both non-relativistic and relativistic cases. The situation is similar to the uniformly accelerated case in flat space. It was shown in Ref. [35] that if we consider the problem in a non-inertial (Rindler) reference frame: while the observers which are co-accelerated with the charge see no radiation, if there is mutual motion between the observer and the charge in the Rindler frame, an energy-flux will be present. It would be interesting to do a systematic study in the lines of Ref. [35], of the classical radiation emitted by charges on arbitrary tra-jectories on the expanding dS. Also it would be interesting to see how our results change if we consider proper Dirac electrons.
One more thing is worth noting. It is a pleasing fact that out of all 1st order processes the one studied here is the only one that falls off as an inverse power of µ as we go towards the flat space limit. As we have also signaled in Ref. [36], the probabilities for all other 1st order processes are exponentially suppressed as e −α(θ)µ , including the one depicted in Fig.(1b). This is linked to the fact that the process in Fig.(1a), that forms the object of this study, is the only one that has a classical analogue.