A blueprint for a simultaneous test of quantum mechanics and general relativity in a space-based quantum optics experiment

In this paper we propose an experiment designed to observe a general-relativistic effect on single photon interference. The experiment consists of a folded Mach-Zehnder interferometer, with the arms distributed between a single Earth orbiter and a ground station. By compensating for other degrees of freedom and the motion of the orbiter, this setup aims to detect the influence of general relativistic time dilation on a spatially superposed single photon. The proposal details a payload to measure the required effect, along with an extensive feasibility analysis given current technological capabilities.


Introduction
Since the beginning of the 20th century, quantum mechanics and general relativity have been the theoretical foundations of modern physics. However, while their predictive power has enabled us to understand the universe at both very small and very large scales, their validity has only been explored separately. Regarding general relativity, historic tests include measuring the perihelion of Mercury [1] and the Pound-Rebka experiment to quantify gravitational redshift [2]. Even with the precision attainable with contemporary experiments, general relativity remains robust; for example, estimates of the Shapiro delay (the slowing of light as it passes by a massive body) have been acutely measured via radio links with the Cassini spacecraft [3]. Additionally, tests of the (weak) equivalence principle now concur with general relativity with accuracy 10 −13 [4], and will be bounded by 10 −15 with the space-based MICROSCOPE experiment [5]. This resilience to experimental analysis is also shared by quantum mechanics. The most stringent tests of quantum theory are precision tests of quantum electrodynamics (QED); these have now been verified to an accuracy of 10 −12 using a one-electron quantum cyclotron [6].
While quantum field theory, which is currently the best description we have of the quantum world, is fundamentally incompatible with general relativity, descriptions of quantum field theories on static but curved background spacetimes are internally consistent and produce novel predictions [7] (such as Hawking radiation [8]). However, any full description of gravity itself as a quantum field theory is bound to fail due to a lack of renormalisability. Additionally, quantum field theory and general arXiv:1611.01327v1 [astro-ph.IM] 4 Nov 2016 relativity have contrasting descriptions of physical parameters, the most notable of which is the nature of time; this "problem of time" is an unresolved hurdle in producing a canonically quantised theory of gravity [9].
Without exception, all of the experiments mentioned above are tests of either general relativity or quantum mechanics, but not the interplay of both. As such, their outcomes are consistent with classical dynamics evolving on a curved background (in the former case), or of a quantum field theory evolving in flat space (in the latter case). A controlled experiment exploring quantum dynamics on a curved background is, to the authors' knowledge, yet to be performed. To this end, a theoretical proposal by Zych et al. [10] outlined the premise of a controllable experiment on a system that is both highly quantum, and occurring on a non-trivial background spacetime. Such an experiment would provide vital evidence to aid the successful unification of quantum mechanics and general relativity.
The crux of the proposal in [10] is to perform interferometry of single quanta, but to orient the interferometer in such a way that quantum states in each arm traverse a different gravitational equipotential. Historically, such an idea is not new; the analysis of gravitational effects on matter interferometry was first explored by Collela, Overhauser and Werner ("COW") in 1975 [11]. While the COW experiment did observe a gravitational influence on the output of the interferometer, their data is adequately described by analysing a Mach-Zehnder interferometer placed in a Newtonian gravitational field, rather than a truly relativistic spacetime. The key physical extension of the COW experiment present in [10] is to construct an interferometer large enough that general relativistic effects become apparent in the output of the interferometer, and to let the position of the photon in the interferometer serve as the local clock. Similar proposals exploring gravitational effects on quantum optical systems consider interferometry with a precessing polarisation serving as the clock [12], or explore the degradation of entanglement due to gravity [13].
The objective of this paper is to present a full design of a space-based quantum optics experiment capable of carrying out the experiment proposed by Zych et al. [10], along with a discussion of the practical feasibility of such an experiment. The design consists of a "folded" Mach-Zehnder interferometer distributed between an orbiting satellite (containing a single photon source, a spool of single-mode optical fibre, and a transmitting telescope) and a ground-station (equipped with a receiving telescope, an identical spool of fibre and a bank of single photon detectors). A stripped-down schematic of the proposed design is shown in Figure 1. Moreover, we contend that overcoming the engineering requirements for this proposal would be of significant benefit to those interested in developing quantum communications networks on a global scale.
The paper breakdown is as follows: in Section 2, we extend the analysis present in [10] to derive the predicted interferometric signal given that the satellite is in motion. In Section 3, we give an extensive feasibility analysis of the experiment. This includes a proposed setup and operational parameters for key optical components in Section 3.1, followed by analysis of interferometric stability, systematic error, loss and noise in Sections 3.2, 3.3, 3.4 and 3.5, respectively. All of these elements are combined with a statistical analysis in Section 4, in order to derive the operational duration necessary to test the hypotheses in [10] to a prespecified confidence. Finally, we outline orbit and mission requirements in Section 5 and highlight operational risks to be mitigated in Section 6.

