Search

In a recent work we have proposed an original analytic expression for the partition function of the quartic oscillator. This partition function, which has a simple and compact form with {\it no adjustable parameters}, reproduces some key mathematical properties of the exact partition function and provides free energies accurate to a few percent over a wide range of temperatures and coupling constants. In this work, we present the derivation of the energy spectrum of this model. We also generalize our previous study limited to the quartic oscillator to the case of a general anharmonic oscillator. Numerical application for a potential of the form $V(x)=\frac{\omega^2}{2} x^2 + g x^{2m}$ show that the energy levels are obtained with a relative error of about a few percent, a precision which we consider to be quite satisfactory given the simplicity of the model, the absence of adjustable parameters, and the negligible computational cost., Comment: 7 pages. Revised version

In this paper we examine the numerical approximation of the limiting invariant measure associated with Feynman-Kac formulae. These are expressed in a discrete time formulation and are associated with a Markov chain and a potential function. The typical application considered here is the computation of eigenvalues associated with non-negative operators as found, for example, in physics or particle simulation of rare-events. We focus on a novel \emph{lagged} approximation of this invariant measure, based upon the introduction of a ratio of time-averaged Feynman-Kac marginals associated with a positive operator iterated $l \in\mathbb{N}$ times; a lagged Feynman-Kac formula. This estimator and its approximation using Diffusion Monte Carlo (DMC) have been extensively employed in the physics literature. In short, DMC is an iterative algorithm involving $N\in\mathbb{N}$ particles or walkers simulated in parallel, that undergo sampling and resampling operations. In this work, it is shown that for the DMC approximation of the lagged Feynman-Kac formula, one has an almost sure characterization of the $\mathbb{L}_1$-error as the time parameter (iteration) goes to infinity and this is at most of $\mathcal{O}(\exp\{-\kappa l\}/N)$, for $\kappa>0$. In addition a non-asymptotic in time, and time uniform $\mathbb{L}_1-$bound is proved which is $\mathcal{O}(l/\sqrt{N})$. We also prove a novel central limit theorem to give a characterization of the exact asymptotic in time variance. This analysis demonstrates that the strategy used in physics, namely, to run DMC with $N$ and $l$ small and, for long time enough, is mathematically justified. Our results also suggest how one should choose $N$ and $l$ in practice. We emphasize that these results are not restricted to physical applications; they have broad relevance to the general problem of particle simulation of the Feynman-Kac formula., Comment: 22 pages, 3 figures

Recently, a new distributed implementation of the full configuration interaction (FCI) method has been reported [Gao et al. J. Chem Theory Comput. 2024, 20, 1185]. Thanks to a hybrid parallelization scheme, the authors were able to compute the exact energy of propane (\ce{C3H8}) in the minimal basis STO-3G. This formidable task involves handling an active space of 26 electrons in 23 orbitals or a Hilbert space of \SI{1.3d12} determinants. This is, by far, the largest FCI calculation reported to date. Here, we illustrate how, from a general point of view, selected configuration interaction (SCI) can achieve microhartree accuracy at a fraction of the computational and memory cost, via a sparse exploration of the FCI space. The present SCI calculations are performed with the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) algorithm, as implemented in a determinant-driven way in the \textsc{quantum package} software. The present study reinforces the common wisdom that among the exponentially large number of determinants in the FCI space, only a tiny fraction of them significantly contribute to the energy. More importantly, it demonstrates the feasibility of achieving comparable accuracy using more reasonable and sustainable computational resources, hence reducing the ever-growing carbon footprint of computational chemistry., Comment: 5 pages, 2 figures (supporting information available)

The Diffusion Monte Carlo method with constant number of walkers, also called Stochastic Reconfiguration as well as Sequential Monte Carlo, is a widely used Monte Carlo methodology for computing the ground-state energy and wave function of quantum systems. In this study, we present the first mathematically rigorous analysis of this class of stochastic methods on non necessarily compact state spaces, including linear diffusions evolving in quadratic absorbing potentials, yielding what seems to be the first result of this type for this class of models. We present a novel and general mathematical framework with easily checked Lyapunov stability conditions that ensure the uniform-in-time convergence of Diffusion Monte Carlo estimates towards the top of the spectrum of Schr\"odinger operators. For transient free evolutions, we also present a divergence blow up of the estimates w.r.t. the time horizon even when the asymptotic fluctuation variances are uniformly bounded. We also illustrate the impact of these results in the context of generalized coupled quantum harmonic oscillators with non necessarily reversible nor stable diffusive particle and a quadratic energy absorbing well associated with a semi-definite positive matrix force.

