This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Accelerated-Cherenkov radiation and signatures of radiation reaction

Morgan H. Lynch morgan.lynch@technion.ac.il Department of Electrical Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel    Eliahu Cohen eliahu.cohen@biu.ac.il Physics Department, Centre for Research in Photonics, University of Ottawa, Advanced Research Complex, 25 Templeton, Ottawa, Ontario, Canada K1N 6N5 Faculty of Engineering and the Institute of Nanotechnology and Advanced Materials, Bar Ilan University, Ramat Gan 5290002, Israel    Yaron Hadad yaronhadad@gmail.com Department of Mathematics, University of Arizona, Tucson, AZ, 85721, USA    Ido Kaminer kaminer@technion.ac.il Department of Electrical Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel
(August 16, 2025)
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.

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 a¯\bar{a}, implied that Minkowski vacuum was comprised of a thermalized bath of particles at temperature T=a¯2π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 =c=kB=1\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, A^μ(x)\hat{A}^{\mu}(x), using the current interaction for QED peskin , S^I=d4xj^μ(x)A^μ(x)\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, j^μ(x)=uμq^(τ)δ3(xxtr)\hat{j}_{\mu}(x)=u_{\mu}\hat{q}(\tau)\delta^{3}(x-x_{tr}). Here we defined the Heisenberg evolved charge monopole moment operator q^=eiH^τq^(0)eiH^τ\hat{q}=e^{i\hat{H}\tau}\hat{q}(0)e^{-i\hat{H}\tau} where q^(0)\hat{q}(0) is defined as q^(0)|Ei=|Ef\hat{q}(0)\ket{E_{i}}=\ket{E_{f}} with EiE_{i} and EfE_{f} the initial energy and final energy of the electron moving along the trajectory, xtrx_{tr}, of the current with four velocity uμu_{\mu}. The energy gap of the detector is ΔE=EfEi\Delta E=E_{f}-E_{i} and the charge qq of the electron is given by q=|Ef|q^(0)|Ei|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; 𝒜=i𝐤|Ef|S^I|Ei|0.\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

Γ=q2𝑑ξeiΔEξ/γVμν[x,x]Gμν[x,x].\displaystyle\Gamma=q^{2}\int d\xi e^{-i\Delta E\xi/\gamma}V_{\mu\nu}[x^{\prime},x]G^{\mu\nu}[x^{\prime},x]. (1)

Here we are integrating over the laboratory time ξ\xi. The “velocity tensor” Vμν=vμvνV_{\mu\nu}=v_{\mu}v_{\nu} determines the trajectory of the electron with vμ=(1,0,0,β)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,

dΓdω=αβ[1cos2(θc)].\frac{d\Gamma}{d\omega}=\alpha\beta\left[1-\cos^{2}{(\theta_{c})}\right]. (2)

We note the prefactor is nothing more than the fine structure constant α=q24π\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(θc)=1nβ+ΔEnβωγ\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 Ei=m2+ω2E_{i}=\sqrt{m^{2}+\omega^{2}} and a final state electron energy comprised of a bare electron and recoil momentum k=nωk=n\omega (where n is the frequency-dependent index of refraction) produced by the photon emission, Ef=m2+(nω)2E_{f}=\sqrt{m^{2}+(n\omega)^{2}}, see Fig. 1 for an illustrative example. The recoil momentum, nωn\omega, is characteristic of the non-trivial dispersion relation of the photon. As such, we can determine the energy difference to be

ΔE=(n21)ω22m.\displaystyle\Delta E=\frac{(n^{2}-1)\omega^{2}}{2m}. (3)
Refer to caption
Figure 1: The conceptual foundation for an initial dressed electron subject to recoil after photon emission. This assumption is used to define the energy gap of the Unruh-DeWitt detector. The photon energy contributes to the electrons mass and the non trivial dispersion produced by the indexed medium produces the recoil momentum when emitted. The kinetic energy produced by this recoil defines our energy gap Eq. (3). This phenomenon can also be viewed as a consequence of mass renormalization in an optical medium.

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

cosθc=1nβ+(n21)ω2mγnβ.\cos{\theta_{c}}=\frac{1}{n\beta}+\frac{(n^{2}-1)\omega}{2m\gamma n\beta}. (4)

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,

dΓdω=αβ[11n2β2(1+(n21)ωmγ+(n21)2ω24m2γ2)].\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]. (5)

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

Γ=q2𝑑ξeiΔEξUμν[x,x]Gμν[x,x].\displaystyle\Gamma=q^{2}\int d\xi e^{-i\Delta E\xi}U_{\mu\nu}[x^{\prime},x]G^{\mu\nu}[x^{\prime},x]. (6)

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μν=uμuνU_{\mu\nu}=u_{\mu}u_{\nu}. Recalling that under proper acceleration a¯\bar{a}, the four-velocities at proper time τ\tau will be given uμ=(cosh(a¯τ),0,0,sinh(a¯τ))u^{\mu}=(\cosh{(\bar{a}\tau)},0,0,\sinh{(\bar{a}\tau)}) winitzki . Computation of the emission spectrum, parametrized by the lab frame acceleration aa, with a¯=aγ3\bar{a}=a\gamma^{3}, then yields the result,

