arxiv: 2605.10311 · v2 · submitted 2026-05-11 · ⚛️ physics.chem-ph

Recognition: 2 theorem links

· Lean Theorem

Do Water Molecules Always Stabilize Resonances? Microhydration Effects on Thymine Shape Resonances

Sujan Mandal , Jishnu Narayanan S J , Ankita Gogoi , Madhubani Mukherjee , Idan Haritan , Achintya Kumar Dutta

Authors on Pith no claims yet

Pith reviewed 2026-05-14 21:40 UTC · model grok-4.3

classification ⚛️ physics.chem-ph

keywords thymineshape resonancesmicrohydrationresonance lifetimesnucleobaseselectron attachmentPadé approachDLPNO-EA-EOM-CCSD

0 comments

The pith

Water molecules stabilize thymine's lowest shape resonance, tripling its lifetime from 39 fs to 110 fs in a trihydrated cluster.

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

The paper examines microhydration effects on the three low-lying π* shape resonances of thymine through high-level electronic structure calculations. It finds that the first two resonances shift downward in energy with markedly longer lifetimes as water molecules are added, while the third resonance displays more variable behavior tied to specific geometries. The observed changes result from a competition among hydrogen bonding, electrostatic interactions, microsolvation-induced structural distortions, and contributions from diffuse basis functions on the water molecules. Ghost-atom tests separate genuine physical stabilization from basis-set artifacts. These results indicate that resonance properties in nucleobases are sensitive to the precise local solvation environment rather than following a uniform stabilization pattern.

Core claim

Calculations on isolated thymine and its microhydrated clusters show that the 1π* and 2π* resonances undergo systematic stabilization with significant lifetime increases upon hydration, whereas the 3π* resonance behaves more complexly; the lifetime of the lowest resonance specifically rises from 39 fs to 110 fs in the thymine(H2O)3 cluster. Ghost-atom calculations demonstrate that diffuse basis functions centered on nearby water molecules contribute to apparent stabilization, while explicit inclusion of water produces genuine physical effects. Resonance positions and widths further depend strongly on the local hydrogen-bonding arrangement across different monohydrated conformers.

What carries the argument

Resonance via Padé approach combined with DLPNO-EA-EOM-CCSD calculations performed on thymine and multiple microhydrated clusters, together with ghost-atom basis tests and geometry-specific analysis.

If this is right

Resonance lifetimes in microhydrated nucleobases can increase substantially compared with the gas-phase values.
The local hydrogen-bonding geometry around thymine strongly controls both resonance energy and lifetime.
Diffuse basis functions supplied by nearby water molecules contribute measurably to apparent resonance shifts.
Explicit inclusion of water molecules produces genuine physical stabilization beyond pure basis-set effects.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Longer-lived resonances in hydrated DNA bases could raise the probability of electron-driven damage pathways in aqueous biological settings.
Conformer-dependent resonance behavior implies that different hydration sites on biomolecules may exhibit distinct attachment probabilities.
Parallel microhydration calculations on other nucleobases could reveal whether the observed lifetime extension is a general feature of hydrated DNA constituents.

Load-bearing premise

The DLPNO-EA-EOM-CCSD method paired with the Resonance via Padé approach gives accurate resonance positions and widths for both the isolated molecule and the hydrated clusters without large errors from approximations or basis-set limitations.

What would settle it

An experimental determination of the resonance lifetime for thymine microhydrated with three water molecules that yields a value substantially different from 110 fs would directly test the reported stabilization.

read the original abstract