In this work an approximate analytic expression for the quantum partition function of the quartic oscillator described by the potential $V(x) = \frac{1}{2} \omega^2 x^2 + g x^4$ is presented. Using a path integral formalism, the exact partition function is approximated by the partition function of a harmonic oscillator with an effective frequency depending both on the temperature and coupling constant $g$. By invoking a Principle of Minimal Sensitivity (PMS) of the path integral to the effective frequency, we derive a mathematically well-defined analytic formula for the partition function. Quite remarkably, the formula reproduces qualitatively and quantitatively the key features of the exact partition function. The free energy is accurate to a few percent over the entire range of temperatures and coupling strengths $g$. Both the harmonic ($g\rightarrow 0$) and classical (high-temperature) limits are exactly recovered. The divergence of the power series of the ground-state energy at weak coupling, characterized by a factorial growth of the perturbational energies, is reproduced as well as the functional form of the strong-coupling expansion along with accurate coefficients. Explicit accurate expressions for the ground- and first-excited state energies, $E_0(g)$ and $E_1(g)$ are also presented., Comment: 15 pages, 4 figures. Minor revisions

We report ground- and excited-state dipole moments and oscillator strengths (computed in different ``gauges'' or representations) of full configuration interaction (FCI) quality using the selected configuration interaction method known as \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI). Thanks to a set encompassing 35 ground- and excited-state properties computed in 11 small molecules, the present near-FCI estimates allow us to assess the accuracy of high-order coupled-cluster (CC) calculations including up to quadruple excitations. In particular, we show that incrementing the excitation degree of the CC expansion (from CCSD to CCSDT or from CCSDT to CCSDTQ) reduces the average error with respect to the near-FCI reference values by approximately one order of magnitude., Comment: 12 pages, 8 figures (supporting information available)

