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. (JCAP 11:013, 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, it being in agreement with the Larmor formula. The authors considered two distinct cases 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 difa e-mail: robert.blaga@e-uvt.ro b e-mail: sergiu.busuioc@e-uvt.ro fer. 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 timedependent 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. it 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 per-formed 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 non-comoving 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 Sect. 2 we gather the basics of sQED on de Sitter spacetime. In Sect. 3, 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 Sect. 4 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) are the conformal momenta which are related to the physical momenta by p = 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, 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, which 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 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 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 it equals 2 Note that momentum conservation constrains the momenta p + k = p in the first, and k + q + q = 0 in the second process.
In Ref. [19] we numerically analyzed 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, which 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 the following 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 ofh, 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 ensures 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 ξ = 1 + z 2 + ln z 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) cannot 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 non-relativistic approximation, while p phys m represents an ultrarelativistic 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 resulting angular distribution of the 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 Eq. (24) to remain valid, the 1/μ corrections must remain small as compared to the leading term. We require thus that p / μ (Fig. 2). 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 cutoff, 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 non-relativistic 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 nonlocality is implicit in the result. On the other hand, in Ref. [4], the authors do not find the aforementioned non-locality for the case of a charge moving in an electromagnetic planewave background. The difference is that this calculation is performed using the Schwinger-Keldysh (in-in) 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 in-in formalism with 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 ultrarelativistic motion of the source (z 1), we obtain The energy radiated under a unit solid angle and in a 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 − 2 p 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: Frequency distribution of the radiated energy in the ultrarelativistic limit for p = 1. For large frequencies the radiation falls off as 1/k 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 nonrelativistic 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 toward energy conservation is p p , k → 0. On the other hand in the ultra-relativistic 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 falls off 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 analyzed 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 nonrelativistic 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 Eqs. (25) and (36) we find a vanishing energy in both the non-relativistic and the 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 trajectories 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 toward 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, which forms the object of this study, is the only one that has a classical analog.