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

Pierre-François Loos

Pierre-François Loos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phcond-mat.str-elphysics.comp-phcond-mat.othercond-mat.mtrl-scinucl-thcond-mat.quant-gasmath-phmath.MPcond-mat.mes-hallphysics.atom-phphysics.pop-phquant-ph

Catalog footprint

What is connected

61works
13topics
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

61 published item(s)

preprint2022arXiv

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events, but also to provide reference values for vertical transition energies, hence allowing to benchmark lower-order models. In this category, both CASPT2 and NEVPT2 are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches the present work relies on highly-accurate ($\pm 0.03$ eV) \emph{aug}-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 $n \to π^*$, 119 $π\to π^*$, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially-contracted (PC) and strongly-contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of $0.11$ and $0.13$ eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focussed on single-reference wave function methods (\textit{J.~Chem. Theory Comput.} \textbf{14}, 4360 (2018); \emph{ibid.}, \textbf{16}, 1711 (2020)), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Benchmarking CASPT3 Vertical Excitation Energies

Based on 280 reference vertical transition energies of various natures (singlet, triplet, valence, Rydberg, $n\toπ^*$, $π\toπ^*$, and double excitations) extracted from the QUEST database, we assess the accuracy of third-order multireference perturbation theory, CASPT3, in the context of molecular excited states. When one applies the disputable ionization-potential-electron-affinity (IPEA) shift, we show that CASPT3 provides a similar accuracy as its second-order counterpart, CASPT2, with the same mean absolute error of $0.11$ eV. However, as already reported, we also observe that the accuracy of CASPT3 is almost insensitive to the IPEA shift, irrespective of the transition type and system size, with a small reduction of the mean absolute error to $0.09$ eV when the IPEA shift is switched off.

0 citations0 reviewsNo code linkedOpen access
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

Hierarchy Configuration Interaction: Combining Seniority Number and Excitation Degree

We propose a novel partitioning of the Hilbert space, hierarchy configuration interaction (hCI), where the excitation degree (with respect to a given reference determinant) and the seniority number (i.e., the number of unpaired electrons) are combined in a single hierarchy parameter. The key appealing feature of hCI is that each hierarchy level accounts for all classes of determinants whose number share the same scaling with system size. By surveying the dissociation of multiple molecular systems, we found that the overall performance of hCI usually exceeds or, at least, parallels that of excitation-based CI. For higher orders of hCI and excitation-based CI, the additional computational burden related to orbital optimization usually do not compensate the marginal improvements compared with results obtained with Hartree-Fock orbitals. The exception is orbital-optimized CI with single excitations, a minimally correlated model displaying the qualitatively correct description of single bond breaking, at a very modest computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Static and Dynamic Bethe-Salpeter Equations in the $T$-Matrix Approximation

While the well-established $GW$ approximation corresponds to a resummation of the direct ring diagrams and is particularly well suited for weakly-correlated systems, the $T$-matrix approximation does sum ladder diagrams up to infinity and is supposedly more appropriate in the presence of strong correlation. Here, we derive and implement, for the first time, the static and dynamic Bethe-Salpeter equations when one considers $T$-matrix quasiparticle energies as well as a $T$-matrix-based kernel. The performance of the static scheme and its perturbative dynamical correction are assessed by computing the neutral excited states of molecular systems. Comparison with more conventional schemes as well as other wave function methods are also reported. Our results suggest that the $T$-matrix-based formalism performs best in few-electron systems where the electron density remains low.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Transient Uniform Electron Gases

The uniform electron gas (UEG), a hypothetical system with finite homogenous electron density composed by an infinite number of electrons in a box of infinite volume, is the practical pillar of density-functional theory (DFT) and the foundation of the most acclaimed approximation of DFT, the local-density approximation (LDA). In the last thirty years, the knowledge of analytical parametrizations of the infinite UEG (IUEG) exchange-correlation energy has allowed researchers to perform a countless number of approximate electronic structure calculations for atoms, molecules, and solids. Recently, it has been shown that the traditional concept of the IUEG is not the unique example of UEGs, and systems, in their lowest-energy state, consisting of electrons that are confined to the surface of a sphere provide a new family of UEGs with more customizable properties. Here, we show that, some of the excited states associated with these systems can be classified as transient UEGs (TUEGs) as their electron density is only homogenous for very specific values of the radius of the sphere even though the electronic wave function is not rotationally invariant. Concrete examples are provided in the case of two-electron systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Unphysical Discontinuities, Intruder States and Regularization in $GW$ Methods

