Geodesic Interferometer: Concepts & Applications

• Geodesic interferometers are systems that transduce geometric quantities into measurable phase shifts, delays, or displacements across optical, matter-wave, radio, and orbital implementations.
• Various designs like DFMI optical sensors, Michelson quantum-geometry probes, and GNSS–VLBI setups employ tailored calibration chains to overcome technical noise and achieve sub-picometer precision.
• These instruments support applications from satellite geodesy and gravity gradient measurements to tests of emergent quantum geometry and formation-flying dynamics in space.

Searching arXiv for the cited papers and closely related work on geodesic interferometry. A geodesic interferometer is an interferometric system in which the measured phase, delay, or displacement is tied to a geodesic quantity: the separation of nearly free test masses in geodesy, the correlated path-length response of a Michelson interferometer to hypothesized quantum-geometrical fluctuations, the shortest geodesic used to close a non-cyclic trajectory on the Bloch sphere, or the geometric delay associated with a terrestrial or orbital baseline. The arXiv literature uses the term across optical, matter-wave, radio, and formation-flying contexts rather than for a single canonical instrument (Isleif et al., 2019, Hogan, 2012, Zhou et al., 2019, Skeens et al., 2023, Ito, 2023).

1. Terminological scope and operational meanings

In the cited literature, “geodesic interferometer” has several precise meanings. In compact deep frequency modulation interferometry (DFMI), the device is an optical displacement and tilt sensor proposed for optical gradiometers for satellite geodesy and as a dimensional sensor for ground-based gravity experiments. In Hogan’s Michelson-style formulation, the instrument probes whether emergent quantum geometry produces spatially coherent, transverse position indeterminacy between world lines. In the SU(2) matter-wave experiment, the relevant geodesic is the shortest geodesic on the Bloch sphere connecting the endpoints of a non-cyclic evolution. In GNSS–VLBI work, the interferometer measures geometric delay on a terrestrial baseline for geodetic applications including local tie measurements. In formation-flying studies, the geodesic content is embedded in the orbital dynamics of separated spacecraft whose relative motion must be held within a small-disturbance environment (Isleif et al., 2019, Hogan, 2012, Zhou et al., 2019, Skeens et al., 2023, Ito, 2023).

Context	Geodesic quantity	Primary observable
DFMI optical gradiometer	Relative displacement on a measurement baseline	Displacement and tilt
Michelson quantum-geometry probe	Transverse position indeterminacy between world lines	Correlated displacement noise
SU(2) matter-wave interferometer	Shortest geodesic on the Bloch sphere	Non-cyclic geometric phase
GNSS–VLBI geodetic interferometer	Geometric delay on baseline $\mathbf{B}$	Group delay, phase delay, delay rate
Formation-flying orbital interferometer	Relative orbit around a chief spacecraft	Control acceleration and $\Delta V$ budget

This range of usage suggests that the unifying feature is not a single hardware topology but the encoding of geometry into an interferometric observable.

2. Compact optical implementations for geodesy and gravity gradients

A compact optical realization is the DFMI sensor built as a single-component, prism-shaped fused-silica optic with an approximately $25\,\mathrm{mm}$ base and less than $1\,\mathrm{in}^3$ volume. A polarization-maintaining fiber collimator injects a single Gaussian beam into the prism at an angle of incidence of about $4.1^\circ$. The optic generates a reference arm that propagates entirely inside the glass and a measurement arm that exits the prism, reflects off a nearby test mass, and re-enters the prism. The two beams recombine at the second internal beam-splitter face and leave via complementary “direct” and “transmitted” ports that are phase-shifted by $\pi$. A single photodiode per port, or a quadrant photodiode for tilt, converts the interference signal into a photocurrent, which is amplified by a low-noise transimpedance stage, digitized at at least $250\,\mathrm{kHz}$ per channel, and processed by a software phasemeter (Isleif et al., 2019).