The sampling of the configuration space in diffusion Monte Carlo (DMC) is done using walkers moving randomly. In a previous work on the Hubbard model [\href{https://doi.org/10.1103/PhysRevB.60.2299}{Assaraf et al.~Phys.~Rev.~B \textbf{60}, 2299 (1999)}], it was shown that the probability for a walker to stay a certain amount of time in the same state obeys a Poisson law and that the on-state dynamics can be integrated out exactly, leading to an effective dynamics connecting only different states. Here, we extend this idea to the general case of a walker trapped within domains of arbitrary shape and size. The equations of the resulting effective stochastic dynamics are derived. The larger the average (trapping) time spent by the walker within the domains, the greater the reduction in statistical fluctuations. A numerical application to the Hubbard model is presented. Although this work presents the method for finite linear spaces, it can be generalized without fundamental difficulties to continuous configuration spaces., Comment: 14 pages, 4 figures

We study beryllium dihydride (BeH$_2$) and acetylene (C$_2$H$_2$) molecules using real-space diffusion Monte Carlo (DMC) method. The molecules serve as perhaps the simplest prototypes that illustrate the difficulties with biases in the fixed-node DMC calculations that might appear with the use of effective core potentials (ECPs) or other nonlocal operators. This is especially relevant for the recently introduced correlation consistent ECPs (ccECPs) for $2s2p$ elements. Corresponding ccECPs exhibit deeper potential functions due to higher fidelity to all-electron counterparts, which could lead to larger local energy fluctuations. We point out that the difficulties stem from issues that are straightforward to address by upgrades of basis sets, use of T-moves for nonlocal terms, inclusion of a few configurations into the trial function and similar. The resulting accuracy corresponds to the ccECP target (chemical accuracy) and it is in consistent agreement with independent correlated calculations. Further possibilities for upgrading the reliability of the DMC algorithm and considerations for better adapted and more robust Jastrow factors are discussed as well., Comment: 13 pages, 5 figures

Following our recent work on the benzene molecule [\href{https://doi.org/10.1063/5.0027617}{J.~Chem.~Phys.~\textbf{153}, 176101 (2020)}], itself motivated by the blind challenge of Eriksen \textit{et al.} [\href{https://doi.org/10.1021/acs.jpclett.0c02621}{J.~Phys.~Chem.~Lett.~\textbf{11}, 8922 (2020)}] on the same system, we report accurate full configuration interaction (FCI) frozen-core correlation energy estimates for twelve five- and six-membered ring molecules in the standard correlation-consistent double-$\zeta$ Dunning basis set (cc-pVDZ). Our FCI correlation energy estimates, with estimated error smaller than 1 millihartree, are based on energetically optimized-orbital selected configuration interaction (SCI) calculations performed with the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) 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{\o}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), are also investigated., Comment: 13 pages, 5 figures

The pair coupled cluster doubles (pCCD) method (where the excitation manifold is restricted to electron pairs) has a series of interesting features. Among others, it provides ground-state energies very close to what is obtained with doubly-occupied configuration interaction (DOCI), but with polynomial cost (compared with the exponential cost of the latter). Here, we address whether this similarity holds for excited states, by exploring the symmetric dissociation of the linear \ce{H4} molecule. When ground-state Hartree-Fock (HF) orbitals are employed, pCCD and DOCI excited-state energies do not match, a feature that is assigned to the poor HF reference. In contrast, by optimizing the orbitals at the pCCD level (oo-pCCD) specifically for each excited state, the discrepancies between pCCD and DOCI decrease by one or two orders of magnitude. Therefore, the pCCD and DOCI methodologies still provide comparable energies for excited states, but only if suitable, state-specific orbitals are adopted. We also assessed whether a pCCD approach could be used to directly target doubly-excited states, without having to resort to the equation-of-motion (EOM) formalism. In our $\Delta$oo-pCCD model, excitation energies were extracted from the energy difference between separate oo-pCCD calculations for the ground state and the targeted excited state. For a set comprising the doubly-excited states of \ce{CH+}, \ce{BH}, nitroxyl, nitrosomethane, and formaldehyde, we found that $\Delta$oo-pCCD provides quite accurate excitation energies, with root mean square deviations (with respect to full configuration interaction results) lower than CC3 and comparable to EOM-CCSDT, two methods with much higher computational cost., Comment: 12 pages, 4 figures

We describe our efforts of the past few years to create a large set of more than 500 highly-accurate vertical excitation energies of various natures ($\pi \to \pi^*$, $n \to \pi^*$, double excitation, Rydberg, singlet, doublet, triplet, etc) in small- and medium-sized molecules. These values have been obtained using an incremental strategy which consists in combining high-order coupled cluster and selected configuration interaction calculations using increasingly large diffuse basis sets in order to reach high accuracy. One of the key aspect of the so-called QUEST database of vertical excitations is that it does not rely on any experimental values, avoiding potential biases inherently linked to experiments and facilitating theoretical cross comparisons. Following this composite protocol, we have been able to produce theoretical best estimate (TBEs) with the aug-cc-pVTZ basis set for each of these transitions, as well as basis set corrected TBEs (i.e., near the complete basis set limit) for some of them. The TBEs/aug-cc-pVTZ have been employed to benchmark a large number of (lower-order) wave function methods such as CIS(D), ADC(2), CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, ADC(3), CC3, NEVPT2, and others (including spin-scaled variants). In order to gather the huge amount of data produced during the QUEST project, we have created a website [https://lcpq.github.io/QUESTDB_website] where one can easily test and compare the accuracy of a given method with respect to various variables such as the molecule size or its family, the nature of the excited states, the type of basis set, etc. We hope that the present review will provide a useful summary of our effort so far and foster new developments around excited-state methods., Comment: 44 pages, 5 figures, 4 Tables, supplementary information available at https://doi.org/10.5281/zenodo.4297012

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for \textit{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) must be used. 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.

Quantum Monte Carlo (QMC) methods are some of the most accurate methods for simulating correlated electronic systems. We investigate the compatibility, strengths and weaknesses of two such methods, namely, diffusion Monte Carlo (DMC) and auxiliary-field quantum Monte Carlo (AFQMC). The multi-determinant trial wave functions employed in both approaches are generated using the configuration interaction using a perturbative selection made iteratively (CIPSI) technique. Complete basis set full configuration interaction (CBS-FCI) energies estimated with CIPSI are used as a reference in this comparative study between DMC and AFQMC. By focusing on a set of canonical finite size solid state systems, we show that both QMC methods can be made to systematically converge towards the same energy once basis set effects and systematic biases have been removed. AFQMC shows a much smaller dependence on the trial wavefunction than DMC while simultaneously exhibiting a much larger basis set dependence. We outline some of the remaining challenges and opportunities for improving these approaches.

A Monte Carlo method for evaluating multi-center two-electron-repulsion integrals over any type of orbitals (Slater, Sturmian, finite-range, numerical, etc.) is presented. The approach is based on a simple and universal (orbital-independent) gaussian sampling of the two-electron configuration space and on the use of efficient zero-variance Monte Carlo estimators. Quite remarkably, it is shown that the high level of accuracy required on two-electron integrals to make Hartree-Fock (HF) and configuration interaction (CI) calculations feasible can be achieved. A first zero-variance estimator is built by introducing a gaussian approximation of the orbitals and by evaluating the two-electron integrals using a correlated sampling scheme for the difference between exact and approximate orbitals. A second one is based on the introduction of a general coordinate transformation. The price to pay for this simple and general Monte Carlo scheme is the high computational cost required. However, we argue that the great simplicity of the algorithm, its embarrassingly parallel nature, its ideal adaptation to modern computational platforms and, most importantly, the possibility of using more compact and physically meaningful basis sets make nevertheless the method attractive. HF and near full CI (FCI) calculations using Slater-type orbitals (STO) are reported for Be, CH4 and [H2N(CH)NH2]^+ (a simple model of cyanine). To the best of our knowledge, our largest FCI calculation involving 18 active electrons distributed among 90 orbitals for the cyanine molecule, is the most extensive molecular calculation performed so far using pure STO orbitals (no gaussian approximation, even for the challenging four-center two-electron integrals)., Comment: 11 pages, 2 figures

Due to their diverse nature, the faithful description of excited states within electronic structure theory methods remains one of the grand challenges of modern theoretical chemistry. Quantum Monte Carlo (QMC) methods have been applied very successfully to ground state properties but still remain generally less effective than other non-stochastic methods for electronically excited states. Nonetheless, we have recently reported accurate excitation energies for small organic molecules at the fixed-node diffusion Monte Carlo (FN-DMC) within a Jastrow-free QMC protocol relying on a deterministic and systematic construction of nodal surfaces using the selected configuration interaction (sCI) algorithm known as CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively). Albeit highly accurate, these all-electron calculations are computationally expensive due to the presence of core electrons. One very popular approach to remove these chemically-inert electrons from the QMC simulation is to introduce pseudopotentials (also known as effective core potentials). Taking the water molecule as an example, we investigate the influence of Burkatzki-Filippi-Dolg (BFD) pseudopotentials and their associated basis sets on vertical excitation energies obtained with sCI and FN-DMC methods. Although these pseudopotentials are known to be relatively safe for ground state properties, we evidence that special care may be required if one strives for highly accurate vertical transition energies. Indeed, comparing all-electron and valence-only calculations, we show that using pseudopotentials with the associated basis sets can induce differences of the order of 0.05 eV on the excitation energies. Fortunately, a reasonable estimate of this shift can be estimated at the sCI level., Comment: 7 pages, 3 figures (SI available)

\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 describe a method for imposing the correct electron-nucleus (e-n) cusp in molecular orbitals expanded as a linear combination of (cuspless) Gaussian basis functions. Enforcing the e-n cusp in trial wave functions is an important asset in quantum Monte Carlo calculations as it significantly reduces the variance of the local energy during the Monte Carlo sampling. In the method presented here, the Gaussian basis set is augmented with a small number of Slater basis functions. Note that, unlike other e-n cusp correction schemes, the presence of the Slater function is not limited to the vicinity of the nuclei. Both the coefficients of these cuspless Gaussian and cusp-correcting Slater basis functions may be self-consistently optimized by diagonalization of an orbital-dependent effective Fock operator. Illustrative examples are reported for atoms (\ce{H}, \ce{He} and \ce{Ne}) as well as for a small molecular system (\ce{BeH2}). For the simple case of the \ce{He} atom, we observe that, with respect to the cuspless version, the variance is reduced by one order of magnitude by applying our cusp-corrected scheme., Comment: 23 pages, 5 figures

Excited states exhibiting double excitation character are notoriously difficult to model using conventional single-reference methods, such as adiabatic time-dependent density-functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC). In the present work, we provide accurate reference excitation energies for transitions involving a substantial amount of double excitation using a series of increasingly large diffuse-containing atomic basis sets. Our set gathers 20 vertical transitions from 14 small- and medium-size molecule. Depending on the size of the molecule, selected configuration interaction (sCI) and/or multiconfigurational calculations are performed in order to obtain reliable estimates of the vertical transition energies. In addition, coupled cluster approaches including at least contributions from iterative triples are assessed. Our results clearly evidence that the error in CC methods is intimately related to the amount of double excitation character of the transition. For "pure" double excitations (i.e. for transitions which do not mix with single excitations), the error in CC3 can easily reach 1 eV, while it goes down to few tenths of an eV for more common transitions (like in \emph{trans}-butadiene) involving a significant amount of singles. As expected, CC approaches including quadruples yield highly accurate results for any type of transitions. The quality of the excitation energies obtained with multiconfigurational methods is harder to predict. We have found that the overall accuracy of these methods is highly dependent of both the system and the selected active space. The inclusion of the $\sigma$ and $\sigma^*$ orbitals in the active space, even for transitions involving mostly $\pi$ and $\pi^*$ orbitals, is mandatory in order to reach high accuracy. A theoretical best estimate (TBE) is reported for each transition., Comment: 16 pages, 3 figures

Striving to define very accurate vertical transition energies, we perform both high-level coupled cluster (CC) calculations (up to CCSDTQP) and selected configuration interaction (sCI) calculations (up to several millions of determinants) for 18 small compounds. By systematically increasing the order of the CC expansion, the number of determinants in the CI expansion as well as the size of the one-electron basis set, we have been able to reach near full CI (FCI) quality transition energies. These calculations are carried out on CC3/aug-cc-pVTZ geometries, using a series of increasingly large atomic basis sets systematically including diffuse functions. In this way, we define a list of 110 transition energies for states of various characters to be used as references for further calculations. Benchmark transition energies are provided at the aug-pVTZ level as well as with additional basis set corrections, in order to obtain results close to the complete basis set limit. These reference data are used to benchmark a series of twelve excited-state wave function methods accounting for double and triple contributions, namely ADC(2), ADC(3), CIS(D), CIS(D$_\infty$), CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, CC3, CCSDT and CCSDTQ. It turns out that CCSDTQ yields a negligible difference with the extrapolated CI values with a mean absolute error as small as $0.01$ eV, whereas the coupled cluster approaches including iterative triples are also very accurate (mean absolute error of $0.03$ eV). Consequently, CCSDT-3 and CC3 can be used to define reliable benchmarks. This observation does not hold for ADC(3) that delivers quite large errors for this set of small compounds, with a clear tendency to overcorrect its second-order version, ADC(2). Finally, we discuss the possibility to use basis set extrapolation approaches so as to tackle more easily larger compounds., Comment: 57 pages, 2 figures

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 CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) 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 coworkers. In particular, the multi-state algorithm removes the storage bottleneck of the effective Hamiltonian via a low-rank factorization of the dressing matrix., Comment: 6 pages, 1 figure. Supplementary material available

Quantum Monte Carlo (QMC) is a stochastic method which 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-$\zeta$) basis sets., Comment: 8 pages, 3 figures

In diffusion Monte Carlo (DMC) methods, the nodes (or zeroes) of the trial wave function dictate the magnitude of the fixed-node (FN) error. Within standard DMC implementations, they emanate from short multideterminant expansions, \textit{stochastically} optimized in the presence of a Jastrow factor. Here, following a recent proposal, we follow an alternative route by considering the nodes of selected configuration interaction (sCI) expansions built with the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm. In contrast to standard implementations, these nodes can be \textit{systematically} and \textit{deterministically} improved by increasing the size of the sCI expansion. The present methodology is used to investigate the properties of the transition metal sulfide molecule FeS. This apparently simple molecule has been shown to be particularly challenging for electronic structure theory methods due to the proximity of two low-energy quintet electronic states of different spatial symmetry. In particular, we show that, at the triple-zeta basis set level, all sCI results --- including those extrapolated at the full CI (FCI) limit --- disagree with experiment, yielding an electronic ground state of $^{5}\Sigma^+$ symmetry. Performing FN-DMC simulation with sCI nodes, we show that the correct $^{5}\Delta$ ground state is obtained if sufficiently large expansions are used. Moreover, we show that one can systematically get accurate potential energy surfaces and reproduce the experimental dissociation energy as well as other spectroscopic constants., Comment: 8 pages, 2 figure and 4 tables

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 MC 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 $\sim 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 F$_2$ and Cr$_2$ molecules. Even for the largest case corresponding to the Cr$_2$ 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 \times 10^7$ determinants., Comment: 8 pages, 5 figures

Several aspects of the recently proposed DMC-CIPSI approach consisting in using selected Configuration Interaction (SCI) approaches such as CIPSI (Configuration Interaction using a Perturbative Selection done Iteratively) to build accurate nodes for diffusion Monte Carlo (DMC) calculations are presented and discussed. The main ideas are illustrated with a number of calculations for diatomics molecules and for the benchmark G1 set.

All-electron Fixed-node Diffusion Monte Carlo (FN-DMC) calculations for the nonrelativistic ground-state energy of the water molecule at equilibrium geometry are presented. The determinantal part of the trial wavefunction is obtained from a perturbatively selected Configuration Interaction calculation (CIPSI method) including up to about 1.4 million of determinants. Calculations are made using the cc-pCV$n$Z family of basis sets, with $n$=2 to 5. In contrast with most QMC works no re-optimization of the determinantal part in presence of a Jastrow is performed. For the largest cc-pCV5Z basis set the lowest upper bound for the ground-state energy reported so far of -76.43744(18) is obtained. The fixed-node energy is found to decrease regularly as a function of the cardinal number $n$ and the Complete Basis Set limit (CBS) associated with {\it exact nodes} is easily extracted. The resulting energy of -76.43894(12) -in perfect agreement with the best experimentally derived value- is the most accurate theoretical estimate reported so far. We emphasize that employing selected CI nodes of increasing quality in a given family of basis sets may represent a simple, deterministic, reproducible, and systematic way of controlling the fixed-node error in DMC., Comment: 5 pages, 1 figure

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 $N_{\rm det}^\sigma$ ($\sigma=\uparrow,\downarrow$) with a non-negligible weight in the expansion is of order ${\cal O}(\sqrt{N_{\rm det}})$. We show that a careful implementation of the calculation of the $N_{\rm det}$-dependent contributions can make this step negligible enough so that in practice the algorithm scales as the total number of unique spin-specific determinants, $\; N_{\rm det}^\uparrow + N_{\rm det}^\downarrow$, over a wide range of total number of determinants (here, $N_{\rm det}$ 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 $\sim$ 400 compared to a single-determinant calculation are shown to be feasible., Comment: 9 pages, 3 figures

The potential energy curve of the F$_2$ molecule is calculated with Fixed-Node Diffusion Monte Carlo (FN-DMC) using Configuration Interaction (CI)-type trial wavefunctions. To keep the number of determinants reasonable (the first and second derivatives of the trial wavefunction need to be calculated at each step of FN-DMC), the CI expansion is restricted to those determinants that contribute the most to the total energy. The selection of the determinants is made using the so-called CIPSI approach (Configuration Interaction using a Perturbative Selection made Iteratively). Quite remarkably, the nodes of CIPSI wavefunctions are found to be systematically improved when increasing the number of selected determinants. To reduce the non-parallelism error of the potential energy curve a scheme based on the use of a $R$-dependent number of determinants is introduced. Numerical results show that improved FN-DMC energy curves for the F$_2$ molecule are obtained when employing CIPSI trial wavefunctions. Using the Dunning's cc-pVDZ basis set the FN-DMC energy curve is of a quality similar to that obtained with FCI/cc-pVQZ. A key advantage of using selected CI in FN-DMC is the possibility of improving nodes in a systematic and automatic way without resorting to a preliminary multi-parameter stochastic optimization of the trial wavefunction performed at the Variational Monte Carlo level as usually done in FN-DMC., Comment: 16 pages, 15 figures

We present a comparative study of the spatial distribution of the spin density (SD) of the ground state of CuCl2 using Density Functional Theory (DFT), quantum Monte Carlo (QMC), and post-Hartree-Fock wavefunction theory (WFT). A number of studies have shown that an accurate description of the electronic structure of the lowest-lying states of this molecule is particularly challenging due to the interplay between the strong dynamical correlation effects in the 3d shell of the copper atom and the delocalization of the 3d hole over the chlorine atoms. It is shown here that qualitatively different results for SD are obtained from these various quantum-chemical approaches. At the DFT level, the spin density distribution is directly related to the amount of Hartree-Fock exchange introduced in hybrid functionals. At the QMC level, Fixed-node Diffusion Monte Carlo (FN-DMC) results for SD are strongly dependent on the nodal structure of the trial wavefunction employed (here, Hartree-Fock or Kohn-Sham with a particular amount of HF exchange) : in the case of this open-shell system, the 3N -dimensional nodes are mainly determined by the 3-dimensional nodes of the singly occupied molecular orbital. Regarding wavefunction approaches, HF and CASSCF lead to strongly localized spin density on the copper atom, in sharp contrast with DFT. To get a more reliable description and shed some light on the connections between the various theoretical descriptions, Full CI-type (FCI) calculations are performed. To make them feasible for this case a perturbatively selected CI approach generating multi-determinantal expansions of reasonable size and a small tractable basis set are employed. Although semi-quantitative, these near-FCI calculations allow to clarify how the spin density distribution evolves upon inclusion of dynamic correlation effects. A plausible scenario about the nature of the SD is proposed., Comment: 13 pages, 12 Figures

Various strategies to implement efficiently QMC simulations for large chemical systems are presented. These include: i.) the introduction of an efficient algorithm to calculate the computationally expensive Slater matrices. This novel scheme is based on the use of the highly localized character of atomic Gaussian basis functions (not the molecular orbitals as usually done), ii.) the possibility of keeping the memory footprint minimal, iii.) the important enhancement of single-core performance when efficient optimization tools are employed, and iv.) the definition of a universal, dynamic, fault-tolerant, and load-balanced computational framework adapted to all kinds of computational platforms (massively parallel machines, clusters, or distributed grids). These strategies have been implemented in the QMC=Chem code developed at Toulouse and illustrated with numerical applications on small peptides of increasing sizes (158, 434, 1056 and 1731 electrons). Using 10k-80k computing cores of the Curie machine (GENCI-TGCC-CEA, France) QMC=Chem has been shown to be capable of running at the petascale level, thus demonstrating that for this machine a large part of the peak performance can be achieved. Implementation of large-scale QMC simulations for future exascale platforms with a comparable level of efficiency is expected to be feasible.