By recasting the non-linear frequency-dependent $GW$ quasiparticle equation into a linear eigenvalue problem, we explain the appearance of multiple solutions and unphysical discontinuities in various physical quantities computed within the $GW$ approximation. Considering the $GW$ self-energy as an effective Hamiltonian, it is shown that these issues are key signatures of strong correlation in the $(N\pm1)$-electron states and can be directly related to the intruder state problem. A simple and efficient regularization procedure inspired by the similarity renormalization group is proposed to avoid such issues and speed up convergence of partially self-consistent $GW$ calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Mountaineering Strategy to Excited States: Highly-Accurate Energies and Benchmarks for Bicyclic Systems

Pursuing our efforts to define highly-accurate estimates of the relative energies of excited states in organic molecules, we investigate, with coupled-cluster methods including iterative triples (CC3 and CCSDT), the vertical excitation energies of 10 bicyclic molecules (azulene, benzoxadiazole, benzothiadiazole, diketopyrrolopyrrole, fuofuran, phthalazine, pyrrolopyrrole, quinoxaline, tetrathiafulvalene, and thienothiophene). In total, we provide \emph{aug}-cc-pVTZ reference vertical excitation energies for 91 excited states of these relatively large systems. We use these reference values to benchmark various wave function methods, i.e., CIS(D), EOM-MP2, CC2, CCSD, STEOM-CCSD, CCSD(T)(a)*, CCSDR(3), CCSDT-3, ADC(2), ADC(2.5), ADC(3), as well as some spin-scaled variants of both CC2 and ADC(2). These results are compared to those obtained previously on smaller molecules. It turns out that while the accuracy of some methods is almost unaffected by system size, e.g., CIS(D) and CC3, the performance of others can significantly deteriorate as the systems grow, e.g., EOM-MP2 and CCSD, whereas others, e.g., ADC(2) and CC2, become more accurate for larger derivatives.

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

A weight-dependent local correlation density-functional approximation for ensembles

We report a local, weight-dependent correlation density-functional approximation that incorporates information about both ground and excited states in the context of density-functional theory for ensembles (eDFT). This density-functional approximation for ensembles is specially designed for the computation of single and double excitations within Gross--Oliveira--Kohn (GOK) DFT (i.e., eDFT for neutral excitations), and can be seen as a natural extension of the ubiquitous local-density approximation in the context of ensembles. The resulting density-functional approximation, based on both finite and infinite uniform electron gas models, automatically incorporates the infamous derivative discontinuity contributions to the excitation energies through its explicit ensemble weight dependence. Its accuracy is illustrated by computing single and double excitations in one-dimensional many-electron systems in the weak, intermediate and strong correlation regimes. Although the present weight-dependent functional has been specifically designed for one-dimensional systems, the methodology proposed here is general, i.e., directly applicable to the construction of weight-dependent functionals for realistic three-dimensional systems, such as molecules and solids.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamical Correction to the Bethe-Salpeter Equation Beyond the Plasmon-Pole Approximation

The Bethe-Salpeter equation (BSE) formalism is a computationally affordable method for the calculation of accurate optical excitation energies in molecular systems. Similar to the ubiquitous adiabatic approximation of time-dependent density-functional theory, the static approximation, which substitutes a dynamical (i.e., frequency-dependent) kernel by its static limit, is usually enforced in most implementations of the BSE formalism. Here, going beyond the static approximation, we compute the dynamical correction of the electron-hole screening for molecular excitation energies thanks to a renormalized first-order perturbative correction to the static BSE excitation energies. The present dynamical correction goes beyond the plasmon-pole approximation as the dynamical screening of the Coulomb interaction is computed exactly within the random-phase approximation. Our calculations are benchmarked against high-level (coupled-cluster) calculations, allowing to assess the clear improvement brought by the dynamical correction for both singlet and triplet optical transitions.

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 Bethe-Salpeter Equation Formalism: From Physics to Chemistry

The Bethe-Salpeter equation (BSE) formalism is steadily asserting itself as a new efficient and accurate tool in the ensemble of computational methods available to chemists in order to predict optical excitations in molecular systems. In particular, the combination of the so-called $GW$ approximation, giving access to reliable ionization energies and electron affinities, and the BSE formalism, able to model UV/Vis spectra, has shown to provide accurate singlet excitation energies with a typical error of $0.1$--$0.3$ eV. With a similar computational cost as time-dependent density-functional theory (TD-DFT), BSE is able to provide an accuracy on par with the most accurate global and range-separated hybrid functionals without the unsettling choice of the exchange-correlation functional, resolving further known issues (\textit{e.g.}, charge-transfer excitations). In this \textit{Perspective} article, we provide a historical overview of BSE, with a particular focus on its condensed-matter roots. We also propose a critical review of its strengths and weaknesses in different chemical situations.

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
preprint2020arXiv