Scientific background
While the experiment proposed by Zych et al. in [10] calculates the general relativistic effect on the fringe visibility of the interferometer output signal, we must also introduce special relativistic corrections due to the orbital velocity of the satellite (Section 2.1). This full expression for the time dilation is then used to derive an expected output for the Mach-Zehnder interferometer in a regime where relativistic corrections are important (Section 2.2).

Relativistic time dilation
General relativity predicts that a clock in a lower gravitational potential ticks slower than a clock in a higher one. This means that two clocks initially synchronised and then moved to points of differing gravitational potential accumulate an increasing time delay. This gravitational time dilation is observed daily in global navigation satellite systems, where one also has to account for special relativistic time dilation due to the motion of the positioning satellites. The total time delay can be derived directly in a general relativistic framework starting from the spacetime metric, using the weak field approximation for the Earth's gravitational potential and assuming the speed of a satellite to be much smaller than the speed of light. Figure 2 The total time delay between the fibre clocks, as a function of height, for a fixed fibre length and fibre refractive index. Magenta lines mark lowest and highest altitudes at which measurements could be performed (see feasibility analysis). In this region the time delay is dominated by the Shapiro delay. The choice of fibre length and refractive index is discussed in Section 3. The altitudes and orbital parameters are discussed in Section 5.
It is important to note that the time dilation is independent of the physical realisation of the clock. For example, we could consider an optical fibre loop and use the position of a photon travelling in the fibre as a clock. If two of these fibre-clocks were placed at different heights in the Earth's gravitational field, the one at the lower position would tick slower, meaning that the photon would take more time to travel the entire length of the fibre. Such a time delay is known as the Shapiro delay.
Considering the simplified experimental setup described in Figure 1, we deduce that the time spent by a photon in the satellite loop differs from the time spent in the ground loop, as measured by an observer on the ground. The accumulated time delay, i.e. the difference in photon arrival times, can be obtained by generalising the calculation described in [10] (see Appendix A). For a satellite on an elliptical orbit around the Earth, the time delay is given by where l is the proper length of the satellite fibre, l + dl is the proper length of the ground fibre, n is the fibre's refractive index, G is the gravitational constant, M is the mass of the Earth (5.97219 ×10 24 kg), R ⊕ is the radius of the Earth (6378100 m), the geoidal potential W 0 = L G c 2 where L G = 6.969290134 × 10 −10 by definition, c is the speed of light in vacuum, h is the altitude of the satellite with respect to the Earth's surface and a is the semi-major axis of the orbit. In order Figure 3 The expected interferometer signal vs. physical path length difference of the two interferometer arms, with the parameters as described in the text, at three different heights of the satellite. The signal was calculated by numerically evaluating Eq. (2) with a frequency-dependent refractive index as described in [15].
to compute the gravitational potential experienced by the satellite, we assumed the Earth is a perfect sphere. We can safely neglect the corrections due to the irregular shape of the Earth, both on the gravitational potential and on the satellite's orbit, since they are orders of magnitude smaller than the first order approximation [14]. A plot of this expression, in the case dl = 0, is shown in Figure 2. Note that the general relativistic and special relativistic time dilations act in opposition, such that at low satellite altitudes the time delay is positive and dominated by the special relativistic dilation, whereas at high altitudes it is negative and dominated by the general relativistic dilation. Therefore care must be taken to launch the satellite into a sufficiently high orbit for the general relativistic time dilation to be the dominant effect. Figure 2 demonstrates that the Shapiro delay for a satellite at the orbital heights considered in this proposal is of the order 100 fs.