dΓdω\displaystyle\frac{d\Gamma}{d\omega} =\displaystyle= α2βγ2[(2γ21)[1cos2(θc)][1cos2θc+(aγ2nβω)2]2+4cosθc(aγ2nβω)].\displaystyle\frac{\alpha}{2\beta\gamma^{2}}\left[\left(2\gamma^{2}-1\right)\left[1-\cos^{2}{(\theta_{c})}\right]-\sqrt{\left[1-\cos^{2}{\theta_{c}}+\left(\frac{a\gamma^{2}}{n\beta\omega}\right)^{2}\right]^{2}+4\cos{\theta_{c}}\left(\frac{a\gamma^{2}}{n\beta\omega}\right)}\right]. (7)

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θc=1nβ+ΔEnβωγ\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,

cos(θa)=cosθc2+(aγ2nβω)2.\cos{(\theta_{a})}=\sqrt{\cos{\theta_{c}}^{2}+\left(\frac{a\gamma^{2}}{n\beta\omega}\right)^{2}}. (8)

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.

[Uncaptioned image]
Refer to caption
Figure 2: The structure of Cherenkov spectra with and without acceleration and recoil. We assume a constant index n=1.5n=1.5, β=.7\beta=.7, and with an array of experimentally realizable accelerations. In (a) we see the overall large scale structure and in (b) we have zoomed in on the resonance peak. The acceleration scale 1.8×10181.8\times 10^{18} m/s2 corresponds to a laser field amplitude of 1 GeV/m, an acceleration energy scale of ωa=5.5×104\omega_{a}=5.5\times 10^{-4} eV, and a cutoff frequency, see Eq. (10), of ωc=1.7×103\omega_{c}=1.7\times 10^{-3} eV. Note that both acceleration and recoil are required to see the peak in the emission frequency.

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

ωc=aγ2nβ211cos2(θc).\displaystyle\omega_{c}=\frac{a\gamma}{\sqrt{2}n\beta^{2}}\frac{1}{\sqrt{1-\cos^{2}{(\theta_{c})}}}. (9)

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

ωc\displaystyle\omega_{c} =\displaystyle= ωa1(nβ)21.\displaystyle\omega_{a}\frac{1}{\sqrt{(n\beta)^{2}-1}}. (10)

Here we defined the characteristic energy scale set by the particle’s acceleration, ωa=aγ2β\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 ωa\omega_{a} in the Cherenkov regime. As with the Unruh effect unruh1 , the energy scale ωa\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 ωa=5.5×104\omega_{a}=5.5\times 10^{-4} eV used as the characteristic energy scale in our figures is based off an acceleration of, a1.8×1020a\sim 1.8\times 10^{20} m/s2, that is indeed experimentally realizable in common laser systems with an amplitude of 1\sim 1 GV/m. Moreover, table top plasma wakefield accelerators are even capable of producing accelerations, a10×1022a\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 ΔE=0\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 ΔE=ω2(n21)2m\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,

ωc\displaystyle\omega_{c} =\displaystyle= ωa(nβ)2(1+ωc(n21)2mγ)2\displaystyle\frac{\omega_{a}}{\sqrt{\left(n\beta\right)^{2}-\left(1+\frac{\omega_{c}(n^{2}-1)}{2m\gamma}\right)^{2}}} (11)

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,

ωc\displaystyle\omega_{c^{\prime}} =\displaystyle= 2mγ(nβ1)n21.\displaystyle\frac{2m\gamma(n\beta-1)}{n^{2}-1}. (12)

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, ωa2ω3=(n21)2γm[1+ω(n21)2mγ]\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,

ωp=[2mγωa2n21]1/3.\omega_{p}=\left[\frac{2m\gamma\omega_{a}^{2}}{n^{2}-1}\right]^{1/3}. (13)

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 ωa\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.5n=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.

[Uncaptioned image]
Refer to caption
Figure 3: The structure of Cherenkov spectral cutoffs and peak frequency emitted as a function of β\beta near the Cherenkov threshold. The shaded region determines the allowed frequencies of light emitted. The higher frequencies are bounded by the recoil correction while the lower frequencies are bounded by the acceleration energy scale ωa\omega_{a}. These bounds converge to the peak frequency at threshold. The asymptotes for the low and high frequency cutoffs are given by Eqs. (10) and (12), respectively. The peak frequency is given by Eq. (13) and the solid curve (both blue and red) of the spectral cutoff is given by Eq. (11). Here, we use the experimentally realizable parameters of ωa=5.5×104\omega_{a}=5.5\times 10^{-4} eV (a), ωa=11×104\omega_{a}=11\times 10^{-4} eV (b), and index of refraction n=1.5n=1.5.

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/β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. n1/β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 β0.60.7\beta\sim 0.6-0.7 and with laser field amplitudes currently in regular use 1\sim 1 GV/m, we should then see a signal between .11.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.

Refer to caption
Figure 4: The structure of the peak frequency as a function of the particle velocity and local acceleration. Here, the index of refraction is tuned to n1βn\sim\frac{1}{\beta} so we are always near resonance. Note that all velocity and acceleration profiles used are currently accessible experimentally and produce a peak frequency that is also well within current detection capabilities.

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θc=1nβ+ΔEnβωγ\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);

ωc\displaystyle\omega_{c} =\displaystyle= ΔEγ±(ΔEγ)2+((nβ)21)((ΔEγ)2+ωa2)n2β21.\displaystyle\frac{\frac{\Delta E}{\gamma}\pm\sqrt{\left(\frac{\Delta E}{\gamma}\right)^{2}+\left((n\beta)^{2}-1\right)\left(\left(\frac{\Delta E}{\gamma}\right)^{2}+\omega_{a}^{2}\right)}}{n^{2}\beta^{2}-1}. (14)

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 ωa\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(θc)=1nβ\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.

