Projected Complex Absorbing Potential Multireference Configuration Interaction Approach for Shape and Feshbach Resonances

Anion resonances are formed as metastable intermediates in low-energy electron-induced reactions. Due to the finite lifetimes of resonances, applying standard Hermitian formalism for their characterization presents a vexing problem for computational chemists. Numerous modifications to conventional quantum chemical methods have enabled satisfactory characterization of resonances, but specific issues remain, especially in describing two-particle one-hole (2p-1h) resonances. An accurate description of these resonances and their coupling with single-particle resonances requires a multireference approach. We propose a projected complex absorbing potential (CAP) implementation within the multireference configuration interaction (MRCI) framework to characterize single-particle and 2p-1h resonances. As a first application, we use the projected-CAP-MRCI approach to characterize and benchmark the ²Π_g shape resonance in N₂^-. We test its performance as a function of the size of the subspace and other parameters, and we compute the complex potential energy surface of the ²Π_g shape resonance to show that a smooth curve is obtained. One key benefit of MRCI is that it can describe Feshbach resonances (most common examples of 2p-1h resonances) at the same footing as shape resonances. Therefore, it is uniquely positioned to describe mixing between the different channels. To test these additional capabilities, we compute Feshbach resonances in H₂O^- and anions of dicyanoethylene isomers. We find that CAP-MRCI can efficiently capture the mixing between the Feshbach and shape resonances in dicyanoethylene isomers, which has significant consequences for their lifetimes.

Complex energies and transition dipoles for shape-type resonances of uracil anion from stabilization curves via Padé

Absorption of slow moving electrons by neutral ground state nucleobases has been known to produce resonance metastable states. There are indications that such metastable states may play a key role in DNA/RNA damage. Therefore, herein, we present an ab initio non-Hermitian investigation of the resonance positions and decay rates for the low lying shape-type states of the uracil anion. In addition, we calculate the complex transition dipoles between these resonance states. We employ the resonance via Padé (RVP) method to calculate these complex properties from real stabilization curves by analytical dilation into the complex plane. This method has already been successfully applied to many small molecular systems, and herein, we present the first application of RVP to a medium-sized system. The presented resonance energies are optimized with respect to the size of the basis set and compared with previous theoretical studies and experimental findings. Complex transition dipoles between the shape-type resonances are computed using the optimal basis set. The ability to calculate ab initio energies and lifetimes of biologically relevant systems paves the way for studying reactions of such systems in which autoionization takes place, while the ability to also calculate their complex transition dipoles opens the door for studying photo-induced dynamics of such biological molecules.

Modeling the Ultrafast Electron Attachment Dynamics of Solvated Uracil

Electron attachment to DNA by low energy electrons can lead to DNA damage, so a fundamental understanding of how electrons interact with the components of nucleic acids in solution is an open challenge. In solution, low energy electrons can generate presolvated electrons, e_pre^-, which are efficiently scavanged by pyrimidine nucleobases to form transient negative ions, able to relax to either stable valence bound anions or undergo dissociative electron detachment or transfer to other parts of DNA/RNA leading to strand breakages. In order to understand the initial electron attachment dynamics, this paper presents a joint molecular dynamics and high-level electronic structure study into the behavior of the electronic states of the solvated uracil anion. Both the valence π* and nonvalence e_pre^- states of the solvated uracil system are studied, and the effect of the solvent environment and the geometric structure of the uracil core are uncoupled to gain insight into the physical origin of the stabilization of the solvated uracil anion. Solvent reorganization is found to play a dominant role followed by relaxation of the uracil core.

Electron-induced origins of prebiotic building blocks of sugars: mechanism of self-reactions of a methanol anion dimer

Reactions of methanol dimers in interstellar medium driven by low energy irradiation may lead to prebiotic precursors.