{"id":"bf54b254-5541-478d-a310-c1011416c2ce","arxiv_id":"2608.05905","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":20,"one_line_summary":"This paper reports high-precision J=1-0 rotational frequencies for 172, 174, and 176YbO and a three-state Morse-coupling model that reproduces the perturbed ground-state vibrational levels through v=8.","lead":"New microwave measurements of three ytterbium oxide isotopologues give sharply improved rotational constants and large isotope-dependent corrections, and a coupled-state model reproduces the molecule's irregular vibrational spacings. The results provide cleaner deperturbed potentials for low-lying states of YbO, which are useful for ongoing precision-measurement and theory work on ytterbium-containing molecules.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Morse wavefunctions used for overlaps are not eigenfunctions of the energy formula Eq. (1); the cubic anharmonic term is large at high v, so the size of the induced error in the fitted couplings and deperturbed potentials needs to be quantified before the central claim is accepted.","rationale":"The reader's weakest assumption is the Morse/cubic mismatch, and that is exactly the most load-bearing internal modeling assumption. The paper's central assertion depends on reliable overlap integrals between Morse wavefunctions for three states, but Eq. (1) adds a cubic term whose magnitude grows rapidly with v and is not contained in the Morse Hamiltonian. This is not a stylistic objection: it is a concrete inconsistency between the operator whose eigenvalues are fitted and the basis in which the coupling matrix is evaluated. The quantitative impact is testable with a standard numerical vibrational calculation, and the outcome determines whether the fitted H_E values and deperturbed potentials should be treated as quantitative or merely indicative. The FTMW rotational analysis and multi-isotope BOB work are independent and appear sound; the concern is restricted to the deperturbation model. Since the paper already presents the model as approximate and the reader has issued a CONDITIONAL verdict, no verdict change is needed, but the missing quantification should be supplied or explicitly stated as a limitation in a revision.","tokens_in":15930,"tokens_out":7483,"duration_ms":74655,"concrete_test":"Recompute the central overlaps and refit the couplings using numerical vibrational wavefunctions for the X state that are eigenfunctions of a potential reproducing the fitted Eq. (1) eigenvalues (e.g., Numerov integration of an RKR/inverted-perturbation potential through v=8), while keeping B1 and D1 as specified. Compare the resulting F_XB1(4,2), F_XD1(4,0), fitted H_E values, and X-state spacings with the Morse-based results. If any key overlap or coupling shifts by more than ~20%, the Morse assumption is load-bearing and the deperturbed potentials need quantitative error estimates; if all shifts are small, the acknowledged inconsistency is benign.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The deperturbation model's diagonal energies are generated by Eq. (1), which contains a cubic anharmonic term ωe ye (v+1/2)^3 for the X state. The off-diagonal matrix elements and all Franck-Condon factors use Eqs. (2)-(3) with Morse oscillator wavefunctions that are eigenfunctions of a Hamiltonian containing only the quadratic anharmonic term. The authors explicitly acknowledge this in Section III.C ('Morse wavefunctions (which refer to an anharmonic Hamiltonian including only the quadratic term) are assumed'), but never quantify the mismatch. For the optimized X parameters (ωe=701 cm^-1, ωexe=10.2 cm^-1, ωeye=0.34 cm^-1) the cubic term is ~209 cm^-1 at v=8, about 7% of the vibrational term value, so the wavefunction shapes entering the key X(v=4)-B1(v=2) and X(v=4)-D1(v=0) overlaps can be significantly distorted. Because the observed level shifts depend on products H_E × F, an unquantified bias in F propagates directly into the fitted electronic couplings and the reported deperturbed potentials. The RMSE of 8.4 cm^-1 and the ground-state spacing residuals therefore cannot, by themselves, establish that the extracted potentials are reliable. No error bars or stability analysis are given for the 10 manually fitted parameters, compounding the issue.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports Fourier transform microwave measurements of the N = 1←0 rotational transition for 172YbO, 174YbO, and 176YbO, together with a multi-isotopologue fit that includes Born-Oppenheimer breakdown terms. The rotational constants are combined with near-infrared combination differences from Melville et al. to improve the B and D constants. The second half of the paper presents a deperturbation model of the X 1Σ+, B1, and D1 states using coupled Morse potentials with electronic coupling matrix elements; the model is fitted to vibrational levels through X v = 8 and validated against fluorescence branching ratios. The authors conclude that the X-state spacing anomalies around v = 4 arise from near-degeneracy with B1(v = 2) and D1(v = 0).","tokens_in":16454,"tokens_out":4685,"duration_ms":41949,"significance":"If the deperturbation analysis is reliable, it provides the first quantitative account of the long-noted 3000 cm−1 anomaly in YbO and yields deperturbed potentials useful for the wider YbX family. The rotational analysis is a solid contribution: the multi-isotopologue FTMW fit is self-consistent, the combination-difference approach greatly reduces the uncertainties in B, and the BOB terms are physically plausible. The paper is also honest about its limitations, explicitly acknowledging the Morse-wavefunction mismatch and the lack of independent D1 potential information. However, the deperturbation model's central claims currently rest on approximations whose impact is not quantified, and the manual fit lacks parameter uncertainties.","major_comments":[{"comment":"The paper fits eigenvalues of Eq. (1), which includes the cubic anharmonic term ω_e y_e (v+1/2)^3, but computes vibrational overlaps and Franck-Condon factors with Morse wavefunctions that are eigenfunctions of a Hamiltonian containing only the quadratic anharmonic term. The text acknowledges this ('Morse wavefunctions ... are assumed') but never quantifies the resulting error. For the optimized X-state parameters (ω_e = 701 cm−1, ω_exe = 10.2 cm−1, ω_eye = 0.34 cm−1), the cubic term contributes roughly 209 cm−1 at v = 8, about 7% of the vibrational term value; this is large enough that the shapes of the high-v wavefunctions, and hence the overlaps entering the X(v = 4)–B1(v = 2) and X(v = 4)–D1(v = 0) couplings, may be significantly biased. Because the observed level shifts depend on products H_E × F, an unquantified bias in F propagates directly into the fitted electronic couplings and the reported deperturbed potentials. I request a quantitative estimate of this mismatch, for example by comparing the Morse-overlap results with overlaps obtained from numerical vibrational eigenfunctions of the full potential implied by Eq. (1), or by a sensitivity test on ω_eye.","section":"Section III.C, Eq. (1)"},{"comment":"The ten free parameters are optimized through a 'manual refinement process' and the paper reports no parameter uncertainties, no correlation matrix, and no stability analysis. The RMSE of 8.4 cm−1 is two orders of magnitude smaller than the level spacings, but this alone does not establish that the extracted potentials are reliable, particularly because several fixed inputs (Δr_B1, Δr_D1, ω_e,D1 scaled from theory, and ω_exe,D1 = ω_exe,B1) are taken from ab initio calculations without a reported sensitivity study. The authors should report at least one-dimensional confidence intervals from the fit, for example by scanning each parameter, and should show how the RMSE and the deperturbed potentials change when the theoretical constraints are varied within their stated uncertainties.","section":"Section III.C, Table V"},{"comment":"The model is fitted to roughly 14 observed levels with 10 free parameters, and the fit through X v = 8 excludes higher electronic states; the authors state that the B1 residuals are dominated by the neglected Ω = 0− and Ω = 2 states, but they do not quantify whether those states also affect the X-state levels through v = 8. Since the central conclusion concerns the X-state spacings, a test of the sensitivity of the X(v = 4) shifts to the inclusion of a fourth state, or at least an estimate of the expected magnitude of the neglected couplings, is needed to support the claim that the three-state model is sufficient.","section":"Section III.C, Table IV"}],"minor_comments":[{"comment":"The phrase 'To the extent that (2) is an accurate model' should read 'To the extent that Eq. (2) is an accurate model'.","section":"Section III.C"},{"comment":"The sentence 'equilibrium rotational constants B e and the corresponding bond distances are are derived' contains a duplicated 'are'; it should read 'are derived'.","section":"Table III"},{"comment":"The text contains 'the the 18 parameters originally in the model'; the duplicated article should be removed.","section":"Section III.C"},{"comment":"The reference in the text to 'Fig.III D' should be 'Fig. 4'.","section":"Section III.D"},{"comment":"The Franck-Condon validation uses an assumed upper-state r_e = 1.780 Å, and the authors note that the branching ratios are governed primarily by this choice of r_e; a short sensitivity statement showing the range of r_e that is consistent with the observed 80:10:1 and 30:1 ratios would strengthen the claimed independent validation.","section":"Section III.D"},{"comment":"The paper does not include a data availability statement or a link to the Python code used for the deperturbation fit; providing the code and input files would improve reproducibility.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The rotational and BOB portion is strong and publishable nearly as is. The deperturbation section is interesting and within the journal's scope, but the lack of uncertainty quantification for the manual fit and the unquantified wavefunction mismatch are serious enough to require substantial revision. I would not demand a full multi-state coupled-channel calculation, but the requested sensitivity analyses and error estimates are essential. Reference 36 is cited with an SSRN DOI and 'in press'; the editor should verify that this citation is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this one: the rotational/BOB half is strong and citable, and the deperturbation half is an interesting but under-supported model. The FTMW J=1-0 measurements for 172, 174, 176YbO are clean, the multi-isotopologue SPFIT analysis is sensible, and the improved B, D, and bond distance are useful. The large BOB terms (~50 kHz) are consistent with the known dense manifold of low-lying states, and the authors are appropriately cautious in not over-interpreting the field-shift vs. mass-dependent split. That part of the paper is a solid contribution.\n\nThe deperturbation model is where I part ways with the authors' confidence. The central idea—three mutually interacting Morse potentials for X, B1, and D1—is a reasonable extension of the YbF and PbO work, and the claim that a single perturbing state cannot explain the X-state spacings is well supported. But the model has a load-bearing internal inconsistency: Eq. (1) adds a cubic anharmonic term for the X state, while the wavefunctions used for all overlaps and Franck-Condon factors are pure Morse functions that belong to a Hamiltonian without that cubic term. The authors explicitly acknowledge this but never quantify the error. The stress-test estimate is sobering: at v=8 the cubic term is ~209 cm^-1, about 7% of the vibrational term value. That directly distorts the FCFs that multiply the fitted coupling terms, so the reported H_E values and deperturbed potentials inherit an unquantified bias.\n\nOn top of that, the fit is manual: 10 parameters adjusted to roughly 14 levels, with no parameter uncertainties and no stability analysis. The D1 potential is pinned by theory, and the Franck-Condon validation, while independent, rides on an assumed upper-state r_e. The RMSE of 8.4 cm^-1 is dominated by the B1 state, which the authors attribute to omitted states—fine, but it underscores that the model is phenomenological.\n\nBottom line: the rotational/BOB results are worth publishing and citing. The deperturbation model is a plausible qualitative explanation of the X-state anomaly, but the deperturbed potentials and couplings should not be treated as quantitative until the Morse/cubic mismatch is fixed or bounded, and ideally with error bars from a least-squares fit. I'd send this to a serious referee, with a demand for major revision on the deperturbation section. I'd cite the rotational constants and BOB terms; I'd cite the deperturbation only with a caveat.","headline":"Solid new FTMW rotational data and BOB terms for YbO; the three-state deperturbation model is a plausible first pass but has a load-bearing Morse/cubic inconsistency that needs to be quantified before the deperturbed potentials are trusted.","tokens_in":16944,"tokens_out":1831,"would_cite":true,"duration_ms":19875,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that YbO's known irregularities in ground-state vibrational spacings, particularly the v = 4 anomaly near 3000 cm⁻¹, are caused by two nearly degenerate excited-state levels, and that a three-state effective Hamiltonian…","keywords":["ytterbium oxide","rotational spectroscopy","Fourier transform microwave","Born-Oppenheimer breakdown","de-perturbation","Morse potential","Franck-Condon factor","vibrational perturbations"],"falsifier":"Measure the $D_1(v=1)$ and $D_1(v=2)$ vibrational levels directly with high-resolution near-infrared or laser spectroscopy; the de-perturbed $D_1$ potential from the model predicts their positions, and a deviation of more than a few wavenumbers would falsify that potential. A second, calculation-based check is to recompute the 4518 Å fluorescence branching ratios using numerical (non-Morse) eigenfunctions of the fitted potentials to see whether the 80:10.4:0.7 and 30:1.1 predictions survive the change in wavefunction model.","tokens_in":15723,"feed_emoji":"📡","tokens_out":9872,"duration_ms":82715,"temperature":0.7,"pith_summary":"This paper reports high-precision Fourier transform microwave measurements of the $J = 1 \\leftarrow 0$ rotational transition in three even-mass isotopologues of ytterbium oxide, YbO, and combines them with near-infrared data to sharpen the ground-state rotational and centrifugal distortion constants. The main explanatory claim is that the ground state's irregular vibrational spacings, especially the abnormally small $v = 4 \\rightarrow 3$ interval near 3000 cm$^{-1}$, cannot be produced by a single perturbing state. The authors build a three-state effective Hamiltonian in which the $X\\,^1\\Sigma^+$ ground state mixes with the $B_1$ and $D_1$ excited states through electronic coupling matrix elements, with Morse potentials generating the vibrational overlaps. Diagonalizing this matrix reproduces the measured $X$-state levels through $v = 8$ with an RMS deviation of 8.4 cm$^{-1}$, yielding de-perturbed potentials for all three states. If correct, the model explains the long-noted anomaly and provides de-perturbed potentials for studying other ytterbium-containing molecules.","feed_headline":"Two electronic states explain YbO's ground-state anomaly","feed_subtitle":"Microwave measurements plus a three-state Morse model reproduce the v = 4 spacing oddity within experimental error.","key_machinery":"The central object is an effective Hamiltonian matrix built in the basis of vibrational levels of the three electronic states $X$, $B_1$, and $D_1$. Diagonal terms are Morse-based vibrational energies with an added cubic anharmonic term, $\\omega_e y_e (v+1/2)^3$, needed for the higher ground-state levels; off-diagonal terms are the product of a state-independent electronic coupling $H^E_{s s'}$ and a vibrational overlap $F_{s s'} = \\langle v_{s} | v_{s'} \\rangle$ computed numerically from Morse wavefunctions. The model carries the argument because diagonalizing this matrix simultaneously reproduces the perturbed $X$-state spacings, locates the two perturbing states, and yields wavefunctions whose Franck-Condon factors match the observed fluorescence intensities.","core_discovery":"The central discovery is that the perturbed ground-state spacings of YbO, and in particular the small $v = 4 \\rightarrow 3$ interval, are a three-state resonance: $X(v=4)$ sits nearly degenerate with $B_1(v=2)$ and $D_1(v=0)$, and both excited levels plus a $B_1$–$D_1$ coupling are needed to reproduce the data. After fixing equilibrium bond-length shifts from theory and scaling the $D_1$ vibrational frequency, the model has 10 free parameters and fits all measured levels through $v=8$ for $X$, through $v=3$ for $B_1$, and the origin of $D_1$ with an RMS deviation of 8.4 cm$^{-1}$, two orders of magnitude smaller than the level spacing. The authors state that the model \"reproduces the observed perturbations in the ground state levels to well within the measurement uncertainties.\" An independent check is that the de-perturbed Morse wavefunctions predict fluorescence branching ratios of 80:10.4:0.7 to the $X$-state progression and 30:1.1 to the $B_1$ progression, matching the measured 80:10:1 and 30:1.","pith_inferences":["The use of Morse wavefunctions for overlaps while fitting a cubic-term Hamiltonian is the least-tested piece of the model; a numerical potential fitted to the same levels could show whether the Franck-Condon validation is sensitive to this inconsistency.","Because the $D_1$ vibrational parameters are fixed by theory rather than data, the $D_1(v=1)$ prediction is the model's sharpest testable consequence; a miss there would not necessarily invalidate the $X$-state perturbation picture.","The authors note that the $B_1$-state residuals dominate the RMS error and suspect missing couplings to the lowest $\\Omega = 0^-$ and $\\Omega = 2$ states; adding those states to the matrix could shift the extracted $B_1$ origin and slightly alter the $X$–$B_1$ coupling."],"forward_implications":["The de-perturbed $X$-state potential ($\\omega_e = 701$, $\\omega_e x_e = 10.2$, $\\omega_e y_e = 0.34$ cm$^{-1}$) can be used in place of effective constants for modeling YbO spectra and computing molecular properties.","The $B_1$ and $D_1$ potentials predict the positions of unobserved vibrational levels, particularly $D_1(v=1,2)$, offering concrete targets for new infrared and optical searches.","The large Born-Oppenheimer breakdown terms of 45–50 kHz measured for the heavier isotopologues are consistent with the strong electronic mixing and will need to be included in any isotopologue-independent model of YbO.","The three-state coupled-Morse approach extends the de-perturbation method previously applied to YbF and PbO, providing a template for disentangling the overlapping f-hole states that complicate other lanthanide oxides."],"supporting_citations":[{"why":"Identified the low-lying states of YbO and first noted the anomalous $v=4$ ground-state spacing that the present work explains.","marker":"[4]"},{"why":"Reported the low-lying excited-state energies and the 4518 Å fluorescence intensity ratios used to validate the model's Franck-Condon factors.","marker":"[6]"},{"why":"Supplies the relativistic DFT equilibrium bond-length shifts and $D_1$ vibrational frequency ratio used to fix parameters in the model.","marker":"[8]"},{"why":"Near-infrared combination differences that, combined with the new FTMW data, determine the improved $B$ and $D$ constants.","marker":"[12]"},{"why":"The YbF de-perturbation study whose coupled-Morse-state modeling approach is extended here.","marker":"[43]"},{"why":"The recent PbO de-perturbation analysis used as a template for handling multiple mutually perturbing low-lying states.","marker":"[36]"},{"why":"The least-squares fitting program used for the multiple-isotopologue rotational analysis.","marker":"[47]"},{"why":"The Python Morse-oscillator code used to construct vibrational wavefunctions and compute overlap integrals.","marker":"[55]"}],"fun_headline_variants":["Three-state resonance fixes YbO ground-state spacings","YbO's v=4 oddity traced to B1 and D1 states","Microwave data plus Morse model explain YbO perturbations","De-perturbed YbO potentials match all measured levels","YbO ground-state anomaly solved with three-state model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes Morse-shaped vibrational wavefunctions for the overlaps even though the fitted energy levels come from a Hamiltonian with an extra cubic anharmonic term beyond the Morse form, and the size of the resulting mismatch is never quantified.","fun_headline_variants_meta":{"raw":{"variants":["Three-state resonance fixes YbO ground-state spacings","YbO's v=4 oddity traced to B1 and D1 states","Microwave data plus Morse model explain YbO perturbations","De-perturbed YbO potentials match all measured levels","YbO ground-state anomaly solved with three-state model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000205,"raw_usage":{"total_tokens":1416,"prompt_tokens":988,"completion_tokens":428,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":604,"completion_tokens_details":{"reasoning_tokens":340}},"tokens_in":604,"tokens_out":428,"duration_ms":4679,"temperature":1.0,"reasoning_tokens":340,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T21:23:08.400873+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $D_1(v=1)$ and $D_1(v=2)$ vibrational levels directly with high-resolution near-infrared or laser spectroscopy; the de-perturbed $D_1$ potential from the model predicts their positions, and a deviation of more than a few wavenumbers would falsify that potential. A second, calculation-based check is to recompute the 4518 Å fluorescence branching ratios using numerical (non-Morse) eigenfunctions of the fitted potentials to see whether the 80:10.4:0.7 and 30:1.1 predictions survive the change in wavefunction model.","supporting_citations":[{"cited_title":"and Tarbutt, Michael R","cited_arxiv_id":null,"evidence_quote":"Reported the low-lying excited-state energies and the 4518 Å fluorescence intensity ratios used to validate the model's Franck-Condon factors."},{"cited_title":"and Goodridge, Damian M","cited_arxiv_id":null,"evidence_quote":"The recent PbO de-perturbation analysis used as a template for handling multiple mutually perturbing low-lying states."},{"cited_title":"and Schaller, S","cited_arxiv_id":null,"evidence_quote":"The least-squares fitting program used for the multiple-isotopologue rotational analysis."}],"review_version":1}