Neutrino oscillations in the Unruh radiation
Massimo Blasone, Gaetano Lambiase, Giuseppe Gaetano Luciano and, Luciano Petruzziello

TL;DR
This paper demonstrates that neutrino oscillations occur in the Unruh thermal bath and that flavor neutrinos are essential as asymptotic states, extending the understanding of quantum field theory in accelerated frames.
Contribution
It extends previous work on Unruh radiation by incorporating neutrino mixing and shows that neutrino oscillations naturally occur in the accelerated frame's thermal bath.
Findings
Neutrino oscillations occur in the Unruh thermal bath.
Flavor neutrinos are necessary as asymptotic states.
The analysis confirms the consistency of neutrino mixing with the Unruh effect.
Abstract
The study of the decay of an accelerated proton recently provided a theoretical proof of the Unruh effect. On the basis of general covariance of Quantum Field Theory, indeed, it was found that the decay rates in the inertial and comoving frames do coincide only when the thermal nature of the accelerated vacuum is taken into account. Such an analysis was then extended to the case with neutrino mixing. In this Letter, we show that, by further embedding neutrino oscillations in the above framework, general covariance necessarily entails the use of flavor neutrinos as asymptotic states, as well as the occurrence of neutrino oscillations in the Unruh thermal bath.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Neutrino oscillations in the Unruh radiation
Massimo [email protected], Gaetano [email protected], Giuseppe Gaetano [email protected] and Luciano [email protected]
1Dipartimento di Fisica, Università di Salerno, Via Giovanni Paolo II, 132 I-84084 Fisciano (SA), Italy.
2INFN, Sezione di Napoli, Gruppo collegato di Salerno, Italy.
Abstract
The study of the decay of an accelerated proton recently provided a “theoretical proof” of the Unruh effect. On the basis of general covariance of Quantum Field Theory, indeed, it was found that the decay rates in the inertial and comoving frames do coincide only when the thermal nature of the accelerated vacuum is taken into account. Such an analysis was then extended to the case with neutrino mixing. In this Letter, we show that, by further embedding neutrino oscillations in the above framework, general covariance necessarily entails the use of flavor neutrinos as asymptotic states, as well as the occurrence of neutrino oscillations in the Unruh thermal bath.
pacs:
13.30.–a, 04.62.+v, 14.20.Dh, 95.30.Cq, 14.60.Pq
Introduction. In the history of Physics, the adoption of principles has revealed to be a formidable investigation tool. Although intimately related to the phenomenological realm from which they stem, once elevated to the status of postulates, physical principles act as lighthouses for the development of a consistent theoretical apparatus. Paradigmatic examples are the principle of conservation of energy, which led for instance to the discovery of the neutrino, and the principle of constancy of speed of light, at the basis of the construction of Special Relativity.
Recently, general covariance has been advocated to exhibit that the internal consistency of quantum field theory (QFT) unavoidably requires the existence of the Unruh effect, also known as Fulling-Davies-Unruh effect Fulling ; Davies ; Unruh:1976db . Indeed, starting from the statement that acceleration can influence even the proper lifetime of stable particles Muller:1997rt , in a series of remarkable papers Matsas:1999jx it was shown that the tree-level decay rate of an accelerated proton via inverse -decay is frame-independent only when the thermal nature of the vacuum for a non-inertial observer is considered.
The fact that a theoretical requirement leads to the specific form of the ground state for an accelerated observer should be regarded as a considerable result, especially in view of the perplexities which have been sometimes raised about the physical significance of the Unruh effect scettici . Such a skepticism is enhanced by the lack of direct evidences of this phenomenon, as it also happens for the case of the Hawking radiation Hawking . In fact, at present the most likely arena for (indirect) experimental tests of these effects is given by analogue gravity Barcelo .
In the above cited studies on the accelerated proton decay, the emitted neutrino was initially treated as massless and only in Ref. suzuki as a massive particle. In these works, however, neutrino mixing was not taken into account. This was done for the first time in Ref. Aluw , where a discrepancy between the proton decay rate in the inertial and comoving frames was claimed to arise. Subsequently, it was proved that general covariance does indeed hold in the above analysis: this has been achieved by employing either flavor Blasone:2018czm or mass Cozzella:2018qew eigenstates for neutrinos, thus leaving an essential ambiguity on the very nature of asymptotic neutrino states.
In this Letter we show that, due to the occurrence of neutrino oscillations in the problem at hand, general covariance leads to the conclusion that flavor eigenstates are fundamental and neutrinos belonging to the Unruh radiation must oscillate. Throughout the paper, we use natural units and the Minkowski metric with the mostly negative signature.
Let us start by setting the framework. Following the approach and the notation of Refs. Matsas:1999jx , the proton and the neutron can be viewed as unexcited and excited states of the nucleon, respectively. In addition, we assume to deal with particles that are energetic enough to possess a well-defined trajectory. In these conditions, it is possible to employ the Fermi theory of current-current interaction, where we consider a quantum leptonic and a classical hadronic current , with . Here555We assume the acceleration to occur along the -direction. With this choice, the Rindler coordinates are related with the Minkowski coordinates as follows: , , with and untouched., is the nucleon’s proper time (with being the Rindler time coordinate), its proper acceleration and represents the spatial Rindler coordinate which denotes the world line of the uniformly accelerated particle. The four-velocity of the nucleon is given by , in Rindler and Minkowski coordinates, respectively. In accordance with Refs. Matsas:1999jx ; Birrell , the Hermitian monopole is given by where is the nucleon Hamiltonian and is used to reconstruct the effective Fermi constant
Assuming to deal with a simplified two-flavor model, the interaction of the charged leptons and neutrinos with the nucleon current is described by the Fermi action
[TABLE]
where and are the gamma matrices in Dirac representation (e.g., see Ref. Itzykson ). In Eq. (1), neutrino fields with definite flavors are related to the ones with definite masses by
[TABLE]
In what follows, however, we shall focus on the processes involving only the production of a positron. The case in which an anti-muon appears in the final state can be treated analogously and in an independent way.
Inertial frame. In the inertial frame, the process we consider is the decay of an accelerated proton via inverse -decay
[TABLE]
which is pictorially represented in Fig. 1.
The fermion field is expanded as
[TABLE]
where is the polarization, and is the spinor defined as in Ref. Blasone:2018czm .
The tree-level transition amplitude for the process (i) reads Blasone:2018czm
[TABLE]
where for simplicity we have omitted the k- and -dependence of the lepton states, we have assumed equal momenta and polarizations for neutrinos with definite masses and
[TABLE]
with .
Note that the asymptotic flavor neutrino state in Eq. (Neutrino oscillations in the Unruh radiation) has been expressed in terms of the corresponding mass states by means of Pontecorvo mixing transformations Bilenky:1978nj
[TABLE]
The differential transition probability is defined as
[TABLE]
and the transition rate is given by , with
[TABLE]
being the nucleon total proper time.
In Ref. Blasone:2018czm , it has been shown that the transition rate for the process Eq. (8) is
[TABLE]
Here, we have introduced the shorthand notation
[TABLE]
where and
[TABLE]
The aim of the calculations contained in Ref. Blasone:2018czm is to exhibit the equality between Eq. (10) and its counterpart in the accelerated frame (see below), which guarantees the validity of the principle of general covariance.
At this point, it must be emphasized that in the above calculation an infinite proper time interval is considered, which allows for the electron neutrino produced in the proton decay to oscillate. Thus, we should take into account not only the process contemplated in Eq. (3), but also the following one:
[TABLE]
The above relation must be intended in the sense of Fig. 1: although it is true that the lepton charge must necessarily be conserved in the vertex (at tree-level), as soon as the outgoing neutrino is produced, flavor oscillations will inevitably take place. We remark that this process has not been included in the analysis of Ref. Blasone:2018czm , without however affecting the validity of the results there contained.
The transition amplitude for the process (ii) is indeed non-vanishing and is given by
[TABLE]
In terms of , the quantity of Eq. (Neutrino oscillations in the Unruh radiation) associated to the process (ii) leads to the following transition rate:
[TABLE]
where all the contributions in the r.h.s. of the above expression have already been introduced in Eqs. (11) and (Neutrino oscillations in the Unruh radiation). We notice that the above transition rate is proportional to , thus showing that it is originated by interference.
Finally, observe that
[TABLE]
Accelerated frame. From the point of view of an observer comoving with the proton, the only possible way to render the particle’s decay viable is to suppose the existence of a thermal bath of electrons, neutrinos and the corresponding antiparticles Blasone:2018czm . In such conditions, recalling that here we consider the electron decay channel only, three processes can be taken into account to match the inertial process Eq. (3):
[TABLE]
In order to calculate the proton’s decay rate in the accelerated frame, we need the expansion for the fermion fields in Rindler coordinates, that is
[TABLE]
where , is the Rindler frequency which does not satisfy any dispersion relation, and is defined as in Ref. Blasone:2018czm .
By way of illustration, let us consider the scattering (iii) of Eq. (17a) (see also Fig. 2); similar calculations can be carried out for the processes in Eqs. (17b), (17c). We obtain Blasone:2018czm
[TABLE]
where we have assumed equal frequencies and polarizations for different neutrinos with definite masses, whereas the quantity denoted with is equal to
[TABLE]
with . The spinor component related to the neutrino field contains the information on the mass of , and by means of the current hypothesis this is the only difference between the functions with different indexes.
The sum of the transition rates for the three processes in Eqs. (17) yields Blasone:2018czm
[TABLE]
where
[TABLE]
and . The functions and in Eqs. (21) and (22), respectively, are defined by
[TABLE]
and
[TABLE]
where and we have used the shorthand notation .
In Ref. Blasone:2018czm it was shown that, by taking into account the three interactions in Eq. (17) of the proton with the particles of the thermal bath, the quantity matches the corresponding decay rate evaluated in the inertial frame. In particular, it is possible to exhibit Blasone:2018czm that for , whereas and are equal to each other only up to a first-order expansion in the parameter . In such an approximation, Pontecorvo states Eq. (7) can be identified with the exact flavor neutrino states, defined as eigenstates of flavor charges Blasone:1995zc ; Blasone:2001qa .
On the other hand, we have seen above that an additional contribution to the proton decay rate has to be considered, which in the inertial frame is represented by the process (ii) in Fig. 1. Guided by the principle of general covariance, we now seek the corresponding processes in the comoving frame which should lead to the same result.
To this aim, we consider the three following contributions as potential candidates for the non-inertial counterpart of the decay Eq. (13):
[TABLE]
which are depicted in Fig. 2 in the last three diagrams. Note that, while the process (vi) is of the same type of (ii) since it entails an oscillation of the emitted (electron) neutrino, the processes (vii) and (viii) are essentially due to the oscillation of an (muon) antineutrino that is already present in the Unruh thermal bath.
In order to legitimate the validity of our assumption, we need to perform the same calculations that lead to the decay rate of Eq. (21). The outcome of this procedure turns out to be
[TABLE]
By virtue of the aforesaid observations contained in detail in Ref. Blasone:2018czm which allow us to state that , it is possible to infer that such an equivalence holds also between the decay rates of Eqs. (15) and (26).
Moreover, if we compute the total comoving decay rate which includes neutrino oscillations, we deduce that
[TABLE]
By comparing this with the total inertial decay rate of Eq. (16), we find that
[TABLE]
which means that such a result does not depend on the quantities and , whose treatment would require additional computational effort Blasone:2018czm .
Remarkably, Eq. (28) not only involves a generalization of the analysis of the accelerated proton decay to the case in which the produced neutrino oscillates, but it also unambiguously corroborates our guess of selecting the processes in Eqs. (25) as the counterpart for the decay (ii) in the inertial frame. Hence, the requirement of the principle of general covariance clearly results in the necessity of having an Unruh thermal bath containing flavor neutrinos which do oscillate.
Conclusions. In this Letter, we have extended the formalism firstly introduced in Refs. Matsas:1999jx and then developed in Refs. Aluw ; Blasone:2018czm ; Cozzella:2018qew to the case in which neutrino oscillations are taken into account.
In these works, the principle of general covariance has been exploited in the analysis of accelerated proton decay to probe the internal consistency of the theory and to shed light on new aspects of this problem. In such a perspective, here we have further studied this subject by considering flavor oscillations for the emitted neutrinos. By enforcing general covariance, we have then found that the Unruh radiation “seen” by the accelerated proton must involve oscillating neutrinos. This is a novel feature which emerges in a natural way in the present approach.
A further interesting observation that can be deduced from our analysis is related to the identities Eqs. (16) and (27) that are true for the inertial and the comoving frame, respectively. For this purpose, we recall that the decay rates appearing in the aforementioned equations have been computed by employing neutrino flavor states as asymptotic states. However, we note that similar relations also hold for the quantities and calculated in Ref. Cozzella:2018qew using neutrino mass eigenstates as fundamental objects. We then have
[TABLE]
where the equality has to be regarded both in the inertial and the comoving frames. Such an equation constitutes a consistency check for the correctness of the calculations in Refs. Blasone:2018czm and Cozzella:2018qew . The physical meaning of Eq. (29) can be understood by considering the charges for mixed neutrino fields as derived from Noether’s theorem Blasone:2001qa . Indeed, by denoting with
[TABLE]
the conserved charges for the neutrino fields with definite masses and with
[TABLE]
the (time-dependent) flavor charges, one can see that Blasone:2001qa
[TABLE]
where represents the total charge. The above relation can be interpreted as the conservation of the total lepton number. On the one hand, this can be viewed as the sum of two separately conserved family lepton numbers, when no mixing is present; on the other hand, the same conserved number is obtained by the sum of non-conserved flavor charges, which indeed are associated to neutrino oscillations.
Apart from its relevance in the context of neutrino mixing and oscillations, we stress that the Unruh effect is an excellent benchmark for both testing well-established predictions and pointing out novel effects in fundamental physics, as it combines such wide domains as general relativity, quantum mechanics and thermodynamics. For instance, in Refs. nonthermal it has been shown that the Unruh spectrum may exhibit exotic non-thermal corrections even within the standard QFT, thus emphasizing how such a framework provides an active forge of still unexplored scenarios. A similar non-thermal behavior has been obtained in Refs. Luciano:2019mrz , where Planck scale effects on the Unruh bath have been derived in the context of the Generalized Uncertainty Principle gup . Further features of the Unruh effect may be addressed in connection with entanglement properties for accelerated observers, whose implications have been investigated also in the context of black hole physics Hwang:2010ib . In particular, in Ref. FuentesSchuller:2004xp , it has been proved that entanglement turns out to be an observer-dependent quantity in non-inertial frames, due to the Unruh radiation. The question thus arises as to how this setting is modified in the presence of mixed neutrinos, particularly in view of the discussion of Refs. blasent .
As a final remark, we stress that in the current Letter we have made use of the simplest framework of neutrino mixing among two generations. The extension to three flavors is in principle straightforward, and represents one of the future directions of our investigation. We envisage that the presence of CP violation may introduce interesting additional features which would enrich the non-trivial structure of the Unruh radiation.
References
(2)
S. A. Fulling, Phys. Rev. D 7, 2850 (1973).
(3)
P. C. W. Davies, J. Phys. A 8, 609 (1975).
(4)
W. G. Unruh, Phys. Rev. D 14, 870 (1976).
(5)
R. Muller, Phys. Rev. D 56, 953 (1997).
(6)
G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 59, 094004 (1999); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. D 63, 014010 (2000); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. Lett. 87, 151301 (2001).
(7)
V. A. Belinskii, B. M. Karnakov, V. D. Mur and N. B. Narozhny, Pis’ma Zh. Eksp. Teor. Fiz. 65, 861 [JETP Lett. 65, 902 (1997)]; A. M. Fedotov, V. D. Mur, N. B. Narozhny, V. A. Belinskii and B. M. Karnakov, Phys. Lett. A 254, 126 (1999); N. B. Narozhny, A. M. Fedotov, B. M. Karnakov, V. D. Mur and V. A. Belinskii, Phys. Rev. D 65, 025004 (2001); N. B. Narozhny, A. M. Fedotov, B. M. Karnakov, V. D. Mur and V. A. Belinskii, Phys. Rev. D 70, 048702 (2004).
(8)
S. W. Hawking, Nature 248, 30 (1974); S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) Erratum: [Commun. Math. Phys. 46, 206 (1976)].
(9)
C. Barcelo, S. Liberati and M. Visser, Living Rev. Rel. 14, 3 (2011); J. Drori, Y. Rosenberg, D. Bermudez, Y. Silberberg and U. Leonhardt, Phys. Rev. Lett. 122, 010404 (2019).
(10)
H. Suzuki and K. Yamada, Phys. Rev. D 67, 065002 (2003).
(11)
D. V. Ahluwalia, L. Labun and G. Torrieri, Eur. Phys. J. A 52, 189 (2016).
(12)
M. Blasone, G. Lambiase, G. G. Luciano and L. Petruzziello, Phys. Rev. D 97, 105008 (2018).
(13)
G. Cozzella, S. A. Fulling, A. G. S. Landulfo, G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 97, 105022 (2018).
(14)
N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space, (Cambridge University Press, Cambridge, 1982).
(15)
C. Itzykson and J. B. Zuber, Quantum Field Theory (McGrawHill, New York, 1980).
(16)
S. M. Bilenky and B. Pontecorvo, Phys. Rept. 41, 225 (1978).
(17)
M. Blasone and G. Vitiello, Annals Phys. 244 283 (1995).
(18)
M. Blasone, P. Jizba and G. Vitiello, Phys. Lett. B 517 471 (2001).
(19)
R. Carballo-Rubio, L. J. Garay, E. Martín-Martínez and J. De Ramón, arXiv:1804.00685 [quant-ph]; J. Marino, A. Noto and R. Passante, Phys. Rev. Lett. 113, no. 2, 020403 (2014); M. Blasone, G. Lambiase and G. G. Luciano, Phys. Rev. D 96, 025023 (2017); M. Blasone, G. Lambiase and G. G. Luciano, J. Phys. Conf. Ser. 956, 012021 (2018).
(20)
F. Scardigli, M. Blasone, G. Luciano and R. Casadio, Eur. Phys. J. C 78, 728 (2018); G. G. Luciano and L. Petruzziello, arXiv:1902.07059 [hep-th].
(21)
H. S. Snyder, Phys. Rev. 71, 38 (1947); D. Amati, M. Ciafaloni and G. Veneziano, Phys. Lett. B 197, 81 (1987); M. Maggiore, Phys. Lett. B 319, 83 (1993); A. Kempf, G. Mangano and R. B. Mann, Phys. Rev. D 52, 1108 (1995); F. Scardigli, Phys. Lett. B 452, 39 (1999).
(22)
M. R. Hwang, D. K. Park and E. Jung, Phys. Rev. A 83, 012111 (2011); M. Shamirzai, B. N. Esfahani and M. Soltani, Int. J. Theor. Phys. 51, 787 (2012). S. Khan, N. A. Khan and M. K. Khan, Commun. Theor. Phys. 61, 281 (2014).
(23)
I. Fuentes-Schuller and R. B. Mann, Phys. Rev. Lett. 95, 120404 (2005).
(24)
M. Blasone, F. Dell’Anno, S. De Siena and F. Illuminati, EPL 85, 50002 (2009); M. Blasone, F. Dell’Anno, S. De Siena, M. Di Mauro and F. Illuminati, Phys. Rev. D 77, 096002 (2008).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1)
- 2(2) S. A. Fulling, Phys. Rev. D 7 , 2850 (1973). (3) P. C. W. Davies, J. Phys. A 8 , 609 (1975). (4) W. G. Unruh, Phys. Rev. D 14 , 870 (1976). (5) R. Muller, Phys. Rev. D 56 , 953 (1997). (6) G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 59 , 094004 (1999); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. D 63 , 014010 (2000); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. Lett. 87 , 151301 (2001). (7) V. A. Belinskii, B. M. Karnakov, V. D
- 3(3) P. C. W. Davies, J. Phys. A 8 , 609 (1975). (4) W. G. Unruh, Phys. Rev. D 14 , 870 (1976). (5) R. Muller, Phys. Rev. D 56 , 953 (1997). (6) G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 59 , 094004 (1999); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. D 63 , 014010 (2000); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. Lett. 87 , 151301 (2001). (7) V. A. Belinskii, B. M. Karnakov, V. D. Mur and N. B. Narozhny, Pis’ma Zh. Eksp. Teor. Fiz. 65
- 4(4) W. G. Unruh, Phys. Rev. D 14 , 870 (1976). (5) R. Muller, Phys. Rev. D 56 , 953 (1997). (6) G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 59 , 094004 (1999); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. D 63 , 014010 (2000); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. Lett. 87 , 151301 (2001). (7) V. A. Belinskii, B. M. Karnakov, V. D. Mur and N. B. Narozhny, Pis’ma Zh. Eksp. Teor. Fiz. 65 , 861 [JETP Lett. 65 , 902 (1997)]; A. M. Fedotov, V.
- 5(5) R. Muller, Phys. Rev. D 56 , 953 (1997). (6) G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 59 , 094004 (1999); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. D 63 , 014010 (2000); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. Lett. 87 , 151301 (2001). (7) V. A. Belinskii, B. M. Karnakov, V. D. Mur and N. B. Narozhny, Pis’ma Zh. Eksp. Teor. Fiz. 65 , 861 [JETP Lett. 65 , 902 (1997)]; A. M. Fedotov, V. D. Mur, N. B. Narozhny, V. A. Belinskii and B. M. Karna
- 6(6) G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 59 , 094004 (1999); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. D 63 , 014010 (2000); D. A. T. Vanzella and G. E. A. Matsas, Phys. Rev. Lett. 87 , 151301 (2001). (7) V. A. Belinskii, B. M. Karnakov, V. D. Mur and N. B. Narozhny, Pis’ma Zh. Eksp. Teor. Fiz. 65 , 861 [JETP Lett. 65 , 902 (1997)]; A. M. Fedotov, V. D. Mur, N. B. Narozhny, V. A. Belinskii and B. M. Karnakov, Phys. Lett. A 254 , 126 (1999); N. B. Narozhny,
- 7(7) V. A. Belinskii, B. M. Karnakov, V. D. Mur and N. B. Narozhny, Pis’ma Zh. Eksp. Teor. Fiz. 65 , 861 [JETP Lett. 65 , 902 (1997)]; A. M. Fedotov, V. D. Mur, N. B. Narozhny, V. A. Belinskii and B. M. Karnakov, Phys. Lett. A 254 , 126 (1999); N. B. Narozhny, A. M. Fedotov, B. M. Karnakov, V. D. Mur and V. A. Belinskii, Phys. Rev. D 65 , 025004 (2001); N. B. Narozhny, A. M. Fedotov, B. M. Karnakov, V. D. Mur and V. A. Belinskii, Phys. Rev. D 70 , 048702 (2004). (8) S. W. Hawking,
- 8(8) S. W. Hawking, Nature 248 , 30 (1974); S. W. Hawking, Commun. Math. Phys. 43 , 199 (1975) Erratum: [Commun. Math. Phys. 46 , 206 (1976)]. (9) C. Barcelo, S. Liberati and M. Visser, Living Rev. Rel. 14 , 3 (2011); J. Drori, Y. Rosenberg, D. Bermudez, Y. Silberberg and U. Leonhardt, Phys. Rev. Lett. 122 , 010404 (2019). (10) H. Suzuki and K. Yamada, Phys. Rev. D 67 , 065002 (2003). (11) D. V. Ahluwalia, L. Labun and G. Torrieri, Eur. Phys. J. A 52 , 189 (2016
