
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

71203
10.1038/s41598-024-71203-x
Article
Non-reciprocity in photon polarization based on direction of polarizer under gravitational fields
Noh Hansol 12
Alsing Paul M. 3
Miller Warner A. 3
Ahn Doyeol dahn@uos.ac.kr

13
1 https://ror.org/05en5nh73 grid.267134.5 0000 0000 8597 6969 Department of Electrical and Computer Engineering and Center for Quantum Information Processing, University of Seoul, 163 Seoulsiripdae-Ro, Dongdaemun-Gu, Seoul, 02504 Republic of Korea
2 https://ror.org/04h9pn542 grid.31501.36 0000 0004 0470 5905 Department of Electrical and Computer Engineering, Seoul National University Seoul, Seoul, 08826 Republic of Korea
3 https://ror.org/05p8w6387 grid.255951.f 0000 0004 0377 5792 Department of Physics, Florida Atlantic University, 777 Glades Road, Boca Raton, FL 33431-0991 USA
5 9 2024
5 9 2024
2024
14 208014 4 2024
26 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Unification of gravity with quantum mechanics is still a terra incognita. Photon polarization measurements offer a unique window for probing the interaction between these two fundamental forces. We have revealed that non-reciprocity in the photon polarization angle can arise by tailoring the quantization axis, which corresponds to the direction of polarizer. Due to this non-reciprocity, the measured polarization angle can become ten times larger than that of gravitationally induced frame rotation in both near-Earth and black hole environments. To verify this finding, we propose an astronomical interferometer composed of satellites with the tailored quantization axis, challenging the conventional view of their triviality in closed paths of a photon. Notably, this non-reciprocity can extend to any rotation in the polarization plane, irrespective of the origins, all of which can dictate polarization rotation. Our findings could offer new opportunities for testing fundamental principles in physics.

Subject terms

Quantum information
General relativity and gravity
Theoretical physics
http://dx.doi.org/10.13039/501100003725 National Research Foundation of Korea NRF-2023R1A2C1003570 Ahn Doyeol http://dx.doi.org/10.13039/501100014188 Ministry of Science and ICT, South Korea RS-2023-00225385 Ahn Doyeol http://dx.doi.org/10.13039/100000181 Air Force Office of Scientific Research FA2386-21-1-0089 Ahn Doyeol issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The unification of general relativity and quantum mechanics remains one of the most significant challenge in theoretical physics. General relativity governs the large-scale dynamics of the universe and quantum mechanics dictates the behavior of quantum particles at Planck scales. However, the description of quantum particles under gravitational fields has not been fully understood. Seminal works such as Hawking radiation1 and Unruh effect2, have provided some insights into the interaction between quantum particles and classical gravitational field. Nevertheless, complete understanding of the interplay between two pillars of physics is still lacking, leaving ample room for further discoveries3.

Photons, as inherently massless entities, move along null geodesics, positioning them as sensitive probes of spacetime curvature. The essential roles of photons in advanced experimental systems from interferometers4–6 to quantum communication systems7–23, further emphasize the significance of comprehending the nature of photon states under gravitational influence. Moreover, state-of-the-art advancements in high-precision measurements of relativistic effects, such as the Laser Interferometer Space Antenna (LISA)5, and, like Gravity Probe B4, have made probing gravity’s subtle influence on light increasingly feasible. A recent study employing the Micius satellite examined the decoherence effects of entangled photons in Earth’s gravitational field, providing data that challenges the predictions of event formalism and is consistent with the standard quantum formalism24. These technological developments open new arenas for investigating the intricate interplay between general relativity and quantum mechanics, furthering our understanding of the universe’s fundamental nature.

Furthermore, the polarization of a photon provides a unique window for probing various advanced fields such as quantum gravity25,26, potential violations of Einstein equivalence principle27,28, and broader areas of astrophysics29 and cosmology30, as the polarization changes when the photon moves through spacetime. Nevertheless, while the effects of time dilation have been widely explored, the intricacies of the spin (or polarization) states of a photon in spacetime have been underexplored.

In flat spacetime, the exploration of photon polarization states has yielded key insights, notably into the Wigner rotation—an additional relativistic rotation arising from the non-commutative nature of Lorentz transformations31,32—and its associated quantum phase, the Wigner rotation angle (WRA), which becomes observable with spin states of quantum particles and elucidates how a spin state of a quantum particle interacts with the additional rotation. In this context, introducing a standard frame—where the wave vector is aligned to the third axis (called the quantization axis corresponding to the orientation of a polarizer)—is essential for understanding how quantum particles respond to Lorentz transformations between inertial frames32. It is employed to construct irreducible unitary representations of a single Lorentz transformation. In this frame, the representation decomposes into three sequential Lorentz transformations: from the standard frame to the original, an arbitrary Lorentz transformation in the original, and back to the standard frame. From non-commutativity of these three sequential Lorentz transformations, the Wigner rotation angle (WRA) arises naturally.

The effect of Wigner rotation in flat spacetime on quantum information has been widely studied. Previous study has shown that WRA affects entanglement between the spins of a pair of particles through spin-momentum mixing in moving frames33. Furthermore, it has been shown that momentum-dependent phases (WRA) result in variations in entanglement between a pair of polarization-entangled photon beams in moving frames depending on the boost direction34. Wigner rotation is introduced in both studies to study entanglement of quantum particles observed in moving frames and demonstrate that the WRA can lead to frame-dependent quantum effects.

The phenomenon becomes more intricate in curved spacetime. Employing Einstein equivalence principle and tetrad fields—a set of four orthonormal basis vectors that delineate spacetime curvature into locally flat segments, providing a means to describe various local inertial frames (corresponding to observers in different states of motion: e.g. stationary, freely falling, rotating, etc.)—the WRA can be extended into the geometric phase of photon helicity states in curved spacetime. In this framework, spacetime curvature is encoded in the variations in local inertial frames which manifest as local Lorentz transformations35. Naturally, their irreducible unitary representations and the associated Wigner rotation angles are employed in the quantum field description of local Lorentz transformations as a particle moves in curved spacetime.

For an open path of a photon, in Schwarzschild spacetime, there is a consensus that the phase is zero on open paths only with the proper choice of gauge36. Prior research has accentuated the gravitationally induced WRA can have values of O(10−5) degrees, especially for observers in circular orbits of Earth (e.g. satellites). While these findings underscore the observational potential of WRA, especially in Earth’s vicinity37,38, the WRA of a photon traveling along a closed loop has conventionally been considered trivial, interpreted as a Berry phase influenced by the holonomy group of classical polarization vectors39. In Schwarzschild spacetime, the WRA is known to be zero on closed paths regardless of gauge36. For Kerr spacetime, it has been reported that a gauge-invariant result can be obtained using two future-directed trajectories whose layout is similar to the two arms of a Mach-Zehnder interferometer36. However, given that the effect of the Earth’s spinning angular momentum is negligible, measurement via interferometers has been regarded as infeasible, despite the small values of WRA.

In this study, we explore the effects of the choice of quantization axis, which corresponds to the orientation of a polarizer and has not been previously considered, on the geometric phase (WRA) of a photon helicity state, leading to momentum-dependence in the WRA, for both flat and curved spacetime scenarios. Since the WRA depends on the Lorentz transformations rather than the frame, the choice of various possible standard frames can yield different values for the WRA, resulting in momentum-dependence in WRA. We find that the interplay between the tailored choices of the quantization axis and spacetime curvature results in non-reciprocity of the WRAs. This leads to measurable phase differences in astronomical interferometry conducted near Earth. Also, it is found that the conventional view of trivial WRA on a closed path is recovered with the quantization axis chosen to be parallel to the momentum-vector of a photon.

We propose employing two types of astronomical interferometers with a specifically tailored quantization axis to investigate the interactions between quantum spin states and classical gravitational fields. The first type, the combination of Hong–Ou–Mandel (HOM) and Mach-Zehnder interferometer is appreciated for its precision, particularly its sensitivity and its ability to resist noise from single-photon detection. Through monitoring the coincidence counts, the subtle geometric phase variances induced by spacetime curvature should be measurable. In parallel, we consider a Mach-Zehnder interferometer using a nearly monochromatic light source or a single-photon source. This setup is introduced for validating the equivalence between the Wigner rotation of photon states and the transformation of classical electric fields by comparing the results of the case with the single-photon source, which highlights photon state non-locality.