We consider the use in quantum Monte Carlo calculations of two types of valence bond wave functions based on strictly localized active orbitals, namely valence bond self-consistent-field (VBSCF) and breathing-orbital valence bond (BOVB) wave functions. Complemented by a Jastrow factor, these Jastrow-Valence-Bond wave functions are tested by computing the equilibrium well depths of the four diatomic molecules C$_2$, N$_2$, O$_2$, and F$_2$ in both variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC). We show that it is possible to design compact wave functions based on chemical grounds that are capable of describing both static and dynamic electron correlation. These wave functions can be systematically improved by inclusion of valence bond structures corresponding to additional bonding patterns.

We report state-of-the-art quantum Monte Carlo calculations of the singlet $n \to \pi^*$ (CO) vertical excitation energy in the acrolein molecule, extending the recent study of Bouab\c{c}a {\it et al.} [J. Chem. Phys. {\bf 130}, 114107 (2009)]. We investigate the effect of using a Slater basis set instead of a Gaussian basis set, and of using state-average versus state-specific complete-active-space (CAS) wave functions, with or without reoptimization of the coefficients of the configuration state functions (CSFs) and of the orbitals in variational Monte Carlo (VMC). It is found that, with the Slater basis set used here, both state-average and state-specific CAS(6,5) wave functions give an accurate excitation energy in diffusion Monte Carlo (DMC), with or without reoptimization of the CSF and orbital coefficients in the presence of the Jastrow factor. In contrast, the CAS(2,2) wave functions require reoptimization of the CSF and orbital coefficients to give a good DMC excitation energy. Our best estimates of the vertical excitation energy are between 3.86 and 3.89 eV., Comment: 6 pages, 1 figure, 2 tables, to appear in Progress in Theoretical Chemistry and Physics

