Multireference Correlated Oscillator Strengths from Adiabatic Connection Approaches Based on Extended Random Phase Approximation.

We show that accurate oscillator strengths can be obtained from adiabatic connection (AC) approaches based on the extended random phase approximation (ERPA) combined with multireference (complete active space, CAS) wave functions. The oscillator strengths calculated using the perturbation-corrected ERPA transition density matrices, proposed in this work, and the excitation energies calculated with recently introduced AC correlation energy methods, AC0 and AC0D, compete with accuracy in the perturbational CASPT2 approach and require less computational effort. AC0 and AC0D methods scale more favorably with the number of active orbitals than multiconfigurational perturbation approaches like CASPT2 and NEVPT2 thanks to their dependence on reduced density matrices up to the order of 2. Importantly, the newly developed approach for computing correlated transition dipole moments does not entail any additional costs, as all intermediate quantities become available when AC0 energies are being computed. We also test the performance of the recently proposed AC method corrected for the negative-transition contributions to the correlation energy, AC0D, for triplet excitation energies. Similarly, as for the singlet excitations, the correction improves the performance of the AC0 method, particularly for the low-lying excited states.

On-Surface Synthesis and Determination of the Open-Shell Singlet Ground State of Tridecacene.

The character of the electronic structure of acenes has been the subject of longstanding discussion. However, convincing experimental evidence of their open-shell character has so far been missing. Here, we present the on-surface synthesis of tridecacene molecules by thermal annealing of octahydrotridecacene on a Au(111) surface. We characterized the electronic structure of the tridecacene by scanning probe microscopy, which reveals the presence of an inelastic signal at 126 meV. We attribute the inelastic signal to spin excitation from the singlet diradical ground state to the triplet excited state. To rationalize the experimental findings, we carried out many-body ab initio calculations as well as model Hamiltonians to take into account the effect of the metallic substrate. Moreover, we provide a detailed analysis of how the dynamic electron correlation and virtual charge fluctuation between the molecule and metallic surface reduces the singlet-triplet band gap. Thus, this work provides the first experimental confirmation of the magnetic character of tridecacene.

Variational Quantum Eigensolver Boosted by Adiabatic Connection.

In this work, we integrate the variational quantum eigensolver (VQE) with the adiabatic connection (AC) method for efficient simulations of chemical problems on near-term quantum computers. Orbital-optimized VQE methods are employed to capture the strong correlation within an active space, and classical AC corrections recover the dynamical correlation effects comprising electrons outside of the active space. On two challenging strongly correlated problems, namely, the dissociation of N2 and the electronic structure of the tetramethyleneethane biradical, we show that the combined VQE-AC approach enhances the performance of VQE dramatically. Moreover, since the AC corrections do not bring any additional requirements on quantum resources or measurements, they can actually boost the VQE algorithms. Our work paves the way toward quantum simulations of real-life problems on near-term quantum computers.

Role of Spin Polarization and Dynamic Correlation in Singlet-Triplet Gap Inversion of Heptazine Derivatives.

The new generation of proposed light-emitting molecules for organic light-emitting diodes (OLEDs) has raised considerable research interest due to its exceptional feature─a negative singlet-triplet (ST) gap violating Hund's multiplicity rule in the excited S1 and T1 states. We investigate the role of spin polarization in the mechanism of ST gap inversion. Spin polarization is associated with doubly excited determinants of certain types, whose presence in the wave function expansion favors the energy of the singlet state more than that of the triplet. Using a perturbation theory-based model for spin polarization, we propose a simple descriptor for prescreening of candidate molecules with negative ST gaps and prove its usefulness for heptazine-type molecules. Numerical results show that the quantitative effect of spin polarization decreases linearly with the increasing highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) exchange integral. Comparison of single- and multireference coupled-cluster predictions of ST gaps shows that the former methods provide good accuracy by correctly balancing the effects of doubly excited determinants and dynamic correlation. We also show that accurate ST gaps may be obtained using a complete active space model supplemented with dynamic correlation from multireference adiabatic connection theory.

Efficient Calculation of the Dispersion Energy for Multireference Systems with Cholesky Decomposition: Application to Excited-State Interactions.

Accurate and efficient prediction of dispersion interactions in excited-state complexes poses a challenge due to the complex nature of electron correlation effects that need to be simultaneously considered. We propose an algorithm for computing the dispersion energy in nondegenerate ground- or excited-state complexes with arbitrary spin. The algorithm scales with the fifth power of the system size due to employing Cholesky decomposition of Coulomb integrals and a recently developed recursive formula for density response functions of the monomers. As a numerical illustration, we apply the new algorithm in the framework of multiconfigurational symmetry adapted perturbation theory, SAPT(MC), to study interactions in dimers with localized excitons. The SAPT(MC) analysis reveals that the dispersion energy may be the main force stabilizing excited-state dimers.