Weight Dependence of Local Exchange-Correlation Functionals in Ensemble Density-Functional Theory: Double Excitations in Two-Electron Systems

Gross--Oliveira--Kohn (GOK) ensemble density-functional theory (GOK-DFT) is a time-\textit{independent} extension of density-functional theory (DFT) which allows to compute excited-state energies via the derivatives of the ensemble energy with respect to the ensemble weights. Contrary to the time-dependent version of DFT (TD-DFT), double excitations can be easily computed within GOK-DFT. However, to take full advantage of this formalism, one must have access to a \textit{weight-dependent} exchange-correlation functional in order to model the infamous ensemble derivative contribution to the excitation energies. In the present article, we discuss the construction of first-rung (i.e., local) weight-dependent exchange-correlation density-functional approximations for two-electron atomic and molecular systems (He and H$_2$) specifically designed for the computation of double excitations within GOK-DFT. In the spirit of optimally-tuned range-separated hybrid functionals, a two-step system-dependent procedure is proposed to obtain accurate energies associated with double excitations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

$\mathcal{PT}$-Symmetry in Hartree-Fock Theory

$\mathcal{PT}$-symmetry --- invariance with respect to combined space reflection $\mathcal{P}$ and time reversal $\mathcal{T}$ --- provides a weaker condition than (Dirac) Hermiticity for ensuring a real energy spectrum of a general non-Hermitian Hamiltonian. $\mathcal{PT}$-symmetric Hamiltonians therefore form an intermediate class between Hermitian and non-Hermitian Hamiltonians. In this work, we derive the conditions for $\mathcal{PT}$-symmetry in the context of electronic structure theory, and specifically, within the Hartree-Fock (HF) approximation. We show that the HF orbitals are symmetric with respect to the $\mathcal{PT}$ operator \textit{if and only if} the effective Fock Hamiltonian is $\mathcal{PT}$-symmetric, and \textit{vice versa}. By extension, if an optimal self-consistent solution is invariant under $\mathcal{PT}$, then its eigenvalues and corresponding HF energy must be real. Moreover, we demonstrate how one can construct explicitly $\mathcal{PT}$-symmetric Slater determinants by forming $\mathcal{PT}$ doublets (i.e. pairing each occupied orbital with its $\mathcal{PT}$-transformed analogue), allowing $\mathcal{PT}$-symmetry to be conserved throughout the self-consistent process. Finally, considering the \ce{H2} molecule as an illustrative example, we observe $\mathcal{PT}$-symmetry in the HF energy landscape and find that the symmetry-broken unrestricted HF wave functions (i.e. diradical configurations) are $\mathcal{PT}$-symmetric, while the symmetry-broken restricted HF wave functions (i.e. ionic configurations) break $\mathcal{PT}$-symmetry.

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

A Wigner molecule at extremely low densities: a numerically exact study

In this work we investigate Wigner localization at very low densities by means of the exact diagonalization of the Hamiltonian. This yields numerically exact results. In particular, we study a quasi-one-dimensional system of two electrons that are confined to a ring by three-dimensional gaussians placed along the ring perimeter. To characterize the Wigner localization we study several appropriate observables, namely the two-body reduced density matrix, the localization tensor and the particle-hole entropy. We show that the localization tensor is the most promising quantity to study Wigner localization since it accurately captures the transition from the delocalized to the localized state and it can be applied to systems of all sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Capturing static and dynamic correlation with $Δ\text{NO}$-MP2 and $Δ\text{NO}$-CCSD

The $Δ\text{NO}$ method for static correlation is combined with second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster singles and doubles (CCSD) to account for dynamic correlation. The MP2 and CCSD expressions are adapted from finite-temperature CCSD, which includes orbital occupancies and vacancies, and expanded orbital summations. Correlation is partitioned with the aid of damping factors incorporated into the MP2 and CCSD residual equations. Potential energy curves for a selection of diatomics are in good agreement with extrapolated full configuration interaction results (exFCI), and on par with conventional multireference approaches.

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

Complex Adiabatic Connection: a Hidden Non-Hermitian Path from Ground to Excited States

