Source author record

Emmanuel Giner

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

17 published item(s)

preprint2022arXiv

Basis-set correction based on density-functional theory: Rigorous framework for a one-dimensional model

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

preprint2022arXiv

Basis-set correction for coupled-cluster estimation of dipole moments

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

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.

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.

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.

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.

preprint2020arXiv

Relativistic short-range exchange energy functionals beyond the local-density approximation

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

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.

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.

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.

preprint2016arXiv

Quantum Monte Carlo with reoptimized perturbatively selected configuration-interaction wave functions

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

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.

preprint2015arXiv

A useful underestimate for the convergence of integral functionals

This article deals with the lower compactness property of a sequence of integrands and the use of this key notion in various domains: convergence theory, optimal control, non-smooth analysis. First about the interchange of the weak epi-limit and the symbol of integration for a sequence of integral functionals. These functionals are defined on a topological space $(\mathcal{X, T})$ where $\mathcal{X}$ is a subset of measurable functions and the $\mathcal{T}$ convergence is stronger than or equal to the convergence in the Bitting sense. Given a sequence $(f_{n})_n$ of integrands, if the integrand $f$ is the weak lower sequential epi-limit of the integrands $f_n$ one of the main results of this article asserts that under the Ioffe's criterion, the $\mathcal{T}$-lower sequential epi-limit of the sequence of integral functionals at the point $x$ is bounded below by the value of the integral functional associated to the Fenchel-Moreau-Rockafellar biconjugate of $f$ at the point $x$. Then the strong-weak semicontinuity (respectively the subdifferentiability) are studied in relation with the Ioffe's criterion. This permits with original proofs to give new conditions for the strong-weak lower semicontinuity at a given point, and to obtain necessary and sufficient conditions for the Fréchet and the (weak)Hadamard subdifferentiability of integral functionals on general spaces, particularly on Lebesgue spaces.

preprint2015arXiv

Efficiency of a Multi-Reference Coupled Cluster method

The multi-reference Coupled Cluster method first proposed by Meller et al (J. Chem. Phys. 1996) has been implemented and tested. Guess values of the amplitudes of the single and double excitations (the ${\hat T}$ operator) on the top of the references are extracted from the knowledge of the coefficients of the Multi Reference Singles and Doubles Configuration Interaction (MRSDCI) matrix. The multiple parentage problem is solved by scaling these amplitudes on the interaction between the references and the Singles and Doubles. Then one proceeds to a dressing of the MRSDCI matrix under the effect of the Triples and Quadruples, the coefficients of which are estimated from the action of ${\hat T}^2$. This dressing follows the logics of the intermediate effective Hamiltonian formalism. The dressed MRSDCI matrix is diagonalized and the process is iterated to convergence. The method is tested on a series of benchmark systems from Complete Active Spaces (CAS) involving 2 or 4 active electrons up to bond breakings. The comparison with Full Configuration Interaction (FCI) results shows that the errors are of the order of a few milli-hartree, five times smaller than those of the CASSDCI. The method is totally uncontracted, parallelizable, and extremely flexible since it may be applied to selected MR and/or selected SDCI. Some potential generalizations are briefly discussed.

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.

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.

preprint2013arXiv

An efficient implementation of Slater-Condon rules

Slater-Condon rules are at the heart of any quantum chemistry method as they allow to simplify $3N$-dimensional integrals as sums of 3- or 6-dimensional integrals. In this paper, we propose an efficient implementation of those rules in order to identify very rapidly which integrals are involved in a matrix element expressed in the determinant basis set. This implementation takes advantage of the bit manipulation instructions on x86 architectures that were introduced in 2008 with the SSE4.2 instruction set. Finding which spin-orbitals are involved in the calculation of a matrix element doesn't depend on the number of electrons of the system.