The signal chain uses deep frequency modulation of an external-cavity diode laser with carrier frequency $f_0 \simeq 281\,\mathrm{THz}$ at $\lambda = 1064\,\mathrm{nm}$, sinusoidal frequency modulation at $f_m \simeq 0.8\,\mathrm{kHz}$, and deviation $\Delta V$0. The modulated light is split into a reference interferometer with arm imbalance $\Delta V$1 for laser-noise stabilization and multiple DFMI sensors with $\Delta V$2 per prism. Readout is performed by a real-time Levenberg–Marquardt least-squares fit of a Bessel-series expansion of the first $\Delta V$3 harmonics, with $\Delta V$4. The extracted fit parameters are the interferometric phase $\Delta V$5, contrast $\Delta V$6, modulation index $\Delta V$7, and FM phase $\Delta V$8.

The time-domain photocurrent model is

$\Delta V$9

with demodulation obtained by multiplying $25\,\mathrm{mm}$0 by $25\,\mathrm{mm}$1 and $25\,\mathrm{mm}$2, low-pass filtering, and fitting the complex amplitudes $25\,\mathrm{mm}$3. The interferometric phase for a pair of fields $25\,\mathrm{mm}$4 is written as

$25\,\mathrm{mm}$5

The phase-noise and displacement-noise analysis is explicit. The displacement noise is decomposed as

$25\,\mathrm{mm}$6

Typical contributors reported in the experiment include laser frequency noise after stabilization of about $25\,\mathrm{mm}$7 at $25\,\mathrm{mm}$8, giving $25\,\mathrm{mm}$9; shot noise at $1\,\mathrm{in}^3$0 optical power of about $1\,\mathrm{in}^3$1; thermal path-length fluctuations of about $1\,\mathrm{in}^3$2 at $1\,\mathrm{in}^3$3 with $1\,\mathrm{in}^3$4 character; parasitic beam-jitter coupling with coupling factor about $1\,\mathrm{in}^3$5 and projected contribution about $1\,\mathrm{in}^3$6; and electronic digitization noise and non-linearities of about $1\,\mathrm{in}^3$7 above $1\,\mathrm{in}^3$8. Measured performance reaches displacement noise $1\,\mathrm{in}^3$9 for $4.1^\circ$0, with a floor of about $4.1^\circ$1 at $4.1^\circ$2, and tilt noise $4.1^\circ$3 for $4.1^\circ$4.

For geodesy, the gradient resolution is related to displacement noise by

$4.1^\circ$5

so that for $4.1^\circ$6 and $4.1^\circ$7, one obtains $4.1^\circ$8, where $4.1^\circ$9. Because $\pi$0 is baseline-independent for the prism geometry, the gravity-gradient noise scales as $\pi$1. The same source gives expected values of $\pi$2 for satellite geodesy with $\pi$3 and about $\pi$4 for ground-based gradiometers with $\pi$5 (Isleif et al., 2019).

The satellite-integration problem is correspondingly stringent. Alignment tolerances require prism wedge perpendicularity below $\pi$6 and beam-incidence alignment within $\pi$7 to retain more than $\pi$8 heterodyne contrast and suppress tilt-to-length coupling. Thermal control requires sub-mK stability of the prism mount and bench with coefficient of thermal expansion around $\pi$9 to keep $250\,\mathrm{kHz}$0. The proposed implementation uses a single ECDL, PM-fiber distribution, at least $250\,\mathrm{kHz}$1 per sensor, flight-qualified ADC/DAC with at least 18-bit dynamic range, a radiation-hard FPGA/CPU for real-time phasemetry, fused-silica prisms on a Clearceram bench, vibration isolation, and multi-layer thermal shielding. Radiation tolerance, outgassing, FM-deviation drift, phasemeter calibration drift, and redundancy remain explicit flight-qualification challenges.

3. Michelson geodesic interferometers and emergent quantum geometry

A distinct usage arises in proposals to test whether classical geometry is only an approximate macroscopic behavior of a quantum system at the Planck scale. In this framework, the mean 4-position of a macroscopic body is represented by operators $250\,\mathrm{kHz}$2 whose noncommutativity is postulated as

$250\,\mathrm{kHz}$3

where $250\,\mathrm{kHz}$4 is a fundamental length scale, $250\,\mathrm{kHz}$5 is the expectation position four-vector, $250\,\mathrm{kHz}$6 is the dimensionless 4-velocity, and $250\,\mathrm{kHz}$7 is the Levi-Civita tensor. In the rest frame and at equal time this reduces to