Processes related to electronically excited states are central in many areas of science, however accurately determining excited-state energies remains a major challenge in theoretical chemistry. Recently, higher energy stationary states of non-linear methods have themselves been proposed as approximations to excited states, although the general understanding of the nature of these solutions remains surprisingly limited. In this Letter, we present an entirely novel approach for exploring and obtaining excited stationary states by exploiting the properties of non-Hermitian Hamiltonians. Our key idea centres on performing analytic continuations of conventional quantum chemistry methods. Considering Hartree--Fock theory as an example, we analytically continue the electron-electron interaction to expose a hidden connectivity of multiple solutions across the complex plane, revealing a close resemblance between Coulson--Fischer points and non-Hermitian degeneracies. Finally, we demonstrate how a ground-state wave function can be morphed naturally into an excited-state wave function by constructing a well-defined complex adiabatic connection.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Evaluating 0-0 Energies with Theoretical Tools: a Short Review

For a given electronic excited state, the 0-0 energy ($T_0$ or $T_{00}$) is the simplest property allowing straightforward and physically-sound comparisons between theory and (accurate) experiment. However, the computation of 0-0 energies with \emph{ab initio} approaches requires determining both the structure and the vibrational frequencies of the excited state, which limits the quality of the theoretical models that can be considered in practice. This explains why only a rather limited, yet constantly increasing, number of works have been devoted to the determination of this property. In this contribution, we review these efforts with a focus on benchmark studies carried out for both gas phase and solvated compounds. Over the years, not only as the size of the molecules increased, but the refinement of the theoretical tools has followed the same trend. Though the results obtained in these benchmarks significantly depend on both the details of the protocol and the nature of the excited states, one can now roughly estimate, in the case of valence transitions, the overall accuracy of theoretical schemes as follows: $1$ eV for CIS, $0.2$--$0.3$ eV for CIS(D), $0.2$--$0.4$ eV for TD-DFT when one employs hybrid functionals, $0.1$--$0.2$ eV for ADC(2) and CC2, and $0.04$ eV for CC3, the latter approach being the only one delivering chemical accuracy on a near-systematic basis.

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

Is ADC(3) as Accurate as CC3 for Valence and Rydberg Transition Energies?

The search for new \emph{ab initio} models rapidly delivering accurate excited state energies and properties is one of the most active research lines of theoretical chemistry. Along with these methodological developments, the performances of known methods are constantly reassessed thanks to the emergence of new benchmark values. In this Letter, we show that, in contrast to previous claims, the third-order algebraic diagrammatic construction, ADC(3), does not yield transition energies of the same quality as the third-order coupled cluster method, CC3. There is indeed a significant difference in terms of accuracy between the two approaches, as we clearly and unambiguously demonstrate here thanks to extensive comparisons with several hundreds high-quality vertical transition energies obtained with FCI, CCSDTQ, and CCSDT. Direct comparisons with experimental 0-0 energies of small- and medium-size organic molecules support the same conclusion, which holds for both valence and Rydberg transitions, as well as singlet and triplet states. In regards of these results, we introduce a composite approach that we named ADC(2.5) which consists in averaging the ADC(2) and ADC(3) excitation energies. Although ADC(2.5) does not match the CC3 accuracy, it significantly improves the ADC(3) results, especially for vertical energies. We hope that the present contribution will stimulate further developments and, in particular, improvements of the ADC-type methods which have the indisputable advantage of being computationally lighter than their equivalent-order CC variants.

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

Many-electron integrals over gaussian basis functions. I. Recurrence relations for three-electron integrals

Explicitly-correlated F12 methods are becoming the first choice for high-accuracy molecular orbital calculations, and can often achieve chemical accuracy with relatively small gaussian basis sets. In most calculations, the many three- and four-electron integrals that formally appear in the theory are avoided through judicious use of resolutions of the identity (RI). However, in order not to jeopardize the intrinsic accuracy of the F12 wave function, the associated RI auxiliary basis set must be large. Here, inspired by the Head-Gordon-Pople (HGP) and PRISM algorithms for two-electron integrals, we present an algorithm to compute directly three-electron integrals over gaussian basis functions and a very general class of three-electron operators, without invoking RI approximations. A general methodology to derive vertical, transfer and horizontal recurrence relations is also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The uniform electron gas

The uniform electron gas or UEG (also known as jellium) is one of the most fundamental models in condensed-matter physics and the cornerstone of the most popular approximation --- the local-density approximation --- within density-functional theory. In this article, we provide a detailed review on the energetics of the UEG at high, intermediate and low densities, and in one, two and three dimensions. We also report the best quantum Monte Carlo and symmetry-broken Hartree-Fock calculations available in the literature for the UEG and discuss the phase diagrams of jellium.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nodal surfaces and interdimensional degeneracies

