BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Anthony Scemama

Anthony Scemama contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phphysics.comp-phcond-mat.str-elcond-mat.mtrl-sciquant-phcond-mat.otherphysics.atom-phphysics.class-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

18 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
preprint2022arXiv

Double excitation energies from quantum Monte Carlo using state-specific energy optimization

We show that recently developed quantum Monte Carlo methods, which provide accurate vertical transition energies for single excitations, also successfully treat double excitations. We study the double excitations in medium-sized molecules, some of which are challenging for high level coupled-cluster calculations to model accurately. Our fixed-node diffusion Monte Carlo excitation energies are in very good agreement with reliable benchmarks, when available, and provide accurate predictions for excitation energies of difficult systems where reference values are lacking.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimization of large determinant expansions in quantum Monte Carlo

We present a new method for the optimization of large configuration interaction (CI) expansions in the quantum Monte Carlo (QMC) framework. The central idea here is to replace the non-orthogonal variational optimization of CI coefficients performed in usual QMC calculations by an orthogonal non-Hermitian optimization thanks to the so-called transcorrelated (TC) framework, the two methods yielding the same results in the limit of a complete basis set. By rewriting the TC equations as an effective self-consistent Hermitian problem, our approach requires the sampling of a single quantity per Slater determinant, leading to minimal memory requirements in the QMC code. Using analytical quantities obtained from both the TC framework and the usual CI-type calculations, we also propose improved estimators which reduce the statistical fluctuations of the sampled quantities by more than an order of magnitude. We demonstrate the efficiency of this method on wave functions containing $10^5-10^6$ Slater determinants, using effective core potentials or all-electron calculations. In all the cases, a sub-milliHartree convergence is reached within only two or three iterations of optimization.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The effect of uncertainty on building blocks in molecules

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

What is the number of electrons in a spatial domain?

We like to attribute a number of electrons to spatial domains (atoms, bonds, ...). However, as a rule, the number of electrons in a spatial domain is not a sharp number. We thus study probabilities for having any number of electrons (between 0 and the total number of electrons in the system) in a given spatial domain. We show that by choosing a domain that maximizes a chosen probability (or is close to it), one obtains higher probabilities for chemically relevant regions. The probability to have a given electronic arrangement, - for example, by attributing a number of electrons to an atomic shell - can be low. It remains so even in the "best" case, i.e, if the spatial domain is chosen to maximize the chosen probability. In other words, the number of electrons in a spatial region significantly fluctuates. The freedom of choosing the number of electrons we are interested in shows that a "chemical" question is not always well-posed. We show it using as an example the KrF2 molecule.

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
preprint2021arXiv

Tailoring CIPSI expansions for QMC calculations of electronic excitations: the case study of thiophene

The perturbatively selected configuration interaction scheme (CIPSI) is particularly effective in constructing determinantal expansions for quantum Monte Carlo (QMC) simulations with Jastrow-Slater wave functions: fast and smooth convergence of ground-state properties, as well as balanced descriptions of ground- and excited-states of different symmetries have been reported. In particular, accurate excitation energies have been obtained by the pivotal requirement of using CIPSI expansions with similar second-order perturbation corrections for each state, that is, similar estimated errors with respect to the full configuration interaction limit. Here we elaborate on the CIPSI selection criterion for excited states of the same symmetry as the ground state, generating expansions from a common orbital set. Using these expansions in QMC as determinantal components of Jastrow-Slater wave functions, we compute the lowest, bright excited state of thiophene, which is challenging due to its significant multireference character. The resulting vertical excitation energies are within 0.05~eV of the best theoretical estimates, already with expansions of only a few thousand determinants. Furthermore, we relax the ground- and excited-state structures following the corresponding root in variational Monte Carlo and obtain bond lengths which are accurate to better than 0.01~Å. Therefore, while the full treatment at the CIPSI level of this system would be quite demanding, in QMC we can compute high-quality excitation energies and excited-state structural parameters building on affordable CIPSI expansions with relatively few, well chosen determinants.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Density-Based Basis-Set Incompleteness Correction for GW Methods

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Mountaineering Strategy to Excited States: Highly-Accurate Energies and Benchmarks for Exotic Molecules and Radicals

