Software for the frontiers of quantum chemistry:An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design

Complex Excited State Polarizabilities in the ADC/ISR Framework

We present the derivation and implementation of complex, frequency-dependent polarizabilities for excited states using the algebraic-diagrammatic construction for the polarization propagator (ADC) and its intermediate state representation (ISR). Based on the complex polarizability we evaluate C6 dispersion coefficients for excited states. The methodology is implemented up to third order in perturbation theory in the Python-driven adcc toolkit for the development and application of ADC methods. We exemplify the approach using small model systems and compare it to results from coupled-cluster theory and from experiments.</div

adcc: A versatile toolkit for rapid development of algebraic-diagrammatic construction methods

30 pages, 7 figuresInternational audienceADC-connect (adcc) is a hybrid python/C++ module for performing excited state calculations based on the algebraic-diagrammatic construction scheme for the polarisation propagator (ADC). Key design goal is to restrict adcc to this single purpose and facilitate connection to external packages, e.g., for obtaining the Hartree-Fock references, plotting spectra, or modelling solvents. Interfaces to four self-consistent field codes have already been implemented, namely pyscf, psi4, molsturm, and veloxchem. The computational workflow, including the numerical solvers, are implemented in python, whereas the working equations and other expensive expressions are done in C++. This equips adcc with adequate speed, making it a flexible toolkit for both rapid development of ADC-based computational spectroscopy methods as well as unusual computational workflows. This is demonstrated by three examples. Presently, ADC methods up to third order in perturbation theory are available in adcc, including the respective core-valence separation and spin-flip variants. Both restricted or unrestricted Hartree-Fock references can be employed

Analytical Gradients for Core-Excited States in the Algebraic Diagrammatic Construction (ADC) Framework

Here we present a derivation of the analytical expressions required to determine nuclear gradients for core-excited states at the core-valence separated algebraic diagrammatic construction (CVS-ADC) theory level. Analytical gradients up to and including the extended CVS-ADC(2)-x order have been derived and implemented into a Python module, adc_gradient. The gradients were used to determine core-excited state optimized geometries and relaxed potential energy surfaces for the water, formic acid, and benzne molecules. <br /

Magnetic circular dichroism within the Algebraic Diagrammatic Construction scheme of the polarisation propagator up to third order

We present an implementation of the B term of Magnetic Circular Dichroism 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 using the Coupled-Cluster Singles and Approximate Doubles (CC2) method, literature TD-DFT results as well as with available experimental data

Characterizing Bonding Patterns in Diradicals and Triradicals by Density-Based Wave Function Analysis: A Uniform Approach

Density-based wave function analysis enables unambiguous comparisons of the electronic structure computed by different methods and removes ambiguity of orbital choices. We use this tool to investigate the performance of different spin-flip methods for several prototypical diradicals and triradicals. In contrast to previous calibration studies that focused on energy gaps between high- and low spin-states, we focus on the properties of the underlying wave functions, such as the number of effectively unpaired electrons. Comparison of different density functional and wave function theory results provides insight into the performance of the different methods when applied to strongly correlated systems such as polyradicals. We show that canonical molecular orbitals for species like large copper-containing diradicals fail to correctly represent the underlying electronic structure due to highly non-Koopmans character, while density-based analysis of the same wave function delivers a clear picture of the bonding pattern

Deciphering the Electronic Transitions of Thiophene-Based Donor-Acceptor-Donor Pentameric Ligands Utilized for Multimodal Fluorescence Microscopy of Protein Aggregates

Anionic pentameric thiophene acetates can be used for fluorescence detection and diagnosis of protein amyloid aggregates. Replacing the central thiophene unit by benzothiadiazole (BTD) or quinoxaline (QX) leads to large emission shifts and basic spectral features have been reported [Chem. Eur. J. 2015, 21, 15133-13137]. Here we present new detailed experimental results of solvent effects, time-resolved fluorescence and examples employing multi-photon microscopy and lifetime imaging. Quantum chemical response calculations elucidate how the introduction of the BTD/QX groups changes the electronic states and emissions. The dramatic red-shift follows an increased conjugation and quinoid character of the pi-electrons of the thiophene backbone. An efficient charge transfer in the excited states S-1 and S-2 compared to the all-thiophene analogue makes these more sensitive to the polarity and quenching by the solvent. Taken together, the results guide in the interpretation of images of stained Alzheimer disease brain sections employing advanced fluorescence microscopy and lifetime imaging, and can aid in optimizing future fluorescent ligand development.Funding Agencies|Swedish Research CouncilSwedish Research Council [2016-00748, 2018-4343]; Swedish e-Science Research Centre, SeRC</p

Benchmarking Post-Hartree–Fock Methods To Describe the Nonlinear Optical Properties of Polymethines: An Investigation of the Accuracy of Algebraic Diagrammatic Construction (ADC) Approaches

Third-order nonlinear optical (NLO) properties of polymethine dyes have been widely studied for applications such as all-optical switching. However, the limited accuracy of the current computational methodologies has prevented a comprehensive understanding of the nature of the lowest excited states and their influence on the molecular optical and NLO properties. Here, attention is paid to the lowest excited-state energies and their energetic ratio, as these characteristics impact the figure-of-merit for all-optical switching. For a series of model polymethines, we compare several algebraic diagrammatic construction (ADC) schemes for the polarization propagator with approximate second-order coupled cluster (CC2) theory, the widely used INDO/MRDCI approach and the symmetry-adapted cluster configuration interaction (SAC-CI) algorithm incorporating singles and doubles linked excitation operators (SAC-CI SD-R). We focus in particular on the ground-to-excited state transition dipole moments and the corresponding state dipole moments, since these quantities are found to be of utmost importance for an effective description of the third-order polarizability γ and two-photon absorption spectra. A sum-overstates expression has been used, which is found to quickly converge. While ADC­(3/2) has been found to be the most appropriate method to calculate these properties, CC2 performs poorly