Exploring the role of mean-field potentials and short-range wave function behavior in the adiabatic connection.

In this article, we explore the construction of Hamiltonians with long-range interactions and their corrections using the short-range behavior of the wave function. A key aspect of our investigation is the examination of the one-particle potential, kept constant in our previous work, and the effects of its optimization on the adiabatic connection. Our methodology involves the use of a parameter-dependent potential dependent on a single parameter to facilitate practical computations. We analyze the energy errors and densities in a two-electron system (harmonium) under various conditions, employing different confinement potentials and interaction parameters. The study reveals that while the mean-field potential improves the expectation value of the physical Hamiltonian, it does not necessarily improve the energy of the system within the bounds of chemical accuracy. We also delve into the impact of density variations in adiabatic connections, challenging the common assumption that a mean field improves results. Our findings indicate that as long as energy errors remain within chemical accuracy, the mean field does not significantly outperform a bare potential. This observation is attributed to the effectiveness of corrections based on the short-range behavior of the wave function, a universal characteristic that diminishes the distinction between using a mean field or not.

Modified Expression for the Hamiltonian Expectation Value Exploiting the Short-Range Behavior of the Wave Function.

The expectation value of the Hamiltonian using a model wave function is widely used to estimate the eigenvalues of electronic Hamiltonians. We explore here a modified formula for models based on a long-range interaction. It scales differently the singlet and triplet components of the repulsion between electrons not present in the model (its short-range part). The scaling factors depend uniquely on the parameter used in defining the model interaction and are constructed using only exact properties. We show results for the ground states and low-lying excited states of Harmonium with two to six electrons. We obtain important improvements for the estimation of the exact energy, not only over the model energy but also over the expectation value of the Hamiltonian.

Transcorrelated selected configuration interaction in a bi-orthonormal basis and with a cheap three-body correlation factor.

In this work, we develop a mathematical framework for a selected configuration interaction (SCI) algorithm within a bi-orthogonal basis for transcorrelated (TC) calculations. The bi-orthogonal basis used here serves as the equivalent of the standard Hartree-Fock (HF) orbitals. However, within the context of TC, it leads to distinct orbitals for the left and right vectors. Our findings indicate that the use of such a bi-orthogonal basis allows for a proper definition of the frozen core approximation. In contrast, the use of HF orbitals results in bad error cancellations for ionization potentials and atomization energies (AE). Compared to HF orbitals, the optimized bi-orthogonal basis significantly reduces the positive part of the second-order energy (PT2), thereby facilitating the use of standard extrapolation techniques of hermitian SCI. While we did not observe a significant improvement in the convergence of the SCI algorithm, this is largely due to the use in this work of a simple three-body correlation factor introduced in a recent study. This correlation factor, which depends only on atomic parameters, eliminates the need for re-optimization of the correlation factor for molecular systems, making its use straightforward and user-friendly. Despite the simplicity of this correlation factor, we were able to achieve accurate results on the AE of a series of 14 molecules on a triple-zeta basis. We also successfully broke a double bond until the full dissociation limit while maintaining the size consistency property. This work thus demonstrates the potential of the BiO-TC-SCI approach in handling complex molecular systems.

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.

Reference Energies for Cyclobutadiene: Automerization and Excited States.

Cyclobutadiene is a well-known playground for theoretical chemists and is particularly suitable to test ground- and excited-state methods. Indeed, due to its high spatial symmetry, especially at the D4h square geometry but also in the D2h rectangular arrangement, the ground and excited states of cyclobutadiene exhibit multiconfigurational characters and single-reference methods, such as standard adiabatic time-dependent density-functional theory (TD-DFT) or standard equation-of-motion coupled cluster (EOM-CC), are notoriously known to struggle in such situations. In this work, using a large panel of methods and basis sets, we provide an extensive computational study of the automerization barrier (defined as the difference between the square and rectangular ground-state energies) and the vertical excitation energies at D2h and D4h equilibrium structures. In particular, selected configuration interaction (SCI), multireference perturbation theory (CASSCF, CASPT2, and NEVPT2), and coupled-cluster (CCSD, CC3, CCSDT, CC4, and CCSDTQ) calculations are performed. The spin-flip formalism, which is known to provide a qualitatively correct description of these diradical states, is also tested within TD-DFT (combined with numerous exchange-correlation functionals) and the algebraic diagrammatic construction [ADC(2)-s, ADC(2)-x, and ADC(3)]. A theoretical best estimate is defined for the automerization barrier and for each vertical transition energy.