It is noteworthy that the non-reciprocity arising from the choice of quantization axis can extend beyond the WRA induced by local inertial frame variations. Establishing a standard frame aligned with the chosen quantization axis introduces a sequence of three Lorentz transformations including any arbitrary rotation regardless of its origin: starting from the standard frame, applying the rotation, and then returning to the original frame. Also, it is numerically shown that WRA can be ten times larger than the frame’s rotation angle in both near-Earth and black hole environments by harnessing the non-reciprocity. This suggests that the observed non-reciprocity could be utilized in measuring subtle polarization rotations under varying gravitational conditions, opening doors for further exploration in how light and electromagnetic phenomena behave in different cosmic settings.

Results

Theoretical background

In the study of photon helicity states in spacetime, the standard frame is indispensable, wherein the third axis is designated as the quantization axis. Within this framework, the photon undergoes three sequential Lorentz transformations32; starting from the Lorentz transformation aligning the photon’s wave vector in the standard frame to the transmission direction of a photon, a local Lorentz transformation induced by spacetime curvature, and finally returning to alignment with the quantization axis for comparison with the initial state. The non-commutative nature of the three Lorentz transformations gives rise to the additional rotation, called Wigner Rotation Angle (WRA), underscoring the significant role the quantization axis selection plays in determining the WRA.

Transformation between two local inertial frames, through which a photon passes, is described by an infinitesimal variation of the local orthonormal bases, called tetrads ea^μx, along an infinitesimal displacement of the photon, which, spanning the local inertial frames, satisfy the orthonormal constraints,1 gμνx=ηa^b^eμa^xeνb^x;a^andμ=0,1,2,3,

with ηa^b^=-1,1,1,1. Considering Einstein equivalence principle, the variation of local inertial frames can be viewed as a local Lorentz transformation35,40–43. When a photon state moves along a null geodesic in the geometric optics limit44, photon four-momentum components measured in local frames are transformed under an infinitesimal change of tetrads such that2 ka^(x)→ka^′(x)≡ka^(x)+δka^(x)=δa^b^+λa^b^(x)dξkb^(x)=Λa^b^(x)kb^(x),

where λa^b^(x)=∇kea^ν(x)eνb^(x). The infinitesimal variation rate λa^b^(x) is antisymmetric35,40 and thus can be naturally interpreted as an infinitesimal local Lorentz transformation, Λa^b^(x)=δa^b^+λa^b^(x) ref.35,40,42,43 with Kronecker delta δa^b^. A unitary representation of arbitrary Lorentz transformation, Λ, for a photon state with the helicity, σ, is described by323 U(Λ)k^,σ=Λk0^k0^∑σ′Dσ′σW(Λ,k^)k^Λ,σ′,

where the W(Λ,k^) and DW represent the Wigner’s little group element and its irreducible representation, respectively. Hereafter, the bolded letters indicate the spatial components of the four vector: k^ refers to the spatial components of the wave vector in the local frame and k^Λ signifies the Lorentz-transformed spatial components of the wave vector in the local frame. The necessity of Wigner’s little group in the quantization process lies in its role in steering transformations of particles’ internal degrees of freedom under Lorentz transformations. To formulate a detailed description of the system, it is essential to account for transformations possibly altering the particles’ spin states while preserving momentum invariance. The Wigner’s little group is defined as the subgroup of Poincaré group, which leaves the particle’s four-momentum invariant in a specific reference frame—either the rest frame for a massive particle or the standard frame for a massless one. Accordingly, for photons, the little group W is defined to leave the wave vector invariant in the standard frame conventionally defined as having the photon wave vector aligned along the third axis such that the photon 4-momentum becomes k^std=1,0,0,1. As such, the WRA in Eq. (3) is defined as W=L-1Λk^ΛLk^. Here, Lk^ represents Lorenz transformation from k^std to k^. The non-commutativity of these sequential Lorentz transformations naturally leads to the WRA.

By translating the infinitesimal variations of tetrads into local Lorentz transformations between local inertial frames, the phase of a photon state evolves by an infinitesimal Wigner rotation angle dψ (IWRA) whose ratio with respect to an affine parameter ξ is given by374 dψdξ=λ2^1^+n1^1+n3^-λ2^0^+λ3^2^+n2^1+n3^λ1^0^+λ1^3^=ψ~infgeodetic+ψ~infresidual,

where ni^ is defined as ki^/|ki^| with i^=1,2,and 3. The first term, λ2^1^, in Eq. (4) is the rotation in the 1^-2^ plane about the local third axis, 3^-axis, and the residual terms in the square bracket represents a momentum-dependent phase which we designate as ψ~inf.residual. The first (classical) term corresponds to the rotation of the polarization about the 3^-axis which should be disentangled from other rotations when considering WRA induced by the non-commutativity of local Lorentz transformation. Also, by setting the local third axis to be orthogonal to the observers’ planes, the contribution from the geodetic precession can be isolated in the first term λ2^1^=ψ~inf.geodetic. The total Wigner rotation angle Ψ can be obtained via a time-ordered integration of the IWRAs over the geodesic trajectory xξ of the photon such that5 eiψ(Λ,n→)=eiψgeodetic(Λ,n→)+iψresidual(Λ,n→)=Texpi∫ψ~infgeodetic(Λ(x(ξ)),ni^(ξ))dξ+Texpi∫ψ~infresidual(Λ(x(ξ)),ni^(ξ))dξ,

where T represents the time-ordering operator.

WRA for two open paths of a photon: earth-satellite communication and a spinning black hole

We explore the dependency of the Wigner Rotation Angle (WRA) on the quantization axis by examining two distinct scenarios depicted in Fig. 1: ground-station-to-satellite communication (Fig. 1A), and orbiting observers in the equatorial plane of a rotating black hole (Fig. 1B). Figure 1A illustrates a satellite in the equatorial (yellow line) and polar (light blue line) orbits around Earth. In Fig. 1B, yellow circles represent observers in an equatorial orbit around the black hole such as materials in an accretion disk, stars, or spacecrafts. The black hole scenario is chosen to study WRA in a context where such effects from the dependence are more pronounced. The three sequential Lorentz transformations in these scenarios can be interpreted as follows: prior to photon emission, the polarization or phase of the helicity state is measured in the standard frame where the wave vector is aligned with the local third axis eφ, called the quantization axis (Fig. 1C) and rotated back to the original wave vector (Fig. 1D). Then, as the photons propagate, photons experience local Lorentz transformations induced by the gravitational field. Upon reaching the receiver, the polarization and phase are measured again in the standard frame (Fig. 1E). For brevity, we focus on the case where photon lies in the equatorial plane only. Throughout this paper, a photon field on a curved spacetime is assumed to have a spinor structure32,45 and the (− +  + +) metric signature is used, and the hatted letter represents local flat spacetime. Details of tetrads and photon’s trajectory in a Kerr spacetime can be found in the Method section.Fig. 1 Ground-station-to-satellite (A), and orbiting observers in an equatorial plane of a rotating black hole (B). Yellow and light blue lines in A represent an equatorial and polar orbit, respectively. Yellow circles in B represent observers orbiting around the black hole; we consider the case where photons are emitted from 4.5 times the Schwarzschild radius (rs) away from the black hole. For both cases, various ratios rkφ to kr are considered. As depicted in (C-E), wave vectors of photons are aligned with the local third axis, eφ, i.e., transformed to the standard frame. Before sending a photon, the polarization (i.e. the phase of helicity state) is measured in the standard frame (C). Then, photons are rotated due to the WRA induced by gravity as depicted in (D). At the receiver, polarization and phase are again measured in the standard frame. While both WRAs for the equatorial and polar-orbit cases have the dependence on the impact factors of photon trajectories, the case of observers in the equatorial orbits near the Earth does not give discernible difference due to the weak gravity. For the cases of Earth-satellite-in-polar-orbits and observers in the equatorial plane of M87* black hole, WRAs are plotted as in (F) and (G), respectively. The original image of Earth used in A is obtained from Yeongkwang Kim47.

To isolate contribution from the geodetic precession, another relativistic effect, we intentionally choose the quantization axis to be orthogonal to the orbit plane37. We then examine how the Wigner Rotation Angle (WRA) changes with different orientations of the quantization axis relative to the wave vector. For this, we consider various ratios of rkφ to kr—approximately the impact factors bph of photon’s trajectories42,46—for both cases; for the earth-satellite case, the photon is sent off with the different ratios rkφ /k at the radius of Earth, and for the black hole (BH), photons are sent from an observer at a distance of 4.5 times the Schwarzschild radius from the BH with the various ratios rkφ /k.

