Algebraic Diagrammatic Construction Scheme Research Articles

Efficient and Parallel Implementation of Real and Complex Response Functions Employing the Second-Order Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator.

We present the implementation of an efficient matrix-folded formalism for the evaluation of complex response functions and the calculation of transition properties at the level of the second-order algebraic-diagrammatic construction (ADC(2)) scheme. The underlying algorithms, in combination with the adopted hybrid MPI/OpenMP parallelization strategy, enabled calculations of the UV/vis spectra of a guanine oligomer series ranging up to 1032 contracted basis functions, thereby utilizing vast computational resources from up to 32,768 CPU cores. Further analysis of the convergence behavior of the involved iterative subspace algorithms revealed the superiority of a frequency-separated treatment of response equations even for a large spectral window, including 101 frequencies. We demonstrate the applicability to general quantum mechanical operators by the first reported electronic circular dichroism spectrum calculated with a complex polarization propagator approach at the ADC(2) level of theory.

Cotton-Mouton Effect Using Algebraic Diagrammatic Construction Schemes in the Intermediate State Representation Formalism.

The Cotton-Mouton effect is theoretically investigated for a selected set of molecules by using a novel computational methodology based on algebraic diagrammatic construction (ADC) schemes in the intermediate state representation (ISR) formulation. Therefore, the electronic contributions to the frequency-dependent polarizabilities and, for the first time, to the magnetizabilities as well as mixed electric and magnetic hypermagnetizabilities have been computed in the ADC/ISR framework. In addition to calculation of the Cotton-Mouton constant and the birefringence, the gauge origin dependence of the computed tensors and the applied methodology are thoroughly investigated. The new ADC/ISR methodology, employing the recently presented responsefun package, is applied to a test set of Ne and small molecules (H2, HF, O2, CO2, and benzene) and compared to data from the experiment as well as other ab initio methods. The presented theoretical ab initio ADC/ISR approach is a substantial extension of the available computational methods for the investigation of complex nonlinear properties, however, with a gauge origin dependence inherent to the method that decreases with increasing perturbation order.

Calibrating TDDFT Calculations of the X-ray Emission Spectrum of Liquid Water: The Effects of Hartree-Fock Exchange.

The structure and dynamics of liquid water continue to be debated, with insight provided by, among others, X-ray emission spectroscopy (XES), which shows a split in the high-energy 1b1 feature. This split is yet to be reproduced by theory, and it remains unclear if these difficulties are related to inaccuracies in dynamics simulations, spectrum calculations, or both. We investigate the performance of different methods for calculating XES of liquid water, focusing on the ability of time-dependent density functional theory (TDDFT) to reproduce reference spectra obtained by high-level coupled cluster and algebraic-diagrammatic construction scheme calculations. A metric for evaluating the agreement between theoretical spectra termed the integrated absolute difference (IAD), which considers the integral of shifted difference spectra, is introduced and used to investigate the performance of different exchange-correlation functionals. We find that computed spectra of symmetric and asymmetric model water structures are strongly and differently influenced by the amount of Hartree-Fock exchange, with best agreement to reference spectra for ∼40-50%. Lower percentages tend to yield high density of contributing states, resulting in too broad features. The method introduced here is useful also for other spectrum calculations, in particular where the performance for ensembles of structures are evaluated.

Theoretical insights into ultrafast excited state proton transfer coupled with twisted intramolecular charge transfer mechanism in 7-(2ʹ-pyridyl)indole: Effect of hydrogen bonding

