Search

In recent years there has been a rapid growth in the development and application of new stochastic methods in electronic structure. These methods are quite diverse, from many-body wave function techniques in real space or determinant space to being used to sum perturbative expansions. This growth has been spurred by the more favorable scaling with the number of electrons and often better parallelization over large numbers of central processing unit (CPU) cores or graphical processing units (GPUs) than for high-end non-stochastic wave function based methods. This special issue of the Journal of Chemical Physics includes 33 papers that describe recent developments and applications in this area. As seen from the articles in the issue, stochastic electronic structure methods are applicable to both molecules and solids and can accurately describe systems with strong electron correlation. This issue was motivated, in part, by the 2019 Telluride Science Research Center workshop on Stochastic Electronic Structure Methods that we organized. Below we briefly describe each of the papers in the special issue, dividing the papers into six subtopics. [ABSTRACT FROM AUTHOR]

In Paper I, the performances of pre-screening (PS), extended PS (EPS), and cumulant (CU) approximations to the fourth-order density matrix were examined in the context of second-order N-electron valence state perturbation theory (NEVPT2). It has been found that the CU, PS, and even EPS approximations with loose thresholds may introduce intruder states. In the present work, the origin of these "false intruder" states introduced by approximated density matrices is discussed. Canonical NEVPT2 implementations employ a rank reduction trick. By analyzing its residual error, we find that the omission of the rank reduction leads to a more stable multireference perturbation theory for incomplete active space reference wave functions. Such a full rank (FR)-NEVPT2 formulation is equivalent to the conventional NEVPT2 method for the complete active space self-consistent field/complete active space configuration interaction reference wave function. A major drawback of the FR-NEVPT2 formulation is the necessity of the fifth-order density matrix. To avoid the construction of the high-order density matrices, the combination of the FR-NEVPT2 with the CU approximation is studied. However, we find that the CU approximation remains problematic as it still introduces intruder states. The question of how to robustly and efficiently perform internally contracted multireference perturbation theories with approximate densities remains a challenging field of investigation. [ABSTRACT FROM AUTHOR]

The present paper reformulates and improves a previously proposed determinant-based second-order multi-reference perturbative formalism. Through a rather simple modification of the energy denominators, this formalism takes into account the interactions between the model space determinants, which are repeated in outer space. The method has been shown to be size-consistent when the model space is a complete active space, which is a severe limit. It is shown here that the completeness of the model space is not necessary to keep this property, provided that the zero-order function satisfies some conditions. For instance, size consistency may be obtained from truncated complete active spaces. It may even be satisfied from Singles and Doubles Configuration Interactions, provided that a coupled electron pair approximation is used in the definition of the model space wave function. The physical content of the method is illustrated by a series of model problems, showing its robustness. A major benefit of the fact that the perturbers are single determinants is the possibility to revise with full flexibility the model-space component of the wave function, i.e., to treat the feedback effect of the dynamic correlation on the valence component of the wave function. [ABSTRACT FROM AUTHOR]

In this paper, the Hessian matrix of valence bond (VB) self-consistent field (VBSCF) energy with respect to orbitals are evaluated by applying the nonorthogonal orbital based N-body reduced density matrices, which was presented in Paper I. To this end, an automatic formula/code generator (AFCG) is developed; with which the matrix elements between internally contracted excited configurations of VB wave function and the corresponding codes are generated automatically. Compared to the tedious manual formula deducing and implementing, AFCG is much more convenient and efficient, and enables us to avoid troublesome debugging. With the help of AFCG, the Hessian-based Newton-Raphson algorithm is implemented for the VBSCF orbital optimization. Test calculations indicate that the Newton-Raphson algorithm converges quadratically and has much better convergence behavior than the gradient-based LBFGS algorithms. Furthermore, a combined approach with LBFGS and Newton-Raphson algorithms is applied to reduce the total CPU time of the calculation. [ABSTRACT FROM AUTHOR]

Allowing triplet components of individual geminals, spin-contaminated strongly orthogonal geminal wave functions may emerge, which can be ameliorated by spin-projection techniques. Of the latter, half-projection was previously shown to be useful, offering a compromise between the amount of remaining spin-contamination and the violation of size consistency generated by projection. This paper investigates how a half-projected spin-contaminated geminal wave function can be improved by multi-configuration perturbation theory to incorporate dynamical correlation effects. [ABSTRACT FROM AUTHOR]

The single-reference coupled-cluster method has proven very effective in the ab initio description of atomic and molecular systems, but its successful application is limited to states dominated by a single Slater determinant, which is used as the reference. In cases where several determinants are important in the wave function expansion, i.e., we have to deal with nondynamic correlation effects, a multi-reference version of the coupled-cluster method is required. The multi-reference coupled-cluster approaches are based on the effective Hamiltonian formulation providing a two-step procedure, in which dynamic correlation effects can be efficiently evaluated by the wave operator, while nondynamic correlation contributions are given by diagonalization of the effective Hamiltonian in the final step. There are two classical multi-reference coupled-cluster formulations. In this paper, the focus is on the so-called Fock-space coupled-cluster method in its basic version with one- and two-particle operators in the exponent. Computational schemes using this truncation of the cluster operator have been successfully applied in calculations in one- and two-valence sectors of the Fock space. In this paper, we show that the approach can be easily extended and effectively employed in the three-valence sector calculations. [ABSTRACT FROM AUTHOR]

This paper combines the valence bond block diabatization approach with the idea of orbital breathing. With highly compact wave functions, the breathing orbital valence bond (BOVB) method is applied to investigate several atomic and molecular properties, including the electron affinity of F, the adiabatic and diabatic potential energy curves and the dipole moment curves of the two lowest-lying 1Σ+ states, the electronic coupling curve and the crossing distance of the two diabatic states, and the spectroscopic constants of the ground states for LiF. The configuration selection scheme proposed in this work is quite general, requiring only the selection of several de-excitation and excitation orbitals in a sense like the restricted active space self-consistent field method. Practically, this is also the first time that BOVB results are extrapolated to complete basis set limit. Armed with the chemical intuition provided by valence bond theory, the classic but challenging covalent-ionic interaction in the title molecule is not only conceptually interpreted but also accurately computed. [ABSTRACT FROM AUTHOR]