The aim of this paper is to shed light on the topology and properties of the nodes (i.e. the zeros of the wave function) in electronic systems. Using the "electrons on a sphere" model, we study the nodes of two-, three- and four-electron systems in various ferromagnetic configurations ($sp$, $p^2$, $sd$, $pd$, $p^3$, $sp^2$ and $sp^3$). In some particular cases ($sp$, $p^2$, $sd$, $pd$ and $p^3$), we rigorously prove that the non-interacting wave function has the same nodes as the exact (yet unknown) wave function. The number of atomic and molecular systems for which the exact nodes are known analytically is very limited and we show here that this peculiar feature can be attributed to interdimensional degeneracies. Although we have not been able to prove it rigorously, we conjecture that the nodes of the non-interacting wave function for the $sp^3$ configuration are exact.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Three-electron coalescence points in two and three dimensions

The form of the wave function at three-electron coalescence points is examined for several spin states using an alternative method to the usual Fock expansion. We find that, in two- and three-dimensional systems, the non-analytical nature of the wave function is characterized by the appearance of logarithmic terms, reminiscent of those that appear as both electrons approach the nucleus of the helium atom. The explicit form of these singularities is given in terms of the interelectronic distances for a doublet and two quartet states of three electrons in a harmonic well.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Uniform electron gases: III. Low-density gases on three-dimensional spheres

By combining variational Monte Carlo (VMC) and complete-basis-set limit Hartree-Fock (HF) calculations, we have obtained near-exact correlation energies for low-density same-spin electrons on a three-dimensional sphere (3-sphere), i.e.~the surface of a four-dimensional ball. In the VMC calculations, we compare the efficacies of two types of one-electron basis functions for these strongly correlated systems, and analyze the energy convergence with respect to the quality of the Jastrow factor. The HF calculations employ spherical Gaussian functions (SGFs) which are the curved-space analogs of cartesian Gaussian functions. At low densities, the electrons become relatively localized into Wigner crystals, and the natural SGF centers are found by solving the Thomson problem (i.e. the minimum-energy arrangement of $n$ point charges) on the 3-sphere for various values of $n$. We have found 11 special values of $n$ whose Thomson sites are equivalent. Three of these are the vertices of four-dimensional Platonic solids --- the hyper-tetrahedron ($n=5$), the hyper-octahedron ($n=8$) and the 24-cell ($n=24$) --- and a fourth is a highly symmetric structure ($n=13$) which has not previously been reported. By calculating the harmonic frequencies of the electrons around their equilibrium positions, we also find the first-order vibrational corrections to the Thomson energy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Basis functions for electronic structure calculations on spheres

We introduce a new basis function (the spherical gaussian) for electronic structure calculations on spheres of any dimension $D$. We find \alert{general} expressions for the one- and two-electron integrals and propose an efficient computational algorithm incorporating the Cauchy-Schwarz bound. Using numerical calculations for the $D = 2$ case, we show that spherical gaussians are more efficient than spherical harmonics when the electrons are strongly localized.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Chemistry in One Dimension

We report benchmark results for one-dimensional (1D) atomic and molecular systems interacting via the Coulomb operator $|x|^{-1}$. Using various wavefunction-type approaches, such as Hartree-Fock theory, second- and third-order Møller-Plesset perturbation theory and explicitly correlated calculations, we study the ground state of atoms with up to ten electrons as well as small diatomic and triatomic molecules containing up to two electrons. A detailed analysis of the 1D helium-like ions is given and the expression of the high-density correlation energy is reported. We report the total energies, ionization energies, electron affinities and other interesting properties of the many-electron 1D atoms and, based on these results, we construct the 1D analog of Mendeleev's periodic table. We find that the 1D periodic table contains only two groups: the alkali metals and the noble gases. We also calculate the dissociation curves of various 1D diatomics and study the chemical bond in H$_2^+$, HeH$^{2+}$, He$_2^{3+}$, H$_2$, HeH$^+$ and He$_2^{2+}$. We find that, unlike their 3D counterparts, 1D molecules are primarily bound by one-electron bonds. Finally, we study the chemistry of H$_3^+$ and we discuss the stability of the 1D polymer resulting from an infinite chain of hydrogen atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalized local-density approximation and one-dimensional finite uniform electron gases

We explicitly build a generalized local-density approximation (GLDA) correlation functional based on one-dimensional (1D) uniform electron gases (UEGs). The fundamental parameters of the GLDA \textemdash a generalization of the widely-known local-density approximation (LDA) used in density-functional theory (DFT) \textemdash are the electronic density $ρ$ and a newly-defined two-electron local parameter called the hole curvature $η$. The UEGs considered in this study are finite versions of the conventional infinite homogeneous electron gas and consist of $n$ electrons on a infinitely thin wire with periodic boundary conditions. We perform a comprehensive study of these finite UEGs at high, intermediate and low densities using perturbation theory and quantum Monte Carlo calculations. We show that the present GLDA functional yields accurate estimates of the correlation energy for both weakly and strongly correlated one-dimensional systems and can be easily generalized to higher-dimensional systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Uniform Electron Gases. II. The Generalized Local Density Approximation in One Dimension