The ultrafast proton transfer (PT) process in the excited state of 7-(2′-pyridyl)indole (abbreviated as 7PyIn) has been elucidated by on-the-fly dynamics simulations using the algebraic-diagrammatic construction scheme of second order (ADC(2)) method under the resolution of identity (RI) approximation. The tautomers of 7PyIn alone and associated with an acetonitrile molecule in the 7PyIn-CH3CN complex are most likely found after the excited state intramolecular proton transfer (ESIntraPT) within an ultrafast time scale between 40 and 60 fs. While a methanol molecule in the complex of 7PyIn-CH3OH could establish both ESIntraPT and excited-state intermolecular proton transfer (ESInterPT). ESInterPT of 7PyIn-CH3OH possibly occurs through a concerted process via the intermolecular hydrogen bonding bridge of the methanol molecule almost twice faster than that of ESIntraPT. This slower circumstance of the ESIntraPT process found in 7PyIn-CH3OH may be caused by the competition between formations of the intramolecular and intermolecular hydrogen bonds between 7PyIn and methanol. In addition, after the accomplishment of excited state proton transfer (ESPT), the initiation of conical intersection could be taking place by the structural twist of 7PyIn. The geometries, energies, proton probabilities, and provided mechanisms via the dynamics simulations of 7PyIn and its solvent clusters exhibiting ESIntraPT and ESInterPT are explored to provide the dynamic pictures of the ESPT processes at the molecular level.

Algebraic Diagrammatic Construction Schemes Employing the Intermediate State Formalism: Theory, Capabilities, and Interpretation.

Algebraic diagrammatic construction (ADC) schemes represent a family of ab initio methods for the calculation of excited electronic states and electron-detached and -attached states. All ADC methods have been demonstrated to possess great potential for molecular applications, e.g., for the calculation of absorption or photoelectron spectra or electron attachment processes. ADC originates from Green's function or propagator theory; however, most recent ADC developments heavily rely on the intermediate state representation or effective Liouvillian formalisms, which comprise new ADC methods and computational schemes for high-order properties. The different approaches for the calculation of excitation energies, ionization potentials, and electron affinities are intimately related, and they provide a coherent description of these quantities at equivalent levels of theory and with comparable errors. Most quantum chemical program packages contain ADC methods; however, the most complete ADC suite of methods can be found in the recent release of Q-Chem.

Consistent third-order one-particle transition and excited-state properties within the algebraic-diagrammatic construction scheme for the polarization propagator.

The intermediate state representation (ISR) formalism allows for the straightforward calculation of excited state properties and state-to-state transition moments using the algebraic-diagrammatic construction (ADC) scheme for the polarization propagator. Here, the derivation and implementation of the ISR in third-order perturbation theory for the one-particle operator are presented, enabling, for the first time, the calculation of consistent third-order ADC [ADC(3)] properties. The accuracy of ADC(3) properties is evaluated with respect to high-level reference data and compared to the previously used ADC(2) and ADC(3/2) schemes. Oscillator strengths and excited state dipole moments are computed, and typical response properties are considered: dipole polarizabilities, first-order hyperpolarizabilities, and two-photon absorption strengths. The consistent third-order treatment of the ISR leads to an accuracy similar to that of the mixed-order ADC(3/2) method; the individual performance, however, depends on the property and molecule under investigation. ADC(3) produces slightly improved results in the case of oscillator strengths and two-photon absorption strengths, while excited state dipole moments, dipole polarizabilities, and first-order hyperpolarizabilities exhibit similar accuracy at ADC(3) and ADC(3/2) levels. Taking the significant increase of central processing unit time and memory requirements of the consistent ADC(3) approach into account, the mixed-order ADC(3/2) scheme offers a better compromise between accuracy and efficiency for the properties considered.

Tuning UV Pump X-ray Probe Spectroscopy on the Nitrogen K Edge Reveals the Radiationless Relaxation of Pyrazine: Ab Initio Simulations Using the Quasiclassical Doorway-Window Approximation.

Transient absorption UV pump X-ray probe spectroscopy has been established as a versatile technique for the exploration of ultrafast photoinduced dynamics in valence-excited states. In this work, an ab initio theoretical framework for the simulation of time-resolved UV pump X-ray probe spectra is presented. The method is based on the description of the radiation-matter interaction in the classical doorway-window approximation and a surface-hopping algorithm for the nonadiabatic nuclear excited-state dynamics. Using the second-order algebraic-diagrammatic construction scheme for excited states, UV pump X-ray probe signals were simulated for the carbon and nitrogen K edges of pyrazine, assuming a duration of 5 fs of the UV pump and X-ray probe pulses. It is predicted that spectra measured at the nitrogen K edge carry much richer information about the ultrafast nonadiabatic dynamics in the valence-excited states of pyrazine than those measured at the carbon K edge.