We present a methodology for analyzing chemical bonds embedded in the electronic wave function of molecules, especially in terms of spin correlations or so-called "local spin." In this paper, based on biorthogonal second quantization, the spin correlation functions of molecules are naturally introduced, which enables us to extract local singlet and local triplet elements from the wave function. We also clarify the relationship between these spin correlations and traditional chemical concepts, i.e., resonance structures. Several chemical reactions, including the intramolecular radical cyclization and the formation of preoxetane, are demonstrated to verify the analysis method numerically. [ABSTRACT FROM AUTHOR]

C4n cyclacenes exhibit strong bond-alternation in their equilibrium geometry. In the two equivalent geometries, the system keeps an essentially closed-shell character. The two energy minima are separated by a transition state suppressing the bond-alternation, where the wave function is strongly diradical. This paper discusses the physical factors involved in this energy difference and possible evaluations of the barrier height. The barrier given as the energy difference between the restricted density functional theory (DFT)/B3LYP for the equilibrium and the broken symmetry DFT/B3LYP of the transition state is either negative or small, in contradiction with the most reliable Wave Function Theory calculations. The minimal (two electrons in two molecular orbitals) Complete Active Space self-consistent field (CASSCF) overestimates the barrier, and the subsequent second-order perturbation cancels it. Due to the collective character of the spin-polarization effect, it is necessary to perform a full π CASSCF + second-order perturbation to reach a reasonable value of the barrier, but this type of treatment cannot be applied to large molecules. DFT procedures treating on an equal foot the closed-shell and open-shell geometries have been explored, such as Mixed-Reference Spin-Flip Time-dependent-DFT and a new spin-decontamination proposal, namely, DFT-dressed configuration interaction, but the results still depend on the density functional. M06-2X without or with spin-decontamination gives the best agreement with the accurate wave function results. [ABSTRACT FROM AUTHOR]

With a He– H 3 + interaction potential obtained from advanced electronic structure calculations, we computed the vibration-rotation-tunneling (VRT) states of this complex for total angular momenta J from 0 to 9, both for the vibrational ground state and for the twofold degenerate v2 = 1 excited state of H 3 + . The potential has three equivalent global minima with depth De = 455.3 cm−1 for He in the plane of H 3 + , three equatorial saddle points that separate these minima with barriers of 159.5 cm−1, and two axial saddle points with energies of 243.1 cm−1 above the minima. The dissociation energies calculated for the complexes of He with ortho- H 3 + (o H 3 + ) and para- H 3 + (p H 3 + ) are D0 = 234.5 and 236.3 cm−1, respectively. Wave function plots of the VRT states show that they may be characterized as weakly hindered internal rotor states, delocalized over the three minima in the potential and with considerable amplitude at the barriers. Most of them are dominated by the jk = 10 and 11 rotational ground states of o H 3 + and p H 3 + , with the intermolecular stretching mode excited up to v = 4 inclusive. However, we also found excited internal rotor states: 33 in He–o H 3 + , and 22 and 21 in He–p H 3 + . The VRT levels and wave functions were used to calculate the frequencies and line strengths of all allowed v2 = 0 → 1 rovibrational transitions in the complex. Theoretical spectra generated with these results are compared with the experimental spectra in Paper II [Salomon et al., J. Chem. Phys. 156, 144308 (2022)] and are extremely helpful in assigning these spectra. This comparison shows that the theoretical energy levels and spectra agree very well with the measured ones, which confirms the high accuracy of our ab initio He– H 3 + interaction potential and of the ensuing calculations of the VRT states. [ABSTRACT FROM AUTHOR]

Downfolding coupled cluster techniques have recently been introduced into quantum chemistry as a tool for the dimensionality reduction of the many-body quantum problem. As opposed to earlier formulations in physics and chemistry based on the concept of effective Hamiltonians, the appearance of the downfolded Hamiltonians is a natural consequence of the single-reference exponential parameterization of the wave function. In this paper, we discuss the impact of higher-order terms originating in double commutators. In analogy to previous studies, we consider the case when only one- and two-body interactions are included in the downfolded Hamiltonians. We demonstrate the efficiency of the many-body expansions involving single and double commutators for the unitary extension of the downfolded Hamiltonians on the example of the beryllium atom, and bond-breaking processes in the Li2 and H2O molecules. For the H2O system, we also analyze energies obtained with downfolding procedures as functions of the active space size. [ABSTRACT FROM AUTHOR]

The chemical formula theory (CFT) proposed in Paper I of this series [H. Nakatsuji et al., J. Chem. Phys. 149, 114105 (2018)] is a simple variational electronic structure theory for atoms and molecules. The CFT constructs simple, conceptually useful wave functions for the ground and excited states, simultaneously, from the ground and excited states of the constituent atoms, reflecting the spirits of the chemical formulas. The CFT wave functions are also designed to be used as the initial wave functions of the free complement (FC) theory, that is, the exact theory producing the exact wave functions of the Schrödinger accuracy. This combined theory is referred to as the FC-CFT. We aim to construct an exact wave function theory that is useful not only quantitatively but also conceptually. This paper shows the atomic applications of the CFT and the FC-CFT. For simplicity, we choose the small atoms, Be and Li, and perform variational calculations to essentially exact levels. For these elements, a simple Hylleraas CI type formulation is known to be potentially highly accurate: we realize it with the CFT and the FC-CFT. Even from the CFT levels, the excitation energies to the Rydberg excited states were calculated satisfactorily. Then, with increasing the order of the FC theory in the FC-CFT, all the absolute energies and the excitation energies of the Be and Li atoms were improved uniformly and reached rapidly to the essentially exact levels in order 3 or 4 with moderately small calculational labors. [ABSTRACT FROM AUTHOR]

The quest for "chemical accuracy" is becoming more and more demanded in the field of structure and kinetics of molecules at solid surfaces. In this paper, as an example, we focus on the barrier for hydrogen diffusion on a α-Al2O3(0001) surface, aiming for a couple cluster singles, doubles, and perturbative triples [CCSD(T)]-level benchmark. We employ the density functional theory (DFT) optimized minimum and transition state structures reported by Heiden, Usvyat, and Saalfrank [J. Phys. Chem. C 123, 6675 (2019)]. The barrier is first evaluated at the periodic Hartree–Fock and local Møller–Plesset second-order perturbation (MP2) level of theory. The possible sources of errors are then analyzed, which includes basis set incompleteness error, frozen core, density fitting, local approximation errors, as well as the MP2 method error. Using periodic and embedded fragment models, corrections to these errors are evaluated. In particular, two corrections are found to be non-negligible (both from the chemical accuracy perspective and at the scale of the barrier value of 0.72 eV): the correction to the frozen core-approximation of 0.06 eV and the CCSD(T) correction of 0.07 eV. Our correlated wave function results are compared to barriers obtained from DFT. Among the tested DFT functionals, the best performing for this barrier is B3LYP-D3. [ABSTRACT FROM AUTHOR]

