Field-capable gravimeters and gradiometers are currently used to measure gravity and gravity gradients for numerous applications including surveys of underground mineral and natural resources, predictions of earthquakes, and global climate research such as monitoring Earth's icecaps and water tables. Primary use comes from oil and mining companies, the defense industry, large-scale government-funded projects such as the National Geodetic Survey (NGS), academic researchers focused on studies of fundamental Earth's properties, and government agencies such as the United States Geological Survey (USGS).
The current standard for a high performance “absolute” fieldable gravimeter with high accuracy and low drift is a sensor that is based on measuring the position of a mechanical mirror in free-fall. However, this fieldable instrument suffers from a combination of high cost, high power consumption, frequent recalibration, lack of robustness, and long survey times. For example, the falling mirror wears out after repeated use and requires periodic replacement, limiting sensor utility for long-term monitoring. A second competing technology is a superconducting gravimeter. However, this “relative” instrument compares gravity values at different locations and is susceptible to measurement drift. Further the superconducting sensors require cryogenic cooling, limiting their fieldability; cryogenic cooling consumes kilowatts of power and adds substantial size and weight to the sensor.
Various embodiments of the invention are disclosed in the following detailed description and the accompanying drawings.
The invention can be implemented in numerous ways, including as a process; an apparatus; a system; a composition of matter; a computer program product embodied on a computer readable storage medium; and/or a processor, such as a processor configured to execute instructions stored on and/or provided by a memory coupled to the processor. In this specification, these implementations, or any other form that the invention may take, may be referred to as techniques. In general, the order of the steps of disclosed processes may be altered within the scope of the invention. Unless stated otherwise, a component such as a processor or a memory described as being configured to perform a task may be implemented as a general component that is temporarily configured to perform the task at a given time or a specific component that is manufactured to perform the task. As used herein, the term ‘processor’ refers to one or more devices, circuits, and/or processing cores configured to process data, such as computer program instructions.
A detailed description of one or more embodiments of the invention is provided below along with accompanying figures that illustrate the principles of the invention. The invention is described in connection with such embodiments, but the invention is not limited to any embodiment. The scope of the invention is limited only by the claims and the invention encompasses numerous alternatives, modifications and equivalents. Numerous specific details are set forth in the following description in order to provide a thorough understanding of the invention. These details are provided for the purpose of example and the invention may be practiced according to the claims without some or all of these specific details. For the purpose of clarity, technical material that is known in the technical fields related to the invention has not been described in detail so that the invention is not unnecessarily obscured.
A system for gravity measurement is disclosed. The system comprises one or more atom sources to provide three ensembles of atoms. The system comprises two or more laser beams to cool or interrogate the three ensembles of atoms. The system comprises a polarizing beamsplitter and a retro-reflection prism assembly in a racetrack configuration to route the two or more laser beams in opposing directions around a loop topology intersecting the three ensembles of atoms with appropriate polarizations chosen for cooling or interferometer interrogation, wherein the three ensembles of atoms are positioned coaxially when interrogated.
In some embodiments, the cold-atom technology behind the disclosed gravity measurement system resolves and improves on the shortcomings of the current technology. Atoms are ideal test masses that fall in vacuum and do not couple to the external environment. Laser beams used to measure the position of the cold-atom cloud act as extremely precise rulers, accurate to a fraction of the wavelength of light. The disclosed sensor offers significant performance improvement over current state-of-the-art (SOA) fieldable gravimeters at reduced cost, size, and power consumption.
In some embodiments, the disclosed system comprises a high performance, highly compact, low power cold-atom gravity (CAG) sensor to measure minute variations in local gravity. This device is based on novel cold-atom technology that has demonstrated superior performance over SOA sensors. Cold-atom sensor technology is robust and fieldable in noisy locations without expensive site preparation or active isolation systems. The CAG sensor comprises a sensor head, the laser system, and the control electronics. The highly compact sensor head includes a small vacuum cell, beam routing optics, and a vacuum pump. It has no moving parts, providing long term stability while minimizing power consumption. Laser cooled atoms serve as identical test-masses. Atom clouds are continuously regenerated, avoiding the wear and tear issues of competing technologies. The laser system, which is based on low-power semiconductor technology, generates appropriate beams to manipulate the atoms inside the vacuum cell. The laser beams serve as a very accurate ruler that measures the falling atom cloud's position. The compact electronics control the laser system and record and digitize the gravity measurements. The three main subsystems are attached to a common field-deployable platform for simplified operation. Exceptional CAG sensor precision and intrinsic absolute calibration reduces field survey times. The CAG sensor also delivers excellent long-term stability, necessary for applications such as monitoring of underground movement of oil or gas, polar ice sheet melting, or seismic activity.
In some embodiments, the CAG Sensor has no moving parts, reducing its sensitivity to thermal variation and ground noise. With no moving parts to wear out, the sensor can operate continuously for long-term monitoring missions. The disclosed sensor leverages existing technologies to reduce the sensor volume and power consumption. A novel configuration of the sensor substantially suppresses well-known instrument sensitivity to high frequency vibration noise, enabling low cost field deployment. The sensor also incorporates a simplified optical layout to reduce sensor complexity and cost.
In some embodiments, a light pulse atom (LPA) sensor is an atom analog of an optical Mach-Zehnder interferometer. Lasers cool a cloud of atoms, slowing them from a few hundred m/s to a few cm/s without requiring any cryogens. The atoms are launched in vacuum with no coupling to the external environment. A series of interferometer pulses then act as beam splitters (π/2 pulses) and a mirror (it pulse), splitting and recombining the atomic wavepacket. In some embodiments, optical two-photon Raman transitions are used as the interferometer pulses, enabling unique capabilities for precision sensing. For example, electro-optic devices can dynamically shift the optical phase fronts to attain high dynamic range, high sensitivities, and eliminate certain systematic effects. After each interaction with the light pulses, the photon-to atom momentum transfer induces a phase shift φi=ki·xi, where xi is the atom position, and k is the effective wave vector of the transition. The inertial measurement sensitivity comes from the dependence of atomic trajectories on the local acceleration. The net phase difference between two interferometer arms is Δφ=k·αT2, where T is the time between interferometer pulses and α could be gravitational acceleration. The laser thus acts as a high precision ruler that measures the atom's location to a fraction of a wavelength of light. In some embodiments, a 780 nm laser frequency stabilized to <1 kHz, can measure the atom's position to ˜1:1012. Following the interferometer sequence, resonant fluorescence detection measures the excitation probability, Pe, of an atom from a ground to excited state, where Pe=[1−cos(Δφ+φ)]/2, where Δφ=φ1−2φ2+φ3 for the simplest case of a 3-pulse sequence, and φ is an arbitrary phase.
In some embodiments, the phase of one or more interferometer laser pulses is uniformly adjusted to control the phase of the atoms exiting the interferometer, for example, to servo the interferometer phase to a particular value. In some embodiments, the phase of one or more interferometer laser pulses is spatially modulated to create a spatial phase variation across the ensembles of atoms and facilitates quadrature detection and or tracking of the interferometer phase in the presence of platform vibration induced acceleration noise. In some embodiments, detection is done with a camera for spatial resolution. In some embodiments, detection is done using a photodiode.
In some embodiments, an atomic gravity gradiometer includes two or more gravimeters common-mode coupled to substantially reduce sensitivity to various noise sources and systematic drifts, such as vibrations or laser intensity fluctuations that would otherwise decrease measurement sensitivity and increase the noise floor. For example, two or more atomic fountains may share the same interferometer pulses. The inhomogeneous gravity field induces slightly different phase shifts, ΔφUPPER and ΔφLOWER between the upper and lower clouds separated by L. The differential gradiometer signal is Γ=(ΔφUPPER−ΔφLOWER)/(|k|T2L). The Earth's contribution to this signal (˜3,080 E) is static and can be subtracted by calibrating the apparatus prior to a measurement. In some embodiments, the distance from the target mass to the gradiometer is comparable to the atom cloud separation. Note that Γ equals the gravity gradient (the second derivative of the scalar potential) only in the limit where the atom cloud separation is much smaller than the distance to the mass.
In some embodiments, the gradiometer performance is ultimately limited by the signal to noise ratio (SNR). The acceleration sensitivity is Δα=((SNR)keffT2√R)−1, where R is total atom cloud launch rate. Tis limited by the thermal spread of the atomic cloud, and is typically less than a few hundred milliseconds.
In some embodiments, the gradiometer includes the following: (1) The sensor operates in the large momentum transfer Bragg regime, which improves the gradiometer sensitivity by more than an order of magnitude; and (2) an advanced atom cooling trap reduces the initial atom temperature to a few hundred nanokelvin, more than an order of magnitude better than existing 3D magneto-optical traps with polarization gradient cooling.
In some embodiments, in the Bragg scattering regime, there is no change in the internal energy state of the atom. Atoms scatter off the standing wave created by two counter-propagating light fields with the same polarization. The scattering satisfies the Bragg condition: λL sin(θn)=NλdB. Here, λL is the standing wave's wavelength, θn is the scatter angle for Nth order scatter and λdB is the de Broglie wavelength of the atom. In the atomic interferometer, the detuning is varied, ω12, between the two fields to select the desired higher order scattering (N). The detuning equals the atom recoil frequency, ω12=N2k2/(2m), where m is the atomic mass. For the N=40 scattering order in rubidium, ω12=2π×6 MHz. In the gradiometer configuration, the atoms launch collinearly with the light fields. The momentum transfer is along the light wave vector and θn=0.
In some embodiments, the main obstacle to attaining high-order momentum transfer is the nonzero momentum spread of the atoms. As the momentum transfer order increases, the π/2−π−π/2 pulse sequence transfers more of the atoms to other nearby momentum states. This degrades and eventually destroys interferometer contrast. In some embodiments, operation is in the sequential Bragg regime, sequentially applying several π/2 and π pulses. Each pulse imparts only a few k of momentum to the atoms, ensuring a cleaner overall population transfer from the ground to the excited momentum states.
In some embodiments, for large momentum transfer Bragg interferometers using rubidium atoms, the atomic cloud should be cooled to 100s of nanokelvin to reduce its momentum width. These temperatures, below the recoil limit, are readily achieved via three-dimensional delta-kick cooling. The atoms are allowed to expand in the spherically harmonic potential and then switch off the magnetic field after an appropriate time, ton. The cooling process is then less sensitive to potential anharmonicity and the precise timing of the magnetic field turn-off. In some embodiments, the atoms cool to approximately 100 nK scale temperatures in 50 ms, which does not significantly extend the measurement cycle time. The number of atoms in the interferometer sequence remain effectively the same as without the delta-kick cooling stage. Compared to other cooling techniques, such as evaporative cooling, delta-kick cooling is faster and requires lower driving currents. In some embodiments, the delta-kick cooling is implemented using a specially configured time-orbiting potential (TOP) trap.
In some embodiments, resonant fluorescence detection measures the excitation probability Pe that an atom ends up in a particular ground state after the interferometer sequence: Pe=[1−cos(ΔΦ+φL)]/2. ΔΦ=φ1−2φ2+φ3 is the net phase difference between two paths of the interferometer, and φL is the laser arbitrary phase. φi={right arrow over (k)}i·{right arrow over (x)}i is the phase shift the atoms acquire from Raman pulse i when the atoms are at position {right arrow over (x)}i, and {right arrow over (k)}i is the effective wave vector corresponding to the two-photon Raman transition for pulse i. Evaluation of this expression for uniform accelerations or rotations leads to the inertial measurement sensitivity ΔΦ={right arrow over (k)}·{right arrow over (g)}T2, where {right arrow over (k)} is the average Raman wave vector, {right arrow over (g)} is the local gravity, and T is the interferometer interrogation time. Hence, gravity can be expressed as:
In some embodiments, when a sensor comprises two or more interferometers, the measured phase for each interferometer n is ΔΦn. The gravity gradient
between any two interferometers, u and l linearly separated by distance L is then:
Derivative of the gravity gradient between three equally spaced interferometers upper, middle, and lower: u, m, l, where distance from u and m is L, and the distance from m and l is L, and the distance between u and l is 2L is then
In terms of interferometer phase, the derivative of the gravity gradient is then:
In some embodiments, a gradiometer comprising three gravimeters all sharing a common pair of Raman beams has strong noise suppression properties compared to the traditional two-gravimeter design. In fact, the two main obstacles to fielding a two-gravimeter gradiometer is its susceptibility to laser phase noise and to instrument rotations. Both of these obstacles are overcome in a three-gravimeter design. The noise suppression properties to both of these effects are explained below
In some embodiments, the influence of laser phase noise can be written as
where g(t) is the interferometer sensitivity function, φA(t) is the phase of the master Raman laser, and td is the time delay between the two laser frequencies interacting with the atom cloud. The path difference between lasers A and B determines td.
In some embodiments, the frequency difference between the two lasers A, B is assumed to be perfectly stable. For three equally spaced gravimeters comprising a single gradiometer, it can be shown that the following quantity is approximately zero in appropriate limits and therefore is invariant to phase noise:
Where, Δ is the extra time delay at gravimeters II and III as compared to gravimeter I, and v(t)=dφ/dt is the laser frequency. Simplifying, gives the following:
If it is assumed that Δ is small compared to t−td or that v(t) is sufficiently slowly varying on the time scale of 2Δ (for a gradiometer with 1 m baseline separation, this implies that the frequency noise power spectral density rolls off at frequencies above 300 MHz), the first order Taylor expansion is then v1(t−td+Δ)=v1(t−td)+Δ*v1′(t−td), and v1(t−td+2Δ)=v1(t−td)+2Δ*v1′(t−td).
2δΦII−δΦI−δΦIII=∫−∞∞dt g(t)[Δ*v1′(t−td)−Δ*v1′(t−td)]=0
In some embodiments, gradiometer signal sensitivity to platform rotations is examined by analyzing analytically and numerically the signal for a gravimeter, a gradiometer comprised of two vertically offset gravimeters, and a gradiometer comprised of three vertically offset gravimeters. The signal from the three instruments is examined in the presence of platform rotations due to (1) nearby mass and (2) earth's gravity field. The analytical error model for the three devices is derived by solving the equation of motion for an atom {right arrow over (r)}(t) launched in the rotating frame of reference fixed on the Earth's surface:
Here, {right arrow over (R)} is the distance from the center of the earth to the location of atom launch, {right arrow over (Ω)} is the Earth's rotation, {right arrow over (g)} is the constant gravity acceleration on Earth's surface, and Tij is the gravity gradient tensor. The coordinate system has x, y, z facing North-East-Down respectively, fixed on the Earth's surface. Equation (1) can be solved analytically if the gravity gradient term is treated as a perturbation. Ignoring the gravity gradient term, Eq. (1) is then solved with appropriate initial conditions, and the interferometer phase in the absence of gravity gradients is calculated as
φ0={right arrow over (k)}·[{right arrow over (r)}(t0)−2{right arrow over (r)}(t0+T)+{right arrow over (r)}(t0+2T)],
where T is the delay between successive interferometer pulses and t0 is the time of the first Raman pulse. The contribution of the gravity gradient to the total interferometer phase is then:
Here, Tzz is the zz component of the gravity gradient, and vrec is the atom recoil velocity. The perturbative approach is accurate at the μrad phase level. The total interferometer phase is then φTotal=φ0+φgg. For a three cell gravimeter, the signal is φGRAD=2φM−φU−φL, where φM is the middle gravimeter.
In some embodiments, an analysis of the contribution of platform rotations to the measurement phase noise of two cell and three cell gradiometers using both analytic and numerical approaches shows that the effects of rotations are small for both gradiometer configurations for reasonably small rotation angles during the measurement interval (e.g., 1 mrad). For higher magnitude rotational noise (e.g., 10 mrad), the phase error for a two cell gradiometer approaches 0.1 mrad, or 10% of a typical gradiometer noise floor. The rotational noise is completely suppressed for a 3 cell gradiometer in the earth's gravity field. In the presence of a nearby mass, the noise contribution of platform rotation is small for both three and two cell gradiometer configurations.
In some embodiments, atoms are launched without sub-Doppler cooling.
Although the foregoing embodiments have been described in some detail for purposes of clarity of understanding, the invention is not limited to the details provided. There are many alternative ways of implementing the invention. The disclosed embodiments are illustrative and not restrictive.
This application claims priority to U.S. Provisional Patent Application No. 62/095,671 entitled GRADIOMETER CONFIGURATION INVARIANT TO LASER PHASE NOISE AND SENSOR ROTATIONS filed Dec. 22, 2014 which is incorporated herein by reference for all purposes.
This invention was made with Government support under DE-AC52-07NA27344 awarded by the United States Department of Energy. The Government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
62095671 | Dec 2014 | US |