In this work we present a detailed study of the Fermion Monte Carlo algorithm (FMC), a recently proposed stochastic method for calculating fermionic ground-state energies [M.H. Kalos and F. Pederiva, Phys. Rev. Lett. vol. 85, 3547 (2000)]. A proof that the FMC method is an exact method is given. In this work the stability of the method is related to the difference between the lowest (bosonic-type) eigenvalue of the FMC diffusion operator and the exact fermi energy. It is shown that within a FMC framework the lowest eigenvalue of the new diffusion operator is no longer the bosonic ground-state eigenvalue as in standard exact Diffusion Monte Carlo (DMC) schemes but a modified value which is strictly greater. Accordingly, FMC can be viewed as an exact DMC method built from a correlated diffusion process having a reduced Bose-Fermi gap. As a consequence, the FMC method is more stable than any transient method (or nodal release-type approaches). We illustrate the various ideas presented in this work with calculations performed on a very simple model having only nine states but a full sign problem. Already for this toy model it is clearly seen that FMC calculations are inherently uncontrolled., Comment: 49 pages with 4 postscript figures

An alternative Monte Carlo estimator for the one-body density rho(r) is presented. This estimator has a simple form and can be readily used in any type of Monte Carlo simulation. Comparisons with the usual regularization of the delta-function on a grid show that the statistical errors are greatly reduced. Furthermore, our expression allows accurate calculations of the density at any point in space, even in the regions never visited during the Monte Carlo simulation. The method is illustrated with the computation of accurate Variational Monte Carlo electronic densities for the Helium atom (1D curve) and for the water dimer (3D grid containing up to 51x51x51=132651 points)., Comment: 12 pages with 3 postscript figures