The effect of uncertainty on building blocks in molecules.

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.

Extension of selected configuration interaction for transcorrelated methods.

In this work, we present an extension of popular selected configuration interaction (SCI) algorithms to the Transcorrelated (TC) framework. Although we used in this work the recently introduced one-parameter correlation factor [E. Giner, J. Chem. Phys. 154, 084119 (2021)], the theory presented here is valid for any correlation factor. Thanks to the formalization of the non-Hermitian TC eigenvalue problem as a search of stationary points for a specific functional depending on both left- and right-functions, we obtain a general framework, allowing for different choices for both the selection criterion in SCI and the second order perturbative correction to the energy. After numerical investigations on different second-row atomic and molecular systems in increasingly large basis sets, we found that taking into account the non-Hermitian character of the TC Hamiltonian in the selection criterion is mandatory to obtain a fast convergence of the TC energy. In addition, selection criteria based on either the first order coefficient or the second order energy lead to significantly different convergence rates, which is typically not the case in the usual Hermitian SCI. Regarding the convergence of the total second order perturbation energy, we find that the quality of the left-function used in the equations strongly affects the quality of the results. Within the near-optimal algorithm proposed here, we find that the SCI expansion in the TC framework converges faster than the usual SCI in terms of both the basis set and the number of Slater determinants.

Accurate full configuration interaction correlation energy estimates for five- and six-membered rings.

Following our recent work on the benzene molecule [P.-F. Loos, Y. Damour, and A. Scemama, J. Chem. Phys. 153, 176101 (2020)], motivated by the blind challenge of Eriksen et al. [J. Phys. Chem. Lett. 11, 8922 (2020)] on the same system, we report accurate full configuration interaction (FCI) frozen-core correlation energy estimates for 12 five- and six-membered ring molecules (cyclopentadiene, furan, imidazole, pyrrole, thiophene, benzene, pyrazine, pyridazine, pyridine, pyrimidine, s-tetrazine, and s-triazine) in the standard correlation-consistent double-ζ Dunning basis set (cc-pVDZ). Our FCI correlation energy estimates, with an estimated error smaller than 1 millihartree, are based on energetically optimized-orbital selected configuration interaction calculations performed with the configuration interaction using a perturbative selection made iteratively algorithm. Having at our disposal these accurate reference energies, the respective performance and convergence properties of several popular and widely used families of single-reference quantum chemistry methods are investigated. In particular, we study the convergence properties of (i) the Møller-Plesset perturbation series up to fifth-order (MP2, MP3, MP4, and MP5), (ii) the iterative approximate coupled-cluster series CC2, CC3, and CC4, and (iii) the coupled-cluster series CCSD, CCSDT, and CCSDTQ. The performance of the ground-state gold standard CCSD(T) as well as the completely renormalized CC model, CR-CC(2,3), is also investigated. We show that MP4 provides an interesting accuracy/cost ratio, while MP5 systematically worsens the correlation energy estimates. In addition, CC3 outperforms CCSD(T) and CR-CC(2,3), as well as its more expensive parent CCSDT. A similar trend is observed for the methods including quadruple excitations, where the CC4 model is shown to be slightly more accurate than CCSDTQ, both methods providing correlation energies within 2 millihartree of the FCI limit.

The performance of CIPSI on the ground state electronic energy of benzene.