$250\,\mathrm{kHz}$8

The resulting uncertainty relation implies $250\,\mathrm{kHz}$9, so transverse position indeterminacy grows as the square root of macroscopic separation rather than remaining at the Planck length (Hogan, 2012).

The interferometric translation of this hypothesis is formulated for a Michelson geometry with orthogonal arms of length $f_0 \simeq 281\,\mathrm{THz}$0. A length fluctuation $f_0 \simeq 281\,\mathrm{THz}$1 in arm $f_0 \simeq 281\,\mathrm{THz}$2 produces the optical phase shift

$f_0 \simeq 281\,\mathrm{THz}$3

and the dark-port response depends on the differential phase $f_0 \simeq 281\,\mathrm{THz}$4. The beamsplitter is then modeled as undergoing a random walk in the plane transverse to the two arms, with coherence time $f_0 \simeq 281\,\mathrm{THz}$5, so the differential arm-length picks up a correlated jitter of order $f_0 \simeq 281\,\mathrm{THz}$6.

The predicted one-sided displacement power spectral density for a simple Michelson interferometer is

$f_0 \simeq 281\,\mathrm{THz}$7

For $f_0 \simeq 281\,\mathrm{THz}$8, the spectrum scales as $f_0 \simeq 281\,\mathrm{THz}$9; for $\lambda = 1064\,\mathrm{nm}$0, it tends to a constant of order $\lambda = 1064\,\mathrm{nm}$1. The corresponding phase-noise spectrum is

$\lambda = 1064\,\mathrm{nm}$2

Because photon-shot noise, thermal noise, and seismic noise are uncorrelated between independent interferometers, the proposed discriminator is cross-correlation between co-located Michelsons. If two nominally identical interferometers are placed with beamsplitters separated by much less than $\lambda = 1064\,\mathrm{nm}$3, their dark-port signals should share the same emergent geometry fluctuations for lags up to about $\lambda = 1064\,\mathrm{nm}$4,

$\lambda = 1064\,\mathrm{nm}$5

The Fermilab Planck-precision implementation uses two $\lambda = 1064\,\mathrm{nm}$6-arm Michelsons separated by about $\lambda = 1064\,\mathrm{nm}$7, a $\lambda = 1064\,\mathrm{nm}$8–$\lambda = 1064\,\mathrm{nm}$9 continuous laser at $f_m \simeq 0.8\,\mathrm{kHz}$0, and a bandwidth of $f_m \simeq 0.8\,\mathrm{kHz}$1 to $f_m \simeq 0.8\,\mathrm{kHz}$2 to cover the first few fringes of $f_m \simeq 0.8\,\mathrm{kHz}$3. The critical frequency is about $f_m \simeq 0.8\,\mathrm{kHz}$4 and the first zero is near $f_m \simeq 0.8\,\mathrm{kHz}$5. The target displacement sensitivity is better than $f_m \simeq 0.8\,\mathrm{kHz}$6 and the integration time is of order $f_m \simeq 0.8\,\mathrm{kHz}$7 per full-band cross-spectrum. Commissioning is described as complete and initial runs ongoing; no statistically significant correlated signal had yet been reported, and current upper limits already constrained certain variants of the noncommutative models at the tens-of-percent level (Hogan, 2012).

The significance of this program lies in its separation from ordinary gravitational-wave phenomenology: the proposed fluctuations are described as not metric fluctuations and not gravitational waves, but a different class of macroscopic deviations from classicality.

4. Matter-wave geodesic-rule interferometry and gravitational red-shift

In a third meaning, the “geodesic” in geodesic interferometer is the shortest geodesic on the Bloch sphere required to close a non-cyclic quantum trajectory. For a normalized state $f_m \simeq 0.8\,\mathrm{kHz}$8 with $f_m \simeq 0.8\,\mathrm{kHz}$9, the Mukunda–Simon definition gives

$\Delta V$00

$\Delta V$01

and the gauge-invariant geometric phase

$\Delta V$02

