Quantum communications and quantum metrology in the spacetime of a rotating planet Jan Kohlrus,1 David Edward Bruschi,2, 3 Jorma Louko,1 and Ivette Fuentes1, 4 1School of Mathematical Sciences, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom 2Racah Institute of Physics and Quantum Information Science Centre, the Hebrew University of Jerusalem, Givat Ram, 91904 Jerusalem, Israel 3York Centre for Quantum Technologies, Department of Physics, University of York, Heslington, YO10 5DD York, United Kingdom 4Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Wien, Austria We study how quantum systems that propagate in the spacetime of a rotating planet are affected by the curved background. Spacetime curvature affects wavepackets of photons propagating from Earth to a satellite, and the changes in the wavepacket encode the parameters of the spacetime. This allows us to evaluate quantitatively how quantum communications are affected by the curved spacetime background of the Earth and to achieve precise measurements of Earth's Schwarzschild radius and equatorial angular velocity. We then provide a comparison with the state of the art in parameter estimation obtained through classical means. Satellite to satellite communications and future directions are also discussed. I. INTRODUCTION of the parameter in a single-shot run [4]. Typical applications of quantum metrology range from Quantum communications is a rapidly growing field phase estimation in quantum optics to estimating the which promises several technical improvements to cur- gravitational potential with Bose Einstein Condensates rent classical communications. One example is the use of (BECs) [5, 6]. However, when estimating relativistic pa- quantum cryptography which would make communica- rameters, gravity usually appears as an external poten- tions more secure, thanks to more robust protocols than tial, or a phase modification, which does not overcome the classical ones [1]. More fundamental aspects can also the inherent inconsistency between quantum physics and be studied using quantum communications. For exam- relativity [7]. Recently, this gap has been bridged and ple, the interplay between quantum physics and relativ- quantum field theory in curved spacetime has been em- ity can be probed through quantum communications be- ployed as the core framework to compute the ultimate tween moving observers and within schemes in a curved bounds on ultra-precise measurements of relativistic pa- spacetime background [2]. The results of the measure- rameters. In particular, it was shown that it is possible ments can then be compared with the predictions ob- to use the shifting induced on the frequency distribution tained by theories that were developed in the overlap of of single photons ascending the gravitational potential of quantum physics and relativity, the most well known and a static planet to estimate with great precision the dis- understood of these being Quantum Field Theory (QFT) tance between a user based on Earth and one on a satel- in curved spacetime [3]. lite [2, 8]. In this case, gravity isn't affecting the quantum Knowing quantitatively how quantum communications state as the simple addition of a phase. The effects due are affected by the curved spacetime background would to curved spacetime can therefore not be explained by a enable to compensate undesirable relativistic effects in simple ad hoc implementation of proper time in a classical future quantum technologies. Precise values of the nec- quantum mechanics scheme. Furthermore, it was shown essary corrections in such quantum communication se- that these effects can have potentially high impact on tups can only be obtained with an accurate knowledge of specific types of quantum key distribution (QKD) proto- arXiv:1511.04256v6 [quant-ph] 8 Aug 2017 the spacetime parameters. We thus need to employ tech- cols [2]. This direction has the potential of leading to- niques from the field of quantum metrology, which aims wards the development of new relativistic and quantum at exploiting quantum resources, such as entanglement, technologies aimed at testing the predictions of quantum to estimate physical parameters [4]. Within a standard field theory in curved spacetime in space-based experi- estimation protocol, an input quantum state undergoes ments with satellites. a transformation that encodes the parameter to be esti- In this work we extend the analysis carried out in pre- mated. The resulting state of this transformation is then vious works which investigated quantum estimation tech- compared, by means of the fidelity, to a neighbouring niques in scenarios where photons are exchanged between state which is infinitesimally close in terms of the param- Earth and a satellite [8]. There, the Earth was assumed eter. One can define a distance between these two states as static and the effects on the propagation of the photons that is directly related to the Quantum Fisher Informa- depend only on the Schwarzschild radius of the Earth. tion (QFI), which in turn is directly related to the max- Here we consider a rotating planet, and we model the imum precision one can obtain in an estimation scheme. metric outside the mass distribution by the well known A final measurement provides an estimation of the value Kerr metric [9]. The transformation induced by the cur- 2 vature on the traveling photon reduces to a beam-splitter, coordinates (t; r; φ) reads [9]: a well known linear transformation in quantum optics 2M 1 [10]. We can therefore restrict ourselves to Gaussian ds2 = − 1 − dt2 + dr2 states and employ the powerful covariance matrix formal- r ∆ ism that allows to achieve analytical insight in scenarios 2Ma2 4Ma + r2 + a2 + dφ2 − dt dφ, (1) that involve Gaussian states and linear unitary transfor- r r mations [11, 12]. In particular, we seek out the effects of 2M a2 ∆ = 1 − + : (2) rotation on previously employed entanglement-swapping r r2 protocols [2, 8]. We find the error bound on the equatorial angular ve- For clarity, we will consider the rotating planet to be the locity of the Earth and compare it with that achieved Earth, with mass M, radius rA, angular momentum J and Kerr parameter (i.e., angular momentum per unit with cutting edge technology. The rate of improve- J ment of quantum optical technologies and the rapid in- mass) a = M . crease of the control over quantum systems suggest that A photon is sent radially by Alice from a laboratory in the near future our scheme might provide a reliable on Earth's equator to Bob who is in a satellite circularly way to outperform current technologies based on classi- orbiting at radius rB in the equatorial plane of the Kerr cal means. spacetime. A schematic representation of the setup can be found in Fig. 1. The paper is organised as follows. In Section II, we present the process of exchanging photons between Earth and a satellite, we characterise and model the system, and we give the mathematical formalism that is going to be relevant for the general relativistic calculations that will follow. In Section III, we derive the expression of the frequency shift for the photon travelling through the Kerr spacetime. Section IV consists of the relativistic quan- tum metrology calculations. It introduces the relevant perturbative quantities that are affecting the states, and derives the Quantum Fisher Information (QFI) for the system studied, and hence the estimated error bounds for the spacetime parameters. Section V introduces the satellite to satellite scheme and the related precision es- timations are computed in the same fashion as in the Earth to satellite case. Finally, Section VI briefly dis- cusses how the effects computed in this work can affect FIG. 1: Alice on Earth sends a photon (localised around a simple QKD protocol, specifically comparing the mag- the straight line) to Bob in the satellite. The photon ex- periences the effects of the curvature of spacetime along the nitude of the effect with what has been found in [2]. whole path while propagating, which can be seen in the pic- Throughout the whole paper we employ geometrical ture by the progressive tightening and redshifting of the Gaus- units G = 1 = c. Relevant constants are restored when sian wavepacket. The final effect is a nonlocal and cumulative needed for the sake of clarity. Vectors and matrices effect due to travel along the whole path. are denoted in bold characters. Vectors are written us- ing the usual differential geometry notation [13], namely t r θ φ t r θ φ X = (X ;X ;X ;X ) = X @t +X @r +X @θ +X @φ. Einstein's summation convention is assumed on repeated B. Wave packet characterisation Greek indices. A and B indices denote evaluations at Alice's and Bob's events respectively. In this work we employ photons which are geomet- rically radial, namely with vanishing angular velocities _ _ φγ = 0 = θγ . We will see that in Kerr space-time such II. INTRODUCTION TO THE FORMALISM photons have a non trivial angular momentum constant of motion, contrarily to the Schwarzschild case. The evo- lution of the quantum field is thus a 1+1 dimensional A. Description of the experiment problem. Such a photon can be modelled by a wave packet of frequency distribution F!0 (!) of monochro- In this work we consider a spherical planet that rotates matic plane waves with frequency ! and peaked at !0. slowly. The Kerr metric can be used, to good approxi- The annihilation operator associated to this photon by mation, to model the spacetime background around the an observer (infinitely) far from Earth is: rotating planet [9]. Our work will be constrained to the π Z +1 equatorial plane θ = 2 to be able to work with simple an- −i!t a!0 (t) = d! e F!0 (!) a!: (3) alytical formulas. The reduced metric in Boyer-Lindquist 0 3 The canonical bosonic commutation relations with the amount of shifting being quantified by the new [a (t); ay (t)] = 1 for the bosonic operator (3) at parameter δ defined by1: !0 !0 any instant of time follow directly from those for the y 0 δ = f − 1: (9) sharp frequency operators [a!; a!0 ] = δ(! − ! ) and from the normalisation of the frequency distribution function R 1 2 In the following, we will derive the explicit formula for f F!0 (!) i.e., 0 d!jF!0 (!)j = 1.