In this paper, we report on a correctly scaling novel coupled cluster singles and doubles (CCSD) implementation for arbitrary high-spin open-shell states. The chosen cluster operator is completely spin-free, i.e., employs spatial substitutions only. It is composed of our recently developed Löwdin-type operators [N. Herrmann and M. Hanrath, J. Chem. Phys. 153, 164114 (2020)], which ensure (1) spin completeness and (2) spin adaption, i.e., spin purity of the CC wave function. In contrast to the proof-of-concept matrix-representation-based implementation presented there, the present implementation relies on second quantization and factorized tensor contractions. The generated singles and doubles operators are embedded in an equation generation engine. In the latter, Wick's theorem is used to derive prefactors arising from spin integration directly from the spin-free full contraction patterns. The obtained Wick terms composed of products of Kronecker deltas are represented by special non-antisymmetrized Goldstone diagrams. Identical (redundant) diagrams are identified by solving the underlying graph isomorphism problem. All non-redundant graphs are then automatically translated to locally—one term at a time—factorized tensor contractions. Finally, the spin-adapted and spin-complete (SASC) CCS and CCSD variants are applied to a set of small molecular test systems. Both correlation energies and amplitude norms hint toward a reasonable convergence of the SASC-CCSD method for a Baker–Campbell–Hausdorff series truncation of order four. In comparison to spin orbital CCSD, SASC-CCSD leads to slightly improved correlation energies with differences of up to 1.292mEH (1.10% with respect to full configuration identification) for quintet CH2 in the cc-pVDZ basis set. [ABSTRACT FROM AUTHOR]

Quantum computation promises to provide substantial speedups in many practical applications with a particularly exciting one being the simulation of quantum many-body systems. Adiabatic state preparation (ASP) is one way that quantum computers could recreate and simulate the ground state of a physical system. In this paper, we explore a novel approach for classically simulating the time dynamics of ASP with high accuracy and with only modest computational resources via an adaptive sampling configuration interaction scheme for truncating the Hilbert space to only the most important determinants. We verify that this truncation introduces negligible error and use this new approach to simulate ASP for sets of small molecular systems and Hubbard models. Furthermore, we examine two approaches to speeding up ASP when performed on quantum hardware: (i) using the complete active space configuration interaction (CASCI) wave function instead of the Hartree–Fock initial state and (ii) a nonlinear interpolation between the initial and target Hamiltonians. We find that starting with a CASCI wave function with a limited active space yields substantial speedups for many of the systems examined, while nonlinear interpolation does not. In additional, we observe interesting trends in the minimum gap location (based on the initial state) as well as how state preparation time can depend on certain molecular properties, such as the number of valence electrons. Importantly, we find that the required state preparation times do not show an immediate exponential wall that would preclude an efficient run of ASP on actual hardware. [ABSTRACT FROM AUTHOR]

The computation of the nuclear quantum dynamics of molecules is challenging, requiring both accuracy and efficiency to be applicable to systems of interest. Recently, theories have been developed for employing time-dependent basis functions (denoted modals) with vibrational coupled cluster theory (TDMVCC). The TDMVCC method was introduced along with a pilot implementation, which illustrated good accuracy in benchmark computations. In this paper, we report an efficient implementation of TDMVCC, covering the case where the wave function and Hamiltonian contain up to two-mode couplings. After a careful regrouping of terms, the wave function can be propagated with a cubic computational scaling with respect to the number of degrees of freedom. We discuss the use of a restricted set of active one-mode basis functions for each mode, as well as two interesting limits: (i) the use of a full active basis where the variational modal determination amounts essentially to the variational determination of a time-dependent reference state for the cluster expansion; and (ii) the use of a single function as an active basis for some degrees of freedom. The latter case defines a hybrid TDMVCC/TDH (time-dependent Hartree) approach that can obtain even lower computational scaling. The resulting computational scaling for hybrid and full TDMVCC[2] is illustrated for polyaromatic hydrocarbons with up to 264 modes. Finally, computations on the internal vibrational redistribution of benzoic acid (39 modes) are used to show the faster convergence of TDMVCC/TDH hybrid computations towards TDMVCC compared to simple neglect of some degrees of freedom. [ABSTRACT FROM AUTHOR]

Within the well-established optical response function formalism, a new strategy with the central idea of employing the forward-backward stochastic Schrödinger equations in a segmented way to accurately obtain the two-dimensional (2D) electronic spectrum is presented in this paper. Based on the simple excitonically coupled dimer model system, the validity and efficiency of the proposed schemes are demonstrated in detail, along with the comparison against the deterministic hierarchy equations of motion and perturbative second-order time-convolutionless quantum master equations. In addition, an important insight is provided in this paper that the characteristic frequency of the overdamped environment is an extremely crucial factor to regulate the lifetimes of the oscillating signals in 2D electronic spectra and of quantum coherence features of system dynamics. It is worth noting that the proposed scheme benefiting from its stochastic nature and wavefunction framework and many other advantages of substantially reducing the numerical cost has a great potential to systematically investigate various quantum effects observed in realistic large-scale natural and artificial photosynthetic systems. [ABSTRACT FROM AUTHOR]

X-ray photon absorption leads to the creation of highly excited species, which often decay through the Auger process. The theoretical treatment of Auger decay is challenging because of the resonance nature of the initial core-excited or core-ionized states and the continuous nature of the ejected electron. In Paper I [W. Skomorowski and A. I. Krylov, J. Chem. Phys. 154, 084124 (2021)], we have introduced a theoretical framework for computing Auger rates based on the Feshbach–Fano approach and the equation-of-motion coupled-cluster ansätze augmented with core–valence separation. The outgoing Auger electron is described with a continuum orbital. We considered two approximate descriptions—a plane wave and a Coulomb wave with an effective charge. Here, we use the developed methodology to calculate Auger transition rates in core-ionized and core-excited benchmark systems (Ne, H2O, CH4, and CO2). Comparison with the available experimental spectra shows that the proposed computational scheme provides reliable ab initio predictions of the Auger spectra. The reliability, cost efficiency, and robust computational setup of this methodology offer advantages in applications to a large variety of systems. [ABSTRACT FROM AUTHOR]