[Uncaptioned image]
Refer to caption
Figure 5: The frequency-angle structure of the Cherenkov cone as a function of acceleration for a constant index of n=1.5n=1.5 and with β=.7\beta=.7, see Eq. (8). Note the effect of acceleration will be to fill in the Cherenkov cone with the lower energy end of the emission spectrum and thus produce a “Cherenkov circle” rather than the standard Cherenkov ring.

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.

References

  • (1) W. G. Unruh, Phys. Rev. D 14, 870 (1976).
  • (2) B. S. DeWitt, General Relativity an Einstein centenary survey, S. W. Hawking, W. Israel, Cambridge University Press (Cambridge, 1979).
  • (3) J. Louko and A. Satz, Class. Quant. Grav. 25, 055012 (2008).
  • (4) D. E. Bruschi, et al., Phys. Rev. A 82, 042332 (2010).
  • (5) B. L. Hu, S.-Y. Lin, and J. Louko, Class. Quant. Grav. 29, 224005 (2012).
  • (6) D. E. Bruschi, A. R. Lee, and I. Fuentes, J. Phys. A: Math. Theor. 46, 165303 (2013).
  • (7) A. R. Lee and I. Fuentes, Phys. Rev. D 89, 085041 (2014).
  • (8) P. A. M. Dirac, Proc. Roy. Soc. Lond. A 167, 148 (1938).
  • (9) A. G. R. Thomas, et al., Phys. Rev. X 2, 041004 (2012).
  • (10) J. M. Cole, et al., Phys. Rev. X 8, 011020 (2018).
  • (11) K. Poder, et al., Phys. Rev. X 8, 031004 (2018).
  • (12) T. N. Wistisen, et al., Nat. Commun. 9, 795 (2018).
  • (13) Y. Y. Chen, et al., Phys. Rev. Lett. 121, 074801 (2018).
  • (14) A. Di Piazza, Lett. in Math. Phys. 83, 305 (2008).
  • (15) Y. Hadad, et al., Phys. Rev. D 82, 096012 (2010).
  • (16) D. A. Burton and A. Noble, Contemp. Phys. 55, 110 (2014).
  • (17) A. Ilderton and G. Torgrimsson, Phys. Rev. D 88, 025021 (2013).
  • (18) P. A. Cerenkov, Phys. Rev. 52, 378 (1937).
  • (19) J. D. Jackson, Classical Electrodynamics (Wiley, 1998).
  • (20) N. D. Birrell and P. C. W. Davies, Quantum Field Theory in Curved Space (Cambridge University Press, Cambridge, 1982).
  • (21) V. F. Mukhanov and S. Winitzki, Introduction to Quantum Effects in Gravity (Cambridge University Press, Cambridge, 2010).
  • (22) L. Parker, The Creation of Particles by the Expanding Universe, Ph.D. thesis, Harvard University, 1966; Nature (London) 261, 20 (1976).
  • (23) S. W. Hawking, Nature (London) 248, 30 (1974); Commun. Math. Phys. 43, 199 (1975).
  • (24) R. Muller, Phys. Rev. D 56, 953 (1997).
  • (25) G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 59, 094004 (1999).
  • (26) D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. D 63, 014010 (2000)
  • (27) D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. Lett. 87, 151301 (2001).
  • (28) M. H. Lynch, Phys. Rev. D 90, 024049 (2014).
  • (29) M. H. Lynch, Phys. Rev. D 92, 024019 (2015).
  • (30) I. M. Frank and I.E. Tamm, Compt. Rend. Acad. Sci. 14, (1937).
  • (31) I. Kaminer, et al., Phys. Rev. X 6, 011006 (2016).
  • (32) W. D. Kimura, et al., Phys. Rev. Lett. 74, 546 (1995).
  • (33) Y. Hu, et al., Sci. Rep. 7, 8695 (2017).
  • (34) M. Henstridge, et al., Science 26, 439 (2018).
  • (35) M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory (Westview Press, 1995).
  • (36) V. L. Ginzburg, Applications of Electrodynamics in Theoretical Physics and Astrophysics (CRC Press, 2017).
  • (37) X. Shi, et al., Nature Physics 14, 1001 (2018).
  • (38) V. N. Tsytovich, J. Exptl. Theoret. Phys. 42, 457 (1962).
  • (39) A. Sokolov, Dokl. Akad. Nauk SSSR 28, 415 (1940).
  • (40) J. B. Rosenzweig, et al., Nucl. Instrum. Meth. Phys. Res., Sect. A 653, 98 (2011).
  • (41) E. Martin-Martinez, M. Montero, and M. del Rey, Phys. Rev. D 87, 064038 (2013).
  • (42) F. J. Garcia de Abajo, Rev. Mod. Phys. 82, 209 (2010).
  • (43) V. Ginis, et al., Phys. Rev. Lett. 113, 167402 (2014).
  • (44) X. Lin et. al, Nat. Phys. 14, 816 (2018).
  • (45) V. L. Ginzburg and V. P. Frolov, Sov. Phys. Usp. 30, 1073 (1987).
  • (46) V. P. Frolov and V. L. Ginzburg, Phys. Lett. A 116, 423 (1986).
  • (47) L. C. B. Crispino, A. Higuchi, and G. E. A. Matsas, Rev. Mod. Phys. 80, 787 (2008).
  • (48) M. H. Lynch, et al., arXiv:1903.00043 [gr-qc], (2019).

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, A^μ(x)\hat{A}^{\mu}(x), using the current interaction for QED peskin ,

