Accelerated-Cherenkov radiation and signatures of radiation reaction

TL;DR

This paper investigates how acceleration and recoil influence Cherenkov radiation emitted by a charged particle in an optical medium, revealing a resonance at the superluminal threshold that could enable precise tests of radiation reaction effects.

Contribution

It introduces a theoretical framework combining acceleration and recoil effects in Cherenkov radiation, highlighting a resonance phenomenon useful for testing radiation reaction.

Findings

Recoil and acceleration produce a sharp resonance in Cherenkov emission.

A low-frequency spectral cutoff appears due to acceleration.

The resonance peak can be used for high-precision radiation reaction tests.

Abstract

In this manuscript we examine an accelerated charged particle moving through an optical medium, and explore the emission of accelerated-Cherenkov radiation. The particle's reaction to acceleration creates a low-frequency spectral cutoff in the Cherenkov emission that has a sharp resonance at the superluminal threshold. Moreover, the effect of recoil on the radiation is incorporated kinematically through the use of an Unruh-DeWitt detector by setting an energy gap, i.e., the change in electron energy, to the recoil energy of the emitted photon. The simultaneous presence of recoil and acceleration conspire to produce a localized resonance peak in the emission. These theoretical considerations could be used to construct high precision tests of radiation reaction using Cherenkov emission under acceleration.

Equations

$$\Gamma =q^2\int d\xi e^{-i\Delta E\xi /\gamma }{V}_{\mu \nu }[x^{\prime },x]G^{\mu \nu }[x^{\prime },x].$$

$$\frac{d\Gamma }{d\omega }=\alpha \beta \left[1-{\cos }^2({\theta }_c)\right].$$

$$\Delta E=\frac{(n^2-1){\omega }^2}{2m}.$$

$$\cos {\theta }_c=\frac{1}{n\beta }+\frac{(n^2-1)\omega }{2m\gamma n\beta }.$$

$$\frac{d\Gamma }{d\omega }=\alpha \beta \left[1-\frac{1}{n^2{\beta }^2}\left(1+\frac{(n^2-1)\omega }{m\gamma }+\frac{(n^2-1{)}^2{\omega }^2}{4m^2{\gamma }^2}\right)\right].$$

$$\Gamma =q^2\int d\xi e^{-i\Delta E\xi }{U}_{\mu \nu }[x^{\prime },x]G^{\mu \nu }[x^{\prime },x].$$

$$\frac{d\Gamma }{d\omega }$$

$$\cos ({\theta }_a)=\sqrt{\cos {{\theta }_c}^2+{\left(\frac{a{\gamma }^2}{n\beta \omega }\right)}^2}.$$

$${\omega }_c=\frac{a\gamma }{\sqrt{2}n{\beta }^2}\frac{1}{\sqrt{1-{\cos }^2({\theta }_c)}}.$$

$${\omega }_c$$

$${\omega }_{c^{\prime }}$$

$${\omega }_p={\left[\frac{2m\gamma {\omega }_a^2}{n^2-1}\right]}^{1/3}.$$

$${\omega }_c$$

$${\hat{S}}_I=\int d^{4}x{\hat{j}}_{\mu }(x){\hat{A}}^{\mu }(x).$$

$$\hat{\psi }(x,t)$$

$$\hat{\bar{\psi }}(x,t)$$

$${\hat{j}}_{\mu }(x)=\bar{u}(p){\gamma }_{\mu }u(p){\varphi }_{E_f}^{\ast }(x,t){\hat{a}}_p^{†}{\hat{a}}_p{\varphi }_{E_i}(x,t).$$

$${\hat{j}}_{\mu }(x)=u_{\mu }{g}^{\ast }(x)g(x)e^{i{E}_f\tau }{\hat{a}}_p^{†}{\hat{a}}_p{e}^{-i{E}_i\tau }.$$

$$e^{i{E}_f\tau }{\hat{a}}_p^{†}{\hat{a}}_p{e}^{-i{E}_i\tau }$$

$${\hat{j}}_{\mu }(x)=u_{\mu }\hat{q}(\tau ){\delta }^3(x-x_{tr}).$$

$$\mathcal{A}=i⟨\mathbf{k}|\otimes ⟨E_f|{\hat{S}}_I|E_i⟩\otimes |0⟩.$$

$${\hat{A}}^{\mu }(x)$$

$$\mathcal{P}$$

$$\displaystyle\Big{|}\bra{E_{f}}\hat{j}_{\mu}(x)\ket{E_{i}}\Big{|}^{2}$$

$$\int d^{3}k|⟨\mathbf{k}|{\hat{A}}^{\mu }(x)|0⟩{|}^2=⟨0|{\hat{A}}^{†\nu }(x^{\prime }){\hat{A}}^{\mu }(x)|0⟩.$$

$$⟨0|{\hat{A}}^{†\nu }(x^{\prime }){\hat{A}}^{\mu }(x)|0⟩$$

$$\mathcal{P}$$

$$\Gamma =q^2\frac{1}{(2\pi {)}^3}\frac{1}{2}\int d\xi \int \frac{d^{3}k}{\omega }\frac{u_{\mu }(x)u_{\nu }(x^{\prime }){\sum }_{ij}{ϵ}_i^{\mu }{ϵ}_j^{{\nu }^{\prime }}}{u_0(x)u_0(x^{\prime })}{e}^{-i(\frac{\Delta E}{\gamma }\xi -\mathbf{k}\cdot \Delta {\mathbf{x}}_{tr}(\xi )+\omega \xi )}.$$

$$\Gamma$$

Full text

I Introduction

