Source author record

Peter M. W. Gill

Peter M. W. Gill 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

29works
10topics
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

29 published item(s)

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.

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.

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.

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.

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.

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.

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.

preprint2013arXiv

Distribution of $r_{12} \cdot p_{12}$ in quantum systems

We introduce the two-particle probability density $X(x)$ of $x=\bm{r}_{12}\cdot\bm{p}_{12}=\left(\bm{r}_1-\bm{r}_2\right) \cdot \left(\bm{p}_1-\bm{p}_2\right)$. We show how to derive $X(x)$, which we call the Posmom intracule, from the many-particle wavefunction. We contrast it with the Dot intracule [Y. A. Bernard, D. L. Crittenden, P. M. W. Gill, Phys. Chem. Chem. Phys., 10, 3447 (2008)] which can be derived from the Wigner distribution and show the relationships between the Posmom intracule and the one-particle Posmom density [Y. A. Bernard, D. L. Crittenden, P. M. W .Gill, J.Phys. Chem.A, 114, 11984 (2010)]. To illustrate the usefulness of $X(x)$, we construct and discuss it for a number of two-electron systems.

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.

preprint2013arXiv

Resolutions of the Coulomb operator: VI. Computation of auxiliary integrals

We discuss the efficient computation of the auxiliary integrals that arise when resolutions of two-electron operators (specifically, the Coulomb and long-range Ewald operators) are employed in quantum chemical calculations. We derive a recurrence relation that facilitates the generation of auxiliary integrals for Gaussian basis functions of arbitrary angular momentum and propose a near-optimal algorithm for its use.

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.

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.

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.

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.

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.

preprint2011arXiv

A Remarkable Identity Involving Bessel Functions

We consider a new identity involving integrals and sums of Bessel functions. The identity provides new ways to evaluate integrals of products of two Bessel functions. The identity is remarkably simple and powerful since the summand and integrand are of exactly the same form and the sum converges to the integral relatively fast for most cases. A proof and numerical examples of the identity are discussed.

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

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.

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.

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.

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.

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

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

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.

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.

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.

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.