For a two-level system mapped to $\Delta V$03, if the evolution follows a non-cyclic path $\Delta V$04 from $\Delta V$05 to $\Delta V$06, then

$\Delta V$07

where $\Delta V$08 is the Berry–Pancharatnam one-form and $\Delta V$09 is the shortest geodesic joining $\Delta V$10. Equivalently, the phase is one half of the signed area enclosed by $\Delta V$11 and $\Delta V$12 (Zhou et al., 2019).

The experimental realization is a spatial SU(2) matter-wave interferometer using a Bose–Einstein condensate of $\Delta V$13 prepared in $\Delta V$14. On-chip RF and magnetic-gradient pulses create two spatially separated wavepackets, and a third RF pulse rotates each local spin by angle $\Delta V$15 about the $\Delta V$16-axis so that both packets occupy the same latitude on the Bloch sphere. Because the packets are separated by a few tens of microns in $\Delta V$17, a subsequent magnetic-gradient pulse imprints a differential phase $\Delta V$18 between the $\Delta V$19 and $\Delta V$20 components, placing the packets at $\Delta V$21 and $\Delta V$22.

After release from the trap, the wavepackets overlap in time of flight, and the total phase is extracted by fitting the density modulation

$\Delta V$23

with fringe period $\Delta V$24. The reported result is an unambiguous confirmation of the geodesic rule. For $\Delta V$25, when $\Delta V$26 is scanned from $\Delta V$27 to $\Delta V$28, the total phase remains rigid in each hemisphere and then exhibits a sudden $\Delta V$29 jump as $\Delta V$30 crosses $\Delta V$31. Fitting $\Delta V$32 gives $\Delta V$33 with precision around $\Delta V$34 and baseline phase $\Delta V$35 to about $\Delta V$36. Subtracting the dynamical phase yields a geometric phase that changes sign at the equator and jumps by $\Delta V$37 when $\Delta V$38, in exact agreement with the half-area prediction.

The same experiment connects the result to the Pancharatnam phase by choosing the north-pole state $\Delta V$39 as a third vertex. It further proposes an application to gravitational red-shift. If the two packets carry identical internal clock states with splitting $\Delta V$40 and sit at heights $\Delta V$41 and $\Delta V$42, then to first order the proper-time difference is $\Delta V$43, leading to

$\Delta V$44

For $\Delta V$45, this becomes

$\Delta V$46

The paper states that with atomic masses around $\Delta V$47, vertical separations $\Delta V$48–$\Delta V$49, and interrogation time $\Delta V$50, one can achieve $\Delta V$51–$\Delta V$52, within current interferometric precision of about $\Delta V$53. The main noise sources are magnetic-gradient instability, vibration-induced phase noise, and atom-number fluctuations (Zhou et al., 2019).

5. Geodetic baseline interferometers with GNSS and VLBI

A geodetic interferometer in the radio domain is realized by pairing a commercial geodetic-quality GNSS antenna with a VLBI radio telescope. The GNSS element uses a modified Topcon CR-G5 antenna whose built-in RF passband board is removed and replaced by two high-pass filters with cutoff around $\Delta V$54, two low-noise amplifiers, a bias-tee, and a second band-pass filter covering roughly $\Delta V$55–$\Delta V$56. The amplified broadband signal enters a High-Rate Tracking Receiver with a $\Delta V$57 ADC, $\Delta V$58-bit in-phase and $\Delta V$59-bit quadrature sampling, digital downconversion, CIC decimation, FIR filtering, $\Delta V$60-bit I/Q quantization, and a polyphase channelizer with up to nine $\Delta V$61 bands. Outputs comprise real-time GNSS observables converted to $\Delta V$62 RINEX and raw baseband I/Q stored in HDF5 and converted to VDIF. The VLBI side at Fort Davis uses the $\Delta V$63 VLBA dish, an L-band cryogenic receiver with $\Delta V$64 and $\Delta V$65, two single-polarization $\Delta V$66 IF bands, a $\Delta V$67 2-bit ADC, a Mark 5 recorder, and a DiFX software correlator (Skeens et al., 2023).

The geometric delay for a source at effectively infinite distance is

$\Delta V$68

