Search

Quantum computing promises a computational advantage over classical methods in electronic-structure calculations, with expected applications in drug design and materials science. Accessing a quantitative description of chemical systems while minimizing quantum resources, such as the number of qubits, is an essential challenge given the limited capabilities of current quantum processors. We provide a shortcut towards quantum computations at chemical accuracy by approaching the complete-basis-set limit (CBS) through integrating density-functional theory into quantum algorithms via density-based basis-set corrections coupled to basis-sets crafted on-the-fly and specifically adapted to a given system/user-defined qubit budget. The approach self-consistently accelerates the basis-set convergence, improving electronic densities, ground-state energies, and first-order properties such as dipole moments. It can also serve as a classical, a posteriori, energy correction to quantum hardware calculations. The strategy is assessed using GPU-accelerated state-vector emulation up to 32 qubits. We converge the ground-state energies of four systems (He, Be, H$_2$, LiH) within chemical accuracy of the CBS full-configuration-interaction reference, while offering a systematic increase of accuracy beyond a double-zeta quality for various molecules up to the H$_8$ hydrogen chain. We also obtain dissociation curves for H$_2$ and LiH that reach the CBS limit whereas for the challenging simulation of the N$_2$ triple-bond breaking, we achieve a near-triple-zeta quality at the cost of a minimal basis-set. This hybrid strategy allows us to obtain quantitative results that would otherwise require brute-force quantum simulations using far more than 100 logical qubits, thereby opening up opportunities to explore real-world chemistry with reasonable computational resources.

We present the first application to real molecular systems of the recently proposed linear-response theory for the density-based basis-set correction method [J. Chem. Phys. 158, 234107 (2023)]. We apply this approach to accelerate the basis-set convergence of excitation energies in the equation-of-motion coupled-cluster singles doubles (EOM-CCSD) method. We use an approximate linear-response framework which neglects the second-order derivative of the basis-set correction density functional and consists in simply adding to the usual Hamiltonian the one-electron potential generated by the first-order derivative of the functional. This additional basis-set correction potential is evaluated at the Hartree-Fock density, leading to a very computationally cheap basis-set correction. We tested this approach over a set of about 30 excitation energies computed for five small molecular systems and found that the excitation energies from the ground state to Rydberg states are the main source of basis-set error. These excitation energies systematically increase when the size of the basis set is increased, suggesting a biased description in favour of the excited state. Despite the simplicity of the present approach, the results obtained with the basis-set corrected EOM-CCSD method are encouraging as they yield to a mean absolute deviation of 0.02 eV for the aug-cc-pVTZ basis set, while it is of 0.04 eV using the standard EOM-CCSD method. This might open the path to an alternative to explicitly correlated approaches to accelerate the basis-set convergence of excitation energies.

This work reports an efficient density-fitting implementation of the density-based basis-set correction (DBBSC) method in the MOLPRO software. This method consists in correcting the energy calculated by a wave-function method with a given basis set by an adapted basis-set correction density functional incorporating the short-range electron correlation effects missing in the basis set, resulting in an accelerated convergence to the complete-basis-set limit. Different basis-set correction density-functional approximations are explored and the complementary-auxiliary-basis-set single-excitation correction is added. The method is tested on a benchmark set of reaction energies at the second-order M{\o}ller-Plesset (MP2) level and a comparison with the explicitly correlated MP2-F12 method is provided. The results show that the DBBSC method greatly accelerates the basis convergence of MP2 reaction energies, without reaching the accuracy of the MP2-F12 method but with a lower computational cost., Comment: 6 pages, 2 figures, 1 table

The basis-set correction method based on density-functional theory consists in correcting the energy calculated by a wave-function method with a given basis set by a density functional. This basis-set correction density functional incorporates the short-range electron correlation effects missing in the basis set. This results in accelerated basis convergences of ground-state energies to the complete-basis-set limit. In this work, we extend the basis-set correction method to a linear-response formalism for calculating excited-state energies. We give the general linear-response equations, as well as the more specific equations for configuration-interaction wave functions. As a proof of concept, we apply this approach to the calculations of excited-state energies in a one-dimensional two-electron model system with harmonic potential and a Dirac-delta electron-electron interaction. The results obtained with full-configuration-interaction wave functions expanded in a basis of Hermite functions and a local-density-approximation basis-set correction functional show that the present approach does not help in accelerating the basis convergence of excitation energies. However, we show that it significantly accelerates basis convergences of excited-state total energies.

We consider a one-dimensional effective quantum electrodynamics (QED) model of the relativistic hydrogen-like atom using delta-potential interactions. We discuss the general exact theory and the Hartree-Fock approximation. The present one-dimensional effective QED model shares the essential physical feature of the three-dimensional theory: the nuclear charge polarizes the vacuum state (creation of electron-positron pairs) which results in a QED Lamb-type shift of the bound-state energy. Yet, this 1D effective QED model eliminates some of the most serious technical difficulties of the three-dimensional theory coming from renormalization. We show how to calculate the vacuum-polarization density at zeroth order in the two-particle interaction and the QED Lamb-type shift of the bound-state energy at first order in the two-particle interaction. The present work may be considered as a step toward the development of a quantum-chemistry effective QED theory of atoms and molecules.

We develop semiclassical approximations for calculating photoabsorption cross sections beyond the continuum threshold in quantum many-body systems. These approximations use the fully quantum-mechanical Wigner function of the ground state and semiclassical expansions only for the part of the cross section depending on the continuum states, thus avoiding the difficult explicit calculation of the continuum states. Even though the approach is general, we test it in electronic-structure theory for the photoionization cross sections of the hydrogen and helium atoms. The results suggest that these semiclassical approximations can be used to obtain good estimates of cross sections at high energy.