Following the recent work of Eriksen et al. [J. Phys. Chem. Lett. 11, 8922 (2020)], we report the performance of the configuration interaction using a perturbative selection made iteratively method on the non-relativistic frozen-core correlation energy of the benzene molecule in the cc-pVDZ basis. Following our usual protocol, we obtain a correlation energy of -863.4 mEh, which agrees with the theoretical estimate of -863 mEh proposed by Eriksen et al. [J. Phys. Chem. Lett. 11, 8922 (2020)] using an extensive array of highly accurate new electronic structure methods.

Taming the fixed-node error in diffusion Monte Carlo via range separation.

By combining density-functional theory (DFT) and wave function theory via the range separation (RS) of the interelectronic Coulomb operator, we obtain accurate fixed-node diffusion Monte Carlo (FN-DMC) energies with compact multi-determinant trial wave functions. In particular, we combine here short-range exchange-correlation functionals with a flavor of selected configuration interaction known as configuration interaction using a perturbative selection made iteratively (CIPSI), a scheme that we label RS-DFT-CIPSI. One of the take-home messages of the present study is that RS-DFT-CIPSI trial wave functions yield lower fixed-node energies with more compact multi-determinant expansions than CIPSI, especially for small basis sets. Indeed, as the CIPSI component of RS-DFT-CIPSI is relieved from describing the short-range part of the correlation hole around the electron-electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional CIPSI calculation. Importantly, by performing various numerical experiments, we evidence that the RS-DFT scheme essentially plays the role of a simple Jastrow factor by mimicking short-range correlation effects, hence avoiding the burden of performing a stochastic optimization. Considering the 55 atomization energies of the Gaussian-1 benchmark set of molecules, we show that using a fixed value of µ = 0.5 bohr-1 provides effective error cancellations as well as compact trial wave functions, making the present method a good candidate for the accurate description of large chemical systems.

A basis-set error correction based on density-functional theory for strongly correlated molecular systems.

We extend to strongly correlated molecular systems the recently introduced basis-set incompleteness correction based on density-functional theory (DFT) [E. Giner et al., J. Chem. Phys. 149, 194301 (2018)]. This basis-set correction relies on a mapping between wave-function calculations in a finite basis set and range-separated DFT (RSDFT) through the definition of an effective non-divergent interaction corresponding to the electron-electron Coulomb interaction projected in the finite basis set. This enables the use of RSDFT-type complementary density functionals to recover the dominant part of the short-range correlation effects missing in this finite basis set. To study both weak and strong correlation regimes, we consider the potential energy curves of the H10, N2, O2, and F2 molecules up to the dissociation limit, and we explore various approximations of complementary functionals fulfilling two key properties: spin-multiplet degeneracy (i.e., independence of the energy with respect to the spin projection Sz) and size consistency. Specifically, we investigate the dependence of the functional on different types of on-top pair densities and spin polarizations. The key result of this study is that the explicit dependence on the on-top pair density allows one to completely remove the dependence on any form of spin polarization without any significant loss of accuracy. Quantitatively, we show that the basis-set correction reaches chemical accuracy on atomization energies with triple-ζ quality basis sets for most of the systems studied here. In addition, the present basis-set incompleteness correction provides smooth potential energy curves along the whole range of internuclear distances.

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids-A case study in diamond.

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater-Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

Acute ischemic stroke treated with mechanical thrombectomy and fungal endocarditis: A case report and systematic review of the literature.

Fungal endocarditis is a rare clinical form of infective endocarditis. The main etiology of FE is Candida albicans but also Candida parapsilosis and the overall mortality is high. We report a case of an acute ischemic stroke treated by mechanical thrombectomy, with the histopathological analysis of the retrieved clot followed by the confirmation of fungal endocarditis. An extensive review of the literature has been proposed and three key points concerning the fungal endocarditis predisposing factors, the relation between thrombolysis and hemorrhagic risk and, finally, the importance of clot analysis have been discussed.