We investigate microhydration effects on the three low-lying {\pi}* shape resonances of thymine using the Resonance via Pad\'e approach in combination with the DLPNO-EA-EOM-CCSD method. For isolated thymine, the calculated resonance positions are benchmarked against projected CAP-EA-EOM-CCSD calculations and compared with available theoretical and experimental data. Upon hydration, the 1{\pi}* and 2{\pi}* resonances undergo systematic stabilization accompanied by significant increases in their lifetimes, whereas the 3{\pi}* resonance exhibits a more complex behavior. In particular, the lifetime of the lowest resonance increases from 39 fs in isolated thymine to 110 fs in the thymine(H2O)3 cluster. Detailed analysis reveals that the observed resonance shifts arise from competing contributions involving hydrogen bonding, electrostatic interactions, microsolvation-induced geometric distortion, and finite-basis-set effects. Ghost-atom calculations demonstrate that diffuse basis functions associated with nearby water molecules contribute appreciably to the apparent stabilization, while explicit inclusion of water molecules leads to genuine physical stabilization of the resonance states. Furthermore, calculations on multiple conformers of the monohydrated cluster show that resonance positions and lifetimes depend strongly on the local hydrogen-bonding arrangement and microsolvation geometry. These findings demonstrate that resonance stabilization in microhydrated nucleobases is governed by a subtle interplay between geometry, basis-set effects, and intermolecular interactions.

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.

Desk Editor's Note private letter to a colleague

Microhydration stabilizes thymine's lowest pi* resonance and lengthens its lifetime from 39 fs to 110 fs, but the hydrated widths rest on unbenchmarked DLPNO-Padé numbers.

read the letter

The paper's central result is that adding three water molecules to thymine stabilizes the lowest shape resonance and increases its lifetime from 39 fs to 110 fs. The authors reach this by running DLPNO-EA-EOM-CCSD calculations, extracting widths with the Padé approach, and then dissecting the shift into hydrogen-bonding, electrostatic, geometric, and basis-set pieces using ghost-atom tests and multiple monohydrate conformers. That decomposition is the genuinely new part; earlier work on isolated thymine did not quantify how microsolvation and nearby diffuse functions compete. The isolated-molecule benchmarks against projected CAP-EA-EOM-CCSD and literature values are clean and give the method credibility. The conformer scan is also useful because it shows the resonance position and width are sensitive to the exact hydrogen-bond geometry. The soft spot is exactly the one the stress-test flagged: there is no equivalent CAP or full-method benchmark for the hydrated cluster widths. The DLPNO locality approximation and the Padé fit are only validated on the isolated case, and the resonance orbital becomes more diffuse once waters are present. That leaves the 110 fs number with an unquantified uncertainty, probably in the 20-30% range. The trend is still believable, but a referee would reasonably ask for at least one cross-check on the trihydrate. This work is for computational chemists who model electron attachment to nucleobases or who need concrete numbers on how small water clusters affect resonance lifetimes. A reader building radiation-chemistry or DNA-damage models will get usable data and a clear breakdown of the contributing effects. It deserves peer review because the isolated benchmarks are solid, the new analysis is reproducible from the methods described, and the central claim is narrow enough that targeted additional calculations can settle the remaining uncertainty.

Referee Report

2 major / 2 minor

Summary. The manuscript investigates microhydration effects on the three low-lying π* shape resonances of thymine using the Resonance via Padé approach combined with DLPNO-EA-EOM-CCSD. For isolated thymine, resonance positions are benchmarked against projected CAP-EA-EOM-CCSD calculations and literature data. The study reports that hydration leads to stabilization of the 1π* and 2π* resonances with increased lifetimes, while the 3π* shows more complex behavior; specifically, the lowest resonance lifetime increases from 39 fs (isolated) to 110 fs in the thymine(H2O)3 cluster. Stabilization arises from competing hydrogen bonding, electrostatic interactions, microsolvation-induced geometric distortions, and basis-set effects, as analyzed via ghost-atom calculations and studies of multiple monohydrated conformers.