We consider the calculations of photoionization spectra and core resonances of open-shell systems using range-separated time-dependent density-functional theory. Specifically, we use the time-dependent range-separated hybrid (TDRSH) scheme, combining a long-range Hartree-Fock (HF) exchange potential and kernel with a short-range potential and kernel from a local density-functional approximation, and the time-dependent locally range-separated hybrid (TDLRSH) scheme, which uses a local range-separation parameter. To efficiently perform the calculations, we formulate a spin-unrestricted linear-response Sternheimer approach in a non-orthogonal B-spline basis set and using appropriate frequency-dependent boundary conditions. We illustrate this approach on the Li atom, which suggests that TDRSH and TDLRSH are adequate simple methods for estimating single-electron photoionization spectra of open-shell systems., Comment: arXiv admin note: substantial text overlap with arXiv:2203.05951

The VB-QMC method is presented in this chapter. It consists of using in quantum Monte Carlo (QMC) approaches with a wave function expressed as a usually short expansion of classical Valence-Bond (VB) structures supplemented by a Jastrow factor to account for dynamical correlation. Two variants exist: the VB-VMC (using variational Monte Carlo) and VB-DMC (using diffusion Monte Carlo) methods. QMC algorithms circumvent the notorious non-orthogonality issue of classical VB approaches, and allow highly efficient calculations on massively parallel machines. Calculation of VB weights and resonance energies are possible at the VB-VMC level, which makes VB-VMC a correlated method retaining all the interpretative capabilities of classical VB methods. Several recent applications are shown to illustrate the potential of this method as a modern alternative to classical VB methods to study ground and excited states of molecules., Comment: book chapter submitted to Comprehensive Computational Chemistry

We propose to compute physical properties by Monte Carlo calculations using conditional expectation values. The latter are obtained on top of the usual Monte Carlo sampling by partitioning the physical space in several subspaces or fragments, and subsampling each fragment (i.e., performing side walks) while freezing the environment. No bias is introduced and a zero-variance principle holds in the limit of separability, i.e. when the fragments are independent. In practice, the usual bottleneck of Monte Carlo calculations -- the scaling of the statistical fluctuations as a function of the number of particles N -- is relieved for extensive observables. We illustrate the method in variational Monte Carlo on the 2D Hubbard model and on metallic hydrogen chains using Jastrow-Slater wave functions. A factor O(N) is gained in numerical efficiency.

We explore the merits of linear-response range-separated time-dependent density-functional theory (TDDFT) for the calculation of photoionization spectra. We consider two variants of range-separated TDDFT, namely the timedependent range-separated hybrid (TDRSH) scheme which uses a global range-separation parameter and the timedependent locally range-separated hybrid (TDLRSH) which uses a local range-separation parameter, and compare with standard time-dependent local-density approximation (TDLDA) and time-dependent Hartree-Fock (TDHF). We show how to calculate photoionization spectra with these methods using the Sternheimer approach formulated in a nonorthogonal B-spline basis set with appropriate frequency-dependent boundary conditions. We illustrate these methods on the photoionization spectrum of the Be atom, focusing in particular on the core resonances. Both the TDRSH and TDLRSH photoionization spectra are found to constitute a large improvement over the TDLDA photoionization spectrum and a more modest improvement over the TDHF photoionization spectrum.

The present work proposes an approach to obtain a basis-set correction based on density-functional theory (DFT) for the computation of molecular properties in wave-function theory (WFT). This approach allows one to accelerate the basis-set convergence of any energy derivative of a non-variational WFT method, generalizing previous works on the DFT-based basis-set correction where either only ground-state energies could be computed with non-variational wave functions [J. Phys. Chem. Lett. 10, 2931 (2019)] or properties where computed as expectation values over variational wave functions [J. Chem. Phys. 155, 044109 (2021)]. This work focuses on the basis-set correction of dipole moments in coupled-cluster with single, double, and perturbative triple excitations (CCSD(T)), which is numerically tested on a set of fourteen molecules with dipole moments covering two orders of magnitude. As the basis-set correction relies only on Hartree-Fock densities, its computational cost is marginal with respect to the one of the CCSD(T) calculations. Statistical analysis of the numerical results shows a clear improvement of the basis convergence of the dipole moment with respect to the usual CCSD(T) calculations.

We reexamine the recently introduced basis-set correction theory based on density-functional theory consisting in correcting the basis-set incompleteness error of wave-function methods using a density functional. We use a one-dimensional model Hamiltonian with delta-potential interactions which has the advantage of making easier to perform a more systematic analysis than for three-dimensional Coulombic systems while keeping the essence of the slow basis convergence problem of wave-function methods. We provide some mathematical details about the theory and propose a new variant of basis-set correction which has the advantage of being suited to the development of an adapted local-density approximation. We show indeed how to develop a local-density approximation for the basis-set correction functional which is automatically adapted to the basis set employed, without resorting to range-separated density-functional theory as in previous works, but using instead a finite uniform electron gas whose electron-electron interaction is projected on the basis set. The work puts the basis-set correction theory on firmer grounds and provides an interesting strategy for the improvement of this approach.