The Unruh-DeWitt detector unruh1 ; dewitt has found a variety of applications across several fields of physics. Its most famous manifestation, and original use, was in the discovery of the Unruh effect unruh1 . There, the transition rate of a uniformly accelerated Unruh-DeWitt detector, with proper acceleration $\bar{a}$, implied that Minkowski vacuum was comprised of a thermalized bath of particles at temperature $T=\frac{\bar{a}}{2\pi}$. In addition to its original use in exploring the particle content of general relativistic space times, it has also played a key role in the exploration of relativistic quantum information louko ; bruschi ; hu ; bruschi1 ; lee . This wide ranging applicability is born out of the fact that it enables an exploration of particle systems independent of the major details, i.e. it only cares about the change in energy. This has the effect of taking an otherwise complicated process and rendering it a two level system like an atom with one excited state or a “qubit” of information. By ignoring the fine details and focusing on the particle’s quantized recoil due to its radiation, the Unruh-Dewitt detector approach may be able to provide a new perspective on a long standing question: the problem of radiation reaction dirac . Put simply, radiation reaction is the effect that photon emission has on a charged particle’s trajectory while it radiates, e.g. a recoil correction due to the emitted photon’s momentum. This is particularly exciting since recent theoretical and experimental works involving high intensity lasers and channeling radiation thomas ; cole ; poder ; wistisen ; chen have reported the first signatures of radiation reaction which appears to have quantum aspects involved in them dipiazza ; yaron ; burton ; iiderton .

Here we develop the Unruh-DeWitt detector formalism for its use in incorporating recoil and apply it to Cherenkov radiation cherenkov . Moreover, we also include acceleration in the analysis to gain insight into the effects of radiation reaction beyond the conventional models that consider the recoil as a result of Larmor radiation jackson . Fundamentally, this method is based on quantum field theory in curved spacetime for analyzing radiation emission davies ; winitzki . By incorporating the techniques used to explore radiation in more exotic spacetimes, we can gain considerable insight into the simpler setting of Minkowski space. To explore vector currents coupled to photons via the QED interaction, we use spacetime trajectories/world lines that are often utilized in quantum field theory in curved spacetime to investigate the emission of radiation parker ; hawking ; unruh1 ; muller ; matsas1 ; matsas2 ; matsas3 ; lynch ; lynch1 and apply them to superluminal velocities with acceleration in an optical medium, revealing a novel accelerated-Cherenkov effect. We derived a generalized Frank-Tamm cherenkov ; frank ; ido formula that simultaneously takes into account the quantum recoil due to a single photon emission as well as the acceleration; the inclusion of both effects yields a novel resonant spectral cutoff at the superluminal threshold. We also demonstrate how to incorporate the recoil, e.g. radiation reaction, by using the recoil momenta of the emitted photon as the change in energy of the electron. The preliminary example of Cherenkov radiation is given to outline the functionality of the technique, how to incorporate the recoil, and is then applied to the case where an electron is accelerating in the superluminal regime. The presence of the acceleration sets an energy scale that sharply cuts off the photon emission. This energy scale is quite tunable due to the presence of a strong resonance as the electron velocity approaches the in-medium speed of light. Starting from the cutoff energy, the emitted photon spectra rapidly climbs and asymptotes to the flat Cherenkov emission spectra. The incorporation of the particle recoil due to the emission of a single photon quanta affects the higher frequencies to produce a monotonically decreasing spectra. The combination of both yields a peak in the Cherenkov emission spectrum which may be measurable experimentally kimura ; yhu ; henstridge .

Section 2 of this paper explores photon emission in QED peskin and its relationship to an Unruh-DeWitt detector with the only assumption being that the electron momentum is much bigger than the emitted photon momentum (i.e. the weak recoil approximation). By propagating the detector along an inertial, i.e. non accelerated, world line moving superluminally through an optical medium, we are then able to analyze Cherenkov emission. The overall computational technique is presented and we show how to take into account the quantum recoil of photon emission by setting the Unruh-DeWitt detector energy to the recoil energy produced by the photon emission. Then, in Sec. III we retrace the computation except we move the detector along a uniformly accelerated world line. The effects of incorporating acceleration, including the particle recoil due to the interaction with the second-quantized photon field, is then explored. Strong signatures of the acceleration can be tuned by a resonance at the superluminal threshold and its combined effect with recoil produces a broad peak in the emission spectra. Finally, a summary of conclusions are presented and we have included all pertinent derivations in the appendix. All calculations are performed using the natural units $\hbar=c=k_{B}=1$.

II Cherenkov Emission via Unruh-DeWitt detectors

The emission of radiation produced by inertial currents offers the easiest setup to analyze recoil since the kinematics are purely Minkowskian. This setting provides a rather simple introduction to both radiation emission as well as the recoil. We begin by examining the emission of a photon in a refractive medium, $\hat{A}^{\mu}(x)$, using the current interaction for QED peskin , $\hat{S}_{I}=\int d^{4}x\hat{j}_{\mu}(x)\hat{A}^{\mu}(x)$. To model the electron current semi-classically we will make use the charged current coupled to an Unruh-DeWitt detector, $\hat{j}_{\mu}(x)=u_{\mu}\hat{q}(\tau)\delta^{3}(x-x_{tr})$. Here we defined the Heisenberg evolved charge monopole moment operator $\hat{q}=e^{i\hat{H}\tau}\hat{q}(0)e^{-i\hat{H}\tau}$ where $\hat{q}(0)$ is defined as $\hat{q}(0)\ket{E_{i}}=\ket{E_{f}}$ with $E_{i}$ and $E_{f}$ the initial energy and final energy of the electron moving along the trajectory, $x_{tr}$, of the current with four velocity $u_{\mu}$. The energy gap of the detector is $\Delta E=E_{f}-E_{i}$ and the charge $q$ of the electron is given by $q=|\bra{E_{f}}\hat{q}(0)\ket{E_{i}}|$ davies ; unruh1 ; matsas1 ; matsas2 ; matsas3 .

