No Arabic abstract
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 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 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 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 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 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.