REVIEW 3 major objections 5 minor 51 references
Rydberg Atom Electric Field Sensors as Linear Time-invariant Systems
T0 review · 3 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read A Rydberg atom electric field sensor can be treated as a linear time-invariant system, yielding an analytic transfer function that reproduces full numerical frequency responses up to two orders of magnitude faster.
desk verdict Useful, honest LTI transfer-function paper for Rydberg sensor engineers, with a repairable gap in the formal perturbation expansion. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the transfer matrix $G(\omega) \equiv (i\omega I - A)^{-1} B \vec{\rho}_{\mathrm{SS}}$, assembled from three pieces: the constant matrix $A$ of the linearized Liouville-von Neumann flow, the sparse matrix $B$ that couples only the Rydberg-Rydberg transition, and the steady-state density vector $\vec{\rho}_{\mathrm{SS}}$ set by the local oscillator and optical fields. The linearization step, in which the RF signal enters only through the first-order term $\Omega_{\mathrm{RF}}^{(1)}(t) B \vec{\rho}_{\mathrm{SS}}$, is what converts the sensor into an LTI system with $C=I$ and $D=0$. Doppler broadening is handled by writing $G(\omega;u)$ with velocity-shifted detunings and integrating over the Maxwell distribution, which avoids hundreds of separate numerical integrations.
What would settle it
Simulate one fixed atomic configuration with both the LTI transfer function and full time-dependent density-matrix integration while steadily increasing the ratio of RF signal amplitude to local oscillator amplitude; the point where the two responses diverge by more than the numerical fitting error marks the linearity ceiling, and if that ceiling lies below the signal strengths used in practice, the LTI model would not cover the intended operating regime.
Extended reading notes
Core claim
The central claim is that the first-order perturbative response of the atom vapor, $\tilde{\vec{\rho}}^{(1)}(\omega) = (i\omega I - A)^{-1} B \vec{\rho}_{\mathrm{SS}} \tilde{\Omega}_{\mathrm{RF}}^{(1)}(\omega)$, is the complete linear response of the sensor, so the probe transmission obeys $\tilde{\rho}_{12}^{(1)}(\omega_{\mathrm{IF}}) = G(\omega_{\mathrm{IF}}) \tilde{\Omega}_{\mathrm{sig}}(\omega_{\mathrm{IF}})$ with $G(\omega_{\mathrm{IF}}) = [G_{1,2}(\omega) - G_{2,1}(\omega)]/2i$. The paper demonstrates that this transfer function, Doppler-averaged over thermal velocity classes, agrees with full numerical integration of the equations of motion across the tested intermediate-frequency band, and that the same impulse response can be convolved with a 16QAM waveform to reconstruct the received constellation, including under noise. This is presented as a new application: casting the atomic receiver as a textbook LTI system with standard state-space matrices, so that receiver engineering tools apply directly to Rydberg sensors.
Load-bearing premise
The approach stands or falls on the assumption that the neglected second-order term $\epsilon^2 \Omega_{\mathrm{RF}}^{(1)}(t) B \vec{\rho}^{(1)}(t)$ is negligible, i.e., that the RF signal is weak enough relative to the local oscillator that the linear response dominates; the paper states this premise but gives no quantitative bound on when it fails.
Editorial extensions
If this is right
- The frequency response of a Rydberg sensor can be computed in roughly a minute instead of hours, making parameter scans over detunings, Rabi frequencies, and decay rates practical.
- Any weak RF waveform can be propagated through the sensor by convolution with the impulse response, so modulated formats such as QAM can be evaluated without solving the full equations of motion.
- Key receiver metrics—instantaneous bandwidth, sensitivity, and linear-regime dynamic range—follow directly from $G(\omega_{\mathrm{IF}})$ and its extensions.
- Doppler averaging becomes part of the analytic transfer function rather than a costly separate average over many numerical runs.
- Because the sensor is now expressed in standard state-space form, established systems-engineering analyses of stability, noise, and saturation apply directly to Rydberg receivers.
Reading between the lines
- The paper leaves implicit that the rational form of $G(\omega)$ means its poles carry the sensor's transient time scales; pole-zero analysis could be used to design equalization filters that compensate the nonlinear phase response noted in the paper.
- Since the paper explicitly sets aside the $\epsilon^2 \Omega_{\mathrm{RF}}^{(1)} B \vec{\rho}^{(1)}$ term, carrying that term to second order would produce predictions of intermodulation products and harmonic distortion, quantifying dynamic range in the nonlinear regime.
- The unquantified linearity ceiling could be turned into a measured specification by sweeping the signal-to-local-oscillator amplitude ratio and recording where the LTI prediction first departs from full numerical integration.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents an analytic linear-time-invariant (LTI) model for the response of a Rydberg-atom electric field sensor operated in a heterodyne configuration. The authors linearize the Liouville-von Neumann equation for a 4+1 level density matrix around the steady-state operating point, derive a state-space representation, and obtain the frequency-domain transfer function G(ω) = (iωI−A)^{-1} B ρ_SS. They validate the predicted probe-transmission response against full time-domain integration using RydIQule for one set of parameters, report a speedup of roughly two orders of magnitude, and demonstrate the approach on a 16QAM waveform.
Significance. The contribution is practically valuable: if the transfer function is accurate, it replaces hours of numerical integration with a minute-scale matrix evaluation and enables fast estimates of sensor bandwidth and waveform distortion. The derivation from Eq. 4 to Eq. 18 is straightforward, and the use of an independent numerical solver (RydIQule) as a benchmark is appropriate. The paper's reproducible parameter list in Appendix B is a strength. The main weaknesses are the incomplete definition of the small expansion parameter in Section II and the limited, amplitude-normalized validation; these are repairable and do not invalidate the underlying approach.
major comments (3)
- [Section II, Eqs. (5)-(8)] The perturbative expansion is not anchored to a physical small parameter. Eq. (7) sets E_RF^(1)(t)=E_LO cos(ω_IF t+φ), so the first-order input has the same amplitude as the zeroth-order LO term, and Eq. (8) contains a second-order term proportional to E_LO that does not vanish as the signal amplitude goes to zero. Consequently, the neglect of the ε^2 term in Eq. (11) is not justified by E_sig/E_LO being small. Please replace E_LO by E_sig in Eq. (7) (and drop or re-derive Eq. (8)) or otherwise define ε explicitly and show that the neglected term is O(ε^2) with ε=E_sig/E_LO.
- [Appendix B, "Transfer function model", Fig. 3] The numerical comparison is amplitude-normalized ("We normalize the amplitude response so that the DC response is set to unity"), which removes exactly the information needed to test the linear scaling of the output with the signal amplitude. I recommend reporting the un-normalized transfer functions or a quantitative relative-error curve, and repeating the comparison for at least two additional ratios E_sig/E_LO (e.g., 0.001 and 0.1) to map the linear regime boundary.
- [Section III, Fig. 3, Table II] The validation is limited to a single parameter set and no quantitative error metric is reported between the LTI and RydIQule transfer functions; the speedup claim rests on one wall-clock comparison (two hours versus a little over a minute). At minimum, state that the speedup is anecdotal, or add a second configuration with a different detuning or Rabi frequency to demonstrate generality.
minor comments (5)
- [Eq. (21)] The definition of f(u) contains an integral and no dependence on u; it should read f(u)=e^{-u^2}/√π. The Doppler-averaging formula in Eq. (24) uses the correct form.
- [Section VI heading] The appendix is titled "Apprendix VI" in the main text; please correct the spelling and numbering.
- [Table II] In Table II, the transit dephasing rate is listed as Γt but the text and Eq. (2) use γt; unify the notation throughout.
- [Fig. 3 caption and text] The phase response is described as "wrapped and shifted vertically," but the vertical shift is unspecified; please provide the unwrapped or shifted data so the reader can reproduce the comparison.
- [Reference [14]] Reference [14] is a Ph.D. thesis without a university or year; please complete the citation.
Circularity Check
No circularity: the LTI transfer function is derived algebraically from the Lindblad equation and checked against an independent numerical integration of the same equation, with no fitted parameters or load-bearing self-citation.
full rationale
The paper's derivation chain is self-contained in the relevant sense. Equations (5)-(13) linearize the Liouville-von Neumann equation by expanding the RF Rabi frequency and density matrix in a formal Taylor parameter ε, and Eq. (18) follows by direct algebra: G(ω) = (iωI - A)^{-1} B ρ_SS. No parameter is fitted to reproduce the later numerical results; A, B and ρ_SS are fixed by the model parameters in Table II. The comparison to RydIQule (Fig. 3) is an independent numerical integration of the same master equation, so it validates the LTI approximation without being a fitted-input-called-prediction. The 16QAM demonstration convolves a synthetic waveform with the derived impulse response; its constellation outputs are not used to adjust the model. Self-citations (e.g., Ref. 5) are contextual and not load-bearing, and no uniqueness theorem or ansatz is imported from the authors' prior work. One substantive correctness concern exists but is not circularity: Eq. (7) writes the first-order field as E_LO cos(ω_IF t) rather than E_sig cos(ω_IF t), leaving ε's relation to the physical signal-to-LO ratio unstated; this weakens the formal justification for dropping the ε^2 term but does not make the result equivalent to its inputs by construction.
Assumptions & free parameters
free parameters (3)
- gamma_t (transit dephasing rate) =
2π * 0.1 MHz
- gamma_dsg (dummy state decay) =
2π * 0.5 MHz
- gamma_dsr (dummy state relaxation) =
2π * 0.1 MHz
assumptions (5)
- domain assumption Rotating wave approximation is valid for the field detunings, decay rates, and Rabi frequencies used.
- domain assumption The atomic vapor obeys a Markovian Lindblad master equation with the specified decay and dephasing operators.
- domain assumption The five-level model with a dummy state captures the relevant physics of the Rydberg sensor.
- domain assumption The incident RF signal is weak enough that first-order perturbation theory is accurate and higher-order terms in Eq. 11 can be neglected.
- domain assumption Doppler averaging can be represented by a one-dimensional Gaussian velocity distribution with fixed temperatures.
invented entities (1)
-
Dummy state |5⟩
Cite this review
Pith. "Pith review of Rydberg Atom Electric Field Sensors as Linear Time-invariant Systems." pith.science (2026). https://pith.science/paper/BWHWYZTC
@misc{pith2026250500159,
author = {Pith},
title = {Pith review of: Rydberg Atom Electric Field Sensors as Linear Time-invariant Systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/BWHWYZTC}},
note = {Machine review of arXiv:2505.00159}
}
read the original abstract
Over the past decade, Rydberg atom electric field sensors have been under investigation as potential alternatives or complements to conventional antenna-based receivers for select applications in RF communications, remote sensing, and precision metrology. To understand the potential utility of these devices for various use cases, it is crucial to develop models that accurately predict key performance metrics such as instantaneous bandwidth and dynamic range. However, existing numerical models require solving a large set of coupled differential equations that is computationally intensive and lengthy to solve. We present an analytic approach that can be used to derive an impulse response function that allows up to two orders-of-magnitude reduction in computation time compared to the full time-dependent integration of the equations of motion. This approach can be used to enable rapid assessments of the Rydberg sensor's response to various waveforms.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
G. O. Young, ``Synthetic structure of industrial plastics,'' in Plastics, 2 nd ed., vol. 3, J. Peters, Ed. New York, NY, USA: McGraw-Hill, 1964, pp. 15--64
work page 1964
-
[2]
Chen, Linear Networks and Systems
W.-K. Chen, Linear Networks and Systems. Belmont, CA, USA: Wadsworth, 1993, pp. 123--135
1993
- [3]
-
[4]
E. P. Wigner, ``Theory of traveling-wave optical laser,'' Phys. Rev., vol. 134, pp. A635--A646, Dec. 1965, DOI 10.1109.XXX.123456
work page 1965
-
[5]
E. H. Miller, ``A note on reflector arrays,'' IEEE Trans. Antennas Propagat., to be published
-
[6]
E. E. Reber, R. L. Michell, and C. J. Carter, ``Oxygen absorption in the earth's atmosphere,'' Aerospace Corp., Los Angeles, CA, USA, Tech. Rep. TR-0200 (4230-46)-3, Nov. 1988
1988
-
[7]
J. H. Davis and J. R. Cogdell, ``Calibration program for the 16-foot antenna,'' Elect. Eng. Res. Lab., Univ. Texas, Austin, TX, USA, Tech. Memo. NGL-006-69-3, Nov. 15, 1987
1987
-
[8]
Transmission Systems for Communications, 3 rd ed., Western Electric Co., Winston-Salem, NC, USA, 1985, pp. 44--60
work page 1985
Show all 51 references
-
[9]
Motorola Semiconductor Data Manual, Motorola Semiconductor Products Inc., Phoenix, AZ, USA, 1989
1989
-
[10]
G. O. Young, ``Synthetic structure of industrial plastics,'' in Plastics, vol. 3, Polymers of Hexadromicon, J. Peters, Ed., 2 nd ed. New York, NY, USA: McGraw-Hill, 1964, pp. 15-64. [Online]. Available: http://www.bookref.com
1964
-
[11]
Kurland and Ralph Lerner, eds., Chicago, IL, USA: Univ
The Founders' Constitution, Philip B. Kurland and Ralph Lerner, eds., Chicago, IL, USA: Univ. Chicago Press, 1987. [Online]. Available: http://press-pubs.uchicago.edu/founders/
1987
-
[12]
ZOmega Terahertz Corp., 2014
The Terahertz Wave eBook. ZOmega Terahertz Corp., 2014. [Online]. Available: http://dl.z-thz.com/eBook/zomega\_ebook\_pdf\_1206\_sr.pdf . Accessed on: May 19, 2014
2014
-
[13]
Kurland and Ralph Lerner, eds., The Founders' Constitution
Philip B. Kurland and Ralph Lerner, eds., The Founders' Constitution. Chicago, IL, USA: Univ. of Chicago Press, 1987, Accessed on: Feb. 28, 2010, [Online] Available: http://press-pubs.uchicago.edu/founders/
1987
-
[14]
J. S. Turner, ``New directions in communications,'' IEEE J. Sel. Areas Commun., vol. 13, no. 1, pp. 11-23, Jan. 1995, DOI 10.1109.XXX.123456
1995
-
[15]
W. P. Risk, G. S. Kino, and H. J. Shaw, ``Fiber-optic frequency shifter using a surface acoustic wave incident at an oblique angle,'' Opt. Lett., vol. 11, no. 2, pp. 115--117, Feb. 1986, DOI 10.1109.XXX.123456
1986
-
[16]
Kopyt et al., ``Electric properties of graphene-based conductive layers from DC up to terahertz range,'' IEEE THz Sci
P. Kopyt et al., ``Electric properties of graphene-based conductive layers from DC up to terahertz range,'' IEEE THz Sci. Technol., to be published. DOI: 10.1109/TTHZ.2016.2544142
2016
-
[17]
Intranets: Internet technologies deployed behind the firewall for corporate productivity
PROCESS Corporation, Boston, MA, USA. Intranets: Internet technologies deployed behind the firewall for corporate productivity. Presented at INET96 Annual Meeting. [Online]. Available: http://home.process.com/Intranets/wp2.htp
-
[18]
R. J. Hijmans and J. van Etten, ``Raster: Geographic analysis and modeling with raster data,'' R Package Version 2.0-12, Jan. 12, 2012. [Online]. Available: http://CRAN.R-project.org/package=raster
2012
-
[19]
Lytera UG, Kirchhain, Germany [Online]
Teralyzer. Lytera UG, Kirchhain, Germany [Online]. Available: http://www.lytera.de/Terahertz\_THz\_Spectroscopy.php?id=home , Accessed on: Jun. 5, 2014
2014
-
[20]
U.S. House. 102 nd Congress, 1 st Session. (1991, Jan. 11). H. Con. Res. 1, Sense of the Congress on Approval of Military Action. [Online]. Available: LEXIS Library: GENFED File: BILLS
1991
-
[21]
Musical toothbrush with mirror, by L.M.R. Brooks. (1992, May 19). Patent D 326 189 [Online]. Available: NEXIS Library: LEXPAT File: DES
1992
-
[22]
D. B. Payne and J. R. Stern, ``Wavelength-switched pas- sively coupled single-mode optical network,'' in Proc. IOOC-ECOC, Boston, MA, USA, 1985, pp. 585--590, DOI 10.1109.XXX.123456
1985
-
[23]
Ebehard and E
D. Ebehard and E. Voges, ``Digital single sideband detection for interferometric sensors,'' presented at the 2 nd Int. Conf. Optical Fiber Sensors, Stuttgart, Germany, Jan. 2-5, 1984
1984
-
[24]
Brandli and M
G. Brandli and M. Dick, ``Alternating current fed power supply,'' U.S. Patent 4 084 217, Nov. 4, 1978
1978
-
[25]
J. O. Williams, ``Narrow-band analyzer,'' Ph.D. dissertation, Dept. Elect. Eng., Harvard Univ., Cambridge, MA, USA, 1993
1993
-
[26]
Kawasaki, ``Parametric study of thermal and chemical nonequilibrium nozzle flow,'' M.S
N. Kawasaki, ``Parametric study of thermal and chemical nonequilibrium nozzle flow,'' M.S. thesis, Dept. Electron. Eng., Osaka Univ., Osaka, Japan, 1993
1993
-
[27]
Harrison, private communication, May 1995
A. Harrison, private communication, May 1995
1995
-
[28]
Smith, ``An approach to graphs of linear forms,'' unpublished
B. Smith, ``An approach to graphs of linear forms,'' unpublished
-
[29]
Brahms, ``Representation error for real numbers in binary computer arithmetic,'' IEEE Computer Group Repository, Paper R-67-85
A. Brahms, ``Representation error for real numbers in binary computer arithmetic,'' IEEE Computer Group Repository, Paper R-67-85
-
[30]
IEEE Criteria for Class IE Electric Systems, IEEE Standard 308, 1969
1969
-
[31]
Letter Symbols for Quantities, ANSI Standard Y10.5-1968
1968
-
[32]
Fardel, M
R. Fardel, M. Nagel, F. Nuesch, T. Lippert, and A. Wokaun, ``Fabrication of organic light emitting diode pixels by laser-assisted forward transfer,'' Appl. Phys. Lett., vol. 91, no. 6, Aug. 2007, Art. no. 061103, DOI. 10.1109.XXX.123456
2007
-
[33]
Zhang and N
J. Zhang and N. Tansu, ``Optical gain and laser characteristics of InGaN quantum wells on ternary InGaN substrates,'' IEEE Photon. J., vol. 5, no. 2, Apr. 2013, Art. no. 2600111, DOI. 10.1109.XXX.123456
2013
-
[34]
Azodolmolky et al., Experimental demonstration of an impairment aware network planning and operation tool for transparent/translucent optical networks,'' J
S. Azodolmolky et al., Experimental demonstration of an impairment aware network planning and operation tool for transparent/translucent optical networks,'' J. Lightw. Technol., vol. 29, no. 4, pp. 439--448, Sep. 2011, DOI. 10.1109.XXX.123456
2011
-
[35]
author author M. T. \ Simons , author A. B. \ Artusio-Glimpse , author C. L. \ Holloway , author E. Imhof , author S. R. \ Jefferts , author R. Wyllie , author B. C. \ Sawyer , \ and\ author T. G. \ Walker ,\ title title Continuous radio-frequency electric-field detection thro...
-
[36]
author author S. M. \ Bohaichuk , author D. Booth , author K. Nickerson , author H. Tai , \ and\ author J. P. \ Shaffer ,\ title title Origins of rydberg-atom electrometer transient response and its impact on radio-frequency pulse sensing , \ 10.1103/PhysRevApplied.18.034030 j...
-
[37]
Schmidt , author S
author author M. Schmidt , author S. Bohaichuk , author V. Venu , author F. Christaller , author C. Liu , author F. Ripka , author H. K\" u bler , \ and\ author J. P. \ Shaffer ,\ title title Rydberg-atom-based radio-frequency sensors: amplitude-regime sensing , \ 10.1364/OE.5...
-
[38]
Holloway , author M
author author C. Holloway , author M. Simons , author A. H. \ Haddab , author J. A. \ Gordon , author D. A. \ Anderson , author G. Raithel , \ and\ author S. Voran ,\ title title A multiple-band rydberg atom-based receiver: Am/fm stereo reception , \ 10.1109/MAP.2020.2976914 j...
-
[39]
author author C. T. \ Fancher , author D. R. \ Scherer , author M. C. S. \ John , \ and\ author B. L. S. \ Marlow ,\ title title Rydberg atom electric field sensors for communications and sensing , \ 10.1109/TQE.2021.3065227 journal journal IEEE Transactions on Quantum Enginee...
2021
-
[40]
author author M. T. \ Simons , author A. B. \ Artusio-Glimpse , author A. K. \ Robinson , author N. Prajapati , \ and\ author C. L. \ Holloway ,\ title title Rydberg atom-based sensors for radio-frequency electric field metrology, sensing, and communications , \ https://doi.or...
-
[41]
author author D. H. \ Meyer , author B. N. \ Miller , author T. Virtanen , author C. M. \ O'Brien , \ and\ author K. C. \ Cox ,\ title title Rydiqule: Graph-based modeling of atomic systems , \ in\ 10.1364/CLEO_FS.2024.FTu4L.4 booktitle CLEO 2024 \ ( publisher Optica Publishin...
2024 doi
-
[42]
William L
author author J. William L. Hallauer ,\ @noop title Introduction to Linear, Time-invariant systems, Dynamic Systems for Students of Engineering \ ( publisher A.T. Still University ,\ year 2016 ) NoStop
2016
-
[43]
author author D. F. \ Delchamps ,\ @noop title State Space and input-output linear systems \ ( publisher Springer-Verlag ,\ year 1988 ) NoStop
1988
-
[44]
author author J. P. \ Hespanha ,\ @noop title Linear Systems Theory \ ( publisher Princeton Press ,\ address Princeton, New Jersey ,\ year 2018 )\ note iSBN13: 9780691179575 NoStop
2018
-
[45]
author author M. T. \ Simons , author A. H. \ Haddab , author J. A. \ Gordon , \ and\ author C. L. \ Holloway ,\ title title A Rydberg atom-based mixer: Measuring the phase of a radio frequency wave , \ 10.1063/1.5088821 journal journal Applied Physics Letters \ volume 114 ,\ ...
-
[46]
Jing , author Y
author author M. Jing , author Y. Hu , author J. Ma , author H. Zhang , author L. Zhang , author L. Xiao , \ and\ author S. Jia ,\ title title Atomic superheterodyne receiver based on microwave-dressed rydberg spectroscopy , \ 10.1038/s41567-020-0918-5 journal journal Nature P...
-
[47]
author author D. H. \ Meyer , author Z. A. \ Castillo , author K. C. \ Cox , \ and\ author P. D. \ Kunz ,\ title title Assessment of rydberg atoms for wideband electric field sensing , \ 10.1088/1361-6455/ab6051 journal journal Journal of Physics B: Atomic, Molecular and Optic...
-
[48]
author author R. C. \ Brown , author B. Kayim , author M. A. \ Viray , author A. R. \ Perry , author B. C. \ Sawyer , \ and\ author R. Wyllie ,\ title Very-high- and ultrahigh-frequency electric-field detection using high angular momentum Rydberg states ,\ 10.1103/PhysRevA.107...
-
[49]
author author J. A. \ Sedlacek ,\ @noop title Microwave and surface electrometry with rydberg atoms , \ ( year 2016 ) NoStop
2016
-
[50]
Vajpayee , author V
author author V. Vajpayee , author V. Becerra , author N. Bausch , author J. Deng , author S. Shimjith , \ and\ author A. J. \ Arul ,\ title title Dynamic modelling, simulation, and control design of a pressurized water-type nuclear power plant , \ https://doi.org/10.1016/j.nu...
-
[51]
author author W. L. \ Luyben ,\ @noop title Chemical reactor design and control \ ( publisher John Wiley & Sons, Inc. ,\ year 2007 ) NoStop
2007
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.