We introduce a generalization (gLDA) of the traditional Local Density Approximation (LDA) within density functional theory. The gLDA uses both the one-electron Seitz radius $\rs$ and a two-electron hole curvature parameter $η$ at each point in space. The gLDA reduces to the LDA when applied to the infinite homogeneous electron gas but, unlike the LDA, is is also exact for finite uniform electron gases on spheres. We present an explicit gLDA functional for the correlation energy of electrons that are confined to a one-dimensional space and compare its accuracy with LDA, second- and third-order Møller-Plesset perturbation energies and exact calculations for a variety of inhomogeneous systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exact wave functions for concentric two-electron systems

We show that the exact solution of the Schrödinger equation for two electrons confined to two distinct concentric rings or spheres can be found in closed form for particular sets of the ring or sphere radii. In the case of two concentric rings, we report exact polynomial and irrational solutions. The same methodology is applied to the case of two concentric spheres for which we report exact polynomial solutions for the ground state and the excited states of $S$ symmetry. For these concentric systems, we show that the exact wave function does not contain terms proportional to the interelectronic distance due to the spatial separation of the electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High-density correlation energy expansion of the one-dimensional uniform electron gas

We show that the expression of the high-density (i.e small-$r_s$) correlation energy per electron for the one-dimensional uniform electron gas can be obtained by conventional perturbation theory and is of the form $\Ec(r_s) = -π^2/360 + 0.00845 r_s + ...$, where $r_s$ is the average radius of an electron. Combining these new results with the low-density correlation energy expansion, we propose a local-density approximation correlation functional, which deviates by a maximum of 0.1 millihartree compared to the benchmark DMC calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Uniform electron gases. I. Electrons on a ring

We introduce a new paradigm for one-dimensional uniform electron gases (UEGs). In this model, $n$ electrons are confined to a ring and interact via a bare Coulomb operator. We use Rayleigh-Schrödinger perturbation theory to show that, in the high-density regime, the ground-state reduced (i.e. per electron) energy can be expanded as $\eps(r_s,n) = \eps_0(n) r_s^{-2} + \eps_1(n) r_s^{-1} + \eps_2(n) +\eps_3(n) r_s + \ldots$, where $r_s$ is the Seitz radius. We use strong-coupling perturbation theory and show that, in the low-density regime, the reduced energy can be expanded as $\eps(r_s,n) = η_0(n) r_s^{-1} + η_1(n) r_s^{-3/2} + η_2(n) r_s^{-2} + \ldots$. We report explicit expressions for $\eps_0(n)$, $\eps_1(n)$, $\eps_2(n)$, $\eps_3(n)$, $η_0(n)$ and $η_1(n)$ and derive the thermodynamic (large-$n$) limits of each of these. Finally, we perform numerical studies of UEGs with $n = 2, 3, \ldots, 10$, using Hylleraas-type and quantum Monte Carlo methods, and combine these with the perturbative results to obtain a picture of the behavior of the new model over the full range of $n$ and $r_s$ values.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Correlation energy of the one-dimensional Coulomb gas

We introduce a new paradigm for finite and infinite strict-one-dimensional uniform electron gases. In this model, $n$ electrons are confined to a ring and interact via a bare Coulomb operator. In the high-density limit (small-$r_s$, where $r_s$ is the Seitz radius), we find that the reduced correlation energy is $\Ec(r_s,n) = \eps^{(2)}(n) + O(r_s)$, and we report explicit expressions for $\eps^{(2)}(n)$. In the thermodynamic (large-$n$) limit of this, we show that $\Ec(r_s) = - π^2/360 + O(r_s)$. In the low-density (large-$r_s$) limit, the system forms a Wigner crystal and we find that $\Ec(r_s) = -[\ln(\sqrt{2π})-3/4] r_s^{-1} + 0.359933 r_s^{-3/2} + O(r_s^{-2})$. Using these results, we propose a correlation functional that interpolates between the high- and low-density limits. The accuracy of the functional for intermediate densities is established by comparison with diffusion Monte Carlo results. Application to a non-uniform system is also reported.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exact wave functions of two-electron quantum rings

We demonstrate that the Schrödinger equation for two electrons on a ring, which is the usual paradigm to model quantum rings, is solvable in closed form for particular values of the radius. We show that both polynomial and irrational solutions can be found for any value of the angular momentum and that the singlet and triplet manifolds, which are degenerate, have distinct geometric phases. We also study the nodal structure associated with these two-electron states.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Harmonically trapped jellium