Experimental-computational approach to investigate elastic properties of struvite.

The research presented in the paper concerns the elastic properties of struvite. The article combines theoretical and experimental research. Experimental studies were carried out on struvite single crystals grown in sodium metasilicate gel by single diffusion. This unique method leads to obtaining crystals of sufficiently large size to conduct, for the first time, experimental measurement of elastic properties of monocrystalline struvite. Using the nanoindentation method, the Ez = 29.1 ± 0.7 GPa value of the component of Young's modulus was determined for a struvite single crystal. In addition, the elastic constants C11, C22, and C33 were determined using micro-Brillouin spectroscopy. Theoretical calculations of the abovementioned properties have been carried out by employing density functional theory methods. Scaling of the theoretical elastic constants leads to obtaining good agreement with the experimental values. Values of the Ex and Ey components of the Young's modulus, not available from the experimental nanoindentation technique, have been determined theoretically as 23 GPa and 27 GPa, respectively. Differences in the values of elastic components and Young's modulus components are related to the layered crystal structure of struvite and directional character of the hydrogen-bonding pattern.

TREXIO: A file format and library for quantum chemistry.

TREXIO is an open-source file format and library developed for the storage and manipulation of data produced by quantum chemistry calculations. It is designed with the goal of providing a reliable and efficient method of storing and exchanging wave function parameters and matrix elements, making it an important tool for researchers in the field of quantum chemistry. In this work, we present an overview of the TREXIO file format and library. The library consists of a front-end implemented in the C programming language and two different back-ends: a text back-end and a binary back-end utilizing the hierarchical data format version 5 library, which enables fast read and write operations. It is compatible with a variety of platforms and has interfaces for Fortran, Python, and OCaml programming languages. In addition, a suite of tools have been developed to facilitate the use of the TREXIO format and library, including converters for popular quantum chemistry codes and utilities for validating and manipulating data stored in TREXIO files. The simplicity, versatility, and ease of use of TREXIO make it a valuable resource for researchers working with quantum chemistry data.

Toward more accurate adiabatic connection approach for multireference wavefunctions.

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the "fixed-RDM" approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

Projection-Based Density Matrix Renormalization Group in Density Functional Theory Embedding.

The density matrix renormalization group (DMRG) method has already proved itself as a very efficient and accurate computational method, which can treat large active spaces and capture the major part of strong correlation. Its application on larger molecules is, however, limited by its own computational scaling as well as demands of methods for treatment of the missing dynamical electron correlation. In this work, we present the first step in the direction of combining DMRG with density functional theory (DFT), one of the most employed quantum chemical methods with favorable scaling, by means of the projection-based wave function (WF)-in-DFT embedding. On two proof-of-concept but important molecular examples, we demonstrate that the developed DMRG-in-DFT approach provides a very accurate description of molecules with a strongly correlated fragment.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Interplay between π-Conjugation and Exchange Magnetism in One-Dimensional Porphyrinoid Polymers.

The synthesis of novel polymeric materials with porphyrinoid compounds as key components of the repeating units attracts widespread interest from several scientific fields in view of their extraordinary variety of functional properties with potential applications in a wide range of highly significant technologies. The vast majority of such polymers present a closed-shell ground state, and, only recently, as the result of improved synthetic strategies, the engineering of open-shell porphyrinoid polymers with spin delocalization along the conjugation length has been achieved. Here, we present a combined strategy toward the fabrication of one-dimensional porphyrinoid-based polymers homocoupled via surface-catalyzed [3 + 3] cycloaromatization of isopropyl substituents on Au(111). Scanning tunneling microscopy and noncontact atomic force microscopy describe the thermal-activated intra- and intermolecular oxidative ring closure reactions as well as the controlled tip-induced hydrogen dissociation from the porphyrinoid units. In addition, scanning tunneling spectroscopy measurements, complemented by computational investigations, reveal the open-shell character, that is, the antiferromagnetic singlet ground state (S = 0) of the formed polymers, characterized by singlet-triplet inelastic excitations observed between spins of adjacent porphyrinoid units. Our approach sheds light on the crucial relevance of the π-conjugation in the correlations between spins, while expanding the on-surface synthesis toolbox and opening avenues toward the synthesis of innovative functional nanomaterials with prospects in carbon-based spintronics.