S^I=d4xj^μ(x)A^μ(x).\hat{S}_{I}=\int d^{4}x\hat{j}_{\mu}(x)\hat{A}^{\mu}(x). (A1)

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)u(s,p) and v(s,p)v(s,p) of spin ss and momentum pp that are created by a^s,p\hat{a}^{\dagger}_{s,p} and b^s,p\hat{b}^{\dagger}_{s,p}, we have the following electron field operators

ψ^(x,t)\displaystyle\hat{\psi}(x,t) =\displaystyle= d3ps[a^s,pu(s,p)ϕp(x,t)+b^s,pv(s,p)χp(x,t)],\displaystyle\int d^{3}p\sum_{s}\left[\hat{a}_{s,p}u(s,p)\phi_{p}(x,t)+\hat{b}^{\dagger}_{s,p}v(s,p)\chi_{p}(x,t)\right],
ψ¯^(x,t)\displaystyle\hat{\bar{\psi}}(x,t) =\displaystyle= d3ps[a^s,pu¯(s,p)ϕp(x,t)+b^s,pv¯(s,p)χp(x,t)].\displaystyle\int d^{3}p\sum_{s}\left[\hat{a}^{\dagger}_{s,p}\bar{u}(s,p)\phi^{\ast}_{p}(x,t)+\hat{b}_{s,p}\bar{v}(s,p)\chi^{\ast}_{p}(x,t)\right]. (A2)

The positive and negative frequency modes are given by ϕp(x,t)\phi_{p}(x,t) and χp(x,t)\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, j^μ=ψ¯^γμψ^\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μ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. 1\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

j^μ(x)=u¯(p)γμu(p)ϕEf(x,t)a^pa^pϕEi(x,t).\displaystyle\hat{j}_{\mu}(x)=\bar{u}(p)\gamma_{\mu}u(p)\phi^{\ast}_{E_{f}}(x,t)\hat{a}^{\dagger}_{p}\hat{a}_{p}\phi_{E_{i}}(x,t). (A3)

Now, in keeping with the weak recoil approximation we make use of the Gordon identity peskin . As such, we have u¯(p)γμu(p)=uμ\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 ϕ(x,τ)=g(x)eiEτ\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

j^μ(x)=uμg(x)g(x)eiEfτa^pa^peiEiτ.\displaystyle\hat{j}_{\mu}(x)=u_{\mu}g^{\ast}(x)g(x)e^{iE_{f}\tau}\hat{a}^{\dagger}_{p}\hat{a}_{p}e^{-iE_{i}\tau}. (A4)

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(x)g(x)=δ(xxtr)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ψf(x)ψi(x)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)=pgf(x)gi(x)eiΔEτF(x)=pg^{\ast}_{f}(x)g_{i}(x)e^{i\Delta E\tau}, where pp is the momentum of the electron and ΔE=EfEi\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

eiEfτa^pa^peiEiτ\displaystyle e^{iE_{f}\tau}\hat{a}^{\dagger}_{p}\hat{a}_{p}e^{-iE_{i}\tau} =\displaystyle= eiH^τq^(0)eiH^τ=q^(τ).\displaystyle e^{i\hat{H}\tau}\hat{q}(0)e^{-i\hat{H}\tau}=\hat{q}(\tau). (A5)

Here we defined the Heisenberg evolved charge monopole moment operator q^=eiH^τq^(0)eiH^τ\hat{q}=e^{i\hat{H}\tau}\hat{q}(0)e^{-i\hat{H}\tau} where q^(0)\hat{q}(0) is defined as q^(0)|Ei=|Ef\hat{q}(0)\ket{E_{i}}=\ket{E_{f}} with EiE_{i} and EfE_{f} the initial energy and final energy of the electron moving along the trajectory, xtrx_{tr}, of the current. The energy gap of the detector is ΔE=EfEi\Delta E=E_{f}-E_{i} and the charge qq of the electron is given by q=|Ef|q^(0)|Ei|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

j^μ(x)=uμq^(τ)δ3(xxtr).\hat{j}_{\mu}(x)=u_{\mu}\hat{q}(\tau)\delta^{3}(x-x_{tr}). (A6)

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;