We discuss the model of a $D$-dimensional confined electron gas in which the particles are trapped by a harmonic potential. In particular, we study the non-interacting kinetic and exchange energies of finite-size inhomogeneous systems, and compare the resulting Thomas-Fermi and Dirac coefficients with various uniform electron gas paradigms. We show that, in the thermodynamic limit, the properties of this model are identical to those of the $D$-dimensional Fermi gas.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Understanding excitons using spherical geometry

Using the spherical geometry, we introduce a novel model to study excitons confined in a three-dimensional space, which offers unparalleled mathematical simplicity while retaining much of the key physics. This new model consists of an exciton trapped on the 3-sphere (i.e. the surface of a four-dimensional ball), and provides a unified treatment of Frenkel and Wannier-Mott excitons. Moreover, we show that one can determine, for particular values of the dielectric constant $ε$, the closed-form expression of the exact wave function. We use the exact wave function of the lowest bound state for $ε=2$ to introduce an intermediate regime which gives satisfactory agreement with \alert{the} exact results for a wide range of $ε$ values.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Uniform electron gases

We show that the traditional concept of the uniform electron gas (UEG) --- a homogeneous system of finite density, consisting of an infinite number of electrons in an infinite volume --- is inadequate to model the UEGs that arise in finite systems. We argue that, in general, a UEG is characterized by at least two parameters, \textit{viz.} the usual one-electron density parameter $ρ$ and a new two-electron parameter $η$. We outline a systematic strategy to determine a new density functional $E(ρ,η)$ across the spectrum of possible $ρ$ and $η$ values.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Correlation energy of anisotropic quantum dots

We study the $D$-dimensional high-density correlation energy $\Ec$ of the singlet ground state of two electrons confined by a harmonic potential with Coulombic repulsion. We allow the harmonic potential to be anisotropic, and examine the behavior of $\Ec$ as a function of the anisotropy $α^{-1}$. In particular, we are interested in the limit where the anisotropy goes to infinity ($α\to0$) and the electrons are restricted to a lower-dimensional space. We show that tuning the value of $α$ from 0 to 1 allows a smooth dimensional interpolation and we demonstrate that the usual model, in which a quantum dot is treated as a two-dimensional system, is inappropriate. Finally, we provide a simple function which reproduces the behavior of $\Ec$ over the entire range of $α$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Correlation energy of the spin-polarized uniform electron gas at high density

The correlation energy per electron in the high-density uniform electron gas can be written as $\Ec(r_s,ζ) = \lam_0(ζ) \ln r_s + \eps_0(ζ) + \lam_1(ζ) \,r_s \ln r_s + O(r_s)$, where $r_s$ is the Seitz radius and $ζ$ is the relative spin polarization. We derive an expression for $\lam_1(ζ)$ which is exact for any $ζ$, including the paramagnetic and ferromagnetic limits, $\lam_1(0)$ and $\lam_1(1)$, and discover that the previously published $\lam_1(1)$ value is incorrect. We trace this error to an integration and limit that do not commute. The spin-resolution of $\lam_1(ζ)$ into contributions of electron pairs is also derived.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Leading-order behavior of the correlation energy in the uniform electron gas

We show that, in the high-density limit, restricted Møller-Plesset (RMP) perturbation theory yields $E_{\text{RMP}}^{(2)} = π^{-2}(1-\ln 2) \ln r_s + O(r_s^0)$ for the correlation energy per electron in the uniform electron gas, where $r_s$ is the Seitz radius. This contradicts an earlier derivation which yielded $E_{\text{RMP}}^{(2)} = O(\ln|\ln r_s|)$. The reason for the discrepancy is explained.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thinking outside the box: the uniform electron gas on a hypersphere

We discuss alternative homogeneous electron gas systems in which a finite number $n$ of electrons are confined to a $D$-dimensional sphere. We derive the first few terms of the high-density ($r_s\to0$, where $r_s$ is the Seitz radius) energy expansions for these systems and show that, in the thermodynamic limit ($n\to\infty$), these terms become identical to those of $D$-dimensional jellium.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Tale of Two Electrons: Correlation at High Density