In the case of superluminal motion through an indexed medium, the emission of radiation accompanied by the transition of an internal energy level is associated with the anomalous Doppler effect ginsburg . This extra degree of freedom will allow us to incorporate electronic excitations to higher energy, i.e. the recoil. With the intent to examine Cherenkov radiation, i.e. photon emission from a charge moving faster than the speed of light in medium, we formulate the following amplitude; $\mathcal{A}=i\bra{\mathbf{k}}\otimes\bra{E_{f}}\hat{S}_{I}\ket{E_{i}}\otimes\ket{0}.$ The square of the this quantity determines the transition probability per momentum of the final state photon. Note, by having a detector we are relieved of integrating over the final electron momentum states. Constructing the response function yields

[TABLE]

Here we are integrating over the laboratory time $\xi$. The “velocity tensor” $V_{\mu\nu}=v_{\mu}v_{\nu}$ determines the trajectory of the electron with $v^{\mu}=(1,0,0,\beta)$ given by the lab frame velocity. Computation of the emission spectrum yields the Frank-Tamm formula frank , see the appendix for details. Hence,

[TABLE]

We note the prefactor is nothing more than the fine structure constant $\alpha=\frac{q^{2}}{4\pi}$. This well know expression characterizes the emission of Cherenkov photons by charged particles in the superluminal regime. Note that it scales as $\beta$ and is independent of the frequency modulo the frequency dependence in the index of refraction. The Cherenkov emission angle also contains a contribution from the energy gap which can be used to explore the anomalous Doppler effect and other phenomena ginsburg ; xihang , $\cos{(\theta_{c})}=\frac{1}{n\beta}+\frac{\Delta E}{n\beta\omega\gamma}$. The energy gap serves to deform the emission angle making it smaller if the energy gap is positive (excitations) or larger if the energy gap is negative (decays). If we are examining the recoil from a photon emission we can also determine the explicit form for the energy gap. Considering an initial “dressed” electron with total energy comprised of the electron rest mass and the emitted photons energy $E_{i}=\sqrt{m^{2}+\omega^{2}}$ and a final state electron energy comprised of a bare electron and recoil momentum $k=n\omega$ (where n is the frequency-dependent index of refraction) produced by the photon emission, $E_{f}=\sqrt{m^{2}+(n\omega)^{2}}$, see Fig. 1 for an illustrative example. The recoil momentum, $n\omega$, is characteristic of the non-trivial dispersion relation of the photon. As such, we can determine the energy difference to be

[TABLE]

Here, the energy gap of the Unruh-DeWitt detector is defined as the difference between the energies of the electron before and after emission. This specific energy gap, Eq. (3), is present whether the photon is emitted at constant velocity in the superluminal regime or by acceleration. The energy gap is a by-product of mass renormalization where the electron’s energy can be enhanced on account of it propagating in a material with dispersion tsytovich and exists independent of whether or not the charged particle is accelerated. In general when in the presence of an optical medium one can always define the energy gap in such a way. There may, in fact, even be other terms which characterize other processes as well but we leave the examination of those to a later work. From this relation we also obtain the fully relativistic quantum corrected Cherenkov emission angle ginsburg . Our Cherenkov angle, including the recoil correction, is then given by

[TABLE]

As such, by defining the energy gap in accordance with mass renormalization, we have derived the Frank-Tamm formula generalized to the exact quantum result that includes the reaction of the particle to its recoil by the photon emission sokolov ; ido . Hence,

[TABLE]

To conclude this section, we have outlined the general procedure for the computation of the emission rate of Cherenkov radiation, including recoil, via the use of an Unruh-DeWitt detector. We outlined the general idea for how to incorporate recoil with the Unruh-DeWitt detector; this procedure is discussed in more detail in the appendix. What is most important is the fact that we are able to successfully incorporate recoil into the emission process. In the following section we will follow the same procedure while placing the current along an accelerated trajectory which is also superluminal. We then explore the emission spectrum of Cherenkov radiation, under acceleration, while simultaneously taking into account recoil and how these inherently quantum mechanical processes alter the emission spectrum.

III Cherenkov emission via accelerated currents

Building upon the formalism of the previous section, let us now endeavor to compute the emission spectra by an accelerated electron current in the Cherenkov regime. The overall prescription will be to utilize the kinematics of Rindler space to determine trajectories, velocities, etc. and incorporate their dynamics into the emission. Recalling the response function, Eq. (1), is given by

[TABLE]

Note that we are now integrating over the proper time, also parametrized by $\xi$, of the electron. To analyze accelerated Cherenkov emission we integrate over the proper time rather than laboratory time since the accelerated trajectories are parametrized this way. We will begin by noting the velocity tensor in now given by $U_{\mu\nu}=u_{\mu}u_{\nu}$. Recalling that under proper acceleration $\bar{a}$, the four-velocities at proper time $\tau$ will be given $u^{\mu}=(\cosh{(\bar{a}\tau)},0,0,\sinh{(\bar{a}\tau)})$ winitzki . Computation of the emission spectrum, parametrized by the lab frame acceleration $a$, with $\bar{a}=a\gamma^{3}$, then yields the result,

[TABLE]

