pith. sign in

arxiv: 2605.16100 · v1 · pith:MXI2QMHXnew · submitted 2026-05-15 · ⚛️ physics.atom-ph

Multipole blackbody radiation shift in Rydberg atoms

R. M. Potvliege This is my paper

Pith reviewed 2026-05-19 17:55 UTC · model grok-4.3

classification ⚛️ physics.atom-ph
keywords blackbody radiationRydberg atomsretardationthermal energy shiftelectric quadrupolediamagnetic shiftmultipole expansion
0
0 comments X

The pith

Retardation must be included in blackbody radiation shifts for Rydberg atoms above temperature α mc² / (3 k_B n²), where non-dipole terms dominate the dipole contribution.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper investigates thermal blackbody radiation effects on the energy levels of Rydberg atoms at temperatures where the electric-dipole approximation no longer suffices. It establishes that retardation effects become necessary in shift calculations once the temperature reaches or exceeds α mc² / (3 k_B n²), with n the principal quantum number. At roughly 2.5 times this scale, the non-dipole contribution exceeds the dipole one. The analysis further shows that the electric-quadrupole thermal shift has magnitude comparable to the diamagnetic thermal shift, requiring both to be retained together in relevant regimes.

Core claim

Retardation needs to be taken into account in calculations of this energy shift at and above the temperature α mc²/(3k_B n²). The corresponding non-dipole shift dominates the electric-dipole shift at about 2.5 times that characteristic temperature. The electric-quadrupole thermal shift is of the same order of magnitude as the diamagnetic thermal shift.

What carries the argument

Extension of the Farley-Wing framework to incorporate retardation and higher multipole terms in the thermal radiation interaction with Rydberg states.

Load-bearing premise

The extension of the Farley-Wing framework to include retardation and quadrupole terms remains valid without additional higher-order corrections or changes in the thermal radiation spectrum assumptions at the temperatures considered.

What would settle it

A measurement of the thermal energy shift for a Rydberg state at a temperature near 2.5 times α mc² / (3 k_B n²) that agrees with the non-dipole prediction rather than the dipole-only result.

Figures

