No Arabic abstract
We present a multi-reference generalization of the algebraic diagrammatic construction theory (ADC) [J. Schirmer, Phys. Rev. A 26, 2395 (1982)] for excited electronic states. The resulting multi-reference ADC approach (MR-ADC) can be efficiently and reliably applied to systems, which exhibit strong electron correlation in the ground or excited electronic states. In contrast to conventional multi-reference perturbation theories, MR-ADC describes electronic transitions involving all orbitals (core, active, and external) and enables efficient computation of spectroscopic properties, such as transition amplitudes and spectral densities. Our derivation of MR-ADC is based on the effective Liouvillean formalism of Mukherjee and Kutzelnigg [D. Mukherjee, W. Kutzelnigg, in Many-Body Methods in Quantum Chemistry (1989), pp. 257--274], which we generalize to multi-determinant reference states. We discuss a general formulation of MR-ADC, perform its perturbative analysis, and present an implementation of the first-order MR-ADC approximation, termed MR-ADC(1), as a first step in defining the MR-ADC hierarchy of methods. We show results of MR-ADC(1) for the excitation energies of the Be atom, an avoided crossing in LiF, doubly excited states in C2, and outline directions for our future developments.
We present an implementation and benchmark of new approximations in multireference algebraic diagrammatic construction theory for simulations of neutral electronic excitations and UV/Vis spectra of strongly correlated molecular systems (MR-ADC). Following our work on the first-order MR-ADC approximation [J. Chem. Phys. 2018, 149, 204113], we report the strict and extended second-order MR-ADC methods (MR-ADC(2) and MR-ADC(2)-X) that combine the description of static and dynamic electron correlation in the ground and excited electronic states without relying on state-averaged reference wavefunctions. We present an extensive benchmark of the new MR-ADC methods for excited states in several small molecules, including the carbon dimer, ethylene, and butadiene. Our results demonstrate that for weakly-correlated electronic states the MR-ADC(2) and MR-ADC(2)-X methods outperform the third-order single-reference ADC approximation and are competitive with the results from equation-of-motion coupled cluster theory. For states with multireference character, the performance of the MR-ADC methods is similar to that of an N-electron valence perturbation theory. In contrast to conventional multireference perturbation theories, the MR-ADC methods have a number of attractive features, such as a straightforward and efficient calculation of excited-state properties and a direct access to excitations outside of the frontier (active) orbitals.
We present a second-order formulation of multi-reference algebraic diagrammatic construction theory [Sokolov, A. Yu. J. Chem. Phys. 2018, 149, 204113] for simulating photoelectron spectra of strongly correlated systems (MR-ADC(2)). The MR-ADC(2) method uses second-order multi-reference perturbation theory (MRPT2) to efficiently obtain ionization energies and intensities for many photoelectron transitions in a single computation. In contrast to conventional MRPT2 methods, MR-ADC(2) provides information about ionization of electrons in all orbitals (i.e., core and active) and allows to compute transition intensities in straightforward and efficient way. Although equations of MR-ADC(2) depend on four-particle reduced density matrices, we demonstrate that computation of these large matrices can be completely avoided without introducing any approximations. The resulting MR-ADC(2) implementation has a lower computational scaling compared to conventional MRPT2 methods. We present results of MR-ADC(2) for photoelectron spectra of small molecules, carbon dimer, and equally-spaced hydrogen chains (H10 and H30) and outline directions for future developments.
We present implementation of second- and third-order algebraic diagrammatic construction theory for efficient and accurate computations of molecular electron affinities (EA), ionization potentials (IP), and densities of states (EA-/IP-ADC(n), n = 2, 3). Our work utilizes the non-Dyson formulation of ADC for the single-particle propagator and reports working equations and benchmark results for the EA-ADC(2) and EA-ADC(3) approximations. We describe two algorithms for solving EA-/IP-ADC equations: (i) conventional algorithm that uses iterative diagonalization techniques to compute low-energy EA, IP, and density of states, and (ii) Greens function algorithm (GF-ADC) that solves a system of linear equations to compute density of states directly for a specified spectral region. To assess accuracy of EA-ADC(2) and EA-ADC(3), we benchmark their performance for a set of atoms, small molecules, and five DNA/RNA nucleobases. As our next step, we demonstrate efficiency of our GF-ADC implementation by computing core-level K-, L-, and M-shell ionization energies of a zinc atom without introducing core-valence separation approximation. Finally, we use EA- and IP-ADC methods to compute band gaps of equally-spaced hydrogen chains Hn with n up to 150, providing their estimates near thermodynamic limit. Our results demonstrate that EA-/IP-ADC(n) (n = 2, 3) methods are efficient and accurate alternatives to widely used electronic structure methods for simulations of electron attachment and ionization properties.
We report a new implementation of multireference algebraic diagrammatic construction theory (MR-ADC) for simulations of electron attachment and ionization in strongly correlated molecular systems (EA/IP-MR-ADC). Following our recent work on IP-MR-ADC [J. Chem. Theory Comput. 2019, 15, 5908], we present the first implementation of the second-order MR-ADC method for electron attachment (EA-MR-ADC(2)), as well as two extended second-order approximations (EA- and IP-MR-ADC(2)-X) that incorporate a partial treatment of third-order electron correlation effects. Introducing a small approximation for the second-order amplitudes of the effective Hamiltonian, our implementation of EA- and IP-MR-ADC(2)-X has a low O(M^5) computational scaling with the basis set size M. Additionally, we describe an efficient algorithm for solving the first-order amplitude equations in MR-ADC and partially-contracted second-order N-electron valence perturbation theory (NEVPT2) that completely avoids computation of the four-particle reduced density matrices without introducing any approximations or imaginary-time propagation. For a benchmark set of eight small molecules, carbon dimer, and a twisted ethylene, we demonstrate that EA- and IP-MR-ADC(2)-X achieve accuracy similar to that of strongly-contracted NEVPT2, while having a lower computational scaling with the active space size and providing efficient access to transition properties.
We present an efficient implementation of the second- and third-order single-reference algebraic diagrammatic construction theory for electron attachment (EA) and ionization (IP) energies and spectra (EA/IP-ADC(n), n = 2, 3). Our new EA/IP-ADC program features spin adaptation for closed-shell systems, density fitting for efficient handling of the two-electron integral tensors, as well as vectorized and parallel implementation of tensor contractions. We demonstrate capabilities of our efficient implementation by applying the EA/IP-ADC(n) (n = 2, 3) methods to compute the photoelectron spectrum of the TEMPO radical, as well as the vertical and adiabatic electron affinities of TEMPO and two DNA base pairs (guanine-cytosine and adenine-thymine). The spectra and electron affinities computed using large diffuse basis sets with up to 1028 molecular orbitals are found to be in a good agreement with the best available results from the experiment and theoretical simulations.