This emission spectrum is the main result of this work and illustrates the nature of accelerated-Cherenkov radiation. Note that it reduces to the standard Frank-Tamm formula, Eq. (2), as the acceleration goes to zero. We also recall the standard Cherenkov angle, as derived in the previous section, is given by $\cos{\theta_{c}}=\frac{1}{n\beta}+\frac{\Delta E}{n\beta\omega\gamma}$. The presence of acceleration then yields the new acceleration-dependent Cherenkov angle,

[TABLE]

Note, not only have we successfully recovered the Frank-Tamm formula with an additional acceleration-dependent term but also have the added functionality of the Unruh-DeWitt and its ability to incorporate recoil ginsburg ; ido . Overall, we see the that presence of acceleration serves to not only change the emission spectrum but also the emission angle. These effects are not merely a byproduct of having a time-dependent velocity but exist independently of the standard Cherenkov effect entirely. One such effect is to produce a low frequency cutoff in the emission spectrum. When this effect is combined with the recoil correction, which produces a high energy cutoff, we will find a resonance or local maximum in the spectrum. Figure 2 below details the general structure of Cherenkov radiation with contributions from both acceleration and recoil. Note that both contributions are required to produce a local maximum. Here, and throughout the manuscript, we assume a constant index of refraction.

Without acceleration this resonance structure vanishes. We see that the recoil correction serves to give the Cherenkov spectrum a monotonically decreasing dependence on the frequency. Consequently, a tilted spectra is, modulo other effects such as material dispersion, a signature of recoil in Cherenkov emission. However, with acceleration present in the system we will see that there is clearly a local maximum in the spectrum which depends explicitly on the acceleration. This local maximum requires both acceleration and recoil to be present in the system and therefore may provide a new scenario to explore the phenomena of radiation reaction. To better characterize this process let us examine the peak and cutoff structure. The subsequent expressions, neglect the last term in the square root of Eq. 7, i.e. the cross term, so as to examine analytical forms of the dominant contributions in both the low and high frequency cutoff limits. As such, for the low energy cutoff, i.e. when the emission rate goes to zero, we find the following condition

[TABLE]

Without an internal energy level, i.e. no anomalous Doppler effect, we have the following cutoff

[TABLE]

Here we defined the characteristic energy scale set by the particle’s acceleration, $\omega_{a}=\frac{a\gamma}{\sqrt{2}\beta}$. Note the presence of a resonance at the Cherenkov condition. Moreover, the overall scale of the cutoff frequency from Eq. (10) is specific to having an accelerated particle, e.g. an electron, with characteristic energy $\omega_{a}$ in the Cherenkov regime. As with the Unruh effect unruh1 , the energy scale $\omega_{a}$ set by acceleration is often vanishingly small. The important point however, is that we have found that Cherenkov emission provides another arena to explore acceleration energy scales which may be more experimentally tenable than in the case of the Unruh effect. As an example, the resonance at the superluminal threshold can potentially be used to make this acceleration scale more accessible. It is important to note that an acceleration energy scale of $\omega_{a}=5.5\times 10^{-4}$ eV used as the characteristic energy scale in our figures is based off an acceleration of, $a\sim 1.8\times 10^{20}$ m/s2, that is indeed experimentally realizable in common laser systems with an amplitude of $\sim 1$ GV/m. Moreover, table top plasma wakefield accelerators are even capable of producing accelerations, $a\sim 10\times 10^{22}$ m/s2 rosenz . These systems may also provide a setting to explore this cutoff structure of Cherenkov emission. Unlike the above case with no internal energy levels, i.e. with $\Delta E=0$, the inclusion of an internal energy gap we will yield a quartic equation that determines the cutoff frequency in the presence of radiation reaction and acceleration. As before, the inclusion of recoil is accomplished by setting the energy gap to $\Delta E=\frac{\omega^{2}(n^{2}-1)}{2m}$ in the Cherenkov angle. The fact that we will have a quartic cutoff is due to the frequency dependence in the recoil correction. As such, the cutoff frequency is determined by the solution to the following equation,

[TABLE]

Moreover, the presence of recoil also produces a high frequency cutoff even in the absence of acceleration. From Eq. (7) we can determine the high frequency cutoff due to recoil (in agreement with ido ) to be,

[TABLE]

We can also compute the peak frequency that is produced by the simultaneous presence of acceleration and recoil. A local maxima in the Cherenkov spectrum in these systems therefore implies both acceleration and recoil. The peak itself is determined by computing the maximum of our generalized Frank-Tamm formula Eq. (7). It resultant expression is quartic in nature, $\frac{\omega_{a}^{2}}{\omega^{3}}=\frac{(n^{2}-1)}{2\gamma m}\left[1+\frac{\omega(n^{2}-1)}{2m\gamma}\right]$, but has an approximate solution given by the following equation,

[TABLE]