Interferometry in the presence of gravity
Given this modification to the expected time dilation, we provide a brief exposition on the theory presented in [10], in order to motivate the design of our experiment. The reader is referred to [10], Sections 2 and 3 for a more thorough treatment.
The output of a Mach-Zehnder interferometer is sensitive only to a difference in the relative optical path length of each arm [16]. If the interferometer were to have two ideally-controlled arms of equal proper length (i.e. dl = 0), the relative optical path would be dependent only on the Shapiro delay ∆τ . However, by introducing a controlled difference in the relative length of the two interferometer arms by actively controlling the length of the fibre spools, one could cancel out ∆τ and recover maximal contrast at the output of the interferometer for any height of the satellite. Measuring the path-length difference that gives maximum contrast at different gravitational potentials would constitute a measurement of the Shapiro delay between the parts of the photon state traversing the upper and lower arms of the interferometer.
It can be shown (see Appendix B) that the probability of the photon being detected at the ± outputs of the Mach-Zehnder interferometer (with a single photon input) are given by: where ν is angular frequency in rads −1 defined in the reference frame of the observer at distance r from the Earth, ∆τ is the time delay as given in Eq. (1), P ± is the probability of the photon being detected at the ± output, and f (ν) is the spectrum of the light source, normalised such that ∞ −∞ |f (ν)| 2 dν = 1. In Figure 3 we numerically evaluated Eq. 2 with a range of different fibre lengths for three different satellite heights. We assumed a Gaussian spectral density for the , with parameters taken from the setup in Section 3: a central frequency ν 0 = 2π ×c/1550 nm and width σ = (100 fs/2π) 2 . We take the satellite fibre to be 60 km long, while the length of the ground station fibre differs by the amount given on the x-axis of the figure. The frequency-dependent speed of light in the fibres was calculated from the refractive index for fused silica given in [15].
In our experiment the satellite is sent on an elliptic orbit and we record the physical path-length difference required to recover full visibility, for a range of satellite heights. From this, the corresponding Shapiro delay as a function of satellite height can be inferred.

Implementation
In order to present a feasible implementation, it must be demonstrated that enough data can be collected within the lifetime of the satellite to confirm, or refute, the hypothesised signal to a pre-specified confidence. In this particular proposal, there are four key issues that must be addressed before feasibility can be demonstrated: stability of the interferometer, loss from the source to the detectors, noise from extraneous and atmospheric photons, and mitigation of any other degrees of freedom besides the Shapiro delay that might introduce distinguishability (in particular, photon dispersion). We first outline a feasible setup in Section 3.1, before presenting an analysis of these key issues.

