Complex absorbing potentials in the framework of electron propagator theory. II. Application to temporary anions

A New Strategy to Optimize Complex Absorbing Potentials for the Computation of Resonance Energies and Widths

Complex absorbing potentials (CAPs) are artificial potentials added to electronic Hamiltonians to make the wave function of metastable electronic states square-integrable. This makes the electronic-structure theory of resonances comparable to that of bound states, thus reducing the complexity of the problem. However, the most often used box and Voronoi CAPs depend on several parameters that have a substantial impact on the numerical results. Among these parameters are the CAP strength and a set of spatial parameters that define the onset of the CAP. It has been a common practice to minimize the perturbation of the resonance states due to the CAP by optimizing the strength parameter while fixing the onset parameters, although the performance of this approach strongly depends on the chosen onset. Here, we introduce a more general approach that allows one to optimize not only the CAP strength but also the spatial parameters. We show that fixing the CAP strength and optimizing the spatial parameters is a reliable way to minimize CAP perturbations. We illustrate the performance of this new approach by computing resonance energies and widths of the temporary anions of dinitrogen, ethylene, and formic acid. This is done at the Hartree-Fock and equation-of-motion coupled-cluster singles and doubles levels of theory using full and projected box and Voronoi CAPs.

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.

Investigation of negative-ion resonances using a subspace-projected multiconfigurational electron propagator perturbed with a complex absorbing potential

The transient negative-ion resonances found in scattering experiments are important intermediates in many chemical processes. These metastable states correspond to the continuum part of the Hamiltonian of the projectile-target composite system. Usual bound-state electronic structure methods are not applicable for these. In this work, we develop a subspace-projection method in connection with an electron propagator (EP) defined in terms of a complete-active-space self-consistent-field initial state. The target Hamiltonian (Ĥ) is perturbed by a complex absorbing potential (CAP) for the analytical continuation of the spectrum of Ĥ to complex eigenvalues associated with the continuum states. The resonance is identified as a pole of the EP, which is stable with respect to variations in the strength of the CAP. The projection into a small subspace reduces the size of the complex matrices to be diagonalized, minimizes the computational cost, and affords some insight into the orbitals that are likely to play some role in the capture of the projectile. Two molecular (Πg2N₂ ^- and ²Π CO^-) and an atomic shaperesonance (²P Be^-) are investigated using this method. The position and width of the resonances are in good agreement with the previously reported values.

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.

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.

Intruder Peak-Free Transient Inner-Shell Spectra Using Real-Time Simulations

Real-time methods are convenient for simulating core-level absorption spectra but suffer from nonphysical intruder peaks when using atom-centered basis sets. In transient absorption spectra, these peaks exhibit highly nonphysical time-dependent modulations in their energies and oscillator strengths. In this paper, we address the origins of these intruder peaks and propose a straightforward and effective solution based on a filtered dipole operator. In combination with real-time time-dependent density functional theory (RT-TDDFT), we demonstrate how to compute intruder-free attosecond transient X-ray absorption spectra for the aminophenol (C₆H₇NO) oxygen and nitrogen K-edges and the α-quartz (SiO₂) silicon L-edge. Without filtering, the computed spectra are qualitatively wrong. This procedure is suitable for both static and transient inner-shell spectroscopy studies and can easily be implemented in a range of real-time methodologies.

Calculation of the Lowest Resonant States of H− and Li by the Complex Absorbing Potential Method

The analysis of the features of the method of complex absorbing potential (CAP) is carried out for a single-channel problem with an explicit parameterization of the scattering matrix. It is shown that there can be several types of CAP trajectories depending on the choice of the initial conditions. In any case, the estimation of the resonance parameters from the position of the optimal trajectory point can lead to a systematic error or an ambiguous result. In special cases, the search for the optimal point can be replaced by the averaging over a closed section of the trajectory. The CAP trajectories constructed in the H− and Li resonance calculations correlate well with the model trajectories, which have a curl around the resonance. The averaging over a closed area of the trajectory leads to better estimates of the energy and width of the resonance in comparison with the technique of searching for the optimal point.

CAP/EA-ADC method for metastable anions: Computational aspects and application to π* resonances of norbornadiene and 1,4-cyclohexadiene