Significance. If the central quantitative results hold, the work provides valuable insights into how explicit water molecules modulate temporary anion resonances in nucleobases, with implications for understanding radiation-induced DNA damage in aqueous environments. The decomposition of stabilization mechanisms into physical versus basis-set contributions, the geometry dependence shown across conformers, and the benchmarking of the isolated case against CAP-EA-EOM-CCSD are notable strengths. The application of DLPNO-EA-EOM-CCSD enables treatment of larger microhydrated systems, offering a practical route to falsifiable predictions for experiment.

major comments (2)

[microhydration results (thymine(H2O)3 cluster)] The headline quantitative claim (lifetime of the lowest resonance increasing from 39 fs to 110 fs in thymine(H2O)3) is load-bearing for the central conclusion but rests on the unverified accuracy of DLPNO-EA-EOM-CCSD widths (via Padé) for the hydrated cluster. While isolated thymine positions are benchmarked against projected CAP-EA-EOM-CCSD, no equivalent full-method comparison, convergence data, or error quantification is supplied for the microhydrated systems where the resonance orbital is more diffuse.
[abstract and results sections] The abstract and results report specific lifetime and position values without accompanying error bars, basis-set extrapolation details, or sensitivity analysis of the Padé approximant for the hydrated cases; this leaves the ~20-30% accuracy needed for the lifetime change claim unestablished.

minor comments (2)

[methods] Clarify in the methods whether the same basis set and DLPNO thresholds were used uniformly for isolated and hydrated systems, and add a brief note on any observed sensitivity to these choices.
[figures] Ensure figure captions for resonance spectra or lifetime plots explicitly reference the Padé fitting procedure and any convergence checks performed.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful reading of our manuscript and for recognizing the potential value of our findings on microhydration effects in nucleobase resonances. We address each major comment point by point below, providing clarifications on our methodological choices and outlining specific revisions to strengthen the presentation of uncertainties and validation.

read point-by-point responses

Referee: [microhydration results (thymine(H2O)3 cluster)] The headline quantitative claim (lifetime of the lowest resonance increasing from 39 fs to 110 fs in thymine(H2O)3) is load-bearing for the central conclusion but rests on the unverified accuracy of DLPNO-EA-EOM-CCSD widths (via Padé) for the hydrated cluster. While isolated thymine positions are benchmarked against projected CAP-EA-EOM-CCSD, no equivalent full-method comparison, convergence data, or error quantification is supplied for the microhydrated systems where the resonance orbital is more diffuse.

Authors: We agree that a direct comparison of DLPNO-EA-EOM-CCSD against projected CAP-EA-EOM-CCSD for the thymine(H2O)3 cluster would provide the strongest possible validation of the resonance widths. Such calculations remain computationally prohibitive at present due to the scaling of the CAP implementation and the increased system size. The DLPNO approximation is designed to recover canonical EA-EOM-CCSD results to high accuracy for valence and diffuse states when appropriate thresholds are used, and our isolated-thymine benchmarks already demonstrate close agreement for both positions and widths. In the revised manuscript we will add explicit convergence tests with respect to DLPNO pair-natural-orbital thresholds for the hydrated cluster, together with a Padé-approximant sensitivity analysis that quantifies the uncertainty in the extracted widths. We will also expand the discussion of expected error propagation from the isolated benchmarks to the microhydrated case. These additions will be placed in the results section and will not alter the central conclusions. revision: partial
Referee: [abstract and results sections] The abstract and results report specific lifetime and position values without accompanying error bars, basis-set extrapolation details, or sensitivity analysis of the Padé approximant for the hydrated cases; this leaves the ~20-30% accuracy needed for the lifetime change claim unestablished.

Authors: We accept that the absence of explicit uncertainty estimates weakens the quantitative presentation. In the revision we will supply error bars on all reported resonance positions and lifetimes for both isolated and microhydrated systems, derived from the variation of the Padé approximant across different fitting windows and from basis-set sensitivity tests performed on the monohydrated conformers. A short paragraph discussing the magnitude of these uncertainties relative to the observed lifetime increase (39 fs to 110 fs) will be added to the results section, and the abstract will be updated to indicate that the reported changes exceed the estimated methodological uncertainties. Basis-set extrapolation details already present for the isolated molecule will be extended to the hydrated clusters where additional calculations are feasible. revision: yes

