Efficient EOM-CC-based Protocol for the Calculation of Electron Affinity of Solvated Nucleobases: Uracil as a Case Study

Electron Attachment to the Nucleobase Uracil in Diethylene Glycol: The Signature of a Doorway

The cellular environment plays a significant role in low energy electron-mediated radiation damage to genetic materials. In this study, we have modeled the effect of the bulk medium on electron attachment to nucleobases in diethylene glycol (DEG) using uracil as a test case, in accordance with recent experimental work on the observation of dissociative quasi-free electron attachment to nucleoside via excited anion radical in solution (in DEG). Our EOM-CCSD-based quantum mechanical/molecular mechanical (QM/MM) simulations indicate that the electron scavenging by uracil in DEG is much slower than that observed in the aqueous medium due to its viscosity. This work also establishes that a doorway mechanism exists in uracil microsolvated and bulk solvated with DEG, with the dipole-bound state and solvent-bound state acting as doorway states, respectively.

Toward a numerically efficient description of bulk-solvated anionic states

We investigate the vertical electron attachment energy (VAE) of 1-methyl-4-nitroimidazole, a model radiosensitizer, employing quantum mechanics/molecular mechanics (QM/MM) and QM/polarized continuum (QM/PCM) solvation models. We considered the solvent-excluded surface (QM/PCM-SES) and Van der Waals (QM/PCM-VDW) cavities within the PCM framework, the electrostatic embedding QM/MM (EE-QM/MM) model, and the self-consistent sequential QM/MM polarizable electrostatic embedding (scPEE-S-QM/MM) model. Due to slow VAE convergence concerning the number of QM solvent molecules, full QM calculations prove inefficient. Ensemble averages in these calculations do not align with VAEs computed for the representative solute-solvent configuration. QM/MM and QM/PCM calculations show agreement with each other for sufficiently large QM regions, although the QM/PCM-VDW model exhibits artifacts linked to the cavity. QM/MM models demonstrate good agreement between ensemble averages and VAEs calculated with the representative configuration. Notably, the VAE computed with the scPEE-S-QM/MM model achieves faster convergence concerning the number of QM water molecules compared to the EE-QM/MM model, attributed to enhanced efficiency from MM charge polarization in the scPEE-S-QM/MM approach. This emphasizes the importance of QM/classical models with accurate solute-solvent and solvent-solvent mutual polarization for obtaining converged VAEs at a reasonable computational cost. The full-QM approach is very inefficient, while the microsolvation model is inaccurate. Computational savings in QM/MM models result from electrostatic embedding and the representative configuration, with the scPEE-S-QM/MM approach emerging as an efficient tool for describing bulk-solvated anions within the QM/MM framework. Its potential extends to improving transient anion state descriptions in biomolecules and radiosensitizers, especially given the frequent employment of microsolvation models.

Detangling Solvatochromic Effects by the Effective Fragment Potential Method

Understanding molecular interactions in complex systems opens avenues for the efficient design of new materials with target properties. Energy decomposition methods provide a means to obtain a detailed picture of intermolecular interactions. This work introduces a molecular modeling approach for decomposing the solvatochromic shifts of the electronic excited states into the contributions of the individual molecular fragments of the environment surrounding the chromophore. The developed approach is implemented for the QM/EFP (quantum mechanics/effective fragment potential) model that provides a rigorous first-principles-based description of the electronic states of the chromophores in complex polarizable environments. On the example of two model systems, water pentamer and hydrated uracil, we show how the decomposition of the solvatochromic shifts into the contributions of individual solvent water molecules provides a detailed picture of the intermolecular interactions in the ground and excited states of these systems. The analysis also demonstrates the nonadditivity of solute-solvent interactions and the significant contribution of solute polarization to the total values of solvatochromic shifts.

Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level

The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.

Secondary Electron Attachment-Induced Radiation Damage to Genetic Materials

Reactions of radiation-produced secondary electrons (SEs) with biomacromolecules (e.g., DNA) are considered one of the primary causes of radiation-induced cell death. In this Review, we summarize the latest developments in the modeling of SE attachment-induced radiation damage. The initial attachment of electrons to genetic materials has traditionally been attributed to the temporary bound or resonance states. Recent studies have, however, indicated an alternative possibility with two steps. First, the dipole-bound states act as a doorway for electron capture. Subsequently, the electron gets transferred to the valence-bound state, in which the electron is localized on the nucleobase. The transfer from the dipole-bound to valence-bound state happens through a mixing of electronic and nuclear degrees of freedom. In the presence of aqueous media, the water-bound states act as the doorway state, which is similar to that of the presolvated electron. Electron transfer from the initial doorway state to the nucleobase-bound state in the presence of bulk aqueous media happens on an ultrafast time scale, and it can account for the decrease in DNA strand breaks in aqueous environments. Analyses of the theoretically obtained results along with experimental data have also been discussed.

Electron Attachment to DNA: The Protective Role of Amino Acids

We have studied the effect of amino acids on the electron attachment properties of a DNA nucleobase, with cytosine as a model system. The equation of motion coupled cluster theory with an extended basis set has been used to simulate the electron-attached state of the DNA model system. Arginine, alanine, lysine, and glycine are the four amino acids considered to investigate their role in electron attachment to a DNA nucleobase. The electron attachment to cytosine in all the four cytosine-amino acid gas-phase dimer complexes follows a doorway mechanism, where the electron gets transferred from the initial dipole-bound doorway state to the final nucleobase-bound state through the mixing of electronic and nuclear degrees of freedom. When cytosine is bulk-solvated with glycine, the glycine-bound state acts as the doorway state, where the initial electron density is localized on the bulk amino acid and away from the nucleobase, thus leading to the physical shielding of the nucleobase from the incoming electron. At the same time, the presence of amino acids can increase the stability of the nucleobase-bound anionic state, which can suppress the sugar-phosphate bond rupture caused by dissociative electron attachment to DNA.