Optical setup and components
A detailed illustration of the interferometer setup, with all the key optical components, can be seen in Figure 4.
The choice of operational parameters is ultimately a compromise between what is feasible to implement and what is necessary to recover the desired data. As such, there is no single optimal set of components; relaxing the constraints on their performance may be acceptable if they become more durable and the satellite can operate for longer and take more data, for example (  dispersion, noise and random fluctuations in signal by integrating counts over long time periods; hence, durability of components is critical). Therefore, the following choices of components and parameters are not intended to be a fixed specification for the payload, but to demonstrate one setup that could feasibly produce meaningful results.
We highlight the following optical elements, whose operational parameters are crucial to the success of the experiment.
Single photon source: An off-the-shelf, mode-locked, pulsed laser with operational wavelength of 1550 nm and 1 GHz repetition rate is attenuated to a mean photon number of 0.1 photons per pulse. While an attenuated laser is a straightforward option, true, heralded single photon sources have been space-qualified and operated on the QUESS satellite [17]. However, we are not aware of a true, heralded source of single photons that meets the requirements of repetition rate and pulse width presented here. The reason for this particular level of attenuation is the tradeoff between lowering multiphoton emission and maintaining a high number of total counts -this attenuation reduces the probability of multiphoton emission to 5% and results in a single photon rate of 100 MHz, which matches the resolution of the single photon detectors. The chosen wavelength guarantees high atmospheric transmission and the utilisation of well-developed telecommunication technology. The functional dependence of the fringe contrast on the Shapiro delay ∆τ and the pulse width √ σ (see [10], Eq. 13) demands the utilisation of ultra-short (<1 ps) single photon pulses to measure a noticeable effect (we note single photon sources with a width as short as 65 fs have been demonstrated in [18]). We take √ σ = 150 fs. Classical reference laser: A 1300 nm multi-purpose laser, operating in continuous wave mode, is led alongside the single photons through the interferometer. Its purpose is to provide reference data to estimate phase fluctuations and systematic errors. For precise corrections, an operational wavelength as close as possible to the single photon source is required, but overlap with the bandwidth of the single photon pulses has to be avoided. The chosen wavelength is a reasonable compromise. Frequency stability of both lasers is paramount; see Section 3.2. Additionally, a fraction of the incident laser power is separated at the ground station and used for wavefront reconstruction and compensation of polarisation changes due to satellite movement.
Fibres: We expect the delay fibres must be stabilised to relative length changes of 10 −10 for observing the predicted interference effects, which includes a thermal stabilisation of ±10 −5 K. Moreover, the fibre length must change dynamically to recover full contrast of the interference fringes. Active length corrections are carried out by a piezo fibre-stretcher. Fibres of length l = 60 km and refractive index n = 1.5 (glass) are assumed.
Transmission telescope: The onboard emitting telescope is used to focus the beam from the two sources. The aperture needs to be large enough to emit a strong signal but is limited by size and weight; a reasonable choice is a 60 cm aperture, 1 m long, 6 µrad field of view, Cassegrain reflector telescope. Material choices, such as Beryllium mirrors, would also minimise weight.
Single photon detectors: Since the goal of the detectors is to receive single photons that have travelled astronomical distances, it is extremely beneficial to minimise dark counts and maximise efficiency. We consider single-nanowire single photon detectors (SNSPDs), to benefit from the superior dark count over conventional avalanche photodiodes. While these detectors aren't as commonplace, the technology has already been established in a similar setting in the LLCD mission (NASA) [19].
Timing: Emission of the single photon pulses are tagged with timestamps by a Rb-clock (chosen primarily for its small size, weight and commercial availability). Synchronisation with an identical clock on the ground will be used to exclude background noise at the data post-processing stage and to track the path of the satellite. Precision timing is also necessary to modulate the action of time-gate filters in front of the detectors.

Random fluctuation and stability
As the effect we wish to measure is a minute change in optical path length in the interferometer due to the Shapiro delay, exposing the interferometer to other sources of instability can swamp the desired signal. We find that optical path and phase changes in the atmosphere due to temperature fluctuations are negligible compared to fluctuations in the length of fibre; given the thermal properties of fused silica we calculate necessary temperature stabilisation of the fibres of less than 10 −5 K to ensure a relative path length stability of 10 −10 . Passive insulation and active heat distribution can mitigate a large fraction of the thermal instability, but both the satellite and the ground station include a feedback loop of a frequency stabilised reference laser in combination with a piezo fibre stretcher, to allow for fibre noise reduction. Also, some thermal fluctuation can be erased in post-processing by referring to data from the reference laser.
In addition, continued operation of the experiment requires protection of the reference laser from frequency fluctuations and drift. Given the magnitude of the path length change, we estimate a necessary relative frequency stability of the reference laser less than 10 −11 . Long term accuracy of the reference laser within required precision is most feasible via frequency comb stabilisation [20]. This method will increase further complications and costs, and could be circumvented by improvements in the area of stabilisation by using atomic absorption lines, which today are close to reaching comparable relative accuracies [21].

