- The paper demonstrates LO-free recovery of an RF signal’s phase and amplitude by using a static DC field to Stark-mix a near-degenerate Rydberg pair and close a phase-sensitive optical loop.
- The exact harmonic phase law makes the first demodulated probe harmonic rotate directly with the RF phase, while amplitude estimation requires an injective response map and follows RMSE laws set by harmonic SNR and local sensitivity.
- The analysis shows that phase and amplitude have different optimal bias angles—approximately 0.49 and 0.60 radians, with a balanced value near 0.56 radians—and that bias nonuniformity sharply reduces coherent gain and phase performance.
Overview and problem statement
Rydberg-atom receivers based on ladder electromagnetically induced transparency (EIT) are well established for RF field measurement, but recovering the carrier phase of a single received RF tone has historically required an auxiliary reference: either an RF local oscillator (LO) injected into the atoms (as in atomic-mixer, superheterodyne, and phase-modulation schemes) or additional time-varying optical/RF reference structure. The paper by Katkov and Zlatanov proposes a third route: a static DC bias field that Stark-mixes a near-degenerate upper Rydberg pair, activating an otherwise absent optical pathway and closing a phase-sensitive loop within a receiver driven only by the probe field, the coupling field, and the received signal itself. "LO-free" here means no auxiliary RF LO inside the vapor cell; ordinary synchronous demodulation outside the cell is still assumed.
The received signal is AScos(ωSt+ΦS) with unknown AS and ΦS at a known or separately tracked carrier ωS. The paper's contributions are: (i) the DC-Stark loop-closure mechanism; (ii) an exact harmonic phase law for the periodic steady state under uniform bias; (iii) phase and amplitude estimators from a demodulated probe harmonic; and (iv) high-SNR RMSE laws, mixing-angle design criteria, and a quasistatic treatment of bias nonuniformity.
Reduced model and loop-closure mechanism
The constructive model uses four states {∣1⟩,∣2⟩,∣3⟩,∣4⟩} with probe on ∣1⟩↔∣2⟩, coupling on ∣2⟩↔∣3⟩, and the received RF coupling the bare pair ∣3⟩↔∣4⟩. A static bias Ez mixes only the near-degenerate upper pair, producing Stark states with mixing angle
tan(2θ)=β=ℏΔ342Ez∣μ34z∣,
and effective couplings AS0, AS1, and AS2. The leg AS3 vanishes at AS4 and is the bias-enabled pathway that closes the loop. A design trade-off is visible immediately: increasing AS5 strengthens the loop-closing leg but weakens the effective RF coupling through AS6 — this competition generates the distinct phase- and amplitude-optimal angles derived later. The authors are explicit that the trigonometric forms belong to this minimal isolated-pair model; the general theory requires only well-defined effective couplings.
After a rotating-wave approximation and a gauge transformation placing the factor AS7 on the bias-enabled leg, the reduced Hamiltonian takes Floquet form with Fourier components only at AS8 and AS9. The construction neglects additional nearby Stark states, magnetic-sublevel structure, Doppler averaging, transport, and the diagonal oscillatory term from the transformed signal dipole; these omissions shift quantitative response maps but do not alter the structural phase result, provided the phase dependence remains confined to ΦS0 and the dissipator stays time independent.
Exact harmonic phase law
Because the unknown phase enters the Hamiltonian only through ΦS1, changing ΦS2 is equivalent to a time shift of the drive. Under uniqueness of the periodic steady state (PSS), each Fourier coefficient of the density operator obeys
ΦS3
This law is exact — detunings, linewidths, and dissipation modify the coefficients ΦS4 but not the multiplicative factor. For the probe coherence, harmonic magnitudes are independent of ΦS5 while phases shift linearly with slope ΦS6, so the full PSS can be computed once at ΦS7 via nearest-neighbor-coupled harmonic-balance equations (a block-tridiagonal Liouville-space system after truncation to ΦS8).
Estimators and high-SNR error laws
Phase is recovered from the argument of a demodulated harmonic relative to a calibrated reference:
ΦS9
The first harmonic (ωS0) is emphasized as branch-unambiguous; higher harmonics carry intrinsic modulo-ωS1 ambiguity. Amplitude recovery inverts the response map ωS2, which must be injective on the operating branch — otherwise inversion is ambiguous without prior information or a second observable.
Under additive circular complex Gaussian noise on the demodulated phasor, linearization yields two distinct laws:
ωS3
where ωS4 is the logarithmic sensitivity. Phase accuracy depends only on harmonic SNR; amplitude accuracy additionally depends on the local response slope, packaged into an effective amplitude SNR ωS5. These are local, high-SNR approximations that neglect phase wrapping, inverse-map ambiguity, and curvature effects.
Four mixing-angle criteria follow for ωS6: the phase-optimal angle maximizes the first-harmonic magnitude; the amplitude-optimal angle maximizes ωS7; a weighted joint-design angle minimizes a weighted sum of squared errors; and a weight-free balanced angle minimizes the maximum fractional degradation relative to the individual optima. In the perturbative regime all angular dependence collapses to ωS8, giving the seed rule ωS9 ({∣1⟩,∣2⟩,∣3⟩,∣4⟩}0) — explicitly flagged as an initialization heuristic, not a design theorem.
For {∣1⟩,∣2⟩,∣3⟩,∣4⟩}1 treated quasistatically, every local first-harmonic response carries the same factor {∣1⟩,∣2⟩,∣3⟩,∣4⟩}2, so any linear spatial average preserves the exact phase law while altering the complex coefficient. Two mechanisms act: spatial variation of the local mixing angle, and inhomogeneous RF detuning via the bias-dependent dressed splitting. The coherent-gain factor {∣1⟩,∣2⟩,∣3⟩,∣4⟩}3 reduces the effective phase SNR as {∣1⟩,∣2⟩,∣3⟩,∣4⟩}4, so the phase-RMSE penalty is exactly {∣1⟩,∣2⟩,∣3⟩,∣4⟩}5. Amplitude recovery is affected more subtly: averaging reshapes the effective response-map slope {∣1⟩,∣2⟩,∣3⟩,∣4⟩}6, and moderate nonuniformity can reduce coherent gain while simultaneously steepening the local branch — so amplitude behavior is not determined by {∣1⟩,∣2⟩,∣3⟩,∣4⟩}7 alone.
Numerical results
Production calculations truncate harmonics at {∣1⟩,∣2⟩,∣3⟩,∣4⟩}8 (validated against {∣1⟩,∣2⟩,∣3⟩,∣4⟩}9: relative error below ∣1⟩↔∣2⟩0 at the nominal point by ∣1⟩↔∣2⟩1, below ∣1⟩↔∣2⟩2 at a stress point by ∣1⟩↔∣2⟩3), with direct BDF time-domain integration confirming the Floquet reconstruction. At the nominal point (∣1⟩↔∣2⟩4, ∣1⟩↔∣2⟩5, ∣1⟩↔∣2⟩6, ∣1⟩↔∣2⟩7, balanced angle ∣1⟩↔∣2⟩8 rad), varying ∣1⟩↔∣2⟩9 over ∣2⟩↔∣3⟩0 rotates the first-harmonic phasor rigidly, with residuals at machine precision — a direct confirmation of the exact phase law. At ∣2⟩↔∣3⟩1 the first harmonic collapses to zero, consistent with the absence of loop closure.
Key numerical findings:
| Quantity |
Value |
| Full-model phase optimum ∣2⟩↔∣3⟩2 |
≈ 0.49 rad |
| Amplitude optimum ∣2⟩↔∣3⟩3 |
≈ 0.60 rad |
| Balanced angle ∣2⟩↔∣3⟩4 |
≈ 0.56 rad |
| Perturbative seed ∣2⟩↔∣3⟩5 |
≈ 0.393 rad |
| Local sensitivity at nominal point |
∣2⟩↔∣3⟩6 |
| Sensitivity under nonuniformity |
∣2⟩↔∣3⟩7 |
| Coherent gain ∣2⟩↔∣3⟩8 for ∣2⟩↔∣3⟩9 |
0.32, 0.22, 0.10 |
Two results deserve emphasis. First, the full-model phase optimum lies clearly right of the perturbative seed, confirming that ∣3⟩↔∣4⟩0 is genuinely non-optimal outside the weak-coupling regime — the authors state plainly that it should be used only for initialization. Second, Monte-Carlo RMSE curves (∣3⟩↔∣4⟩1 trials per point) match the high-SNR laws closely; under nonuniform bias, phase curves collapse when replotted against ∣3⟩↔∣4⟩2, and amplitude curves collapse only after the additional slope correction ∣3⟩↔∣4⟩3, validating the two-factor amplitude theory. Notably, the coherent-gain values are severe — even 1% relative nonuniformity costs roughly a factor of three in coherent gain in this parameter set — underscoring the authors' own conclusion that bias homogeneity and calibration are part of the receiver design problem, not secondary corrections.
Limitations and open questions
The framework is theoretical throughout, built on a deliberately minimal four-level model. It relies on the isolated-pair Stark approximation, neglects bias-induced mixing with lower ladder states and magnetic-sublevel structure, omits Doppler and transit effects, treats nonuniformity as static rather than fluctuating, discards the diagonal oscillatory dipole term, and assumes the carrier frequency is known or tracked — it is not a theory of blind carrier acquisition. Harmonic magnitudes, response maps, prefactors, and optimal angles are model dependent and will shift under manifold-level descriptions. No species-specific level selection or experimental realization is provided; the coherent-gain losses reported under even modest nonuniformity raise the concrete open question of how closely a practical vapor-cell geometry with realistic electrode fields can approach the uniform-bias benchmark.
Conclusion
The paper establishes that a static DC bias, by Stark-mixing a near-degenerate Rydberg pair, suffices to close a phase-sensitive reception loop for a single RF tone without any RF LO in the atoms. The resulting periodic steady state obeys an exact harmonic phase law, enabling direct phase estimation and injective-map amplitude recovery, with explicit high-SNR RMSE laws and four distinct bias-controlled operating-point criteria. Spatial bias nonuniformity preserves the phase law but degrades performance through coherent-gain loss and slope reshaping, quantified by collapse laws confirmed numerically. The work provides a minimal analytical foundation whose translation to specific atomic species and device geometries remains the principal outstanding step.