𝒜=i𝐤|Ef|S^I|Ei|0.\mathcal{A}=i\bra{\mathbf{k}}\otimes\bra{E_{f}}\hat{S}_{I}\ket{E_{i}}\otimes\ket{0}. (A7)

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, d𝒫d3k=d3pfδ3(𝐩f𝐩r)|𝒜|2=|𝒜|2.\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 prp_{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 nn, we have

A^μ(x)\displaystyle\hat{A}^{\mu}(x) =\displaystyle= d3k(2π)3/2iϵiμn22ω[a^kei(𝐤𝐱ωt)+a^kei(𝐤𝐱ωt)]\displaystyle\int\frac{d^{3}k}{(2\pi)^{3/2}}\frac{\sum_{i}\epsilon_{i}^{\mu}}{\sqrt{n^{2}2\omega}}\left[\hat{a}_{k}e^{i(\mathbf{k}\cdot\mathbf{x}-\omega t)}+\hat{a}^{\dagger}_{k}e^{-i(\mathbf{k}\cdot\mathbf{x}-\omega t)}\right]

Here we have a dispersion relation which relates the frequency to momenta via nω=kn\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

𝒫\displaystyle\mathcal{P} =\displaystyle= d3k|𝐤|Ef|d4xj^μ(x)A^μ(x)|Ei|0|2\displaystyle\int d^{3}k\Big{|}\bra{\mathbf{k}}\otimes\bra{E_{f}}\int d^{4}x\hat{j}_{\mu}(x)\hat{A}^{\mu}(x)\ket{E_{i}}\otimes\ket{0}\Big{|}^{2} (A9)
=\displaystyle= d3kd4xd4x|Ef|j^μ(x)|Ei|2|𝐤|A^μ(x)|0|2.\displaystyle\int d^{3}k\int d^{4}x\int d^{4}x^{\prime}\Big{|}\bra{E_{f}}\hat{j}_{\mu}(x)\ket{E_{i}}|^{2}|\bra{\mathbf{k}}\hat{A}^{\mu}(x)\ket{0}\Big{|}^{2}.

Evaluation of the square of the electron current matrix element yields

|Ef|j^μ(x)|Ei|2\displaystyle\Big{|}\bra{E_{f}}\hat{j}_{\mu}(x)\ket{E_{i}}\Big{|}^{2} =\displaystyle= |Ef|uμ(x)eiH^τq^(0)eiH^τδ3(xx(t))|Ei|2\displaystyle\Big{|}\bra{E_{f}}u_{\mu}(x)e^{i\hat{H}\tau}\hat{q}(0)e^{-i\hat{H}\tau}\delta^{3}(x-x(t))\ket{E_{i}}\Big{|}^{2} (A10)
=\displaystyle= q2uμ(x)uν(x)δ3(xxtr(t))δ3(xxtr(t))eiΔE(ττ)\displaystyle q^{2}u_{\mu}(x)u_{\nu}(x^{\prime})\delta^{3}(x-x_{tr}(t))\delta^{3}(x^{\prime}-x^{\prime}_{tr}(t^{\prime}))e^{-i\Delta E(\tau^{\prime}-\tau)}

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,

d3k|𝐤|A^μ(x)|0|2=0|A^ν(x)A^μ(x)|0.\displaystyle\int d^{3}k|\bra{\mathbf{k}}\hat{A}^{\mu}(x)\ket{0}|^{2}=\bra{0}\hat{A}^{\dagger\nu}(x^{\prime})\hat{A}^{\mu}(x)\ket{0}. (A11)

Here we have factored out the completeness relation, 𝑑k|kk|=1\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

0|A^ν(x)A^μ(x)|0\displaystyle\bra{0}\hat{A}^{\dagger\nu}(x^{\prime})\hat{A}^{\mu}(x)\ket{0} =\displaystyle= 0|d3k(2π)3/2iϵiνn22ω[a^kei(𝐤𝐱ωt)+a^kei(𝐤𝐱ωt)]\displaystyle\bra{0}\int\frac{d^{3}k^{\prime}}{(2\pi)^{3/2}}\frac{\sum_{i}\epsilon_{i}^{{}^{\prime}\nu}}{\sqrt{n^{\prime 2}2\omega^{\prime}}}\left[\hat{a}_{k^{\prime}}e^{i(\mathbf{k^{\prime}}\cdot\mathbf{x^{\prime}}-\omega^{\prime}t^{\prime})}+\hat{a}^{\dagger}_{k^{\prime}}e^{-i(\mathbf{k^{\prime}}\cdot\mathbf{x^{\prime}}-\omega^{\prime}t^{\prime})}\right] (A12)
×\displaystyle\;\;\;\;\;\times d3k(2π)3/2iϵiμn22ω[a^kei(𝐤𝐱ωt)+a^kei(𝐤𝐱ωt)]|0\displaystyle\int\frac{d^{3}k}{(2\pi)^{3/2}}\frac{\sum_{i}\epsilon_{i}^{\mu}}{\sqrt{n^{2}2\omega}}\left[\hat{a}_{k}e^{i(\mathbf{k}\cdot\mathbf{x}-\omega t)}+\hat{a}^{\dagger}_{k}e^{-i(\mathbf{k}\cdot\mathbf{x}-\omega t)}\right]\ket{0}
=\displaystyle= 1(2π)312d3kn(ω)ωijϵiμϵjνei(𝐤Δ𝐱ω(tt)).\displaystyle\frac{1}{(2\pi)^{3}}\frac{1}{2}\int\frac{d^{3}k}{n(\omega)\omega}\sum_{ij}\epsilon_{i}^{\mu}\epsilon_{j}^{\nu^{\prime}}e^{i(\mathbf{k}\cdot\Delta\mathbf{x}-\omega(t^{\prime}-t))}.

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

𝒫\displaystyle\mathcal{P} =\displaystyle= d4xd4x|Ef|j^μ(x)|Ei|2|𝐤|A^μ(x)|0|2\displaystyle\int d^{4}x\int d^{4}x^{\prime}\Big{|}\bra{E_{f}}\hat{j}_{\mu}(x)\ket{E_{i}}\Big{|}^{2}\Big{|}\bra{\mathbf{k}}\hat{A}^{\mu}(x)\ket{0}\Big{|}^{2} (A13)
=\displaystyle= q21(2π)312d4xd4xd3kn(ω)ωuμ(x)uν(x)ijϵiμϵjν\displaystyle q^{2}\frac{1}{(2\pi)^{3}}\frac{1}{2}\int d^{4}x\int d^{4}x^{\prime}\int\frac{d^{3}k}{n(\omega)\omega}u_{\mu}(x)u_{\nu}(x^{\prime})\sum_{ij}\epsilon_{i}^{\mu}\epsilon_{j}^{\nu^{\prime}}
×\displaystyle\;\;\;\;\;\times δ3(xxtr(t))δ3(xxtr(t))eiΔE(ττ)+i(𝐤Δ𝐱ω(tt))\displaystyle\delta^{3}(x-x_{tr}(t))\delta^{3}(x^{\prime}-x_{tr}^{\prime}(t^{\prime}))e^{-i\Delta E(\tau^{\prime}-\tau)+i(\mathbf{k}\cdot\Delta\mathbf{x}-\omega(t^{\prime}-t))}
=\displaystyle= q21(2π)312𝑑t𝑑td3kn(ω)ωuμ(x)uν(x)ijϵiμϵjνu0(x)u0(x)ei(ΔE(ττ)𝐤Δ𝐱tr+ω(tt)).\displaystyle q^{2}\frac{1}{(2\pi)^{3}}\frac{1}{2}\int dtdt^{\prime}\int\frac{d^{3}k}{n(\omega)\omega}\frac{u_{\mu}(x)u_{\nu}(x^{\prime})\sum_{ij}\epsilon_{i}^{\mu}\epsilon_{j}^{\nu^{\prime}}}{u_{0}(x)u_{0}(x^{\prime})}e^{-i(\Delta E(\tau^{\prime}-\tau)-\mathbf{k}\cdot\Delta\mathbf{x}_{tr}+\omega(t^{\prime}-t))}.

Note, in the last line we changed our integration variable from proper time to coordinate time via dt=u0dτdt=u_{0}d\tau. Then, by defining the difference and average laboratory times, namely ξ=tt\xi=t^{\prime}-t and η=(t+t)/2\eta=(t^{\prime}+t)/2 respectively, we can decouple the temporal integrations lynch ; lynch1 . Note, this parameterization has a unit Jacobian, i.e. dtdt=dξdη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 Δx=x(t)x(t)=x(tt)=x(ξ)\Delta x=x(t^{\prime})-x(t)=x(t^{\prime}-t)=x(\xi). As a consequence we will have Δxtr=Δxtr(ξ)\Delta x_{tr}=\Delta x_{tr}(\xi). Moreover, with a trajectory that is constant in velocity we will have ϵiμ\epsilon_{i}^{\mu} and uμ(x)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 d𝒫dη=Γ\frac{d\mathcal{P}}{d\eta}=\Gamma. As such,

Γ=q21(2π)312𝑑ξd3kωuμ(x)uν(x)ijϵiμϵjνu0(x)u0(x)ei(ΔEγξ𝐤Δ𝐱tr(ξ)+ωξ).\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}\epsilon_{i}^{\mu}\epsilon_{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)}. (A14)

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=uμ(x)uν(x)i,jϵiμϵjνu0(x)u0(x)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 zz-axis aligned along the direction of the velocity, we have