The semistochastic heat-bath configuration interaction (SHCI) method is a selected configuration interaction plus perturbation theory method that has provided near-full configuration interaction (FCI) levels of accuracy for many systems with both single- and multi-reference character. However, obtaining accurate energies in the complete basis set limit is hindered by the slow convergence of the FCI energy with respect to basis size. Here we show that the recently developed basis-set correction method based on range-separated density-functional theory can be used to significantly speed up basis-set convergence in SHCI calculations. In particular, we study two such schemes that differ in the functional used, and apply them to transition metal atoms and monoxides to obtain total, ionization, and dissociation energies well converged to the complete-basis-set limit within chemical accuracy., Comment: 5 figures, supplemental material

This work provides a self-consistent extension of the recently proposed density-based basis-set correction method for wave-function electronic-structure calculations [J. Chem. Phys. 149, 194301 (2018)]. In contrast to the previously used approximation where the basis-set correction density functional was a posteriori added to the energy from a wave-function calculation, here the energy minimization is performed including the basis-set correction. Compared to the non-self-consistent approximation, this allows one to lower the total energy and change the wave function under the effect of the basis-set correction. This work addresses two main questions: i) What is the change in total energy compared to the non-self-consistent approximation, and ii) can we obtain better properties, namely dipole moments, with the basis-set corrected wave functions? We implement the present formalism with two different basis-set correction functionals and test it on different molecular systems. The main results of the study are that i) the total energy lowering obtained by the self-consistent approach is extremely small, which justifies the use of the non-self-consistent approximation, and ii) the dipole moments obtained from the basis-set corrected wave functions are improved, being already close to their complete-basis-set values with triple-zeta basis sets. Thus, the present study further confirms the soundness of the density-based basis-set correction scheme.

In this chapter, we provide a review of ground-state Kohn-Sham density-functional theory of electronic systems and some of its extensions, we present exact expressions and constraints for the exchange and correlation density functionals, and we discuss the main families of approximations for the exchange-correlation energy: semilocal approximations, single-determinant hybrid approximations, multideterminant hybrid approximations, dispersion-corrected approximations, as well as orbital-dependent exchange-correlation density functionals. The chapter aims at providing both a consistent bird's-eye view of the field and a detailed description of some of the most used approximations. It is intended to be readable by chemists/physicists and applied mathematicians., Comment: Chapter in the book "Density Functional Theory" edited by Eric Canc\`es and Gero Friesecke

A relativistic density-functional theory based on a Fock-space effective quantum-electrodynamics (QED) Hamiltonian using the Coulomb or Coulomb-Breit two-particle interaction is developed. This effective QED theory properly includes the effects of vacuum polarization through the creation of electron-positron pairs but does not include explicitly the photon degrees of freedom. It is thus a more tractable alternative to full QED for atomic and molecular calculations. Using the constrained-search formalism, a Kohn-Sham scheme is formulated in a quite similar way to non-relativistic density-functional theory, and some exact properties of the involved density functionals are studied, namely charge-conjugation symmetry and uniform coordinate scaling. The usual no-pair Kohn-Sham scheme is obtained as a well-defined approximation to this relativistic density-functional theory., Comment: version 2, with a few improvements, to appear in SciPost Chemistry

We construct the complementary short-range correlation relativistic local-density-approximation functional to be used in relativistic range-separated density-functional theory based on a Dirac-Coulomb Hamiltonian in the no-pair approximation. For this, we perform relativistic random-phase-approximation calculations of the correlation energy of the relativistic homogeneous electron gas with a modified electron-electron interaction, we study the high-density behavior, and fit the results to a parametrized expression. The obtained functional should eventually be useful for electronic-structure calculations of strongly correlated systems containing heavy elements.

The recently developed semistochastic heat-bath configuration interaction (SHCI) method is a systematically improvable selected configuration interaction plus perturbation theory method capable of giving essentially exact energies for larger systems than is possible with other such methods. We compute SHCI atomization energies for 55 molecules which have been used as a test set in prior studies because their atomization energies are known from experiment. Basis sets from cc-pVDZ to cc-pV5Z are used, totaling up to 500 orbitals and a Hilbert space of $10^{32}$ Slater determinants for the largest molecules. For each basis, an extrapolated energy well within chemical accuracy (1 kcal/mol or 1.6 mHa/mol) of the exact energy for that basis is computed using only a tiny fraction of the entire Hilbert space. We also use our almost exact energies to benchmark coupled-cluster [CCSD(T)] energies. The energies are extrapolated to the complete basis set limit and compared to the experimental atomization energies. The extrapolations are done both without and with a basis-set correction based on density-functional theory. The mean absolute deviations from experiment for these extrapolations are 0.46 kcal/mol and 0.51 kcal/mol, respectively. Orbital optimization methods used to obtain improved convergence of the SHCI energies are also discussed., Comment: 6 figures, 4 supplemental files

We develop relativistic short-range exchange energy functionals for four-component relativistic range-separated density-functional theory using a Dirac-Coulomb Hamiltonian in the no-pair approximation. We show how to improve the short-range local-density approximation exchange functional for large range-separation parameters by using the on-top exchange pair density as a new variable. We also develop a relativistic short-range generalized-gradient approximation exchange functional which further increases the accuracy for small range-separation parameters. Tests on the helium, beryllium, neon, and argon isoelectronic series up to high nuclear charges show that this latter functional gives exchange energies with a maximal relative percentage error of 3 %. The development of this exchange functional represents a step forward for the application of four-component relativistic range-separated density-functional theory to chemical compounds with heavy elements., Comment: Journal of Chemical Physics, American Institute of Physics, In press

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-zeta quality basis sets for most of the systems studied here. Also, the present basis-set incompleteness correction provides smooth potential energy curves along the whole range of internuclear distances.