In our analysis of satellites in polar orbits, we noted that the signs of the ratios of rkφ to kr significantly influence the WRA, as illustrated in Fig. 1F; WRAs for the negative ratios have ten times larger values than those for the positive. Also, for equatorial orbits, while WRA variations are not discernible for the Earth case, in black hole orbit scenarios, the impact of sign change is apparent for WRAs, indicating a complex interplay with the sign of the impact parameter bph on both orbits in general. The WRAs calculated for observers in equatorial plane of M87* black hole are presented in Fig. 1G.

The observed discrepancy between the two different orbits near Earth is attributed to the momentum-dependence of WRA, i.e., the dependence on the relative direction of the quantization axis compared to the wave vector. Central to this dependence is the choice of the quantization axis, which determines the two Lorentz transformations—the initial transformation from the standard frame to the original frame and the subsequent return transformation to the standard frame after the local Lorentz transformation induced by spacetime curvature. Such intricacies hold significant implications for the polar orbit case, where both wave vectors and quantization axis lie in the same plane. In such a case, when a local Lorentz transformation is applied, WRA changes depending on the sign of z^-component of the wave vector, n3^, while the direction of the z^-component wave vector n3^ alone without any Lorentz transformation does not affect the Wigner rotation angle.

For the specific case illustrated in Fig. 2, as an example, we consider the case where photon momentum lies in y^z^-plane of a local frame. In the absence of any rotation, WRA is null as illustrated in Fig.2A–E. Also, we consider the situation where the local frame undergoes the rotation Ry^Δϕ about the y^-axis by the angle Δϕ, the corresponding rotation R-1k^′—aligning the rotated wave vector k^′=Ry^Δϕk^ to the standard wave vector k^std=1,0,0,1—is directly influenced by the relative orientation of the quantization axis (z^-axis, marked by purple arrows in Fig. 2) with respect to the wave vector k^′, as per the definition of the rotation R-1k^′. This interrelation leads to asymmetry in WRA depending on the relative orientation of the quantization axis as well as Ry^Δϕ, as depicted in Figs. 2J. For infinitesimal y^-axis rotation angle Δϕ, the corresponding WRA is Arg1-Δϕn2^/0.5n2^Δϕ+i1+n3^.Fig. 2 Dependence of WRA on the choice of quantization axis. Using the z^-axis as the quantization axis and the photons on the y^z^-plane as an example. Wave vector k^ and polarization vector p^ are illustrated by red and blue arrows, respectively (A). In the absence of any Lorentz-induced rotation, when the wave vector’s z^-component aligns with the quantization axis, the z^-component matches this axis in the standard frame (B) by rotating the frame about k^×z^, represented by the green arrows (along x^ in A). (C and D) present the case where the direction of z^-component of wave vector is opposite to the quantization axis. Regardless of the direction of the z^-component, the polarization angle (with respect to the z^-axis) is consistent in both cases as shown in (E). However, when the system is rotated by spatial rotation or boost of the frame, z^-component direction of the wave vector, resulting in varied polarization angles in the standard frame. As an example, under a frame rotation about y^-axis, (F-J) show the polarization angle varies in the standard frame depending on the direction of z^-component of wave vector but does not equate to the frame rotation angle due to the non-commutativity. A general observation as depicted in J, is that polarization angle in the standard frame depends on the relative direction of z^-component of wave vector compared to the quantization axis (here the z^-axis), as shown in J, which could lead to asymmetry in WRA. Here, k′^ and p^′ are represent the wave vector and polarization vector after the Lorentz transformation of k^ and p^, respectively.

Violation of local-time reversal symmetry in WRA

In Fig. 3, we illustrate scenarios for the study of local time reversal symmetry of the WRA, where photons are transmitted between satellites in polar orbits. The red shapes represent the wave vector k and 4-velocity of each satellite et^, while blue shapes indicate those of the corresponding time-reversed case (Fig. 3A), and the scenario with photon transmission in the opposite direction -kϕ in the azimuthal plane (Fig. 3B). When dealing with local-time reversal symmetry, the special relativity and equivalence principle govern the behavior of the local wave vectors, represented as k^. Under local-time reversal operator, special relativity dictates that the spatial components of the momentum of the particles reverses direction, while the energy (time component) remains unchanged. In accordance with the equivalence principle, this leads to sign flip of the local spatial components of momentum as shown in Fig. 3A. Also, by the definition of the local-time reversal operator, the sign of timelike component of the tetrad fields should change because time reversal symmetry corresponds to reversing the direction of time, while keeping spatial directions unchanged. These conditions can be achieved by flipping the signs of affine parameter ξ→ξ′=-ξ and the proper time τ→τ′=-τ such that6 eμt^,eμi^=dxν/dτgμν,ei^νgμν→ξ→ξ′τ→τ′dxν/dτ′gμν,ei^νgμν=-dxν/dτgμν,ei^νgμν=-eμt^,eμi^,k^,k^=k^,dxμ/dξeμi^→ξ→ξ′τ→τ′k^,dxμ/dξ′eμi^=k^,-dxμ/dξeμi^=k^,-k^,∇k=dxμ/dξ∇μ→ξ→ξ′τ→τ′dxμ/dξ′∇μ=-dxμ/dξ∇μ=-∇k,

with i^=1,2,and3. Under these transformations, the signs of local infinitesimal rotations, λj^i^=(∇kei^μ)eμj^37,42, are flipped while those of local infinitesimal boosts, λj^t^=(∇ket^μ)eμj^=∇kdxμ/dτeμj^ref.37,42, remain unchanged as shown in Fig. 3C. The corresponding time reversed infinitesimal Wigner rotation angle rate becomes7 dψk^dξ→T^dψ-k^dξ′=-dψ-k^dξ=-λ2^1^+-n1^1-n3^-λ2^0^-λ3^2^+-n2^1-n3^λ1^0^-λ1^3^,

Fig. 3 Local time reversal symmetry violation and WRAs of wave vectors with opposite azimuthal component. Polar orbits of satellites are considered whose 4-velocity vectors are orthogonal to the equatorial plane as shown in (A and B). Under local-time reversal symmetry, photons sent between satellites with wave vector k (depicted with the red arrow in A), the signs of local spatial components should be flipped, while the local time-component (local frequency), remain unchanged, as dictated by the special relativity and equivalence principle. The local-time reversed vectors are depicted with blue in A. Additionally, since by definition local-time reversal symmetry implies a sign flip of time component of tetrad ei^, with no sign change in the same spatial components ei^, the corresponding wave vector in Boyer-Lindquist coordinate bases, local boosts, and rotations should be transformed as shown in (C). When photons are sent with the opposite azimuthal wave vector component, (depicted by the blue arrow in B), the corresponding transformation of the Christoffel symbols of the spherically symmetric Schwarzschild metric are given in (D), where α^is 0, 1, and 2. The original image of Earth used in Fig. 3 is obtained from Yeongkwang Kim47.

(compare with Eq. (4)). Thus, the time-reversed total Wigner rotation angle can be written in terms of the infinitesimal local Lorentz transformation λb^a^ and unit vector ni^ along the direction of photon’s momentum k^ as8 ψk^ξ=∫0ξdψk^dξdξ→T^∫0-ξdψ-k^dξ′dξ′=-∫0ξ-dψ-k^dξdξ=∫0ξλ2^1^dξ-∫0ξn1^1-n3^λ2^0^+λ3^2^+n2^1-n3^-λ1^0^+λ1^3^dξ.

Then, the WRA difference between a path of photons and the time-reversed one represented with red and blue arrows in Fig. 3A, respectively, is non-zero in general such that9 Δψ=∫photon′spathdψk^-T^∫photon′spathdψk^=∫0ξdψk^dξdξ+∫0ξ-dψ-k^dξdξ=∫xa^0xa^ξ∂a^ψk^dxa^-∫xa^0xa^ξ∂a^ψ-k^dxa^=∫0ξn1^λ3^2^1-n3^2+n2^λ1^3^1-n3^2dξ+∫0ξ2n1^n3^λ2^0^1-n3^2-2n2^n3^λ1^0^1-n3^2dξ.

For further verification of local-time reversal symmetry breakdown, in the Supplementary Information (SI), we derive the violation of time reversal symmetry in Euler-Lagrangian framework. In both Eq. (8) and the result of Euler-Lagrangian mechanics, the violation comes from the effect of the choice of quantization axis yielding non-zero components for n1^ and n2^, leading to non-zero WRA difference between the path and the time-revsersed. This symmetry breakdown is caused by the difference between Loretnz transformations Lk^ from k^std to k^ and LTk^ from k^std to Tk^. Here, k^std and Tk^ represent the wave vector of a photon in standard and local-time reversed frame, respectively. Accordingly, when choosing the quantization axis such that the effects of these Lorentz transformations remain in the WRA, the WRA exhibits a localized breakdown in time-reversal symmetry. Accordingly, by selecting a quantization axis parallel to the spatial component of the photon’s wave vector, which avoids retaining the effects of these Lorentz transformations within the WRA, we can ensure that quantum communication systems maintain stability and reliability of quantum states, undisturbed by this localized time-reversal symmetry breakdown.

