REVIEW 3 major objections 5 minor 27 references
The double minimum E(3)$^1\Sigma^{+}_{\mathrm{u}}$ state in Cs$_2$
T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read The paper claims that a potential energy curve built from 6,727 measured transitions reproduces the energies of 5,631 levels of the double-minimum E(3)1Σ+u state of Cs2 with an rms deviation of 0.043 cm-1, while showing that the outer…
desk verdict A substantial, honest Cs2 E-state line list and fit; the PEC is a good data compressor but the barrier position and absolute v numbering are conditional on the theoretical outer well. 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 nearest-energy strategy, in which for each experimental term energy and a given J, the code searches for the closest level calculated from the current potential instead of relying on an assigned vibrational quantum number. This is implemented in a single-channel Fourier grid Hamiltonian (FGH) code, and the potential is refined with the inverted perturbation approach (IPA). The strategy converges only when the initial potential is fairly good, which is why the theoretical PEC of Spies is used as the starting point; the adiabatic single-channel model is justified by the Dressler criterion parameter gamma > 35 for the E(3)1Σ+u and 4 1Σ+u pair.
What would settle it
If a future experiment observed transitions to levels supported by the outer well (for example, from a different electronic state or through a different excitation scheme), and these levels deviated from the predictions of the Table I potential by more than the stated accuracy, the claim of a quantitatively reliable potential would be falsified. Alternatively, a measurement involving another isotopologue of caesium (though none is stable) could fix the absolute vibrational numbering and expose errors in the numbering used here.
Extended reading notes
Core claim
The central claim is that the IPA/FGH potential energy curve in Table I, built from the nearest-energy strategy and starting from the theoretical potential of Spies, reproduces the measured energies of levels in the inner well and above the internal barrier with an rms deviation of 0.043 cm-1, within the experimental accuracy of 0.05 cm-1. The barrier is located near 5.33 Å at about 20195 cm-1. The authors further show that four different refitted potentials (U9, U18, U27, Up) match the data with rms near 0.045 cm-1 while carrying different vibrational numberings, proving that the absolute v numbering and the outer-well shape are not uniquely fixed by the data.
Load-bearing premise
The entire construction rests on the assumption that the initial theoretical potential from Spies is close enough to the true potential that the nearest-energy search converges to the correct correspondence between calculated and observed levels, and that the single-channel adiabatic model accurately describes the state.
Editorial extensions
If this is right
- The term values and potential provide a benchmark for testing ab initio calculations of the double minimum state and its barrier.
- The measured frequencies in the range up to 700 cm-1 above the barrier can be used to plan photoassociation experiments to form ultracold Cs2 molecules.
- The Dunham-type coefficients in Table II allow compact calculation of above-barrier levels for v' = 101-171 and J' = 29-198.
- The demonstration that several PECs fit the data equally well clarifies that the outer well shape and absolute numbering are undetermined without outer-well data.
Reading between the lines
- The ambiguity in the outer well may affect estimates of the long-range behavior of the potential, which could matter for photoassociation rates, so theoretical long-range coefficients should be used with caution until outer-well data are obtained.
- The nearest-energy strategy combined with controlled numerical experiments on the unconstrained part of the potential is a transferable protocol for assessing the reliability of double minimum potential fits; it could be applied to other molecules with a single stable isotopologue.
- If a scheme to populate outer-well levels (e.g., via optical Raman transfer or in a different state) becomes available, the current potential can be tested and refined, potentially fixing the absolute numbering.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. Jastrzebski et al. report a spectroscopic study of the double-minimum E(3)1Σ+u state of Cs2 using polarization labelling spectroscopy. They identify 6727 rotationally resolved E←X transitions, convert them into 5642 term values (including 92 levels from Amiot et al.), and construct a single-channel adiabatic potential energy curve in Table I using the Fourier Grid Hamiltonian method with the nearest-energy strategy and an IPA-type inversion. The potential reproduces 5631 of the 5642 term values with an rms deviation of 0.043 cm-1, covering the inner well and energies up to about 700 cm-1 above the internal barrier. The outer well is fixed to the theoretical PEC of Spies. Section V reports numerical experiments showing that several alternative PECs (U9, U18, U27, Up) fit the same data with similar rms while differing in the outer-well shape and absolute vibrational numbering. The paper concludes that the outer well and parts of the potential above the barrier cannot be determined uniquely from the present data.
Significance. The data set and the fitted PEC provide a quantitative characterization of the inner well and above-barrier region of the E(3)1Σ+u state, which is of genuine interest for photoassociation experiments. The paper's strengths include a large, carefully assigned data set; the use of well-established IPA/FGH methods; an explicit discussion of the non-uniqueness problem; and the stated intention to make supplementary data available. The reproduction of 5631 levels at 0.043 cm-1 rms is a legitimate fit-quality result, and the authors do not market it as an independent prediction. The numerical distortion experiments in Section V are valuable sensitivity checks that probe the uniqueness of the fit. The main caveat is that the derived barrier position and absolute vibrational numbering depend on the theoretically supplied outer well, as the authors themselves demonstrate.
major comments (3)
- [Section V and Table I] The numerical experiments in Section V show that the experimental data do not select a unique PEC: U9, U18, U27 and Up all achieve rms near 0.045 cm-1 yet place the outer-well minimum and parts of the above-barrier potential at different positions and assign different absolute vibrational numbers (e.g., the lowest inner-well level is v=28 in Table I/U9, v=27 in U18/U27, and v=31 in Up). Since the abstract and the Table I caption present the curve as the IPA potential of the E state, and Fig. 2 labels the vibrational axis with the absolute numbering from Table I, the paper overstates the uniqueness of these results. The authors should explicitly state that Table I is one member of an equivalence class of potentials consistent with the data, quantify the spread in the barrier position/height and absolute numbering across the U-family, and adjust the claims in the abstract and introduction accordingly.
- [Section IV, Table I] Table I states that 5631 of 5642 measured levels are reproduced, leaving 11 levels unexplained, but these 11 levels are never identified or discussed. Without knowing whether they are unassigned lines, perturbed levels, or outliers, the reader cannot assess the quality of the fit or the completeness of the data set. The authors should list these levels (at least their v, J and term values or their positions in the supplementary data) and provide a reason for their exclusion from the fit.
- [Table I and Section IV] No uncertainties are quoted for the PEC grid points. The experimental accuracy is 0.05 cm-1 and the fit rms is 0.043 cm-1, so the potential is determined to roughly this precision, but the paper should provide at least the standard errors of the fitted grid points, or of the derived barrier parameters (position and height). Quantifying this uncertainty is particularly important because the numerical experiments demonstrate that different acceptable fits exist and the spread among them is a direct measure of the systematic uncertainty.
minor comments (5)
- [Abstract] The word 'unambiguousness' should be replaced by 'uniqueness' or 'non-uniqueness' for clarity.
- [Section III] The relationship between 6727 spectral lines and 5642 term values should be stated explicitly (e.g., P/R doublets from different lower levels reaching the same upper level), since only a brief mention is given in the text.
- [Figure 3 caption] The caption 'parts of four potentials' should specify which four potentials are shown (Table I, U9, U18, U27) to avoid ambiguity.
- [Equation (1)] The summation limits in Eq. (1) are not defined; the range of m and n should be stated, or the reader should be referred to the table.
- [References] Reference [5] is a 1990 PhD thesis; if more recent theoretical PECs for Cs2 exist, citing one would help the reader assess the quality of the starting potential.
Circularity Check
No significant circularity: the PEC is fitted to the data it reproduces, but the paper does not present this agreement as an independent prediction and explicitly documents the non-uniqueness of the outer well.
full rationale
The paper constructs an IPA/FGH potential that reproduces measured term values with rms 0.043 cm−1, which is a fit statistic rather than an independent prediction; the authors do not frame the agreement as a test of a first-principles result. The ambiguous absolute vibrational numbering and the uncertain outer-well shape are openly discussed in Section V, including numerical experiments (U9, U18, U27, Up) that yield different outer wells and different vibrational numberings with comparable rms. This is an honest statement of non-uniqueness, not a circular derivation. Citations to the authors' earlier methodological work and to Spies' theoretical PEC are used as tools or starting points, not as imported uniqueness theorems or as proof of the final potential. No equation is defined in terms of the quantity it purports to derive, and no fitted parameter is renamed as a prediction. The central physical claims are explicitly conditioned on the unobserved outer well, which is a correctness/robustness caveat rather than circularity.
Assumptions & free parameters
free parameters (3)
- PEC grid point values (Table I) =
28 rotationless grid points, e.g. U(5.33 Å)=20195.4541 cm-1
- Dunham coefficients A_mn (Table II) =
A00=25759.435 cm-1, A10=-227.183 cm-1, etc.
- Distortion amplitude a in numerical experiments =
2, 4, and 6 cm-1 Å-2
assumptions (5)
- domain assumption Single-channel Born-Oppenheimer adiabatic model is valid for the E(3)1Σ+u state.
- ad hoc to paper The initial theoretical PEC of Spies [5] is accurate enough for the nearest-energy assignment.
- domain assumption Ground X-state constants from Amiot et al. [2] can be used to convert transition wavenumbers to term values.
- domain assumption Franck-Condon factors prevent observation of outer-well levels.
- standard math Natural cubic spline interpolation between the grid points of Table I.
Cite this review
Pith. "Pith review of The double minimum E(3)$^1\Sigma^{+}_{\mathrm{u}}$ state in Cs$_2$." pith.science (2026). https://pith.science/paper/MHMGNPEP
@misc{pith2026250105271,
author = {Pith},
title = {Pith review of: The double minimum E(3)$^1\Sigma^+_\mathrmu$ state in Cs$_2$},
year = {2026},
howpublished = {\url{https://pith.science/paper/MHMGNPEP}},
note = {Machine review of arXiv:2501.05271}
}
abstract
The double minimum E$^1\Sigma^{+}_{\mathrm{u}}$ state in caesium dimer was investigated by analysing spectra of the E$^1\Sigma^{+}_{\mathrm{u}}$ $\leftarrow$ X$^1\Sigma^{+}_{\mathrm{g}}$ band system, simplified by polarisation labelling. A total of 6727 rotationally resolved transitions to levels situated in the inner well and above the internal barrier were identified and a potential energy curve allowing to reproduce their energies was constructed using the Fourier grid Hamiltonian and inverted perturbation approach methods. The unambiguousness of the potential curve in view of lack of data related to levels located in the outer well is discussed.
Figures
Reference graph
Works this paper leans on
- [1]
- [2]
-
[3]
M. Pichler, W. C. Stwalley, R. Beuc, and G. Pichler, For- mation of ultracold Cs 2 molecules through the double- minimum Cs 2 31Σ+ u state, Phys. Rev. A69, 013403 (2004)
work page 2004
-
[4]
T. Ban, S. Ter-Avetisyan, R. Beuc, H. Skenderovi´ c, and G. Pichler, Photoassociation of cesium atoms into the double minimum Cs2 31Σ+ u state, Chem. Phys. Lett. 313, 7 /s56 /s49/s48 /s49/s50 /s49/s52 /s49/s46/s57/s56 /s50/s46/s48/s55 /s50/s46/s49/s54 /s85/s40/s82/s41/s120/s49 /s48 /s48 /s48 /s48 /s32/s99/s109 /s45/s49 /s82/s32/s91/s65/s93 /s32/s111/s114/...
work page 1999
-
[5]
The double minimum E(3)$^1\Sigma^{+}_{\mathrm{u}}$ state in Cs$_2$
we constructed a potential energy curve of the E 1Σ+ u state allowing to reproduce energies of the rovibrational levels involved in the observed transitions with accuracy better than 0.05 cm −1. The position and height of the barrier on the potential curve was determined. As a par- ticular benefit, we provide the measured frequencies of transitions betwee...
work page Pith review arXiv 2025
-
[6]
N. Spies, Theoretische Untersuchung von elektronisch an- geregten Zust¨ anden der Molek¨ ule Li2 und Cs2, PhD The- sis, Universit¨ at Kaiserslautern, 1990
work page 1990
- [7]
- [8]
Show all 27 references
-
[9]
Nishimiya, Y
N. Nishimiya, Y. Yasuda, T. Yukiya, and M. Suzuki, Sub- Doppler polarization spectroscopy of the B 1Πu ← X1Σ+ g system for Cs 2 with a titanium sapphire ring laser, J. Mol. Spectrosc. 255, 194–198 (2009)
2009
-
[10]
Pashov, P
A. Pashov, P. Kowalczyk, A. Grochola, J. Szczepkowski, and W. Jastrzebski, Coupled-channels analysis of the (51Σ+ u , 51Πu, 53Πu, 23∆u) complex of electronic states in rubidium dimer, J. Quant. Spectrosc. Radiat. Transfer, 22, 225–232 (2018)
2018
-
[11]
Lefebvre-Brion and R
H. Lefebvre-Brion and R. W. Field, The Spectra and Dynamics of Diatomic Molecules , Elsevier, Amsterdam, 2004
2004
-
[12]
Dressler, in Photophysics and Photochemistry above 6 eV, ed
K. Dressler, in Photophysics and Photochemistry above 6 eV, ed. F. Lahmani, p. 327, Elsevier, Amsterdam, 1985
1985
-
[13]
Amiot and O
C. Amiot and O. Dulieu, The Cs 2 ground electronic state by Fourier transform spectroscopy: Dispersion co- efficients, J. Chem. Phys. 117, 5155–5164 (2002)
2002
-
[14]
Kosman and J
W. Kosman and J. Hinze, Inverse perturbation analy- sis: Improving the accuracy of potential energy curves, J. Mol. Spectrosc. 56, 93–103 (1975)
1975
-
[15]
Vidal and H
C. Vidal and H. Scheingraber, Determination of diatomic molecular constants using an inverted perturbation ap- proach: Application to the A 1Σ+ u –X1Σ+ g system of Mg 2, J. Mol. Spectrosc. 65, 46–64 (1977)
1977
-
[16]
Pashov, W
A. Pashov, W. Jastrz¸ ebski, and P. Kowalczyk, Construc- tion of potential curves for diatomic molecular states by the IPA method, Comp. Phys. Commun. 128, 622–634 (2000)
2000
-
[17]
J. Y. Seto, R. J. Le Roy, J. Verg` es, and C. Amiot, Direct potential fit analysis of the X 1Σ+ g state of Rb 2. Nothing else will do!, J. Chem. Phys. 113, 3067–3076 (2000)
2000
-
[18]
Modified Lennard-Jones oscillator
P. G. Hajigeorgiou and R. J. Le Roy, A “Modified Lennard-Jones oscillator” model for diatom potential functions, J. Chem. Phys. 112, 3949–3957 (2000)
2000
-
[19]
Samuelis, E
C. Samuelis, E. Tiesinga, T. Laue, M. Elbs, H. Kn¨ ockel, and E. Tiemann, Cold atomic collisions studied by molec- ular spectroscopy, Phys. Rev. A63, 012710 (2000)
2000
-
[20]
Pashov, W
A. Pashov, W. Jastrz¸ ebski, W. Ja´ sniecki, V. Bednarska, and P. Kowalczyk, Accurate potential curve of the double minimum 21Σ+ u state of Na2, J. Mol. Spectrosc. 203, 264– 267 (2000)
2000
-
[21]
Jastrz¸ ebski, W
W. Jastrz¸ ebski, W. Ja´ sniecki, P. Kowalczyk, R. Nadyak, and A. Pashov, Spectroscopic investigation of the dou- ble minimum 21Σ+ u state of potassium dimer, Phys. Rev. A62, 042509 (2000)
2000
-
[22]
Grochola, P
A. Grochola, P. Kowalczyk, W. Jastrzebski, and A. Pashov, A regularized inverted perturbation approach method: potential energy curve of the 41Σ+ u state in Na2, J. Chem. Phys. 121, 5754–5760 (2004)
2004
-
[23]
Pashov, W
A. Pashov, W. Jastrz¸ ebski, and P. Kowalczyk, Improved description of the double minimum 6 1Σ+ state of NaK by an IPA potential energy curve, J. Phys. B: At. Mol. Opt. Phys. 33, L611–L614 (2000)
2000
-
[24]
Szczepkowski, A
J. Szczepkowski, A. Grochola, P. Kowalczyk, W. Jas- trzebski, E. A. Pazyuk, A. V. Stolyarov, and A. Pashov, The spin-orbit coupling of the 6 1Σ+ and 4 3Π states in KCs: observation and deperturbation, J. Quant. Spec- trosc. Radiat. Transfer, 239, 106650 (2019)
2019
-
[25]
C. C. Marston and G. G. Balint-Kurti, The Fourier grid Hamiltonian method for bound state eigenvalues and eigenfunctions, J. Chem. Phys. 91, 3571–3576 (1989)
1989
-
[26]
De Boor, A Practical Guide to Splines , Springer, Berlin, 1978
C. De Boor, A Practical Guide to Splines , Springer, Berlin, 1978
1978
-
[27]
Pashov, P
A. Pashov, P. Kowalczyk, and W. Jastrzebski, Absolute vibrational numbering from isotope shifts in fragmentary spectroscopic data, J. Mol. Spectrosc. 347, 48–55 (2018)
2018
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.