The photons of this peak frequency will be the only photons emitted at the Cherenkov threshold. Consequently, by modulating the electrons velocity and acceleration at or near the threshold, one can essentially tune the Cherenkov radiation to a unique frequency that is also emitted into a unique direction, i.e. controlling the electrons velocity can put the system at resonance and controlling the acceleration can tune the peak frequency emitted by setting the energy scale $\omega_{a}$. As such, with a sufficiently intense source of modulated electrons one could create a tunable Cherenkov emitter. In Fig. 3 we examine the structure of the cutoff frequency produced by both the recoil and the acceleration as an electrons velocity approaches Cherenkov threshold from above for an index of refraction of $n=1.5$. The peak resonance that is converged to is given by the solution of Eq. (13) and gives the peak frequency emitted for any superluminal velocity under the assumption of a constant index. As such, we note that by examining the solution to the cubic resonance equation, we can tune the dominant emission frequency by driving both the velocity and acceleration of a Cherenkov system at the superluminal threshold. Moreover, we note that the index of refraction will also depend on frequency. This can also be taken into account by insertion of the frequency dependent index of refraction into the cutoff and peak frequency equation and numerically solving that equation as well. Any vertical slice of Fig. 3 will yield the maximum frequency as well as the high and low frequency cutoffs. The full cutoff curve is determined by Eq. (11). The asymptotic low frequency cutoff produced by the acceleration is determined by Eq. (10). The asymptotic high frequency cutoff produced by the recoil is determined by Eq. (12). The maximum frequency emitted is determined by Eq. (13) and it should be noted that this maximum determines the unique frequency that is emitted at the $\beta$ threshold and can be used to “tune” the emission. This threshold is also shifted due to the presence of both acceleration and recoil. In addition, we can see from Fig. 3 that radiation emission is suppressed at higher energies due to the recoil. This further serves to ensure that the electron, modeled as an Unruh-DeWitt detector, remains pointlike with a wavelength much smaller than that of the emitted photons. If this were not the case then we would have to take into account smearing of the detector and its modification on the emission spectrum; see e.g. eduardo for a QED-based analysis. We should also comment that additional factors such as a strong dependence of the index of refraction on frequency or higher order processes could further modify the spectrum which would undoubtedly give rise to another rich area of investigation.

The ability to tune an electron’s velocity with high precision at or near the Cherenkov threshold would imply the ability to control the peak frequency emitted. It is this peak frequency that the low and high energy cutoffs converge to at the threshold limit and will thus determine the only frequency emitted at resonance. The broadness of this peak will be determined by the control of the electron velocity; a sharp velocity distribution in an electron beam will then yield a sharp Cherenkov emission spectrum. As we drive the electrons with specific velocity and acceleration profiles we will get different peak frequencies. Figure 4 shows the various peak frequencies as a function of velocity and acceleration at resonance; We substitute an index of refraction $n=1/\beta$ everywhere. Note that for initial electron velocities and accelerations that are currently accessible in electron microscopes we would expect to see a peak frequency in, or near, the lower end of the visible spectrum. To be clear, this analysis assumes a relatively constant index of refraction as a function of frequency in the vicinity of the threshold. More specifically, given a beam velocity $\beta$ we must then select a material with index of refraction that will place the beam at or near the threshold, i.e. $n\sim 1/\beta$. Then, we must also require that material has an index that is fairly constant with respect to frequency, near the peak, in the vicinity of the threshold. With these conditions met, it stands to reason that a signal of the peak frequency can be measured. We can see from Fig. 4 that for typical electron velocities in electron microscopes $\beta\sim 0.6-0.7$ and with laser field amplitudes currently in regular use $\sim 1$ GV/m, we should then see a signal between $.1-1$ eV. This signal is well within the measurement capabilities of electron microscopes, e.g. an infrared detector or indirectly via electron energy loss spectroscopy abajo . Moreover, these considerations can further add to the robustness of the Cherenkov effect’s applicability to not only particle identification ginis ; ido2 but also to radiation emission.

For the sake of completeness, we also consider the anomalous Doppler effect ginsburg ; xihang , i.e. an energy gap that does not depend on the emitted photons’ frequency, we then have $\cos{\theta_{c}}=\frac{1}{n\beta}+\frac{\Delta E}{n\beta\omega\gamma}$ and a frequency cutoff given by following roots of Eqn. (9);

[TABLE]

Again, we have a resonance but we see that the cutoff frequency can shift by any internal energy level. We also note the resonances can provide an additional method to probe the anomalous Doppler effect. Note that for both excitations and decays we potentially have two cutoff frequencies. This phenomenon exists both with and without acceleration. In the limit of zero energy gap we arrive back at the Cherenkov cutoff in the presence of acceleration, Eq. (10). The fact we have an acceleration energy scale $\omega_{a}$ is also reminiscent of the Unruh effect unruh1 . There is indeed a close relationship between the anomalous Doppler effect and the Unruh effect frolov1 ; frolov2 that has been explored in detail. The results presented here further stress the similarities between these two processes and outline the details of a new effect due to acceleration that is present in Cherenkov radiation.

Finally, in addition to the novel structure of the Cherenkov emission spectrum in the presence of acceleration, there is also a deformation of the Cherenkov cone. If we are to assume an index of refraction that is independent of frequency, then the Cherenkov radiation will all be emitted into the well defined cone; the angle relative to direction of prorogation is defined by $\cos{(\theta_{c})}=\frac{1}{n\beta}$. When acceleration is present, it allows the Cherenkov cone to fill in with the lower frequency end of the emission spectrum; see Fig. (5). Therefore, for a sufficiently thin sample, an accelerating charged particle in the superluminal regime, would reveal a filled-in circle rather than the standard ring in an imaging sensor.

To conclude this section, we have revealed a new aspect in the theory of Cherenkov emission which incorporates both acceleration as well as recoil. Its high tunability, accompanied by resonance peaks, may be employed as a high precision test of radiation reaction using Cherenkov emission under acceleration since the recoil is easily incorporated. We explored the various ways in which to characterize the radiation reaction, namely the acceleration-dependent cutoff frequency and local maxima in the spectrum. It is interesting to note the fact that we have an energy scale set by the acceleration is thus reminiscent of the Unruh effect matsas4 .

