{"id":"015195b6-7236-4f2e-b043-cc47ed377f75","arxiv_id":"2608.06739","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Simulated Raman spectra of 13C16O2 and 16O13C18O reproduce the main experimental bands, and the apparent discrepancy between two experiments is attributed to rotational structure in the energy levels.","lead":"This paper simulates Raman spectra of two carbon dioxide isotopologues using algebraic vibrational wave functions and polarizability derivatives fitted to the main isotopologue, then compares the results with two independent experiments. A smart generalist might read it because isotope-specific Raman spectra at high temperature matter for gas sensing, combustion diagnostics, and geological carbon isotope analysis.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The rotational-structure explanation for the 570 K discrepancy is not established: the new comparison uses different temperature, pressure, and SNR, and no anisotropic calculation is made.","rationale":"The reader's weakest assumption was the interpretation of the discrepancy, and I partially agree. The deeper issue is that the paper's own comparison is not sufficient to separate the rotational and anisotropic hypotheses. The model itself has real supporting evidence: the 298 K 13CO2 spectrum reproduces relative intensities, and Tables 2 and 3 show many transition moments within experimental uncertainty. I do not lean on the rms inconsistency emphasized by the reader, because the reported rms values appear consistent with a relative-residual definition (Table 3's 0.108 matches the relative RMS of its residuals), so that criticism should be moderated. The 500/650 K experiment is a reasonable idea but is not a controlled test: the claim that acquisition time alone explains the absence or presence of rotational structure is physically implausible as stated, and the full Alvarez window is never fitted quantitatively. The conclusion that anisotropic effects are unimportant is especially undersupported because no anisotropic calculation appears in the paper. The model may still be correct, but the resolution claimed in Sections 5 and 6 needs a full-window quantitative fit and an explicit anisotropic test. This is consistent with the reader's CONDITIONAL verdict, so I would not change it.","tokens_in":22512,"tokens_out":11549,"duration_ms":121847,"concrete_test":"Re-analyze the Alvarez 570 K spectrum over the full 1200-1460 cm^-1 range with the authors' ro-vibrational simulation (Eq. (22) with Eq. (26), using B_nu constants from Refs. [69-72]) and compute a reduced chi-squared or normalized band-by-band residuals. Then repeat the same full-window comparison with an added anisotropic polarizability term, taking anisotropy derivatives from Ref. [68] or from a fresh ab initio calculation. If the isotropic ro-vibrational model fails on bands outside 1375-1385 cm^-1, or if adding the anisotropic term changes the simulated contour by more than a few percent, the rotational-vs-anisotropic conclusion is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sections 5 and 6 conclude that the discrepancy between the Alvarez 570 K spectrum and the authors' 500/650 K spectra is rotational structure, and that anisotropic scattering is comparatively unimportant. This is the load-bearing link between the spectral agreement and the claimed physical resolution. The support is incomplete in three concrete ways. (i) The two experiments differ not only in acquisition time but also in temperature (570 vs 500/650 K), pressure (1 MPa in the new capillary vs unspecified in Ref. [61]), sample composition, and instrumentation; any of these can alter line contours. (ii) The only demonstration of the rotational mechanism is the narrow 1375-1385 cm^-1 window in Figure 6; no full-window fit or residual is reported for the Alvarez 1200-1460 cm^-1 spectrum. (iii) The conclusion about anisotropic scattering is reached without computing or fitting any anisotropic contribution. The new spectra are unpolarized, so both isotropic and anisotropic parts contribute; matching an unpolarized spectrum with an isotropic-only model does not by itself rule out a significant anisotropic component. Finally, the physical claim that a 300 s x 3 acquisition records only 'most probable vibrational transitions' while a 50 min acquisition detects rotational transitions is not a mechanism: longer integration improves SNR but does not gate out rotational transitions. The central claim, as formulated, is therefore not demonstrated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports simulations of the Raman spectra of 12C16O2, 13C16O2, and 16O13C18O using vibrational wave functions from the SU1(2)×U(3)×SU2(2) algebraic model, with mean-polarizability derivatives fitted to 42 experimental transition moments of the principal isotopologue. The simulations are compared with experimental spectra from Álvarez et al. (Ref. [61]) and with new measurements at 500 and 650 K. The authors argue that the discrepancy between the 570 K spectrum of Ref. [61] and their simulated spectra arises from rotational structure rather than from anisotropic scattering, and conclude that the model provides transferable Raman intensities for CO2 isotopologues.","tokens_in":22864,"tokens_out":4548,"duration_ms":38947,"significance":"If the central claim holds, the paper offers a practical route to predicting high‑temperature Raman spectra of CO2 isotopologues from a fit to the principal isotopologue, with applications to isotope‑ratio measurements and combustion diagnostics. The comparison of predicted transition moments for 13C16O2 and 16O13C18O with the independent experimental moments in Tables 2 and 3 is a genuine cross‑check, not a fit to the target data, because the derivatives are determined from the principal isotopologue. The inclusion of diagonal ro‑vibrational energy corrections is a useful methodological improvement for spectral shape. However, the strongest conclusion—that the observed experimental differences are due to rotational structure and that anisotropic contributions are unimportant—is not established by the evidence presented.","major_comments":[{"comment":"The claim that the discrepancy between the 570 K spectrum of Ref. [61] and the authors' spectra is due to rotational structure, rather than to anisotropic scattering, is not supported by the evidence. The two experiments differ in temperature (570 K vs 500/650 K), pressure (1 MPa in the new capillary vs unspecified in Ref. [61]), sample composition, and instrumentation; any of these can alter line contours. The only demonstration of the rotational mechanism is the narrow 1375–1385 cm⁻¹ window in Fig. 6; no full‑window fit or residual is reported for the 1200–1460 cm⁻¹ spectrum of Ref. [61]. Moreover, the statement that a 300 s × 3 acquisition records only the 'most probable events provided by the vibrational transitions' while a 50 min acquisition detects ro‑vibrational transitions is not a mechanism: longer integration improves signal‑to‑noise but does not gate out rotational transitions. Finally, the conclusion that anisotropic contributions are 'not so important' is reached without computing or fitting any anisotropic term; the new spectra are unpolarized, so both isotropic and anisotropic parts contribute, and agreement with an isotropic‑only model does not by itself rule out a significant anisotropic component.","section":"§5, Fig. 6 and §6"},{"comment":"The abstract's claim of 'excellent agreement with experiment' is overstated. The simulation for 12C16O2 at 1780 K shows 'significant discrepancies in the intensities ... in the range 1360–1500 cm⁻¹' as explicitly acknowledged in §5. Tables 2 and 3 also contain residuals that exceed the experimental uncertainty by a large factor: Table 2 band 13 (|2000⟩→|3000⟩) has Δ = −4.19 with σ = 3.00, band 19 (|1200⟩→|1400⟩) has Δ = −3.95 with σ = 0.54, and Table 3 band g has Δ = 2.12 with σ = 0.79. These outliers are not discussed in the text and undermine the statement that the transition moments are 'well determined.'","section":"§5, Fig. 1 and abstract"},{"comment":"The polarizability expansion is truncated at cubic terms with only the five derivatives listed in Table 1, and the sufficiency of this truncation is assumed without justification. In particular, the anisotropic part of the polarizability is omitted entirely, so the model cannot by itself distinguish isotropic from anisotropic contributions. The conclusion in §6 that anisotropic effects are smaller than suggested in Ref. [68] is therefore not a result of the present calculation; at minimum, the authors should explicitly state that this conclusion is conditional on the assumed expansion and on the neglect of anisotropy.","section":"§3, Eq. (23) and Table 1"},{"comment":"The ro‑vibrational stick spectrum in Fig. 7 is computed with a maximum rotational quantum number Jmax = 100, but no justification is given for this cutoff or for the neglect of the J‑dependent Boltzmann factor beyond the diagonal energy correction in Eq. (26). The improvement in Fig. 6 is demonstrated only in a narrow spectral window; a quantitative comparison over the full 1200–1460 cm⁻¹ range, with residuals or an R² measure, is needed to support the claim of a 'remarkable improvement.'","section":"§5, Fig. 7 and Eq. (22)"}],"minor_comments":[{"comment":"Equation (1) of the experimental setup contains apparent encoding errors: '25 ţm' and '−133 řC' should read '25 µm' and '−133 °C'.","section":"§4"},{"comment":"The dynamical group is written as 'SU1(2) × SU(3) × SU2(2)' in §6 but as 'SU1(2) × U(3) × SU2(2)' throughout the rest of the paper; please make the notation consistent.","section":"§6 vs. §1–§3"},{"comment":"The tables state that the rms deviation is calculated using 'the definition given in Eq. (22) of Ref. [60]', but the formula is not reproduced in this manuscript and Ref. [60] is a preprint on SSRN; please provide the explicit definition of the rms residual.","section":"Table 2 and Table 3"},{"comment":"The sentence 'the experimental spectrum of Álvarez et al. display the transition lines' should be 'displays', and the grammatical errors in the third paragraph of §5 ('we readily identify our simulations closely fitted') should be corrected.","section":"§5"},{"comment":"The caption does not specify which simulated spectrum (red or blue) corresponds to which isotopologue in the left display; the labeling should be defined explicitly in the caption.","section":"Figure 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript depends heavily on Ref. [60], which is an SSRN preprint; if that preprint is not yet peer‑reviewed, its status should be clarified. The comparison with Ref. [68] (which includes anisotropic contributions) is cursory; a more quantitative engagement with that work would strengthen the paper. The authors mention that a fit of the derivatives using isotopologue data is 'in progress'; incorporating such a fit, or at least a sensitivity analysis, would make the transferability claim more robust."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe short version: this is a useful but over-claimed paper. The genuinely new pieces are the two new experimental Raman spectra of a 13C16O2/16O13C18O mixture at 500 and 650 K, and the demonstration that adding diagonal rotational corrections to the Raman shift dramatically improves the simulated contour in a narrow 1375-1385 cm^-1 window. The comparison of the model's transition moments to the independent experimental moments from Álvarez et al. is a real predictive test, and the agreement is mostly decent. The target audience—people using Raman for CO2 isotope ratios or high-temperature diagnostics—will find that cross-check informative.\n\nBut the paper needs work before I'd trust its conclusions. The abstract says 'excellent agreement' while Figure 1 shows clear intensity discrepancies for 12C16O2 in the 1360-1500 cm^-1 range; the text later admits 'significant discrepancies.' The reported rms values in Tables 2 and 3 are not the plain rms of the residuals shown—several residuals are multiple sigma (band 19 in Table 2 is -3.95 with a 0.54 uncertainty), and the definition referenced from Ref. [60] doesn't obviously reconcile the numbers. That needs to be spelled out.\n\nThe load-bearing claim—that the discrepancy between the Álvarez 570 K spectrum and the authors' own spectra comes from rotational structure rather than anisotropic scattering—is plausible but not established. The two experiments differ in temperature, pressure, and integration time; the 'rotational structure' demonstration is limited to one narrow window with no quantitative residual for the full spectrum; and the conclusion about anisotropic scattering is reached without computing any anisotropic contribution. Matching an unpolarized spectrum with an isotropic-only model doesn't rule out a significant anisotropic part. Also, the idea that a 300 s x 3 acquisition 'records only pure vibrational transitions' while a 50-minute acquisition detects rotational transitions is not physically precise: longer integration improves SNR, not resolution. The data should be deposited, and the agreement of the new spectra should be quantified.\n\nIn short: the new experiment and the ro-vibrational correction are real contributions, and the transition moment cross-check is worth a referee's time. But the paper overstates its case, and the rotational-structure/anisotropy conclusion needs either a full-spectrum quantitative analysis or a more modest framing. I'd send it to review, with the expectation of major revision.\n\nBest,\n[Your name]","headline":"Useful new data and a plausible rotational correction, but the paper overstates agreement and does not support its conclusion about anisotropic scattering.","tokens_in":23401,"tokens_out":5973,"would_cite":true,"duration_ms":53041,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["33.20.Fb","33.15.Mt"],"model":"deepseek-v4-flash","headline":"The SU1(2) × U(3) × SU2(2) algebraic model of CO2, with polarizability derivatives fitted only to the principal isotopologue, reproduces the high-temperature Raman spectra of 13C16O2 and 16O13C18O, and attributes differences between two…","keywords":["Raman spectroscopy","carbon dioxide isotopologues","algebraic vibrational model","mean polarizability","transition moments","Born-Oppenheimer transferability","ro-vibrational corrections","high-temperature Raman"],"falsifier":"A controlled experiment on the same isotopic mixture at a single temperature, with integration time varied from 300 s to 50 minutes and all other conditions fixed, would settle the claim: if longer acquisition does not grow the rotational structure around 1380.46 $\\mathrm{cm}^{-1}$ and the other discrepant line, the rotational explanation fails. A separate check is a polarized Raman measurement that isolates the anisotropic contribution; if the depolarized component contributes more than the few percent the authors assume, the dismissal of anisotropic scattering would be wrong.","tokens_in":22278,"feed_emoji":"🔬","tokens_out":8621,"duration_ms":80027,"temperature":0.7,"pith_summary":"The paper tries to establish that Raman intensities are transferable across carbon dioxide isotopologues: polarizability derivatives fitted once to the most abundant isotopologue, combined with algebraic vibrational wave functions, are enough to simulate the high-temperature Raman spectra of 13C16O2 and 16O13C18O. The simulations agree with experiment for the transition moments, a 298 K spectrum, and 570 K spectra, with root-mean-square residuals of $0.130$ and $0.108$ in units of $10^{-42}\\,\\mathrm{CV}^{-1}\\mathrm{m}^2$ in the two moment tables. The paper also claims to explain why two independent 570 K measurements look different: one records rotational structure on top of the vibrational bands, the other does not, and this difference comes from rotational energy corrections, not from anisotropic scattering. If correct, one calibrated polarizability surface plus isotope-specific wave functions would be a practical route to Raman-based isotope analysis of CO2 at flame and reservoir temperatures.","feed_headline":"Fitted polarizability predicts Raman spectra of CO2 isotopologues","feed_subtitle":"A model fitted on ordinary CO2 intensities transfers to heavier isotopologues and explains a 570 K spectral discrepancy through rotation.","key_machinery":"The central object is the dynamical group SU1(2) × U(3) × SU2(2), in which the two stretching modes are described by SU(2) ladder operators built from Morse-type anharmonic oscillators and the degenerate bending by the U(3) model with a fixed boson number. Vibrational eigenstates are built in a local basis, projected onto $D_{\\infty h}$ symmetry, and diagonalized with a polyad-preserving Hamiltonian. The intensity machinery is the mean polarizability $\\bar{\\alpha}$ expanded in curvilinear coordinates to cubic order, realized algebraically through canonical and anharmonic mappings; its derivatives, such as $(\\partial\\bar{\\alpha}/\\partial S_g)_0$, are fit parameters determined once from 42 experimental transition moments of the main isotopologue and transferred to other isotopologues under the Born–Oppenheimer approximation. The simulation uses the isotropic differential-cross-section formula and applies the key correction $(B_{\\nu'}-B_\\nu)[J(J+1)-\\ell^2]$ to the Raman shift, so state-dependent rotational constants broaden and shift the computed lines.","core_discovery":"The central claim is that the SU1(2) × U(3) × SU2(2) dynamical-group description of CO2, whose stretching modes are Morse-like SU(2) oscillators and whose bending is a U(3) model, yields vibrational wave functions of spectroscopic quality, with root-mean-square deviations of $0.06$ and $0.07~\\mathrm{cm}^{-1}$ for the two isotopologues. When the mean polarizability is expanded in curvilinear coordinates through cubic terms and its derivatives are fitted to 42 experimental transition moments of the principal isotopologue, the same derivatives produce transition moments for 13C16O2 and 16O13C18O that mostly agree within experimental uncertainty. The visible mismatches at 1275.03 and 1380.46 $\\mathrm{cm}^{-1}$ are then traced to the omission of rotational structure: adding the diagonal term $(B_{\\nu'}-B_\\nu)[J(J+1)-\\ell^2]$ to the Raman shifts, without changing the wave functions, transforms the simulated band shapes and intensities into the measured ones. From this the authors conclude that the difference between the two experimental spectra reflects rotational content, not anisotropic scattering, and that a faithful Raman simulation needs the rotational energy structure at least through diagonal corrections.","pith_inferences":["If the transferability claim holds, the same fit could be used to simulate the Raman spectra of 13C18O2 and 17O13C18O, for which vibrational term values are sparse, making the model a predictive tool rather than an interpolation scheme; the paper does not report those spectra yet.","The acquisition-time explanation implies a directly testable prediction: recording the same gas at the same temperature with integration times from 300 s to 50 minutes should show rotational branches growing in; the paper's two experiments also differ in temperature, so a controlled time series would separate the two variables.","A natural next step is to promote the rotational constants to state-dependent fitted values from a full ro-vibrational effective Hamiltonian; the authors note that line positions would then sit exactly on the experimental peaks.","The method could be combined with spatial-resolution Raman measurements to map carbon isotope ratios in heterogeneous gas samples, an application the paper cites as motivation but does not pursue."],"forward_implications":["The same fitted polarizability derivatives, combined with wave functions for any of the nine isotopologues, allow Raman spectra to be predicted for species with sparse experimental data, including the rare symmetric isotopologues.","Raman simulations intended for quantitative comparison should include diagonal rotational corrections to transition energies even when the wave functions are computed in the rigid-rotor approximation.","Anisotropic scattering, previously suggested as a major source of error in CO2 Raman simulations, is argued to contribute less than about five percent and is not needed to explain the tested spectra.","The two discrepant spectral features are explained by band-specific rotational-constant differences: a wide interval in $B_\\nu$ produces a broad band, a narrow interval a sharp peak, so intensity mismatches can be diagnosed from rotational constants alone.","Transition moments computed for hot bands up to 21,400 $\\mathrm{cm}^{-1}$ extend the energy range over which Raman thermometry calibrations for CO2 could be constructed."],"supporting_citations":[{"why":"Supplies the 42 experimental transition moments of the principal isotopologue used to fit the polarizability derivatives, and the 1780 K benchmark spectrum.","marker":"[57]"},{"why":"Provides the experimental transition moments and the 298 K and 570 K spectra of 13C16O2 and 16O13C18O that the simulations are compared against.","marker":"[61]"},{"why":"Provides the vibrational term values and wave functions for the principal and symmetric isotopologues used in the simulations.","marker":"[54]"},{"why":"Supplies the vibrational energies and wave-function fits for 13C16O2 used in the computed Raman shifts.","marker":"[55]"},{"why":"Supplies the vibrational description and term values for the asymmetric isotopologue 16O13C18O.","marker":"[56]"},{"why":"Establishes the earlier algebraic Raman simulation and the polyad choice used here.","marker":"[59]"},{"why":"Describes the fitting method for polarizability derivatives and the extrapolated polarizability surfaces on which the transfer relies.","marker":"[60]"},{"why":"Is the recent simulation including anisotropic scattering that the paper argues overweights anisotropy.","marker":"[68]"},{"why":"Provide observed ro-vibrational term values used for the rotational corrections to the Raman shifts.","marker":"[69, 70]"}],"fun_headline_variants":["CO2 isotopologue Raman spectra reproduced by algebraic model","Rotational correction explains CO2 isotopologue Raman mismatch","Model predicts hot CO2 isotopologue Raman lines accurately","Dynamical-group polarizability fits CO2 isotopologue Raman","Raman discrepancy resolved by rotation in CO2 isotopologues"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The explanation of the experimental discrepancy assumes that the only relevant difference between the earlier 570 K spectrum and the authors' 500 K and 650 K spectra is acquisition time and therefore how much rotational structure is recorded; if temperature, pressure, sample composition, or calibration differ between the measurements, the conclusion that rotational energy corrections rather than anisotropic scattering produce the line shapes would not follow.","fun_headline_variants_meta":{"raw":{"variants":["CO2 isotopologue Raman spectra reproduced by algebraic model","Rotational correction explains CO2 isotopologue Raman mismatch","Model predicts hot CO2 isotopologue Raman lines accurately","Dynamical-group polarizability fits CO2 isotopologue Raman","Raman discrepancy resolved by rotation in CO2 isotopologues"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000268,"raw_usage":{"total_tokens":1662,"prompt_tokens":1034,"completion_tokens":628,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":650,"completion_tokens_details":{"reasoning_tokens":544}},"tokens_in":650,"tokens_out":628,"duration_ms":6349,"temperature":1.0,"reasoning_tokens":544,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:35:44.283173+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A controlled experiment on the same isotopic mixture at a single temperature, with integration time varied from 300 s to 50 minutes and all other conditions fixed, would settle the claim: if longer acquisition does not grow the rotational structure around 1380.46 $\\mathrm{cm}^{-1}$ and the other discrepant line, the rotational explanation fails. A separate check is a polarized Raman measurement that isolates the anisotropic contribution; if the depolarized component contributes more than the few percent the authors assume, the dismissal of anisotropic scattering would be wrong.","supporting_citations":[{"cited_title":"Álvarez, G","cited_arxiv_id":null,"evidence_quote":"Supplies the 42 experimental transition moments of the principal isotopologue used to fit the polarizability derivatives, and the 1780 K benchmark spectrum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental transition moments and the 298 K and 570 K spectra of 13C16O2 and 16O13C18O that the simulations are compared against."},{"cited_title":"Bermudez-Montaña, R","cited_arxiv_id":null,"evidence_quote":"Provides the vibrational term values and wave functions for the principal and symmetric isotopologues used in the simulations."},{"cited_title":"Bermúdez-Montaña, M","cited_arxiv_id":null,"evidence_quote":"Supplies the vibrational energies and wave-function fits for 13C16O2 used in the computed Raman shifts."},{"cited_title":"Bermúdez-Montaña, M","cited_arxiv_id":null,"evidence_quote":"Supplies the vibrational description and term values for the asymmetric isotopologue 16O13C18O."},{"cited_title":"BermúdezMontaña, M","cited_arxiv_id":null,"evidence_quote":"Establishes the earlier algebraic Raman simulation and the polyad choice used here."},{"cited_title":"Suárez, C","cited_arxiv_id":null,"evidence_quote":"Describes the fitting method for polarizability derivatives and the extrapolated polarizability surfaces on which the transfer relies."}],"review_version":1}