It should be also highlighted that the local space-inversion operator not only reverses the direction of momentum but also alters the sign of helicity. Consequently, considering the non-reciprocity of the Wigner Rotation Angle (WRA) that maintains the helicity, we focus on the local time-reversal operator instead (See SI).

Asymmetry in WRA under different sign of the momentum component along the quantization axis

For the case of an opposite azimuthal component of wave vector, the corresponding transformation of each component of wave vector, tetrads, and local Lorentz transformation is given in Fig. 3d. In Schwarzschild spacetime (J = 0), only the λi^3^, with î = 1 and 2, are dependent on the azimuthal component kϕ and is given as10 -kϕ1-rsrFi^1-3rsrθ-π2,

where Fi^x=cosxfori^=1sinxfori^=2. Accordingly, the WRA difference ΔΨ between two paths, represented with red (path 1) and blue arrows (path 2) in Fig. 3b, is as follows:11 Δψ=∫path1dψ-∫path2dψ=∫0ξn1^λ3^2^1-n3^2+n2^λ1^3^1-n3^2dξ+∫0ξ2n1^n3^λ2^0^1-n3^2-2n2^n3^λ1^0^1-n3^2dξ.

The relative WRA in Eq. (11) has the same form as that of the case of the local time-reversal symmetry violation in Eq. (9). Ascribed to the asymmetry illustrated in Fig. 2 of the main text, the sign of n3^ affects both scenarios depicted in Fig. 3 through the denominator 1-n3^2 of Eqs. (9), (11). As the value of n3^ increases, the denominator 1-n3^2 decreases, and hence the value of the equations can be enhanced in the case of polar Earth orbits up to measurable orders. This also explains why observers in polar and equatorial orbits of Earth exhibit different WRA dependences based on the sign of azimuthal component wave vector kϕ of photons in the equatorial plane, where in the former case, the quantization axis (azimuthal direction) was set to be orthogonal to orbit plane, while for latter, the quantization axis was set to the polar (or zenith) direction. Since we set the photon trajectories to lie in the equatorial plane, the asymmetry induced by the form of the denominator in Eqs. (9), (11) only affects the case of polar orbits. Moreover, by choosing the quantization axis along the wave vector, for which n1^=n2^=0, the relative WRA of Eqs. (9), (11) become zero, consistent with the result from the previous work39. This means the effect intertwined with local Lorentz transformations induces the asymmetry and violation of local-time reversal symmetry.

Detecting WRA in an astronomical interferometer

The relative WRA induced by the effect from choice of the quantization axis can be measured near Earth in both quantum source and classical light source (Fig. 4). For the former, the combination of a Hong-Ou-Mandel (HOM) and Mach–Zehnder interferometer consisting of four satellites, each on polar orbits, can be utilized.; At Alice’s local frame, two indistinguishable photon helicity states, 1,1aI,bI, are launched, passing through a 50:50 beam splitter, after which a π/2-phase shift on one of the output ports dI is applied (Fig. 4B). In this setup, we select a two-photon state with helicity σ as an input state of the interferometer to observe the WRA Ψ in the coincidence rate. For the helicity state, WRA manifests as a phase shift of σΨ rather than a rotation angle of the polarization. This WRA-induced phase difference between the two paths could lead to the variation of the coincidence rate.Fig. 4 Astronomical Mach-Zehnder interferometry setup designed to observe WRA (depicted in A). For a quantum source scenario, two indistinguishable photons enter two input ports aI and bI , respectively, of the first beam splitter (BS). In front of one output port bI, π/2-phase shift is applied as shown in (B). In David’s frame, photons are detected with single photon detectors (SPDs) after being aligned along the quantization axis depicted in (C). With Hong-Ou-Mandel (HOM) effect, photons emit from each output port with ½ possibility with rkφ/kr = tanα . For a classical source one, light passes through as in (D) and is sent to David with rkφ/kr = tanα. The relative phase difference is measured in David’s frame, assuming that Allice is at the altitude of 300 km, as shown in (E). WRA from Bob and Charlie to David are calculated as integrating infinitesimal WRA along photon’s geodesics from image1 and 2 of Alice to David and subtract the part integrated along images to Bob and Charlie depicted with green lines. The original image of Earth used in Fig. 4 is obtained from Yeongkwang Kim47.

Using the creation operators of the two input aI† and bI† and two output ports cI† and dI† of the 50:50 beam splitter, the two-photon helicity states at each in input port can be written in terms of output-port operators as12 1,1aI,bI=aI†bI†0,0=12cI†+dI†cI†-dI†0,0=12cI†20,0-12dI†20,0.=122,0cI,dI-120,2cI,dI

Correspondingly, the state after the phase shift on an output port dI becomes 2,0cI,dI-eiπ/20,2cI,dI/2.

Following this, the photons from each output port are sent from Alice to Bob to David, and from Alice to Charlie to David with rkϕ/kr=tanαas shown in Fig. 4A), which is approximately the impact parameter (bph) of a photon’s trajectory. Then, considering the relative Wigner rotation angle ΔΨ between the two arms of the interferometer, the photon states passing through output ports cII and dII of the second beam splitter (Fig. 4C) are as follows:13 eσiΔψ2cI†20,0-i12dI†20,0=14eσiΔψcII†2+2cII†dII†+dII†20,0-i4cII†2-2cII†dII†+dII†20,0=-iieσiΔψ+14cII†2+dII†2+-1+ieσiΔψ2cII†dII†0,0=-iieσiΔψ+1222,0cIIdII+0,2cIIdII+-1+ieσiΔψ21,1cIIdII,

where the relative WRA Δψ has the same form of Eq. (9). The relative WRA Δψchanges the linear combination of two indistinguishable photon states in Eq. (13), leading to the coincidence rate at David’s frame, (1-sin(σΔψ))/2, depending on both the helicity and relative WRA.

In addition, we consider a Mach-Zehnder interferometer with a nearly monochromatic light source or a single-photon source. In Alice’s local frame the classical light or a single photon enters through one of the input ports of the first beam splitter (as illustrated in the Fig. 4D) and then in David’s frame, the relative phase difference is measured after the light or the photon passes through the second beam splitter after aligning wave vectors to the quantization axis. By comparing the results from two different sources, it is possible to validate the equivalence between the Wigner rotation of photon states with the transformation of classical electric fields (proven in the Method)34,37,42. Moreover, this Mach-Zehnder interferometer could underscore the interplay of the WRA and non-local properties of a photon when a single-photon source is employed.

The relative WRA differences between the two arms of an interferometer near the Earth have a small order of 10−4 degrees compared to a few degrees (See Fig.S3D) for an interferometer around M87*. Considering the precision of LIGO and LISA where a few tens of attometer and picometer precision are required with near-infrared light respectively, O(10−4) degrees of the relative WRA (shown in Fig. 4E) should be measurable. Moreover, while the spinning angular momentum of M87* black hole would lead to a noticeable difference in the WRA, on the O(10–2) degrees, but that of Earth does not significantly affect the WRA to measurable order for both equatorial and polar orbits, as supported by the results of the Gravity Probe B (See SI).

Discussion

Polarization measurements of photons offer a room a window to unveil a variety of phenomena, ranging from the quantum gravity26 to the possibility of violations of the Einstein Equivalence Principle27. The polarization of a photon undergoes alterations as it interacts with the curvature of spacetime during its passage. In this work, we uncover a new facet: the non-reciprocity that emerges from the strategic choice of a quantization axis. Numerical calculations show that in different environments, near-Earth to near a black hole, the WRA can be measured as ten times larger than the rotation angle of the frame by exploiting this non-reciprocity (See SI). This finding suggests that the non-reciprocity could be harnessed to detect subtle polarization rotations under varying gravitational conditions.

This non-reciprocity extends its influence beyond the usual scope of the Wigner rotation angle (WRA), associated with frame variations in spacetime. When a standard frame is established with the chosen quantization axis, there would be a sequence of three Lorentz transformations including a polarization rotation, regardless of its origins: transforming the standard frame to the original, applying the polarization rotation around the photon’s wave vector, and returning to the standard frame. While in this work we only consider polarization rotations in classical gravitational fields, following the Einstein Equivalence Principle, polarization rotation can be induced by various factors such as quantum gravity and violations of the Einstein Equivalence Principle—If the polarization rotations are considered as an SO(2) rotation on the polarization plane in each formalism, we could still utilize the non-reciprocity by introducing a standard frame defined with a tailored direction of the quantization axis.