Aiming at completing the sets of FCI-quality transition energies that we recently developed (\textit{J.~Chem.~Theory Comput.} \textbf{14} (2018) 4360--4379, \textit{ibid.}~\textbf{15} (2019) 1939--1956, and \textit{ibid.}~\textbf{16} (2020) 1711--1741), we provide, in the present contribution, ultra-accurate vertical excitation energies for a series of "exotic" closed-shell molecules containing F, Cl, P, and Si atoms and small radicals, such as CON and its variants, that were not considered to date in such investigations. This represents a total of 81 high-quality transitions obtained with a series of diffuse-containing basis sets of various sizes. For the exotic compounds, these transitions are used to perform benchmarks with a vast array of lower-level models, $\textit{i.e.}$ CIS(D), EOM-MP2, (SOS/SCS)-CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, (SOS-)ADC(2), and ADC(3). Additional comparisons are made with literature data. For the open-shell compounds, we have compared the performances of both the unrestricted and restricted open-shell CCSD and CC3 formalisms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Mountaineering Strategy to Excited States: Highly-Accurate Energies and Benchmarks for Medium Size Molecules

Following our previous work focussing on compounds containing up to 3 non-hydrogen atoms [\emph{J. Chem. Theory Comput.} {\bfseries 14} (2018) 4360--4379], we present here highly-accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms. To obtain these energies, we use equation-of-motion coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ), selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational $n$-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/\emph{aug}-cc-pVQZ vertical excitation energies as well as CC3/\emph{aug}-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of $\pm 0.04$ eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse and accurate set of more than 200 highly-accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, $n \rightarrow π^*$, $π\rightarrow π^*$, \ldots). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, as well as NEVPT2. The results of these benchmarks are compared to the available literature data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pros and Cons of the Bethe-Salpeter Formalism for Ground-State Energies

The combination of the many-body Green's function $GW$ approximation and the Bethe-Salpeter equation (BSE) formalism has shown to be a promising alternative to time-dependent density-functional theory (TD-DFT) for computing vertical transition energies and oscillator strengths in molecular systems. The BSE formalism can also be employed to compute ground-state correlation energies thanks to the adiabatic-connection fluctuation-dissipation theorem (ACFDT). Here, we study the topology of the ground-state potential energy surfaces (PES) of several diatomic molecules near their equilibrium bond length. Thanks to comparisons with state-of-art computational approaches (CC3), we show that ACFDT@BSE is surprisingly accurate, and can even compete with lower-order coupled cluster methods (CC2 and CCSD) in terms of total energies and equilibrium bond distances for the considered systems. However, we sometimes observe unphysical irregularities on the ground-state PES in relation with difficulties in the identification of a few $GW$ quasiparticle energies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Quest For Highly Accurate Excitation Energies: A Computational Perspective

We provide an overview of the successive steps that made possible to obtain increasingly accurate excitation energies with computational chemistry tools, eventually leading to chemically accurate vertical transition energies for small- and medium-size molecules. First, we describe the evolution of \textit{ab initio} methods employed to define benchmark values, with originally Roos' CASPT2 method, then the CC3 method as in the renowned Thiel set, and more recently the resurgence of selected configuration interaction methods. The latter method has been able to deliver consistently, for both single and double excitations, highly accurate excitation energies for small molecules, as well as medium-size molecules with compact basis sets. Second, we describe how these high-level methods and the creation of representative benchmark sets of excitation energies have allowed to assess fairly and accurately the performance of computationally lighter methods. We conclude by discussing the future theoretical and technological developments in the field.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Density-Based Basis-Set Correction For Wave Function Theory

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} $μ$, the key ingredient here is a range-separation \textit{function} $μ(\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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Chemically Accurate Excitation Energies With Small Basis Sets

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-$ζ$ 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.

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