Scaled opposite-spin atomic-orbital based algebraic diagrammatic construction scheme for the polarization propagator with asymptotic linear-scaling effort: Theory and implementation.

A novel local approach for the quantum-chemical computation of excited states is presented, where the concept of the atomic-orbital formulation of the second-order Møller-Plesset energy expression is extended to the second-order algebraic diagrammatic construction scheme by virtue of the Laplace transform. The scaled opposite-spin second-order algebraic diagrammatic construction method with Cholesky decomposed densities and density-fitting, or CDD-DF-SOS-ADC(2) for short, exploits the sparsity of the two-electron repulsion integrals, the atomic ground-state density matrix, and the atomic transition density matrix to drastically reduce the computational effort. By using a local density-fitting approximation, it is shown that asymptotically linear scaling can be achieved for linear carboxylic acids. For electron-dense systems, sub-cubic scaling can be achieved if the excitation is local, and hence the transition density is sparse. Furthermore, the memory footprint and accuracy of the CDD-DF-SOS-ADC(2) method are explored in detail.

Perturbation theoretical approaches to strong light-matter coupling in ground and excited electronic states for the description of molecular polaritons.

Quantum chemical methods for the description of molecular polaritonic states in the strong coupling regime based on the Pauli-Fierz Hamiltonian are introduced. Based on a quantum electrodynamic Hartree-Fock (QED-HF) reference, a QED Møller-Plesset perturbation theory of second order for the electronic ground state and a second order quantum electrodynamic algebraic diagrammatic construction scheme for the polarization propagator [QED-ADC(2)] for excited electronic states have been derived, implemented, and tested for polaritons in hydrogen fluoride. Analogous approaches based on a standard non-polaritonic HF reference are also presented and thoroughly compared, both algebraically and numerically, to those based on the QED-HF reference. Furthermore, a promising route to approximate QED-ADC methods based on a unitary transformation of the algebraic expression into a restricted state space is outlined showing excellent agreement in second order with QED-ADC(2). All presented novel methods are compared to and tested against other existing abinitio approaches, mostly QED coupled cluster theory, including single and double excitations, and show qualitative agreement at a reduced computational effort.

The Multistate Quantum Monte Carlo Algebraic Diagrammatic Construction Method.

A multistate formulation of the recently developed quantum Monte Carlo (QMC) algebraic diagrammatic construction (ADC) method, QMCADC, is presented. QMCADC solves the Hermitian eigenvalue problem of the second-order ADC scheme for the polarization propagator stochastically by combining ADC schemes with projector quantum Monte Carlo (PQMC). It allows for massively parallel distributed computing and exploits the sparsity of the effective ADC matrix, thereby relaxing memory and processing requirements of ADC methods significantly. Here, the theory and implementation of the multistate variant of QMCADC are described, and our first proof-of-principle calculations for various molecular systems are shown. Indeed, multistate QMCADC enables sampling of an arbitrary number of low-lying excited states and can reproduce their vertical excitation energies with a marginal controllable error. The performance of multistate QMCADC is examined in terms of state-wise and overall accuracy as well as with respect to the balance in the treatments of excited states relatively to each other. The results are very promising as they show bias and imbalances among excited states to diminish as the number of sampling points increases. Furthermore, the impact of the quality of trial wave functions on the vertical excitation energies is investigated. A black-box approach for the generation of high quality trial wave functions internally is given.

The fourth-order algebraic diagrammatic construction scheme for the polarization propagator.