Γ\displaystyle\Gamma =\displaystyle= q21(2π)212𝑑ξdcosθdωUωnei(ΔEγξωnΔxtr(ξ)cos(θ)+ωξ).\displaystyle q^{2}\frac{1}{(2\pi)^{2}}\frac{1}{2}\int d\xi\int d\cos{\theta}d\omega U\omega ne^{-i(\frac{\Delta E}{\gamma}\xi-\omega n\Delta x_{tr}(\xi)\cos{(\theta)}+\omega\xi)}. (A15)

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

Γ\displaystyle\Gamma =\displaystyle= q21(2π)212𝑑ξdcosθdωUωnei(ΔEγωnβcos(θ)+ω)ξ\displaystyle q^{2}\frac{1}{(2\pi)^{2}}\frac{1}{2}\int d\xi\int d\cos{\theta}d\omega U\omega ne^{-i(\frac{\Delta E}{\gamma}-\omega n\beta\cos{(\theta)}+\omega)\xi} (A16)
=\displaystyle= q214πdcosθdωUωn(ω)δ(ΔEγωnβcos(θ)+ω).\displaystyle q^{2}\frac{1}{4\pi}\int d\cos{\theta}d\omega U\omega n(\omega)\delta(\frac{\Delta E}{\gamma}-\omega n\beta\cos{(\theta)}+\omega).

Let us examine the scalar quantity UU 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

U\displaystyle U =\displaystyle= (iuμu0ϵiμ)(juνu0ϵjν)=β2sin2(θ).\displaystyle\left(\sum_{i}\frac{u_{\mu}}{u_{0}}\epsilon^{\mu}_{i}\right)\left(\sum_{j}\frac{u_{\nu}}{u_{0}}\epsilon^{\nu}_{j}\right)^{\dagger}=\beta^{2}\sin^{2}{(\theta)}. (A17)

Here we see the sin(θ)\sin{(\theta)} rather than cos(θ)\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(π/2θ)=sin(θ)\cos{(\pi/2-\theta)}=\sin{(\theta)} from the dot product. As such, our emission rate now becomes

Γ\displaystyle\Gamma =\displaystyle= q2β24πdcosθdωsin2(θ)ωnδ(ΔEγωnβcos(θ)+ω).\displaystyle q^{2}\frac{\beta^{2}}{4\pi}\int d\cos{\theta}d\omega\sin^{2}{(\theta)}\omega n\delta(\frac{\Delta E}{\gamma}-\omega n\beta\cos{(\theta)}+\omega). (A18)

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,

cosθ=1nβ+ΔEnβωγ.\cos{\theta}=\frac{1}{n\beta}+\frac{\Delta E}{n\beta\omega\gamma}. (A19)

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 Ei=m2+ω2E_{i}=\sqrt{m^{2}+\omega^{2}} and a final state electron energy comprised of a bare electron and recoil momentum k=nωk=n\omega (where n is the frequency-dependent index of refraction) produced by the photon emission, Ef=m2+(nω)2E_{f}=\sqrt{m^{2}+(n\omega)^{2}}, see Fig. 1 for an illustrative example. The recoil momentum, nωn\omega, is characteristic of the non-trivial dispersion relation of the photon. As such, we can determine the energy difference to be