Systematic transmission errors and dispersion
We highlight three systematic errors present in the experiment, mitigation of which are critical: dispersion, both in the fibre and the atmosphere; Doppler shift due to the velocity of the satellite; and changes in optical path length due to the ellipticity of the orbit.
Dispersion is prevalent both in the optical fibre and in the atmosphere. We estimate a requirement for fibre dispersion of < 5 fs/km/nm, ensuring a broadening of the pulse width in the fibres of < 0.5% per km; this is a stringent enough requirement to fix dispersion as the primary hurdle facing this experiment. The current state-of-the-art for dispersion-limited fibre is a factor of ten worse than this: 50 fs/km/nm [22]. Current technologies are capable of dispersion compensation in optical fibres of 0.5 ps/km/nm, which has to be improved by about a factor of ten to make the scientific requirement feasible. However, the utilisation of telecommunication fibre wavelengths ensure ongoing research in that scientific area. For example dispersion-free transmission of 610 fs laser pulses over a 160 km fibre has been demonstrated, but dispersion-free transmission techniques typically suffer from large losses [23] [24]. On the other hand active compensation of atmospheric chromatic dispersion of a 250 fs pulse over a 200 km propagation distance, to a uncertainty of ±10 fs, has been demonstrated [25]. This indicates that control of atmospheric dispersion of comparably short pulses is well within reach, especially considering the pulses in this experiment are travelling vertically, on a much shorter trajectory through the atmosphere.
Doppler shift, conversely, presents much less of a problem. This is straightforward to calculate based on the velocity of the satellite relative to the Earth. For the worst possible case (at apogee), the Doppler redshift ν is calculated using the formula: where c is the speed of light, f p is the frequency of the single photon source and v 0 is the maximum velocity of the satellite at the lowest sampling point. This gives a shift in photon frequency of 2.48 GHz. Considering the bandwidth size and that the photon frequency exceeds 193 THz, this is considered negligible. Finally, the radial motion of the satellite causes constant variation in optical path length. To keep the path lengths of the interferometer equal this variation must be continuously compensated for. The time spent in the fibre loop by the part of the photon state in the upper arm, before it is emitted from the satellite, is given by: where = 60 km is the fibre length and n = 1.5 is the refractive index. Multiplying this quantity by the radial velocity gives the change in path length between parts of the superposition due to the motion of the satellite; this change in path is plotted in Figure 5. However, this effect can be precisely calculated in advance and so either active compensation with optical components or passive compensation at the post-processing stage can be built into the experiment beforehand.

Loss
Given the repetition rate of the laser and strength of attenuation, we can calculate the expected rate of signal photons registering at the detectors. The analogous discussion on noise photons reaching the detectors is deferred to Section 3.5. There are three major sources of signal attenuation: divergence of the beam from the transmitting telescope to the ground; loss due to atmospheric irradiance and interference; and loss within the fibres. First, we examine the effects of atmospheric transmission. A wavelength of 1550nm is deemed suitable as it is standardised for telecoms use, maximises transmittivity in the fibre and is not readily absorbed by the atmosphere. The high frequency of the emitted photon allows us to ignore ionospheric effects on polarisation and is easily distinguishable from auroral activity over the ground stations [26]. Assuming homogeneity, using a variant of the Beer-Lambert Law we can estimate the probability of a single photon passing through the atmosphere as: Where τ 0 is the optical depth of the atmosphere, and η 0 is the angle of incidence of the beam. The optical depth varies over time and depends on myriad dynamic factors such as aerosol content, atmospheric mixing and Raman and Rayleigh scattering. Modelling these effects is essential, so one might consider optical depth readings from MODIS, MISR or future Sentinel satellites. A numerical weather prediction model might then be produced detailing daily optical depth over ground stations. An alternative might be to model turbulence transfer functions as in [27]. As shown in Figure 6, we can allow for an optical depth of 0.5 before the satellite signal is severely attenuated. Combining these factors produces an estimate for the atmospheric transmission loss to be −2.3 dB. As for beam divergence, a downlink direction is chosen as accentuated divergence would then occur only in the last 12 km of travel. This downlink beam divergence can be computed for a photon of wavelength λ by: where w 0 = 30 cm corresponds to the radius of the transmitting telescope on the satellite. Consequently, the beam diameter on the ground at maximum and minimum altitude is given by D ground 2hθ, which implies that 58 m ≤ D ground ≤ 105 m, taking bounding values for the altitude as 18000 m ≤ h ≤ 32000 m.
Atmospheric turbulence must also be considered. As the atmospheric parameters used for post-processing are extracted from data from the reference laser, this turbulence must not radically change in the time between measurement of the reference laser and of the single photon source. Assuming that the gap between reference and single photon measurement is half the repetition rate of the laser, this delay comes to 5 ns. Assuming a wind speed of 10 m s −1 , an air parcel would move by about 50 nm in some direction between the single photon measurement and reference laser measurement. This is considered negligible compared to other atmospheric effects.
Atmospheric transmission loss, turbulence and beam divergence can be compiled into a "link budget" that calculates the full transmission loss from satellite to ground station. We can modify and simplify the Friis Transmission Equation to give the following link budget equation [28]: where here D T is the transmitter diameter, D R is the receiver diameter, λ is the wavelength of the single photon source, h is the satellite altitude, L p is the pointing loss (taken to be 0.63 [29]), and L t the atmospheric transmission loss as calculated above.
Assuming an apogee of 32000 km and a lowest altitude bin of 18000 km, a ground receiver diameter of 3 m, a satellite transmitter diameter of 0.6 m, and superconducting detector efficiencies of about 90% [30]; we calculate the baseline signal gain from the link budget to be -30 dB at perigee, and -35 dB at apogee.
Further factoring in losses from the optical fibre of -15.5 dB and fibre blackening from radiation of -1.7 dB, the link budget achieves a total attenuation of -52.2 dB in the worst case, towards the end of mission lifetime. Of course, to further reduce signal attenuation, one could increase the aperture diameters of the receiver and transmitter. However, this causes a loss in manoeuvrability and a significant cost increase for rapidly decreasing returns.

