Towards Communication in a Curved Spacetime Geometry
Qasem Exirifard∗, Eric Culf, and Ebrahim Karimi∗ Department of Physics, University of Ottawa, 25 Templeton St., Ottawa, Ontario, K1N 6N5 Canada and ∗Corresponding authors: [email protected] and [email protected]
The current race in quantum communication – endeavouring to establish a global quantum network – must account for special and general relativistic effects. The well-studied general relativistic effects include Shapiro time-delay, gravitational lensing, and frame dragging which all are due to how a mass distribution alters geodesics. Here, we report how the curvature of spacetime geometry affects the propagation of information carriers along an arbitrary geodesic. An explicit expression for the distortion onto the carrier wavefunction in terms of the Riemann curvature is obtained. Furthermore, we investigate this distortion for anti de Sitter and Schwarzschild geometries. For instance, the spacetime curvature causes a 0.10 radian phase-shift for commu- nication between Earth and the International Space Station on a monochromatic laser beam and quadrupole astigmatism; can cause a 12.2 % cross-talk between structured modes traversing through the solar system. Our finding shows that this gravitational distortion is significant, and it needs to be either pre- or post-corrected at the sender or receiver to retrieve the information.
Introduction duce all the effects of a curved space-time geometry. In [19], all the multi-polar ` modes are presented only at the level of Photons, electromagnetic waves, are widely used in clas- the equations, however, the upper value of ` = 100 on the sical and quantum communication since they do not possess multi-polar modes is considered to compute the solution. On electric charge or rest mass. However, a photon’s traits, e.g., the surface of the earth, a narrow beam with an initial width of group and phase velocity, wavelength, linear and optical an- 10 cm and a large value of Rayleigh range requires taking into gular momentum, are modified inside or during propagation account the contribution of multi-polar modes up to at-least 9 through a linear or a nonlinear medium. These traits are gov- ` = 10 . So the solutions presented in ref. [19] do not take erned by Maxwell’s equations which are the relativistic quan- into consideration all the multi-polar modes required to cal- tum field theory of the U(1) gauge connection. Understanding culate the effects of the curvature at the vicinity of the Earth. how these optical properties are altered upon propagation is a Here, we present a computationally simple method to calcu- key element for any optical communication network. In opti- late the distortion by the curvature of any spacetime geometry cal communication, the sender and the receiver, namely Alice on any localised wavepacket. and Bob, use one or several internal photonic degrees of free- In flat spacetime geometry, the propagation of each polari- dom, such as wavelength, polarisation, transverse mode, or sation of photon is isomorphic to the propagation of a mass- time-bins, to share information, including a ciphertext and the less scalar field. Based on our understanding of the Einstein secret key to decrypt the ciphertext. The propagation, e.g., equivalence principle, we expect that studying a massless through fibre, air or underwater channels, causes these pho- scalar field theory also captures some features of a photon’s tonic degrees of freedom to be altered, and thus causes un- propagation in curved spacetime geometry. Therefore, in this desired errors on the shared information. Therefore, the al- study, we consider the propagation of a massless scalar in the teration to those degrees of freedom for any communication bulk of the paper. In Supplementary Note 2, we prove that in channel needs to be considered and well-examined. Sharing the Lorentz gauge, each linear polarisation of the photon in a information with a longer range or with moving objects, e.g. curved spacetime geometry gets corrected as if it were a mass- satellites, airplanes, submersibles [1–5], requires the optical less scalar field. It is reported that the Riemann tensor quan- beam to not only traverse through a medium, but in a few tum mechanically alters the wavepacket propagating along a cases, also in the fabric of the spacetime geometry, where gen- geodesic. The alteration operators depend on the geodesic and arXiv:2009.04217v2 [gr-qc] 1 Aug 2021 eral relativistic effects manifest [6–9]. Effects associated to components of the Riemann tensor on the geodesic. The al- the change of the geodesic due to a mass distribution, such as teration is calculated for examples including the space-time Shapiro time-delay [10], gravitational lensing [11], and frame geometry around the Earth and the Sun. dragging [12], are well-studied and observed [13–16]. We explore the propagation of relativistic wavepackets along an arbitrary null geodesic in a general curved space- Results and Discussion time geometry, and show how the curvature of the spacetime geometry distorts the wavepacket as it travels along the null We start by considering a relativistic massless scalar field geodesic. Different methods are used to tackle this study. For ψ := ψ(xµ) that propagates in a curved spacetime geometry instance, the four-dimensional Klein-Gordon equation is ap- xµ = (t, xi) = (t, x, y, z) with an arbitrary Riemann curvature proximated to a simple two-dimensional partial differential tensor Rabcd. The units are chosen such that the speed of light equation by ignoring all the multi-polar modes [17, 18]. In in vacuum and the reduced Planck constant are equal to one, particular, in derivation of eq.(7) from eq.(6) in [17], the sec- i.e. c = 1 and ~ = 1. Let us consider a localized wave function ond term on the left-hand side of eq. (41) of [19] has not been (information carrier) whose size is small compared to the cur- taken into account, therefore, [17, 18] cannot claim to repro- vature of the spacetime geometry. At the leading order, there- 2
a b c
t x+
Bob Alice x+ x γ x−
1 y x
z x2
FIG. 1: Schematics of communication in a general curved spacetime geometry and proper chosen coordinates. a Two parties, namely Alice and Bob, communicate in a general curved spacetime geometry. Alice encodes her message in a sequence of information carriers and sends them to Bob. The information traverses through the spacetime over a geodesic γ. The physical traits of information are distorted by the curvature of the space-time, causing errors to the field wavepacket. b The shown null geodesic γ. Close to γ, at the local coordinates, the metric is approximately pseudo-Euclidean. c The associated Fermi coordinates, the null-geodesic-path is mapped to x+.
fore, the carrier can be treated as a massless point-like particle tion of a point-like particle along the geodesic; the rest are that travels along a null geodesic γ, see Fig. 1 a. We choose quantum. The effective action of a free photon propagating in the local Fermi coordinates [21] along the geodesics in or- curved spacetime geometry coincides to the tree-level action, der to compute the quantum relativistic corrections. The met- therefore, we omit the subscript c. ric’s components in the Fermi coordinates can be expanded The massless scalar (quantum) field ψ obeys the (covari- − 1 1 µν in terms of the components of the Riemann tensor Rabcd and ant) wave equation, ψ = (−g) 2 ∂µ(−g) 2 g ∂νψ = 0. In the its covariant derivatives evaluated on the geodesic, see Fig. 1 Fermi coordinates, gµν can be viewed as a perturbation to the b,c. The expansion of the metric up to quadratic order in the Minkowski metric, inducing expansion series for the inverse transverse coordinates of the geodesic of a massless particle γ and determinant of the metric: gµν = ηµν + εδgµν + O(ε2) and p is given by [20], ln( − det g) = εδg + O(ε2). ε is the systematic perturbation parameter introduced to keep track of the perturbation series, 2 + − a b a¯ b¯ + 2 ds = 2dx dx + δabdx dx − R+a¯+b¯ x x (dx ) (1) which means that all the components of the Riemann tensor in Eq. (1) are multiplied with ε, and ε is treated as an infinitesi- 4 b¯ c¯ + a¯ 1 b¯ c¯ a¯ b¯ mal parameter. At the end of the computation, we set ε = 1. + R+b¯a¯c¯ x x (dx dx ) + Ra¯c¯b¯d¯x x (dx dx ) + ..., 3 3 This technique helps us to systematically perform perturba- √ tions for small curvatures. Utilizing the perturbation gives, where x± = (x3 ± t)/ 2 represent the Dirac light-cone coor- + (0) µν µν 2 dinates [22] in the Fermi coordinates (x is always tangent ψ = ψ + ε(∂µ(δg ∂νψ) + η ∂νδg∂µψ) = O(ε ), (2) to the null geodesic Fig. 1-b), δab is the Kronecker delta, ˆ (0) µν 2 γ˙ is the tangent of the null geodesic and a, b ∈ {1, 2}, and where = η ∂µ∂ν = 2∂−∂+ + ∇⊥ is the d’Alembert oper- a¯ − a 2 2 2 (x ) = (x , x ) and the curvature components are evaluated ator in the flat spacetime geometry, and ∇⊥ = ∂1 + ∂2. The on γ. The tree-level action of a massless scalar field ψ in perturbative nature of Eq. (2) seeks for a series expansion, (0) (1) 2 (0) a general√ curved spacetime geometry is given by S [ψ] = ψ = ψ + εψ + O(ε ). Here, ψ satisfies the scalar wave 1 R 4 µν (0) (0) 2 dx −gg ∂µψ∂νψ, where g is the determinant of the equation in the flat spacetime geometry ψ = 0, and the (1) metric gµν, and µ, ν ∈ {±, a}. Since the tree-level action is perturbed term to the wavefunction, ψ , yields, quadratic in terms of ψ, it is quantum mechanically exact, (0) (1) µν (0) µν (0) which can be verified by looking√ at its generating function, ψ = −∂µ(δg ∂νψ ) + η ∂νδg∂µψ . (3) R R 4 R i.e., Z[J] = Dψe−i(S [ψ]+ d x − det gJψ)/ Dψe−iS [ψ] – Dψ We assume the Fourier expansion in terms of x− variable, represents the integration over all field configurations and J (0) R (0) + a iωx− (0) is the source field. Its exact effective action, as defined by the ψ = dω fω (x , x )e , where fω satisfies the paraxial 2 (0) Legendre transformation of ln Z , coincides with the tree-level equation 2iω∂+ + ∇⊥ fω = 0. This implies that the so- action, i.e., Γ[ψc] = S [ψc]. For a general action, ϕc resembles lutions, given by the paraxial approximation in optics [23, a “classical” field whose action is given by Γ[ψc], while Γ[ψc] 24], are exact. The paraxial equation is isomorphic to the encapsulates all the quantum loop corrections. The exact ef- Schrodinger¨ equation, and its solutions (the transverse and fective action includes both the classical and quantum effects. longitudinal parts) can be expressed in the form of Laguerre- The classical effects are those that can be reproduced by mo- Gauss (LG) modes (with an azimuthally-symmetric intensity
2 Alice (0) (1) − a¯ b¯ 2 (0) ψ = R+a¯+b¯ x x ∂−ψ . (4) ra Bob (1) R (1) + a iωx− (1) Here, ψ = dω fω (x , x )e with fω being the correc- tion to the structure function for the frequency of ω. We have found the solutions to Eq. (4) – see the Supplementary Note FIG. 2: Propagation of wavepacket and associated null geodesics 1 for more detail on the derivation. The solution is, in the Schwarzschild spacetime geometry. a The null geodesics (dashed curves) for beams that propagate very close to the event hori- Z zon - we set the event horizon at 1. b Schematic of wavepacket prop- X − µ iωx ω ω ω (0) agation radially in the Schwarzschild spacetime geometry. Alice and ψ(x ) = dω cp`n(ω)e 1 + ε(O + Q + Q ) f , U N ω r r θ ϕ p,`,n Bob are located at a and b, respectively, while and are polar and azimuthal angles of the standard spherical coordinates. where cp,`,n(ω) are defined based on the initial and boundary ω ω ω conditions. The operators, O , QU and QN encodes how the curvature of the space-time geometry distorts the wavepacket. They are given by, de Sitter and anti de Sitter spacetime geometries: The de Sitter and anti de Sitter space-times are maximally symmetric iω i and the Riemann tensor at any given event in the space-time Oω = − G xa xb + G˜ xa∂b + G˜ ∂a∂b, 2 ab ab ω ab in any coordinates, including the Fermi coordinates, is given 2 ! i ω by R 0 0 = Λ(g 0 g 0 − g 0 g 0 ). The value of Λ determines Qω = 1 − x−2 G − ix−ωxaG + x−G˜ ∂a, µνµ ν µµ νν µν νµ U ω 2 −− −a −a different geometries: Λ > 0 represents the de Sitter space- 2i time geometry; Λ < 0 represents the anti de Sitter spacetime Qω = x−G + G xa + G˜ ∂a, (5) N −− −a ω −a geometry; and Λ = 0 is the Minkowki spacetime geometry. R+−+− = Λ is the only non-zero component for R+a¯+b¯ evalu- ˜ where Gab, G˜ab, and G˜ ab are integrals of the components of the ated on the geodesic, and thus the correction operators, Eq. (5) Riemann tensor R+a¯+b¯ evaluated on the geodesic – see Supple- are, mentary Note 1: 2 ! ω ω i ω −2 + ω + − Z τ Z τ Z τ O = 0, QU = 1 − x Λx , QN = Λx x . (7) ˜ ˜ ˜ ω 2 Ga¯b¯ = dτR+a¯+b¯ , Ga¯b¯ = dτGa¯b¯ , Ga¯b¯ = dτGa¯b¯ , 0 0 0 (6) Let us consider a Gaussian wavepacket with normal distri- where τ is the affine parameter on the geodesic. The distor- bution for ω around ω0 with the width of σ, i.e. ψAlice = − 2 tions provided by (5)is the solution to the exact quantum effec- (0) + 1 2 iω x− − (σx ) f (x , x , x )e 0 e 2 . The validity of the perturbative tive action and cannot be reproduced by motion of a point-like 2 2 solution demands that |Λ| ω0, |Λ| σ and σ ω0. particle along a geodesic. They do not exist in flat spactime ω0 The wavepacket after the propagation is ψBob ' (1 + εQU + − geometry, so they manifest a set of quantum effects in curved ω − (x )2σ2 0 iω0 x (0) − 2 spacetime geometry. A similar approach can be used to find εQN )e f e . ω0 ω0 the wavefunction of a massive scalar particle. The physical QN and QU change the wavepacket amplitude and phase, ω0 − degrees of U(1) gauge fields get corrected as if they were respectively. The maximum of |QN ψAlice| occurs at x = scalar fields, see Supplementary Note 2 and Supplementary ±1/σ. Requiring√ it to be smaller than 1 yields T < τA where Note 3. Supplementary Note 4 presents the operators in the τA = (σ e)/|Λ|. τA is the maximum time that the wavepacket Hilbert space that corresponds to these corrections. The fol- feels the curvature of the spacetime geometry and keeps its amplitude intact. Qω0 alters the phase of the wavepacket. lowing subsections show how these operators distort the phys- U √ ω0 − ical information encoded in wave-packets traveling along cou- The maximum of |QU ψAlice| occurs at x = ± 2/σ. Re- ple of examples of null geodesics in the Solar system and quiring√ it to be smaller than 1 results in T < τϕ where τϕ = around the Earth. (σ τA e)/(2ω0). τϕ is the maximum time that a wavepacket We now study the distortion operators in a couple of space- can feel the curvature of the spacetime geometry and keep its time geometries, including de Sitter and Schwarzschild space- phase intact. For T > τϕ, information about phase is lost at time geometries. We first study the de Sitter and anti de Sitter perturbation. This may point to a “gravitational decoherence”, space-time geometries because their symmetry allows one to and its possible consequence on anti-de Sitter/conformal field immediately write down the components of the Riemann ten- theory correspondence [26] demands attention. τA represents sor in Fermi coordinates evaluated on the geodesic. the amount of time of interaction with the curvature that the
a b wavepacket can keep its amplitude intact. We observe that 1.0 1. 0.6 0.6 1
0.8 0.8 − 1 ◆ - − 2 b )
τ τ . So the phase changes sooner than the change in the 1 r
ϕ A ◆ )
) 0.4 0.4 ) 2 b 1 ) r −
0.6 0.6 ) - 2 - a − 1 r a.u. a.u. 2 ( ( a ( ( 1 amplitude. ✓ ( 0.2 0.2 r ( 2 χ 0 ✓ A 0.4 0.4 σ ω | 4 A | δ 0.0 0
Schwarzschild spacetime geometry: | We choose the stan- δχ 0.2 0.2 | dard spherical coordinates r, θ, ϕ where geodesics are extrema -0.2 -0.2 0.0 0 -4 -2 0 2 4 -4 -2 0 2 4 of, σ x-− σ x-− m m−1 L = − 1 − t˙2 + 1 − r˙2 + r2(θ˙2 + sin2 θϕ˙2) , (8) FIG. 3: Propagation of Gaussian wavepacket along radial direc- r r tion in the Schwarzschild spacetime geometry. a Amplitude (pur- ple) and phase (orange) of the initial Gaussian beam at the sender and m = 2GM• is the Schwarzschild radius, M• is the mass of (Alice) respectively presented by A(a.u.) and χ(a.u.). Their units, the blackhole, and G is the gravitational constant – the units a.u., are chosen such that both are normalized to their maximum val- are such that m = 1. The components of the Riemann tensor ues. b The changes in the amplitude (δA) and the phase (δχ) of the in the Fermi-coordinates adapted to a general null geodesic Gaussian wavepacket around frequency ω0 with width of σ trans- of Schwarzschild spacetime geometry are derived in the Sup- mitted over the radial null geodesic in the Schwarzschild geometry plementary Note 5. Figure 2-a shows several null geodesics between ra and r = rb, the geometry is shown in Fig. 2-b. The units that go very close to a blackhole. We first consider that the are chosen where the Schwarzschild radius and light speed are one. wavepacket propagates along the radial direction, Figure 2- b, where the only non-zero components of the Riemann ten- − 3 sor is R+−+− = 1/r . This is the same component that used to extend the communication alphabet beyond bits, i.e., appeared in the de Sitter spacetime geometry. The radial ω 0 and 1 [27]. The longitudinal and frequency distributions are geodesic has l = 0, and its correction operators are O = 0, σ − 2 assumed to be Gaussian. When is large, the dominant cor- ω i (ωx ) 1 1 ω 1 1 1 − Q = − 1 − 2 − 2 , and Q = − 2 − 2 x , rection operator is calculated to be: U 2ω 2 ra r N 2 ra r where Alice is located at ra. These correction terms do not 2 9.5ia 2 2 contain derivatives of the spatial transverse coordinates. Thus, O = − x − y , (10) zR the Riemann tensor does not affect the spatial transverse pro- √ √ + + file of the wavepacket. This is due to the symmetry, as the where x = 2x1/w(x ) and y = 2x2/w(x ) are dimen- q 2 radial geodesic inherits the static and spherical symmetry of + x+ sionless coordinates, w(x ) = w0 1 + is the beam ra- the background. zR 1 2 Figure 3-a shows the amplitude and phase of a Gaussian dius, zR = 2 ω0w0 is the Rayleigh range, a = ra/l is a time-bin signal. Figure 3-b depicts the alterations in the am- scaling parameter, w0 is the beam radius at Alice’s position plitude and the phase; - see Supplementary Note 6 for more details. This opera- tor contains coordinate parameters x and y, and thus, alters − 2 (σx−) both the amplitude and phase of the transverse modes upon εx 1 1 − 2 δA = − − e 2 + O ε (9) 2 2 propagation. The correction for the solar system, when Al- 2 ra rb ice and Bob are at the mean Earth-Sun distance from the (σx−)2 iω0 1 1 −2 − 2 Sun and the wavepacket passes at l = 2R , and for z = δχ = − σx e 2 + O ε , R 2 2 2 2.8 km × a2 = 1.34 × 109 m, remains perturbative and is given 4σ ra rb (1) µ 2 2 µ by εψ (x ) = i 0.10 x − y ψAlice(x )|x+=T . where rb is the location of Bob. The maximum alteration to The amplitude and phase of the mode ψ(0)(xµ) and the the amplitude and phase of√ the Gaussian wavepacket occur at correction εψ(1)(xµ) for a few Hermite-Gaussian modes are − − x = ±1/σ and x = ± 2/σ, respectively. For a Gaussian shown in Fig. 4. As seen, these alterations on the modes are wavepacket that propagates radially close to the Earth, the considerable. For instance, it causes up to 12.2% crosstalk be- maximum alteration to the amplitude and phase are respec- tween HG0,3 and HG0,1 modes. The crosstalk between mode | A | 1 1 − 1 ω0 m⊕ 1 − 1 M N |hN|Mi|2 tively δ max. = 2σ r2 r2 and δχmax. = 2σ2 r2 r2 , and mode is given by . Therefore, these alter- a b a b ations need to be accounted for when information is encrypted where m⊕ is the Schwarzschild radius of Earth. in the spatial modes. The action of curved spacetime geome- For ν0 = 456 THz and ∆ν = 1 kHz, δχmax. = 0.10 rad. Supplementary Note 8 provides further details on choosing try on the wavepacket is linear. Therefore, a target beam that these values. does not possess the information can be used as a reference As a final example, we examine the weak regime of grav- to monitor the distortion of the information carrier, and an ac- ity when the beam possesses well-defined transverse modes tive system can be employed for compensating the distortion – see Fig. 2-(a). The wavepacket carrying a well-defined in real-time – in conjunction, results in retrieving the origi- transverse mode traverses through space and reaches to the nal information. Moreover, whenever ε ' 1, the higher-order terms of correction need to be considered. For instance, for minimum distance of l to the central mass – here, we as- 2 sume l is large. We now consider a specific wavepacket, zR < 28 km × a , the correction becomes larger than 1 and we (0) − + 1 2 need to take into account higher ε terms. Taking into account a Hermite-Gauss transverse mode fω,p,`,q(x , x , x , x ) = − 2 all the corrections is tantamount to knowing the Riemann ten- iω x− − (σx ) + 1 2 e 0 e 2 HGm,n(x , x , x ) – Hermite-Gauss modes are sor in whole of the spacetime geometry, a piece of knowledge
0 0 0 0
+ 1 2 FIG. 4: Propagation of Hermite-Gaussian modes in a weak gravitational field. Hermite-Gaussian mode HGm,n(x , x , x ) is not shape invariant under propagation in the Schwarzschild spacetime geometry, and both the intensity and phase profiles modify. a Intensity and phase distributions of the first 9 HG modes m, n ∈ {0, 1, 2} prior to the free-space propagation. b The intensity |εψ(1)|2 and phase arg εψ(1) distributions of the corrected term, for the first 9 HG modes m, n ∈ {0, 1, 2}, after the propagation in a weak gravitational field. The rows and columns are associated to m ∈ {0, 1, 2} and n ∈ {0, 1, 2}, respectively. The dimensionless coordinates are used for these plots.
which is not attainable. Thus, we tend to argue that once the places, respectively. perturbation breaks, in addition to known well-studied gravi- tational decoherence [28, 28–32], a decoherence occurs. We observe that, in addition to the known decoherence of a bipar- tite entangled system when each particle traverses through a Conclusion different gravitational field gradient [32, 33], a coherent beam decoheres when different segments of the spatial spread of the We have presented how the curvature of the spacetime ge- wave experience different tidal gravitational field gradients. ometry affects the propagation of an arbitrary wavepacket The phenomenon we are reporting also occurs for geodesics along a general geodesic in a general curved spacetime ge- passing very close to the event horizon – see Supplementary ometry. The effect is beyond classical general relativity, re- Note 7. siding in the same category as Hawking radiation [34]. A Finally, it is noteworthy that photon pairs |ψientangled, e.g. set of linear operators are presented that encode the effect entangled in spatial, frequency or temporal modes, would be of the curvature. The corrections to the information carrier affected by the curved spacetime geometry whenever they are wavepacket are investigated in cases of de Sitter (anti de Sit- shared between two parties, namely Alice and Bob. The final ter) and Schwarzschild spacetime geometries. It has been state of the entangled photon, indeed, is given by applying shown that the corrections accumulate over time and distort the non-local operators U = (1 + ε(O + QU + QN )) onto the the wavepacket. The gravitational distortion, therefore, needs entangled states, (UA ⊗ UB) |ψientangled - here, UA and UB are to be accounted for in quantum communication performed associated with the correction operators at Alice and Bob’s over long distances in a curved spacetime geometry.
References [8] Chou, C.-W., Hume, D. B., Rosenband, T. & Wineland, D. J. Optical clocks and relativity. Science 329, 1630–1633 (2010). [1] Ursin, R. et al. Entanglement-based quantum communication [9] Ashby, N. Relativity in the Global Positioning System. Living over 144 km. Nature Physics 3, 481–486 (2007). Reviews in Relativity 6, 1 (2003). [2] Yin, J. et al. Satellite-based entanglement distribution over [10] Shapiro, I. I. et al. Fourth test of general relativity: Preliminary 1200 kilometers. Science 356, 1140–1144 (2017). results. Physical Review Letters 20, 1265–1269 (1968). [3] Yin, J. et al. Satellite-to-ground entanglement-based quantum [11] Perlick, V. Gravitational Lensing from a Spacetime Perspective. key distribution. Physical Review Letters 119, 200501 (2017). Living Reviews in Relativity 7, 9 (2004). [4] Wengerowsky, S., et. al. 2020. Passively stable distribution [12] Everitt, C. et al. Gravity Probe B: Final Results of a Space of polarisation entanglement over 192 km of deployed optical Experiment to Test General Relativity. Physical Review Letters fibre. Quantum Information 6, (2020) 106, 221101 (2011). [5] Hufnagel, F., Sit, A., Bouchard, F., Zhang, Y., England, D., [13] Will, C. The Confrontation between General Relativity and Ex- Heshami, K., Sussman, B. J., and Karimi, E. Investigation of periment. Living Reviews in Relativity 17, 4 (2011). underwater quantum channels in a 30 meter flume tank using [14] Liu, T., Wu, S. & Cao, S. The influence of the Earth’s curved structured photons. New Journal of Physics,22(2020) spacetime on Gaussian quantum coherence. Laser Physics Let- [6] Einstein, A. Uber¨ den einfluß der schwerkraft auf die ausbre- ters 16 (2019). itung des lichtes. Annalen der Physik 340, 898–908 (1911). [15] Kish, S. P. & Ralph, T. C. Quantum Metrology in the Kerr [7] Pound, R. V. & Rebka, G. A. Apparent weight of photons. Metric. Physical Review D 99, 124015 (2019). Physical Review Letters 4, 337–341 (1960). [16] Pierini, R. Effects of gravity on continuous-variable quantum
key distribution. Physical Review D 98, 125007 (2018). [33] Ralph, T. C., Milburn, G. J., Downes, T. Quantum connectiv- [17] Bruschi, D. E., Ralph, T., Fuentes, I., Jennewein, T. & Razavi, ity of space-time and gravitationally induced decorrelation of M. Spacetime effects on satellite-based quantum communica- entanglement. Physical Review A 79, (2009). tions. Physical Review D 90, 045041 (2014). [34] Hawking, S. W. Black hole explosions? Nature (London) 248, [18] Bruschi, D.E., Datta, A., Ursin, R., Ralph, T. C. & Fuentes, I. 30–31 (1974). Quantum estimation of the Schwarzschild spacetime parame- [35] Manasse, F. K. and Misner, C. W., Fermi Normal Coordinates ters of the Earth. Physical Review D 90, 124001 (2014). and Some Basic Concepts in Differential Geometry, Journal of [19] Jonsson, R. H.,Aruquipa, D. Q., Casals, M.,Kempf,A. and Mathematical Physics 4, 735 (1963), Mart´ın-Mart´ınez, E., Communication through quantum fields [36] Henry, R. C. Kretschmann scalar for a kerr-newman black hole. near a black hole, Physical Review D 101, 125005 (2020). [ Astrophysical Journal 535, 350 (2000). [20] Blau, M., Frank, D. & Weiss, S. Fermi coordinates and Penrose [37] Goto, T. On the Longitudinal and Scalar Photons in Lorentz limits. Classical and Quantum Gravity 23, 3993–4010 (2006). Gauge and Lorentz Condition. Progress of Theoretical Physics, [21] Fermi, E. Sopra i fenomeni che avvengono in vinicinanza di 37 (1967). una linea oraria. Atti R. Accad. Lincei Rend. Cl. Sci. Fis. Mat. Acknowledgments This work was supported by the High Nat. 31, 21–51 (1922). [22] Dirac, P. A. M. Forms of relativistic dynamics. Review of Mod- Throughput and Secure Networks Challenge Program at the ern Physics 21, 392–399 (1949). National Research Council of Canada, the Canada Research [23] Siegman, A. E. Lasers university science books. Mill Valley, Chairs (CRC) and Canada First Research Excellence Fund CA 37, 462–466 (1986). (CFREF) Program, and Joint Centre for Extreme Photonics [24] Bialynicki-Birula, I. & Bialynicka-Birula, Z. Beams of electro- (JCEP). The authors would like to thank Alicia Sit, Benjamin magnetic radiation carrying angular momentum: The riemann– Sussman, Khabat Heshami, Thomas Jennewein, Christoph Si- silberstein vector and the classical–quantum correspondence. mon, Mathias Blau, Ida Zadeh and Loriano Bonora for fruit- Optics Communications 264, 342–351 (2006). ful discussions and thoughtful feedback, and Haorong Wu for [25] Bliokh, K. Y., Bliokh, Y. P., Savel’Ev, S. & Nori, F. Semiclassi- the email correspondence. Q.E. would like to thank Fernando cal dynamics of electron wavepacket states with phase vortices. Quevedo and Atish Dabholkar for the nice hospitality in ICTP Physical Review Letters 99, 190404 (2007). where part of the work was conducted. [26] Maldacena, J. M. The Large N limit of superconformal field theories and supergravity. International Journal of Theoretical Physics 38, 1113–1133 (1999). Author contributions Q.E. and E.K. conceived the idea [27] Krenn, M., Fickler, R., Fink, M., Handsteiner, J., Malik, M., and developed the theoretical framework. E.C. checked and Scheidl, T., Ursin, R. & Zeilinger, A. Communication with confirmed the computation. All authors contributed to the spatially modulated light through turbulent air across Vienna. manuscript preparation. New Journal of Physics 16, 113028 (2014). [28] Pang, B. H., Chen, Y. & Khalili, F. Y. Universal Decoherence Data availability under Gravity: A Perspective through the Equivalence Princi- The authors declare that the data supporting ple. Physical Review Letters 117, 090401 (2016). the findings of this study are available within the paper and its [29] Stefanov, V., Siutsou, I. & Mogilevtsev, D. Gravitational de- supplementary information file. phasing in spontaneous emission of atomic ensembles in timed Dicke states. Physical Review D 101, 044042 (2020). Competing interests The authors declare no competing inter- [30] Bassi, A., Großardt, A. & Ulbricht, H. Gravitational Decoher- ests. ence. Classical Quantum Gravity 34, 193002 (2017). [31] Penrose, R. On gravity’s role in quantum state reduction. Gen- eral Relativity and Gravitation 28, 581–600 (1996). Author Information Correspondence and requests for ma- [32] S. K. Joshi, et al. [Space QUEST topical Team], Space QUEST terials should be addressed to Q.E. ([email protected]) or mission proposal: Experimentally testing decoherence due to E.K. ([email protected]). gravity New Journal of Physics 20, Supplementary Information for Towards Communication in a Curved Spacetime Geometry
Qasem Exirifard, Eric Culf, and Ebrahim Karimi Department of Physics, University of Ottawa, 25 Templeton St., Ottawa, Ontario, K1N 6N5 Canada; [email protected], [email protected]
Supplementary Note 1: Leading order correction term
We consider the Fourier integral representation of ψ(1) in terms of x−: Z (1) (1) + a iωx− ψ = dω fω (x , x )e , (S1)
(1) (0) R (0) + a iωx− where fω is the correction to the structure function for the frequency ω. Substituting (S1) and ψ = dω fω (x , x )e in (0) (1) 2 a¯ b¯ (0) ψ = +ω R+a¯+b¯ x x ψ yields, Z Z iωx− 2 (1) 2 a¯ b¯ (0) iωx− dωe 2iω∂+ + ∇⊥ fω = dωω R+a¯+b¯ x x fω e . (S2)
− iωx− iωx− We can use x e = −i∂ωe , and integration by parts to simplify the r.h.s of the above equation to, Z h 2 (0)00 2 (0)0 a 2 (0) a b i iωx− dω − ω fω R+−+− + 2i ω fω x R+−+a + ω fω x x R+a+b e , (S3) and then obtain,
2 (1) 2 (0) 00 2 (0) 0 a 2 (0) a b 2iω∂+ + ∇⊥ fω = −(ω fω ) R+−+− + 2i(ω fω ) x R+−+a + ω fω x x R+a+b , (S4) where ()0 stands for the derivatives w.r.t. ω. We would like to solve (S4) for all f (0) that solves
2 (0) (∇⊥ + 2iω∂+) f = 0. (S5)
(0) So fω can be expanded in the Hermite-Gaussian basis or the Laguerre-Gaussian basis. In the Hermite-Gaussian basis, one can write
(0) + a X + a fω (x , x ) = cln(ω) fln,ω(x , x ), (S6) l,n where cln(ω) is the coefficient of the expansion in the Hermite-Gaussian basis. This implies
(0) X ∂ω f = fln,ω∂ωcln + cln∂ω fln,ω . (S7) l,n We consider field configurations wherein
| fln,ω∂ωcln| |cln∂ω fln,ω|. (S8) which means that the frequency dependency of the structure function is mainly given by the coefficients of cln. This means that we are not considering white noise, and we are considering a wavepacket that has a sufficiently sharp Gaussian peak in the frequency spectrum. We are considering configuration that have sharp peak around the mean frequency. Because of (S8), we can approximate (S7) to
0 X ∂ω f = fln,ω∂ωcln, (S9) l,n which can be utilized to write
2 0 X 2 (∇⊥ + 2iω∂+)∂ω f = ∂ωcln(∇⊥ + 2iω∂+) fln,ω. l,n 3
Due to (S5) and (S6), it holds,
2 (∇⊥ + 2iω∂+) fln,ω = 0 Therefore,
2 0 (∇⊥ + 2iω∂+)∂ω f = 0. (S10) Utilizing the same method, eq. (S8), implies that all derivatives of f (0) with respect to ω solve (S5). The right-hand-side of (S4), therefore, includes terms that vanish by the same operator which appears in the left-hand-side of (S4). This resembles resonance where the source oscillates with the frequency of the system. In the following we, show that, similar to the resonance, the corrections accumulate over time. To solve (S4), along the geodesic γ withγ ˙ µ = xˆ+, we define:
Z x+ + + Ga¯b¯ = dξ R+a¯+b¯ (ξ ) , (S11a) 0 Z x+ ˜ + + Ga¯b¯ = dξ Ga¯b¯ (ξ ) , (S11b) 0 Z x+ ˜ + ˜ + Ga¯b¯ = dξ Ga¯b¯ (ξ ) , (S11c) 0 where ξ+ is the affine parameter on the geodesic. Next write f (1) by: 00 0 2 (0) 2 (0) ω fω ω fω iω f (1) = f˜(1) − G + xaG − G xa xb f (0) , (S12) ω ω 2iω −− ω −a 2 ab ω where f˜(1) is an arbitrary function of x+ and xa. Utilizing (S12) in (S4) yields:
iω 2 0 2iω∂ + ∇2 f˜(1) = 2iωG xa∂b f (0) + G f (0)∇2 xa xb − G ∂a ω2 f (0) (S13) + ⊥ ω ab ω 2 ab ω ⊥ ω −a ω We notice that the ++ component of the Ricci tensor along the geodesic holds:
µν ab η R+µ+ν = R+a+bη = 0 , (S14)
ab which implies δ Gab = 0 and simplifies (S13) to, 2 0 2iω∂ + ∇2 f˜(1) = 2iωG xa∂b f (0) − G ∂a ω2 f (0) . (S15) + ⊥ ω ab ω ω −a ω Now let i 0 f˜(1) = f˜(1) + G˜ xa∂b f (0) + G˜ ∂a ω2 f (0) , (S16) ω ω ab ω ω2 −a ω
˜(1) + a where f˜ω is an arbitrary function of x and x . Utilizing (S16) in (S15) yields: 2 ˜(1) ˜ a b (0) 2iω∂+ + ∇⊥ fω = −2Gab∂ ∂ fω , (S17) where we have noticed that derivatives of solutions to homogeneous differential equations are also solutions. The solution for f˜(1) is given by, i f˜(1) = G˜ ∂a∂b f (0) . (S18) ω ω ab ω Combining (S12),(S16) and (S18) gives: ! iω i i 00 1 0 i 0 f (1) = − G xa xb + G˜ xa∂b + G˜ ∂a∂b f (0) + G ω2 f (0) + xaG ω2 f (0) + G˜ ∂a ω2 f (0) (S19) ω 2 ab ab ω ab ω 2ω −− ω ω −a ω ω2 −a ω which we also directly have checked that solves (S4). Let us define: iω i Oω = − G xa xb + G˜ xa∂b + G˜ ∂a∂b , (S20) 2 ab ab ω ab i 00 1 0 i 0 Q[ f (0)] = G ω2 f (0) + xaG ω2 f (0) + G˜ ∂a ω2 f (0) . (S21) ω 2ω −− ω ω −a ω ω2 −a ω 4
Therefore,
ψ = ψ(0) + εψ(1) + O(ε2) (S22) Z iωx− (0) ω (0) (0) 2 = dωe fω + εO fω + εQ[ fω ] + O(ε ) ,
Integration by parts can be employed to rewrite (S22) to: Z µ iωx− ω ω ω (0) 2 ψ(x ) = dω e 1 + εO + εQU + εQN fω + O(ε ) , (S23) where ! i ω2 Qω = 1 − x−2 G − ix−ωxaG + x−G˜ ∂a , (S24) U ω 2 −− −a −a 2i Qω = x−G + G xa + G˜ ∂a . (S25) N −− −a ω −a
Supplementary Note 2: Corrections to U(1) gauge field
We would like to compute the correction due to the curvature of the spacetime geometry to an arbitrary structured photon propagating along a general null geodesic in a given curved spacetime geometry. We start from the action of a U(1) gauge field in a general curved spacetime endowed by metric gµν which is given by, Z 1 4 µµ0νν0 S [A ] = − d x g F F 0 0 , (S26) µ 4 µν µ ν 0 0 1 p 0 0 0 0 gµµ νν = − det g(gµµ gνν − gµν gνµ ) , (S27) 2 where Fµν is the field strength of the U(1) connection: Fµν = ∂µAν − ∂νAµ. The functional variation of the action with respect to the gauge field gives the its equation of motion:
µµ0νν0 ∂µ g Fµ0ν0 = 0 . (S28)
We would like to study perturbative solutions for the ε expansion of the metric. To this aim, we write an expansion series for the 0 0 gauge field and gµµ νν :
(0) (1) 2 Aµ = Aµ + εAµ + O(ε ) , (S29) 0 0 0 0 0 0 gµµ νν = g(0)µµ νν + εg(1)µµ νν + O(ε2) . (S30)
Utilizing (S29) and (S30) in (S28) and keeping the zero order yields
(0)µµ0νν0 (0) ∂µ g Fµ0ν0 = 0 → (0) (0)ν ηµν − ∂µ∂ν A = 0 , (S31)
µν (0)ν νλ (0) where from here forward η is utilized to move up or down the indices: A = η Aλ . We choose the Lorentz gauge, (0)ν (0) ∂νA = 0, that simplifies the equation for A to,
(0) (0) 2 (0) Aµ = 2∂+∂− + ∇⊥ Aµ = 0 . (S32) We are interested in the structured photon whose field configuration is given by,
(0) (0) + a iωx− Aµ = fµ (x , x )e , (S33)
(0) (0) (0) and the structure function, f , is very slowly varying: |∂µ fν | ω| fν |. The Lorentz gauge then implies:
(0) f+ = 0 , (S34) (0) a (0) ∂+ f− + ∂ fa = 0 , (S35) 5
(0) (0) that we choose to solve for f− . This leaves fa as the physical modes and maps (S32) to,
(0) (0) 2 (0) Aa = (2iω∂+ + ∇⊥)Aa = 0 . (S36)
2 (0) (0) iωx− In the plane wave approximation where ∇⊥Aa vanishes, it holds, Aa = Cae . Here, Ca are constants encoding the linear polarisation of the photon – this is in accordance with the presentation used in [37], where A1 and A2 are the transverse modes. They can be understood as the distribution of the polarisation of the photon. Utilizing (S29) and (S30) in (S28) and keeping the linear terms in ε yield
(0)µµ0νν0 (1) (1)µµ0νν0 (0) ∂µ g Fµ0ν0 = −∂µ g Fµ0ν0 → (0) (1)µ (1)µµ0νν0 (0) A = −∂ν g Fµ0ν0 . (S37) where it is noted that the zero term in ε vanishes due (S31), and the Lorentz gauge is used too. This the equation of motion for (1) (1)µν (1)µν Aµ . We are interested in the solutions where |∂λg | ω|g |. Therefore, ∂− in the rhs of (S37) does not act on the metric’s components:
(0) (1)µ (1)µµ0νν0 (0) A = −g ∂νFµ0ν0 , (S38)
As the beam’s shape is changing very slowly compared to ω, we should keep only ∂− in the rhs of (S38):
(0) (1)µ 2 (1)−−µα (1)−αµ− (0) A = −ω g − g Aα , (S39)
2 (0) − ω2 (0) wherein ∂−Aα = c2 Aα is utilized. Expressing (S39) in term of the components of the metric and utilizing (S34) returns:
(0) (1) A+ = 0 , (S40) (0) (1) −− (0) Aa = δg Aa , (S41)
(1) (1) (1) the first of which is solved by A+ = 0. Recalling the Lorentz gauge, we note that A− can be expressed in term of Aa so we have spared writing its equation. Utilizing the metric expansion series then gives,
(0) (1) 2 a¯ b¯ (0) Ai = ω R+a¯+b¯ x x Ai , (S42) which is identical to the equation for a scalar field. So its solution can be straight fully derived from (S22). For the zero-order solution of, Z (0) (0) + a iωx− Ai = dωAi (ω, x , x )e + c.c. (S43)
The corrected solution is given by, Z iωx− (0) + a 2 Ai = dωe (1 + εO + εQU + εQN ) Ai (ω, x , x ) + O(ε ) , (S44) where O and Q are defined in (S20) and (S21). This concludes the correction for a U(1) field.
Supplementary Note 3: Corrections to massive scalar field
Alice and Bob can communicate by exchanging wavepackets made of massive particles. We use the Fermi Coordinates adapted to the time-like geodesic, given by Manasse and Misner [35]: ! 2 1 ds2 = +c2dt2 −1 + R xl xm + R xl xmdtdxi + dxidx j δ + R xl xm + O(xl xm xn) , (S45) 0l0m 3 0lim i j 3 il jm where Rµαβν represents the components of the Riemann tensor computed along the time-like geodesics of the packet. The components of the Riemann tensor in (S45) are functions of time, and time is the proper time as measured in the rest frame of the massive wavepacket. We use the Greek alphabet µ from {0, 1, 2, 3} where 0 stands for time. We use Latin indices to show the directions. The generally covariant form of a wave equation of a massive scalar particle ϕ with rest mass of m0 is,
2 µν µν p 2 − m0 ϕ = ∂µ (g ∂νϕ) + ∂νϕg ∂µ ln( − det g) − m0ϕ = 0 . (S46) 6
Note that we set the units such that c = 1 and ~ = 1. Utilizing the ε expansion series for the metric then yields,
(0) 2 µν µν 2 − m0 ϕ = −ε∂µ(δg ∂νϕ) − εη ∂νϕ∂µδg + O(ε ) , (S47)
(0) µν 2 2 where = η ∂µ∂ν = −∂t + |∇| . We would like to solve (S47) in a perturbative fashion. In so doing we consider:
ϕ = f (0)eiωt + ε f (1)eiωt + O(ε2) , (S48) where ω := m0 is understood. Next assume that |∂µ f | ω| f | which gives:
2 ! ~ 2 (0) − ∇ − i~∂t f = 0 , (S49) 2m0 2 ! ~ 2 (1) 2 l m (0) − ∇ − i~∂t f = ω R0l0m(t)x x f . (S50) 2m0 Eq. (S49) is the Schrodinger¨ equation in 3 + 1 dimensions. Notice that (S50) can be directly obtained from considering the (1) (1) Schrodinger¨ equation in the approximation of Newtonian potential associated to g00 . The equation for f can be mapped to a generalization of (S4) to five dimensions where the components of the Riemann tensor in 5D vanishes for the x− transfer direction, and t = x+. This mapping gives: iω i f (1) = − G xa xb + G˜ xa∂b + G˜ ∂a∂b f (0), (S51) 2 ab ab ω ab where
Z T Gab = dτ R0a0b(τ) , (S52) 0 Z T G˜ab = dτ Gab(τ) , (S53) 0 Z T ˜ G˜ ab = dτ G˜ab(τ) , (S54) 0 where τ is an affine parameter and a, b ∈ {1, 2, 3}, and Bob gets the message T seconds (T in the rest frame of the particle) after a Alice writes it. We directly have checked that (S51) solves (S50) where we notice R0a0 = 0. We again notice that if we represent,
(0) iωt ϕAlice = f e , (S55) (0) (1) iωt ϕBob = f + ε f e , (S56) then,
ϕBob = O ϕAlice , (S57) where the spacetime operator O depends on the spacetime geometry and the geodesic: iω i O = 1 − G xa xb + G˜ xa∂b + G˜ ∂a∂b . (S58) 2 ab ab ω ab This concludes the correction for a massive scalar.
Supplementary Note 4: Corrections in the Hilbert space
We consider the Hilbert space for the scalar field ψ. We write an ε expansion series for the state |ψi: E E |ψi = ψ(0) + ε ψ(1) + O(ε2) (S59)
Let us consider: D E D E hx|ψi = x ψ(0) + ε x ψ(1) + O(ε2) = ψ(0) + εψ(1) + O(ε2) , (S60) 7 where ψ(0) is, Z (0) (0) + a iωx− ψ = dω fω (x , x )e , (S61) and ψ(1) is presented in (S1): Z (1) (1) + a iωx− ψ = dω fω (x , x )e . (S62)
(0) We notice that the equation that governs the dynamics of fω (S5), which is also the paraxial approximation of the Helmholtz equation, is isomorphic to the Schrodinger¨ equation. In fact me f f = ~ω maps (S5) to the Schrodinger¨ equation for a free particle + a of mass me f f in 2 + 1 dimensions of (x , x ):
2 ~ 2 (0) (0) − ∇⊥ fω = i~∂x+ fω . (S63) 2meff This isomorphism implies that x+ in the Hilbert space should be treated as a time-variable rather than an operator. The Hamilto- nian associated to x+ is:
2 P 1 2 2 H = = (P1 + P2) , (S64) 2meff 2ω
1 2 where P1 and P2 are the momentum operators in the transverse directions x and x . We choose a basis of the Hamiltonian to (0) represent fω (x):
Z + ! − X iE x ψ(0) = dωeiωx C(0) exp − s f (x1, x2) , (S65) s,ω ~ s s where fs is an eigenstate of the Hamiltonian with energy of Es. This allows us to write,
Z + !† E − X iEs x ψ(0) = dωeiωx C(0) exp − |0i , (S66) s,ω ~ s ω,s
† where aω,s is the creation operator for a mode with frequency ω and structure s. |0i is the vacuum state in the Fermi coordinates for ε = 0 (adapted to the null geodesic). Eq. (S23) presents the value for (S60):
Z + ! − X iE x hx|ψi = dω eiωx (1 + εOω + εQω + εQω ) C(0) exp − s f + O(ε2) , (S67) U N s,ω ~ s s
a a iPa In order to write |ψi, we first write the quantum analog of (S20),(S24) and (S25) by changing x → X and ∂a → ~ :
iω i G˜ Oω = − G XaXb + G˜ (XaPb + PbXa) − ab PaPb , (S68) 2 ab 2~ ab ~2ω ! i ω2 ix− Qω = 1 − (x−)2 G − ix−ωG Xa + G˜ Pa , (S69) U ω 2 −− −a ~ −a 2 Qω = x−G + G Xa − G˜ Pa , (S70) N −− −a ~ω −a
a ω ω ω + ˜ ˜ + where X are operators. Notice that O , QU and QN are functions of x because G, G and G are function of x . The corrected state (S59) then follows from (S23):
Z + X iEs x − s − ~ iωx ω ω ω † 2 |ψi = dωC (ω)e e 1 + ε(O + QU + QN ) as (ω) |0i + O(ε ) . (S71) s
+ ω ω ω + Alice is at x = 0. Notice that O , QU and QN vanish at x = 0 and we get,
|ψAlicei = |ψi |x+=0 , (S72) 8
Bob is at x+, so the state he sees is given by,
|ψBobi = |ψi |x+=T . (S73)
s Let us choose a normal distribution for C (ω) in terms of ω around ω0 with width of σ:
(ω−ω )2 s 1 − 0 s C (ω) = √ e 2σ2 C˜ , (S74) σ 2π where C˜ s is a real number. We assume that σ ω. This allows us to simplify (S71) to:
iHx+ σ ω0 ω0 ω0 − 2 |ψi = (1 + εO + εQ + εQ )e ~ |ψi | + + O ε , ε . (S75) U N x =0 ω We notice that (S75) can be utilized to write: ω ω ω − iHT 2 σ |ψ i = (1 + εO 0 + εQ 0 + εQ 0 )e ~ |ψ i + O ε , ε . (S76) Bob U N Alice ω
ω0 ω0 We highlight that 1 + εQU and 1 + εO are unitary operators:
Oω0 † + Oω0 = 0 , (S77)
ω0 † ω0 QU + QU = 0 , (S78)
ω0 while 1 + εQN is not unitary:
ω0 † ω0 QN + QN , 0 . (S79)
We refer to QU and O as the generators of quasi unitary operators, QN represents the generator of a non-unitary. All of QU , O and QN are generators of non-local operators. Each transverse polarisation of the U(1) gauge field gets corrected by the spacetime geometry as if it were a massless scalar. So (S76) is valid for the U(1) guage field too. This concludes the correction for the quantum realm.
Supplementary Note 5: Propagation in the Schwarzschild spacetime geometry
µ In a spacetime endowed with the metric gµν, a geodesic x (τ) can be obtained from an effective action: Z S = dτL , (S80)
µ ν L = gµν x˙ x˙ , (S81) where τ is an affine parameter. For the Schwarzschild black hole in the standard coordinates, this is,
2 m 2 r˙ 2 ˙2 2 2 L = − 1 − t˙ + m + r θ + sin θϕ˙ , (S82) r 1 − r where m = 2GN M• that M• is the mass of the black hole and c = 1 is understood. We choose the units such that m = 1. Due π ˙ to the spherical symmetry, without loss of generality, we can choose the equatorial plane θ = 2 and θ = 0 to describe any given geodesic at all times. The cyclic variables of ϕ and t lead to invariant quantities: ∂L = 0 → r2ϕ˙ = L , (S83) ∂ϕ ! ∂L 1 = 0 → 1 − t˙ = E . (S84) ∂t r We consider null geodesics reaching the asymptotic infinity and set E = 1. Due to the form of the Lagrangian, its Legendre transformation, which is the Lagrangian itself, is invariant. We consider a null geodesic, and set L = 0 giving, s ! 1 1 |r˙| = γ(r) = 1 − 1 − L2 . (S85) r2 r 9
π The non-zero components of the Riemann tensor on the geodesic (θ = 2 ) in the standard spherical coordinates are:
1 R = − , (S86a) trtr r3 Rθϕθϕ = r (S86b) r − 1 R = R = , (S86c) tθtθ tϕtϕ 2r2 1 R = R = − . (S86d) rθrθ rϕrϕ 2 (r − 1)
The coordinate independent representation of the Riemann tensor, therefore, follows:
1 r − 1 R = − (dt∧dr) ⊗ (dt∧dr) + r (dθ∧dϕ) ⊗ (dθ∧dϕ) + ((dt∧dϕ) ⊗ (dt∧dϕ) + (dt∧dθ) ⊗ (dt∧dθ)) r3 2r2 1 − ((dr∧dϕ) ⊗ (dr∧dϕ) + (dr∧dθ) ⊗ (dr∧dθ)) . (S87) 2(r − 1)
We would like to compute the components of the Riemann tensor in x+, xa¯ coordinates. To this aim, we first consider an infinitesimal displacement δ~x:
dr d~x = f dteˆ + eˆ + rdϕeˆ + rdθeˆ , (S88) t f r ϕ θ wheree ˆt,e ˆr,e ˆϕ ande ˆθ are the normalized unit vectors in t, r, θ, ϕ coordinates, and,
r 1 f = 1 − . (S89) r
The displacement vector in the x+, xa¯ coordinates follows
+ − 1 2 d~x = dx eˆ+ + dx eˆ− + dx eˆ1 + dx eˆ2 , (S90)
+ a¯ wheree ˆ+,e ˆ−,e ˆ1 ande ˆ2 are the normalized unit vectors in the x , x coordinates. We are not interested in explicit reproduction of the gravitational redshift. So instead of settinge ˆ+ equal to the the tangent of the geodesic, we sete ˆ+ proportional to the tangent of the geodesic where the vectors ofe ˆ−,e ˆ+,e ˆ1 ande ˆ2 are:
1 eˆ+ = √ +eˆt + r˙eˆr + r f ϕ˙eˆϕ , (S91a) 2 1 eˆ− = √ −eˆt + r˙eˆr + r f ϕ˙eˆϕ , (S91b) 2 eˆ1 = −r f ϕ˙eˆr + r˙eˆϕ , (S91c)
eˆ2 = eˆθ , (S91d)
We use (S91) in (S88) and equate it to (S90), and obtain:
1 dt = √ dx+ − dx− , (S92a) f 2 f r˙ dr = √ dx+ + dx− − r f 2ϕ˙dx1 , (S92b) 2 f ϕ˙ r˙ dϕ = √ dx+ + dx− + dx1 , (S92c) 2 r 1 dθ = dx2 . (S92d) r 10
+ We would like to compute R+a¯+b¯ . In doing so, we keep only the linear dx terms in the wedge products that appeared in (S87): r f ϕ˙ dt ∧ dr = rdx˙ + ∧ dx− − √ dx+ ∧ dx1 , (S93a) 2 r˙ dt ∧ dϕ = ϕ˙dx+ ∧ dx− + √ dx+ ∧ dx1 , (S93b) f r 2 1 dt ∧ dθ = √ dx+ ∧ dx2 , (S93c) r f 2 f dr ∧ dϕ = √ dx+ ∧ dx1 , (S93d) r 2 f r˙ dr ∧ dθ = √ dx+ ∧ dx2 , (S93e) r 2 f ϕ˙ dθ ∧ dϕ = − √ dx+ ∧ dx2 . (S93f) r 2
Utilizing (S93) in (S87) enables us to extract the components of R+a¯+b¯ :
2 r˙ r − 1 2 R − − = − + ϕ˙ , (S94a) + + r3 2r2 1 r˙2 f 2ϕ˙2 R = − + − , (S94b) +1+1 4r3 4r3 2r 1 r˙2 f 2ϕ˙2 R − , +2+2 = 3 3 + (S94c) 4r√ 4r 2r 3 2 R − = f r˙ϕ˙ , (S94d) + +1 2r2 R+−+2 = R+1+2 = 0 . (S94e)
2 As a consistency check, we observe that (S94) satisfies (S14), and R+−+2 and R+1+2 vanish due to parity symmetry in the x direction: x2 → −x2. We use (S83), (S85) and (S89) to simplify (S94):
3L2(r − 1) 1 R − − = − , (S95a) + + 2r6 r3 3L2(r − 1) R = − , (S95b) +1+1 4r6 3L2(r − 1) R = , (S95c) +2+2 4r6 r 3L 2 R − = 2 − r˙ , (S95d) + +1 2r4 r R+−+2 = R+1+2 = 0 . (S95e)
Eq. (S95) enables us to compute the G coefficients given (S11) for any geodesic.
A. Propagating along a general null geodesic
This section considers that Alice and Bob are at the same far distance ra from the black hole. Alice sends a wavepacket toward Bob. The wavepacket goes near the black hole and is deflected toward Bob, fig. S1. The package takes x+ = T of the affine + T + T parameter to reach Bob. We notice thatr ˙ is negative for x ∈ {0, 2 } and is positive for x ∈ { 2 , T}. To calculate the integration ˜ of the R+a¯+b¯ over the geodesic, we divide integration to intervals of positive and negativer ˙ :
Z T Z T Z r Z r Z r Z r 2 min dr a dr a dr a dr Ga¯b¯ = dτR+a¯+b¯ + dτR+a¯+b¯ = R+a¯+b¯ + R+a¯+b¯ + R+a¯+b¯ + R+a¯+b¯ , (S96) T r˙ r˙ γ(r) γ(r) 0 2 ra rmin rmin rmin where red color is used for the interval wherer ˙ is negative, and blue is used whenr ˙ is positive. In the last line, the integration limit of integrals are written from a smaller value of r to a larger one, andr ˙ = γ(r) is used for blue integral, andr ˙ = −γ(r) is 11
Alice
ra
WavepacketWave packet trajectory/ lighttrajectory beam Event horizon horizon of the black holeof the black hole
ra
rmin
Bob
Supplementary Figure S1: The geometry of communication for an arbitrary geodesic path, where two satellites share a wavepacket. The wavepacket traverses through the spacetime, goes near the black hole, and is deflected toward Bob. used for red integral where γ is defined in (S85). Utilizing (S95) then gives
Z ra dr G−− = 2 R+−+−, (S97a) rmin γ(r) G−1 = 0 , (S97b) Z ra dr G11 = −G22 = 2 R+1+1 . (S97c) rmin γ(r)
˜ where γ(r) is defined in (S85). We also need to calculate Ga¯b¯ in (S11):
Z T Z τ ˜ 0 0 Ga¯b¯ = dτ dτ R+a¯+b¯ (τ ) , (S98) 0 0
where T is the total amount of affine parameter that the light packet takes to reach from Alice to Bob. We would like to change dr the integration over τ to an integration over r by dτ = r˙ . We break the integral to intervals of constant sign forr ˙,
Z T Z τ Z T Z T Z T Z τ ˜ 2 0 2 0 0 Ga¯b¯ = dτ dτ R+a¯+b¯ + dτ dτ R+a¯+b¯ + dτ dτ R+a¯+b¯ , (S99) T T T 0 0 2 0 2 2 where red color is for the interval of negativer ˙ and blue colour is for positiver ˙. Next we change the integral from τ to r through dr dτ = r˙ :
Z rmin dr Z r dr0 Z ra dr Z rmin dr0 Z ra dr Z r dr0 G˜ − a¯b¯ = 0 R+a¯+b¯ 0 R+a¯+b¯ + 0 R+a¯+b¯ , (S100) ra γ(r) ra γ(r ) rmin γ(r) ra γ(r ) rmin γ(r) rmin γ(r )
Z ra dr Z ra dr0 Z ra dr Z ra dr0 Z ra dr Z r dr0 G˜ a¯b¯ = 0 R+a¯+b¯ + 0 R+a¯+b¯ + 0 R+a¯+b¯ . (S101) rmin γ(r) r γ(r ) rmin γ(r) rmin γ(r ) rmin γ(r) rmin γ(r ) Utilizing (S95) then gives:
Z ra Z ra dr 0 R+−+1 G˜−1 = −2 dr , (S102a) rmin γ(r) r r˙ ! Z ra dr G˜11 = −G˜22 = G11 . (S102b) rmin γ(r) ˜ To compute Ga¯b¯ , notice that it is given by
Z T Z τ1 Z τ2 ˜ ˜ Ga¯b¯ = dτ1 dτ2 dτ3Ra¯b¯ (τ3) . (S103) 0 0 0 We divide the domain of integration to intervals that have constant sign ofr ˙,
T T Z Z τ1 Z τ2 Z T Z Z τ2 ˜ 2 2 Ga¯b¯ = dτ1 dτ2 dτ3R+a¯+b¯ (τ3) + dτ1 dτ2 dτ3R+a¯+b¯ (τ3) T 0 0 0 2 0 0 T Z T Z τ1 Z 2 Z T Z τ1 Z τ2 + dτ1 dτ2 dτ3R+a¯+b¯ (τ3) + dτ1 dτ2 dτ3R+a¯+b¯ (τ3) , (S104) T T T T T 2 2 0 2 2 2 where the red color is associated to intervals wherer ˙ is negative, blue is associated to the intervals whereinr ˙ is positive. Now we change the integral from over τ to r:
Z rmin Z r1 Z r2 Z ra Z rmin Z r2 ˜ dr1 dr2 dr3 dr1 dr2 dr3 Ga¯b¯ = − R+a¯+b¯ + R+a¯+b¯ ra γ(r1) ra γ(r2) ra γ(r3) rmin γ(r1) ra γ(r2) ra γ(r3) Z ra Z r1 Z rmin Z ra Z r1 Z r2 dr1 dr2 dr3 dr1 dr2 dr3 − R+a¯+b¯ + R+a¯+b¯ , (S105) rmin γ(r1) rmin γ(r2) ra γ(r3) rmin γ(r1) rmin γ(r2) rmin γ(r3) wherer ˙ = γ(r) is used for the blue integrals,r ˙ = −γ(r) is used for red integrals, and γ is defined in (S85). We interchange the domain of all the integrals, such that all the integrals start from the smaller value of r:
Z ra Z ra Z ra Z ra Z ra Z ra ˜ dr1 dr2 dr3 dr1 dr2 dr3 Ga¯b¯ = R+a¯+b¯ + R+a¯+b¯ rmin γ(r1) r1 γ(r2) r2 γ(r3) rmin γ(r1) rmin γ(r2) r2 γ(r3) Z ra Z r1 Z ra Z ra Z r1 Z r2 dr1 dr2 dr3 dr1 dr2 dr3 + R+a¯+b¯ + R+a¯+b¯ . (S106) rmin γ(r1) rmin γ(r2) rmin γ(r3) rmin γ(r1) rmin γ(r2) rmin γ(r3)
Noticing that R+1+1 does not depend onr ˙, and after simplification we then obtain:
Z ra !2 Z ra Z ra ! ˜ ˜ dr dr1 dr2 G˜ 11 = −G˜ 22 = G11 − G11 rmimγ(r) rmin γ(r1) r1 γ(r2) Z ra Z ra Z ra dr1 dr2 dr3 + 2 R+1+1(r3) . (S107) rminγ(r1) r1 γ(r2) r2 γ(r3) Notice that in the simplification we haven’t tried to get the minimum number of integrations. We have ordered the boundary of all integrals from a lower value of r to ra. Eq. (S95), (S97), (S102) and (S107) identify all the G-coefficients of O (S20), QU (S24) and QN (S25) at Bob’s position for a general geodesic.
Supplementary Note 6: Deflection by the Sun
This section considers the communication between Alice and Bob in the weak regime of gravity of the Schwarzschild space- time geometry, see Fig. S2. Alice sends a wavepacket/beam toward Bob. The wavepacket reaches a distance of rmin =l from the 13
Bob
Alice ra l ra
Supplementary Figure S2: Alice and Bob are in the weak gravity regime at distance of ra from the central mass. Alice sends a pulse beam (contains information) toward Bob. The wavepacket reaches distance of l to the central mass. l is large and the wavepackage is very slightly deflected. Its trajectory remains almost a straight line. central mass and l is large. So the packet is very slightly deflected. Its trajectory remains almost a straight line. In so doing, we treat the Schwarzschild radius as a small parameter, and approximate γ defined in (S85) to r l2 γ = 1 − , (S108) r2 and simplify (S97), (S102) and (S107) to, √ a2 − 1 l2G = , (S109a) −− a3 l2G = 0 , (S109b) −1 √ a2 − 1 l2G = −l2G = (1 + 2a2) , (S109c) 11 22 2a3 p (a − 1) 2(a2 − 1) lG˜ = − , (S109d) −1 2a 1 lG˜ = −lG˜ = (a2 − 1)(1 + 2a2) , (S109e) 11 22 2a3 2 3 2 (a − 1) 2 (1 + a ) G˜ = −G˜ = − , (S109f) 11 22 2a3
where a = ra/l is a dimensionless parameter a. Utilizing (S109) in (S20), (S24) and (S25) yields: iω i Oω = − G (x2 − x2) + G˜ (x ∂1 − x ∂2) + G˜ (∂1)2 − (∂2)2 (S110a) 2 11 1 2 11 1 2 ω 11 ! i ω2 Qω = 1 − (x−)2 G + x−G˜ ∂1 , (S110b) U ω 2 −− −1 2i Qω = x−G + G˜ ∂1 . (S110c) N −− ω −1 Employing them in (S23) results in the corrected wave: Z µ iωx− ω ω ω (0) 2 ψ(x ) = dω e 1 + εO +εQU + εQN fω + O(ε ) . (S111)
We assume that Alice prepares a frequency Gaussian wavepacket with the central frequency ω0 and width σ:
2 1 (ω−ω0) (0) − 2 (0) + 1 2 fω = √ e 2σ f (x , x , x ) , (S112) σ 2π
(0) where f satisfies the wave equation and σ ω0. This assumption allows us to replace ω in the parentheses of the integral of (S111) with ω0, and perform the integration over ω to obtain:
− 2 ! − (σx ) ω ω εσ µ iω0 x − 2 ω0 0 0 (0) 2 ψ(x ) = e e 1 + ε O + QU + QN f + O ε , . (S113) ω0
1 2 (0) + a Alice sends a beam with a Hermite-Gaussian distribution in the transverse directions of x and x , note that fω (x , x ):= (0) + fω,m,n,q = HGq(x ) HGm,n(x1, x2). To represent such a solution, let w0 be the initial width of the beam, zR be the Rayleigh range, and w(x+) be the width of the packet at x+ s + !2 + x w(x ) = w0 1 + , (S114) zR 1 z = ω w2 . (S115) R 2 0 0 Define δ(x+) as the Gouy phase: ! x+ δ(x+) = (m + n + 1) arctan . (S116) zR
For Hermite-Gaussian modes of m, n around the mean value of ω0 with the width of σ: √ √ ! 2 2 (σx−)2 µ iω x− − w0 2x1 2x2 ω0(x1 + x2) + ψ x C e 0 e 2 − iδ x , Alice( ) = m,n + Hm + Hn + exp + + ( ) (S117a) w(x ) w(x ) w(x ) 2zR − 2ix where Hn(.) are Hermite polynomials, and m, n are positive integer values that define the transverse modes. Figure 4-(a) shows the intensity profile of Hermite-Gauss modes for first 9 modes, m, n ∈ {0, 1, 2}. Note that (S117) represents exact solutions for beams propagating in Minkowski spacetime geometry. We express (S110) at x+ = T in terms of x and y:
2iG˜ O = −iω w2(T)G (x2 − y2) + G˜ (x∂ − y∂ ) + 11 ∂2 − ∂2 , (S118a) 0 11 11 x y ω w2(T) x y √ 0 i ω2 2G˜ Q = (1 − 0 (x−)2)G + −1 x−∂ , (S118b) U ω 2 −− w(T) x 0 √ − 2 2i QN = x G−− + G˜−1∂x , (S118c) ω0w(T) √ √ + + where x = 2x1/w(x ) and y = 2x2/w(x ) are dimensionless coordinates, and √ T = 2l a2 − 1, (S119) is the time at which the wavepacket reaches to Bob – for simplicity, we have dropped the superscript of ω0. Thus, the wavepacket at T is, µ 2 ψBob(x ) = ((1 + ε (O + QU + QN )) ψAlice) |x+=T + O ε , εσ . (S120) We would like to study the correction for the solar system, when Alice and Bob are at the mean-Earth-Sun distance from the sun and the wavepacket passes at l = 2R , Fig. S3. We assume that zR T. So we set:
m = m = 2953 meters , (S121) 11 ra = 1.51 × 10 meters , (S122) 9 l = 2R = 1.39 × 10 meters , (S123) a = 218.6 . (S124) This allows us to approximate: T w(T) = w0 , (S125) zR iT ∂x ≡ − x , (S126) zR iT ∂y ≡ − y , (S127) zR 15
Alice Bob
l ra ra
Supplementary Figure S3: Alice uses the large telescopes on Earth, and sends a Hermite Gaussian signal toward Bob such that the width of the signal remains smaller than the Sun’s radius at its closest encounter to the Sun. where
T = 2la, (S128)
For the case that we are interested in, a is large. So we approximate (S109) for large a, where the dominant term is given by O:
19ia2 O = − (x2 − y2) . (S129) 2zR
zR We observe that, for w0 a , that holds true for a large value of zR, the width of the package near the Sun remains much smaller than the Sun’s radius. We assume that σ is not very sharp, so we ignore terms proportional to x−. We assume that a is large. So O gives the largest contribution. We, therefore, keep only O and write:
2 ! 9.5m a 2 2 ψBob = 1 − i (x − y ) ψAlice|x+=T (S130) zR 2 ! 28 km × a 2 2 = 1 − i (x − y ) ψAlice|x+=T (S131) zR 9 ! 1.12 × 10 meters 2 2 = 1 − i (x − y ) ψAlice|x+=T . (S132) zR
2 For zR < 28 km × a , the correction becomes larger than 1 and we need to take into account higher ε terms. In order to keep the correction perturbative, we should choose either a sufficiently large value of zR or a null geodesic with smaller value of a. For 2 9 zR = 2.8 km × a = 1.34 × 10 m, the correction remains small and the wavepacket at Bob’s position is,
µ 2 2 µ ψBob(x ) = (1 − 0.10i(x − y ))ψAlice(x )|x+=T . (S133)
The correction term εψ(1), thus is,
(1) µ 2 2 µ εψ (x ) = 0.10i(x − y )ψAlice(x )|x+=T . (S134)
We have depicted the amplitude and phase of the εψ(1) in Fig. 4. We highlight that (S130) can be extended to the communication between Alice and Bob around the Earth as well,
2 ! µ 9.5m⊕a 2 2 µ ψBob(x ) = 1 − i (x − y ) ψAlice(x )|x+=T (S135) zR 2 ! 8.4 cm × a 2 2 µ = 1 − i (x − y ) ψAlice(x )|x+=T . (S136) zR
We, however, see that for sufficiently large value of zR that satisfies all the constraints, the magnitude of the correction is far smaller than what we have predicted for the correction around the Sun.
104
1000
a 100 r
10
1 5 10 50 100
rmin
Supplementary Figure S4: The shaded area shows the regime where ra and rmin for which |G˜11| remains smaller than 1. The unit of length is chosen equal to the Schwarzschild radius.
Supplementary Note 7: On Gravitational Decoherence
The consistency of the perturbation requires that all the ε terms in (S20), (S25) and (S25) remain smaller than 1. For a wavepacket with a frequency width of σ and beam width of W, we can estimate the magnitude of each term by replacing x with 1 − 1 ˜ W, ∂ with W and x with σ . In particular, estimating the contribution of G11 in (S102b) to be smaller than 1, yields,
2 T , (S137) |G11| where,
Z ra dr T = 2 , (S138) rmin γ(r) is the time that the packet has felt the curved spacetime geometry during its journey. G11 is a non-zero small parameter presented in (S97c). So (S137) identifies the maximum amount of time that the wavepacket can retain its phase and coherence intact. For |G | m rmin ra and large l, (S109c) will prove that 11 = l2 and simplifies (S137) to,
l2 cT , (S139) Max m which shows that the curvature of the spacetime geometry ultimately changes the phase of the wavepacket moving on any null geodesic. For T > TMax, Bob would need to take into account all the ε terms to recover the phase shift of the wavepacket that Alice has sent. This means that Bob must calculate the corrections up to an infinite derivatives of the Riemann tensor. Taking into account all the derivatives of the Riemann tensor in a smooth geometry, is tantamount of knowing the exact value for the Riemann tensor in whole of the spacetime geometry, a piece of knowledge that is not practically possible to gain. So it is tentative to argue that for T > TMax, Bob practically does not have any chance to compute the effect of the correction due to the curvature of the spacetime geometry and recover what Alice has sent. We refer to this phenomenon as practical gravitational decoherence. For small values of l, the null geodesic goes very close to the event horizon of the black hole. We have numerically performed the integration and calculated G˜11 for small values of rmin. Fig. S4 shows the regime of ra and rmin where G˜11 remains smaller than 1. We observe that when the wavepacket passes very close to the event horizon of the black hole, G˜11 becomes larger than one at around the black hole. This implies that a close encounter with the event horizon of the black hole completely changes phase of the wavepacket. The curvature of the spacetime geometry near the event horizon modifies the phase of the wavepacket to such a degree that the perturbation cannot be used to predict what far observers see. This could be interpreted as gravitational decoherence in the perturbative regime, and implies that Bob and Alice cannot easily communicate over a null geodesic that passes very close to the event horizon of the black hole, where the black hole adds lots of noise, and reducing the noise would be a computationally extensive or impossible task. 17
Supplementary Note 8: On detecting the phase distortion near the Earth
We assume that Alice sends a Gaussian wavepacket around ω = ω0 with the width of σ towards Bob. In other words, Alice at ra produces a wavepacket in the form of
− − (σx )2 (0) + 1 2 iω0 x − ΨAlice = f (x , x , x )e e 2 , (S140) and sends it towards Bob at rb. The wavepacket moves along a radial null geodesic and reaches Bob. A radial null geodesic that represents an outgoing beam. It is described by l = 0 in (S95). The only none-zero components of R+a¯+b¯ for an outgoing beam, therefore, is, 1 R = − . (S141) +−+− r3 + Equation (S85) implies that radial geodesic holdsr ˙ = 0 which is solved to x = r − ra where Alice is located at ra. The only non-zero components of (S11) contributing to the correction is, ! Z dτ Z r dr 1 1 1 G − − − − . −− = 3 = 3 = 2 2 (S142) r ra rr˙ 2 ra r We can utilize (S142) in (S20), (S24) and (S25), and obtain,
Oω = 0 , (S143a) 2 ! ! i (ωx−) 1 1 Qω − − − , U = 1 2 2 (S143b) 2ω 2 ra r ! 1 1 1 Qω − − x− . N = 2 2 (S143c) 2 ra r The transverse structure of the wavepacket is not affected, as there is no derivative with respect to spatial transverse directions present in (S143). This is due to high amount of symmetries present, the background is static and spherical, and the null geodesic respects the symmetries. In order to use the perturbative solution of the previous section, we notice that the curvature of the spacetime geometry during its journey has its largest value at Alice’s position, where the Kretschmann invariant [36] is, 12 K = R Rµνλη = . (S144) µνλη 6 ra We use the square root of the Kretschmann invariant as a measure to estimate the curvature of the spacetime geometry, i. e., 4 K ω0. We, therefore, impose, √ 2 3 ω2r3 , (S145a) √ 0 a 2 3 2 3 σ ra , (S145b) then σ ω0 reduces (S145) to √ 2 3 2 3 σ ra , (S146a) σ ω0 . (S146b)
ω ω Since σ ω, QU and QN at Bob’s position can be approximated to:
− 2 ! i (ω0 x ) 1 1 Qω = − 1 − − , (S147a) U ω 2 2 2 0 2 ra rb x− 1 1 Qω = − − , (S147b) N 2 2 2 ra rb here, ω is approximated to ω0 in (S143). Therefore, Bob receives the following wavepacket,
− 2 ΨBob ω ω iω x− − (σx ) 2 = 1 + εQ 0 + εQ 0 e 0 e 2 + O(ε ) , (S148) f (0) U N 18
ω ω where QN and QU alter the amplitude and phase of the wavepacket, respectively. The change in the amplitude of the wavepacket is given by,
− 2 − − 2 ω − (σx ) 2 εx 1 1 − (σx ) 2 δA = εQ 0 e 2 + O(ε ) = − − e 2 + O(ε ) . (S149) N 2 2 2 ra rb
− 1 The maximum change in the wavepacket amplitude occurs at x = ± σ , and the maximum change yields 1 1 1 1 |δAmax.| = − √ . (S150) σ 2 2 2 ra rb e The consistency of the perturbation requires δA < 1, and thus, results in,
2 2 1 ra 1 − < σ2r3 . (S151) er 2 a 4 a rb
We notice that (S151) is already satisfied by (S145) since 1 ra. The imaginary part of the correction adds to the phase of the wavepackage, which is,
− 2 − 2 ω − (σx ) 2 iω0 1 1 − 2 − (σx ) 2 δχ = εQ 0 e 2 + O(ε ) ≈ − (σx ) e 2 + O(ε ) . (S152) U 2 2 2 4σ ra rb √ − 2 The maximum change in the wavepacket phase occurs at x = ± σ , which is, ω0m⊕ 1 1 δχ = − . (S153) max. 2 2 2 2eσ ra rb
2GM⊕ where m⊕ is the Schwarzschild radius of Earth (m⊕ = c2 ) – here, c is recovered. Since the spacetime geometry is static, and considering that the equation we have solved are linear, there is no need to assume that the length of the wavepacket is small in the x− coordinate. We can extend the correction for small values of σ. This can be easily understood as considering segment of a large wavepacket in x− direction, solving the equation for each segment and adding the solutions for all segments. Let us assume that Alice is on the surface of the Earth with ra = 6400 km while Bob is in the International Space Station with rb = 6800 km. Equation (S153) then gives: ν δχ = (2.2 × 10−10Hz) 0 , (S154) max. ∆ν2 where the laser line-width ∆ν = 2πσ, and ν0 = 2πω0 are used. Note that the size of the wavepacket in time should not be larger than the distance between Alice and Bob, i.e., 400 km, or equivalently
∆ν > 749 Hz . (S155)
In deriving (4) from (3), it is assumed that the metric inside the wavepacket can be approximated to a constant metric. Equation (S95) gives the ε corrections to the metric, so they can be approximated to constant parameter inside the wavepacket. For radial null geodesic L = 0, only R+−+− is non-vanishing. For communication over radial null geodesic between rb and ra (ra < rb), the difference in the component of the metric gµν inside the wavepacket to be much than one gives
∆gµν ∆R+−+− = 3 . (S156) gµν R+−+− Utilizing the (S95), and noting that the linewidth of ∆ν corresponds to a wavepacket with size of c/∆ν in time direction, then implies 3c ∆ν. (S157) ra
Therefore, for ra = 6400 km, we get the lower bound of
∆ν 149 Hz (S158) 19
There exists commercial portable continuous wave (cw) lasers with the ultra narrow linewidth of 1 Hz for λ = 657 nm (https://www.menlosystems.com/products/ultrastable-lasers/ors-cubic/). However, in our calculation, we assume Alice gener- ates a very good coherent Gaussian beam in the form of (S140) with linewidth of 1 kHz. Notice that ∆ν = 1 kHz satisfies both (S155) and (S158). For example, Eq. (S154) for ν0 = 456 THz (λ0 = 657 nm), and ∆ν = 1 kHz gives a phase alteration of
δχmax. = 0.10 , (S159) which is of a measurable magnitude. Although measuring phase of a wavepacket with a high-precision is not a big challenge, in order to close this section we suggest the following strategy. We consider that Alice and Bob have the same equipment to produce wavepackets. Alice first sends a wavepacket in the form of (S140) toward Bob. Bob receives the wavepacket, measures its amplitude over time and identifies σ. Bob generates the wavepacket that Alice has sent, identifying as a local oscillator. We call this as the “standard” wavepacket. Alice prepares the second wavepacket in the form (S140) and sends it to Bob. Alice sends the second wavepacket τA seconds after sending the first one. Bob in his frame would receive the second wavepacket τB seconds after getting the first package. τA and τB are not the same due to gravitational and relativistic redshift. Bob, however, knows when Alice will send the second wavepacket and can compute when he would receive the wavepacket. Bob generates a standard wavepacket and measure its phase difference with the second wavepacket. The phase difference is a function of time, i.e., δχexp = δχexp(t). Equation (S152) maps the quadratic term in t to the distortion due to the curvature of the space-time geometry. So Bob fits χexp to the best quadratic polynomial in t:
2 δχexp = (a + δa)(σt) + (b + δb)(σt) + (c + δc) , (S160) and finds the numbers associated to a, b, c, and their errors δa, δb and δc. Notice that b encodes un-accounted Doppler or gravitational redshift while a is due to (S152), and is given by: ω0m⊕ 1 1 a = − , (S161) 2 2 2 4σ ra rb The magnitude of (S159) shows that a can be measured. However, the atmosphere between ISS and the ground base can introduce an error of such a magnitude that would not allow us to measure a between Earth and the ISS. The experiment, therefore, should be performed between two satellites which are in different latitudes. We can substitute ra and rb in (S161) with the latitudes of the two satellites, and (S161) describes the correction when they communicate over a radial null geodesic. For example for radial communication between satellites with ra = 6600 km and rb = 7000 km, (S153) is given by: ν δχ = (1.98 × 10−10Hz) 0 , (S162) max. ∆ν2 that ν0 = 456 THz, and ∆ν = 1 kHz results
δχmax. = 0.09 , (S163) which is of a measurable magnitude and is free of atmosphere turbulence as there exists no air between the two satellites.