We review our recent progress in the determination of the high-density correlation energy $\Ec$ in two-electron systems. Several two-electron systems are considered, such as the well known helium-like ions (helium), and the Hooke's law atom (hookium). We also present results regarding two electrons on the surface of a sphere (spherium), and two electrons trapped in a spherical box (ballium). We also show that, in the large-dimension limit, the high-density correlation energy of two opposite-spin electrons interacting {\em via} a Coulomb potential is given by $\Ec \sim -1/(8D^2)$ for any radial external potential $V(r)$, where $D$ is the dimensionality of the space. This result explains the similarity of $\Ec$ in the previous two-electron systems for D=3.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Correlation energy of two electrons in a ball

We study the ground-state correlation energy $E_{\rm c}$ of two electrons of opposite spin confined within a $D$-dimensional ball ($D \ge 2$) of radius $R$. In the high-density regime, we report accurate results for the exact and restricted Hartree-Fock energy, using a Hylleraas-type expansion for the former and a simple polynomial basis set for the latter. By investigating the exact limiting correlation energy $E_{\rm c}^{(0)} = \lim_{R \to 0} \Ec$ for various values of $D$, we test our recent conjecture [J. Chem. Phys. {\bf 131} (2009) 241101] that, in the large-$D$ limit, $E_{\rm c}^{(0)} \sim -δ^2/8$ for any spherically-symmetric confining external potential, where $δ=1/(D-1)$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Correlation energy of two electrons in the high-density limit

We consider the high-density-limit correlation energy $\Ec$ in $D \ge 2$ dimensions for the $^1S$ ground states of three two-electron systems: helium (in which the electrons move in a Coulombic field), spherium (in which they move on the surface of a sphere), and hookium (in which they move in a quadratic potential). We find that the $\Ec$ values are strikingly similar, depending strongly on $D$ but only weakly on the external potential. We conjecture that, for large $D$, the limiting correlation energy $\Ec \sim -δ^2/8$ in any confining external potential, where $δ= 1/(D-1)$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ground state of two electrons on a sphere

We have performed a comprehensive study of the singlet ground state of two electrons on the surface of a sphere of radius $R$. We have used electronic structure models ranging from restricted and unrestricted Hartree-Fock theory to explicitly correlated treatments, the last of which lead to near-exact wavefunctions and energies for any value of $R$. Møller-Plesset energy corrections (up to fifth-order) are also considered, as well as the asymptotic solution in the large-$R$ regime.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ground state of two electrons on concentric spheres

We extend our analysis of two electrons on a sphere [Phys. Rev. A {\bf 79}, 062517 (2009); Phys. Rev. Lett. {\bf 103}, 123008 (2009)] to electrons on concentric spheres with different radii. The strengths and weaknesses of several electronic structure models are analyzed, ranging from the mean-field approximation (restricted and unrestricted Hartree-Fock solutions) to configuration interaction expansion, leading to near-exact wave functions and energies. The Møller-Plesset energy corrections (up to third-order) and the asymptotic expansion for the large-spheres regime are also considered. We also study the position intracules derived from approximate and exact wave functions. We find evidence for the existence of a long-range Coulomb hole in the large-spheres regime, and infer that unrestricted Hartree-Fock theory over-localizes the electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hooke's law correlation in two-electron systems

We study the properties of the Hooke's law correlation energy ($\Ec$), defined as the correlation energy when two electrons interact {\em via} a harmonic potential in a $D$-dimensional space. More precisely, we investigate the $^1S$ ground state properties of two model systems: the Moshinsky atom (in which the electrons move in a quadratic potential) and the spherium model (in which they move on the surface of a sphere). A comparison with their Coulombic counterparts is made, which highlights the main differences of the $\Ec$ in both the weakly and strongly correlated limits. Moreover, we show that the Schrödinger equation of the spherium model is exactly solvable for two values of the dimension ($D = 1 \text{and} 3$), and that the exact wave function is based on Mathieu functions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Invariance of the correlation energy at high density and large dimension in two-electron systems

We prove that, in the large-dimension limit, the high-density correlation energy $\Ec$ of two opposite-spin electrons confined in a $D$-dimensional space and interacting {\em via} a Coulomb potential is given by $\Ec \sim -1/(8D^2)$ for any radial confining potential $V(r)$. This result explains the observed similarity of $\Ec$ in a variety of two-electron systems in three-dimensional space.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Two electrons on a hypersphere: a quasi-exactly solvable model

We show that the exact wave function for two electrons, interacting through a Coulomb potential but constrained to remain on the surface of a $\mathcal{D}$-sphere ($\mathcal{D} \ge 1$), is a polynomial in the interelectronic distance $u$ for a countably infinite set of values of the radius $R$. A selection of these radii, and the associated energies, are reported for ground and excited states on the singlet and triplet manifolds. We conclude that the $\mathcal{D}=3$ model bears the greatest similarity to normal physical systems.

0 citations0 reviewsNo code linkedOpen access