Noise
A feasibility case must also ensure a signal-to-noise ratio (SNR) great enough to produce enough meaningful data for statistical analysis. It must be stressed that the signal received on Earth is fixed by the link budget above, however since we are trying to detect a single photon from a plethora of solar and planetary photon noise, optimising the SNR is crucial.
The effects of noise on space to ground quantum channels has previously been explored [31]. The noise power received by the ground telescope (P b ) can be expressed as where H b is the brightness of the sky in units of W m −2 sr −1 µm −1 , Ω fov is the field of view, ∆ν is the bandwidth and ∆t d is the detection time. Typical sky brightness for quantum cryptography applications is discussed in [31], using data from [32]. However, the data therein must be modified in light of the experiment proposed here. Firstly, the data presented in [31,32] is for a frequency band just below 1550 nm, which is subject to much less noise than at 1550 nm itself. Conversely, the primary source of noise photons at 1550 nm is hydroxyl airglow, the strength of which is strongly dependent on ambient temperature. The data from [32] assumes a receiving station in the tropics at Mauna Kea, whereas we propose an arctic station at Svalbard. Utilising instead polar sky brightness data from [33], we take a sky brightness of 2×10 −5 W m −2 sr −1 µm −1 at our operational frequency. The received signal band is filtered to a bandwidth of 15 nm, corresponding to twice the bandwidth of the single photon pulses (twice the full width at half maximum of the Lorentzian pulse). The field of view of the receiving telescope is assumed to be 10 µrad, but the effective field of view is further reduced by a factor of 10 with a variable iris diaphragm [34]. Noise power is further reduced with a 50 ps time gate filter leading to an available detection time of 50 ms per second. Reflections from the satellite or its black body radiation do not significantly contribute to the background noise [31]. Besides received background photons, total noise power also depends on the dark count rate of the single photon detectors. Superconducting detectors have negligible intrinsic dark count rate but as a worst case estimate the detector system dark count rate is assumed to be 1 kHz (although a large fraction of these counts are neglected due to time gating) [30]. The rate of detected signal photons is calculated as 590 Hz (using the repetition rate of the laser, attenuation and loss). Combining this figure with the expected noise from ambient photons and system dark count gives an SNR ≈ 9.0 =9.6 dB.