It has been shown that specific violations of the Einstein Equivalence Principle (EEP) would lead to a rotation in the plane of polarization of radiation from distant radio sources28. In the effective field theory approach within quantum gravity, Lorentz- and CPT-violating dispersion relations can arise, leading to different group velocities for different polarization states26. This results in a polarization rotation during propagation, as described by a modified dispersion relation. These additional rotations in the polarization plane can be examined using the standard frame with a tailored direction of the quantization axis, thus allowing us to utilize the non-reciprocity. Moreover, validating the non-reciprocity induced from the choice of quantization axis could be feasible on an optical table. The non-reciprocity itself can also arise from inherent non-commutativity of rotations with no need for a boost as depicted in Fig. 2, which can be effectively realized using polarizers and mirrors by effectively simulating rotations. Thus, optical table experiments could demonstrate and validate the non-reciprocity in a controlled laboratory environment.

Method

All the data in the GitHub link have been obtained using a Mathematica notebook (Wolfram Research Europe Ltd, Long Hanborough, UK).

Geodesics in Kerr spacetime

The Kerr metric is applied to model spacetimes considered in this work with spin angular momentum (J/Mc = a) of a gravitating object, which is given by4814 ds2=-1-rsrΣc2dt2+ΣΔdr2+Σdθ2+r2+a2+rsra2Σsin2θsin2θdϕ2-2rsrasin2θΣcdtdϕ,

where Σ≡r2+a2cos2θ, Δ≡r2-rsr+a2, and rs≡2GM/c2 is the Schwarzschild radius. G and c are the gravitational constant and the speed of light. J and M are the angular momentum and the mass of the gravitating object, respectively.

The geodesics of a particle in Kerr spacetime is described by4815 dtdξdrdξdθdξdϕdξ=1ΔΣEr2+a22-Δa2sin2θ-2MraΦ±1ΣEr2+a22-aΦ2-ΔK+δ1r2±1ΣK-δ1a2cos2θ-aEsinθ-Φ/sinθ21ΔΣ2MraE+Σ-2MrΦ/sin2θ.

Here, K=C+(Φ-aE)2. The parameters (ξ, δ1) are (affine parameter, 0) or (proper time, 1) for null or time-like geodesics, respectively. E, φ, and C are the conserved quantities, corresponding to energy at infinity, axial angular momentum of an orbiting particle, and Carter’s (fourth) constant, defined as C=dθ/dξ2+cos2θa2δ12-E2+Φcosecθ248–50, respectively.

Co-moving tetrads with a satellite, parallel transported along its geodesic

To describe local inertial frames of satellites, we obtain co-moving and non-spinning tetrads {ei^=(ei^t,ei^r,ei^θ,ei^ϕ)|i^=0,1,2,and3} by transforming Marck’s {λi^=(λi^0,λi^1,λi^2,λi^3)|i^=0,1,2,and3}48 from Carter’s symmetric bases (ω0,ω1,ω2,ω3) back to the Boyer-Lindquist coordinate ones (dt,dr,dθ,dϕ). Co-moving condition ensures that the reference frame moves along with the satellite, providing a local rest frame for measurements, and non-spinning condition is essential to eliminate any arbitrary rotational effects within the reference frame. The Marck’s solution satisfies the two conditions by obtaining the tetrads parallel-transported along a test particle’s geodesics49, specifically those of satellites in our case.

In details, the Marck’s tetrads are defined as follows49:16 λ0^λ1^λ2^λ3^=10000cosΨ0-sinΨ00100sinΨ0cosψΔ/Σe0^t-asin2θe0^ϕΣ/Δe0^rΣe0^θae0^t-r2+a2e0^ϕsinθe0^t/ΣαΣ/KΔre0^rα1/KΣΔrEr2+a2-aΦβ1/KΣacosθaEsinθ-Φ/sinθ-βΣ/Kacosθe0^θΣ/KΔacosθe0^r1/KΣΔacosθEr2+a2-aΦ-1/KΣraEsinθ-Φ/sinθΣ/Kre0^θα1/ΣΔEr2+a2-aΦ1/KΣΔacosθEr2+a2-aΦβΣe0^θβ1/ΣaEsinθ-Φ/sinθ.

The parameter Ψ is introduced for the parallel-transport condition of spacelike component of tetrads along the geodesics of observers, i.e., satellites. The parameter Ψ can be obtained by integrating the following equation49,17 dΨdξ=K1/2ΣEr2+a2-aΦr2+K+aΦ-aEsin2θK-a2cos2θ,

we solve the Eq. (17) with respect to (r, θ, φ) from (the radius of Earth, Π/2, −Π) to (an altitude of 36000 km, π, π) with the intervals (100 km, Π/10, Π/20) for the case of near-Earth and for the black hole, we solve the equation from (4.5rs, Π/2,−Π) to (4.5rs + 2 × 1011 km, Π, Π) with the intervals (2 × 109 km, Π/10, Π/20).

For equatorial orbits, as r and θ should be constant, Eq. (17) can be rewritten as18 dΨdϕ=K1/2ΣEr2+a2-aΦr2+K+aΦ-aEsin2θK-a2cos2θdξdϕ.

In the case of polar orbits, if Ψ depends on r, θ, and φ, dΨ /dξ can be rewritten as (∂Ψ(r,θ,ϕ)/∂r)(dr/dξ)+(∂Ψ(r,θ,ϕ)/∂ϕ)(dϕ/dξ)+(∂Ψ(r,θ,ϕ)/∂θ)(dθ/dξ) using the chain rule. Since the four-velocity vector is independent of φ (as shown in Eq. (15)), the second term (∂Ψ(r,θ,ϕ)/∂ϕ)(dϕ/dξ) cannot be generally zero. Given that Eq. (17) are also independent of φ, this leads to a contradiction. Therefore, we conclude Ψ is independent of φ. Accordingly, considering, as the case of equatorial orbit, r is constant on the orbits, the following equation dictates the parameter Ψ of the tetrads parallel-tranported on polar orbits:19 dΨdθ=K1/2ΣEr2+a2-aΦr2+K+aΦ-aEsin2θK-a2cos2θdξdθ.

Then, transforming the Marck’s from Carter’s symmetric bases to the Boyer-Lindquist coordinate ones, we obtain the tetrads eî used in this paper such that20 ei^tei^rei^θei^ϕ=a2+r2ΔΣ00-asinθΣ0ΔΣ00aΔΣ00-cscθΣ00-1Σ0λi^(0)λi^(1)λi^(2)λi^(3).

The tetrads (et^,er^,eθ^,eϕ^) obtained from Eq. (20) are asymptotically parallel to the unit vectors of global coordinates (∂t,∂r,∂θ,∂ϕ) as r goes to infinity. This transformation also ensures the tetrads obtained from Marck’s become equivalent to those used to study of WRA in Schwarzschild spacetime in our previous work37 by setting spin angular momentum (J/Mc = a) of gravitating object as zero.

Also, we additionally rotate the tetrads about the local first axis to choose the local third axis as that of geodetic precession, compensating the frame-dragging effect. The corresponding rotation angle Ξ is obtained by21 -cosΞe2^θ+sinΞe3^θ=0.

This approach is used to ensure a consistent comparison with our previous work37 on the Schwarzschild metric, In that previous work, the local third axis, the quantization axis, was chosen to be the geodetic precession axis to isolate the geodetic precession.

Wigner rotation and the complementarity of classical and quantum theories: equivalence between Wigner rotation of a photon and classical polarization rotation

The polarization vector, denoted as ε(p̂, σ =  ± 1), is defined with R(p̂) rotating the wave vector k̂std = (1,0,0,1) in the standard frame into the direction of arbitrary momentum p̂ such that3222 εp^,σ=±1=Rp^2εk^std,σ=±1=Rp^201±i0.

Here, σ represents the helicity. Throughout this paper, the hatted letter represents local or global flat spacetime. If the polarization vector is treated as a four-vector, it can be directly subjected to a Lorentz transformation ʌ and be rewritten as follows:23 Λεσ=±1p^=LΛp^L-1Λp^ΛLp^L-1p^εp^,σ=±1=LΛp^WΛ,p^εk^std,σ=±1=LΛp^Sα,βRzψΛ,p^/p^εk^std,σ=±1=LΛp^eiσψΛ,p^/p^εk^std,σ=±1+α+iβ2k^stdk^std=RΛp^eiσψΛ,p^/p^εk^std,σ=±1+α+iβ2k^stdLΛp^k^std=eiσψΛ,p^/p^εΛp^,σ=±1+α+iβ2k^stdΛp^,