Until today, perturbation-theoretical consistent algebraic diagrammatic construction (ADC) schemes for the polarization propagator had been derived and implemented up to third order. They have turned out to be versatile and reliable ab initio single-reference methods for the quantum chemical investigation of electronic transitions as well as excited-state properties. Here we present, for the first time, the derivation of consistent fourth-order ADC(4) schemes exploiting novel techniques of automated equation and code generation. The accuracies of the resulting ADC(4) excitation energies have been benchmarked against recent high-level, near exact reference data. The mean absolute error for singly and doubly excited states turns out to be smaller than 0.1 and 0.5eV, respectively. These developments open also new avenues toward highly accurate ADC methods for electron-detached and attached states.

Magnetic circular dichroism within the algebraic diagrammatic construction scheme of the polarization propagator up to third order.

We present an implementation of the B term of Magnetic Circular Dichroism (MCD) within the Algebraic Diagrammatic Construction (ADC) scheme of the polarization propagator and its Intermediate State Representation. As illustrative results, the MCD spectra of the ADC variants ADC(2), ADC(2)-x, and ADC(3) of the molecular systems uracil, 2-thiouracil, 4-thiouracil, purine, hypoxanthine 1,4-naphthoquinone, 9,10-anthraquinone, and 1-naphthylamine are computed and compared with results obtained by using the Resolution-of-Identity Coupled-Cluster Singles and Approximate Doubles method, with literature Time-Dependent Density Functional Theory results, and with available experimental data.

Exploring the accuracy and usefulness of semi-empirically scaled ADC schemes by blending second and third order terms.

Different approaches to mixed-order algebraic-diagrammatic construction (ADC) schemes are investigated. The performance of two different strategies for scaling third-order contributions to the ADC secular matrix is evaluated. Both considered schemes employ a single tuning parameter and conserve general properties inherent to all ADC methods, such as hermiticity and size-consistency. The first approach, scaled-matrix ADC[(2) + x(3)], scales all contributions first occurring in ADC(3) equally and leads to an improvement of the accuracy of excitation energies compared to ADC(3) for x = 0.4-0.5. However, with respect to excited state dipole moments, this method provides lower accuracy than ADC(3). The second scaling approach, MP[(1) + x(2)] - ISR(3), scales the second order contributions of the ground-state wavefunction and derives a rigorous ADC scheme via the intermediate state representation formalism. Although the error in excitation energies is not improved, this method provides insight into the relevance of the individual terms of the ADC(3) matrix and indicates that the MP(2) wavefunction is, indeed, the optimal reference wavefunction for deriving a third-order single-reference ADC scheme.

Theoretical analysis and comparison of unitary coupled-cluster and algebraic-diagrammatic construction methods for ionization.

This article describes a novel approach for the calculation of ionization potentials (IPs), or, more generally, electron-detachment energies, based on a unitary coupled-cluster (UCC) parameterization of the ground-state wave function. Explicit working equations for a scheme referred to as IP-UCC3 are given, providing electron-detachment energies and spectroscopic amplitudes of electron-detached states dominated by one-hole excitations correct through third order. In the derivation, an expansion of the UCC transformed Hamiltonian involving Bernoulli numbers as expansion coefficients is employed. Both the secular matrix and the effective transition moments are shown to be essentially equivalent to the strict third-order algebraic-diagrammatic construction scheme for the electron propagator (IP-ADC). Interestingly, due to the Bernoulli expansion, neglecting triple substitutions in the UCC expansion manifold does not affect the third-order consistency of the IP-UCC effective transition moments. Finally, the equivalence between ADC and UCC excited-state schemes is shown to not hold in fourth or higher order due to a different treatment of the correlated excited-state basis.

Quantum Monte Carlo formulation of the second order algebraic diagrammatic construction: Toward a massively parallel correlated excited state method.

