The Clusterization Technique: A Systematic Search for the Resonance Energies Obtained via Padé

Shape resonance induced electron attachment to cytosine: The effect of aqueous media

We have investigated the impact of microsolvation on shape-type resonance states of nucleobases, taking cytosine as a case study. To characterize the resonance position and decay width of the metastable states, we employed the newly developed DLPNO-based EA-EOM-CCSD method in conjunction with the resonance via Padé (RVP) method. Our calculations show that the presence of water molecules causes a redshift in the resonance position and an increase in the lifetime for the three lowest-lying resonance states of cytosine. Furthermore, there are some indications that the lowest resonance state in isolated cytosine may get converted to a bound state in the presence of an aqueous environment. The obtained results are extremely sensitive to the basis set used for the calculations.

Use of bound state methods to calculate partial and total widths of shape resonances

In this work we study the 2Π resonances of a two-site model system designed to mimic a smooth transition from the 2Πg temporary anion of N2 to the 2Π temporary anion of CO. The model system possesses the advantage that scattering and bound state (L2) methods can be directly compared without obfuscating electron-correlation effects. Specifically, we compare resonance parameters obtained with the complex Kohn variational (CKV) method with those from stabilization, complex absorbing potential, and regularized analytical continuation calculations. The CKV calculations provide p-wave and d-wave widths, the sum of which provides a good approximation of the total width. Then we demonstrate that the width obtained with modified bound state methods depends on the basis set employed: It can be the total width, a partial width, or an ill-defined sum of partial widths. Provided the basis set is chosen appropriately, widths from bound state methods agree well with the CKV results.

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.

Theory of electronic resonances: fundamental aspects and recent advances

Electronic resonances are states that are unstable towards loss of electrons. They play critical roles in high-energy environments across chemistry, physics, and biology but are also relevant to processes under ambient conditions that involve unbound electrons. This feature article focuses on complex-variable techniques such as complex scaling and complex absorbing potentials that afford a treatment of electronic resonances in terms of discrete square-integrable eigenstates of non-Hermitian Hamiltonians with complex energy. Fundamental aspects of these techniques as well as their integration into molecular electronic-structure theory are discussed and an overview of some recent developments is given: analytic gradient theory for electronic resonances, the application of rank-reduction techniques and quantum embedding to them, as well as approaches for evaluating partial decay widths.

Analysis of Stabilization and Extrapolation Methods for Determining Energies and Lifetimes of Metastable Electronic States

Two methods that make use of standard electronic structure tools, the stabilization and extrapolation methods, are discussed with an eye toward pointing out their relative strengths and weaknesses and for improving their applications. In the former, whether to utilize energy data from only one or from both branches of an avoided crossing between the quasi-bound and pseudo-continuum states is one issue that is focused on. Another is the decision of where along the stabilization plot's branches (i.e., far from or close to the avoided crossing) to create data points for optimal performance given a reasonable (10^-5-10^-7 eV) precision in the electronic energy. A third issue is how many parameters to use in fitting energy data to the (one or two) branches of the stabilization plot. In extrapolation methods, one uses energy data computed when the metastable state's energy has been rendered stable by the application of an external potential, which thus produces a one-branch function. The main issues in implementing this method are the functional form for how the energy E depends on the strength of the external potential especially as the energy evolves from the bound-state region toward the unbound region and how to choose data points so that energy values of a reasonable precision are capable of determining the parameters in the formula that produces the metastable state's energy E and half-width Γ/2 (inversely related to the state's lifetime). In addition to explaining, critiquing, and comparing these two methods, several suggestions are offered for their further testing and improvements.

Variational Solutions for Resonances by a Finite-Difference Grid Method

We demonstrate that the finite difference grid method (FDM) can be simply modified to satisfy the variational principle and enable calculations of both real and complex poles of the scattering matrix. These complex poles are known as resonances and provide the energies and inverse lifetimes of the system under study (e.g., molecules) in metastable states. This approach allows incorporating finite grid methods in the study of resonance phenomena in chemistry. Possible applications include the calculation of electronic autoionization resonances which occur when ionization takes place as the bond lengths of the molecule are varied. Alternatively, the method can be applied to calculate nuclear predissociation resonances which are associated with activated complexes with finite lifetimes.

Role of Overlap between the Discrete State and Pseudocontinuum States in Stabilization Calculations of Metastable States