Hypothesis testing
Once an experimental profile has been outlined, we can perform a statistical analysis to extract the confidence with which we can confirm, or refute, the hypothesis in Eq. 2 (and likewise, the theoretical analysis presented in [10]).
We assume data taken in k altitude bins, uniformly spaced in h from h min = 18000 km to h max = 32000 km, with n counts per bin. Then, we take normally distributed experimental parameters with mean as given in the specification in Section 3.1, and variance as fixed by the discussion on stability in Section 3.2. We assume that a fraction of the photons are lost according to the link budget in Section 3.4, and that there is an ambient background of noise photons as per Section 3.5. We then perform Monte Carlo simulations for a range of number of altitude bins, k, and counts per altitude bin, n, to generate artificial data sets given this initial prescription.
To extract a statistical significance from these simulations, we perform a Kolmogorov-Smirnov statistical significance test between the simulated data set and the mathematical formulation in Eq. 2 and [10]. The test is devised to extract the goodness of fit for a data set to an arbitrary functional hypothesis. Using this test, we can predict the number of single photons that must be sent in order to test the hypothesis to a specified confidence.
The consequential data from the simulation and statistical test is shown in Figure 7, for k = 500 altitude bins. Here, we can see that the number of received, true counts per altitude bin must exceed ≈ 6.5 × 10 5 counts in order to test the hypothesis to a 5σ confidence. Given the figure for total loss of −52.2 dB, this amounts to a necessary emission total of ≈ 10 11 single photons per altitude bin. Given the repetition rate of the laser and the degree of attenuation, we expect raw emission rates of 100 MHz and thus reach the required emission total in each altitude bin within ∼ 1000 seconds = 16.7 minutes of continuous operation.

Mission Design
In order to reach the desired confidence defined in Section 4, data must be taken over a sufficiently large difference in gravitational potential. An orbit with a perigee of 700 km and an apogee of 32000 km gives access a relativistic time delay of 150 fs, which is both large enough to be resolved by the detectors and still gives a reasonable count rate at the apogee. During the measurement procedure, every other light source introduces noise. The ambient sunlight is indeed strong enough to wash out the signal from the satellite. Therefore, measurements must be performed when the ground station is not illuminated by the Sun. Additionally, in order to have the maximum number of measurements per orbit and to minimise noise photons from airglow, the ground station needs to be placed as close as possible to the North Pole.
These considerations lead to the selection of the ground station located in Svalbard, Norway. By considering the maximum optical path that gives a valid measurement and the movement capabilities of the ground telescope, a connection cone of 45°around the zenith can be defined. This indeed allows satellite access times for up to 7 hours per orbit. This ground station is ideal for the mission, though it introduces a constraint on orbital inclination. In fact, only a polar orbit can provide the maximum visibility time from the ground.
The range of orbital heights chosen to send single photons is from 18000 km to 32000 km (to minimise special relativistic effects, as discussed in Section 2.1). Given the number of altitude bins and length of measurement windows, we estimate that a mission lifetime of 1.5 years will completely satisfy the previously stated requirements. The mission profile is thus composed of three different phases. The first, starting immediately after launch, is a 6-month commissioning period used to calibrate the orbit and the measurement system with laser ranging and radio communication. During the satellite motion, due to the perturbation from the oblateness of the Earth and the high eccentricity of the orbit, the orbital apse line rotates clockwise with a rate of 77.1°/year. With an initial argument of perigee of ω = 350°and a launch during Svalbard's spring season, after the commissioning phase ω 310°, meaning that the apogee is contained in the ground station visibility cone. This situation is represented by the thicker line in Figure 8. The main operational phase thus begins in winter and lasts for approximately 6 months. It is then followed by a second mission phase in Svalbard's summer, during which the first results are examined and the orbit is further calibrated. Finally, the last 6-month phase provides additional measurements. At the end of the 1.5-year lifetime, a 55 ms −1 thruster burn lowers the perigee to approximately 200 km. This manoeuvre eventually leads to a further lowering of the apogee and then to a controlled de-orbit into the atmosphere.
The total satellite mass is approximated at 400 kg. The small spacecraft size allows using the VEGA launcher in order to bring the satellite into a parking orbit with a perigee of 700 km and an apogee of 20000 km. The final orbit is then reached using an on board bi-propellant propulsion system through a perigee burn with a ∆V of 323 ms −1 .

Risk Analysis
Here we discuss risk control and mitigation during the mission.