with Sα,β=1+α2+β22αβ-α2+β22α10-αβ01-βα2+β22αβ1-α2+β22 ref.32. The Wigner’s little group W(ʌ,p̂), defined as L-1(ʌp̂)ʌL(ʌ), which is a subgroup of the Lorentz group, leaves the wave vector in the standard frame invariant and can be decomposed into S(α,β)Rz(ψ(ʌ, p̂/|p̂|))32 where S(α,β) is a subgroup isomorphic to the translation of a Euclidean plane and Rz(ψ) represents the rotation about z-axis by ψ. The rotation angle ψ(ʌ, p̂/|p̂|) corresponds to the Wigner rotation angle. Here, we use the decomposition of the Lorentz transformation L(p̂), which maps the wave vector k̂std in the standard frame to a wave vector p̂, into R(p̂)Bz(|p̂|) along with the invariance of the polarization vector under a boost Bz(|p̂|) in the z-direction in the standard frame51. As the polarization vector is essentially a three-dimensional vector in the spatial part of the spacetime, the corresponding representation of Lorentz transformation U(ʌ) can be defined as24 UΛεΛp^,σ=±1=Λεp^,σ=±1-Λεp^,σ=±10^Λp^0^Λp^=eiσψΛ,p^/p^εΛp^,σ=±1.

Here, it is used that the time component of Lorentz transformed Λε0^ corresponds to +i2k^stdΛp^0^, since polarization vector ε has no time component.

In classical description, if the light can be described as monochromatic and circular polarized wave and then its corresponding potential four vector φ is the function of inner product between momentum p̂ and position x̂—φ=ε(p̂, σ =  ± 1)ϕ(p̂∙x̂)—the electric fields have the same relation of Eq. (24) under the Lorentz transformation Λ. It implies that a Lorentz transformation of classical light, ensemble of photons, can have the Wigner rotation angle within the framework of Coulomb gauge fixing. Here, ε(p̂, σ =  ± 1) is the four vector defined as Eq. (22) and ϕ(p̂∙x̂) is the scalar function. In details, the form of gauge-independent electromagnetic field derived from a Lorentz transformed potential field is as follows34:25 Ei^x^=-Λεp^,σ=±1i^∂0^ϕp^·x^-Λεp^,σ=±10^∂i^ϕp^·x^∝Λεp^,σ=±1i^Λp^0^-Λεp^,σ=±10^Λp^i^∝Λεp^,σ=±1i^-Λεp^,σ=±10^Λp^0^Λp^i^,

where Eî(x̂) is î-th component of electric vector field. It is worth mentioning, while the Wigner rotation can be described within the classical framework of electromagnetic waves, classical description cannot account for a super-positioned state of a single photon—a quantum mechanical phenomenon—but only the average effect of Lorentz transformations on a photon ensemble. In other words, classical description does not distinguish whether each photon itself becomes super-positioned or photons are statistically in one of two different spin states. Accordingly, the correspondence between classical and quantum descriptions of Wigner rotation of a photon showcases the deep connection between the two frameworks and further illustrates the complementary nature of classical and quantum theories in explaining physical phenomena within their respective domains of applicability rather than negating the other.

Furthermore, to derive the equivalent relation, it is imperative that the potential field of light be a four-vector, but not the electric or magnetic fields, and that the light exhibits circular polarization—the unique discretized form of polarization corresponding to the quantized helicity. This result underscores the primacy of the Wigner rotation description for the properties of the electric field under a Lorentz transformation since the Wigner rotation is rooted in intrinsic properties of the Wigner’s little group depending on a particle’s spin and the relative direction of a photon’s path compared to the quantization axis of the spin32,52. This perspective could furnish a more foundational approach to understanding the dynamics of spacetime transformations.

This Wigner rotation on a quantum state should be differentiated from that in the context of Thomas precession observed in spinning particles or rotating macroscopic gyroscopes undergoing curvilinear motion53,54. When an observer in an accelerated frame of reference carries out a sequence of infinitesimal Lorentz transformations (boosts), an additional rotation manifests due to the non-commutativity of these transformations (a well-known special relativistic effect). Meanwhile, in the quantum domain, while Wigner rotation itself is induced from the non-commutativity of boosts as Thomas precession, its application to quantum states results in super-positioned spin states of a photon with different phases, called Wigner rotation angles,—distinguished from those which can be obtained from polarizer rotation. The Wigner rotation angle becomes observable via spin of particles, i.e., Wigner rotation induced from the little group does not affect spinless particles. Thus, the distinctiveness of the Wigner rotation angle is characterized by the introduction of superposition spin states without any rotation of the polarizer.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71203-x.

Acknowledgements

This work was supported by Korean Ministry of Science and ICT through the National Information Agency (NIA). D.A. was also supported by Korea National Research Foundation (NRF) grant No. NRF-2023R1A2C1003570, ICT R&D program of MSIT/IITP 2021-0-01810, RS-2023-00225385, RS-2024-00422330, AFOSR grant FA2386-21-1-0089, and AFOSR grant FA2386-22-1-4052. P.M.A. would like to acknowledge support from the Air Force Office of Scientific Research (AFOSR). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the Air Force Research Laboratory.

Author contributions

H.N. introduced the HOM interferometer concept for measuring WRA differences, conducted numerical simulations, and developed a theoretical framework to demonstrate the non-reciprocity and violation of local-time reversal symmetry in WRA. P.M.A., W.M., and D.A. reviewed and validated the paper’s arguments, suggesting an exploration of WRA in Kerr spacetime. All the authors contributed to the manuscript’s composition.

Funding

National Research Foundation of Korea,NRF-2023R1A2C1003570,Ministry of Science and ICT,South Korea,RS-2023-00225385,Air Force Office of Scientific Research,FA2386-21-1-0089.

Data availability

The data supporting the plots and other findings presented in this study are openly accessible on our GitHub repository at the following link: https://github.com/gksshcks/WRA_For_different_impact_parameters.

Code availability