where $\Delta V$69. For a GNSS satellite at finite position $\Delta V$70,

$\Delta V$71

Including instrumental and clock terms, the total model delay is

$\Delta V$72

The cross-correlation function is

$\Delta V$73

and fringe fitting searches for maxima of a channelized, delay-rate-corrected coherent sum. From the phase, group delay, and delay rate, one obtains the baseline vector through weighted least squares,

$\Delta V$74

A central technical element is the Precise Point Positioning extension method. For carrier phase $\Delta V$75, a Kalman filter estimates receiver position, clock bias, wet troposphere, and ambiguities, with typical residual clock-bias uncertainty of $\Delta V$76. The clock correction is converted to a phase correction

$\Delta V$77

The reported effect is substantial: with rubidium clocks, raw coherent integration saturated at about $\Delta V$78 because of decorrelation, whereas after PPP-based phase correction no signal-to-noise loss was seen up to at least $\Delta V$79.

The experimental results include a strong interferometric response with signal-to-noise ratio over $\Delta V$80 from GPS and Galileo satellites, and detections of natural radio sources including Galactic supernova remnants and active galactic nuclei as far as one gigaparsec. On the $\Delta V$81 baseline, typical $\Delta V$82 accumulations on GPS gave $\Delta V$83; Galileo BOC signals produced about $\Delta V$84. On $\Delta V$85 baselines, the SNR dropped by a factor of about $\Delta V$86–$\Delta V$87. The residual delay uncertainties per scan are given as approximately $\Delta V$88 from thermal noise at SNR about $\Delta V$89, less than $\Delta V$90 from residual tropospheric mismatch on short baselines, about $\Delta V$91 from post-PPP clock circulation residual, up to about $\Delta V$92 from GNSS-antenna multipath, and less than $\Delta V$93 from digitization and quantization. The resulting net group-delay precision is at most about $\Delta V$94, corresponding to baseline components of at most about $\Delta V$95. With extended broadband and multi-scan stacking, the same work anticipates sub-millimeter repeatability in a full local-tie campaign (Skeens et al., 2023).

This instrument therefore occupies a different branch of the geodesic-interferometer family: its principal observable is not test-mass displacement or geometric phase on $\Delta V$96, but group and phase delay referenced directly to terrestrial reference-frame realization.

6. Formation-flying geodesic interferometers in geocentric orbit

For spaceborne interferometry, the central problem is relative orbital dynamics rather than internal phase extraction alone. In a near-circular geocentric orbit, the chief–deputy relative motion in the Earth-centered inertial frame is modeled by

$\Delta V$97

and after transformation to the LVLH frame and linearization for $\Delta V$98, by the Clohessy–Wiltshire equations

$\Delta V$99

with

$25\,\mathrm{mm}$00

A bounded relative orbit is the general circular orbit

$25\,\mathrm{mm}$01

With perturbations from $25\,\mathrm{mm}$02, $25\,\mathrm{mm}$03, lunisolar gravity, atmospheric drag, solar radiation pressure, and CW nonlinearities, the error dynamics become

$25\,\mathrm{mm}$04

and a propellant-efficient control law is

$25\,\mathrm{mm}$05

so that $25\,\mathrm{mm}$06. The stability criterion adopted is a residual control acceleration

$25\,\mathrm{mm}$07

corresponding to integrated $25\,\mathrm{mm}$08 per year or less (Ito, 2023).

Three candidate geocentric regimes are identified.

Regime	Geometry and orbit	Control budget
High Earth orbit	Triangular laser-interferometric gravitational-wave telescope, $25\,\mathrm{mm}$09 size; $25\,\mathrm{mm}$10, altitude $25\,\mathrm{mm}$11, $25\,\mathrm{mm}$12	$25\,\mathrm{mm}$13, $25\,\mathrm{mm}$14
Medium Earth orbit	Linear astronomical interferometer, baseline $25\,\mathrm{mm}$15; $25\,\mathrm{mm}$16, altitude $25\,\mathrm{mm}$17, $25\,\mathrm{mm}$18	$25\,\mathrm{mm}$19, $25\,\mathrm{mm}$20
Low Earth orbit	Demonstration interferometer, baseline about $25\,\mathrm{mm}$21–$25\,\mathrm{mm}$22; $25\,\mathrm{mm}$23, altitude $25\,\mathrm{mm}$24, sun-synchronous	$25\,\mathrm{mm}$25, $25\,\mathrm{mm}$26