We propose a new algorithm for sampling the $N$-body density $|\Psi({\bf R})|^2/\int_{\mathbb{R}^{3N}} |\Psi|^2$ in the Variational Monte Carlo (VMC) framework. This algorithm is based upon a modified Ricci-Ciccotti discretization of the Langevin dynamics in the phase space $({\bf R},{\bf P})$ improved by a Metropolis acceptation/rejection step. We show through some representative numerical examples (Lithium, Fluorine and Copper atoms, and phenol molecule), that this algorithm is superior to the standard sampling algorithm based on the biased random walk (importance sampling)., Comment: 23 pages

A simple and stable method for computing accurate expectation values of observable with Variational Monte Carlo (VMC) or Diffusion Monte Carlo (DMC) algorithms is presented. The basic idea consists in replacing the usual ``bare'' estimator associated with the observable by an improved or ``renormalized'' estimator. Using this estimator more accurate averages are obtained: Not only the statistical fluctuations are reduced but also the systematic error (bias) associated with the approximate VMC or (fixed-node) DMC probability densities. It is shown that improved estimators obey a Zero-Variance Zero-Bias (ZVZB) property similar to the usual Zero-Variance Zero-Bias property of the energy with the local energy as improved estimator. Using this property improved estimators can be optimized and the resulting accuracy on expectation values may reach the remarkable accuracy obtained for total energies. As an important example, we present the application of our formalism to the computation of forces in molecular systems. Calculations of the entire force curve of the H$_2$,LiH, and Li$_2$ molecules are presented. Spectroscopic constants $R_e$ (equilibrium distance) and $\omega_e$ (harmonic frequency) are also computed. The equilibrium distances are obtained with a relative error smaller than 1%, while the harmonic frequencies are computed with an error of about 10%.