standing simulated objections not resolved

Direct full-method benchmark of DLPNO-EA-EOM-CCSD widths against projected CAP-EA-EOM-CCSD for the thymine(H2O)3 cluster, which exceeds current computational resources.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper's results are obtained by direct application of the DLPNO-EA-EOM-CCSD method with the Resonance via Padé approach to compute resonance positions and widths for both isolated and microhydrated thymine. Isolated thymine positions are benchmarked against independent projected CAP-EA-EOM-CCSD calculations and external literature data. Microhydration effects, including the reported lifetime increase, follow from explicit calculations on the clusters with ghost-atom tests to isolate basis contributions; no parameters are fitted to the target lifetimes or positions, and no equations reduce the outputs to the inputs by construction. The chain is self-contained against external benchmarks with no self-definitional, fitted-prediction, or load-bearing self-citation circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard quantum chemistry approximations for electron correlation and resonance extraction with no new free parameters or invented entities; results depend on the validity of DLPNO-EA-EOM-CCSD for temporary anion states.

axioms (2)

domain assumption The DLPNO approximation in EA-EOM-CCSD accurately captures electron attachment energies and resonance parameters for nucleobases.
Core method invoked without additional justification in the abstract.
domain assumption The Resonance via Padé approach correctly extracts resonance positions and widths from the computed complex energies.
Central technique for obtaining the reported positions and lifetimes.

pith-pipeline@v0.9.0 · 5585 in / 1528 out tokens · 41126 ms · 2026-05-14T21:40:59.867063+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We investigate microhydration effects on the three low-lying π* shape resonances of thymine using the Resonance via Padé approach in combination with the DLPNO-EA-EOM-CCSD method.
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

the lifetime of the lowest resonance increases from 39 fs in isolated thymine to 110 fs in the thymine(H2O)3 cluster

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages

[1]

Introduction High-energy ionizing radiation (such as X-rays or γ-rays) is well known to damage biomolecules such as DNA and RNA. In recent decades, substantial experimental and theoretical efforts have been devoted to understanding the mechanisms of radiation-biomolecule interactions.1–7 Radiation damage to DNA occurs through multiple pathways that are br...

work page
[2]

Computational Details Theoretical simulation of anionic resonances is more challenging than that of bound states because resonances are embedded in the continuum .48 As a result, the corresponding wave functions are not square-integrable and cannot be represented as discrete eigenstates within the conventional Hermitian framework. Non-Hermitian quantum me...

work page
[3]

Results and Discussion 3.1 Isolated thymine Theoretical modeling of temporarily bound anions is challenging due to their strong sensitivity to both the basis set and the level of electronic structure theory. The reliability of the RVP-EA-EOM- DLPNO-CCSD approach ( hereafter referred to as RVP) has been extensively examined in our previous studies.63,65,94...

work page
[4]

The surrounding water molecules do not alter the qualitative nature of the resonance states but modify their positions and lifetim es

Conclusion We have studied the effect of solvation on the shape resonance states of thymine, taking microhydrated thymine as a test case. The surrounding water molecules do not alter the qualitative nature of the resonance states but modify their positions and lifetim es. The addition of a single water molecule does not necessarily guarantee stabilization...

work page doi:10.1016/s0167-6881(01)80023-9 2022
[5]

(13) Huels, M

https://doi.org/10.1126/science.287.5458.1658. (13) Huels, M. A.; Boudaïffa, B.; Cloutier, P.; Hunting, D.; Sanche, L. Single, Double, and Multiple Double Strand Breaks Induced in DNA by 3−100 eV Electrons. J. Am. Chem. Soc. 2003, 125 (15), 4467–4477. https://doi.org/10.1021/ja029527x. (14) Abdoul-Carime, H.; Gohlke, S.; Fischbach, E.; Scheike, J.; Illenb...

