Methods of scanning an object using a pendulum gravimeter are disclosed. The pendulum gravimeter may include an interferometer, such as a sagnac interferometer, to determine a displacement on the pendulum by way of a mirror attached to the pendulum. Scanning of the object may be performed in 1D, 2D, or 3D, and may result in an image of the object. In another aspect, a mass may be tracked while in motion using a pendulum gravimeter to detect the gravitational attraction of the object.
|
12. A method of scanning an object, comprising:
determining, using a sagnac interferometer, a first displacement of a pendulum caused by the object;
moving the object with respect to the pendulum;
measuring a second displacement of the pendulum caused by the object; and
determining a profile of the object based on the displacement of the pendulum.
1. A method of scanning an object, comprising:
providing a pendulum gravimeter comprising:
a pendulum having a mirror, the pendulum configured to have a restoring force for restoring the pendulum to a neutral position; and
an interferometer configured to measure a displacement of the pendulum from the neutral position;
measuring the gravitational field of the object at a first measurement location;
measuring the gravitational field of the object at a second measurement location; and
calculating the mass distribution of the object based on the measured gravitational fields.
2. The method of
measuring a displacement of the pendulum from the neutral position; and
determining the gravitational field based on the measured displacement.
3. The method of
4. The method of
5. The method of
9. The method of
10. The method of
11. The method of
|
This application claims priority to U.S. Provisional Application No. 62/609,314, filed on Dec. 21, 2017, now pending, the disclosure of which is incorporated herein by reference.
This invention was made with government support under contract no. W911NF-15-0296 awarded by the ARMY Research Office. The government has certain rights in the invention.
The present disclosure relates to measuring gravitational attraction.
There is great interest in measuring gravitational fields nearly perpendicular to the Earth's gravitational field. A sensitive gravimetric camera could “see” a person walking around through hundreds of meters of concrete or monitor shipping crates for hazardous materials. Currently, precision gravimetry tends to have relatively low resolution and is used primarily for large scale structures. Gravimetry is used extensively in mapping the earth's local gravity, oil and gas exploration, mining, mapping temporal geological shifts, the determination of Newton's gravitational constant and gravitationally imaging opaque systems. Precision of the order of 1 mGal (1 Gal=1 cm/s2) to 1 μGal are often used for mapping geological variations. Both relative and absolute measurements are employed. A standard in the industry for absolute gravimetry is measuring interference fringes due to the free-fall of a corner cube in one arm of a Mach-Zehnder interferometer with a sensitivity of 100 μGal/Hz1/2. Another competing gravimetric technology employs atomic interferometry achieving a resolution of 100 nGal after two days of integration. At the forefront, superconducting spheres suspended in the field of a superconducting coil have achieved 3 nGal resolution after one month integration or 1 nGal after one year. The current standard has been set by Kasevich's group (Stanford), which has achieved 1 nGal per shot (2 seconds) in their recently built 10 m atomic fountain and is working on a squeezed atom implementation with the potential for 20 dB improvement per shot.
Embodiments of precision gravimeters are disclosed herein. In a first embodiment, a gravimeter is based on a simple pendulum with a few 10's of ms response time. A shot noise limited acceleration sensitivity of 100 fm/s2/Hz1/2 or 10 pGal/Hz1/2 can be achieved with a few mW of low-coherence laser power. The gravimeter is a physical pendulum which is allowed to oscillate in only a single plane. The pendulum and interferometric readout measure relative gravitational fields perpendicular to the Earth's gravitational field using the fringes in the dark port of a displaced Sagnac interferometer operated in the inverse weak value regime. The system is small, lightweight, simple to build, easily portable, inexpensive, and requires no cryostat, or cold atom technology. This gravimeter can be scaled to several detectors to build a gravity camera. In this way, a one kg mass moving around at a few meters can be remotely tracked using only the gravitational attraction of the mass.
In a second embodiment, an gravimeter uses a long-period torsion pendulum. The shot noise limited predicted sensitivities are a function of the restoring torque or alternatively considered as the restoring period. A device according to this embodiment is able to remotely image precision gravimetric density profiles in nontrivial mass distributions (e.g., tumors in-situ). With relatively reasonable device parameters, the shot noise limited sensitivity of the torsion balance embodiment is sufficient to measure the center of mass of a 1 gram mass 10 cm away with picometer resolution. Used in a near-field configuration, the device may be used to, for example, measure the curvature of space due to the presence of 1 Watt of energy
at a distance of a few mm.
For a fuller understanding of the nature and objects of the disclosure, reference should be made to the following detailed description taken in conjunction with the accompanying drawings, in which:
and fixing the center of mass in place, able to rotate only. A nearby mass M creates a torque on this torsion pendulum, causing it to rotate to an equilibrium angle δθ0, which is detected optically. Lengths, angles, and mass labels are shown in the figure.
,x)−g(−
,x), which is plotted here versus x. This is also called the “point spread function” in Fourier imaging techniques. Note it is an odd function, and most sensitive at x=0.
Precision gravimeters are disclosed herein, suitable for studying gravitational fields perpendicular to the earth's gravitational field. The gravimeters may be used for fundamental studies of gravity, high resolution bulk imaging, as a part of a gravity camera, etc. Shot noise limited measurements with 1 mW of laser power predict precision of the order of 10 pGal/Hz1/2.
With reference to
The inferometer 30 includes a laser 32 for providing a laser beam A. A beam splitter 34 of the inferometer 30 is configured to split the laser beam A into a first beam B and a second beam C. Where a Sagnac inferometer is used, the first beam B and second beam C propagate in opposite directions and are each directed around a looped path to impinge on a mirror 22 affixed to the pendulum 20. The first beam B impinges on the mirror 22 at a first location to provide a first reflected beam D, and the second beam C impinges on the mirror 22 at a second location to provide a second reflected beam E. For example, where the pendulum is configured to move in a vertical plane, such as a simple pendulum, the first location and second location are arranged vertically with respect to each other. In embodiments where the pendulum is configured to move in a horizontal plane, such as a torsion pendulum, the first location and second location are arranged horizontally with respect to each other.
A photodetector 36 is configured to receive the first reflected beam D and the second reflected beam E. A controller 40 is in electrical communication with the photodetector 36. The controller 40 is configured to determine a displacement of the pendulum 20 based on interference between the first reflected beam D and the second reflected beam E. The controller 40 may be, for example, a programmable controller, programmed to determined the displacement of the pendulum.
In some embodiments, the pendulum is a torsion pendulum comprising a mass assembly. In some embodiments, the mass assembly is a rigid rod. In some embodiments, the mass assembly is made up of at least two masses rigidly coupled to each other, such as, for example, coupled using a rigid rod. The mass assembly is balanced on a fulcrum and configured to rotate about the fulcrum. The mass is biased to a first rotational position and configured such that a gravitational attraction (e.g., from an object) will rotate the mass assembly away from the first rotational position. For example, the mass assembly may be suspended by a wire, which acts as the fulcrum at the location where the wire is affixed to the mass assembly. In this way, the mass assembly is balanced while suspended by the wire, and able to rotate by twisting the suspension wire. Such twisting (rotational) motion also creates a restorative torque in the wire.
I. Weak Values
Precision gravimetry using weak value amplification techniques is discussed herein. Weak values were born through asking fundamental questions about quantum measurement limits. Unlike expectation values, weak values consider a normalized expectation of an operator (e.g., the Pauli operator Â=|++|−|−
−|) using pre- and post-selected quantum states ψi,f
Because weak values can be much larger than their respective expectation values when ψf|ψi
→0, they have been used to amplify small effects. Weak value amplification has been shown to be exceptionally valuable in suppressing technical noise in precision measurements. While these techniques do not beat the shot noise limit (with some exceptions), they can come close to reaching it because of the dramatically suppressed technical noise. Of particular interest is the recent inverse weak value work where an angular tilt measurement noise floor of 200 frad Hz−1/2 was achieved. Remarkably, this sensitivity was for signals down to 1 Hz, where noise suppression can be incredibly difficult. This tilt corresponds to a displacement of less than a hair's breadth at the distance of the moon in one second of measurement time using only a few milliwatts of laser power. Using these techniques for precision gravimetry, even at the classical optical fundamental limits, they would push gravimetric sensitivity by several orders of magnitude beyond the state-of-the-art.
II. Inverse Weak Value Interferometry
A specialized weak values interferometer employs a laser beam of transverse width σ in a Sagnac interferometer. The laser beam enters a beamsplitter and propagates in opposite directions around a Sagnac interferometer. When the beams recombine, a small relative phase ϕ between the two returning beams causes a small transverse tilt, k, of the mirror attached to the pendulum to be amplified in the (nearly) dark port of the beamsplitter. The weak value limit occurs when ϕ»kσ. For this particular setup, the amplification of the small transverse tilt k shows up in the dark port beam as a spatial shift by kσ2/ϕ. In terms of weak values, the pre-selected state ψi is the field after passing through beamsplitter the first time. The post-selected state ψf is set by a combination of a phase shift and the second pass through the beamsplitter. The tilt of the mirror yields a small amount of which-path information of a photon in the interferometer, which is a weak measurement of σz.
Conversely, for an inverse weak values experiment, the parameters satisfy the inequality ϕ«kσ. In this case, we fix transverse tilt k and use the known k to amplify a small unknown phase. In this latter experimental regime, the interference pattern of the two beams in the dark port is now a bimodal distribution with a dark fringe at the center of the interference pattern for ϕ=0. The dark fringe moves rapidly with small changes in relative phase. These phase shifts are determined by measuring the relative intensity of the left versus right side of the interference pattern via a split detector. The amplification of the phase in this inverse weak value regime is given by the mean shift of the beam in the dark port, ϕ/k, which now is proportional to the inverse weak value Aw−1.
A displaced Sagnac interferometer may be used to measure the relative phase shift ϕ for this inverse weak value regime. In a displaced Sagnac interferometer, two beams propagate with a transverse displacement (albeit parallel) in opposite directions. A small tilt of a mirror inside the interferometer causes a relative phase between the two paths since the path length increases for one path and decreases for the other. At the output port of the interferometer, the beams are brought back together to interfere with each other.
The shot noise limited angular resolution can be understood from a geometric argument. The relative phase between the two paths goes as δϕ=2√{square root over (2)}πLθ/λ, where √{square root over (2)} comes from impinging at 45 degrees, L is the distance between the centers of the beams propagating in opposite directions (see
we find.
where N is the number of detected photons. Using L=1 cm, approximately 3 mW of laser power, and a wavelength of 500 nm, we achieve a shot noise limited angular sensitivity of 30 frad/Hz1/2. The inverse weak value method of readout for the optical phase ϕ can achieve this shot-noise limited sensitivity, up to a factor of √{square root over (π)}/2 associated with the resolution loss on the split detector.
This same inverse weak value setup with a displaced Sagnac is further described below. Further consideration is given to the physical limitations and sensing capabilities when the tilt mirror in the simple pendulum is replaced with a mirror rigidly connected to a torsion pendulum as shown in
Simple Pendulum Gravimeter
An example of a pendulum gravimeter according to an embodiment of the present disclosure is shown in
A displaced Sagnac interferometer operated in the inverse weak value regime may be used to read out the tilt of a critically damped physical pendulum. A laser beam traverses two trajectories in the interferometer. One trajectory strikes the mirror on the pendulum at one height and the other beam strikes the mirror at a different height. The pendulum has a mass m and the height difference between the beams on the mirror is L. The pendulum design is such that it can only swing in a single plane. For convenience in discussion the operation of this embodiment, the pendulum was considered to be acted upon only by Earth's gravity and by the force of a distant remote object of mass M , which induces a force F perpendicular to Earth's gravity on the test mass of the pendulum. In the presence of the external object, the pendulum mass has an angle θ with respect to the vertical. The equilibrium condition is then F cos θ=mg sin θ. The small angle approximation was used to arrive at
Rather than measure the angular displacement θ of the mass, the relative distance the mass shifts Lθ was measured. For normal reflection from the pendulum, the relative phase is given by ϕ=2kLθ, where k is the wavenumber. A force equation
was obtained. For shot noise limited measurements, an approximation
is made, where N is the number of photons measured For relatively modest parameters, this system can measure very sensitive relative changes in gravity. For example, using visible light photons (k≈107 m−1), approximately 1016 photons/second flux (approximately 1 mW of cw laser power) in the interferometer, a relative pendulum height L of 2 cm, an interferometric sensitivity on the order of 10−13 m/s2 or 10 pGal after one second of integration is made. Note that the presently-disclosed pendulum gravimeter differs from several previous pendulum gravimeters, which were based on a Fabry-Perot etalons.
A system according to the present embodiment can be used to measure modest masses at relatively large distances. The force F is simply the gravitational force due to a mass M at a distance d. Solving for d, yields:
Using the same numbers as above with a mass of M=1 kg, yields d≈25 m. In other words, the device is sensitive enough to measure one kilogram mass up to 25 meters away with a shot-noise limited signal-to-noise ratio of 1 after one second of integration.
Achieving the shot noise limited sensitivity for gravimeters is exceedingly difficult owing to significant environmental noise that can couple into the gravimeter. The present disclosure draws upon ideas in weak values where technical noise can be suppressed, in some cases, by several orders of magnitude and phase stabilization can be made to a high degree of precision. A displaced Sagnac interferometer is operated in the inverse weak value domain. For the inverse weak value, ktσ<1, where kt is the transverse momentum shift and σ is the beam width. The operation of the interferometer in this configuration results in an average beam displacement given by ϕ/kt, which can make the displacement a very sensitive function of phase for small kt. However, this large shift does not change the ultimate sensitivity, because it reduces the signal strength at the detector by (ktσ cos(ϕ/2))2.
The displaced Sagnac interferometer allows the beams to sample two different heights of the reflecting mirror while still achieving Sagnac stability and using low coherence sources. The parity of the interferometer is determined by the number of mirrors, in addition to the beam splitter, inside the interferometer. The Sagnac configuration cancels translational vibration noise since both the upper and lower beams are translated in the same way, which preserves the phase stability. Further, noise sources outside the interferometer such as beam jitter or detector jitter are not amplified. The Sagnac technique was found to be advantageous over the Michelson technique, having improved stability by several orders of magnitude.
Unfortunately, weak values cannot suppress all types of noise. If the noise is in the device that is generating the signal, there is no attenuation of the noise. This means that environmental noise (e.g., acoustic, microseismic, laser heating, laser shot-noise torques) may couple into the pendulum and limit the resolution of the detector. Measures may be used to lower such noise and improve the performance of the device. For example, noise may be reduced by using active isolation, having the interferometer in high vacuum, increasing the effective mass of the device against acoustic noise, balancing the readout laser with respect to the oscillation axis, etc.
In the experimental embodiment shown in
To readout the pendulum shifts, the interferometer was aligned very close to the perfect dark port until a bimodal interference pattern with a large dark fringe results. One axis of the piezo actuated mount was then used to place the dark fringe in the middle of the intensity pattern. The lobes of the bimodal intensity pattern were aligned to two halves of a quad detector so that their difference signal could be measured. The rapidly moving fringe then resulted in large changes in the amount of light striking each side of the quad detector. The sensitivity of the device was calibrated using two different techniques. The first technique used readouts of the sum and difference voltages on the quad cell along with the measurement of the beam waist to determine kt and σ, respectively. The second technique was based on the careful characterization of the other piezo-actuated mirror in the interferometer. By applying a calibrated voltage on the other mirror and measuring the balanced signal, the mass-equivalent shift could be calibrated. Then the sensitivity of the device was calibrated with a known mass at a given distance while monitoring the signal voltage. This allowed for calibration of the signal response for the given test mass for that alignment. While the ultimate sensitivity of the inverse weak value and the standard homodyne technique resulted in the same shot noise limit, the inverse weak value was useful, because the small phase shifts that result from very small pendulum tilts can lead to large dark fringe shifts, which are easily determined from the detector.
Preliminary measurements of the exemplary device were performed using a lock-in amplifier. A 1.5 kg mass was placed on a 3 rpm turntable. The center of the turntable was 46 cm from the center of mass of the pendulum. The 1.5 kg test mass was placed 17 cm from the center of the turntable. The experiments were performed early in the morning when the acoustic and seismic noise was at a minimum. The lock-in amplifier was able to find the signal and hold the voltage to 1 part in 10 of the central voltage. The lock-in voltage was close to the simulated value for the given experimental parameters. This would imply that gravitational attraction was measured to a few nanoGal after a few 10's of minutes integration. Further tests are expected to demonstrate linearity in mass and the quadratic dependence on distance. The use of further damping, vacuum, and vibration isolation is expected approach the shot-noise limited sensitivity without the need for the lock-in amplifier.
Torsion Pendulum
In a second embodiment, shown schematically in
The torsion pendulum will allow ultra-precise near-field gravimetric measurements and far-field gravity gradiometry for high resolution imaging.
III. Geometry and Torque Calculation
Suppose there is a system comprising a point mass M that is detected via a torsion pendulum comprising two masses m, connected via a rigid massless rod of length 2. The mass M is located according to
A. Converting Torque Into Angle
In order to detect this small torque, recall that the pendula have a linear restoring torque quantified by the torsion spring constant κ, such that:
τext=−κδθ, (4)
where δθ denotes the angular distance from its equilibrium position. Note that κ can be empirically found by finding the period of the oscillations. When the pendulum swings freely,
τext=Iα=2m2{umlaut over (θ)}=−κδθ, (5)
where I=2 m2 is the moment of inertia. Putting this equation in the form {umlaut over (θ)}+ω02θ=0, the natural frequency of the pendulum is found to be ω0=(1/
)√{square root over (κ/(2m))}, which may be inverted to find κ from a measurement of the frequency, or period of the pendulum
,
Adding in damping of the pendulum brings the dynamics into the form of Eq. (7).
For a vertical simple pendulum subject to Earth's gravitational field, the restoring force is simply the gravitational acceleration from the earth, causing a restoring torque of τ=−gmδθ for small angles. This also gives rise to dynamics of the form Eq. (7) but with a natural frequency given by ω0=√{square root over (g/)}.
The source of the gravitational signal is a mass M near the pendulum, thereby applying an external torque τ. We wish to measure this signal. We assume that τ can be time-dependent in general. Damping of the oscillations will be critical for a quickly responding detector, so we also add a velocity-dependent damping term, to find the equation of motion of a damped/driven oscillator,
{umlaut over (θ)}+2ζω0{dot over (θ)}+ω02θ=τ(t)/I (7)
Here, ω0=√{square root over (κ/I)}, and ζ is a dimensionless damping coefficient. The underdamped case corresponds to 0≤ζ≤1, whereas the overdamped case corresponds to ζ>1. The general solution for τ=0 (the homogeneous solution) can be expressed as
where θ0 and {dot over (θ)}0 are initial conditions for the pendulum, and the terms inside the braces { . . . } describe decay for the overdamped case, and oscillations in the underdamped case.
If τ is fixed in time, for large damping ζ, then the oscillator will converge to its new equilibrium position exponentially in time with a rate ζω0 according to the solution (8). After this time, the angular displacement can be approximated by the fixed point of (7), given by
We wish to design the pendulum to respond sensitively to stimuli from the target objects, but do not want it to oscillate for a long time before returning to a new equilibrium position. There is a trade-off between sensitivity of the measurement and the speed of the response as will be explored in the following sections.
B. Pendulum Model
The pendulum is fixed with respect to its center of mass motion, and is allowed to only rotate about its center of mass in the plane of the figure. This geometry can be analyzed by computing the torque about the middle of the torsion pendulum.
Plane trigonometry dictates that the distances defined in the figure are given by:
r12=(x−)2+d2, r22=(x+
)2+d2, r02=x2+d2. (10)
The gravitational force between mass M and mass mj, according to Newton, is given by:
The torque τ generated on mass j is given by:
τj=Fj cos θj, (12)
where Fj is the magnitude of gravitational force of mass M on mass mj.
For a Simple Pendulum (SP), we only have one of the masses in the pendulum of
The above results are relevant for a single test mass other than the pendulum mass. In the following sections we will use either the earth's gravity as the restoring force (as usual pendulums do), or by orienting the pendulum perpendicular to the earth's field, can also use the restoring force of a rod to obtain longer periods.
For a balanced Torsion Pendulum (TP), we include torques that nearly counterbalance each other (the torque on mass 1 is positive in sign, and the torque on mass 2 is negative in sign). The net torque is given by:
where we have replaced cos θj=d/rj.
The simple pendulum responds to the bare force on the sensing mass, and thus decays as 1/r2 with respect to the test mass distance. The torsion pendulum balances the average force, and thus responds to the gradient of the field across the size of the torsion pendulum. This effect leads to a less sensitive response to objects far away; it may be beneficial since it efficiently screens out far away objects and allows the sensor to focus on nearby objects.
C. Limiting Case
In some experiments, this expression (Eq. (14)) can be further simplified because it is expected that «d, x, r0. The expression for τΣ is proportional to the difference of the functions g(
)−g(−
), where:
Since is a small parameter, one can approximate g(
)−g(−
)≈g′(0)(2
). It is found that:
so that, to a good approximation,
where we approximate θ1≈θ2=θ, and write x=r0 sin θ and d=r0 cos θ. In this limit, the sensor does not respond to the net force, but rather to its gradient, as indicated by the r−3 law.
In this limit, the one-armed device (SP) equilibrates to an angle
while its two-armed counterpart equilibrates to
In both cases we have used Eq. (9). We stress that in both cases, the sensing mass m only appears in the natural frequency, and in the case of the torsion pendulum, the length also drops out, indicating that small sensors work as well as large ones so long as their periods are the same.
These expressions can be applied to make an approximate survey of the sensitivity of the device to different objects. Specifically, we show the best-case angular response to a target M at distance r0 in
TABLE I
Example parameter values
Pendulum mass
m
100
g
Wavelength of light
λ
500
nm
Pendulum length
5
cm
Period of oscillator
500
s
Torsion spring constant
κ
7.9 × 10−8 kg m2/s2
Length between beams on the mirror
L
1
cm
The above analysis may be extended to a continuous mass distribution by replacing the mass M by a differential element dM=ρ(x)dx, where a body with mass per unit distance ρ(x) distributed along the x direction is imagined. In that case, the next torque for such a mass distribution is given by:
τΣ=∫ƒ(x)ρ(x)dx. (20)
In the general case, ƒ(x)=dGm(r1(x)−3−r2(x)−3), whereas in the limiting case, it is given by ƒ=6
2dGmx/(d2+x2)5/2.
IV. Noise Considerations
It is noted that the pendulum will experience several kinds of noise that may be mitigated in order to reach the fundamental limits of angle detection that the system is capable of. This section focuses on three types of noise: thermal noise, measurement heating noise, and quantum noise of the oscillator.
A. Thermal Noise
Contributions of the thermal noise from the surrounding environment can be computed via the equipartition theorem assuming large temperatures. Both the mean kinetic energy and potential energy are given by the thermal energy for one degree of freedom each. In general,
In the limit of high temperatures, the equipartition of potential energy indicates that
1/2κδθ2
=1/2kBT, (22)
which gives the typical rms noise of the torsion pendulum,
δθrms=√{square root over (kBT/τ)}. (23)
The value can be estimated using the values given in Table I, and room temperature kBT=4.1×10−21 J to find, δθrms=2.3×10−7 rad. In order to access below the picoradian regime, it will therefore be necessary to either cool the oscillator, or to time-average the signal for some time. One could also increase the value of κ, but that would also decrease the angular precision.
B. Measurement Heating
As demonstrated in the first section, it is advantageous to apply as much optical power as possible to the interferometer to maximize the precision of the angle measurement. However, because the gravitational sensor is freely moving, it is possible that the sensing laser may drive excitations, effectively heating the torsion pendulum. We will now calculate the effect of this heat, which may put a bound on the sensing power.
The displaced Sagnac geometry (see
from the axis of rotation. We can imagine that photons each carrying energy hc/λ=k0c arrive at the mirror at discrete times, either at
On average, since the intensity of the light on the left and right side of the mirror is the same from the 50/50 beamsplitter, we have
The time-averaged optical power being applied to the pendulum is N
k0c/T, which could be understood as arising from
provided that the duration T is much longer than the typical time between individual photons, which arrive at the tn, and where k0/√{square root over (2)}, where the factor √{square root over (2)} arises because the photons strike the mirror at an angle. The time-averaged torque may be expressed using the same logic as above, i.e.,
or for T»|tn−tn−1|,
Since the number of photons striking on either side of the point of rotation is equal on average, we see that =0; however, there will be fluctuations from the coherent states of light, which depend on
(N+−N−)2
=
N+2−N−2−2N+N−
=
N
. (27)
The term N+N−
=0 because the ± photon beams are uncorrelated. Furthermore, given the geometry and using statistical properties of coherent states,
Therefore the torque fluctuations averaged over a given time interval T are given by
where we have re-written the duration of the experiment in terms of the power δ
=
/κ2, or taking the square root, fluctuations in the coherent states of light lead to
Equating with the shot noise limited value (2), we find that the measurement heating makes a smaller contribution to the noise than shot noise, when the condition
is satisfied. This implies that for the values in Table I, the noise contribution due to measurement heating would overtake the shot noise when the measurement laser power exceeds
C. uantum Noise
An intriguing aspect of the oscillators is the fundamental limitation of sensitivity due to quantum noise. As will be shown, the shot noise limited resolution is approximately equal to the quantum ground state uncertainty of the oscillator when the integration time is approximately equal to the period of the pendulum. From a quantum mechanical perspective, ground state quantum noise limitations are quite interesting in light of the large masses used in these experiment. Such studies may be valuable in probing quantum gravity. On the other hand, this also places fundamental noise limits on the resolution.
To determine the ground state angular uncertainty, we set the mean potential energy of the oscillator to the ground state energy of the oscillator
(1/2)κ(δθ)2
=(1/4)
ω, (31)
which follows from Eq. (21) when kbT«ω. Solving for the angle we obtain
δθrms=√{square root over (ω/2κ.)} (32)
Using values listed in Table I, we find δθrms=2.9×10−15 rad. This resolution can be achieved when the integration time is roughly equal to the period of the pendulum using a few milliwatts of laser power.
V. Limits of Resolution
The results of the previous sections can now be combined to give the sensitivity limits of the SP and TP to forces, which can be translated into either mass or range uncertainty. Using Eq. (9) and the angular uncertainty Eq. (2), we find at the optimally sensitive response point (x=0 so d≈r, θ=0), the resolution of (usual) simple pendulum acceleration a, relative to the gravitational acceleration g of the SP to be
Consequently, the acceleration uncertainty in units of the accelerations due to gravity near the surface of the earth is simply the same as the angular uncertainty.
If instead, we consider a one-armed torsion pendulum with a torsion constant of κ, oriented perpendicular to the gravitational field of the earth, then the period of the oscillation can be much longer. The angular resolution is given by Eq. (18) so the acceleration uncertainty is reduced to
For the parameters in Table I, this reduces the acceleration uncertainty by a factor of 1.6×10−6, leading to 60 zepto g Hz−1/2. Remarkably, the speed is only a thousand times slower, because of the inverse square relationship of Eq. (6). In either geometry, the acceleration is given by a=GM/r02, so the sensitivity of the acceleration to a change in test mass δM at fixed r0, or a change in the distance r0 for fixed test mass M is given by
δa=GδM/r02−2GMδr0/r03, (35)
from which the mass or distance uncertainty is easily found. The response of the one-armed torsion pendulum is plotted in
For a balanced torsion pendulum, a test mass far from the pendulum will respond according to Eqs. (16, 18). Setting θ=π/4 for maximum sensitivity, the angular response to a gravitating body will be
The r0−3 law gives a smaller sensitivity, but also screens off distant objects. This cannot be directly translated into acceleration of a single mass, but gives the response of the detector to the gradient of the gravitational field. The torsion pendulum response is plotted in
We now briefly discuss the angular response of both types of pendula as the test mass is placed at different angles relative to the axis of rotation. The one-armed torsion pendulum has blind spots at θ≈0 and θ≈π, where a target mass applies no torque, and its sensitivity is maximized at θ≈π/2 and θ≈3π/2. The two-armed torsion pendulum has four blind-spots, as illustrated in , and that the smaller we make
, the larger the deflection angle will get (the moment of inertia in the denominator wins out over the greater torque with greater arm length). The
dependence cancels out entirely from the two-armed device, except for its appearance in the natural frequency.
Near Field Measurements
To illustrate the sensitivity of the device, assume that the mass of interest is in the near-field of the pendulum (i.e., a distance of the mass from one pendulum mass is much closer than the separation between the two pendulum masses), the minimum detectable acceleration of a single mass close to the test mass is:
where l is the radius of the pendulum masses, k is the wavenumber, L is the distance between the beams reading out the deflection of the pendulum, N is the number of photons measured during an integrated time interval and T is the period of the pendulum. For a critically damped pendulum, the damping time is of the order of T/2. For r=0.2 m, L=0.1 m, T=1000 s and a 3 mW beam, the shot noise limit for acceleration is 10−22 m/s2. To demonstrate the force sensitivity of the device, such a gravimeter is sufficient to measure the effective mass (m=E/c2) of 1 Watt (m=1017 kg) of energy at a few mm away from the pendulum, assuming all other effects due to the energy could be ameliorated (for example, one can imagine a resistor dissipating 1 W close to the pendulum with an ultra-broadband mirror reflecting the radiation away from the gravimeter in between). Using a pendulum mass of 100 g, the pendulum could measure forces with a resolution of 10 yocto Newtons, which could be used in almost any force measurement from protein folding to the Watt Balance. Interestingly, a ns, visible, single-photon normally reflecting from a mirror, imparts 10−19 N. There may be many limitations on the resolution of the device including the fact that shot noise in the read-out laser will cause torques on the mirrors.
GravCam, Solving the Inverse Problem
More than one of the devices can be arrayed. For example, with four such devices, a three-dimensional determination of the position and mass of an object can be determined assuming a large enough signal-to-noise ratio using the correlation of the independent detectors. Movement of the object may also be determined. Because the damping times are still relatively long, the objects must move relatively slowly, but with a knowledge of the velocity of an object (velocity can be determined from temporal responses of the pendula), additional structure information can be gathered. Further, with a large array of sensors allowed to pivot about their tilt-axis or translated, even static environments can be gravitationally imaged.
The sensitivity of embodiments of the disclosed device is expected to be sufficient to gravitationally map density variations within the body. As a simple proof-of-principle simulation, a uniform smaller sphere of known density, but unknown size and position inside a larger sphere of known radius and uniform density was considered. For simplicity, the position along a single axis was estimated. With two unknowns, size and single-axis position, only two gravity measurements were needed. Two gravity detectors separated by a distance of one meter with the same z=0 cm, y=0 cm coordinate were used. The center of each sphere was assumed to have a z=30 cm and y=0 m, but unknown x for the smaller sphere, albeit within the larger sphere. For the simulation, the density of the larger sphere was assumed to be ρL=1000 kg/m3 (the density of water) and only a slightly higher density ρS=1100 kg/m3 for the smaller sphere. The radius of the larger sphere was 10 cm and the center of the large sphere was half way between the two detectors. A precision of 10 pGal was assumed and using a smaller sphere radius of 3 cm. Solving the inverse problem, the simulation was found to predict the position of the smaller sphere to within 0.5 cm along the x axis and the simulation determined the mass to within 1%. The density problem is easier in Airport scanners where hidden objects typically have much higher density contrasts. In some embodiments, the detectors are moved such that even static environments can be imaged. Several pendula may be used in a gravity camera along with the torsion pendulum discussed later.
In another example, a torsion pendulum may be used to obtain high-resolution images, for example, when using the gravitational effects of both masses of the pendulum (see
Setting this equal to the restoring torque and taking the differential, an advantageous spatial resolution is found to be:
Using the numbers from above, along with a mass of 1 g, a 2 picometer spatial resolution is achievable at 20 cm away. This would indicate that the gravimetric gradiometer may be used for performing high resolution imaging well beyond the diffraction limit deep inside a bulk object.
Suppose the torsion pendulum is used as a sensor to detect an unknown mass density profile ρ(x). Suppose the mass distribution is moving at a constant speed v along the x direction. This will give a time-dependent torque τΣ(t) (proportional to the equilibrium angle δθ0(t)). Then, given this time dependent signal, the mass distribution ρ(x) can be found.
The time-dependent torque is given by:
τΣ(t)=∫−∞∞ƒ(x)ρ(x−vt)dx. (39)
Defining a coordinate y=vt, the central insight is to interpret this integral as a convolution (ƒ*ρ)(y) of two functions ƒ and ρ (where * denotes convolution). By the convolution theorem, the unknown function ρ(x) is given by:
where {tilde over (ƒ)}(k) and {tilde over (τ)}Σ(k) are defined as the Fourier transforms of the functions ƒ(x), which is a known function, and τΣ(y), which is the measured data. It is noted that both the Fourier transform of the measurement data and the inverse Fourier transform to find the solution can be made with numerical fast Fourier transforms. Consequently, this procedure can be done very quickly, once the functional form of Fourier transform of ƒ(x) is known. The analytic form of the function ƒ is calculated below, giving the final step in the solution of the problem.
Fourier Transform of ƒ(x)
The Fourier transform of ƒ (in the general case) is proportional to the difference of the Fourier transform of g(x,)−g(x,−
), where:
The integral is nontrivial because of the branch cut. The result is found as:
where K1(z) is the first modified Bessel function. This function vanishes at k=0 since K1(0) is finite, and decays exponentially at large k. Taking the difference of {tilde over (g)}(k,) at
and −
, as well as restoring the prefactors, the following result is found,
{tilde over (ƒ)}(k)=(−4i)Gm sin(k)|k|
K1(|k|d). (43)
Together with the main result Eq. (40), the problem is solved. It is noted that the velocity v is unimportant, so long as it is slow enough to allow the torsion pendulum to come to equilibrium, and the optics to readout the angle to sufficiently high precision.
An example of a 1-D mass distribution is shown in
Generalization to Two and Three Dimensions
Suppose that a 2D or 3D object is to be gravitationally image. There will be a mass distribution ρ(x,y) or ρ(x,y,z). One may proceed by assuming the torsion pendulum can rotate only about one axis (the z-axis), and consequently only the z-component of the torque is computed. This is found as:
τz={circumflex over (z)}·({right arrow over (r)}i×{right arrow over (F)}i), (44)
where {right arrow over (r)}i is a vector from the origin to mass mi, and {right arrow over (F)}i is the force vector acting on mass i from the differential mass region at point (x,y,z). The net z-component of the torque is found as,
where r′1,2=√{square root over ((x∓)2+y2+z2)} has been modified. Here, the basic effect of the z direction is seen as contained in the mass distribution itself, and the additional fall off in the denominator the further the mass is from the pendulum. In the y direction, the constant value d now can change as y, as one gets closer and further away from the pendulum. The two dimensional case corresponds to fixing the y variable.
The continuum limit can be reached by introducing a mass density ρ, with a net torque of:
In 2D, one only integrates over x and z. The inversion in 2D and 3D works in principle as the 1D case. Focusing here on the 2D case of imaging in a plane, the torque is only mapped out as the mass density is being raster scanned. This gives the torque as a function of the relative coordinates τz,Σ(x′,z′). As in the 1D case, this is the same as Eq. (46), with ρ replaced by ρ(x−x′,z−z′). This is now a two-dimensional convolution. Consequently, ρ(x,z) can be represented with the inverse Fourier transform ,
ρ(x,z)=[{tilde over (τ)}z,Σ({right arrow over (k)})/{tilde over (ƒ)}({right arrow over (k)})], (48)
where the known function {tilde over (ƒ)}=Gm({tilde over (g)}(
)−{tilde over (g)}(−
)) is given by the multi-dimensional Fourier transform,
Finite Sampling in 3D
In the previous sections, a continuous data curve was assumed as the imaged object was displaced in time. It is likely more realistic to consider a set of data points taken at discrete locations around the object in question. One may consider finite data sets from such an array of measurements at discrete points. For simplicity, consider a uniform grid of object points in three dimensions, {xi,yi,zi}, where i=1, . . . , p. These points approximate a three dimensional object by shrinking the mass within a lattice unit cell into a point of mass mi. This mass array is sensed with a detector sampling the gravitational field at a series of discrete points {xjd,yjd,zjd}, with j=1, . . . , q. These points can be taken around the outside of the body that is being imaged, with several detectors, or a single detector that is displaced (or kept stationary, with the body being imaged moved around). Since the gravitational force is proportional to the two attractive masses, and is additive as a vector, a general relation can be derived between the measured quantities θj and the system masses mi, given as a linear relationship,
θ=R·M, (50)
where θ is a vector of dimension q with elements θj, and M is a vector of dimension p with elements mi. The matrix R is a q×p matrix, which encodes both the geometric relationships between detector and system, as well as the type of detector in use (a free particle, a simple pendulum, a torsion pendulum, etc.). For each of those detectors, the measured quantities may be different (acceleration, displacement, an angle, etc.)
If the dimensions of the system and detector degrees of freedom match, then the matrix has a unique inverse (assuming is has a nonzero determinant). Otherwise, the solution is either over or under determined in general. The Moore-Penrose pseudoinverse, R+ can be used to give the best fit solution to this linear equation,
Mbƒ=R+·θ, (51)
where the pseduo-inverse solution minimizes the norm of all solutions M. The pseduo-inverse may be constructed in a simple way with the singular value decomposition, R=UΣV†, where U,V are unitary matricies, and Σ is a diagonal matrix,
R+=VΣ+U†. (52)
The pseduo-inverse of the diagonal matrix Σ, denoted as Σ+, may be found by inverting the nonzero elements of Σ and taking the transpose.
Matrix R is now computed for the special case of a torsion pendulum detector discussed above. Adapting the geometry shown in
Assuming the masses of the torsion pendulum are the same at every point, and therefore set the masses of the pendulum mj=m for all j. The net displacement angle of the pendulum located at position i, for a restoring torsion constant κ, is given by
where the last equality defines the matrix Rji.
In order to switch to another imaging plane, the torque given in Eq. (53) can be modified accordingly. For example, to align the pendulum in the z direction, and sample in the plane of constant y, the numerator would involve yjd−yi, and the denominator's dependence would be in the y degree of freedom. The matrix R′ji obtained from that plane would then be appended to the matrix Rji from the first measurement plane, to obtain a generally rectangular matrix of all the data.
Once the response matrix R is defined, it is then straightforward to change the mass distribution within the region, and to use the same reconstruction matrix, which will still be valid so long as the geometric relationships between the measurement planes and grid points are the same.
Numerical Examples
The problem to solve is a deconvolution problem: Given a known “point spread function,” ƒ({right arrow over (x)}), and a measured torque function versus scan coordinate {right arrow over (x)}′, the mass density profile ρ({right arrow over (x)}) is to be reconstructed where there is added white Gaussian noise ν from the measurement process,
τ=ƒ*ρ+ν (55)
Here * indicates the convolution operation. The simple inverse of the Fourier transform (in the absence of noise) is problematic, since the Fourier transform of the point spread function rapidly goes to zero, so its inverse diverges rapidly, leading to instability in the deconvolution. There exist a suite of deconvolution techniques, the classic one being the Wiener filter, defined as:
Here is the signal-to-noise ratio of the data collection, which cuts off the divergence in inverting the small values of the transform of ƒ. The reconstructed approximation is given by:
92 rc=τ*W (57)
Note that because [ƒ](k=0), the overall scale of ρ cannot be directly reconstructed, so the inverse is defined only up to an overall factor; that is, Cδ(k) can be added to the Fourier transform, or a constant function to ƒ(x).
One having skill in the art will understand that the Wiener filter is the oldest example in what is now a lengthy number of possible methods to help in the deconvolution problem. The present disclosure should not be limited to the disclosed methods—e.g., the non-limiting, exemplary use of the Wiener filter.
To obtain a 2D scan, the device may be raster scanned or multiple devices may be used. In embodiments where a 3D scan is obtained, multiple detectors can be incorporated to work together.
Resolution of the Imaging System
With respect to the sensitivity of such a detector to fine-grain details of an image, in an optical context, this resolution is limited by the Rayleigh Criterion. However, in the gravitational case, there is no known “wavelength,” so in principle there is no limit to the resolution within Newtonian gravity.
In practice, however, there will be limit set by how precisely the detector can measure the angle shift of the torsion pendulum. Consider Δθ as the standard deviation of the measured angle. As can be seen from
The resolution Δx may be found by error propagation,
where one considers when the slope of θ is steepest with respect to x—i.e., x=0.
The equilibrium angle for a point mass M is given by θ(x,,y)=C(g(x,y,
)−g(x,y,31
)), where C=
yGmM/κ. The slope of the function at the origin S0 is computed,
In order to made the uncertainty on Δx as small as possible, the slope S0 is to be as large as possible. Therefore, S0 versus is maximized, keeping y fixed. The maximal slope is found to be when
=√{square root over (2/3y)}≈0.8615y, or y≈1.22
, which corresponds to
where again, y is assumed to be fixed.
Estimating the position uncertainty δx, m and M are assumed to be about 1 kg, and the precision with which the object can be located is considered. Using the same torsion constant as before, κ of 10−3 Nm/rad, let y=1 m and G≈6×10−11 Nm2/(kg)2. For the angular uncertainty, a value found in previous work was used, Δθ=400 frad=4×10−13 rad. Putting this all together, for a 1 kg mass at 1 meter,
Δx≈6 μm, (61)
which shows very good results. If the test mass is reduced to a mass of 1 mg, the resolution decreases by three orders of magnitude to about 6 mm, which is still good.
Torsion Pendulum Rod Embodiment
Considering another geometry where there is just a rod pendulum of mass m and length L. A point mass is considered, located at perpendicular position y from the rod, and lateral position x from the center of mass of the rod.
Integrating the net torque produced by a mass element dm over the length of the rod, the torque is found as,
The moment of inertia is I=(1/12)mL2. It is noted that this moment of inertia of the bar is down by a factor of ⅓ compared to the dumbell (two masses interconnected by a rigid rod) of the same length and total mass (m→2 m, L→2). This gives the period T since τΣ=iα is of the pendulum. Here, the torsion constant κ is expressed, in terms of T, as:
This gives the angle θ=τΣ/κ.
The function Eq. (62), divided by κ, is the new “point spread function” for imaging a mass distribution with the rod imaging system.
Exemplary Methods
In an embodiment, a method of scanning an object includes determining a first displacement of a pendulum caused by the object. The displacement may be determined using a Sagnac inferometer as described herein. The object is moved with respect to the pendulum. For example, either of the object or the pendulum can be moved, or both can be moved. A second displacement of the pendulum (i.e., displacement of the pendulum caused by the object) is measured. A profile of the object is determined. The profile may include, for example, a position (such as, for example, a 1D position), a mass, and/or a mass distribution (e.g., gravimetric density profile). As described above, additional measurements of the pendulum displacement may be made while the object is in various locations.
In another embodiment, a method of scanning an object comprises providing a pendulum gravimeter. The pendulum gravimeter may be of any variety described herein. For example, in an embodiment, the pendulum gravimeter has a pendulum with a mirror. The pendulum is configured to have a restoring force, for restoring the pendulum to a neutral position. The pendulum gravimeter may further include an interferometer, such as, for example, a Sagnac interferometer. The interferometer is configured to measure a displacement of the pendulum from the neutral position. The method includes measuring the gravitational field of the object at a first measurement location and measuring the gravitational field of the object at a second measurement location. Measuring the gravitational field may include measuring a displacement of the pendulum from the neutral position, and determining the gravitational field based on the measured displacement. The mass distribution of the object is calculated based on the measured gravitational fields. The mass distribution of the object may be calculated in any way described herein.
For example, the object may be measured at a plurality of additional measurement locations for a total of q measurement locations. The object may be scanned in three dimensions and the measured gravitational field of the object θ is represented by a matrix of the measurements θj at each measurement location j having coordinates {xjd,yjd,zjd}, where j=1, . . . , q and q is the total number of data points. The mass distribution of the object M may be determined as a matrix of point masses mi over a set of grid points {xi,yi,zi}, where i=1, . . . , p and p is the total number of grid points. The mass distribution M of the object may be determined according to M=R+·θ, where R+ is the pseudoinverse of a predetermined scanner matrix R.
The calculated mass distribution may be used to create, for example, an image of the object. In such a case, scanning the object may be considered as imaging the object, and the gravimeter may be considered to be a “gravimetric camera.”
In some embodiments, the object is measured at the second measurement location using a second pendulum gravimeter. In some embodiments, the object is measured at the second measurement location by moving the pendulum gravimeter (e.g., moving the gravimeter and subsequently obtaining the measurement). In some embodiments, the object is measured at the second measurement location by moving the object with respect to the pendulum gravimeter (e.g., moving the object and obtaining the reading when the gravimeter is at the second measurement location with respect to the object).
Although the present disclosure has been described with respect to one or more particular embodiments, it will be understood that other embodiments of the present disclosure may be made without departing from the spirit and scope of the present disclosure.
Jordan, Andrew N., Howell, John C.
Patent | Priority | Assignee | Title |
Patent | Priority | Assignee | Title |
10591632, | Sep 16 2015 | University of Rochester | Interferometric pendulum gravimeter, and method for same |
2137963, | |||
3011346, | |||
3424006, | |||
3429184, | |||
3672473, | |||
3688584, | |||
3744909, | |||
3926054, | |||
3960344, | Jul 29 1974 | Precision mobile gravity survey and guidance system | |
4333341, | Jun 09 1980 | Determination of the acceleration of gravity | |
4383762, | Feb 14 1980 | ERWIN KAYSER - THREDE GESELLSCHAFT MIT BESCHRANKTER HAFTUNG | Two-beam interferometer for Fourier spectroscopy with rigid pendulum |
4841772, | Dec 03 1987 | MARYLAND, UNIVERSITY OF | Three-axis superconducting gravity gradiometer |
5144150, | Mar 07 1990 | Matsushita Electric Industrial Co., Ltd. | Configuration measuring apparatus |
5203199, | Oct 12 1990 | TELEDYNE INDUSTRIES, INC A CORP OF CALIFORNIA | Controlled acceleration platform |
5204568, | Sep 23 1991 | GWR Instruments | Superconducting bearing for borehole and survey gravimeters |
5351122, | May 22 1992 | MICRO-G SOLUTIONS, INC | Absolute gravity inline measuring apparatus incorporating improved operating features |
5892151, | May 21 1997 | Lockheed Martin Corporation | Differential interferometric ballistic gravity measurement apparatus and method |
6079267, | Nov 06 1998 | The United States of America as represented by the United States | Method and apparatus for measuring gravitational acceleration utilizing a high temperature superconducting bearing |
6494091, | Feb 28 2001 | Apparatus for measuring a gravitational attraction of the earth | |
6615660, | Mar 07 2001 | Lockheed Martin Corporation | Real-time self-compensating gravity gradiometer instrument |
6671057, | Oct 31 2001 | Schlumberger Technology Corporation | Gravity and differential gravity sensor, and system and method for monitoring reservoirs using same |
7401514, | Jun 08 2005 | MAX-PLANCK-GESELLSCHAFT ZUR FORDERUNG DER WISSENCHAFTER E V ; MAX-PLANCK-GESELLSCHAFT ZUR FORDERUNG DER WISSENSCHAFTEN E V | Asymmetric torsion balance gravimeter |
7466396, | Oct 13 2005 | Nikon Corporation | Lithography apparatus and method utilizing pendulum interferometer system |
8141259, | Jun 09 2010 | Schlumberger Technology Corporation | Method of determining the dip of a formation |
8442777, | Oct 15 2012 | KING FAHD UNIVERSITY OF PETROLEUM AND MINERALS | System and method for measuring rolling resistance |
8661894, | Mar 29 2007 | ENI S P A | Microgravimeter for geophysical prospecting |
9291742, | Feb 17 2014 | Micro-g LaCoste, Inc. | Double pendulum gravimeter and method of measuring gravity using the same |
9507048, | Feb 17 2014 | Micro-g LaCoste, Inc. | Pendulum arm flexure and method of construction |
9939551, | Sep 24 2012 | Schlumberger Technology Corporation | Systems, devices and methods for borehole gravimetry |
9952154, | Jun 22 2016 | The Charles Stark Draper Laboratory, Inc | Separated parallel beam generation for atom interferometry |
20020152810, | |||
20030136190, | |||
20110303005, | |||
20150234088, | |||
20170371065, |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Dec 21 2018 | University of Rochester | (assignment on the face of the patent) | / | |||
Nov 25 2019 | JORDAN, ANDREW N | University of Rochester | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 051432 | /0602 | |
Dec 03 2019 | HOWELL, JOHN C | University of Rochester | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 051432 | /0602 |
Date | Maintenance Fee Events |
Dec 21 2018 | BIG: Entity status set to Undiscounted (note the period is included in the code). |
Feb 04 2019 | SMAL: Entity status set to Small. |
Date | Maintenance Schedule |
Apr 20 2024 | 4 years fee payment window open |
Oct 20 2024 | 6 months grace period start (w surcharge) |
Apr 20 2025 | patent expiry (for year 4) |
Apr 20 2027 | 2 years to revive unintentionally abandoned end. (for year 4) |
Apr 20 2028 | 8 years fee payment window open |
Oct 20 2028 | 6 months grace period start (w surcharge) |
Apr 20 2029 | patent expiry (for year 8) |
Apr 20 2031 | 2 years to revive unintentionally abandoned end. (for year 8) |
Apr 20 2032 | 12 years fee payment window open |
Oct 20 2032 | 6 months grace period start (w surcharge) |
Apr 20 2033 | patent expiry (for year 12) |
Apr 20 2035 | 2 years to revive unintentionally abandoned end. (for year 12) |