Radiation effects
The spacecraft's immediate radiation environment, composed of fluxes from the solar wind and galactic cosmic rays, has a non-negligible effect on the performance of the components. During its orbit around Earth, the satellite passes twice through the Van Allen belts. Therefore, suitable shielding is necessary to ensure the correct operation of the payload. Particular precaution has to be taken to shield the spools of optical fibre, due to the effect of Radiation-Induced Attenuation (RIA). Single mode fibres with an RIA of 0.5 dB/km per 9000 rad have been demonstrated in the literature [35], a figure which can be used for a worst-case analysis. Assuming a 20 mm thick spherical layer of aluminium shielding, the fibre receives an estimated Total Ionizing Dose (TID) between 500 and 2000 rad/year, which yields an RIA between -1.7 and -6.7 dB/year. An aluminium shielding of 2 mm is provided for the lasers and the optical bench. This is also ample shielding to guarantee low radiation doses for the remaining optics.
The blackening of fibre cables due to radiation is thus one of the main operational risks. Although this process is slowed by adequate shielding, occurrence at a higher than expected rate would pose a major obstacle to collecting adequate data. Active avoidance of the Van Allen belts by modifying the orbit could be pursued, if needed.

Thermal Control System
During the mission, the satellite undergoes cyclic eclipses; the resulting temperature fluctuations have to be analysed in order to design the thermal control system. The most critical requirement is the temperature stability of the optical bench. Insulation ensures that the fibre temperature remains within 0.5 K of the equilibrium value. Specifically, the satellite is equipped with a Multi Layer Insulation (MLI) aluminium and Kapton coating and heat pipes to convey and redistribute heat. With the appropriate sizing for radiators and heat pipes it is possible to achieve the required thermal stability, although this represents another primary risk for the mission.

Attitude and Orbit Control System
Precise determination of the satellite's orbit and fine pointing are also fundamental.
Firstly, we assumed in Sec 3.4 that the satellite's pointing accuracy is 0.5 µrad, with 0.1 mrad/s maximum slew rate. This strict requirement ensures efficient transmission of the required number of photons to the ground, but is a relevant risk for the mission. Indeed, if the satellite failed on maintaining the required accuracy, a prohibitively large fraction of the signal photons would be lost. The satellite is thus equipped with a tracking telescope to have a first estimate of the attitude; then, two star trackers refine the determination. A system of reaction wheels and vibration dampers allows for a fine control of the satellite orientation.
As mentioned above, orbital determination is of key importance for the mission. Indeed, in order to account for the transmission error due to the radial motion of the satellite, the radial velocity must be known to within ∼ 1 mm/s precision. This is feasible with current radio tracking systems working at either S-or X-band. In fact, the commonly used Ultra-Stable Oscillators often have a Van Allen stability better than 10 −13 over an integration time from 10 s to 1000 s, which is completely adequate to fulfil the requirement.

Conclusion
We have demonstrated herein a proposal capable of probing the interplay between quantum mechanics and general relativity using single photon interferometry -in particular, exploring the dichotomous nature of time in the two theories. Moreover, the apparatus is well within the reach of current quantum optics technology, and there is already historical precedent for space-qualification and launch of similar payloads. Conversely, successful operation of this experiment would provide strong precedent itself for future projects involving commercial quantum communications satellites. SP was supported by the Bristol Quantum Engineering Centre for Doctoral Training, EPSRC grant EP/L015730/1. PK was supported by the UK EPSRC through the Quantum Technology Hub: Networked Quantum Information Technology (grant reference EP/M013243/1). NG was supported by the Groupement Luxembourgeois de l'Aronautique et de l'Espace (grant GLAE-Alpbach). This work was carried out primarily at the Alpbach and Post-Alpbach conferences. The Summer School Alpbach: Sixty highly qualified European science and engineering students converge annually for stimulating ten days of work in the Austrian Alps. Four teams are formed who each design a space mission which are then judged by a jury of experts. Students learn how to approach the design of a satellite mission and explore new and startling ideas supported by experts. The Summer School Alpbach enjoys 40 years of tradition in providing in-depth teaching on different topics of space science and technology, featuring lectures and concentrated working sessions on mission studies in self-organised working groups. The Summer School Alpbach is organised by FFG and co-sponsored by ESA, ISSI and the national space authorities of ESA member and cooperating states. We would like to specifically acknowledge the tutelage of Günter Kargl, Manuela Werner, Siddharth Koduru Joshi, Esteban Castro-Ruiz and Rupert Ursin, whose input was instrumental.
Finally, for an elliptical orbit with semi-major axis a, indicating with h the altitude of the satellite on the Earth's surface, the time delay of Eq. (13) becomes: which is exactly Eq. (1) in the text.