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

James J. Shepherd

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

ResearcherUnclaimed source record
physics.chem-phphysics.comp-phcond-mat.mtrl-scicond-mat.quant-gascond-mat.str-elcond-mat.othercond-mat.softcond-mat.stat-mechSoftware Engineering

Catalog footprint

What is connected

22works
9topics
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

22 published item(s)

preprint2022arXiv

Machine learning for a finite size correction in periodic coupled cluster theory calculations

We introduce a straightforward Gaussian process regression (GPR) model for the transition structure factor of metal periodic coupled cluster singles and doubles (CCSD) calculations. This is inspired by the method introduced by Liao and Grüneis for interpolating over the transition structure factor to obtain a finite size correction for CCSD [J. Chem. Phys. 145, 141102 (2016)], and by our own prior work using the transition structure factor to efficiently converge CCSD for metals to the thermodynamic limit [Nat. Comput. Sci. 1, 801 (2021)]. In our CCSD-FS-GPR method to correct for finite size errors, we fit the structure factor to a 1D function in the momentum transfer, $G$. We then integrate over this function by projecting it onto a k-point mesh to obtain comparisons with extrapolated results. Results are shown for lithium, sodium, and the uniform electron gas.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Piecewise Interaction Picture Density Matrix Quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact $N$-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC method (IP-DMQMC) which overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Effective Hamiltonians for the study of real metals using quantum chemical theories

Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are at the heart of computational materials science. In that respect the coupled cluster hierarchy of methods plays a central role in molecular quantum chemistry because of its systematic improvability and computational efficiency. In this hierarchy, coupled cluster singles and doubles (CCSD) is one of the most important steps in moving towards chemical accuracy and, in recent years, its scope has successfully been expanded to the study of insulating surfaces and solids. Here, we show that CCSD theory can also be applied to real metals. In so doing, we overcome the limitation of needing extremely large supercells to capture long range electronic correlation effects. An effective Hamiltonian can be found using the transition structure factor--a map of electronic excitations from the Hartree--Fock wavefunction--which has fewer finite size effects than conventional periodic boundary conditions. This not only paves the way of applying coupled cluster methods to real metals but also reduces the computational cost by two orders of magnitude compared to previous methods. Our applications to phases of lithium and silicon show a resounding success in reaching the thermodynamic limit, taking the first step towards a truly universal quantum chemical treatment of solids.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Power laws used to extrapolate the coupled cluster correlation energy to the thermodynamic limit

Recent calculations using coupled cluster on solids have raised discussion of using a $N^{-1/3}$ power law to fit the correlation energy when extrapolating to the thermodynamic limit, an approach which differs from the more commonly used $N^{-1}$ power law which is (for example) often used by quantum Monte Carlo methods. In this paper, we present one way to reconcile these viewpoints. Coupled cluster doubles calculations were performed on uniform electron gases reaching system sizes of $922$ electrons for an extremely wide range of densities ($0.1<r_s<100.0$) to study how the correlation energy approaches the thermodynamic limit. The data were corrected for basis set incompleteness error and use a selected twist angle approach to mitigate finite size error from shell filling effects. Analyzing these data, we initially find that a power law of $N^{-1/3}$ appears to fit the data better than a $N^{-1}$ power law in the large system size limit. However, we provide an analysis of the transition structure factor showing that $N^{-1}$ still applies to large system sizes and that the apparent $N^{-1/3}$ power law occurs only at low $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerating convergence to the thermodynamic limit with twist angle selection applied to methods beyond many-body perturbation theory

We recently developed a scheme to use low-cost calculations to find a single twist angle where the couple cluster doubles energy of a single calculation matches the twist-averaged coupled cluster doubles energy in a finite unit cell. We used initiator full configuration interaction quantum Monte Carlo ($i$-FCIQMC) as an example of an exact method beyond coupled cluster doubles theory to show that this selected twist angle approach had comparable accuracy in methods beyond coupled cluster. Further, at least for small system sizes, we show that the same twist angle can also be found by comparing the energy directly (at the level of second-order Moller-Plesset theory) suggesting a route toward twist angle selection which requires minimal modification to existing codes which can perform twist averaging.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

NECI: N-Electron Configuration Interaction with emphasis on state-of-the-art stochastic methods

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24000 CPU cores. In this paper, we describe the core functionalities of NECI and recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs and it is licensed under GPL-3.0.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Accurate exchange-correlation energies for the warm dense electron gas