Further implications on the effect of radiation reaction from charge renormalization would be interesting to investigate as well. Our work showed that the recoil correction, which is consistent with mass renormalization, manifests itself as a change in the electron’s energy due to the recoil momentum imparted on the electron by the photon emission. The specific effect of charge renormalization peskin would not only change the effective charge of the electron but should, in principle, contribute a new term in the energy gap as well. The charge functions as the monopole moment operator in the interaction and any renormalization of the charge may have an explicit presence there. Another potential correction, which is not discussed in this work, is field renormalization. That could provide another promising avenue to investigate sources of radiation reaction. In the end, we have that the utility of the Unruh-DeWitt detector applied to the problem of radiation reaction may be of particular value.

IV conclusions

In this manuscript we have made use of an Unruh-DeWitt detector to examine the effect of recoil and acceleration on Cherenkov radiation. Overall, our investigation has determined the effect of acceleration on both the emission spectrum as well as the emission angle. The combined presence of recoil and acceleration was then used to gain insight into the paradoxical nature of radiation reaction. We found a resonant energy cutoff resulting from the acceleration and tilted spectra due to the recoil which conspire together to produce a local maximum in the emission spectra. We also found that the presence of acceleration serves to fill in the Cherenkov cone with the lower frequency end of the emission spectrum. Having carried out the analysis via the use of an Unruh-DeWitt detector, we hope this work may open the way to expand its utility and also help develop a better understanding for the phenomenon of radiation reaction crystal .

Acknowledgments

We would like to thank Daniel Abutbul and Yarden Sheffer for many stimulating discussions. M.L. was supported in part at the Technion by a Zuckerman fellowship. E.C. was supported by the Canada Research Chairs (CRC) Program.

V Appendix

V.1 Cherenkov Emission via Unruh-DeWitt detectors

The emission of radiation produced by inertial currents offers the easiest setup to analyze recoil since the kinematics are purely Minkowskian. This setting provides a rather simple introduction to both radiation emission as well as the recoil. We begin by examining the emission of a photon in a refractive medium, $\hat{A}^{\mu}(x)$, using the current interaction for QED peskin ,

[TABLE]

To model the electron current semi-classically we will make use of the weak recoil limit. To see how this is accomplished we begin by recalling that for spinors $u(s,p)$ and $v(s,p)$ of spin $s$ and momentum $p$ that are created by $\hat{a}^{\dagger}_{s,p}$ and $\hat{b}^{\dagger}_{s,p}$, we have the following electron field operators

[TABLE]

The positive and negative frequency modes are given by $\phi_{p}(x,t)$ and $\chi_{p}(x,t)$ respectively. Normally these modes are plane waves, however in more general spacetimes our only requirement is that they are positive and negative frequency modes with respect to the particle’s/detector’s proper time. Using these fields, we will formulate the electron current, $\hat{j}_{\mu}=\hat{\bar{\psi}}\gamma_{\mu}\hat{\psi}$. In order to enforce the weak recoil approximation we will ensure that our initial and final Fock states are labeled by their respective energies but have the same momentum, $p_{\mu}$. The legitimacy of this approximation is restricted to the regime where the energy of the emitted radiation is much less than the rest mass of emitting the particle, e.g. $\sim 1$ MeV for an electron. We will also neglect any spin effects which are described elsewhere, see e.g. ido . Focusing strictly on electrons, i.e. no antiparticles, our current then reduces to the following

[TABLE]

Now, in keeping with the weak recoil approximation we make use of the Gordon identity peskin . As such, we have $\bar{u}(p)\gamma_{\mu}u(p)=u_{\mu}$. We now make use of the fact that our positive and negative frequency mode solutions can be separated into their spatial and temporal components via $\phi(x,\tau)=g(x)e^{-iE\tau}$. We have chosen to parametrize our fields via the electron’s proper time, $\tau$, to incorporate the Rindler coordinate chart when analyzing the accelerated case davies . Our current now reduces to

[TABLE]

We note that for sufficiently localized electronic wave functions, e.g. with a wavelength much smaller than the wavelength of emitted radiation, we have $g^{\ast}(x)g(x)=\delta{(x-x_{tr})}$ along the classical trajectory of the electron; which is assumed to be uniform in the weak recoil limit. This can also be seen by examining a smeared Unruh-DeWitt detector, see e.g. eduardo . There, the smearing function which encodes a detector’s spatial extent is given by $F(x)=-i\psi_{f}(x)\nabla\psi_{i}(x)$. Applying this formula for the inertial plane waves that are typically used to model electrons in the Cherenkov effect yields $F(x)=pg^{\ast}_{f}(x)g_{i}(x)e^{i\Delta E\tau}$, where $p$ is the momentum of the electron and $\Delta E=E_{f}-E_{i}$. Then, by inspection of Eq. (A4), we see that our use of a delta function to model the semi-classical trajectory is equivalent to a detector with a sharply defined spatial extent. It would be interesting indeed to look for finite size effects, i.e. smearing, of an electron in the emission spectrum in processes like Cherenkov radiation. Finally, by attaching the time-dependence to the creation and annihilation operators we have

[TABLE]

Here we defined the Heisenberg evolved charge monopole moment operator $\hat{q}=e^{i\hat{H}\tau}\hat{q}(0)e^{-i\hat{H}\tau}$ where $\hat{q}(0)$ is defined as $\hat{q}(0)\ket{E_{i}}=\ket{E_{f}}$ with $E_{i}$ and $E_{f}$ the initial energy and final energy of the electron moving along the trajectory, $x_{tr}$, of the current. The energy gap of the detector is $\Delta E=E_{f}-E_{i}$ and the charge $q$ of the electron is given by $q=|\bra{E_{f}}\hat{q}(0)\ket{E_{i}}|$. As such, we have transformed our fermionic current into a semi-classical charged current coupled to an Unruh-DeWitt detector davies ; unruh1 ; matsas1 ; matsas2 ; matsas3