We extend the range-separated double-hybrid RSH+MP2 method [J. G. Angyan et al., Phys. Rev. A 72, 012510 (2005)], combining long-range HF exchange and MP2 correlation with a short-range density functional, to a fully self-consistent version using the optimized-effective-potential technique in which the orbitals are obtained from a local potential including the long-range HF and MP2 contributions. We test this approach, that we name RS-OEP2, on a set of small closed-shell atoms and molecules. For the commonly used value of the range-separation parameter $\mu=0.5$ bohr$^{-1}$, we find that self-consistency does not seem to bring any improvement for total energies, ionization potentials, and electronic affinities. However, contrary to the non-self-consistent RSH+MP2 method, the present RS-OEP2 method gives a LUMO energy which physically corresponds to a neutral excitation energy and gives local exchange-correlation potentials which are reasonably good approximations to the corresponding Kohn-Sham quantities. At a finer scale, we find that RS-OEP2 gives largely inaccurate correlation potentials and correlated densities, which points to the need of further improvement of this type of range-separated double hybrids., Comment: Journal of Chemical Theory and Computation, American Chemical Society, In press

Similar to other electron correlation methods, many-body perturbation theory methods based on Green functions, such as the so-called $GW$ approximation, suffer from the usual slow convergence of energetic properties with respect to the size of the one-electron basis set. This displeasing feature is due to lack of explicit electron-electron terms modeling the infamous Kato electron-electron cusp and the correlation Coulomb hole around it. Here, we propose a computationally efficient density-based basis set correction based on short-range correlation density functionals which significantly speeds up the convergence of energetics towards the complete basis set limit. The performance of this density-based correction is illustrated by computing the ionization potentials of the twenty smallest atoms and molecules of the GW100 test set at the perturbative $GW$ (or $G_0W_0$) level using increasingly large basis sets. We also compute the ionization potentials of the five canonical nucleobases (adenine, cytosine, thymine, guanine, and uracil) and show that, here again, a significant improvement is obtained., Comment: 11 pages, 2 figures (supporting information available)

By combining extrapolated selected configuration interaction (sCI) energies obtained with the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm with the recently proposed short-range density-functional correction for basis-set incompleteness [Giner et al.,J. Chem. Phys. 2018, 149, 194301], we show that one can get chemically accurate vertical and adiabatic excitation energies with, typically, augmented double-$\zeta$ 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., Comment: 9 pages, 6 figures

We construct range-separated double-hybrid schemes which combine coupled-cluster or random-phase approximations with a density functional based on a two-parameter Coulomb-attenuating-method-like decomposition of the electron-electron interaction. We find that the addition of a fraction of short-range electron-electron interaction in the wave-function part of the calculation is globally beneficial for the range-separated double-hybrid scheme involving a variant of the random-phase approximation with exchange terms. Even though the latter scheme is globally as accurate as the corresponding scheme employing only second-order M{{\o}}ller-Plesset perturbation theory for atomization energies, reaction barrier heights, and weak intermolecular interactions of small molecules, it is more accurate for the more complicated case of the benzene dimer in the stacked configuration. The present range-separated double-hybrid scheme employing a random-phase approximation thus represents a new member in the family of double hybrids with minimal empiricism which could be useful for general chemical applications., Comment: arXiv admin note: text overlap with arXiv:1804.03373

We report a universal density-based basis-set incompleteness correction that can be applied to any wave function method. The present correction, which appropriately vanishes in the complete basis set (CBS) limit, relies on short-range correlation density functionals (with multi-determinant reference) from range-separated density-functional theory (RS-DFT) to estimate the basis-set incompleteness error. Contrary to conventional RS-DFT schemes which require an \textit{ad hoc} range-separation \textit{parameter} $\mu$, the key ingredient here is a range-separation \textit{function} $\mu(\bf{r})$ that automatically adapts to the spatial non-homogeneity of the basis-set incompleteness error. As illustrative examples, we show how this density-based correction allows us to obtain CCSD(T) atomization and correlation energies near the CBS limit for the G2 set of molecules with compact Gaussian basis sets., Comment: 7 pages, 2 figures and 1 table (supporting information available)

We investigate the performance of the range-separated hybrid (RSH) scheme, which combines long-range Hartree-Fock (HF) and a short-range density-functional approximation (DFA), for calculating photoexcitation/photoionization spectra of the H and He atoms, using a B-spline basis set in order to correctly describe the continuum part of the spectra. The study of these simple systems allows us to quantify the influence on the spectra of the errors coming from the short-range exchange-correlation DFA and from the missing long-range correlation in the RSH scheme. We study the differences between using the long-range HF exchange (nonlocal) potential and the long-range exact exchange (local) potential. Contrary to the former, the latter supports a series of Rydberg states and gives reasonable photoexcitation/photoionization spectra, even without applying linear-response theory. The most accurate spectra are obtained with the linear-response time-dependent range-separated hybrid (TDRSH) scheme. In particular, for the He atom at the optimal value of the range-separation parameter, TDRSH gives slightly more accurate photoexcitation and photoioniza-tion spectra than standard linear-response time-dependent HF. More generally, the present work shows the potential of range-separated density-functional theory for calculating linear and nonlinear optical properties involving continuum states.