In the full optimized reaction space and valence-complete active space self-consistent field (vCAS) methods, a set of active orbitals is defined as the union of the valence orbitals on the atoms, all possible configurations involving the active orbitals are generated, and the orbitals and configuration coefficients are self-consistently optimized. Such wave functions have tremendous flexibility, which makes these methods incredibly powerful but can also lead to inconsistencies in the description of the electronic structure of molecules. In this paper, the problems that can arise in vCAS calculations are illustrated by calculations on the BH and BF molecules. BH is well described by the full vCAS wave function, which accounts for molecular dissociation and 2s–2p near-degeneracy in the boron atom. The same is not true for the full vCAS wave function for BF. There is mixing of core and active orbitals at short internuclear distances and swapping of core and active orbitals at large internuclear distances. In addition, the virtual 2π orbitals, which were included in the active space to account for the 2s–2p near degeneracy effect, are used instead to describe radial correlation of the electrons in the F2pπ-like pairs. Although the above changes lead to lower vCAS energies, they lead to higher vCAS+1+2 energies as well as irregularities and/or discontinuities in the potential energy curves. All of the above problems can be addressed by using the spin-coupled generalized valence bond-inspired vCAS wave function for BF, which includes only a subset of the atomic valence orbitals in the active space. [ABSTRACT FROM AUTHOR]

This paper presents the first implementation of the row-orthonormal hyperspherical coordinate formalism for the computation of the vibrational spectrum of a tetratomic system. The wavefunction of Ne4 is expanded on a large basis set of hyperspherical harmonics generated numerically. This method not only provides spectra with reasonable accuracy, but also gives physical insight into the vibrational dynamics of the system. The characteristics of the spectra are related to the symmetry and localization of the wavefunction in configuration space. [ABSTRACT FROM AUTHOR]

Quantum entangled photons provide a sensitive probe of many-body interactions and offer a unique experimental portal for quantifying many-body correlations in a material system. In this paper, we present a theoretical demonstration of how photon–photon entanglement can be generated via interactions between coupled qubits. Here, we develop a model for the scattering of an entangled pair of photons from a molecular dimer. We develop a diagrammatic theory for the scattering matrix and show that one can correlate the von Neumann entropy of the outgoing bi-photon wave function with exciton exchange and repulsion interactions. We conclude by discussing possible experimental scenarios for realizing these ideas. [ABSTRACT FROM AUTHOR]

Fermi's golden rule defines the transition rate between weakly coupled states and can thus be used to describe a multitude of molecular processes including electron-transfer reactions and light-matter interaction. However, it can only be calculated if the wave functions of all internal states are known, which is typically not the case in molecular systems. Marcus theory provides a closed-form expression for the rate constant, which is a classical limit of the golden rule, and indicates the existence of a normal regime and an inverted regime. Semiclassical instanton theory presents a more accurate approximation to the golden-rule rate including nuclear quantum effects such as tunneling, which has so far been applicable to complex anharmonic systems in the normal regime only. In this paper, we extend the instanton method to the inverted regime and study the properties of the periodic orbit, which describes the tunneling mechanism via two imaginary-time trajectories, one of which now travels in negative imaginary time. It is known that tunneling is particularly prevalent in the inverted regime, even at room temperature, and thus, this method is expected to be useful in studying a wide range of molecular transitions occurring in this regime. [ABSTRACT FROM AUTHOR]

The exponential parametrization of the wave function used in the coupled-cluster approaches has proven very successful in the ab initio description of atomic and molecular systems. This concerns first of all the single-reference version of the method that is designed for states dominated by a single Slater determinant. Usually, the coupled-cluster methods with one- and two-body excitation operators in the exponent form the basic computational schemes. The inclusion of three-body effects in the cluster operator to increase the accuracy of the results is numerically expensive, so their approximate evaluation is rather used in practice. In the case of the single-reference coupled-cluster approach, the problem of approximate evaluation of three-body effects in the cluster operator has been well studied, and computational schemes of both noniterative and iterative nature have been proposed. The situation is different in the case of multireference coupled-cluster methods which are required to describe open shell and quasidegenerate states. The multireference approaches in their standard effective Hamiltonian formulations are more complicated and less frequently used in routine calculations; however, one of them, the so-called Fock-space coupled-cluster method, becomes very effective if reformulated within the intermediate Hamiltonian framework. Both the basic version of the method with one- and two-body clusters and the extended one that includes up to three-body operators in the exponent are implemented. The latter approach provides more accurate results, but its relatively high numerical cost limits its applicability. For this reason, going beyond the basic scheme with one- and two-body clusters through an approximate evaluation of the impact of three-body clusters is of great interest. In the paper, we investigate different ways of approximate inclusion of the three-body effects in the Fock-space coupled-cluster method designated for excitation energy calculations. [ABSTRACT FROM AUTHOR]