The quantum Monte Carlo (QMC) algebraic diagrammatic construction (ADC) method is introduced, which solves the eigenvalue problem of the second-order ADC scheme for the polarization propagator stochastically within the framework of QMC methodology allowing for massively parallel computations. As common virtue of the Monte Carlo integration techniques, quantum Monte Carlo algebraic diagrammatic construction (QMCADC) enables exploitation of the sparsity of the effective ADC matrix, and it reduces the memory requirements by storing only a portion of configurations at each iteration. Furthermore, distributing memory and processing loads to different computing nodes enables the use of fast developing parallel computing resources. Here, the theory and implementation of QMCADC is reported and its viability is demonstrated by the first proof-of-principle calculations. The focus lies on the first excited state and the reproduction of the corresponding lowest vertical excitation energy of various molecular systems. QMCADC is shown to be a genuine stochastic solution of the ADC eigenvalue problem, and exact ADC values can be obtained with a marginal controllable error.

Using core-hole reference states for calculating X-ray photoelectron and emission spectra.

For the calculation of core-ionization energies (IEs), X-ray photoelectron spectra (XPS), and X-ray emission spectra (XES), a commonly applied approach is to use non-Aufbau reference states with a core-hole as either final (IE and XPS) or initial (XES) state. However, such reference states can introduce numerical instabilities in post-HF methods, relating to the denominator of the energy corrections involved. This may become arbitrarily close to zero if a negative virtual orbital is present, e.g. a core-hole, leading to near-singularities. The resulting instabilities lead to severe convergence issues of the calculation schemes and, in addition, can strongly affect both energies and intensities, with oscillator strengths seen to reach values up to 4 × 107. For the K-edge we propose freezing the highest-energy virtual orbitals which contribute to any denominator below a threshold of 0.1 Hartree. Stable and reliable spectra are then produced, with minimal influence due to freezing energetically high-lying virtual orbitals (typically removing <5% of the total number of MOs). The developed protocol is here tested for Møller-Plesset perturbation theory and for the algebraic diagrammatic construction scheme for the polarization propagator, and it is also relevant for coupled cluster theory and other related methods.

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 2Πg resonance of the nitrogen anion N2 -. 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.

Ab Initio Excited-State Electronic Circular Dichroism Spectra Exploiting the Third-Order Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator.

Excited-state rotatory strengths are reported for the first time at a correlated ab initio level, here with the algebraic diagrammatic construction scheme of the polarization propagator up to the third order. To demonstrate the capabilities of this computational approach, the gas phase S1 electronic circular dichroism spectra of the bicyclic ketones (1R)-camphor, (1R)-norcamphor, and (1R)-fenchone have been calculated at the ADC(3) level of theory. Furthermore, the solution excited-state spectra of the energetically lowest conformer of R-(+)-1,1'-bi(2-naphthol) have been computed with inclusion of a polarizable continuum model at the ADC(2) level of theory.

Gator: A Python‐driven program for spectroscopy simulations using correlated wave functions

AbstractThe Gator program has been developed for computational spectroscopy and calculations of molecular properties using real and complex propagators at the correlated level of wave function theory. Currently, the focus lies on methods based on the algebraic diagrammatic construction (ADC) scheme up to the third order of perturbation theory. An auxiliary Fock matrix‐driven implementation of the second‐order ADC method for excitation energies has been realized with an underlying hybrid MPI/OpenMP parallelization scheme suitable for execution in high‐performance computing cluster environments. With a modular and object‐oriented program structure written in a Python/C++ layered fashion, Gator additionally enables time‐efficient prototyping of novel scientific approaches, as well as interactive notebook‐driven training of students in quantum chemistry.This article is categorized under: Computer and Information Science &gt; Computer Algorithms and Programming Electronic Structure Theory &gt; Ab Initio Electronic Structure Methods Software &gt; Quantum Chemistry

Deciphering excited state properties utilizing algebraic diagrammatic construction schemes of decreasing order.

Excited state properties are difficult to trace back to the common molecular orbital picture when the excited state wavefunction is a linear combination of two or more Slater determinants. Here, a theoretical methodology is introduced based on the algebraic diagrammatic construction scheme for the polarization propagator (ADC(n)) that allows to make this connection and to eventually derive structure-function relationships. The usefulness of this approach is demonstrated by an analysis of the transition dipole moments of the low-lying 1B3u and 2B3u states of anthracene and (1,4,5,8)-tetraazaanthracene.