\textsc{Quantum Package} is an open-source programming environment for quantum chemistry specially designed for wave function methods. Its main goal is the development of determinant-driven selected configuration interaction (sCI) methods and multi-reference second-order perturbation theory (PT2). The determinant-driven framework allows the programmer to include any arbitrary set of determinants in the reference space, hence providing greater methodological freedoms. The sCI method implemented in \textsc{Quantum Package} is based on the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm which complements the variational sCI energy with a PT2 correction. Additional external plugins have been recently added to perform calculations with multireference coupled cluster theory and range-separated density-functional theory. All the programs are developed with the IRPF90 code generator, which simplifies collaborative work and the development of new features. \textsc{Quantum Package} strives to allow easy implementation and experimentation of new methods, while making parallel computation as simple and efficient as possible on modern supercomputer architectures. Currently, the code enables, routinely, to realize runs on roughly 2\,000 CPU cores, with tens of millions of determinants in the reference space. Moreover, we have been able to push up to 12\,288 cores in order to test its parallel efficiency. In the present manuscript, we also introduce some key new developments: i) a renormalized second-order perturbative correction for efficient extrapolation to the full CI limit, and ii) a stochastic version of the CIPSI selection performed simultaneously to the PT2 calculation at no extra cost., Comment: 21 pages, 8 figures

We introduce an approximation to the short-range correlation energy functional with multide-terminantal reference involved in a variant of range-separated density-functional theory. This approximation is a local functional of the density, the density gradient, and the on-top pair density, which locally interpolates between the standard Perdew-Burke-Ernzerhof correlation functional at vanishing range-separation parameter and the known exact asymptotic expansion at large range-separation parameter. When combined with (selected) configuration-interaction calculations for the long-range wave function, this approximation gives accurate dissociation energy curves of the H2, Li2, and Be2 molecules, and thus appears as a promising way to accurately account for static correlation in range-separated density-functional theory., Comment: The Journal of Chemical Physics, In press

A formally exact Bethe-Salpeter-like equation for the linear-response function is introduced with a kernel which depends only on the one frequency of the applied field. This is in contrast with the standard Bethe-Salpeter equation (BSE) which involves multiple-frequency integrals over the kernel and response functions. From the one-frequency kernel, known approximations are straightforwardly recovered. However, the present formalism lends itself to more powerful approximations. This is demonstrated with the exact analytical solution of the Hubbard molecule. Similarities and differences of the $GW$+BSE approach with the self-consistent random-phase approximation (RPA) is also discussed.

We lay out the extension of range-separated density-functional theory to a four-component relativistic frame-work using a Dirac-Coulomb-Breit Hamiltonian in the no-pair approximation. This formalism combines a wave-function method for the long-range part of the electron-electron interaction with a density(-current) functionalfor the short-range part of the interaction. We construct for this formalism a short-range exchange local-densityapproximation based on calculations on a relativistic homogeneous electron gas with a modified Coulomb-Breitelectron-electron interaction. More specifically, we provide the relativistic short-range Coulomb and Breit ex-change energies per particle of the relativistic homogeneous electron gas in the form of Pad e approximantswhich are systematically improvable to arbitrary accuracy. These quantities, as well as the associated effectiveCoulomb-Breit exchange hole, show the important impact of relativity on short-range exchange effects for highdensities.

The beryllium dimer and trimer are, despite their small number of electrons, excellent systems for assessing electronic-structure computational methods. With reference data provided by multi-reference averaged coupled-pair functional calculations, we assess several variants of range-separated density-functional theory, combining long-range second-order perturbation theory or coupled-cluster theory with a short-range density functional. The results show that i) long-range second-order perturbation theory is not sufficient, ii) long-range coupled-cluster theory gives reasonably accurate potential energy curves, but iii) provided a relatively large value of $\mu=1$ bohr$^{-1}$ for the range-separationparameter is used. The article is dedicated to the memory of J\'anos G. \'Angy\'an., Comment: Theoretical Chemistry Accounts: Theory, Computation, and Modeling, Springer Verlag, In press

A clear understanding of the mechanisms that control the electron dynamics in strong laser field is still a challenge that requires to be interpreted by advanced theory. Development of accurate theoretical and computational methods, able to provide a precise treatment of the fundamental processes generated in the strong field regime, is therefore crucial. A central aspect is the choice of the basis for the wave-function expansion. Accuracy in describing multiphoton processes is strictly related to the intrinsic properties of the basis, such as numerical convergence, computational cost, and representation of the continuum. By explicitly solving the 1D and 3D time-dependent Schr\"odinger equation for H$^{+}_{2}$ in presence of an intense electric field, we explore the numerical performance of using a real-space grid, a B-spline basis, and a Gaussian basis (improved by optimal Gaussian functions for the continuum). We analyze the performance of the three bases for high-harmonic generation and above-threshold ionization for H$^{+}_{2}$. In particular, for high-harmonic generation, the capability of the basis to reproduce the two-center interference and the hyper-Raman phenomena is investigated., Comment: 35 pages, 12 figures

The present work proposes to use density-functional theory (DFT) to correct for the basis-set error of wave-function theory (WFT). One of the key ideas developed here is to define a range-separation parameter which automatically adapts to a given basis set. The derivation of the exact equations are based on the Levy-Lieb formulation of DFT, which helps us to define a complementary functional which corrects uniquely for the basis-set error of WFT. The coupling of DFT and WFT is done through the definition of a real-space representation of the electron-electron Coulomb operator projected in a one-particle basis set. Such an effective interaction has the particularity to coincide with the exact electron-electron interaction in the limit of a complete basis set, and to be finite at the electron-electron coalescence point when the basis set is incomplete. The non-diverging character of the effective interaction allows one to define a mapping with the long-range interaction used in the context of range-separated DFT and to design practical approximations for the unknown complementary functional. Here, a local-density approximation is proposed for both full-configuration-interaction (FCI) and selected configuration-interaction approaches. Our theory is numerically tested to compute total energies and ionization potentials for a series of atomic systems. The results clearly show that the DFT correction drastically improves the basis-set convergence of both the total energies and the energy differences. For instance, a sub kcal/mol accuracy is obtained from the aug-cc-pVTZ basis set with the method proposed here when an aug-cc-pV5Z basis set barely reaches such a level of accuracy at the near FCI level.