In this paper, we analyze new approximations of the Green's function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i(n, m) method, where the excitation level of the inner auxiliary operators (m) used to describe the ionization potential and electron affinity effects in the N − 1 and N + 1 particle spaces is higher than the excitation level (n) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called "n + 1" rule in this category [or the GFCC-i(n, n + 1) method], which states that in order to maintain size-extensivity of the Green's function matrix elements, the excitation level of inner auxiliary operators Xp(ω) and Yq(ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N2 and CO molecules in the inner and outer valence regimes. In comparison with the Green's function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions. [ABSTRACT FROM AUTHOR]

Probabilities to find a chosen number of electrons in flexible domains of space are calculated for highly correlated wave functions. Quantum mechanics can produce higher probabilities for chemically relevant arrangements of electrons in these regions. However, the probability to have a given arrangement, e.g., that corresponding to chemical formulas (bonds or atoms), is low although being often maximal. Like in valence bond theory, it is useful to consider alternative distributions of electrons. Exchanges of electrons should be considered not only between atoms but also between other types of regions, such as those attributed to lone pairs. It is useful to have definitions flexible enough to allow users to find the most relevant representations. We tentatively suggest a tool (the effective number of parties) to help one make the choice. [ABSTRACT FROM AUTHOR]

For a quantum system subject to a time-dependent perturbing field, Dirac's analysis gives the probability of transition to an excited state |k⟩ in terms of the norm square of the entire excited-state coefficient ck(t) in the wave function. By integrating by parts in Dirac's equation for ck(t) at first order, Landau and Lifshitz separated ck(1)(t) into an adiabatic term ak(1)(t) that characterizes the gradual adjustment of the ground state to the perturbation without transitions and a nonadiabatic term bk(1)(t) that depends explicitly on the time derivative of the perturbation at times t′ ≤ t. Landau and Lifshitz stated that the probability of transition in a pulsed perturbation is given by |bk(t)|2, rather than by |ck(t)|2. We use the term "transition probability" to refer to the probability that a true excited-state component is present in the time-evolved wave function, as opposed to a smooth modification of the initial state. In recent work, we have examined the differences between |bk(t)|2 and |ck(t)|2 when a system is perturbed by a harmonic wave in a Gaussian envelope. We showed that significant differences exist when the frequency of the harmonic wave is off-resonance with the transition frequency. In this paper, we consider Gaussian perturbations and pulses that rise via a half Gaussian shoulder to a level plateau and later return to zero via a down-going half Gaussian. While the perturbation is constant, the transition probability |bk(t)|2 does not change. By contrast, |ck(t)|2 continues to oscillate while the perturbation is constant, and its time averaged value differs from |bk(t)|2. We suggest a general type of experiment to prove that the transition probability is given by |bk(t)|2, not |ck(t)|2. We propose a ratio test that does not require accurate knowledge of transition matrix elements or absolute field intensities. [ABSTRACT FROM AUTHOR]

Recently it has been shown how to employ the adiabatic connection (AC) formalism to obtain correlation energy for multireference wavefunctions [K. Pernal, Phys. Rev. Lett. 120, 013001 (2018)]. Approximations to the exact AC formulation have been based on assuming that a one-electron reduced density matrix is constant along the AC path and by employing the extended random phase approximation. In this paper, the importance of these approximations is examined by comparing approximate AC integrands with their exact counterparts obtained for the hydrogen molecule in its ground and excited states. Encouraging results obtained for H2 indicate that AC is a viable and promising approach to a correlation energy problem not only for ground but also for excited states of electronic systems. [ABSTRACT FROM AUTHOR]

Methods which aim at universal applicability must be able to describe both weak and strong electronic correlation with equal facility. Such methods are in short supply. The combination of symmetry projection for strong correlation and coupled cluster theory for weak correlation offers tantalizing promise to account for both on an equal footing. In order to do so, however, the coupled cluster portion of the wave function must be optimized in the presence of the symmetry projection. This paper discusses how this may be accomplished, and shows the importance of doing so for both the Hubbard model Hamiltonian and the molecular Hamiltonian, all with a computational scaling comparable to that of traditional coupled cluster theory. [ABSTRACT FROM AUTHOR]

We present the theory and algorithms for computing analytical energy gradients for explicitly correlated second-order Møller-Plesset perturbation theory (MP2-F12). The main difficulty in F12 gradient theory arises from the large number of two-electron integrals for which effective two-body density matrices and integral derivatives need to be calculated. For efficiency, the density fitting approximation is used for evaluating all two-electron integrals and their derivatives. The accuracies of various previously proposed MP2-F12 approximations [3C, 3C(HY1), 3*C(HY1), and 3*A] are demonstrated by computing equilibrium geometries for a set of molecules containing first- and second-row elements, using double-ζ to quintuple-ζ basis sets. Generally, the convergence of the bond lengths and angles with respect to the basis set size is strongly improved by the F12 treatment, and augmented triple-ζ basis sets are sufficient to closely approach the basis set limit. The results obtained with the different approximations differ only very slightly. This paper is the first step towards analytical gradients for coupled-cluster singles and doubles with perturbative treatment of triple excitations, which will be presented in the second part of this series. [ABSTRACT FROM AUTHOR]

The Coupled-Cluster expansion, truncated after single and double excitations (CCSD), provides accurate and reliable molecular electronicwave functions and energies for many molecular systems around their equilibrium geometries. However, the high computational cost, which is well-known to scale as O(N6) with system size N, has limited its practical application to small systems consisting of not more than approximately 20-30 atoms. To overcome these limitations, low-order scaling approximations to CCSD have been intensively investigated over the past fewyears. In our previouswork, we have shown that by combining the pair natural orbital (PNO) approach and the concept of orbital domains it is possible to achieve fully linear scaling CC implementations (DLPNO-CCSD and DLPNO-CCSD(T)) that recover around 99.9% of the total correlation energy [C. Riplinger et al., J. Chem. Phys. 144, 024109 (2016)]. The production level implementations of the DLPNO-CCSD and DLPNO-CCSD(T) methods were shown to be applicable to realistic systems composed of a few hundred atoms in a routine, black-box fashion on relatively modest hardware. In 2011, a reduced-scaling CCSD approach for high-spin open-shell unrestricted Hartree-Fock reference wave functions was proposed (UHF-LPNOCCSD) [A. Hansen et al., J. Chem. Phys. 135, 214102 (2011)]. After a few years of experience with this method, a few shortcomings of UHF-LPNO-CCSD were noticed that required a redesign of the method, which is the subject of this paper. To this end, we employ the high-spin open-shell variant of the N-electron valence perturbation theory formalism to define the initial guess wave function, and consequently also the open-shell PNOs. The new PNO ansatz properly converges to the closed-shell limit since all truncations and approximations have been made in strict analogy to the closed-shell case. Furthermore, given the fact that the formalism uses a single set of orbitals, only a single PNO integral transformation is necessary, which offers large computational savings. We show that, with the default PNO truncation parameters, approximately 99.9% of the total CCSD correlation energy is recovered for open-shell species, which is comparable to the performance of the method for closedshells. UHF-DLPNO-CCSD shows a linear scaling behavior for closed-shell systems, while linear to quadratic scaling is obtained for open-shell systems. The largest systems we have considered contain more than 500 atoms and feature more than 10 000 basis functions with a triple-ζ quality basis set. [ABSTRACT FROM AUTHOR]

The impact of the antisymmetrization is often addressed as a local property of the many-electron wave function, namely that the wave function should vanish when two electrons with parallel spins are in the same position in space. In this paper, we emphasize that this presentation is unduly restrictive: we illustrate the strong non-local character of the antisymmetrization principle, together with the fact that it is a matter of spin symmetry rather than spin parallelism. To this aim, we focus our attention on the simplest representation of various states of two-electron systems, both in atomic (helium atom) and molecular (H2 and the π system of the ethylene molecule) cases. We discuss the non-local property of the nodal structure of some two-electron wave functions, both using analytical derivations and graphical representations of cuttings of the nodal hypersurfaces. The attention is then focussed on the impact of the antisymmetrization on the maxima of the two-body density, and we show that it introduces strong correlation effects (radial and/or angular) with a non-local character. These correlation effects are analyzed in terms of inflation and depletion zones, which are easily identifiable, thanks to the nodes of the orbitals composing the wave function. Also, we show that the correlation effects induced by the antisymmetrization occur also for anti-parallel spins since all Ms components of a given spin state have the same N-body densities. Finally, we illustrate that these correlation effects occur also for the singlet states, but they have strictly opposite impacts: the inflation zones in the triplet become depletion zones in the singlet and vice versa. [ABSTRACT FROM AUTHOR]

Double-hybrid density functional approximations (DH-DFAs) provide an accurate description of the electronic structure of molecules by semiempirically mixing density functional and wavefunction theory. In this paper, we investigate the properties of the potential used in such approximations. By using the optimized effective potential approach, the consistent Kohn-Sham (KS) potential for a double-hybrid functional (including the second-order perturbational contribution) can be generated. This potential is shown to provide an improved description of orbital energies as vertical ionization potentials (IPs), relative to the perturbation-free KS potential typically used. Based on this observation, we suggest that DH-DFAs should be constructed in such a way that the potential provides accurate orbital energies. As a proof of principle, the B2-PLYP functional is reparameterized to obtain the IP-optimized B2IP-PLYP functional, using a small set of vertical IPs and atomization energies as reference data. This functional is shown to outperform B2-PLYP in a wide range of benchmarks and is en par with the related B2GP-PLYP. In particular, it is shown to be the most reliable choice in electronically difficult and multireference cases. [ABSTRACT FROM AUTHOR]

In this paper, a two-layer scheme is outlined for the coupled coherent states (CCS) method, dubbed two-layer CCS (2L-CCS). The theoretical framework is motivated by that of the multiconfigurational Ehrenfest method, where different dynamical descriptions are used for different subsystems of a quantum mechanical system. This leads to a flexible representation of the wavefunction, making the method particularly suited to the study of composite systems. It was tested on a 20-dimensional asymmetric system-bath tunnelling problem, with results compared to a benchmark calculation, as well as existing CCS, matching-pursuit/split-operator Fourier transform, and configuration interaction expansion methods. The two-layer method was found to lead to improved short and long term propagation over standard CCS, alongside improved numerical efficiency and parallel scalability. These promising results provide impetus for future development of the method for on-the-fly direct dynamics calculations. [ABSTRACT FROM AUTHOR]

One of the challenges facing on-the-fly ab initio semiclassical time evolution is the large expense needed to converge the computation. In this paper, we suggest that a significant saving in computational effort may be achieved by employing a semiclassical initial value representation (SCIVR) of the quantum propagator based on the Heisenberg interaction representation. We formulate and test numerically a modification and simplification of the previous semiclassical interaction representation of Shao and Makri [J. Chem. Phys. 113, 3681 (2000)]. The formulation is based on the wavefunction form of the semiclassical propagation instead of the operator form, and so is simpler and cheaper to implement. The semiclassical interaction representation has the advantage that the phase and prefactor vary relatively slowly as compared to the "standard" SCIVR methods. This improves its convergence properties significantly. Using a one-dimensional model system, the approximation is compared with Herman-Kluk's frozen Gaussian and Heller's thawed Gaussian approximations. The convergence properties of the interaction representation approach are shown to be favorable and indicate that the interaction representation is a viable way of incorporating on-the-fly force field information within a semiclassical framework. [ABSTRACT FROM AUTHOR]

The aim of this paper is to unravel the physical phenomena involved in the calculation of the spin density of the CuCl2 and [CuCl4]2- systems using wave function methods. Various types of wave functions are used here, both variational and perturbative, to analyse the effects impacting the spin density. It is found that the spin density on the chlorine ligands strongly depends on the mixing between two types of valence bond structures. It is demonstrated that the main difficulties found in most of the previous studies based on wave function methods come from the fact that each valence bond structure requires a different set of molecular orbitals and that using a unique set of molecular orbitals in a variational procedure leads to the removal of one of them from the wave function. Starting from these results, a method to compute the spin density at a reasonable computational cost is proposed. [ABSTRACT FROM AUTHOR]

We developed a simple method to calculate exchange repulsion between a quantum mechanical (QM) solute and a molecular mechanical (MM) molecule in the QM/MM approach. In our method, the size parameter in the Buckingham type potential for the QM solute is directly determined in terms of the one-electron wave functions of the solute. The point of the method lies in the introduction of the exchange core function (ECF) defined as a Slater function which mimics the behavior of the exterior electron density at the QM/MM boundary region. In the present paper, the ECF was constructed in terms of the Becke-Roussel (BR) exchange hole function. It was demonstrated that the ECF yielded by the BR procedure can faithfully reproduce the radial behavior of the electron density of a QM solute. The size parameter of the solute as well as the exchange repulsion are, then, obtained using the overlap model without any fitting procedure. To examine the efficiency of the method, it was applied to calculation of the exchange repulsions for minimal QM/MM systems, hydrogen-bonded water dimer, and H3O+-H2O. We found that our approach is able to reproduce the potential energy curves for these systems showing reasonable agreements with those given by accurate full quantum chemical calculations. [ABSTRACT FROM AUTHOR]

Multideterminant wavefunctions, while having a long history in quantum chemistry, are increasingly being used in highly accurate quantum Monte Carlo calculations. Since the accuracy of QMC is ultimately limited by the quality of the trial wavefunction, multi-Slater determinants wavefunctions offer an attractive alternative to Slater-Jastrow and more sophisticated wavefunction ansatz for several reasons. They can be efficiently calculated, straightforwardly optimized, and systematically improved by increasing the number of included determinants. In spite of their potential, however, the convergence properties of multi-Slater determinant wavefunctions with respect to orbital set choice and excited determinant selection are poorly understood, which hinders the application of these wavefunctions to large systems and solids. In this paper, by performing QMC calculations on the equilibrium and stretched carbon dimer, we find that convergence of the recovered correlation energy with respect to number of determinants can depend quite strongly on basis set and determinant selection methods, especially where there is strong correlation. We demonstrate that properly chosen orbital sets and determinant selection techniques from quantum chemistry methods can dramatically reduce the required number of determinants (and thus the computational cost) to reach a given accuracy, which we argue shows clear need for an automatic QMC-only method for selecting determinants and generating optimal orbital sets. [ABSTRACT FROM AUTHOR]

The free-complement (FC) method is a general method for solving the Schrödinger equation (SE): The produced wave function has the potentially exact structure as the solution of the Schrödinger equation. The variables included are determined either by using the variational principle (FC-VP) or by imposing the local Schrödinger equations (FC-LSE) at the chosen set of the sampling points. The latter method, referred to as the local Schrödinger equation (LSE) method, is integral-free and therefore applicable to any atom and molecule. The purpose of this paper is to formulate the basic theories of the LSE method and explain their basic features. First, we formulate three variants of the LSE method, the AB, HS, and HTQ methods, and explain their properties. Then, the natures of the LSE methods are clarified in some detail using the simple examples of the hydrogen atom and the Hooke's atom. Finally, the ideas obtained in this study are applied to solving the SE of the helium atom highly accurately with the FC-LSE method. The results are very encouraging: we could get the world's most accurate energy of the helium atom within the sampling-type methodologies, which is comparable to those obtained with the FC-VP method. Thus, the FC-LSE method is an easy and yet a powerful integral-free method for solving the Schrödinger equation of general atoms and molecules. [ABSTRACT FROM AUTHOR]

In quantum chemistry, Wick's theorem is an important tool to reduce products of fermionic creation and annihilation operators. It is especially useful in computations employing reference states. The original theorem has been generalized to tackle multiconfigurational wave functions or nonorthogonal orbitals. One particular issue of the resulting structure is that the metric and density matrices are intertwined despite their different origin. Here, an alternative, rather general tensorial formulation of Wick's theorem is proposed. The main difference is the separation of the metric—the coefficients at normal-ordered operators become products of an n-electron density matrix element and the Pfaffian of a matrix formed by orbital overlaps. Different properties of the formalism are discussed, including the use of density cumulants, the particle–hole symmetry, and applications to transition density matrices, i.e., the case of different bra and ket reference states. The metric-separated version of Wick's theorem provides a platform for the derivation of various quantum chemical methods, especially those complicated by non-trivial reference states and nonorthogonality issues. [ABSTRACT FROM AUTHOR]

The transition of energy from the 4f to the 5d state is a fundamental element driving various applications, such as phosphors and optoelectronic devices. The positioning of the 4f ground states and the 5d excited states significantly influences this energy shift. In our research, we delve into the placement of these states utilizing a hybrid density functional theory (DFT) combined with spin–orbit coupling (SOC) via the supercell method. Additionally, we scrutinize the transition energy, applying the constrained density functional theory (cDFT) approach in conjunction with the ΔSCF method. Our study illustrates that the synergy of cDFT and SOC generates a discrepancy of about 2% for Ce1 and 4% for Ce2 when comparing the calculated results to experimental data. Moreover, We have determined the positions of the 4f ground states to be 2.73 eV above the Valence Band Maximum (VBM) for Ce1 and 2.70 eV for Ce2. We also note a tight correlation between the 5d levels identified in the experimental data and the theoretical outcomes derived from wave function calculations at the CASPT2 accuracy level. [ABSTRACT FROM AUTHOR]

Quantum-chemical calculations aimed at deriving magnetic coupling constants in exchange-coupled spin clusters commonly utilize a broken-symmetry (BS) approach. This involves calculating several distinct collinear spin configurations, predominantly by density-functional theory. The energies of these configurations are interpreted in terms of the Heisenberg model, H ̃ = ∑ i < j J i j s ̃ i ⋅ s ̃ j , to determine coupling constants Jij for spin pairs. However, this energy-based procedure has inherent limitations, primarily in its inability to provide information on isotropic spin interactions beyond those included in the Heisenberg model. Biquadratic exchange or multi-center terms, for example, are usually inaccessible and hence assumed to be negligible. The present work introduces a novel approach employing BS mean-field solutions, specifically Hartree–Fock wave functions, for the construction of effective spin Hamiltonians. This expanded method facilitates the extraction of a broader range of coupling parameters by considering not only the energies, but also Hamiltonian and overlap elements between different BS states. We demonstrate how comprehensive s = 1 2 Hamiltonians, including multi-center terms, can be straightforwardly constructed from a complete set of BS solutions. The approach is exemplified for small clusters within the context of the half-filled single-band Hubbard model. This allows to contrast the current strategy against exact results, thereby offering an enriched understanding of the spin-Hamiltonian construction from BS solutions. [ABSTRACT FROM AUTHOR]

The Kohn-Sham theory addresses the challenge of representing the kinetic energy by re-quantizing density functional theory at a level of non-interacting electrons. It transforms the many-electron problem into a fictitious non-interacting electron problem, with the many-electron effects concealed within the exchange-correlation (XC) energy, which is expressed in terms of the electron density ρ(r). Unlike the wave function, ρ(r) can be viewed as a classical quantity, and expressing the XC energy in terms of it circumvents the need for correlated wave functions. In this work, we once again employ the re-quantization strategy and determine the XC energy using a local one-particle Schrödinger equation. The ground-state eigenfunction of the corresponding Hamiltonian is a reference point (r) dependent orbital φr,σ(u, σ′) which is subsequently used to generate the XC hole and the XC energy. The spin coordinate is denoted by σ and u is the electron-electron separation. The one-particle equation for φr,σ(u, σ′) includes a local potential vr,σ(u, σ′) that we approximate using two simple physical constraints. We assess the approximation by applying it to the helium iso-electronic series, the homogeneous electron gas, and the dissociation of the hydrogen molecule. [ABSTRACT FROM AUTHOR]

The core valence separation (CVS) approximation is the most employed strategy to prevent the variational collapse of standard wave function optimization when attempting to compute electronic states bearing one or more electronic vacancies in core orbitals. Here, we explore the spurious consequences of this approximation on the properties of the computed core hole states. We especially focus on the less studied case of double core hole (DCH) states, whose spectroscopic interest has recently been rapidly growing. We show that the CVS error leads to a systematic underestimation of DCH energies, a property in stark contrast with the case of single core hole states. We highlight that the CVS error can then be interpreted as an over relaxation effect and design a new correction strategy adapted to these specificities. [ABSTRACT FROM AUTHOR]

A recently initiated variant of density functional theory utilizes a set of spherically symmetric densities instead of the density. The exact functionals are unknown in the new theory akin to the standard density functional theory. In order to test approximate functionals exactly solvable models are introduced. A harmonic molecular ion, the analogue to the hydrogen molecule ion and a harmonic two-electron molecule showing analogy to the hydrogen molecule are proposed. It has been found that the wave function and the density can be given analytically. The exact spherical densities and the effective potentials of the Euler equations also have analytical form. It has been shown that the models can be easily extended to several "nuclei." [ABSTRACT FROM AUTHOR]

We suggest a consistent framework for the embedding of reduced-space correlated vibrational wave functions in a potential of the remaining modes and generalize this concept to arbitrary many subspaces. We present an implementation of this framework for vibrational coupled-cluster theory and response treatments. For C=O stretches of small molecules, we show that the embedded treatment accelerates convergence for enlarging subsets. For the water dimer and trimer as well as a water wire in bacteriorhodopsin, we investigate different partitioning schemes for the embedding approach: In the local partitioning of the vibrations, the modes dominated by motions in the same spatial region are correlated, whereas in the energy-based partitioning, modes of similar fundamental frequencies are correlated. In most cases, we obtain better agreement with superset reference results for the local partitioning than for energy-based partitioning. This work represents an important step toward multi-level methodologies in vibrational-structure theory required for its application to sizable (bio-)molecular systems. [ABSTRACT FROM AUTHOR]

Computing accurate yet efficient approximations to the solutions of the electronic Schrödinger equation has been a paramount challenge of computational chemistry for decades. Quantum Monte Carlo methods are a promising avenue of development as their core algorithm exhibits a number of favorable properties: it is highly parallel and scales favorably with the considered system size, with an accuracy that is limited only by the choice of the wave function Ansatz. The recently introduced machine-learned parametrizations of quantum Monte Carlo Ansätze rely on the efficiency of neural networks as universal function approximators to achieve state of the art accuracy on a variety of molecular systems. With interest in the field growing rapidly, there is a clear need for easy to use, modular, and extendable software libraries facilitating the development and adoption of this new class of methods. In this contribution, the DeepQMC program package is introduced, in an attempt to provide a common framework for future investigations by unifying many of the currently available deep-learning quantum Monte Carlo architectures. Furthermore, the manuscript provides a brief introduction to the methodology of variational quantum Monte Carlo in real space, highlights some technical challenges of optimizing neural network wave functions, and presents example black-box applications of the program package. We thereby intend to make this novel field accessible to a broader class of practitioners from both the quantum chemistry and the machine learning communities. [ABSTRACT FROM AUTHOR]

Electron occupations that arise from pure quantum states are restricted by a stringent set of conditions that are said to generalize the Pauli exclusion principle. These generalized Pauli constraints (GPCs) define the boundary of the set of one-electron reduced density matrices (1-RDMs) that are derivable from at least one N-electron wavefunction. In this paper, we investigate the sparsity of the Slater-determinant representation of the wavefunction that is a necessary, albeit not sufficient, condition for its 1-RDM to lie on the boundary of the set of pure N-representable 1-RDMs or in other words saturate one of the GPCs. The sparse wavefunction, we show, is exact not only for 3 electrons in 6 orbitals but also for 3 electrons in 8 orbitals. For larger numbers of electrons and/or orbitals in the lowest spin state, the exact wavefunction does not generally saturate one of the GPCs, and hence, the sparse representation is typically an approximation. Because the sparsity of the wavefunction is a necessary but not sufficient condition for saturation of one of the GPCs, optimization of the sparse wavefunction Ansatz to minimize the ground-state energy does not necessarily produce a wavefunction whose 1-RDM exactly saturates one of the GPCs. While the sparse Ansatz can be employed with arbitrary orbitals or optimized orbitals, in this paper, we explore the Ansatz with the natural orbitals from full configuration interaction, which yields an upper bound to the ground-state energy that equals the exact energy for a given basis set if the full-configuration-interaction wavefunction saturates the Ansatz’s GPC. With calculations on the boron isoelectronic sequence, the dinitrogen cation N 2 + , hydrogen chains, and cyclic conjugated π systems, we examine the quality of the sparse wavefunction Ansatz from the amount of correlation energy recovered. [ABSTRACT FROM AUTHOR]

In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representative of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and the relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which is industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow it to generate energies close to chemical accuracy compared to full configuration interaction or experimental results. In the field of material science, we study the electronic structure properties of cuprates through the quasi-1D Fermi–Hubbard model derived from CuBr2, where we find that trial wave functions with both significantly larger fidelities and lower energies over a mean-field solution do not necessarily lead to AFQMC results closer to the exact ground state energy. [ABSTRACT FROM AUTHOR]

In recent years, Green's function methods have garnered considerable interest due to their ability to target both charged and neutral excitations. Among them, the well-established GW approximation provides accurate ionization potentials and electron affinities and can be extended to neutral excitations using the Bethe–Salpeter equation (BSE) formalism. Here, we investigate the connections between various Green's function methods and evaluate their performance for charged and neutral excitations. Comparisons with other widely known second-order wave function methods are also reported. Additionally, we calculate the singlet-triplet gap of cycl[3,3,3]azine, a model molecular emitter for thermally activated delayed fluorescence, which has the particularity of having an inverted gap thanks to a substantial contribution from the double excitations. We demonstrate that, within the GW approximation, a second-order BSE kernel with dynamical correction is required to predict this distinctive characteristic. [ABSTRACT FROM AUTHOR]

As in the hierarchy of the Hartree–Fock theory in terms of spin symmetry, the extension of spin functions in a theory of two-electron units or geminal, was developed in this study. A trial wave function is constructed as an antisymmetrized product of geminals, in which singlet and triplet two-electron functions are fully mixed. We present a variational optimization method for this generalized pairing wave function in the strong orthogonality condition. The present method is considered an extension of the antisymmetrized product of strongly orthogonal geminals or perfect pairing generalized valence bond methods, maintaining the compactness of the trial wave function. The obtained broken-symmetry solutions were similar to the unrestricted Hartree–Fock wave functions in terms of spin contamination while giving lower energy due to the inclusion of the electron correlation effect in geminals. The degeneracy of the obtained broken-symmetry solutions in the Sz space is reported for the tested four-electron systems. [ABSTRACT FROM AUTHOR]