[TABLE]

This formalism allows us to examine transitions both up and down in energy. In the case of uniform motion through an indexed medium, transitions up in energy are associated with the anomalous Doppler effect ginsburg . This extra degree of freedom will allow us to incorporate electronic excitations to higher energy, i.e. the recoil. With the intent to examine Cherenkov radiation, i.e. photon emission from a charge moving faster than the speed of light in medium, we formulate the following amplitude;

[TABLE]

The square of the this quantity determines the transition probability per momentum of the final state photon. Note, by having a detector we are relieved of integrating over the final electron momentum states. The final momentum is essentially fixed by the detector gap and is tuned directly to the energy change experienced during the recoil. This is equivalent to a delta function density of states, $\frac{d\mathcal{P}}{d^{3}k}=\int d^{3}p_{f}\delta^{3}{(\mathbf{p}_{f}-\mathbf{p}_{r})}|\mathcal{A}|^{2}=|\mathcal{A}|^{2}.$ The nature of the recoil momentum $p_{r}$ will be explored later on when we define the Unruh-DeWitt detector energy gap. Using the mode decomposition for a massless vector field in a dielectric medium with index of refraction $n$, we have

[TABLE]

Here we have a dispersion relation which relates the frequency to momenta via $n\omega=k$ jackson . This non-trivial dispersion relation is at the heart of the recoil correction to radiation reaction and we will see how it leads to a nonzero recoil momentum by strict conservation of energy arguments; see section IV below. Constructing the probability from Eq. (A7) yields

[TABLE]

Evaluation of the square of the electron current matrix element yields

[TABLE]

Next, we shall evaluate the vector field inner product. For this we will need to integrate over the final state momenta, thereby developing the total emission probability. The result is the Wightman function davies for the photon. Hence,

[TABLE]

Here we have factored out the completeness relation, $\int dk\ket{k}\bra{k}=1$ over momentum eigenstates. The resultant two point function, with vector indices, factors into a scalar two point function weighted by the sum over indexed vector polarizations. Evaluation of this two point function, and being sure to keep polarization sum inside the integral, and without evaluating the integral over the momentum yields

[TABLE]

Combining all the pieces we can formulate the total transition probability. Hence,

[TABLE]

Note, in the last line we changed our integration variable from proper time to coordinate time via $dt=u_{0}d\tau$. Then, by defining the difference and average laboratory times, namely $\xi=t^{\prime}-t$ and $\eta=(t^{\prime}+t)/2$ respectively, we can decouple the temporal integrations lynch ; lynch1 . Note, this parameterization has a unit Jacobian, i.e. $dtdt^{\prime}=d\xi d\eta$. In the presence of trajectories that are linear in time, i.e. unaccelerated, the trajectory will obey the following relation $\Delta x=x(t^{\prime})-x(t)=x(t^{\prime}-t)=x(\xi)$. As a consequence we will have $\Delta x_{tr}=\Delta x_{tr}(\xi)$. Moreover, with a trajectory that is constant in velocity we will have $\epsilon_{i}^{\mu}$ and $u_{\mu}(x)$ independent of time but dependent on the photon emission angle. Noting the resultant expression is therefore independent of the laboratory time $\eta$ we can formulate the transition rate via $\frac{d\mathcal{P}}{d\eta}=\Gamma$. As such,

[TABLE]

We have defined the $\xi$-dependent Lorentz gamma $\gamma$ in the last line to boost the proper time difference to the lab frame. We can then define the scalar quantity $U=\frac{u_{\mu}(x)u_{\nu}(x^{\prime})\sum_{i,j}\epsilon_{i}^{\mu}\epsilon_{j}^{\nu}}{u_{0}(x)u_{0}(x^{\prime})}$. Moving into momentum space spherical coordinates, with the $z$-axis aligned along the direction of the velocity, we have

[TABLE]

Here we used the dispersion relation $\omega=k/n$. Now, for inertial spacetime trajectories under constant velocity with respect to the lab frame we have $\Delta x_{tr}(\xi)=\beta\xi$. As such we can integrate over time along the trajectory yielding

[TABLE]

Let us examine the scalar quantity $U$ in more detail. Since it is the 4-velocity of the charged particle contracted with the polarization vector we can evaluate it using the properties of the polarization. Since the polarization vector has no zero component and will always be orthogonal to the wave vector of the photon, we find

[TABLE]

Here we see the $\sin{(\theta)}$ rather than $\cos{(\theta)}$ comes out of the dot product. This is because the angle $\theta$ lies between the charged particle velocity and momentum of the photon. Since the polarization is everywhere perpendicular to the momentum we get $\cos{(\pi/2-\theta)}=\sin{(\theta)}$ from the dot product. As such, our emission rate now becomes

[TABLE]

Prior to integrating over the angle we note the presence of the delta function which enforces conservation of energy contains the celebrated Cherenkov emission angle deformed by the presence of the energy gap, this is the anomalous Doppler effect ginsburg ; xihang . Hence, the root yields the Cherenkov angle,

[TABLE]

If we are examining a recoil from a photon emission we can find an explicit form for the energy gap, see section IV below for further discussion. Considering an initial “dressed” electron with total energy comprised of the electron rest mass and the emitted photons energy $E_{i}=\sqrt{m^{2}+\omega^{2}}$ and a final state electron energy comprised of a bare electron and recoil momentum $k=n\omega$ (where n is the frequency-dependent index of refraction) produced by the photon emission, $E_{f}=\sqrt{m^{2}+(n\omega)^{2}}$, see Fig. 1 for an illustrative example. The recoil momentum, $n\omega$, is characteristic of the non-trivial dispersion relation of the photon. As such, we can determine the energy difference to be