Density matrix quantum Monte Carlo (DMQMC) is used to sample exact-on-average $N$-body density matrices for uniform electron gas systems of up to 10$^{124}$ matrix elements via a stochastic solution of the Bloch equation. The results of these calculations resolve a current debate over the accuracy of the data used to parametrize finite-temperature density functionals. Exchange-correlation energies calculated using the real-space restricted path-integral formalism and the $k$-space configuration path-integral formalism disagree by up to $\sim$$10$\% at certain reduced temperatures $T/T_F \le 0.5$ and densities $r_s \le 1$. Our calculations confirm the accuracy of the configuration path-integral Monte Carlo results available at high density and bridge the gap to lower densities, providing trustworthy data in the regime typical of planetary interiors and solids subject to laser irradiation. We demonstrate that DMQMC can calculate free energies directly and present exact free energies for $T/T_F \ge 1$ and $r_s \le 2$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence of many-body wavefunction expansions using a plane wave basis in the thermodynamic limit

Basis set incompleteness error and finite size error can manifest concurrently in systems for which the two effects are phenomenologically well-separated in length scale. When this is true, we need not necessarily remove the two sources of error simultaneously. Instead, the errors can be found and remedied in different parts of the basis set. This would be of great benefit to a method such as coupled cluster theory since the combined cost of $n_{occ}^6 n_{virt}^4$ could be separated into $n_{occ}^6$ and $n_{virt}^4$ costs with smaller prefactors. In this Communication, we present analysis on a data set due to Baardsen and coworkers, containing coupled cluster doubles energies for the 2DEG for $r_s=$ 0.5, 1.0 and 2.0 a.u.~at a wide range of basis set sizes and particle numbers. In obtaining complete basis set limit thermodynamic limit results, we find that within a small and removable error the above assertion is correct for this simple system. This approach allows for the combination of methods which separately address finite size effects and basis set incompleteness error.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Using full configuration interaction quantum Monte Carlo in a seniority zero space to investigate the correlation energy equivalence of pair coupled cluster doubles and doubly occupied configuration interaction

Over the past few years pair coupled cluster doubles (pCCD) has shown promise for the description of strong correlation. This promise is related to its apparent ability to match results from doubly occupied configuration interaction (DOCI), even though the latter method has exponential computational cost. Here, by modifying the full configuration interaction quantum Monte Carlo (FCIQMC) algorithm to sample only the seniority zero sector of Hilbert space, we show that the DOCI and pCCD energies are in agreement for a variety of 2D Hubbard models, including for systems well out of reach for conventional configuration interaction algorithms. Our calculations are aided by the sign problem being much reduced in the seniority zero space compared with the full space. We present evidence for this, and then discuss the sign problem in terms of the wave function of the system which appears to have a simplified sign structure.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interaction Picture Density Matrix Quantum Monte Carlo

The recently developed density matrix quantum Monte Carlo (DMQMC) algorithm stochastically samples the N -body thermal density matrix and hence provides access to exact properties of many-particle quantum systems at arbitrary temperatures. We demonstrate that moving to the interaction picture provides substantial benefits when applying DMQMC to interacting fermions. In this first study, we focus on a system of much recent interest: the uniform electron gas in the warm dense regime. The basis set incompleteness error at finite temperature is investigated and extrapolated via a simple Monte Carlo sampling procedure. Finally, we provide benchmark calculations for a four-electron system, comparing our results to previous work where possible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Open-source development experiences in scientific software: the HANDE quantum Monte Carlo project

The HANDE quantum Monte Carlo project offers accessible stochastic algorithms for general use for scientists in the field of quantum chemistry. HANDE is an ambitious and general high-performance code developed by a geographically-dispersed team with a variety of backgrounds in computational science. In the course of preparing a public, open-source release, we have taken this opportunity to step back and look at what we have done and what we hope to do in the future. We pay particular attention to development processes, the approach taken to train students joining the project, and how a flat hierarchical structure aids communication

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coupled Cluster Channels in the Homogeneous Electron Gas

We discuss diagrammatic modifications to the coupled cluster doubles (CCD) equations, wherein different groups of terms out of rings, ladders, crossed-rings and mosaics can be removed to form approximations to the coupled cluster method, of interest due to their similarity with various types of random phase approximations. The finite uniform electron gas is benchmarked for 14- and 54-electron systems at the complete basis set limit over a wide density range and performance of different flavours of CCD are determined. These results confirm that rings generally overcorrelate and ladders generally undercorrelate; mosaics-only CCD yields a result surprisingly close to CCD. We use a recently developed numerical analysis [J. J. Shepherd and A. Grüneis, Phys. Rev. Lett. 110, 226401 (2013)] to study the behaviours of these methods in the thermodynamic limit. We determine that the mosaics, on forming the Brueckner Hamltonian, open a gap in the effective one-particle eigenvalues at the Fermi energy. Numerical evidence is presented which shows that methods based on this renormalisation have convergent energies in the thermodynamic limit including mosaic-only CCD, which is just a renormalised MP2. All other methods including only a single channel, namely ladder-only CCD, ring-only CCD and crossed-ring-only CCD, appear to yield divergent energies; incorporation of mosaic terms prevents this from happening.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Range Separated Brueckner Coupled Cluster Doubles Theory