The second- and third-order algebraic-diagrammatic construction schemes for the electron propagator for studies of electron attachment processes [EA-ADC(2) and EA-ADC(3)] have been extended to include the complex absorbing potential (CAP) method for the treatment of electronic resonances. Theoretical and conceptual aspects of the new CAP/EA-ADC methodology are studied in detail at the example of the well-known ²Π_g resonance of the nitrogen anion N₂ ^-. The methodology is further applied to π* shape resonances, for which ethylene is considered as a prototype. Furthermore, the first many-body treatment of the π₊ ^* and π_- ^* resonances of norbornadiene and 1,4-cyclohexadiene is provided, which have served as model systems for the concept of through-space and through-bond interactions for a long time.

Exploring spin symmetry-breaking effects for static field ionization of atoms: Is there an analog to the Coulson-Fischer point in bond dissociation?

Löwdin's symmetry dilemma is an ubiquitous issue in approximate quantum chemistry. In the context of Hartree-Fock (HF) theory, the use of Slater determinants with some imposed constraints to preserve symmetries of the exact problem may lead to physically unreasonable potential energy surfaces. On the other hand, lifting these constraints leads to the so-called broken symmetry solutions that usually provide better energetics, at the cost of losing information about good quantum numbers that describe the state of the system. This behavior has previously been extensively studied in the context of bond dissociation. This paper studies the behavior of different classes of HF spin polarized solutions (restricted, unrestricted, and generalized) in the context of ionization by strong static electric fields. We find that, for simple two electron systems, unrestricted Hartree-Fock (UHF) is able to provide a qualitatively good description of states involved during the ionization process (neutral, singly ionized, and doubly ionized states), whereas RHF fails to describe the singly ionized state. For more complex systems, even though UHF is able to capture some of the expected characteristics of the ionized states, it is constrained to a single M_s (diabatic) manifold in the energy surface as a function of field intensity. In this case, a better qualitative picture can be painted by using generalized Hartree-Fock as it is able to explore different spin manifolds and follow the lowest solution due to lack of collinearity constraints on the spin quantization axis.

An Electron Propagator Approach Based on a Multiconfigurational Reference State for the Investigation of Negative-Ion Resonances Using a Complex Absorbing Potential Method

Negative-ion resonances are important metastable states that result from the collision between an electron and a neutral target. The course of many chemical processes in nature is often dictated by how an intermediate resonance state falls apart. This article reports on the development of an electron propagator (EP) based on a Hamiltonian Ĥ perturbed by a complex absorbing potential (CAP) and a multiconfigurational self-consistent field (MCSCF) initial state to study these resonances. Perturbation of Ĥ by a CAP makes the resonances amenable to a bound-state method like MCSCF. Resonances stand out among the non-resonant states as persistent complex eigenvalues of the perturbed Ĥ when the strength (η) of the CAP is varied. The MCSCF method gives a reliable and accurate description of the target states, especially when the non-dynamical correlations are dominant. The resonance energies are obtained from the poles of the EP. We propose three variants of our EP depending on how the effect of the CAP is introduced. We find that the computationally most efficient variant is the one in which the reference state of the EP is an unperturbed MCSCF wavefunction and a non-zero CAP is defined only on the virtual orbital subspace of the reference state. The onset of the CAP is carefully optimized in order to minimize the artifacts due to reflections from the CAP. An extrapolation method (based on a Padé approximant) and a de-perturbation method are adopted in order to account for the limitations of finite basis sets used and determine the resonance energy in the limit of η → 0. ²P Be^-, ²Π_g N₂^-, and ²Π CO^- shape resonances are investigated. The position and width of these resonances computed in this study agree well with those reported earlier in the literature.

Combination of a Voronoi-Type Complex Absorbing Potential with the XMS-CASPT2 Method and Pilot Applications