A range-separated double-hybrid (RSDH) scheme which generalizes the usual range-separated hybrids and double hybrids is developed. This scheme consistently uses a two-parameter Coulomb-attenuating-method (CAM)-like decomposition of the electron-electron interaction for both exchange and correlation in order to combine Hartree-Fock exchange and second-order M{{\o}}ller-Plesset (MP2) correlation with a density functional. The RSDH scheme relies on an exact theory which is presented in some detail. Several semi-local approximations are developed for the short-range exchange-correlation density functional involved in this scheme. After finding optimal values for the two parameters of the CAM-like decomposition, the RSDH scheme is shown to have a relatively small basis dependence and to provide atomization energies, reaction barrier heights, and weak intermolecular interactions globally more accurate or comparable to range-separated MP2 or standard MP2. The RSDH scheme represents a new family of double hybrids with minimal empiricism which could be useful for general chemical applications., Comment: Journal of Chemical Physics, American Institute of Physics, In press

A key goal in quantum chemistry methods, whether ab initio or otherwise, is to achieve size consistency. In this manuscript we formulate the related idea of "Casimir-Polder size consistency" that manifests in long-range dispersion energetics. We show that local approximations in time-dependent density functional theory dispersion energy calculations violate the consistency condition because of incorrect treatment of highly non-local "xc kernel" physics, by up to 10% in our tests on closed-shell atoms., Comment: Journal of Chemical Theory and Computation, American Chemical Society, 2017

A G{\"o}rling-Levy (GL)-based perturbation theory along the range-separated adiabatic connection is assessed for the calculation of electronic excitation energies. In comparison with the Rayleigh-Schr{\"o}dinger (RS)-based perturbation theory introduced in a previous work [E. Rebolini, J. Toulouse, A. M. Teale, T. Helgaker, A. Savin, Mol. Phys. 113, 1740 (2015)], this GL-based perturbation theory keeps the ground-state density constant at each order and thus gives the correct ionization energy at each order. Excitation energies up to first order in the perturbation have been calculated numerically for the helium and beryllium atoms and the hydrogen molecule without introducing any density-functional approximations. In comparison with the RS-based perturbation theory, the present GL-based perturbation theory gives much more accurate excitation energies for Rydberg states but similar excitation energies for valence states.

Many chemical systems cannot be described by quantum chemistry methods based on a singlereference wave function. Accurate predictions of energetic and spectroscopic properties require a delicate balance between describing the most important configurations (static correlation) and obtaining dynamical correlation efficiently. The former is most naturally done through a multiconfigurational (MC) wave function, whereas the latter can be done by, e.g., perturbation theory. We have employed a different strategy, namely, a hybrid between multiconfigurational wave functions and density-functional theory (DFT) based on range separation. The method is denoted by MC short-range (sr) DFT and is more efficient than perturbative approaches as it capitalizes on the efficient treatment of the (short-range) dynamical correlation by DFT approximations. In turn, the method also improves DFT with standard approximations through the ability of multiconfigurational wave functions to recover large parts of the static correlation. Until now, our implementation was restricted to closed-shell systems, and to lift this restriction, we present here the generalization of MC-srDFT to open-shell cases. The additional terms required to treat open-shell systems are derived and implemented in the DALTON program. This new method for open-shell systems is illustrated on dioxygen and [Fe(H2O)6]3+., Comment: 37 pages, 3 figures, 4 tables, 1 appendix and 79 references Changes in v2: 1) Appendix B and reference 81 removed 2) Removed dublicated reference and corrected reference 31. 3) Added spin-charge cross terms to GGA (Appendix A). Code changed accordingly and GGA results recalculated. All GGA results are revised -only small modifications observed. Conclusions are unchanged