ΔE(n21)ω22m.\displaystyle\Delta E\approx\frac{(n^{2}-1)\omega^{2}}{2m}. (A20)

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

cosθ=1nβ+(n21)ω2mγnβ.\cos{\theta}=\frac{1}{n\beta}+\frac{(n^{2}-1)\omega}{2m\gamma n\beta}. (A21)

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 θcr\theta_{cr}. We also pick up a factor of ωnβ\omega n\beta in the denominator from the Jacobian. As such we obtain

Γ\displaystyle\Gamma =\displaystyle= q2β24πdcosθdωsin2(θ)ωnδ(cos(θ)cos(θcr))ωnβ\displaystyle q^{2}\frac{\beta^{2}}{4\pi}\int d\cos{\theta}d\omega\sin^{2}{(\theta)}\omega n\frac{\delta(\cos{(\theta)}-\cos{(\theta_{cr})})}{\omega n\beta} (A22)
=\displaystyle= q2β4π𝑑ω(1cos2(θcr)).\displaystyle q^{2}\frac{\beta}{4\pi}\int d\omega\left(1-\cos^{2}{(\theta_{cr})}\right).

We note the prefactor is nothing more than the fine structure constant α=q24π\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,

dΓdω=αβ[11n2β2(1+(n21)ωmγ+(n21)2ω24m2γ2)].\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]. (A23)

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)

𝒫=q21(2π)312𝑑τ𝑑τd3kωσμσνUμνei(ΔEΔτ𝐤Δ𝐱tr+ωΔt).\mathcal{P}=q^{2}\frac{1}{(2\pi)^{3}}\frac{1}{2}\int d\tau d\tau^{\prime}\int\frac{d^{3}k}{\omega}\sigma^{\mu}\sigma^{\dagger\nu}U_{\mu\nu}e^{-i(\Delta E\Delta\tau-\mathbf{k}\cdot\Delta\mathbf{x}_{tr}+\omega\Delta t)}. (A24)