The Mathematica codes that support the data used in this paper and other finding are available from the corresponding author upon reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Hawking SW Particle creation by black holes Commun. Math. Phys. 1975 43 199 220 10.1007/BF02345020
Hawking, S. W. Particle creation by black holes. Commun. Math. Phys. 43, 199–220. 10.1007/BF02345020 (1975).10.1007/BF02345020
2. Unruh WG Notes on black-hole evaporation Phys. Rev. D 1976 14 870 10.1103/PhysRevD.14.870
Unruh, W. G. Notes on black-hole evaporation. Phys. Rev. D 14, 870. 10.1103/PhysRevD.14.870 (1976).10.1103/PhysRevD.14.870
3. Murugan J Weltman A Ellis GF Foundations of Space and Time: Reflections on quantum Gravity 2012 Cambridge University Press
Murugan, J., Weltman, A. & Ellis, G. F. Foundations of Space and Time: Reflections on quantum Gravity (Cambridge University Press, 2012).
4. Everitt CWF Gravity probe B: Final results of a space experiment to test general relativity Phys. Rev. Lett. 2011 106 221101 10.1103/PhysRevLett.106.221101 21702590
Everitt, C. W. F. et al. Gravity probe B: Final results of a space experiment to test general relativity. Phys. Rev. Lett. 106, 221101. 10.1103/PhysRevLett.106.221101 (2011).21702590 10.1103/PhysRevLett.106.221101
5. Amaro-Seoane P Laser interferometer space antenna arXiv preprint arXiv:1702.00786 2017 10.48550/arXiv.1702.00786
Amaro-Seoane, P. et al. Laser interferometer space antenna. arXiv preprint arXiv:1702.0078610.48550/arXiv.1702.00786 (2017).10.48550/arXiv.1702.00786
6. Abbott BP Observation of gravitational waves from a binary black hole merger Phys. Rev. Lett. 2016 116 061102 10.1103/PhysRevLett.116.061102 26918975
Abbott, B. P. et al. Observation of gravitational waves from a binary black hole merger. Phys. Rev. Lett. 116, 061102. 10.1103/PhysRevLett.116.061102 (2016).26918975 10.1103/PhysRevLett.116.061102
7. Yin J Entanglement-based secure quantum cryptography over 1120 kilometres Nature 2020 582 501 505 10.1038/s41586-020-2401-y 32541968
Yin, J. et al. Entanglement-based secure quantum cryptography over 1120 kilometres. Nature 582, 501–505. 10.1038/s41586-020-2401-y (2020).32541968 10.1038/s41586-020-2401-y
8. Yin J Satellite-based entanglement distribution over 1200 kilometers Science 2017 356 1140 10.1126/science.aan3211 28619937
Yin, J. et al. Satellite-based entanglement distribution over 1200 kilometers. Science 356, 1140. 10.1126/science.aan3211 (2017).28619937 10.1126/science.aan3211
9. Ren J-G Ground-to-satellite quantum teleportation Nature 2017 549 70 73 10.1038/nature23675 28825708
Ren, J.-G. et al. Ground-to-satellite quantum teleportation. Nature 549, 70–73. 10.1038/nature23675 (2017).28825708 10.1038/nature23675
10. Liao S-K Satellite-to-ground quantum key distribution Nature 2017 549 43 47 10.1038/nature23655 28825707
Liao, S.-K. et al. Satellite-to-ground quantum key distribution. Nature 549, 43–47. 10.1038/nature23655 (2017).28825707 10.1038/nature23655
11. Liao S-K Satellite-relayed intercontinental quantum network Phys. Rev. Lett. 2018 120 030501 10.1103/PhysRevLett.120.030501 29400544
Liao, S.-K. et al. Satellite-relayed intercontinental quantum network. Phys. Rev. Lett. 120, 030501. 10.1103/PhysRevLett.120.030501 (2018).29400544 10.1103/PhysRevLett.120.030501
12. Buttler WT Free-space quantum-key distribution Phys. Rev. A 1998 57 2379 2382 10.1103/PhysRevA.57.2379
Buttler, W. T. et al. Free-space quantum-key distribution. Phys. Rev. A 57, 2379–2382. 10.1103/PhysRevA.57.2379 (1998).10.1103/PhysRevA.57.2379
13. Gibson G Free-space information transfer using light beams carrying orbital angular momentum Opt. Express 2004 12 5448 5456 10.1364/OPEX.12.005448 19484105
Gibson, G. et al. Free-space information transfer using light beams carrying orbital angular momentum. Opt. Express 12, 5448–5456. 10.1364/OPEX.12.005448 (2004).19484105 10.1364/OPEX.12.005448
14. Schmitt-Manderbach T Experimental demonstration of free-space decoy-state quantum key distribution over 144 km Phys. Rev. Lett. 2007 98 010504 10.1103/PhysRevLett.98.010504 17358463
Schmitt-Manderbach, T. et al. Experimental demonstration of free-space decoy-state quantum key distribution over 144 km. Phys. Rev. Lett. 98, 010504. 10.1103/PhysRevLett.98.010504 (2007).17358463 10.1103/PhysRevLett.98.010504
15. Aolita L Walborn SP Quantum communication without alignment using multiple-Qubit single-photon states Phys. Rev. Lett. 2007 98 100501 10.1103/PhysRevLett.98.100501 17358518
Aolita, L. & Walborn, S. P. Quantum communication without alignment using multiple-Qubit single-photon states. Phys. Rev. Lett. 98, 100501. 10.1103/PhysRevLett.98.100501 (2007).17358518 10.1103/PhysRevLett.98.100501
16. Laing A Scarani V Rarity JG O’Brien JL Reference-frame-independent quantum key distribution Phys. Rev. A 2010 82 012304 10.1103/PhysRevA.82.012304
Laing, A., Scarani, V., Rarity, J. G. & O’Brien, J. L. Reference-frame-independent quantum key distribution. Phys. Rev. A 82, 012304. 10.1103/PhysRevA.82.012304 (2010).10.1103/PhysRevA.82.012304
17. Zhang P Reference-frame-independent quantum-key-distribution server with a telecom tether for an on-chip client Phys. Rev. Lett. 2014 112 130501 10.1103/PhysRevLett.112.130501 24745397
Zhang, P. et al. Reference-frame-independent quantum-key-distribution server with a telecom tether for an on-chip client. Phys. Rev. Lett. 112, 130501. 10.1103/PhysRevLett.112.130501 (2014).24745397 10.1103/PhysRevLett.112.130501
18. Erven C Studying free-space transmission statistics and improving free-space quantum key distribution in the turbulent atmosphere New J. Phys. 2012 14 123018 10.1088/1367-2630/14/12/123018
Erven, C. et al. Studying free-space transmission statistics and improving free-space quantum key distribution in the turbulent atmosphere. New J. Phys. 14, 123018. 10.1088/1367-2630/14/12/123018 (2012).10.1088/1367-2630/14/12/123018
19. Vallone G Free-space quantum key distribution by rotation-invariant twisted photons Phys. Rev. Lett. 2014 113 060503 10.1103/PhysRevLett.113.060503 25148310
Vallone, G. et al. Free-space quantum key distribution by rotation-invariant twisted photons. Phys. Rev. Lett. 113, 060503. 10.1103/PhysRevLett.113.060503 (2014).25148310 10.1103/PhysRevLett.113.060503
20. Perdigues Armengol JM Quantum communications at ESA: Towards a space experiment on the ISS Acta Astronaut. 2008 63 165 178 10.1016/j.actaastro.2007.12.039
Perdigues Armengol, J. M. et al. Quantum communications at ESA: Towards a space experiment on the ISS. Acta Astronaut. 63, 165–178. 10.1016/j.actaastro.2007.12.039 (2008).10.1016/j.actaastro.2007.12.039
21. Ursin R Space-quest, experiments with quantum entanglement in space Europhys. News 2009 40 26 29 10.1051/epn/2009503
Ursin, R. et al. Space-quest, experiments with quantum entanglement in space. Europhys. News 40, 26–29. 10.1051/epn/2009503 (2009).10.1051/epn/2009503
22. Rideout D Fundamental quantum optics experiments conceivable with satellites—reaching relativistic distances and velocities Class. Quantum Gravity 2012 29 224011 10.1088/0264-9381/29/22/224011
Rideout, D. et al. Fundamental quantum optics experiments conceivable with satellites—reaching relativistic distances and velocities. Class. Quantum Gravity 29, 224011. 10.1088/0264-9381/29/22/224011 (2012).10.1088/0264-9381/29/22/224011
23. Bruschi DE Ralph TC Fuentes I Jennewein T Razavi M Spacetime effects on satellite-based quantum communications Phys. Rev. D 2014 90 045041 10.1103/PhysRevD.90.045041
Bruschi, D. E., Ralph, T. C., Fuentes, I., Jennewein, T. & Razavi, M. Spacetime effects on satellite-based quantum communications. Phys. Rev. D 90, 045041. 10.1103/PhysRevD.90.045041 (2014).10.1103/PhysRevD.90.045041
24. Xu P Satellite testing of a gravitationally induced quantum decoherence model Science 2019 366 132 135 10.1126/science.aay5820 31604316
Xu, P. et al. Satellite testing of a gravitationally induced quantum decoherence model. Science 366, 132–135. 10.1126/science.aay5820 (2019).31604316 10.1126/science.aay5820
25. Myers RC Pospelov M Ultraviolet modifications of dispersion relations in effective field theory Phys. Rev. Lett. 2003 90 211601 10.1103/PhysRevLett.90.211601 12786546
Myers, R. C. & Pospelov, M. Ultraviolet modifications of dispersion relations in effective field theory. Phys. Rev. Lett. 90, 211601. 10.1103/PhysRevLett.90.211601 (2003).12786546 10.1103/PhysRevLett.90.211601
26. Toma K Strict limit on CPT violation from polarization of γ-ray bursts Phys. Rev. Lett. 2012 109 241104 10.1103/PhysRevLett.109.241104 23368301
Toma, K. et al. Strict limit on CPT violation from polarization of γ-ray bursts. Phys. Rev. Lett. 109, 241104. 10.1103/PhysRevLett.109.241104 (2012).23368301 10.1103/PhysRevLett.109.241104
27. Li C Zhao H Cai Y-F New test on the Einstein equivalence principle through the photon ring of black holes Phys. Rev. D 2021 104 064027 10.1103/PhysRevD.104.064027
Li, C., Zhao, H. & Cai, Y.-F. New test on the Einstein equivalence principle through the photon ring of black holes. Phys. Rev. D 104, 064027. 10.1103/PhysRevD.104.064027 (2021).10.1103/PhysRevD.104.064027
28. Carroll SM Field GB Einstein equivalence principle and the polarization of radio galaxies Phys. Rev. D 1991 43 3789 10.1103/PhysRevD.43.3789
Carroll, S. M. & Field, G. B. Einstein equivalence principle and the polarization of radio galaxies. Phys. Rev. D 43, 3789. 10.1103/PhysRevD.43.3789 (1991).10.1103/PhysRevD.43.3789
29. Tinbergen J Astronomical polarimetry Astron. Polarim. 2005 10.1017/CBO9780511525100.005
Tinbergen, J. Astronomical polarimetry. Astron. Polarim.10.1017/CBO9780511525100.005 (2005).10.1017/CBO9780511525100.005
30. Toma K Statistical properties of gamma-ray burst polarization Astrophys. J. 2009 698 1042 1053 10.1088/0004-637x/698/2/1042
Toma, K. et al. Statistical properties of gamma-ray burst polarization. Astrophys. J. 698, 1042–1053. 10.1088/0004-637x/698/2/1042 (2009).10.1088/0004-637x/698/2/1042
31. Vigoureux J Calculations of the Wigner angle Eur. J. Phys. 2001 22 149 10.1088/0143-0807/22/2/307
Vigoureux, J. Calculations of the Wigner angle. Eur. J. Phys. 22, 149. 10.1088/0143-0807/22/2/307 (2001).10.1088/0143-0807/22/2/307
32. Weinberg S The Quantum Theory of Fields 1995 Cambridge University Press
Weinberg, S. The Quantum Theory of Fields (Cambridge University Press, 1995).
33. Gingrich RM Adami C Quantum entanglement of moving bodies Phys. Rev. Lett. 2002 89 270402 10.1103/PhysRevLett.89.270402 12513186
Gingrich, R. M. & Adami, C. Quantum entanglement of moving bodies. Phys. Rev. Lett. 89, 270402. 10.1103/PhysRevLett.89.270402 (2002).12513186 10.1103/PhysRevLett.89.270402
34. Gingrich RM Bergou AJ Adami C Entangled light in moving frames Phys. Rev. A 2003 68 042102 10.1103/PhysRevA.68.042102
Gingrich, R. M., Bergou, A. J. & Adami, C. Entangled light in moving frames. Phys. Rev. A 68, 042102. 10.1103/PhysRevA.68.042102 (2003).10.1103/PhysRevA.68.042102
35. Terashima H Ueda M Einstein-Podolsky-Rosen correlation in a gravitational field Phys. Rev. A 2004 69 032113 10.1103/PhysRevA.69.032113
Terashima, H. & Ueda, M. Einstein-Podolsky-Rosen correlation in a gravitational field. Phys. Rev. A 69, 032113. 10.1103/PhysRevA.69.032113 (2004).10.1103/PhysRevA.69.032113
36. Brodutch A Demarie TF Terno DR Photon polarization and geometric phase in general relativity Phys. Rev. D 2011 84 104043 10.1103/PhysRevD.84.104043
Brodutch, A., Demarie, T. F. & Terno, D. R. Photon polarization and geometric phase in general relativity. Phys. Rev. D 84, 104043. 10.1103/PhysRevD.84.104043 (2011).10.1103/PhysRevD.84.104043
37. Noh H Alsing PM Ahn D Miller WA Park N Quantum mechanical rotation of a photon polarization by Earth’s gravitational field npj Quantum Inf. 2021 10.1038/s41534-021-00471-6
Noh, H., Alsing, P. M., Ahn, D., Miller, W. A. & Park, N. Quantum mechanical rotation of a photon polarization by Earth’s gravitational field. npj Quantum Inf.10.1038/s41534-021-00471-6 (2021).10.1038/s41534-021-00471-6
38. Dahal PK Terno DR Polarization rotation and near-Earth quantum communications Phys. Rev. A 2021 10.1103/PhysRevA.104.042610
Dahal, P. K. & Terno, D. R. Polarization rotation and near-Earth quantum communications. Phys. Rev. A10.1103/PhysRevA.104.042610 (2021).10.1103/PhysRevA.104.042610
39. Lindner NH Peres A Terno DR Wigner’s little group and Berry’s phase for massless particles J. Phys. A Math. Gen. 2003 36 L449 L454 10.1088/0305-4470/36/29/101
Lindner, N. H., Peres, A. & Terno, D. R. Wigner’s little group and Berry’s phase for massless particles. J. Phys. A Math. Gen. 36, L449–L454. 10.1088/0305-4470/36/29/101 (2003).10.1088/0305-4470/36/29/101
40. Nakahara M Geometry, Topology and Physics 2003 CRC Press
Nakahara, M. Geometry, Topology and Physics (CRC Press, 2003).
41. Schwinger J Energy and momentum density in field theory Phys. Rev. 1963 130 800 805 10.1103/PhysRev.130.800
Schwinger, J. Energy and momentum density in field theory. Phys. Rev. 130, 800–805. 10.1103/PhysRev.130.800 (1963).10.1103/PhysRev.130.800
42. Alsing PM Stephenson GJ Jr The Wigner rotation for photons in an arbitrary gravitational field arXiv preprint arXiv:0902.1399 2009 10.48550/arXiv.0902.1399
Alsing, P. M. & Stephenson, G. J. Jr. The Wigner rotation for photons in an arbitrary gravitational field. arXiv preprint arXiv:0902.139910.48550/arXiv.0902.1399 (2009).10.48550/arXiv.0902.1399
43. Alsing PM Stephenson GJ Jr Kilian P Spin-induced non-geodesic motion, gyroscopic precession, Wigner rotation and EPR correlations of massive spin 1/2 particles in a gravitational field arXiv preprint arXiv:0902.1396 2009 10.48550/arXiv.0902.1396
Alsing, P. M., Stephenson, G. J. Jr. & Kilian, P. Spin-induced non-geodesic motion, gyroscopic precession, Wigner rotation and EPR correlations of massive spin 1/2 particles in a gravitational field. arXiv preprint arXiv:0902.139610.48550/arXiv.0902.1396 (2009).10.48550/arXiv.0902.1396
44. Misner CW Thorne KS Wheeler JA Gravitation 1973 Macmillan
Misner, C. W., Thorne, K. S. & Wheeler, J. A. Gravitation (Macmillan, 1973).
45. Birrell ND Davies PCW Quantum Fields in Curved Space 1984 Cambridge University Press
Birrell, N. D. & Davies, P. C. W. Quantum Fields in Curved Space (Cambridge University Press, 1984).
46. Hartle JB Gravity: An Introduction to Einstein’s General Relativity 2003 American Association of Physics Teachers
Hartle, J. B. Gravity: An Introduction to Einstein’s General Relativity (American Association of Physics Teachers, 2003).
47. Yeongkwang Kim. https://gofile.me/6VGLQ/0D1FbsC4f
48. Chandrasekhar S Thorne KS The Mathematical Theory of Black Holes 1985 American Association of Physics Teachers
Chandrasekhar, S. & Thorne, K. S. The Mathematical Theory of Black Holes (American Association of Physics Teachers, 1985).
49. Marck J-A Solution to the equations of parallel transport in Kerr geometry; tidal tensor Proc. R Soc. Lond. A 1983 385 431 438 10.1098/rspa.1983.0021
Marck, J.-A. Solution to the equations of parallel transport in Kerr geometry; tidal tensor. Proc. R Soc. Lond. A 385, 431–438. 10.1098/rspa.1983.0021 (1983).10.1098/rspa.1983.0021
50. Thorne KS Misner CW Wheeler JA Gravitation 2000 Freeman San Francisco
Thorne, K. S., Misner, C. W. & Wheeler, J. A. Gravitation (Freeman San Francisco, 2000).
51. Alsing PM Milburn GJ On entanglement and Lorentz transformations Quantum Info. Comput. 2002 2 487 512
Alsing, P. M. & Milburn, G. J. On entanglement and Lorentz transformations. Quantum Info. Comput. 2, 487–512 (2002).
52. Saldanha PL Vedral V Physical interpretation of the Wigner rotations and its implications for relativistic quantum information New J. Phys. 2012 10.1088/1367-2630/14/2/023041
Saldanha, P. L. & Vedral, V. Physical interpretation of the Wigner rotations and its implications for relativistic quantum information. New J. Phys.10.1088/1367-2630/14/2/023041 (2012).10.1088/1367-2630/14/2/023041
53. Baskal S Kim YS Wigner rotations in laser cavities Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 2002 66 026604 10.1103/PhysRevE.66.026604 12241308
Baskal, S. & Kim, Y. S. Wigner rotations in laser cavities. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 66, 026604. 10.1103/PhysRevE.66.026604 (2002).12241308 10.1103/PhysRevE.66.026604
54. Thomas LHI The kinematics of an electron with an axis Lond. Edinb. Dublin Philos. Mag. J. Sci. 1927 3 1 22 10.1080/14786440108564170
Thomas, L. H. I. The kinematics of an electron with an axis. Lond. Edinb. Dublin Philos. Mag. J. Sci. 3, 1–22. 10.1080/14786440108564170 (1927).10.1080/14786440108564170