An efficient implementation of time-dependent linear-response theory for strongly orthogonal geminal wave function models.

We present an implementation of time-dependent linear-response equations for strongly orthogonal geminal wave function models: the time-dependent generalized valence bond (TD-GVB) perfect-pairing theory and the antisymmetrized product of strongly orthogonal geminals. The geminal wave functions are optimized using a restricted-step second-order algorithm suitable for handling many geminals, and the linear-response equations are solved in an efficient way using a direct iterative approach. The wave function optimization algorithm features an original scheme to create initial orbitals for the geminal functions in a black-box fashion. The implementation is employed to examine the accuracy of the geminal linear response for singlet excitation energies of small and medium-sized molecules. In systems dominated by dynamic correlation, geminal models constitute only a minor improvement with respect to time-dependent Hartree-Fock. Compared to the linear-response complete active space self-consistent field, TD-GVB either misses or gives large errors for states dominated by double excitations.

Efficient Adiabatic Connection Approach for Strongly Correlated Systems: Application to Singlet-Triplet Gaps of Biradicals.

Strong electron correlation can be captured with multireference wave function methods, but an accurate description of the electronic structure requires accounting for the dynamic correlation, which they miss. In this work, a new approach for the correlation energy based on the adiabatic connection (AC) is proposed. The ACn method accounts for terms up to order n in the coupling constant, and it is size-consistent and free from instabilities. It employs the multireference random phase approximation and the Cholesky decomposition technique, leading to a computational cost growing with the fifth power of the system size. Because of the dependence on only one- and two-electron reduced density matrices, ACn is more efficient than existing ab initio multireference dynamic correlation methods. ACn affords excellent results for singlet-triplet gaps of challenging organic biradicals. The development presented in this work opens new perspectives for accurate calculations of systems with dozens of strongly correlated electrons.

Dispersion Interactions in Exciton-Localized States. Theory and Applications to π-π* and n-π* Excited States.

We address the problem of intermolecular interaction energy calculations in molecular complexes with localized excitons. Our focus is on the correct representation of the dispersion energy. We derive an extended Casimir-Polder formula for direct computation of this contribution through second order in the intermolecular interaction operator V̂. An alternative formula, accurate to infinite order in V̂, is derived within the framework of the adiabatic connection (AC) theory. We also propose a new parametrization of the VV10 nonlocal correlation density functional, so that it corrects the CASSCF energy for the dispersion contribution and can be applied to excited-state complexes. A numerical investigation is carried out for benzene, pyridine, and peptide complexes with the local exciton corresponding to the lowest π-π* or n- π* states. The extended Casimir-Polder formula is implemented in the framework of multiconfigurational symmetry-adapted perturbation theory, SAPT(MC). A SAPT(MC) analysis shows that the creation of a localized exciton affects mostly the electrostatic component of the interaction energy of investigated complexes. Nevertheless, the changes in Pauli repulsion and dispersion energies cannot be neglected. We verify the performance of several perturbation- and AC-based methods. Best results are obtained with a range-separated variant of an approximate AC approach employing extended random phase approximation and CASSCF wave functions.

Dispersion Interactions between Molecules in and out of Equilibrium Geometry: Visualization and Analysis.

The London dispersion interactions between systems undergoing bond breaking, twisting, or compression are not well studied due to the scarcity and the high computational cost of methods being able to describe both the dynamic correlation and the multireference character of the system. Recently developed methods based on the Generalized Valence Bond wave function, such as EERPA-GVB and SAPT(GVB) (SAPT = symmetry-adapted perturbation theory), allow one to accurately compute and analyze noncovalent interactions between multireference systems. Here, we augment this analysis by introducing a local indicator for dispersion interactions inspired by Mata and Wuttke's Dispersion Interaction Density [ J. Comput. Chem. 2017, 38, 15-23] applied on top of an EERPA-GVB computation. Using a few model systems, we show what insights into the nature and evolution of the dispersion interaction during bond breaking and twisting such an approach is able to offer. The new indicator can be used at a minimal cost additional to an EERPA-GVB computation and can be complemented by an energy decomposition employing the SAPT(GVB) method. We explain the physics behind the initial increase, followed by a decrease in the interaction of linear molecules upon bond stretching. Namely, the elongation of covalent bonds leads to the enhancement of attractive dispersion interactions. For even larger bond lengths, this effect is canceled by the increase of the repulsive exchange forces resulting in a suppression of the interaction and finally leading to repulsion between monomers.

