BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michel Caffarel

Michel Caffarel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phphysics.comp-phcond-mat.str-elcond-mat.otherphysics.atom-ph

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Diffusion Monte Carlo using domains in configuration space

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

QUESTDB: a database of highly-accurate excitation energies for the electronic structure community

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 ($π\to π^*$, $n \to π^*$, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Influence of Pseudopotentials on Excitation Energies From Selected Configuration Interaction and Diffusion Monte Carlo

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Quantum Package 2.0: An Open-Source Determinant-Driven Suite of Programs

\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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Self-Consistent Electron-Nucleus Cusp Correction for Molecular Orbitals

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Toward an improved control of the fixed-node error in quantum Monte Carlo: The case of the water molecule

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fixed-Node Diffusion Monte Carlo potential energy curve of the fluorine molecule F2 using selected configuration interaction trial wavefunctions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spin density distribution in open-shell transition metal systems: A comparative post-Hartree-Fock, Density Functional Theory and quantum Monte Carlo study of the CuCl2 molecule

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum Monte Carlo with Jastrow-valence-bond wave functions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum Monte Carlo calculations of electronic excitation energies: the case of the singlet $n \to π^*$ (CO) transition in acrolein

We report state-of-the-art quantum Monte Carlo calculations of the singlet $n \to π^*$ (CO) vertical excitation energy in the acrolein molecule, extending the recent study of Bouabç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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

An efficient sampling algorithm for Variational Monte Carlo

We propose a new algorithm for sampling the $N$-body density $|Ψ({\bf R})|^2/\int_{\mathbb{R}^{3N}} |Ψ|^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).

0 citations0 reviewsNo code linkedOpen access