work page doi:10.1126/science.287.5458.1658 2003
[6]

(27) Arumainayagam, C

https://doi.org/10.1038/s41467-018-08005-z. (27) Arumainayagam, C. R.; Lee, H.-L.; Nelson, R. B.; Haines, D. R.; Gunawardane, R. P. Low- Energy Electron -Induced Reactions in Condensed Matter. Surf. Sci. Rep. 2010, 65, 1 –44. https://doi.org/10.1016/j.surfrep.2009.09.001. (28) Kumari, B.; Huwaidi, A.; Robert, G.; Cloutier, P.; Bass, A. D.; Sanche, L.; Wag...

work page doi:10.1038/s41467-018-08005-z 2010
[7]

(41) Davis, D.; Sajeev, Y

https://doi.org/10.1021/acs.jpca.5c04881. (41) Davis, D.; Sajeev, Y. A Hitherto Unknown Stability of DNA Basepairs. Chem. Commun. 2020, 56 (93), 14625–14628. https://doi.org/10.1039/D0CC06641A. (42) Anstöter, C. S.; DelloStritto, M.; Klein, M. L.; Matsika, S. Modeling the Ultrafast Electron Attachment Dynamics of Solvated Uracil. J. Phys. Chem. A 2021, 12...

work page doi:10.1021/acs.jpca.5c04881 2020
[8]

(57) Fennimore, M

https://doi.org/10.1021/acs.jpclett.6b01601. (57) Fennimore, M. A.; Matsika, S. Core-Excited and Shape Resonances of Uracil. Phys. Chem. Chem. Phys. 2016, 18 (44), 30536–30545. https://doi.org/10.1039/C6CP05342D. (58) Winstead, C.; McKoy, V. Low -Energy Electron Collisions with Gas -Phase Uracil. J. Chem. Phys. 2006, 125 (17), 174304. https://doi.org/10.1...

work page doi:10.1021/acs.jpclett.6b01601 2016
[9]

(90) Nooijen, M.; Bartlett, R

https://doi.org/10.1063/1.1378323. (90) Nooijen, M.; Bartlett, R. J. Equation of Motion Coupled Cluster Method for Electron Attachment. J. Chem. Phys. 1995, 102 (9), 3629–3647. https://doi.org/10.1063/1.468592. (91) Dutta, A. K.; Saitow, M.; Demoulin, B.; Neese, F.; Izsák, R. A Domain -Based Local Pair Natural Orbital Implementation of the Equation of Mot...

work page doi:10.1063/1.1378323 1995
[10]

(96) Neese, F

https://doi.org/10.1039/C9CP06869D. (96) Neese, F. Software Update: The ORCA Program System —Version 5.0. WIREs Comput. Mol. Sci. 2022, 12 (5), e1606. https://doi.org/10.1002/wcms.1606. (97) Automatic-RVP: GitHub repository. https://github.com/haritan/RVP (accessed 202 6-05- 10). (98) Zuev, D.; Jagau, T.-C.; Bravaya, K. B.; Epifanovsky, E.; Shao, Y.; Sund...

work page doi:10.1039/c9cp06869d 2022
[11]

(102) Arora, S.; Narayanan S J, J.; Haritan, I.; Adhikary, A.; Dutta, A

https://doi.org/10.1007/s12039-025-02451-1. (102) Arora, S.; Narayanan S J, J.; Haritan, I.; Adhikary, A.; Dutta, A. K. Effect of Protein Environment on the Shape Resonances of RNA Pyrimidine Nucleobases: Insights from a Model System. J. Chem. Phys. 2025, 163 (13), 134103. https://doi.org/10.1063/5.0288514. (103) Aflatooni, K.; Gallup, G. A.; Burrow, P. D...

work page doi:10.1007/s12039-025-02451-1 2025