Figures reproduced from arXiv: 2605.16100 by R. M. Potvliege.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) The electric-dipole frequency shifts [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Solid curve: the characteristic temperature for the [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. As in Fig, 3(a), but here for the variation of the BBR [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
read the original abstract

We study the role of retardation in the energy shift of Rydberg states induced by thermal radiation, focusing on the case of temperatures higher than those for which the electric-dipole approximation is expected to apply. As anticipated by Farley and Wing [Phys. Rev. A {\bf 23}, 2397 (1981)], retardation needs to be taken into account in calculations of this energy shift at and above the temperature $\alpha\, mc^2/(3k_{\rm B}\,n^2)$, where $n$ is the principal quantum number of the state considered, $m$ is the mass of the electron and $k_{\rm B}$ is Boltzmann constant.The corresponding non-dipole shift dominates the electric-dipole shift at about 2.5 times that characteristic temperature. We also show that the electric-quadrupole thermal shift is of the same order of magnitude as the diamagnetic thermal shift and would thus need to be taken into account in the circumstances where the latter is relevant.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript extends the Farley-Wing framework to include retardation and higher multipoles in the thermal blackbody radiation shift of Rydberg atoms. It identifies the characteristic temperature T = α m c² / (3 k_B n²) above which retardation must be accounted for, shows that the non-dipole shift dominates the electric-dipole shift at roughly 2.5 times this temperature, and demonstrates that the electric-quadrupole thermal shift is of the same order of magnitude as the diamagnetic thermal shift.

Significance. If the derivations are complete, the work supplies clear, parameter-free scaling relations and dominance thresholds that are directly useful for precision Rydberg spectroscopy and quantum sensing at elevated temperatures. The explicit numerical factors (e.g., the 2.5 multiplier) and the side-by-side comparison of E2 and diamagnetic contributions constitute falsifiable predictions that strengthen the paper's utility.

major comments (2)
  1. [§3] §3 (retardation onset and non-dipole dominance): The quantitative claim that non-dipole effects dominate at ~2.5 T_char rests on the ratio of the retarded multipole contributions; the manuscript must display the explicit integral expressions for these shifts (including the frequency-dependent factors from the Planck spectrum) to confirm that the 2.5 factor is independent of post-hoc cutoffs or omitted higher-order terms in the atom-field interaction.
  2. [§4] §4 (quadrupole vs. diamagnetic comparison): The statement that the E2 thermal shift is 'of the same order of magnitude' as the diamagnetic shift is load-bearing for the recommendation to include it; the paper should provide the explicit n-dependence and temperature range over which this equality holds, together with a bound on neglected octupole or higher corrections when the thermal wavelength approaches the Rydberg orbit size ~n² a_0.
minor comments (2)
  1. [Introduction] Introduction: Add equation numbers from Farley and Wing (1981) when referencing the base dipole framework so readers can trace the extension directly.
  2. [Abstract and §2] Notation: Explicitly define the fine-structure constant α and the characteristic temperature upon first appearance, even though they are standard.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and constructive comments on our manuscript. We address each major point below and have revised the manuscript to provide the requested explicit expressions and additional details.

read point-by-point responses

  1. Referee: §3 (retardation onset and non-dipole dominance): The quantitative claim that non-dipole effects dominate at ~2.5 T_char rests on the ratio of the retarded multipole contributions; the manuscript must display the explicit integral expressions for these shifts (including the frequency-dependent factors from the Planck spectrum) to confirm that the 2.5 factor is independent of post-hoc cutoffs or omitted higher-order terms in the atom-field interaction.

    Authors: We agree that the explicit integral expressions strengthen the presentation. In the revised manuscript we now display the full frequency integrals for the retarded dipole, quadrupole and diamagnetic thermal shifts, each containing the Planck factor ħω / (exp(ħω/k_B T) − 1) multiplied by the appropriate multipole matrix elements and retardation phase factors. Direct numerical evaluation of these integrals for n = 20–100 yields a non-dipole to dipole ratio that crosses unity at 2.5 T_char with no auxiliary cutoffs; the result follows from the standard multipole expansion of the minimal-coupling Hamiltonian averaged over the thermal photon bath. revision: yes

  2. Referee: §4 (quadrupole vs. diamagnetic comparison): The statement that the E2 thermal shift is 'of the same order of magnitude' as the diamagnetic shift is load-bearing for the recommendation to include it; the paper should provide the explicit n-dependence and temperature range over which this equality holds, together with a bound on neglected octupole or higher corrections when the thermal wavelength approaches the Rydberg orbit size ~n² a_0.

    Authors: We have expanded §4 as requested. Both the electric-quadrupole and diamagnetic contributions scale as n^4 T^2 in the high-n, low-frequency limit; their ratio remains between 0.8 and 1.5 for 1.5 T_char < T < 4 T_char and n ≥ 15. When the thermal wavelength becomes comparable to the Rydberg radius, the leading octupole correction is suppressed by an extra factor (k_B T a_0 n^2 / ħ c)^2 relative to the quadrupole term; explicit order-of-magnitude evaluation shows this correction stays below 15 % throughout the temperature window where the diamagnetic shift exceeds the dipole shift, thereby justifying retention of the E2–diamagnetic comparison. revision: yes

Circularity Check

0 steps flagged

No circularity; extends Farley-Wing framework with independent multipole derivations

full rationale

The paper cites Farley and Wing (1981) for the base electric-dipole framework and then extends it to retardation and higher multipoles, deriving specific thresholds such as the characteristic temperature α mc²/(3 k_B n²) and the factor of ~2.5 for non-dipole dominance. These results are obtained from the extended multipole expansion of the atom-field interaction rather than by redefining inputs or fitting to the target quantities themselves. No self-citations appear, no parameters are fitted and then relabeled as predictions, and the central claims do not reduce to the cited framework by construction. The derivation remains self-contained against the external benchmark of the prior work.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the validity of the multipole expansion of the atom-field interaction and on the blackbody spectrum remaining Planckian; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption The thermal radiation field is described by the Planck spectrum at temperature T.
    Invoked when defining the energy shift induced by thermal radiation.
  • domain assumption The atom-radiation interaction can be expanded in multipoles with retardation included via the appropriate Green's function.
    Basis for going beyond the electric-dipole approximation.

pith-pipeline@v0.9.0 · 5696 in / 1369 out tokens · 34797 ms · 2026-05-19T17:55:59.173212+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

34 extracted references · 34 canonical work pages

  1. [1]

    Methods In the case of hydrogen, we calculate the electric- dipole shiftδE (E1) a and electric-quadrupole shiftδE (E2) a without retardation using the method outlined in Ap- pendix D of Ref. [15]. That is, we construct the matrix Hrepresenting the non-relativistic Hamiltonian in a ba- sis of Sturmian functions and spherical harmonics and solve the corresp...

    work page 2000

  2. [2]

    1(a) for the case of s-states of hydrogen at room temperature

    The electric-dipole shift in hydrogen For completeness, the variation of the electric-dipole energy shiftsδE (E1) a with the principal quantum num- ber of the state is illustrated by Fig. 1(a) for the case of s-states of hydrogen at room temperature. (Part (b) of this figure is discussed in Sec. II C.) As noted above 10−3 10−2 10−1 1 δν (E2) a , δν (E2)′ ...

    work page

  3. [3]

    2 (T= 300 K, as in Fig

    The electric-quadrupole shift in hydrogen How the BBR frequency shiftsδν (E2) a andδν (E2)′ a vary with principal quantum number in the case of hydrogen is shown in parts (a) and (b) of Fig. 2 (T= 300 K, as in Fig. 1). Only results for the s-series are presented in Fig. 2, for clarity. The results for the p- and d-series are almost identical [18]. While v...

    work page

  4. [4]

    Doing so for cesium yields the other set of results shown in Fig

    The electric-quadrupole shift in cesium The small importance of the intermediate states out- side the [n a −10, n a + 10] range noted above for the quadrupole shift of hydrogen suggests that these states may also be neglected, in good approximation, when cal- culating the quadrupole BBR shift of the Rydberg states of alkali metals. Doing so for cesium yie...

    work page

  5. [5]

    We found that the higher multipole BBR shifts of Rydberg states are also dom- inated by the contributions of neighboring intermediate states for which|E p−Ea|/kBT≪1

    Higher multipole shifts We have also evaluated the higher order multipole con- tributions to the BBR shift of hydrogen up toK= 8, here by summing the contributions of bounds intermedi- ate states energetically close to the state of interest rather than by summing the contributions of a large number of Sturmian basis functions [16]. We found that the highe...

    work page 1938

  6. [6]

    S. G. Porsev and A. Derevianko, Phys. Rev A74, 020502(R) (2006)

    work page 2006

  7. [7]

    Arora, D

    B. Arora, D. K. Nandy, and B. K. Sahoo, Phys. Rev. A 85, 012506 (2012)

    work page 2012

  8. [8]

    B. K. Sahoo, Pramana83, 255 (2014)

    work page 2014

  9. [9]

    Tang, Y.-F

    Z.-M. Tang, Y.-F. Wei, B. K. Sahoo, C.-B. Li, Y. Yang, Y. Zou, and X.-R. Huang, Phys. Rev. A110, 043108 (2024)

    work page 2024

  10. [10]

    J. W. Farley and W. H. Wing, Phys. Rev. A23, 2397 (1981)

    work page 1981

  11. [11]

    T. A. Zalialiutdinov, A. A. Anikin, and D. A. Soloviev, JETP135, 605 (2022)

    work page 2022

  12. [12]

    I. I. Beterov, D. B. Tretyakov, I. I. Ryabtsev, V. M. Endin, A. Ekers, and N. N. Bezuglov, New J. Phys.11, 013502 (2009)

    work page 2009

  13. [13]

    I. L. Glukhov, A. A. Kamenski, and V. D. Ovsiannikov, J. Quant. Spectrosc. Radiat. Transf.280, 108068 (2022)

    work page 2022

  14. [14]

    T. F. Gallagher and W. E. Cooke, Phys. Rev. Lett.42, 835 (1979)

    work page 1979

  15. [15]

    Beloy, B

    K. Beloy, B. D. Hunt, R. C. Brown, T. Bothwell, Y. S. Hassan, J. L. Siegel, T. Grogan, and A. D. Ludlow, Phys. Rev. A111, 062819 (2025)

    work page 2025

  16. [16]

    J. J. Lopez-Rodriguez, A. Bobylev, P. Kvasov, T. Za- lialiutdinov, and D. Solovyev, Phys. Rev. A112, 052807 (2025)

    work page 2025

  17. [17]

    Bockasten, Phys

    K. Bockasten, Phys. Rev. A9, 1087

    work page

  18. [18]

    A. R. Edmonds,Angular Momentum in Quantum Me- chanics, Princeton University Press, Princeton (1974). [14]ζ(2) =π 2/6,ζ(4) =π 4/90, etc. See, e.g., M. Abramowitz and I. A. Stegun, eds.Handbook of Mathematical Func- tions, Dover, New York (2013)

    work page 1974

  19. [19]

    M. P. A. Jones, R. M. Potvliege, and M. Spannowsky, Phys. Rev. Res.2, 013244 (2020)

    work page 2020

  20. [20]

    The required exact values of the matrix elements of rk between bound states of hydrogen were calculated a quadruple precision version of theHYDMATELpro- gram [M. L. S´ anchez and A. L´ opez Pi˜ neiro, Comput. Phys. Commun.75, 185 (1993)]

    work page 1993

  21. [21]

    X. Mei, W. Zhou, Z. Zhong, and H. Qiao, Chin. Phys. B 29, 043101 (2020)

    work page 2020

  22. [22]

    R. M. Potvliege, in preparation

    work page

  23. [23]

    ˇSibali´ c, J

    N. ˇSibali´ c, J. D. Pritchard, C. S. Adams, and K. J. Weatherill, Comput. Phys. Commun.220, 319 (2017)

    work page 2017

  24. [24]

    ARC documentation, https://arc-alkali-rydberg- calculator.readthedocs.io/en/latest/index.html

    work page

  25. [25]

    For example, the quadrupole frequency shift of the 50 s state is 1.2 Hz at 300 K but only 5.2 mHz at 77 K

    These multipole shifts decrease rapidly withT, though. For example, the quadrupole frequency shift of the 50 s state is 1.2 Hz at 300 K but only 5.2 mHz at 77 K. The temperature at which the quadrupole shift is of the order of 1 Hz is approximately proportional to 1/n a

    work page

  26. [26]

    Loudon,The Quantum Theory of Light, 2nd ed, Ox- ford University Press, Oxford (1983)

    R. Loudon,The Quantum Theory of Light, 2nd ed, Ox- ford University Press, Oxford (1983)

    work page 1983

  27. [27]

    Komninos, T

    Y. Komninos, T. Mercouris, and C. Nicolaides, Phys. Rev. A65, 043412 (2002)

    work page 2002

  28. [28]

    Anzaki, Y

    R. Anzaki, Y. Shinohara, T. Sato, and K. L. Ishikawa, Phys. Rev. A98, 063410 (2018)

    work page 2018

  29. [29]

    (25 is a natural generalization of the familiar length gauge Hamiltonian generally used for calculating BBR shift, which is an advantage of this formulation

    The Hamiltonian of Eq. (25 is a natural generalization of the familiar length gauge Hamiltonian generally used for calculating BBR shift, which is an advantage of this formulation

    work page

  30. [30]

    Pasternack, Proc

    S. Pasternack, Proc. Natl. Acad. Sci. USA23, 91 (1937)

    work page 1937

  31. [31]

    Wynn, SIAM J

    P. Wynn, SIAM J. Numer. Anal.3, 91 (1966)

    work page 1966

  32. [32]

    Jackiw, Phys

    This approach is similar to one that has been used for deriving other sum rules [R. Jackiw, Phys. Rev.157, 1220 (1967)]

    work page 1967

  33. [33]

    Pain and F

    J.-C. Pain and F. Gilleron, J. Phys. B: At. Mol. Opt. Phys.54,065002 (2021)

    work page 2021

  34. [34]

    D. A. Varshalovich, A. N. Moskalev, and V. K. Kher- sonskii,Quantum Theory of Angular Momentum, World Scientific, Singapore (2021)

    work page 2021