Using a {\it collective} Monte Carlo algorithm we study the low-temperature and long-distance properties of two systems of two-dimensional classical tops. Both systems have the same spin-wave dynamics (low-temperature behavior) as a large class of Heisenberg frustrated spin systems. They are constructed so that to differ only by their topological properties. The spin-stiffnesses for the two systems of tops are calculated for different temperatures and different sizes of the sample. This allows to investigate the role of topological defects in frustrated spin systems. Comparisons with Renormalization Group results based on a Non Linear Sigma model approach and with the predictions of some simple phenomenological model taking into account the topological excitations are done., Comment: RevTex, 25 pages, 14 figures, Minor changes, final version. To appear in Phys.Rev.B

We present a general approach to greatly increase at little cost the efficiency of Monte Carlo algorithms. To each observable to be computed we associate a renormalized observable (improved estimator) having the same average but a different variance. By writing down the zero-variance condition a fundamental equation determining the optimal choice for the renormalized observable is derived (zero-variance principle for each observable separately). We show, with several examples including classical and quantum Monte Carlo calculations, that the method can be very powerful., Comment: 9 pages, Latex, to appear in Phys. Rev. Lett

We investigate the metal-insulator transition of the one-dimensional SU(N) Hubbard model for repulsive interaction. Using the bosonization approach a Mott transition in the charge sector at half-filling (k_F=\pi/Na_0) is conjectured for N > 2. Expressions for the charge and spin velocities as well as for the Luttinger liquid parameters and some correlation functions are given. The theoretical predictions are compared with numerical results obtained with an improved zero-temperature quantum Monte Carlo approach. The method used is a generalized Green's function Monte Carlo scheme in which the stochastic time evolution is partially integrated out. Very accurate results for the gaps, velocities, and Luttinger liquid parameters as a function of the Coulomb interaction U are given for the cases N=3 and N=4. Our results strongly support the existence of a Mott-Hubbard transition at a {\it non-zero} value of the Coulomb interaction. We find $U_c \sim 2.2$ for N=3 and $U_c \sim 2.8$ for N=4., Comment: 22 pages, 9 Figs