Chemically accurate excitation energies with small basis sets.

By combining extrapolated selected configuration interaction (sCI) energies obtained with the Configuration Interaction using a Perturbative Selection made Iteratively algorithm with the recently proposed short-range density-functional correction for basis-set incompleteness [E. Giner et al., J. Chem. Phys. 149, 194301 (2018)], we show that one can get chemically accurate vertical and adiabatic excitation energies with, typically, augmented double-ζ basis sets. We illustrate the present approach on various types of excited states (valence, Rydberg, and double excitations) in several small organic molecules (methylene, water, ammonia, carbon dimer, and ethylene). The present study clearly evidences that special care has to be taken with very diffuse excited states where the present correction does not catch the radial incompleteness of the one-electron basis set.

Selected configuration interaction dressed by perturbation.

Selected configuration interaction (sCI) methods including second-order perturbative corrections provide near full CI (FCI) quality energies with only a small fraction of the determinants of the FCI space. Here, we introduce both a state-specific and a multi-state sCI method based on the configuration interaction using a perturbative selection made iteratively (CIPSI) algorithm. The present method revises the reference (internal) space under the effect of its interaction with the outer space via the construction of an effective Hamiltonian, following the shifted-Bk philosophy of Davidson and co-workers. In particular, the multi-state algorithm removes the storage bottleneck of the effective Hamiltonian via a low-rank factorization of the dressing matrix. Illustrative examples are reported for the state-specific and multi-state versions.

Excitation energies from diffusion Monte Carlo using selected configuration interaction nodes.

Quantum Monte Carlo (QMC) is a stochastic method that has been particularly successful for ground-state electronic structure calculations but mostly unexplored for the computation of excited-state energies. Here, we show that within a Jastrow-free QMC protocol relying on a deterministic and systematic construction of nodal surfaces using selected configuration interaction (sCI) expansions, one is able to obtain accurate excitation energies at the fixed-node diffusion Monte Carlo (FN-DMC) level. This evidences that the fixed-node errors in the ground and excited states obtained with sCI wave functions cancel out to a large extent. Our procedure is tested on two small organic molecules (water and formaldehyde) for which we report all-electron FN-DMC calculations. For both the singlet and triplet manifolds, accurate vertical excitation energies are obtained with relatively compact multideterminant expansions built with small (typically double-ζ) basis sets.

Hybrid stochastic-deterministic calculation of the second-order perturbative contribution of multireference perturbation theory.

A hybrid stochastic-deterministic approach for computing the second-order perturbative contribution E(2) within multireference perturbation theory (MRPT) is presented. The idea at the heart of our hybrid scheme-based on a reformulation of E(2) as a sum of elementary contributions associated with each determinant of the MR wave function-is to split E(2) into a stochastic and a deterministic part. During the simulation, the stochastic part is gradually reduced by dynamically increasing the deterministic part until one reaches the desired accuracy. In sharp contrast with a purely stochastic Monte Carlo scheme where the error decreases indefinitely as t-1/2 (where t is the computational time), the statistical error in our hybrid algorithm displays a polynomial decay ∼t-n with n = 3-4 in the examples considered here. If desired, the calculation can be carried on until the stochastic part entirely vanishes. In that case, the exact result is obtained with no error bar and no noticeable computational overhead compared to the fully deterministic calculation. The method is illustrated on the F2 and Cr2 molecules. Even for the largest case corresponding to the Cr2 molecule treated with the cc-pVQZ basis set, very accurate results are obtained for E(2) for an active space of (28e, 176o) and a MR wave function including up to 2×107 determinants.

A Jeziorski-Monkhorst fully uncontracted multi-reference perturbative treatment. I. Principles, second-order versions, and tests on ground state potential energy curves.

