{"id":"8b6770d3-c25b-4e79-8643-b30ce2ae6568","arxiv_id":"2505.07361","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"New experimental and simulated valence photoelectron spectra for 2-, 4-, and 2,4-dithiouracil assign the first ionization bands to sulfur and oxygen based orbitals.","lead":"Researchers recorded sharp gas-phase electron spectra of three sulfur-substituted uracil variants and simulated them with coupled-cluster quantum chemistry. The data give the first complete valence photoemission reference for these molecules, helping future ultrafast studies of their photochemistry.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Abstract's 'well reproduced up to 12 eV' is contradicted by the paper's own report of largely overestimated D3 intensity for 2,4-dTU within that range.","rationale":"The reader identified the time-independent double-harmonic adiabatic Hessian approximation as the weakest assumption. That approximation is indeed the root cause of the intensity overestimation, but the most load-bearing concern is the resulting mismatch between the abstract's claim and the paper's own reported results. The authors explicitly acknowledge the failure for higher states and the overestimated D3 intensity, yet the abstract still claims 'well reproduced' for bands up to 12 eV, which includes the problematic D3 band. This is an internal inconsistency, not merely a limitation shared with the field. The concern does not invalidate the experimental data or the qualitative orbital assignment, both of which appear well supported. It does require tempering the abstract and adding quantitative error estimates, so the CONDITIONAL verdict from the reader remains appropriate. Other potential concerns, such as tautomer composition or energy-dependent photoionization cross-sections, are less directly tied to the paper's central claim and are not necessary to resolve for the main conclusions.","tokens_in":14194,"tokens_out":4826,"duration_ms":44547,"concrete_test":"Compute the integrated calculated and experimental spectral intensities for 2,4-dTU in the 9-11 eV binding energy range. If the calculated D3 contribution exceeds the experimental band area by more than a factor of two, the 'well reproduced' claim is quantitatively false. A complementary check is to fit the experimental spectrum allowing D3 intensity to scale freely and report the best-fit scaling factor; if it is significantly below 1, the overestimation is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in the abstract that 'ionization bands up to 12 eV show strong vibrational progressions and are well reproduced by the calculations' is not supported by the authors' own results. In the 'Vibronic transitions' section, they state: 'The D3 intensity for the 2,4-dTU is largely overestimated compared to the experimental spectrum.' For 2,4-dTU, D3 has a vertical ionization potential of 9.137 eV (Table 4), well within the 12 eV range claimed. Thus the blanket 'well reproduced' statement is internally inconsistent. The paper also notes that D4 and D7 intensities are strongly overestimated (D7 for 2,4-dTU), and that optimizations of higher cationic states hit crossing regions, so the double-harmonic adiabatic Hessian approach fails for higher states. The authors themselves conclude that 'this approach fails to reproduce the bands for higher cationic states.' The abstract's unqualified claim overstates the agreement, particularly for 2,4-dTU where the D3 band is a major feature around 9-10 eV. This matters because the abstract is the primary statement of the paper's contribution and currently misrepresents the level of agreement between theory and experiment.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports gas-phase valence photoelectron spectra of 2-thiouracil, 4-thiouracil, and 2,4-dithiouracil measured at 100 eV photon energy, including what the authors state are the first spectra of 4-TU and 2,4-dTU. The experimental bands are assigned with EOM-IP-CCSD vertical ionization energies, Dyson orbitals, and time-independent double-harmonic adiabatic Hessian vibronic simulations. The central claims are that the first ionization potentials lie near 8.4-8.7 eV, about 1 eV below uracil; that the bands up to 12 eV show strong vibrational progressions and are well reproduced by the calculations; and that these bands arise from sulfur- and oxygen-localized orbitals, with nonadiabatic couplings important at higher binding energies.","tokens_in":14367,"tokens_out":3152,"duration_ms":29619,"significance":"If the results stand, the paper provides valuable high-resolution reference spectra for three thionated nucleobases, with consistent experimental-theoretical orbital assignments across the series. The experimental first IPs have small uncertainties, the S/O assignment is corroborated by the systematic absence of the oxygen band in 2,4-dTU, and the data are openly deposited. The EOM-IP-CCSD calculations and Dyson orbital analysis are a strength, as is the explicit identification of where the double-harmonic adiabatic Hessian approximation breaks down. The main weakness is that the abstract overstates the level of agreement with theory, and one numerical claim in the abstract is not consistent with the paper's own Table 1.","major_comments":[{"comment":"The abstract states that the first ionization potential is between 8.4 and 8.7 eV for all three thiouracils, but Table 1 reports 8.74 ± 0.03 eV for 2-TU. This is outside the stated range. Please revise the range or the wording (e.g., '8.4-8.8 eV' or 'around 8.4-8.7 eV') so the abstract matches the reported data.","section":"Abstract vs. Table 1"},{"comment":"The abstract's claim that 'ionization bands up to 12 eV show strong vibrational progressions and are well reproduced by the calculations' is internally inconsistent with the manuscript's own results. In the 'Vibronic transitions' section the authors state that the D3 intensity for 2,4-dTU is largely overestimated compared with experiment, and D3 has a vertical IP of 9.137 eV (Table 4), well inside the 12 eV range. The D7 intensity is also described as strongly overestimated. The abstract should be qualified, for example by restricting the 'well reproduced' statement to the first two bands of 2-TU and 4-TU, and by noting the overestimation in 2,4-dTU.","section":"Abstract vs. §Vibronic transitions and Fig. 2"},{"comment":"The reported agreement is achieved after applying per-molecule energy shifts of 0.05, 0.10, -0.02, and -0.07 eV to the theoretical spectra (Fig. 1 caption), in addition to a hand-chosen Lorentzian width of 0.02 eV. Because the 'well reproduced up to 12 eV' claim depends on these alignment choices, the abstract and the results discussion should explicitly disclose that the theoretical spectra were shifted to match experiment. Without this disclosure, the claim overstates the parameter-free character of the comparison.","section":"Computational Details and Fig. 1 caption"}],"minor_comments":[{"comment":"In the equation for I(E), the linewidth parameter appears to be missing in the rendered text; it should read γ = 0.02 eV rather than ' = 0.02 eV'.","section":"Methods, Lorentzian broadening equation"},{"comment":"The sentence 'This was also observed in the D2 optimizations for 2,4-TU' should read '2,4-dTU' for consistency with the rest of the manuscript.","section":"§Vibronic transitions"},{"comment":"The phrase 'the undertaken geometry optimizations for higher energy cationic excited states usually hit intersection regions' is awkward; consider rephrasing to 'the geometry optimizations of higher cationic states usually reached intersection regions'.","section":"§Vibronic transitions"},{"comment":"The manuscript alternates between '2,4-TU' and '2,4-dTU'; please standardize on '2,4-dTU'.","section":"Throughout"},{"comment":"The conclusion says the approach 'fails to reproduce the bands for higher cationic states'; this is appropriately hedged, but the abstract should be brought in line with this statement.","section":"Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The experimental work is timely and the orbital assignment is convincing, but the abstract contains a numerical inconsistency with Table 1 and an overstatement of theory-experiment agreement that the authors' own results contradict. These are fixable by revised wording, so major revision rather than rejection seems appropriate. I would also encourage the editor to ask the authors to state the applied energy shifts in the main text and not only in the figure caption."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the useful part. The paper reports the first gas-phase valence photoelectron spectra of 4-thiouracil and 2,4-dithiouracil, measured at 30 meV resolution, and gives the complete thiouracil series. The first ionization energies are credible, with uncertainties under 0.03 eV, and the assignment of the low bands to sulfur- versus oxygen-localized orbitals is convincing. The identification of a C=O stretch progression in the second band of 2-TU and 4-TU is a concrete result. The calculations are at a respectable level: EOM-IP-CCSD, Dyson orbital norms, and a double-harmonic Franck-Condon simulation with explicit mode assignments. Derived data are on Zenodo.\n\nThe soft spot is the abstract. Claiming bands up to 12 eV are 'well reproduced' is not supported by the authors' own results. They state that the D3 intensity for 2,4-dTU is largely overestimated and D7 is strongly overestimated; D3 lies at roughly 9.1 eV vertical, inside the 12 eV range. The conclusion's more cautious phrase 'good match for the first two ionization bands' is accurate. The per-molecule energy shifts (0.02–0.10 eV) and the Lorentzian width are reasonable, but the paper does not discuss the sensitivity of the comparison to them, and there are no error bars on the theoretical curves. The double-harmonic approximation is explicitly said to fail for higher cationic states, so the abstract's scope should be limited to the first two bands per molecule.\n\nMinor: tautomers are not discussed. At 120 °C, enol/thiol forms are conceivable; prior 2-TU work makes the thione assumption plausible, but a sentence would help.\n\nOverall, this is a solid reference paper for people doing time-resolved photoelectron spectroscopy on thionucleobases. The new spectra are the contribution; the theory supports the low-band assignment but not the higher-energy structure. The abstract overreach is fixable in revision. I'd send it to peer review and ask for a revised abstract and a short discussion of the alignment shifts.","headline":"New reference spectra for two thiouracils with a sound S/O assignment, but the abstract overstates theory–experiment agreement up to 12 eV.","tokens_in":15013,"tokens_out":3157,"would_cite":true,"duration_ms":28008,"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":"Valence photoelectron spectra of all three thiouracils show a first ionization potential of 8.4-8.7 eV, about 1 eV below uracil, with bands up to 12 eV reproduced by EOM-IP-CCSD calculations.","keywords":["valence photoelectron spectroscopy","thiouracils","ionization potential","vibrational progression","EOM-IP-CCSD","Dyson orbitals","adiabatic Hessian approximation","nucleobases"],"falsifier":"Measure the high-resolution valence photoelectron spectrum of 2,4-dithiouracil in the 10-11 eV binding-energy range: the paper's assignment predicts no oxygen-localized band there, so seeing a structured band would contradict the central orbital assignment. Alternatively, comparing the measured vibrational spacing in 4-thiouracil's second band with the calculated ~1836 cm-1 mode would test whether the harmonic model holds.","tokens_in":13938,"feed_emoji":"🧪","tokens_out":9921,"duration_ms":85384,"temperature":0.7,"pith_summary":"Thiouracils are thionated nucleobases of medical interest, but their gas-phase valence photoelectron reference spectra were incomplete. This paper fills that gap for the full series—2-thiouracil, 4-thiouracil, and 2,4-dithiouracil—and establishes that replacing uracil's oxygen with sulfur lowers the first ionization potential to between 8.4 and 8.7 eV, about 1 eV below uracil. The vibrationally resolved bands up to 12 eV are reproduced by equation-of-motion coupled-cluster calculations combined with the double-harmonic adiabatic Hessian approximation, allowing each band to be assigned to ionization of sulfur- or oxygen-localized orbitals. The paper also shows that the higher-energy bands, above about 12 eV, cannot be captured by this single-state harmonic picture because the cationic states cross.","feed_headline":"Thiouracils ionize at 8.4-8.7 eV, about 1 eV below uracil","feed_subtitle":"First gas-phase valence spectra of 4-thiouracil and 2,4-dithiouracil pin down sulfur- and oxygen-band assignments.","key_machinery":"The argument is carried by pairing EOM-IP-CCSD (equation-of-motion coupled cluster for ionization potentials) with the time-independent double-harmonic adiabatic Hessian approach. EOM-IP-CCSD computes the ionization energies and cationic states directly, while the adiabatic Hessian approach computes Franck-Condon factors from the harmonic normal modes of the neutral and cation states, including Duschinsky rotations and temperature effects. The resulting stick spectra are broadened and compared with the experimental spectra, and Dyson orbital norms estimate relative band intensities. Normal-mode analysis then identifies the specific vibrations that generate the long progressions, such as mode 26, an NH-bending plus C=O-stretching plus C-C-stretching mode near 1800-2100 cm-1, which dominates the oxygen-ionization band in 2- and 4-thiouracil.","core_discovery":"The paper's central claim is that every thiouracil ionizes at 8.4-8.7 eV—measured as 8.74 ± 0.03 eV for 2-thiouracil, 8.42 ± 0.01 eV for 4-thiouracil, and 8.45 ± 0.01 eV for 2,4-dithiouracil—roughly 1 eV below uracil, and that this lowering is a general consequence of thionation. The spectra of 4-thiouracil and 2,4-dithiouracil are reported for the first time. EOM-IP-CCSD/cc-pVTZ calculations reproduce the vibrational progressions of the first two bands and show that the lowest band comes from ionizing sulfur-dominated orbitals, while the band near 10-11 eV in the two singly thionated molecules comes from oxygen-dominated orbitals; 2,4-dithiouracil lacks this band because it has no oxygen. The same calculations show that Koopmans' orbital-energy estimates overshoot the ionization energies by up to about 2.3 eV and that the Dyson orbitals for the low-lying transitions are close to canonical Hartree-Fock orbitals, with norms above 0.8. For binding energies above about 12 eV, the paper argues that nonadiabatic couplings between cationic states must be included, since geometry optimizations of higher states repeatedly land in crossing regions.","pith_inferences":["The measured ~1 eV drop in ionization energy, which the paper does not discuss in redox terms, implies thiouracils are markedly easier to oxidize than uracil, a property that could matter in charge-transfer contexts.","The oxygen-localized band near 10-11 eV could serve as a spectroscopic fingerprint for distinguishing singly from doubly thionated uracils, an application the paper does not propose.","A natural next step, not taken here, would be to build a two-state vibronic-coupling model for the crossing regions and test whether it reproduces the broad 12-16 eV envelope in the experimental spectra."],"forward_implications":["The first ionization potentials of all three thiouracils fall between 8.4 and 8.7 eV, so thionation lowers uracil's ionization energy by roughly 1 eV regardless of which oxygen is replaced.","The structured band near 10-11 eV in 2- and 4-thiouracil comes from ionizing an oxygen-localized orbital, and 2,4-dithiouracil shows no such band, confirming the assignment.","For binding energies up to about 12 eV, EOM-IP-CCSD with the double-harmonic adiabatic Hessian approximation reproduces the measured vibrational progressions, enabling mode assignments such as the C=O-stretch-dominated progression in the oxygen band.","Simple Koopmans' theorem estimates overestimate the measured ionization energies by up to about 2.3 eV, so a correlated treatment is necessary even for assigning the lowest bands.","Above about 12 eV, the single-state harmonic approximation breaks down because cationic states cross, so interpreting those bands requires a nonadiabatic or multistate treatment."],"supporting_citations":[{"why":"Previous simulated photoelectron spectra of the three thiouracils; the paper compares its orbital assignments and Dyson norms against them.","marker":"[18]"},{"why":"Defines the equation-of-motion coupled-cluster method for ionization potentials used to compute cation states and energies.","marker":"[22,23]"},{"why":"Supplies the time-independent double-harmonic adiabatic Hessian method used for the vibrationally resolved spectra.","marker":"[32]"},{"why":"Provides Dyson orbital norms used to estimate transition intensities and to justify the one-electron picture.","marker":"[33]"},{"why":"Provides the high-resolution uracil valence photoelectron spectrum used as the reference for the thionation shift.","marker":"[34]"},{"why":"Earlier measurements of the 2-thiouracil first ionization energy with which the new value agrees.","marker":"[10,35]"},{"why":"Earlier uracil ionization energies used to quantify the ~1 eV shift upon thionation.","marker":"[36,37]"}],"fun_headline_variants":["Thiouracils ionize ~1 eV below uracil","Sulfur substitution lowers thiouracil ionization by 1 eV","First thiouracil spectra reveal 1 eV ionization drop","Thiouracil ionization: ~1 eV shift from uracil"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculations assume that the low-lying cation states responsible for the bands below 12 eV are each a single harmonic well, well separated from other electronic states; if that assumption is false, the calculated vibrational progressions and band assignments lose their footing.","fun_headline_variants_meta":{"raw":{"variants":["Thiouracils ionize ~1 eV below uracil","Sulfur substitution lowers thiouracil ionization by 1 eV","First thiouracil spectra reveal 1 eV ionization drop","Thiouracil ionization: ~1 eV shift from uracil"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000743,"raw_usage":{"total_tokens":3361,"prompt_tokens":1036,"completion_tokens":2325,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":652,"completion_tokens_details":{"reasoning_tokens":2255}},"tokens_in":652,"tokens_out":2325,"duration_ms":16110,"temperature":1.0,"reasoning_tokens":2255,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:18:22.023190+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the high-resolution valence photoelectron spectrum of 2,4-dithiouracil in the 10-11 eV binding-energy range: the paper's assignment predicts no oxygen-localized band there, so seeing a structured band would contradict the central orbital assignment. Alternatively, comparing the measured vibrational spacing in 4-thiouracil's second band with the calculated ~1836 cm-1 mode would test whether the harmonic model holds.","supporting_citations":[],"review_version":1}