The Hubbard model is exactly solved for two particles with opposite spins on d-dimensional hypercubes. It is shown that the spectrum can be separated into two parts: a trivial (U-independent) part resulting from symmetries of hypercubes and a non-trivial part expressed as a single-impurity problem on a set of finite chains of size $d^\prime+1$ ($d^\prime \le d$). The exact expression for the one-particle Green's function is given. Finally, we discuss the extension of these results to standard hypercubic lattices with periodic boundary conditions., Comment: 5 pages, Latex, with 2 eps figures

The low-energy physics of the one-dimensional Pair-Hopping (PH) and attractive Hubbard models are expected to be similar. Based on numerical calculations on small chains, several authors have recently challenged this idea and predicted the existence of a phase transition at half-filling and finite positive coupling for the pair-hopping model. We re-examine the controversy by making systematic comparisons between numerical results obtained for the PH and attractive Hubbard models. To do so, we have calculated the Luttinger parameters (spin and charge velocities, stiffnesses, etc...) of the two models using both the Density Matrix Renormalization Group method for large systems and Lancz\'os calculations with twisted boundary conditions for smaller systems. Although most of our results confirm that both models are very similar we have found some important differences in the spin properties for the small sizes considered by previous numerical studies (6-12 sites). However, we show that these differences disappear at larger sizes (14-42 sites) when sufficiently accurate eigenstates are considered. Accordingly, our results strongly suggest that the ground-state phase transition previously found for small systems is a finite size artefact. Interpreting our results within the framework of the Luttinger liquid theory, we discuss the origin of the apparent contradiction between the predictions of the perturbative Renormalization group approach and numerical calculations at small sizes., Comment: 23 pages, latex, 8 postscript figures

Using a collective-mode Monte Carlo method (the Wolff-Swendsen-Wang algorithm), we compute the spin-stiffness of the two-dimensional classical Heisenberg model. We show that it is the relevant physical quantity to investigate the behaviour of the model in the very low temperature range inaccessible to previous studies based on correlation length and susceptibility calculations., Comment: 6 pages, latex, 3 postscript figures appended, DIM preprint 93-30

We present a powerful method for calculating the thermodynamic properties of the Hubbard model in infinite dimensions, using an exact diagonalization of an Anderson model with a finite number of sites. At finite temperatures, the explicit diagonalization of the Anderson Hamiltonian allows the calculation of Green's functions with a resolution far superior to that of Quantum Monte Carlo calculations. At zero temperature, the Lancz\`os method is used and yields the essentially exact zero-temperature solution of the model, except in a region of very small frequencies. Numerical results for the half-filled case in the paramagnetic phase (quasi-particle weight, self-energy, and also real-frequency spectral densities) are presented., Comment: 7 pages, revtex, 4 postscript figures appended, LPQ preprint (1993)

We apply an exact diagonalization method to the the infinite-D two-band Hubbard model. The method is essentially exact for the calculation of thermodynamic properties for all but the smallest frequencies and yields a resolution unavailable in Monte Carlo calculations. We establish the instability of the normal state with respect to singlet superconductivity at small frequencies at small doping, the regime of relevance for High-$T_c$ superconductors. We also present evidence for the existence of an instability towards triplet superconductivity in the large doping regime $n \sim 2$., Comment: 7 pages, revtex, 3 postscript figures appended, LPSENS preprint 93-17