Recent Advances in the Study of Negative-Ion Resonances Using Multiconfigurational Propagator and a Complex Absorbing Potential

The transient resonances are a challenge to bound state quantum mechanics. These states lie in the continuum part of the spectrum of the Hamiltonian. For this, one has to treat a continuum problem due to electron-molecule scattering and the many-electron correlation problem simultaneously. Moreover, the description of a resonance requires a wavefunction that bridges the part that resembles a bound state with another that resembles a continuum state such that the continuity of the wavefunction and its first derivative with respect to the distance between the incoming projectile and the target is maintained. A review of the recent advances in the theoretical investigation of the negative-ion resonances (NIR) is presented. The NIRs are ubiquitous in nature. They result from the scattering of electrons off of an atomic or molecular target. They are important for numerous chemical processes in upper atmosphere, space and even biological systems. A contextual background of the existing theoretical methods as well as the newly-developed multiconfigurational propagator tools based on a complex absorbing potential are discussed.

Low-Energy Shape Resonances of a Nucleobase in Water

When high-energy radiation passes through aqueous material, low-energy electrons are produced which cause DNA damage. Electronic states of anionic nucleobases have been suggested as an entrance channel to capture the electron. However, identifying these electronic resonances have been restricted to gas-phase electron-nucleobase studies and offer limited insight into the resonances available within the aqueous environment of DNA. Here, resonance and detachment energies of the micro-hydrated uracil pyrimidine nucleobase anion are determined by two-dimensional photoelectron spectroscopy and are shown to extrapolate linearly with cluster size. This extrapolation allows the corresponding resonance and detachment energies to be determined for uracil in aqueous solution as well as the reorganization energy associated with electron capture. Two shape resonances are clearly identified that can capture low-energy electrons and subsequently form the radical anion by solvent stabilization and internal conversion to the ground electronic state. The resonances and their dynamics probed here are the nucleobase-centered doorway states for low-energy electron capture and damage in DNA.

Computation of Oxidation Potentials of Solvated Nucleobases by Static and Dynamic Multilayer Approaches

The determination of the redox properties of nucleobases is of paramount importance to get insight into the charge-transfer processes in which they are involved, such as those occurring in DNA-inspired biosensors. Although many theoretical and experimental studies have been conducted, the value of the one-electron oxidation potentials of nucleobases is not well-defined. Moreover, the most appropriate theoretical protocol to model the redox properties has not been established yet. In this work, we have implemented and evaluated different static and dynamic approaches to compute the one-electron oxidation potentials of solvated nucleobases. In the static framework, two thermodynamic cycles have been tested to assess their accuracy against the direct determination of oxidation potentials from the adiabatic ionization energies. Then, the introduction of vibrational sampling, the effect of implicit and explicit solvation models, and the application of the Marcus theory have been analyzed through dynamic methods. The results revealed that the static direct determination provides more accurate results than thermodynamic cycles. Moreover, the effect of sampling has not shown to be relevant, and the results are improved within the dynamic framework when the Marcus theory is applied, especially in explicit solvent, with respect to the direct approach. Finally, the presence of different tautomers in water does not affect significantly the one-electron oxidation potentials.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Probing Interfacial Effects on Ionization Energies: The Surprising Banality of Anion-Water Hydrogen Bonding at the Air/Water Interface

Liquid microjet photoelectron spectroscopy is an increasingly common technique to measure vertical ionization energies (VIEs) of aqueous solutes, but the interpretation of these experiments is subject to questions regarding sensitivity to bulk versus interfacial solvation environments. We have computed aqueous-phase VIEs for a set of inorganic anions, using a combination of molecular dynamics simulations and electronic structure calculations, with results that are in excellent agreement with experiment regardless of whether the simulation data are restricted to ions at the air/water interface or to those in bulk aqueous solution. Although the computed VIEs are sensitive to ion-water hydrogen bonding, we find that the short-range solvation structure is sufficiently similar in both environments that it proves impossible to discriminate between the two on the basis of the VIE, a conclusion that has important implications for the interpretation of liquid-phase photoelectron spectroscopy. More generally, analysis of the simulation data suggests that the surface activity of soft anions is largely a second or third solvation shell effect, arising from disruption of water-water hydrogen bonds and not from significant changes in first-shell anion-water hydrogen bonding.

Electron Attachment to Cytosine: The Role of Water

We present an EOM-CCSD-based quantum mechanical/molecular mechanical (QM/MM) study on the electron attachment process to solvated cytosine. The electron attachment in the bulk solvated cytosine occurs through a doorway mechanism, where the initial electron is localized on water. The electron is subsequently transferred to cytosine by the mixing of electronic and nuclear degrees of freedom, which occurs on an ultrafast time scale. The bulk water environment stabilizes the cytosine-bound anion by an extensive hydrogen-bond network and drastically enhances the electron transfer rate from that observed in the gas phase. Microhydration studies cannot reproduce the effect of the bulk water environment on the electron attachment process, and one needs to include a large number of water molecules in the calculation to obtain converged results. The predicted adiabatic electron affinity and electron transfer rate obtained from our QM/MM calculations are consistent with the available experimental results.