The high-orbit case uses a general circular orbit with $25\,\mathrm{mm}$27, $25\,\mathrm{mm}$28, and fixed arm length $25\,\mathrm{mm}$29. Its disturbance breakdown is reported as about $25\,\mathrm{mm}$30 from $25\,\mathrm{mm}$31, $25\,\mathrm{mm}$32 from $25\,\mathrm{mm}$33, $25\,\mathrm{mm}$34 from $25\,\mathrm{mm}$35, and $25\,\mathrm{mm}$36 from $25\,\mathrm{mm}$37, with one $25\,\mathrm{mm}$38-day observation season per year under a Sun-avoidance angle of at least $25\,\mathrm{mm}$39. The medium-Earth linear array benefits from $25\,\mathrm{mm}$40-driven nodal precession of about $25\,\mathrm{mm}$41, yielding declination coverage from $25\,\mathrm{mm}$42 to $25\,\mathrm{mm}$43 over five years and eclipse durations up to $25\,\mathrm{mm}$44 per orbit near solstices. The low-Earth demonstration case has larger control cost, with drag and $25\,\mathrm{mm}$45 dominating and visibility below about $25\,\mathrm{mm}$46, but is explicitly identified as suitable for experimental purposes.

The comparison with Sun–Earth $25\,\mathrm{mm}$47 is also quantitative. Beyond-Earth orbits offer $25\,\mathrm{mm}$48, continuous illumination, and fixed geometry, whereas geocentric orbits provide proven GNSS-based formation autonomy, routine launch and transfer, and faster operational cadence. The paper’s conclusion is that geocentric space contains “sweet spots” for formation-flying interferometry from LEO through HEO, with HEO at $25\,\mathrm{mm}$49 approaching an $25\,\mathrm{mm}$50-like disturbance environment while retaining lower mission cost and complexity (Ito, 2023).

7. Comparative interpretation and recurring misconceptions

The cited literature shows that a geodesic interferometer is not synonymous with a single Michelson layout. In one usage it is a compact optical sensor for geodesy; in another it is a Michelson interferometer designed to test emergent quantum geometry; in another it is a matter-wave interferometer whose phase is determined by a shortest geodesic on the Bloch sphere; in yet another it is a GNSS–VLBI instrument for local ties and terrestrial reference frames; and in orbital studies it denotes a formation whose relative motion is maintained in a small-disturbance geocentric environment (Isleif et al., 2019, Hogan, 2012, Zhou et al., 2019, Skeens et al., 2023, Ito, 2023).

A second misconception is to treat all geodesic interferometers as measuring the same physical object. The observables are in fact heterogeneous: sub-picometer displacement and sub-$25\,\mathrm{mm}$51 tilt in DFMI; correlated broad-band displacement noise with the spectral form $25\,\mathrm{mm}$52 in the Michelson quantum-geometry proposal; hemisphere-dependent sign changes and $25\,\mathrm{mm}$53 jumps in SU(2) geometric phase; picosecond-scale group delays in GNSS–VLBI; and control acceleration and annual $25\,\mathrm{mm}$54 budgets in formation flying.

A plausible implication is that the term is best understood functionally rather than architecturally. Across all implementations, geometry enters the interferometric output through a calibration chain that must dominate technical noise: laser frequency noise, thermal path-length fluctuations, parasitic beam jitter, and electronic digitization in DFMI; shot noise and environmental decorrelation in cross-correlated Michelsons; magnetic-gradient instability and vibration-induced phase noise in matter-wave experiments; clock drift, multipath, and quantization in GNSS–VLBI; and $25\,\mathrm{mm}$55, drag, solar radiation pressure, and lunisolar perturbations in space formations. The common scientific objective is therefore not a uniform hardware standard, but a controlled transduction of geometrical structure into a measurable interferometric phase, delay, or displacement.