We propose a method for obtaining effective lifetimes of scattering electronic states for avoiding the artificially confinement of the wave function due to the use of incomplete basis sets in time-dependent electronic-structure calculations of atoms and molecules. In this method, using a fitting procedure, the lifetimes are extracted from the spatial asymptotic decay of the approximate scattering wave functions obtained with a given basis set. The method is based on a rigorous analysis of the complex-energy solutions of the Schr{\"o}dinger equation. It gives lifetimes adapted to any given basis set without using any empirical parameters. The method can be considered as an ab initio version of the heuristic lifetime model of Klinkusch et al. [J. Chem. Phys. 131, 114304 (2009)]. The method is validated on the H and He atoms using Gaussian-type basis sets for calculation of high-harmonic-generation spectra.

We present the extension of variational Monte Carlo (VMC) to the calculation of electronic excitation energies and oscillator strengths using time-dependent linear-response theory. By exploiting the analogy existing between the linear method for wave-function optimisation and the generalised eigenvalue equation of linear-response theory, we formulate the equations of linear-response VMC (LR-VMC). This LR-VMC approach involves the first-and second-order derivatives of the wave function with respect to the parameters. We perform first tests of the LR-VMC method within the Tamm-Dancoff approximation using single-determinant Jastrow-Slater wave functions with different Slater basis sets on some singlet and triplet excitations of the beryllium atom. Comparison with reference experimental data and with configuration-interaction-singles (CIS) results shows that LR-VMC generally outperforms CIS for excitation energies and is thus a promising approach for calculating electronic excited-state properties of atoms and molecules., Comment: Advances in Quantum Chemistry, 2017, Novel Electronic Structure Theory: General Innovations and Strongly Correlated Systems

We introduce an orbital-optimized double-hybrid (DH) scheme using the optimized-effective-potential (OEP) method. The orbitals are optimized using a local potential corresponding to the complete exchange-correlation energy expression including the second-order M{{\o}}ller-Plesset (MP2) correlation contribution. We have implemented a one-parameter version of this OEP-based self-consistent DH scheme using the BLYP density-functional approximation and compared it to the corresponding non-self-consistent DH scheme for calculations on a few closed-shell atoms and molecules. While the OEP-based self-consistency does not provide any improvement for the calculations of ground-state total energies and ionization potentials, it does improve the accuracy of electron affinities and restores the meaning of the LUMO orbital energy as being connected to a neutral excitation energy. Moreover, the OEP-based self-consistent DH scheme provides reasonably accurate exchange-correlation potentials and correlated densities., Comment: in The Journal of Chemical Physics, American Institute of Physics, 2016

We investigate fractional-charge and fractional-spin errors in range-separated density-functional theory. Specifically, we consider the range-separated hybrid (RSH) method which combines long-range Hartree-Fock (HF) exchange with a short-range semilocal exchange-correlation density functional, and the RSH+MP2 method which adds long-range second-order M{{\o}}ller-Plesset (MP2) correlation. Results on atoms and molecules show that the fractional-charge errors obtained in RSH are much smaller than in the standard Kohn-Sham (KS) scheme applied with semilocal or hybrid approximations, and also generally smaller than in the standard HF method. The RSH+MP2 method tends to have smaller fractional-charge errors than standard MP2 for the most diffuse systems, but larger fractional-charge errors for the more compact systems. Even though the individual contributions to the fractional-spin errors in the H atom coming from the short-range exchange and correlation density-functional approximations are smaller than the corresponding contributions for the full-range exchange and correlation density-functional approximations, RSH gives fractional-spin errors that are larger than in the standard KS scheme and only slightly smaller than in standard HF. Adding long-range MP2 correlation only leads to infinite fractional-spin errors. This work clarifies the successes and limitations of range-separated density-functional theory approaches for eliminating self-interaction and static-correlation errors.

We explore the computation of high-harmonic generation spectra by means of Gaussian basis sets in approaches propagating the time-dependent Schr{\"o}dinger equation. We investigate the efficiency of Gaussian functions specifically designed for the description of the continuum proposed by Kaufmann et al. [J. Phys. B 22, 2223 (1989)]. We assess the range of applicability of this approach by studying the hydrogen atom, i.e. the simplest atom for which "exact" calculations on a grid can be performed. We notably study the effect of increasing the basis set cardinal number, the number of diffuse basis functions, and the number of Gaussian pseudo-continuum basis functions for various laser parameters. Our results show that the latter significantly improve the description of the low-lying continuum states, and provide a satisfactory agreement with grid calculationsfor laser wavelengths $\lambda$0 = 800 and 1064 nm. The Kaufmann continuum functions therefore appear as a promising way of constructing Gaussian basis sets for studying molecular electron dynamics in strong laser fields using time-dependent quantum-chemistry approaches.

We explore the use in quantum Monte Carlo (QMC) of trial wave functions consisting of a Jastrow factor multiplied by a truncated configuration-interaction (CI) expansion in Slater determinants obtained from a CI perturbatively selected iteratively (CIPSI) calculation. In the CIPSI algorithm, the CI expansion is iteratively enlarged by selecting the best determinants using perturbation theory, which provides an optimal and automatic way of constructing truncated CI expansions approaching the full CI limit. We perform a systematic study of variational Monte Carlo (VMC) and fixed-node diffusion Monte Carlo (DMC) total energies of first-row atoms from B to Ne with different levels of optimization of the parameters (Jastrow parameters, coefficients of the determinants, and orbital parameters) in these trial wave functions. The results show that the reoptimization of the coefficients of the determinants in VMC (together with the Jastrow factor) leads to an important lowering of both VMC and DMC total energies, and to their monotonic convergence with the number of determinants. In addition, we show that the reoptimization of the orbitals is also important in both VMC and DMC for the Be atom when using a large basis set. These reoptimized Jastrow-CIPSI wave functions appear as promising, systematically improvable trial wave functions for QMC calculations., Comment: in Molecular Physics, Taylor \& Francis, 2016

We present a range-separated linear-response time-dependent density-functional theory (TDDFT) which combines a density-functional approximation for the short-range response kernel and a frequency-dependent second-order Bethe-Salpeter approximation for the long-range response kernel. This approach goes beyond the adiabatic approximation usually used in linear-response TDDFT and aims at improving the accuracy of calculations of electronic excitation energies of molecular systems. A detailed derivation of the frequency-dependent second-order Bethe-Salpeter correlation kernel is given using many-body Green-function theory. Preliminary tests of this range-separated TDDFT method are presented for the calculation of excitation energies of the He and Be atoms and small molecules (H2, N2, CO2, H2CO, and C2H4). The results suggest that the addition of the long-range second-order Bethe-Salpeter correlation kernel overall slightly improves the excitation energies., Comment: in The Journal of Chemical Physics, American Institute of Physics, 2016

We provide a pedagogical introduction to the two main variants of real-space quantum Monte Carlo methods for electronic-structure calculations: variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC). Assuming no prior knowledge on the subject, we review in depth the Metropolis-Hastings algorithm used in VMC for sampling the square of an approximate wave function, discussing details important for applications to electronic systems. We also review in detail the more sophisticated DMC algorithm within the fixed-node approximation, introduced to avoid the infamous Fermionic sign problem, which allows one to sample a more accurate approximation to the ground-state wave function. Throughout this review, we discuss the statistical methods used for evaluating expectation values and statistical uncertainties. In particular, we show how to estimate nonlinear functions of expectation values and their statistical uncertainties., Comment: Advances in Quantum Chemistry, 2015, Electron Correlation in Molecules -- ab initio Beyond Gaussian Quantum Chemistry, pp.0000

Quantum chemistry methods exploiting density-functional approximations for short-range electron-electron interactions and second-order M{{\o}}ller-Plesset (MP2) perturbation theory for long-range electron-electron interactions have been implemented for periodic systems using Gaussian-type basis functions and the local correlation framework. The performance of these range-separated double hybrids has been benchmarked on a significant set of systems including rare-gas, molecular, ionic, and covalent crystals. The use of spin-component-scaled MP2 for the long-range part has been tested as well. The results show that the value of $\mu$ = 0.5 bohr^{--1} for the range-separation parameter usually used for molecular systems is also a reasonable choice for solids. Overall, these range-separated double hybrids provide a good accuracy for binding energies using basis sets of moderate sizes such as cc-pVDZ and aug-cc-pVDZ.

We consider several spin-unrestricted random-phase approximation (RPA) variants for calculating correlation energies, with and without range separation, and test them on datasets of atomization energies and reaction barrier heights. We show that range separation greatly improves the accuracy of all RPA variants for these properties. Moreover, we show that a RPA variant with exchange, hereafter referred to as RPAx-SO2, first proposed by Sz-abo and Ostlund [A. Szabo and N. S. Ostlund, J. Chem. Phys. 67, 4351 (1977)] in a spin-restricted closed-shell formalism, and extended here to a spin-unrestricted formalism , provides on average the most accurate range-separated RPA variant for atomization energies and reaction barrier heights. Since this range-separated RPAx-SO2 method had already been shown to be among the most accurate range-separated RPA variants for weak intermolecular interactions [J. Toulouse, W. Zhu, A. Savin, G. Jansen, and J. G. {\'A}ngy{\'a}n, J. Chem. Phys. 135, 084119 (2011)], this works confirms range-separated RPAx-SO2 as a promising method for general chemical applications.

In this paper, an alternative method to range-separated linear-response time-dependent density-functional theory and perturbation theory is proposed to improve the estimation of the energies of a physical system from the energies of a partially interacting system. Starting from the analysis of the Taylor expansion of the energies of the partially interacting system around the physical system, we use an extrapolation scheme to improve the estimation of the energies of the physical system at an intermediate point of the range-separated or linear adiabatic connection where either the electron--electron interaction is scaled or only the long-range part of the Coulomb interaction is included. The extrapolation scheme is first applied to the range-separated energies of the helium and beryllium atoms and of the hydrogen molecule at its equilibrium and stretched geometries. It improves significantly the convergence rate of the energies toward their exact limit with respect to the range-separation parameter. The range-separated extrapolation scheme is compared with a similar approach for the linear adiabatic connection, highlighting the relative strengths and weaknesses of each approach., Comment: to appear in Phys. Rev. A

We present a study of the variation of total energies and excitationenergies along a range-separated adiabatic connection. This connectionlinks the non-interacting Kohn-Sham electronic system to the physicalinteracting system by progressively switching on theelectron-electron interactions whilst simultaneously adjusting aone-electron effective potential so as to keep the ground-statedensity constant. The interactions are introduced in arange-dependent manner, first introducing predominantly long-range,and then all-range, interactions as the physical system is approached,as opposed to the conventional adiabatic connection where theinteractions are introduced by globally scaling the standard Coulomb interaction.Reference data are reported for the He and Be atoms and the H2molecule, obtained by calculating the short-range effective potentialat the full configuration-interaction level using Lieb'sLegendre-transform approach. As the strength of the electron-electroninteractions increases, the excitation energies, calculated for thepartially interacting systems along the adiabatic connection, offerincreasingly accurate approximations to the exact excitation energies.Importantly, the excitation energies calculated at an intermediatepoint of the adiabatic connection are much better approximations tothe exact excitation energies than are the corresponding Kohn-Shamexcitation energies. This is particularly evident in situationsinvolving strong static correlation effects and states with multipleexcitation character, such as the dissociating H2 molecule. Theseresults highlight the utility of long-range interacting referencesystems as a starting point for the calculation of excitation energiesand are of interest for developing and analyzing practical approximaterange-separated density-functional methodologies.

We explore the possibility of calculating electronic excited states by using perturbation theory along a range-separated adiabatic connection. Starting from the energies of a partially interacting Hamiltonian, a first-order correction is defined with two variants of perturbation theory: a straight-forward perturbation theory, and an extension of the G{\"o}rling--Levy one that has the advantage of keeping the ground-state density constant at each order in the perturbation. Only the first, simpler, variant is tested here on the helium and beryllium atoms and on the dihydrogene molecule. The first-order correction within this perturbation theory improves significantly the total ground-and excited-state energies of the different systems. However, the excitation energies are mostly deterio-rated with respect to the zeroth-order ones, which may be explained by the fact that the ionization energy is no longer correct for all interaction strengths. The second variant of the perturbation theory should improve these results but has not been tested yet along the range-separated adiabatic connection.