Precision Measurement of the Newtonian Gravitational Constant Using Cold Atoms G. Rosi1, F. Sorrentino1, L. Cacciapuoti2, M. Prevedelli 3 & G. M. Tino1 1Dipartimento di Fisica e Astronomia and LENS, Universita` di Firenze - INFN Sezione di Firenze, Via Sansone 1, 50019 Sesto Fiorentino, Italy 2European Space Agency, Keplerlaan 1 - P.O. Box 299, 2200 AG Noordwijk ZH, The Netherlands 3Dipartimento di Fisica, Universita` di Bologna, Via Irnerio 46, 40126, Bologna, Italy About 300 experiments tried to determine the value of the Newtonian gravitational constant G to date but large discrepancies in the results prevent from knowing its value precisely1. The weakness of the gravitational interaction and the impossibility of shielding the effects of gravity make it very difficult to measure G keeping systematic effects under control. Most of the experiments performed so far were based on the torsion pendulum or torsion balance scheme as in the experiment by Cavendish2 in 1798 and in all cases they used macroscopic masses. Here we report on the first precise determination of G using laser-cooled atoms and quantum interferometry to probe gravity. We obtain the value G= 6.67191(99) × 10−11 m3 kg−1 s−2 with a relative uncertainty of 150 ppm. Our value is at 1.5 combined standard deviations from the current recommended value of the Committee on Data for Science and Technology (CODATA)3. Such a conceptually different experiment is important to identify the so far elusive systematic errors thus improving the confidence in the value of G. There is no definitive relationship indeed between G and the other fundamental constants and no theoretical pre- 1 diction for its value to test the experimental results. Improving the knowledge of G has not only a pure metrological interest, but is also important for the key role that this fundamental constant plays in theories of gravitation, cosmology, particle physics, astrophysics, and geophysical models. The basic idea of our experiment is to use an atom interferometer as gravity sensor and a well-characterized mass as source of a gravitational field. From the precise measurement of the atoms’ acceleration produced by the source mass and from the knowledge of the mass distribution, it is possible to extract the value of the gravitational constant G using the well-known formula F(r) = −G M1M2ˆr (1) r2 Atom interferometers4,5 are new tools for experimental gravitation as, for example, in precision measurements of gravity acceleration6 and gravity gradients7, as gyroscopes based on the Sagnac effect8, for testing the 1/r2 law9, general relativity10 and quantum gravity models11, and for applications in geophysics12. Proof-of-principle experiments to measure G using atom interferometry were reported13–15. Ongoing studies show that future experiments in space will allow to take full advantage of the potential sensitivity of atom interferometers for fundamental physics tests16. The possibility of using atom interferometry for gravitational waves detection is being studied17. Since the problem in the determination of G depends on the presence of unidentified systematic errors, our experiment was designed with a double-differential configuration in order to be as immune as possible from such effects: the atomic sensor was a double interferometer in a gravity- 2 gradiometer configuration, in order to subtract common-mode spurious signals, and we used two sets of well-characterized tungsten masses as source of the gravitational field that were placed in two different positions in order to modulate the relevant gravitational signal. An additional cancellation of common-mode spurious effects was obtained by reversing the direction of the two-photon recoil used to split and recombine the wave packets in the interferometer18. Efforts were devoted to the control of systematics related to atomic trajectories, positioning of the atoms, and effects due to stray fields. The high density of tungsten was instrumental to maximize the signal and to compensate the Earth’s gravity gradient in the atom interferometers regions thus reducing the sensitivity to the vertical position and size of the atomic probes. The atom interferometer is realized using light pulses to stimulate 87Rb atoms on the twophoton Raman transition between the hyperfine levels F =1 and F =2 of the ground state19. The light field is generated by two counter-propagating laser beams with wave vectors k1 and k2 ≃ −k1 aligned along the vertical direction. The gravity gradiometer consists of two vertically separated atom interferometers operated in differential mode. Two atomic clouds launched along the vertical direction are simultaneously interrogated by the same π/2-π-π/2 pulse sequence. The difference of the phase shifts detected at the output of each interferometer provides a direct measurement of the differential acceleration induced by gravity on the two atomic samples. In this way, any spurious acceleration induced by vibrations or seismic noise on the common reference frame identified by the vertical Raman beams is efficiently rejected. Fig. 1 shows a schematic of the experiment. The atom interferometer apparatus and the 3 source masses assembly are described in detail elsewhere20,21. In the vacuum chamber at the bottom of the apparatus, a magneto-optical trap (MOT) collects ∼ 109 rubidium atoms. After turning the MOT magnetic field off, the atoms are launched vertically along the symmetry axis of the vacuum tube by using the moving molasses technique. During the launch sequence, atoms are laser cooled to a temperature of ∼ 4 µK. We juggle two atomic samples to have them reach the apogees of their ballistic trajectories at about 60 cm and 90 cm above the MOT with a vertical separation of 328 mm. The atoms are velocity selected and prepared in the magnetic-field-insensitive F =1, mF =0 hyperfine state with a combination of three π Raman pulses and two resonant laser pulses used to blow away the atoms occupying the wrong hyperfine state. The interferometers are realized at the center of the vertical tube shown in Fig. 1. In this region, surrounded by two cylindrical magnetic shields, a uniform magnetic field of 29 µT along the vertical direction defines the quantization axis. Here, atoms are subjected to the Raman three-pulse interferometer sequence. The central π pulse occurs about 6 ms before the atoms reach the apogees of their trajectories. At the end of their ballistic flight, the population of the ground state is measured by selectively exciting the atoms in both hyperfine levels of the ground state and detecting the light-induced fluorescence emission. We typically detect 105 atoms on each rubidium sample at the end of the interferometer sequence. Each measurement takes 1.9 s. The information on the relative phase shift between the two atom interferometers is obtained from a Lissajous plot of the signal of one interferometer vs the other. Experimental points distribute along an ellipse. The differential phase shift is extracted from the eccentricity and the rotation angle of the ellipse fitting the data22. The instrument sensitivity for 4 differential acceleration is 3 × 10−9 g for 1 s of integration. The source mass is composed of 24 tungsten alloy cylinders, for a total mass of about 516 kg. Each cylinder is machined to a diameter of 99.90 mm and a height of 150.11 mm. They are positioned on two titanium platforms and distributed with hexagonal symmetry around the vertical axis of the tube (see Fig. 1). The cylinders’ centers lie along two circles with nominal radii 2R √ and 2R 3 respectively, where R is the radius of a single cylinder. The vertical positioning of the two platforms is ensured by precision screws synchronously driven by stepper motors and by an optical readout system. The reproducibility of the positioning was verified with a laser tracker to be within 1 µm20. With respect to the position of the apogee of the lower atomic cloud, the centers of the lower and upper set of cylinders lie at a vertical distance of 40 mm and 261 mm respectively for the C configuration and -74 mm and 377 mm respectively for the F configuration. The value of the Newtonian gravitational constant was obtained from a series of gravity gradient measurements performed by periodically changing the vertical position of the source masses between configuration C and F . Fig. 2 shows the data used for the determination of G. Data were collected in 100 hours during one week in July 2013. Each phase measurement was obtained by fitting a 360-point scan of the atom interference fringes to an ellipse. The modulation of the differential phase shift produced by the source mass is well visible and can be resolved with a signal-to-noise ratio of 1000 after about one hour. The resulting value of the differential phase shift is 0.547870(63) rad from which the G value was obtained. The measured differential phase shift is produced for 97.0% by the cylinders, for 2.8% by the cylinders supports20, and for 0.2% by 5 the additional moving masses (translation stages, optical rulers, screws). A budget of the uncertainty sources affecting the value of G is presented in Table 1. Positioning errors account for uncertainties in the positions of the 24 tungsten cylinders along the radial and vertical direction, both in configuration C and F . Density inhomogeneities of the source masses were measured by cutting and weighing a spare cylinder20 and modelled in the data analysis. The precise knowledge of the atomic trajectories is of key importance to analyze the experimental results and derive the value of G. The velocity of the atomic clouds and their position at the time of the first interferometer pulse were calibrated by time-of-flight measurements and by detecting the atoms when they crossed a horizontal light sheet while moving upward and downward. The Earth rotation affects the atom interferometer signal because of the transverse velocity distribution of the atoms in the fountain. Following the method demonstrated for a single interferometer23,24, we implemented a tip-tilt scheme for the mirror retroreflecting the Raman beams in our double interferometer. The analysis to extract the value of G from the data involved the following steps: (1) calculation of the gravitational potential produced by the source masses, (2) calculation of the phase shift for single atom trajectories, (3) Monte Carlo simulation of the atomic cloud, (4) calculation of the corrections for the effects not included in the Monte Carlo simulation which are listed in Table 1. After an analysis of the error sources affecting our measurement, we obtain the value G = 6.67191(77)(65)×10−11 m3 kg−1 s−2. The statistical and systematic errors, reported in parenthesis as one standard deviation, lead to a combined relative uncertainty of 150 ppm. In Fig. 3, this 6 result is compared with the values of recent experiments and CODATA adjustments. Our value, obtained with a method completely different from the ones used before, is at 1.5 combined standard deviations from the current CODATA value G = 6.67384(80) × 10−11 m3 kg−1 s−2. Such conceptually different experiments will be important to identify errors that produced the still unsolved discrepancies amongst previous measurements. Given the relevance of the gravitational constant in several fields ranging from cosmology to particle physics and in the absence of a complete theory linking gravity to other forces, high precision measurements based on different methods are crucial to improve the confidence in the value of this fundamental constant. It is possible to envisage a new generation of experiments to measure G using ultracold atoms confined in optical traps for their precise positioning. The remaining major contribution to the systematic uncertainty in our experiment derives indeed from the positioning of the atoms with respect to the source mass. The choice of a different atom might also reduce other sources of uncertainties. For example, Sr atoms can be rapidly cooled down to Bose Einstein condensation25, they are virtually immune to magnetic fields and collisional effects in their ground state, and the possibility of precise positioning and gravity measurements using optical lattices was demonstrated26,27. Higher sensitivity atom interferometers with larger splitting of the atomic wave-packet would enable the use of a smaller highly homogeneous source mass such as gold or eventually silicon crystals. We foresee prospects for pushing the measurement accuracy in the determination of the gravitational constant to better than 10 ppm using atoms. 7 Methods Summary The Raman lasers frequencies ν1 and ν2 match the resonance condition ν2 − ν1 ≃ 6.8 GHz, corresponding to the separation between the F =1 and F =2 levels of 87Rb ground state. During absorption and stimulated emission processes, the resonant light field ex- changes with atoms a total momentum of h¯k = h¯k1 − h¯k2, coupling the two states |F = 1, p and |F = 2, p + h¯k , where p is the initial atomic momentum. The interferometer is composed of a π/2-π-π/2 sequence of Raman pulses separated by a time T = 160 ms. The pulses split, redirect and recombine the wave-packets producing atom interference. At the end of the sequence, the probability of detecting atoms in |F = 2, p + h¯k is P2 = 1 2 [1 − cos φ], where φ is the phase difference accumulated along the interferometer arms. In a uniform gravity field, atoms experi- ence a phase shift φ = k · gT 2 that can be measured thus providing direct information on the local acceleration due to gravity. Accurate control of noise sources and systematic phase shifts is crucial to optimize the sensitivity and the long-term stability of gravity gradient measurements28. Active loops have been implemented to stabilize the optical intensity of cooling, Raman and probe laser beams and the Raman mirror tilt. In this way, we reach a quantum-projection-noise-limited sensi- tivity to differential accelerations of 3 × 10−9 g after 1 s of integration time and an accuracy on the G measurement at the level of 100 ppm. The measured phase shift is compared to the theoretical value obtained by evaluating the action integral along the trajectories given by the solution of the Lagrange equation: φ = 1 h¯ Ldt. (2) Following a perturbative approach29, we separate the Lagrangian L into L0 = p2/(2m) − mgz, m being the atomic mass, and L1 accounting for Earth gravity gradient and source masses contribu- 8 tion. A Monte Carlo simulation is used to evaluate the measurement results and to determine the sensitivity of the final phase angle to the relevant parameters. To this purpose, the key parameters entering the simulation are varied and their derivatives calculated to estimate the uncertainty on the G measurement. 1. Quinn, T. Measuring big G. Nature 408, 919–921 (2000). 2. Cavendish, H. Experiments to determine the density of the Earth. Philosophical Transactions of the Royal Society of London 88, 469–526 (1798). 3. Mohr, P. J., Taylor, B. N. & Newell, D. B. CODATA recommended values of the fundamental physical constants: 2010. Rev. Mod. Phys. 84, 1527–1605 (2012). 4. Cronin, A. D., Schmiedmayer, J. & Pritchard, D. E. Optics and interferometry with atoms and molecules. Rev. Mod. Phys. 81, 1051–1129 (2009). 5. Tino, G. M. & Kasevich, M. A. (eds.). Atom Interferometry, Proceedings of the International School of Physics ”Enrico Fermi”, Course CLXXXVIII, Varenna 2013 (Societa` Italiana di Fisica and IOS Press). 6. Peters, A., Chung, K. Y. & Chu, S. Measurement of gravitational acceleration by dropping atoms. Nature 400, 849–852 (1999). 7. McGuirk, J. M., Foster, G. T., Fixler, J. B., Snadden, M. J. & Kasevich, M. A. Sensitive absolute-gravity gradiometry using atom interferometry. Phys. Rev. A 65, 033608 (2002). 9 8. Gustavson, T. L., Landragin, A. & Kasevich, M. A. Rotation sensing with a dual atom interferometer Sagnac gyroscope. Class. Quantum Grav. 17, 2385–2398 (2000). 9. Ferrari, G., Poli, N., Sorrentino, F. & Tino, G. M. Long-Lived Bloch Oscillations with Bosonic Sr Atoms and Application to Gravity Measurement at the Micrometer Scale. Phys. Rev. Lett. 97, 060402 (2006). 10. Dimopoulos, S., Graham, P., Hogan, J. & Kasevich, M. Testing General Relativity with Atom interferometry. Phys. Rev. Lett. 98, 111102 (2007). 11. Amelino-Camelia, G., La¨mmerzahl, C., Mercati, F. & Tino, G. M. Constraining the energymomentum dispersion relation with Planck-scale sensitivity using cold atoms. Phys. Rev. Lett. 103, 171302 (2009). 12. de Angelis, M. et al. Precision gravimetry with atomic sensors. Meas. Sci. Technol. 20, 022001 (2009). 13. Bertoldi, A. et al. Atom interferometry gravity-gradiometer for the determination of the Newtonian gravitational constant G . Eur. Phys. J. D 40, 271–279 (2006). 14. Fixler, J. B., Foster, G. T., McGuirk, J. M. & Kasevich, M. Atom Interferometer Measurement of the Newtonian Constant of Gravity. Science 315, 74–77 (2007). 15. Lamporesi, G., Bertoldi, A., Cacciapuoti, L., Prevedelli, M. & Tino, G. M. Determination of the Newtonian gravitational constant using atom interferometry. Phys. Rev. Lett. 100, 050801 (2008). 10 16. Tino, G. M. et al. Precision Gravity Tests with Atom Interferometry in Space. Nucl. Phys. B (Proc. Suppl.) 243-244, 203–217 (2013). 17. Special issue on “Gravitational waves detection with atom interferometry”. In Tino, G. M., Vetrano, F. & La¨mmerzahl, C. (eds.) Gen. Relativ. Gravit., vol. 43, 1901–2088 (2011). 18. Louchet-Chauvet, A. et al. The influence of transverse motion within an atomic gravimeter. New J. Phys. 13, 065025 (2011). 19. Kasevich, M. & Chu, S. Atomic interferometry using stimulated Raman transitions. Phys. Rev. Lett. 67, 181–184 (1991). 20. Lamporesi, G. et al. Source masses and positioning system for an accurate measurement of G. Rev. Sci. Instrum. 78, 075109 (2007). 21. Sorrentino, F. et al. Sensitive gravity-gradiometry with atom interferometry: progress towards an improved determination of the gravitational constant. New J. Phys. 12, 095009 (2010). 22. Foster, G. T., Fixler, J. B., McGuirk, J. M. & Kasevich, M. A. Novel Method of Phase Extraction Between Coupled Atom Interferometers Using Ellipse-Specific Fitting. Opt. Lett. 27, 951–953 (2002). 23. Hogan, J., Johnson, D. & Kasevich, M. Light-pulse atom interferometry. In Arimondo, E., Ertmer, W., Schleich, W. & Rasel, E. (eds.) Proceedings of the International School of Physics ”Enrico Fermi”, Course CLXVIII, 411–447 (SIF, Bologna and IOS Press, Amsterdam, 2009). 11 24. Lan, S.-Y., Kuan, P.-C., Estey, B., Haslinger, P. & Mu¨ller, H. Influence of the Coriolis force in atom interferometry. Phys. Rev. Lett. 108, 090402 (2012). 25. Stellmer, S., Pasquiou, B., Grimm, R. & Schreck, F. Laser cooling to quantum degeneracy. Phys. Rev. Lett. 110, 263003 (2013). 26. Sorrentino, F. et al. Quantum sensor for atom-surface interactions below 10µm. Phys. Rev. A 79, 013409 (2009). 27. Poli, N. et al. Precision measurement of gravity with cold atoms in an optical lattice and comparison with a classical gravimeter. Phys. Rev. Lett. 106, 038501 (2011). 28. Sorrentino, F. et al. Sensitivity limits of a Raman atom interferometer as a gravity gradiometer. Phys. Rev. A 89, 023607 (2014). 29. Storey, P. & Cohen-Tannoudji, C. The Feynman path integral approach to atomic interferometry. A tutorial. J. Phys. II France 4, 1999–2027 (1994). Acknowledgements G.M.T. acknowledges seminal discussions with M.A. Kasevich, and J. Faller and useful suggestions by A. Peters in the initial phase of the experiment. We are grateful to A. Cecchetti and B. Dulach of INFN-LNF for the design of the source mass support and to A. Peuto, A. Malengo, and S. Pettorruso of INRIM for density tests on W mass. We thank D. Wiersma for a critical reading of the manuscript. This work was supported by INFN (MAGIA experiment). Author Contributions G. M. T. conceived the experiment, supervised it and wrote the manuscript. G. R., F. S. and L. C. performed the experiment. M. P. contributed to the experiment and analysed the data. 12 Author Information The authors declare that they have no competing financial interests. Correspondence and requests for materials should be addressed to G. M. T. (guglielmo.tino@fi.infn.it). Methods Experimental setup Our interferometer uses Raman pulses to drive 87Rb atoms on the two-photon transition between the hyperfine levels |1 ≡ |F = 1, mF = 0 and |2 ≡ |F = 2, mF = 0 of the ground state19. The light field is generated by two counter-propagating laser beams with wave vectors k1 and k2 ≃ −k1 aligned along the vertical direction. This configuration is obtained by retro-reflecting the Raman lasers by a mirror placed on top of the vertical tube (see Fig. 1). The laser frequencies ν1 and ν2 match the resonance condition ν2 − ν1 = ν0 + νAC + νrec, where hν0 is the energy corresponding to the |1 → |2 transition, hνAC accounts for the AC Stark effect and hνrec is the energy shift due to the two-photon recoil. As a consequence, the transition modifies both the internal energy and the total momentum of the atom, coupling the two atomic states |1, p and |2, p + h¯k , where k = k1 − k2 is the effective wave vector of the Raman transition. The h¯k momentum transfer is responsible for the spatial separation of the atomic wave-packet 13 along the two physical arms of the interferometer. During the interferometer, a sequence of three Raman pulses provides the atom-optical elements of a typical Mach-Zehnder configuration: at t = 0 a π/2 beam splitter pulse prepares the atoms, initially in |1, p , in an equal and coherent superposition of |1, p and |2, p + h¯k ; after a time T = 160 ms a π pulse plays the role of a mirror, swapping state |1, p with |2, p + h¯k and redirecting the atoms towards the output ports of the interferometer; at t = 2T a final π/2 pulse recombines the atomic wave-packets in the two output ports where they interfere. When interacting with the light field, the phase of the Raman lasers is imprinted on the atomic wave-function, sampling the atomic motion at the three interaction events of the π/2-π-π/2 sequence. This information can be read through the effect of matter-waves interference by detecting the atomic population in the F = 1 and F = 2 hyperfine states. If φ is the phase difference accumulated by the atoms along the two interferometer arms, the probability of detecting atoms in |2 can be expressed as P2 = 1 2 [1 − cos φ]. In the presence of a gravity field, atoms experience a phase shift φ = k · gT 2. Measuring φ is therefore equivalent to a measurement of the local acceleration due to the gravity field along the direction of the effective wave vector k. Our instrument is a gravity gradiometer consisting of two Raman-pulse interferometers simultaneously probing two clouds of laser-cooled 87Rb atoms aligned along the vertical axis and separated by a distance of 328 mm 21,28. The two atomic clouds are interrogated by the same Raman lasers with the same interferometric pulse sequence. The difference of the phase shifts at each interferometer provides a direct measurement of the differential acceleration induced by gravity on the two atomic samples. This configuration allows a very efficient rejection of common mode noise sources including spurious accelerations at the instrument platform induced by mechanical 14 vibrations or seismic noise, g variations due to tidal effects and phase noise introduced by the Raman lasers. In addition, the difference of the measurements performed in configurations C and F (see Fig. 1) efficiently rejects long-term systematic drifts which do not depend on the distribution of the source masses. They include wavefront distortions, magnetic field perturbations, Coriolis acceleration and light shift. Controlling sources of statistical noise and systematic phase shifts on the gravity gradient measurement is crucial for the experiment28. Inertial effects induced by Earth’s rotation are an important source of both noise and systematic errors. Gravity gradient measurements are indeed sensitive to the Coriolis acceleration when differential velocities along the East-West direction are present between the two atomic samples and when the launch direction itself is not stable. Precise control on the effects of the Coriolis acceleration is achieved by applying a uniform rotation rate to the retro-reflecting Raman mirror by means of PZT actuators during the atom interferometry sequence23,24. The optimal rotation rate, counteracting the local projection of the Earth rotation rate on the horizontal plane, is determined to better than 2 µrad/s after maximizing the interferometers contrast and minimizing the noise on the differential phase measurements. In such conditions, the residual shift amounts to 1.5·∆θ, where ∆θ is the variation in the launch direction expressed in rad when switching from configuration C to F and viceversa. As a consequence, a G measurement down to 100 ppm level requires a control on the launch direction to better than 37 µrad. We can control ∆θ to better than 8 µrad. The velocity distribution of the atoms due to the finite temperature of the sample (∼ 4 µK) also contributes a systematic shift on the G measurement. In our experiment, the Raman mirror tilt is controlled to better than 100 nrad bringing the corresponding error 15 on the G measurement to negligible levels. Magnetic fields can perturb atom interferometry measurements both through their mechanical action on the atomic trajectories and via the second-order Zeeman effect. A 103 shielding factor to external magnetic fields is achieved in the interferometer volume by passively isolating the vertical tube through a system of two cylindric µ-metal layers. In addition, systematic shifts that slowly change over a typical time scale of a few seconds and that do not depend on the direction of the effective wave vector k are efficiently cancelled by reversing k with a frequency shift of the Raman lasers at each measurement cycle and taking the average phase shift18. They include the second-order Zeeman effect as well as the first-order light shift. Finally, the main parameters of the experiment have been actively stabilized to optimize the sensitivity and the long-term stability of gravity gradient measurements: the optical intensity of cooling, Raman and probe laser beams are regulated by acting on the RF power driving the acousto-optical modulators; the Raman mirror tilt is controlled by a piezo tip/tilt system. In this way, we obtain a sensitivity to differential accelerations of 3 × 10−9 g for an integration time of 1 s, in agreement with the calculated quantum projection noise limit for 2 × 105 atoms28. Monte Carlo simulation and measurement systematics Gravity gradient measurements are compared to a numerical simulation that evaluates the gravitational potential produced by a given configuration of the source masses, implements the calculation of the phase shift experienced by single atoms at the two interferometers and finally runs a Monte Carlo simulation on the atomic trajectories by varying initial position and velocity out of the density and velocity distribution of the atomic cloud best fitting the experimental density profiles. The gravitational potential generated by the source masses is computed analytically using a multipole expansion. The formula has been ver- 16 ified by calculating the potential of a single cylindrical mass along the classical path of the atoms during the interferometric sequence and comparing it against the numerical evaluation of the exact solution. The multipole expansion has been found to introduce negligible error while drastically reducing the computation time. The phase shift φ for an atom is obtained as 1/h¯ times the classical action, evaluated along the integration paths given by the solution of the Lagrange equation (Eq. 2). We use a perturbative method29 separating the Lagrangian L into L0 = p2/(2m) − mgz, m being the atomic mass, and L1 accounting for the contributions of the Earth’s gravity gradient and of the source masses. The phase shift due to L1 is then calculated by integration over the unperturbed path given by the solution of the Lagrange equation for L0. This approximation gives a negligible error in the G determination since the ratio |L1/L0| computed along the unperturbed paths is lower than 10−7 in our experiment. The validity of the perturbative approach was additionally checked by comparison with the exact solution for the case of a uniform gradient30. The interaction processes between the atom and the Raman laser beams during the interferometry sequence are modelled in the code to account for the finite duration of the Raman pulses. At detection, we evaluate the convolution effects due to both the finite width of the detection beams and the finite bandwidth of the photodiodes collecting the fluorescence light at the F =1 and F =2 channels. Crucial for the Monte Carlo simulation is the spatial and momentum distribution used to describe the atomic sample and its evolution through the velocity selection pulses and the atom interferometer pulse-sequence. Different techniques have been used to characterize the density and velocity profile of the atomic clouds. In particular, time-of-flight measurements provide the required information along the vertical direction, while the profiles in the horizontal plane are measured both by performing CCD 17 imaging and by scanning the atomic cloud using the vertical Raman laser beams with a reduced diameter of 4 mm. The experimental data show that the spatial and velocity profiles of atomic clouds at launch, right after release from the magneto-optical trap, are well described by Gaussian distributions with spatial deviation ∼ 3 mm, velocity deviation ∼ 16 mm/s for the lower cloud, and ∼ 22 mm/s for the upper cloud. The velocity selection during the interferometer reduces the vertical velocity deviation to ∼ 3 mm/s and the horizontal velocity deviation to ∼ 6 mm/s. Average transverse velocities are below 1.5 mm/s. The Monte Carlo simulation is used not only to evaluate the results of our measurements, but also to determine the sensitivity of the final phase angle to the relevant parameters and finally estimate the uncertainty on the G measurement. Derivatives are computed with respect to 28 parameters including the positions of the source masses, positions and velocities of the upper and lower clouds, their time-of-flight at detection, the Gaussian widths of the detection lasers and their position with respect to the atomic clouds, the interferometry sequence duration 2T , the effective wave vector k of the Raman lasers, the gravitational acceleration and its gradient. From an analysis of the results, the spread of the atomic sample along the vertical direction appears to be responsible for the dominant contribution, requiring a measurement of the Gaussian width for both the upper and lower clouds to better than 0.1 mm. The determination of the atomic trajectories requires a good knowledge of the local gravitational acceleration and of the gravity gradient. Gravity acceleration was measured by using an absolute gravimeter to an uncertainty that contributes with a negligible error to the G measurement. A recent measurement of the gravity gradient performed with our atom interferometer gave 3.11 × 10−6 s−2, with an effect of less than 18 2 ppm on the determination of the Newtonian gravitational constant. The G measurement has been corrected for the effect of the finite density of the volume of air that replaces the source masses when they are moved from the configuration C to F and viceversa. The corresponding correction on the G value amounts to (+60 ± 6) ppm. Finally, a correction has been applied to account for the non-uniform density of the tungsten cylinders20; the axial density inhomogeneity of the cylinders gives the main contribution corresponding to a correction on the G value of (+90 ± 20) ppm. The uncertainty budget of our measurement of the Newtonian gravitational constant is detailed in Table 1. Our result is reported in Fig. 3 and compared with the CODATA value3 and recent determinations for G31–43. 30. Antoine, C. & Borde´, C. Exact phase shifts for atom interferometry. Physics Letters A 306, 277–284 (2003). 31. Luther, G. G. & Towler, W. R. Redetermination of the Newtonian gravitational constant G. Phys. Rev. Lett. 48, 121–123 (1982). 32. Karagioz, O. & Izmailov, V. Measurement of the gravitational constant with a torsion balance. Measurement Techniques 39, 979–987 (1996). 33. Bagley, C. H. & Luther, G. G. Preliminary results of a determination of the Newtonian constant of gravitation: A test of the Kuroda hypothesis. Phys. Rev. Lett. 78, 3047–3050 (1997). 34. Gundlach, J. H. & Merkowitz, S. M. Measurement of Newton’s constant using a torsion balance with angular acceleration feedback. Phys. Rev. Lett. 85, 2869–2872 (2000). 19 35. Quinn, T. J., Speake, C. C., Richman, S. J., Davis, R. S. & Picard, A. A new determination of G using two methods. Phys. Rev. Lett. 87, 111101 (2001). 36. Kleinevoss, U. Bestimmung der Newtonschen Gravitationskonstanten G, Ph.D. thesis (University of Wuppertal). (2002). 37. Armstrong, T. R. & Fitzgerald, M. P. New measurements of G using the measurement standards laboratory torsion balance. Phys. Rev. Lett. 91, 201101 (2003). 38. Hu, Z.-K., Guo, J.-Q. & Luo, J. Correction of source mass effects in the HUST-99 measurement of G. Phys. Rev. D 71, 127505 (2005). 39. Schlamminger, S. et al. Measurement of Newton’s gravitational constant. Phys. Rev. D 74, 082001 (2006). 40. Luo, J. et al. Determination of the Newtonian gravitational constant G with time-of-swing method. Phys. Rev. Lett. 102, 240801 (2009). 41. Tu, L.-C. et al. New determination of the gravitational constant G with time-of-swing method. Phys. Rev. D 82, 022001 (2010). 42. Parks, H. V. & Faller, J. E. Simple pendulum determination of the gravitational constant. Phys. Rev. Lett. 105, 110801 (2010). 43. Quinn, T., Parks, H., Speake, C. & Davis, R. Improved determination of G using two methods. Phys. Rev. Lett. 111, 101102 (2013). 20 Figure 1 Schematic of the experiment showing the Rb atom interferometer operated as a gravity gradiometer and the W masses used as source of the gravitational field. For the measurement of G, the position of the source masses is alternated between configurations F and C. Plots of gravity acceleration produced along the symmetry axis by the source masses are also shown; a constant value for Earth’s gravity was subtracted. The spatial regions of the upper and lower atom interferometers are indicated by the thick lines. The vertical acceleration plots show the effect of source masses to cancel the local gravity gradient at the position of the atomic apogees. Figure 2 Experimental data. a) Typical Lissajous figures obtained by plotting the output signal of the upper atom interferometer vs the lower one for the two configurations of the source masses. b) Modulation of the differential phase shift for the two configurations of source masses for a given direction of the Raman beams k vector. Each phase measurement is obtained by fitting a 360-point scan of the atom interference fringes to an ellipse. The error bars, not visible at this scale, are given by the standard error of the least-squares fit to the ellipse. c) Results of the measurements to determine G. Each point is obtained by averaging the signals recorded for each of the two directions of the Raman k vector (see Methods). The data acquisition for each point took about one hour. These data were recorded in different days spanning over one week in July 2013. The error bars are given by the combined error on the angles for the four ellipses. d) Histogram of the data in c). 21 Figure 3 Result of this experiment for the Newtonian gravitational constant G compared with the values obtained in previous experiments and with the recent CODATA adjustments. Only the experiments considered for the current CODATA 2010 value are included and the successive BIPM-13 result. For details on the experiments and their identification with the acronyms used in the figure, see ref. [3] and the additional references in the Methods section. 22 Table 1: Effects, relative corrections and uncertainties considered for the determination of G in this experiment. Uncertainties are quoted as one standard deviation. The third column contains the corrections we applied due to the effects not included in the Monte Carlo simulation. The bias and systematic error from ellipse fitting are evaluated by a numerical simulation on synthetic data. Other effects include: cylinders mass, cylinders vertical position, gravity gradient, gravity acceleration, Raman mirror tilt, Raman k vector, timing. 23 Phase shift (rad) Upper interferometer signal Interferometer phase (rad) 0.8 a) 1.9 0.7 1.8 0.6 1.7 0.5 0.4 0.3 0.3 0.4 0.5 0.6 0.7 0.8 Lower interferometer signal 0.550 1.6 1.5 1.4 0 4 8 12 Time (hours) c) 0.549 0.548 0.547 0.546 0.545 0 20 40 60 80 100 0 10 20 30 Measurement # Counts b) d) 40 6.665 6.670 6.675 G (10-11m3kg-1s-2) NIST-82 torsion balance TR&D-96 torsion balance LANL-97 torsion balance CODATA 1998 UWash-00 torsion balance BIPM-01 torsion balance UWup-02 simple pendulum CODATA 2002 MSL-03 torsion balance HUST-05 torsion balance UZur-06 beam balance CODATA 2006 HUST-09 torsion balance JILA-10 simple pendulum CODATA 2010 BIPM-13 torsion balance THIS WORK atom interferometry 6.680 Effect Uncertainty Air density Apogee time Atomic clouds horizontal size Atomic clouds vertical size Atomic clouds horizontal position Atomic clouds vertical position Atoms launch direction change C/F Cylinders density inhomogeneity Cylinders radial position Ellipse fitting Size of detection region Support platforms mass Translation stages position Other effects Systematic uncertainty Statistical uncertainty Total 10 % 30 µs 0.5 mm 0.1 mm 1 mm 0.1 mm 8 µrad 10−4 10 µm 1 mm 10 g 0.5 mm Correction to G (ppm) 60 91 -13 <2 137 Relative uncertainty ∆G/G (ppm) 6 6 24 56 37 5 36 18 38 4 13 5 6 1 92 116 148