Note that we are back to integrating over the proper time of the electron and thus reabsorbed our u0u_{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 a¯\bar{a}, the four-velocities at proper time τ\tau will be given uμ=(cosh(a¯τ),0,0,sinh(a¯τ))u^{\mu}=(\cosh{(\bar{a}\tau)},0,0,\sinh{(\bar{a}\tau)}) winitzki . Hence,

U\displaystyle U =\displaystyle= (iuμϵiμ)(juνϵjν)=sinh(a¯τ)sinh(a¯τ)|sin(θ)|2.\displaystyle\left(\sum_{i}u_{\mu}\epsilon^{\mu}_{i}\right)\left(\sum_{j}u_{\nu}\epsilon^{\nu}_{j}\right)^{\dagger}=\sinh{(\bar{a}\tau)}\sinh{(\bar{a}\tau^{\prime})}|\sin{(\theta)}|^{2}. (A25)

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

sinh(a¯τ)sinh(a¯τ)=12[2cosh2(a¯η)1cosh(a¯ξ)].\sinh{(\bar{a}\tau)}\sinh{(\bar{a}\tau^{\prime})}=\frac{1}{2}\left[2\cosh^{2}{(\bar{a}\eta^{\prime})}-1-\cosh{(\bar{a}\xi^{\prime})}\right]. (A26)

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

Γη\displaystyle\Gamma_{\eta^{\prime}} =\displaystyle= q21(2π)314𝑑ξd3kn2(ω)ω[2cosh2(a¯η)1cosh(a¯ξ)]|sin(θ)|2ei(ΔEξ𝐤Δ𝐱tr+ωΔt).\displaystyle q^{2}\frac{1}{(2\pi)^{3}}\frac{1}{4}\int d\xi^{\prime}\int\frac{d^{3}k}{n^{2}(\omega)\omega}\left[2\cosh^{2}{(\bar{a}\eta^{\prime})}-1-\cosh{(\bar{a}\xi^{\prime})}\right]|\sin{(\theta)}|^{2}e^{-i(\Delta E\xi^{\prime}-\mathbf{k}\cdot\Delta\mathbf{x}_{tr}+\omega\Delta t)}. (A27)

Note that with respect to the proper time η\eta^{\prime}, the Lorentz gamma and velocity are given by γ=cosh(a¯η)\gamma=\cosh{(\bar{a}\eta^{\prime})} and β=tanh(a¯η)\beta=\tanh{(\bar{a}\eta^{\prime})} respectively. The trajectory is determined by Δ𝐱tr=βΔt\Delta\mathbf{x}_{tr}=\beta\Delta t. The time difference is given by Δt=2a¯sinh(a¯ξ/2)γ\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 Δt=ξγ\Delta t=\xi^{\prime}\gamma. Now, to examine the momentum integrations we move to spherical coordinates with the momentum aligned along the zz-axis. This will yield our emission spectrum dΓηdω\frac{d\Gamma_{\eta^{\prime}}}{d\omega}. Hence

dΓηdω\displaystyle\frac{d\Gamma_{\eta^{\prime}}}{d\omega} =\displaystyle= q21(2π)214𝑑ξdcos(θ)nω[2γ21cosh(a¯ξ)]|sin(θ)|2ei(ΔEnωcos(θ)βγ+ωγ)ξ.\displaystyle q^{2}\frac{1}{(2\pi)^{2}}\frac{1}{4}\int d\xi^{\prime}d\cos{(\theta)}n\omega\left[2\gamma^{2}-1-\cosh{(\bar{a}\xi^{\prime})}\right]|\sin{(\theta)}|^{2}e^{-i(\Delta E-n\omega\cos{(\theta)}\beta\gamma+\omega\gamma)\xi^{\prime}}. (A28)

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,

dΓηdω\displaystyle\frac{d\Gamma_{\eta^{\prime}}}{d\omega} =\displaystyle= q21(2π)214𝑑ξdcos(θ)nω[2γ2112ea¯ξ12ea¯ξ]|sin(θ)|2ei(ΔEnωcos(θ)βγ+ωγ)ξ\displaystyle q^{2}\frac{1}{(2\pi)^{2}}\frac{1}{4}\int d\xi^{\prime}d\cos{(\theta)}n\omega\left[2\gamma^{2}-1-\frac{1}{2}e^{\bar{a}\xi^{\prime}}-\frac{1}{2}e^{-\bar{a}\xi^{\prime}}\right]|\sin{(\theta)}|^{2}e^{-i(\Delta E-n\omega\cos{(\theta)}\beta\gamma+\omega\gamma)\xi^{\prime}} (A29)
=\displaystyle= q21(2π)14dcos(θ)|sin(θ)|2nω[(2γ21)δ(ΔEnωcos(θ)βγ+ωγ)\displaystyle q^{2}\frac{1}{(2\pi)}\frac{1}{4}\int d\cos{(\theta)}|\sin{(\theta)}|^{2}n\omega\Big{[}\left(2\gamma^{2}-1\right)\delta{(\Delta E-n\omega\cos{(\theta)}\beta\gamma+\omega\gamma)}
12δ(ΔEnωcos(θ)βγ+ωγ+ia¯)12δ(ΔEnωcos(θ)βγ+ωγia¯)].\displaystyle\;\;\;\;-\frac{1}{2}\delta{(\Delta E-n\omega\cos{(\theta)}\beta\gamma+\omega\gamma+i\bar{a})}-\frac{1}{2}\delta{(\Delta E-n\omega\cos{(\theta)}\beta\gamma+\omega\gamma-i\bar{a})}\Big{]}.

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

cosθc\displaystyle\cos{\theta_{c}} =\displaystyle= 1nβ+ΔEnβωγ\displaystyle\frac{1}{n\beta}+\frac{\Delta E}{n\beta\omega\gamma}
cosθ±a\displaystyle\cos{\theta_{\pm a}} =\displaystyle= cosθc±ia¯nβωγ.\displaystyle\cos{\theta_{c}}\pm\frac{i\bar{a}}{n\beta\omega\gamma}. (A30)

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

dΓηdω\displaystyle\frac{d\Gamma_{\eta}}{d\omega} =\displaystyle= α2βγ[(2γ21)[1cos2(θc)]1(cosθc+ia¯nβωγ)21(cosθcia¯nβωγ)2]\displaystyle\frac{\alpha}{2\beta\gamma}\left[\left(2\gamma^{2}-1\right)\left[1-\cos^{2}{(\theta_{c})}\right]-\sqrt{1-\left(\cos{\theta_{c}}+\frac{i\bar{a}}{n\beta\omega\gamma}\right)^{2}}\sqrt{1-\left(\cos{\theta_{c}}-\frac{i\bar{a}}{n\beta\omega\gamma}\right)^{2}}\right] (A31)
=\displaystyle= α2βγ[(2γ21)[1cos2(θc)][1cos2θc+(a¯nβωγ)2]2+4cosθc(a¯nβωγ)].\displaystyle\frac{\alpha}{2\beta\gamma}\left[\left(2\gamma^{2}-1\right)\left[1-\cos^{2}{(\theta_{c})}\right]-\sqrt{\left[1-\cos^{2}{\theta_{c}}+\left(\frac{\bar{a}}{n\beta\omega\gamma}\right)^{2}\right]^{2}+4\cos{\theta_{c}}\left(\frac{\bar{a}}{n\beta\omega\gamma}\right)}\right].

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

dΓdω\displaystyle\frac{d\Gamma}{d\omega} =\displaystyle= α2βγ2[(2γ21)[1cos2(θc)][1cos2θc+(a¯nβωγ)2]2+4cosθc(a¯nβωγ)].\displaystyle\frac{\alpha}{2\beta\gamma^{2}}\left[\left(2\gamma^{2}-1\right)\left[1-\cos^{2}{(\theta_{c})}\right]-\sqrt{\left[1-\cos^{2}{\theta_{c}}+\left(\frac{\bar{a}}{n\beta\omega\gamma}\right)^{2}\right]^{2}+4\cos{\theta_{c}}\left(\frac{\bar{a}}{n\beta\omega\gamma}\right)}\right]. (A32)

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, aa, and proper acceleration is given by a¯=aγ3\bar{a}=a\gamma^{3} we have

dΓdω\displaystyle\frac{d\Gamma}{d\omega} =\displaystyle= α2βγ2[(2γ21)[1cos2(θc)][1cos2θc+(aγ2nβω)2]2+4cosθc(aγ2nβω)].\displaystyle\frac{\alpha}{2\beta\gamma^{2}}\left[\left(2\gamma^{2}-1\right)\left[1-\cos^{2}{(\theta_{c})}\right]-\sqrt{\left[1-\cos^{2}{\theta_{c}}+\left(\frac{a\gamma^{2}}{n\beta\omega}\right)^{2}\right]^{2}+4\cos{\theta_{c}}\left(\frac{a\gamma^{2}}{n\beta\omega}\right)}\right]. (A33)

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.