We introduce a range-separation approximation to coupled cluster doubles (CCD) theory that successfully overcomes limitations of regular CCD when applied to the uniform electron gas. We combine the short-range ladder channel with the long-range ring channel in the presence of a Bruckner renormalized one-body interaction and obtain ground-state energies with an accuracy of 0.001 a.u./electron across a wide range of density regimes. Our scheme is particularly useful in the low-density and strongly-correlated regimes, where regular CCD has serious drawbacks. Moreover, we cure the infamous overcorrelation of approaches based on ring diagrams (i.e. the particle-hole random phase approximation). Our energies are further shown to have appropriate basis set and thermodynamic limit convergence, and overall this scheme promises energetic properties for realistic periodic and extended systems which existing methods do not possess.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The sign problem in full configuration interaction quantum Monte Carlo: Linear and sub-linear representation regimes for the exact wave function

We investigate the sign problem for full configuration interaction quantum Monte Carlo (FCIQMC), a stochastic algorithm for finding the ground state solution of the Schrödinger equation with substantially reduced computational cost compared with exact diagonalisation. We find $k$-space Hubbard models for which the solution is yielded with storage that grows sub-linearly in the size of the many-body Hilbert space, in spite of using a wave function that is simply linear combination of states. The FCIQMC algorithm is able to find this sub-linear scaling regime without bias and with only a choice of Hamiltonian basis. By means of a demonstration we solve for the energy of a 70-site half-filled system (with a space of $10^{38}$ determinants) in 250 core hours, substantially quicker than the $\sim$10$^{36}$ core hours that would be required by exact diagonalisation. This is the largest space that has been sampled in an unbiased fashion. The challenge for the recently-developed FCIQMC method is made clear: expand the sub-linear scaling regime whilst retaining exact on average accuracy. This result rationalizes the success of the initiator adaptation (i-FCIQMC) and offers clues to improve it. We argue that our results changes the landscape for development of FCIQMC and related methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy landscapes, structural topologies and rearrangement mechanisms in clusters of dipolar particles

Clusters of spherical particles with isotropic attraction favour compact structures that maximise the number of energetically optimal nearest-neighbour interactions. In contrast, dipolar interactions lead to the formation of chains with a low coordination number. When both isotropic and dipolar interactions are present, the competition between them can lead to intricate knot, link and coil structures. Here, we investigate how these structures may self-organise and interconvert in clusters bound by the Stockmayer potential (Lennard-Jones plus point dipole). We map out the low-lying region of the energy landscape using disconnectivity graphs to follow how it evolves as the strength of the dipolar interactions increases. From comprehensive surveys of isomerisation pathways, we identify a number of rearrangement mechanisms that allow the topology of chain-like structures to interconvert.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Explicitly correlated plane waves: Accelerating convergence in periodic wavefunction expansions

We present an investigation into the use of an explicitly correlated plane wave basis for periodic wavefunction expansions at the level of second-order Møller-Plesset perturbation theory (MP2). The convergence of the electronic correlation energy with respect to the one-electron basis set is investigated and compared to conventional MP2 theory in a finite homogeneous electron gas model. In addition to the widely used Slater-type geminal correlation factor, we also derive and investigate a novel correlation factor that we term Yukawa-Coulomb. The Yukawa-Coulomb correlation factor is motivated by analytic results for two electrons in a box and allows for a further improved convergence of the correlation energies with respect to the employed basis set. We find the combination of the infinitely delocalized plane waves and local short-ranged geminals provides a complementary, and rapidly convergent basis for the description of periodic wavefunctions. We hope that this approach will expand the scope of discrete wavefunction expansions in periodic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Many-body quantum chemistry for the electron gas: convergent perturbative theories

We investigate the accuracy of a number of wavefunction based methods at the heart of quantum chemistry for metallic systems. Using Hartree-Fock as a reference, perturbative (Møller-Plesset, MP) and coupled cluster (CC) theories are used to study the uniform electron gas model. Our findings suggest that non-perturbative coupled cluster theories are acceptable for modelling electronic interactions in metals whilst perturbative coupled cluster theories are not. Using screened interactions, we propose a simple modification to the widely-used coupled-cluster singles and doubles plus perturbative triples method (CCSD(T)) that lifts the divergent behaviour and is shown to give very accurate correlation energies for the homogeneous electron gas.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence of many-body wavefunction expansions using a plane wave basis: from the homogeneous electron gas to the solid state