In a diabatic picture metastable states subject to decay by electron detachment can be viewed as arising from the coupling between a discrete state and a continuum. In treating such states with bound-state quantum chemical methods, the continuum is discretized. In this study, we elucidate the role of overlap in this interaction in the application of the stabilization method to temporary anion states. This is accomplished by use of a minimalist stabilization calculation on the lowest energy l=2 (D) resonance of the finite spherical well potential using two basis functions, one describing the diabatic discrete state and the other a diabatic discretized continuum state. We show that even such a simple treatment predicts a complex resonance energy in good agreement with the exact result. If the energy of the discrete state is assumed to be constant, which is tantamount to orthogonalizing the discretized continuum state to the discrete state, it is demonstrated that the square of the off-diagonal coupling has a maximum close to the crossing point of the orthogonalized diabatic curves and that the curvature in the coupling is responsible for the complex stationary point associated with the resonance. Moreover, this curvature is a consequence of the overlap between the two diabatic states.

Uniform vs Partial Scaling within Resonances via Padé Based on the Similarities to Other Non-Hermitian Methods: Illustration for the Beryllium 1s22p3s State

Resonance via Padé (RVP) is an efficient method for calculating autoionization resonance states. It is based on the stabilization technique in which the basis set is scaled. The scaling can be uniform (i.e., all basis functions are scaled) or partial. Herein, we compare the two RVP scaling schemes for calculating an autoionization eigenvalue; moreover, the effect of freezing the core electrons is intertwined within this comparison. In order to study the different behavior of the RVP schemes, we associate each RVP scaling scheme with a complex contour of integration. Similarities between RVP and other non-Hermitian methods emerge from the generated contours, which suggest that RVP introduces similar outgoing boundary conditions as the complex scaling (CS), complex basis function (CBF), and reflection-free complex absorbing potential (RF-CAP) methods. A uniform-RVP contour, unlike a partial one, immediately penetrates the complex plane and influences the interaction region. Hence, uniform scaling within RVP destroys the description of the core electrons, as well as the description of the reference state, and yields less reliable results than partial scaling. The 1s²2p3s ¹P autoionization state of Be, at the equation-of-motion coupled-cluster level, is used as our case study model.

Quantum Effects Dominating the Interatomic Coulombic Decay of an Extreme System

LiHe is an extreme open-shell system. It is among the weakest bound systems known, and its mean interatomic distance extends dramatically into the classical forbidden region. Upon 1s → 2p excitation of He, interatomic Coulombic decay (ICD) takes place in which the electronically excited helium atom relaxes and transfers its excess energy to ionize the neighboring lithium atom. A substantial part of the decay is found to be to the dissociation continuum producing Li⁺ and He atoms. The distribution of the kinetic energy released by the ICD products is found to be highly oscillatory. Its analysis reveals that quantum phase shifts between the decaying states and the dissociating final states are controlling this ICD reaction. The semiclassical reflection principle, which commonly explains ICD reactions, fails. The process is expected to be amenable to experiment.

Ab Initio Complex Transition Dipoles between Autoionizing Resonance States from Real Stabilization Graphs

Electronic transition dipoles are crucial for investigating light-matter interactions. Transition dipoles between metastable (autoionizing resonance) states become complex within non-Hermitian formalism, analogous to the resonance energies. Herein, we put forward a robust method for evaluating complex transition dipoles based on real ab initio stabilization calculations. The complex transition dipoles are obtained by analytical continuation via the Padé approximant and are identified as stationary solutions in the complex plane. The capability of the new approach is demonstrated for several transition dipoles of the doubly excited helium resonance states, for which exact values are available for comparison. Nevertheless, the method presented here has no inherent limitation and is suitable for polyatomic systems.

Evidence of Nonrigidity Effects in the Description of Low-Energy Anisotropic Molecular Collisions of Hydrogen Molecules with Excited Metastable Helium Atoms

Cold collisions serve as a sensitive probe of the interaction potential. In the recent study of Klein et al. (Nature Phys.2017, 13, 35–38), the one-parameter scaling of the interaction potential was necessary to obtain agreement between theoretical and observed patterns of the orbiting resonances for excited metastable helium atoms colliding with hydrogen molecules. Here, we show that the effect of nonrigidity of the H₂ molecule on the resonant structure, absent in the previous study, is critical to predict the correct positions of the resonances in that case. We have complemented the theoretical description of the interaction potential and revised reaction rate coefficients by proper inclusion of the flexibility of the molecule. The calculated reaction rate coefficients are in remarkable agreement with the experimental data without empirical adjustment of the interaction potential. We have shown that even state-of-the-art calculations of the interaction energy cannot ensure agreement with the experiment if such an important physical effect as flexibility of the interacting molecule is neglected. Our findings about the significance of the nonrigidity effects can be especially crucial in cold chemistry, where the quantum nature of molecules is pronounced.