[TABLE]

From this relation we also obtain the fully relativistic quantum corrected Cherenkov emission angle as computed by Ginzburg ginsburg . We also note this is consistent with mass renormalization where the electron’s energy contains losses on account of it propagating in a material with dispersion tsytovich . These losses manifest themselves as a frequency dependent change in the electron’s mass. Our Cherenkov angle, including the recoil correction, is then given by

[TABLE]

Next, to integrate over the angles we note the root of the delta function, either given as above in the case of recoil or from Eq. (A19) for the more general anomalous Doppler effect. For angles which satisfy this root we define $\theta_{cr}$. We also pick up a factor of $\omega n\beta$ in the denominator from the Jacobian. As such we obtain

[TABLE]

We note the prefactor is nothing more than the fine structure constant $\alpha=\frac{q^{2}}{4\pi}$. Moreover, we have now derived the Frank-Tamm formula generalized to the exact quantum result that includes the reaction of the particle to its recoil by the photon emission sokolov ; ido . Hence,

[TABLE]

To conclude this section, we have outlined the general procedure for the computation of the emission rate of Cherenkov radiation, including recoil, via the use of an Unruh-DeWitt detector. We outlined the general idea for how to incorporate recoil with the Unruh-DeWitt detector; this procedure is discussed in more detail in Section IV. What is most important is the fact that we are able to successfully incorporate recoil into the emission process. In the following section we will follow the same procedure while placing the current along an accelerated trajectory which is also superluminal. We then explore the emission spectrum of Cherenkov radiation, under acceleration, while simultaneously taking into account recoil and how these inherently quantum mechanical processes alter the emission spectrum.

V.2 Cherenkov emission via accelerated currents

Building upon the formalism of the previous section, let us now endeavor to compute the emission spectra by an accelerated electron current in the Cherenkov regime. The overall prescription will be to utilize the kinematics of Rindler space to determine trajectories, velocities, etc. and incorporate their dynamics into the emission. Recalling the photon emission probability, Eq. (A13)

[TABLE]

Note that we are back to integrating over the proper time of the electron and thus reabsorbed our $u_{0}$ components. To analyze accelerated Cherenkov emission we integrate over the proper time rather than laboratory time since the accelerated trajectories are parametrized this way. We will begin by examining the polarization vectors that are contracted with our velocity tensor. Recalling that under proper acceleration $\bar{a}$, the four-velocities at proper time $\tau$ will be given $u^{\mu}=(\cosh{(\bar{a}\tau)},0,0,\sinh{(\bar{a}\tau)})$ winitzki . Hence,

[TABLE]

Note, we have made explicit the magnitude squared of the sin term since our roots in the forthcoming delta functions will be complex in nature. From here we utilize the same change of variables to the difference and average proper times, $\xi^{\prime}$ and $\eta^{\prime}$. Moreover we will make use of the hyperbolic double angle formulas to obtain

[TABLE]

This analysis is simplified using the time-dependent formalism presented in lynch1 . Combining all the above pieces we can now formulate the emission rate, $\Gamma_{\eta^{\prime}}=\frac{d\mathcal{P}}{d\eta^{\prime}}$. Thus

[TABLE]

Note that with respect to the proper time $\eta^{\prime}$, the Lorentz gamma and velocity are given by $\gamma=\cosh{(\bar{a}\eta^{\prime})}$ and $\beta=\tanh{(\bar{a}\eta^{\prime})}$ respectively. The trajectory is determined by $\Delta\mathbf{x}_{tr}=\beta\Delta t$. The time difference is given by $\Delta t=\frac{2}{\bar{a}}\sinh{(\bar{a}\xi^{\prime}/2)}\gamma$ for uniform acceleration and we will make use of the short time approximation such that $\Delta t=\xi^{\prime}\gamma$. Now, to examine the momentum integrations we move to spherical coordinates with the momentum aligned along the $z$-axis. This will yield our emission spectrum $\frac{d\Gamma_{\eta^{\prime}}}{d\omega}$. Hence

[TABLE]

The integration over proper time can now be evaluated by the use of delta functions. Note, the first term in the square bracket is the standard Cherenkov term and the hyperbolic cosine encodes the acceleration into the emission. Conversion of the acceleration dependent term into exponents and integration yields,

[TABLE]

Converting the delta functions for the subsequent $\cos(\theta)$ integration yields a Jacobian of $n\beta\omega\gamma$ for each delta function as well for each of the following three roots;

[TABLE]

Recalling that $\sin{(\theta)}=\sqrt{1-\cos^{2}{(\theta)}}$. Evaluation of the angular integral then yields

[TABLE]

Recalling that the emission rate is parametrized by the electron’s proper time, we need to boost back into the lab frame via $\frac{d\Gamma_{\eta^{\prime}}}{d\omega}=\gamma\frac{d\Gamma}{\omega}$. We then have

[TABLE]

Care must be taken here since we have the proper acceleration and we need to boost back into the lab frame to have a well formulated emission rate. Recalling the relationship between the laboratory acceleration, $a$, and proper acceleration is given by $\bar{a}=a\gamma^{3}$ we have

[TABLE]

Note, we have successfully recovered the Frank-Tamm formula but now with an additional acceleration-dependent term. The above expressions shows the acceleration dependence in Cherenkov emission and, with the energy gap, can include the standard case, anomalous Doppler effect, and the recoil correction ginsburg ; ido . Here, and throughout the manuscript, we assume a constant index of refraction.