Density Matrix Renormalization Group with Dynamical Correlation via Adiabatic Connection.

The quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as one of the methods of choice for calculations of strongly correlated molecular systems. Despite its great ability to capture strong electronic correlation in large active spaces, it is not suitable for computations of dynamical electron correlation. In this work, we present a new approach to the electronic structure problem of strongly correlated molecules, in which DMRG is responsible for a proper description of the strong correlation, whereas dynamical correlation is computed via the recently developed adiabatic connection (AC) technique which requires only up to two-body active space reduced density matrices. We report the encouraging results of this approach on typical candidates for DMRG computations, namely, n-acenes (n = 2 → 7), Fe(II)-porphyrin, and the Fe3S4 cluster.

Symmetry-Adapted Perturbation Theory Based on Multiconfigurational Wave Function Description of Monomers.

We present a formulation of the multiconfigurational (MC) wave function symmetry-adapted perturbation theory (SAPT). The method is applicable to noncovalent interactions between monomers which require a multiconfigurational description, in particular when the interacting system is strongly correlated or in an electronically excited state. SAPT(MC) is based on one- and two-particle reduced density matrices of the monomers and assumes the single-exchange approximation for the exchange energy contributions. Second-order terms are expressed through response properties from extended random phase approximation (ERPA). The dispersion components of SAPT(MC) have been introduced in our previous works [Hapka, M. J. Chem. Theory Comput. 2019, 15, 1016-1027; Hapka, M. J. Chem. Theory Comput. 2019, 15, 6712-6723]. SAPT(MC) is applied either with generalized valence bond perfect pairing (GVB) or with complete active space self-consistent field (CASSCF) treatment of the monomers. We discuss two model multireference systems: the H2 ··· H2 dimer in out-of-equilibrium geometries and interaction between the argon atom and excited state of ethylene. Using the C2H4* ··· Ar complex as an example, we examine second-order terms arising from negative transitions in the linear response function of an excited monomer. We demonstrate that the negative-transition terms must be accounted for to ensure qualitative prediction of induction and dispersion energies and develop a procedure allowing for their computation. Factors limiting the accuracy of SAPT(MC) are discussed in comparison with other second-order SAPT schemes on a data set of small single-reference dimers.

Excited states in the adiabatic connection fluctuation-dissipation theory: Recovering missing correlation energy from the negative part of the density response spectrum.

The adiabatic connection (AC) theory offers an alternative to the perturbation theory methods for computing correlation energy in the multireference wavefunction framework. We show that the AC correlation energy formula can be expressed in terms of the density linear response function as a sum of components related to positive and negative parts of the transition energy spectrum. Consequently, generalization of the adiabatic connection fluctuation-dissipation theory to electronically excited states is obtained. The component of the linear response function related to the negative-transition energy enters the correlation energy expression with an opposite sign to that of the positive-transition part and is non-negligible in the description of excited states. To illustrate this, we analyze the approximate AC model in which the linear response function is obtained in the extended random phase approximation (ERPA). We demonstrate that AC can be successfully combined with the ERPA for excited states, provided that the negative-excitation component of the response function is rigorously accounted for. The resulting AC0D model, an extension of the AC0 scheme introduced in our earlier works, is applied to a benchmark set of singlet excitation energies of organic molecules. AC0D constitutes a significant improvement over AC0 by bringing the excitation energies of the lowest excited states to a satisfactory agreement with theoretical best estimates, which parallels or even exceeds the accuracy of the n-electron valence state perturbation theory method. For higher excitations, AC0D is less reliable due to the gradual deterioration of the underlying ERPA linear response.

Extension of an Atom-Atom Dispersion Function to Halogen Bonds and Its Use for Rational Design of Drugs and Biocatalysts.

A dispersion function Das in the form of a damped atom-atom asymptotic expansion fitted to ab initio dispersion energies from symmetry-adapted perturbation theory was improved and extended to systems containing heavier halogen atoms. To illustrate its performance, the revised Das function was implemented in the multipole first-order electrostatic and second-order dispersion (MED) scoring model. The extension has allowed applications to a much larger set of biocomplexes than it was possible with the original Das. A reasonable correlation between MED and experimentally determined inhibitory activities was achieved in a number of test cases, including structures featuring nonphysically shortened intermonomer distances, which constitute a particular challenge for binding strength predictions. Since the MED model is also computationally efficient, it can be used for reliable and rapid assessment of the ligand affinity or multidimensional scanning of amino acid side-chain conformations in the process of rational design of novel drugs or biocatalysts.