Electronic resonances are metastable (N + 1) electron states, in other words, discrete states embedded in an electronic continuum. While great progress has been made for certain types of resonances-for example, temporary anions created by attaching one excess electron to a closed shell neutral-resonances in general remain a great challenge of quantum chemistry because a successful description of the decay requires a balanced description of the bound and continuum aspect of the resonance. Here, a smoothed Voronoi complex absorbing potential (CAP) is combined with the XMS-CASPT2 method, which enables us to address the balance challenge by appropriate choice of the CAS space. To reduce the computational cost, the method is implemented in the projected scheme. In this pilot application, three temporary anions serve as benchmarks: the π* resonance state of formaldehyde; the π* and σ* resonance states of chloroethene as functions of the C-Cl bond dissociation coordinate; and the ⁴Π_u and ²Π_u resonance states of N₂^-. The convergence of the CAP/XMS-CASPT2 results has been systematically examined with respect to the size of the active space. Resonance parameters predicted by the CAP/XMS-CASPT2 method agree well with CAP/SAC-CI results (deviations of about 0.15 eV); however, as expected, CAP/XMS-CASPT2 has clear advantages in the bond dissociation region. The advantages of CAP/XMS-CASPT2 are further demonstrated in the calculations of ⁴Π_u and ²Π_u resonance states of N₂^- including their ³Σ_u⁺ and ³Δ_u parent states. Three of the involved states (²Π_u, ³Σ_u⁺, and ³Δ_u) possess multireference character, and CAP/XMS-CASPT2 can easily describe these states with a relatively modest active space.

Simple DFT Scheme for Estimating Negative Electron Affinities

A simple density functional theory (DFT) scheme is proposed for estimating negative vertical electron affinities of neutral systems, based on a consideration of the integer discontinuity and density scaling homogeneity. The key feature is the derivation of two system-dependent exchange-correlation functionals, one appropriate for the electron deficient side of the integer and one appropriate for the electron abundant side. The electron affinity is evaluated as a linear combination of frontier orbital energies from self-consistent Kohn-Sham calculations on the neutral system using these functionals. For two assessments comprising a total of 43 molecules, the scheme provides electron affinities that are in good agreement with experimental values and which are an improvement over those from the DFT method of Tozer and De Proft [ J. Phys. Chem. A 2005 , 109 , 8923 ]. The scheme is trivial to implement in any Kohn-Sham program, and the computational cost is that of a series of generalized gradient approximation Kohn-Sham calculations. More generally, the study provides a prescription for performing low-cost, self-consistent Kohn-Sham calculations that yield frontier orbital energies that approximately satisfy the appropriate Koopmans conditions, without the need for exact exchange.

Resonance Energies and Lifetimes from the Analytic Continuation of the Coupling Constant Method: Robust Algorithms and a Critical Analysis

The energy of a metastable state can be computed by adding an artificial stabilizing potential to the Hamiltonian, increasing the stabilization until the metastable state is turned into a bound one, and then further increasing the stabilization until enough bound-state data have been collected so that these can be extrapolated back to vanishing stabilization. The lifetime of the metastable state can be obtained from the same data, but only if the extrapolation is performed by analytic continuation. This extrapolation method is called analytic continuation of the coupling constant (ACCC). Here we introduce preconditioning schemes for two of the three established extrapolation algorithms and critically compare results from all three extrapolation schemes in a variety of situations: As examples for resonance states serve the π* temporary anions of ethylene and formaldehyde as well as a model potential, which provides a case where input data with full numeric precision are available. In the data collection step, three different stabilizing potentials are employed, a Coulomb potential, a short-range Coulomb potential, and a soft-box Voronoi potential. Effects of different orders of the extrapolating Padé approximant are investigated, and last, the energy range of input data for the extrapolation is studied. Moreover, all ACCC results are compared to resonance parameters that have been independently obtained with the same theoretical method, but with a different continuum approach-complex scaling for the model and complex absorbing potentials for the temporary anions.

Ab initio calculation of the cross sections for electron impact vibrational excitation of CO via the (2)Π shape resonance

The stabilization method is used to calculate the complex potential energy curve of the (2)Π state of CO(-) as a function of bond length, with the refinement that separate potentials are determined for p-wave and d-wave attachment and detachment of the excess electron. Using the resulting complex potentials, absolute vibrational excitation cross sections are calculated as a function of electron energy and scattering angle. The calculated cross sections agree well with experiment.

Anions of xenon clusters bound by long-range electron correlations

In contrast with the single atom, atomic van der Waals clusters can form stable anions where the excess electron is bound due to long-range correlations with the electrons of the cluster. We report on extensive all-electron many-body ab initio studies on Xe clusters. Three-dimensional, planar, and linear structures of the clusters are investigated and compared. In particular, we find that the minimal number of Xe atoms in the cluster required to form a stable anion is 5 independently of the dimensionality of the cluster. We provide electron affinities for clusters made of 5, 6, and 7 atoms in all dimensions and find that the planar clusters form the most stable anions. The Dyson orbitals of the excess electrons are computed and analyzed.