The present paper introduces a new multi-reference perturbation approach developed at second order, based on a Jeziorski-Mokhorst expansion using individual Slater determinants as perturbers. Thanks to this choice of perturbers, an effective Hamiltonian may be built, allowing for the dressing of the Hamiltonian matrix within the reference space, assumed here to be a CAS-CI. Such a formulation accounts then for the coupling between the static and dynamic correlation effects. With our new definition of zeroth-order energies, these two approaches are strictly size-extensive provided that local orbitals are used, as numerically illustrated here and formally demonstrated in the Appendix. Also, the present formalism allows for the factorization of all double excitation operators, just as in internally contracted approaches, strongly reducing the computational cost of these two approaches with respect to other determinant-based perturbation theories. The accuracy of these methods has been investigated on ground-state potential curves up to full dissociation limits for a set of six molecules involving single, double, and triple bond breaking together with an excited state calculation. The spectroscopic constants obtained with the present methods are found to be in very good agreement with the full configuration interaction results. As the present formalism does not use any parameter or numerically unstable operation, the curves obtained with the two methods are smooth all along the dissociation path.

Alternative definition of excitation amplitudes in multi-reference state-specific coupled cluster.

A central difficulty of state-specific Multi-Reference Coupled Cluster (MR-CC) in the multi-exponential Jeziorski-Monkhorst formalism concerns the definition of the amplitudes of the single and double excitation operators appearing in the exponential wave operators. If the reference space is a complete active space (CAS), the number of these amplitudes is larger than the number of singly and doubly excited determinants on which one may project the eigenequation, and one must impose additional conditions. The present work first defines a state-specific reference-independent operator T∼^m which acting on the CAS component of the wave function |Ψ0m⟩ maximizes the overlap between (1+T∼^m)|Ψ0m⟩ and the eigenvector of the CAS-SD (Singles and Doubles) Configuration Interaction (CI) matrix |ΨCAS-SDm⟩. This operator may be used to generate approximate coefficients of the triples and quadruples, and a dressing of the CAS-SD CI matrix, according to the intermediate Hamiltonian formalism. The process may be iterated to convergence. As a refinement towards a strict coupled cluster formalism, one may exploit reference-independent amplitudes provided by (1+T∼^m)|Ψ0m⟩ to define a reference-dependent operator T^m by fitting the eigenvector of the (dressed) CAS-SD CI matrix. The two variants, which are internally uncontracted, give rather similar results. The new MR-CC version has been tested on the ground state potential energy curves of 6 molecules (up to triple-bond breaking) and two excited states. The non-parallelism error with respect to the full-CI curves is of the order of 1 mEh.

Quantum Monte Carlo with very large multideterminant wavefunctions.

An algorithm to compute efficiently the first two derivatives of (very) large multideterminant wavefunctions for quantum Monte Carlo calculations is presented. The calculation of determinants and their derivatives is performed using the Sherman-Morrison formula for updating the inverse Slater matrix. An improved implementation based on the reduction of the number of column substitutions and on a very efficient implementation of the calculation of the scalar products involved is presented. It is emphasized that multideterminant expansions contain in general a large number of identical spin-specific determinants: for typical configuration interaction-type wavefunctions the number of unique spin-specific determinants Ndetσ ( σ=↑,↓) with a non-negligible weight in the expansion is of order O(Ndet). We show that a careful implementation of the calculation of the Ndet -dependent contributions can make this step negligible enough so that in practice the algorithm scales as the total number of unique spin-specific determinants,  Ndet↑+Ndet↓, over a wide range of total number of determinants (here, Ndet up to about one million), thus greatly reducing the total computational cost. Finally, a new truncation scheme for the multideterminant expansion is proposed so that larger expansions can be considered without increasing the computational time. The algorithm is illustrated with all-electron fixed-node diffusion Monte Carlo calculations of the total energy of the chlorine atom. Calculations using a trial wavefunction including about 750,000 determinants with a computational increase of ∼400 compared to a single-determinant calculation are shown to be feasible. © 2016 Wiley Periodicals, Inc.