Using the finite simulation-cell homogeneous electron gas (HEG) as a model, we investigate the convergence of the correlation energy to the complete basis set (CBS) limit in methods utilising plane-wave wavefunction expansions. Simple analytic and numerical results from second-order Møller-Plesset theory (MP2) suggest a 1/M decay of the basis-set incompleteness error where M is the number of plane waves used in the calculation, allowing for straightforward extrapolation to the CBS limit. As we shall show, the choice of basis set truncation when constructing many-electron wavefunctions is far from obvious, and here we propose several alternatives based on the momentum transfer vector, which greatly improve the rate of convergence. This is demonstrated for a variety of wavefunction methods, from MP2 to coupled-cluster doubles theory (CCD) and the random-phase approximation plus second-order screened exchange (RPA+SOSEX). Finite basis-set energies are presented for these methods and compared with exact benchmarks. A transformation can map the orbitals of a general solid state system onto the HEG plane wave basis and thereby allow application of these methods to more realistic physical problems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Correlation Energy Divergences in Metallic Systems

We numerically examine divergences of the total energy in metallic systems of approximate many-body theories using Hartree--Fock as a reference, including perturbative (M\oller-Plesset, MP), coupled cluster (CC) and configuration interaction (CI) approaches. Controlling for finite size effects and basis set incompleteness error by comparison with energies from the random phase approximation (RPA), we suggest convincingly that non-perturbative coupled cluster theories are acceptable for modelling electronic interactions in metals whilst perturbative coupled cluster theories are not. Data are provided from the RPA with which it is possible to test other approximate finite-basis methods for divergences with only modest computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Emergence of Critical Phenomena in Full Configuration Interaction Quantum Monte Carlo

There has been recent literature discussion on the origin and severity of the `sign problem' in full configuration interaction quantum Monte Carlo (FCIQMC) and its `initiator' adaptation (i-FCIQMC), methods of interest and potential because they allow for exact (FCI) ground-state solutions to be obtained often at a much reduced computational cost. In this study we aim to use a simple order parameter, describing the `sign structure' of the stochastic wavefunction representation, to empirically characterise the fundamentally different collective behaviour of the walker population in both methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Investigation of the Full Configuration Interaction Quantum Monte Carlo Method Using Homogeneous Electron Gas Models

Using the homogeneous electron gas (HEG) as a model, we investigate the sources of error in the `initiator' adaptation to Full Configuration Interaction Quantum Monte Carlo (i-FCIQMC), with a view to accelerating convergence. In particular we find that the fixed shift phase, where the walker number is allowed to grow slowly, can be used to effectively assess stochastic and initiator error. Using this approach we provide simple explanations for the internal parameters of an i-FCIQMC simulation. We exploit the consistent basis sets and adjustable correlation strength of the HEG to analyze properties of the algorithm, and present finite basis benchmark energies for N=14 over a range of densities $0.5 \leq r_s \leq 5.0$ a.u. A \emph{single-point extrapolation} scheme is introduced to produce complete basis energies for 14, 38 and 54 electrons. It is empirically found that, in the weakly correlated regime, the computational cost scales linearly with the plane wave basis set size, which is justifiable on physical grounds. We expect the fixed shift strategy to reduce the computational cost of many \iFCIQMC calculations of weakly correlated systems. In addition, we provide benchmarks for the electron gas, to be used by other quantum chemical methods in exploring periodic solid state systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Full Configuration Interaction Perspective on the Homogeneous Electron Gas

Highly accurate results for the homogeneous electron gas (HEG) have only been achieved to date within a diffusion Monte Carlo (DMC) framework. Here, we introduce a newly developed stochastic technique, Full Configuration Interaction Quantum Monte Carlo (FCIQMC), which samples the exact wavefunction expanded in plane wave Slater determinants. Despite the introduction of a basis set incompleteness error, we obtain a finite-basis energy which is significantly, and variationally lower than any previously published work for the 54-electron HEG at $r_s$ = 0.5 a.u., in a Hilbert space of $10^{108}$ Slater determinants. At this value of $r_s$, as well as of 1.0 a.u., we remove the remaining basis set incompleteness error by extrapolation, yielding results comparable or better than state-of-the-art DMC backflow energies. In doing so, we demonstrate that it is possible to yield highly accurate results with the FCIQMC method in sizable periodic systems.

0 citations0 reviewsNo code linkedOpen access