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

George H. Booth

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

ResearcherUnclaimed source record
cond-mat.str-elphysics.comp-phphysics.chem-phquant-phcond-mat.mtrl-scicond-mat.quant-gasphysics.atom-ph

Catalog footprint

What is connected

26works
7topics
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

26 published item(s)

preprint2022arXiv

The quantum Gaussian process state: A kernel-inspired state with quantum support data

We introduce the quantum Gaussian process state, motivated via a statistical inference for the wave function supported by a data set of unentangled product states. We show that this condenses down to a compact and expressive parametric form, with a variational flexibility shown to be competitive or surpassing established alternatives. The connections of the state to its roots as a Bayesian inference machine as well as matrix product states, also allow for efficient deterministic training of global states from small training data with enhanced generalization, including on application to frustrated spin physics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fully Algebraic and Self-consistent Effective Dynamics in a Static Quantum Embedding

Quantum embedding approaches involve the self-consistent optimization of a local fragment of a strongly correlated system, entangled with the wider environment. The `energy-weighted' density matrix embedding theory (EwDMET) was established recently as a way to systematically control the resolution of the fragment-environment coupling, and allow for true quantum fluctuations over this boundary to be self-consistently optimized within a fully static framework. In this work, we reformulate the algorithm to ensure that EwDMET can be considered equivalent to an optimal and rigorous truncation of the self-consistent dynamics of dynamical mean-field theory (DMFT). A practical limitation of these quantum embedding approaches is often a numerical fitting of a self-consistent object defining the quantum effects. However, we show here that in this formulation, all numerical fitting steps can be entirely circumvented, via an effective Dyson equation in the space of truncated dynamics. This provides a robust and analytic self-consistency for the method, and an ability to systematically and rigorously converge to DMFT from a static, wave function perspective. We demonstrate that this improved approach can solve the correlated dynamics and phase transitions of the Bethe lattice Hubbard model in infinite dimensions, as well as one- and two-dimensional Hubbard models where we clearly show the benefits of this rapidly convergent basis for correlation-driven fluctuations. This systematically truncated description of the effective dynamics of the problem also allows access to quantities such as Fermi liquid parameters and renormalized dynamics, and demonstrates a numerically efficient, systematic convergence to the zero-temperature dynamical mean-field theory limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A wave function perspective and efficient truncation of renormalised second-order perturbation theory

We present an approach to renormalized second-order Green's function perturbation theory (GF2) which avoids all dependency on continuous variables, grids or explicit Green's functions, and is instead formulated entirely in terms of static quantities and wave functions. Correlation effects from MP2 diagrams are iteratively incorporated to modify the underlying spectrum of excitations by coupling the physical system to fictitious auxiliary degrees of freedom, allowing for the single-particle orbitals to delocalize into this additional space. The overall approach is shown to be rigorously $\mathcal{O}[N^5]$, after an appropriate compression of this auxiliary space. This is achieved via a novel scheme which ensures that a desired number of moments of the underlying occupied and virtual spectra are conserved in the compression, allowing a rapid and systematically improvable convergence to the limit of the effective dynamical resolution. The approach is found to then allow for the qualitative description of stronger correlation effects, avoiding the divergences of MP2, as well as its orbital-optimized version. On application to the G1 test set, we find that modifications to only up to the third spectral moment of the underlying spectrum from which the double excitations are built is required for accurate energetics, even in strongly correlated regimes. This is beyond simple self-consistent changes to the density matrix of the system, but far from requiring a description of the full dynamics of the frequency-dependent self-energy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Driven Imposters: Controlling Expectations in Many-Body Systems

We present a framework to control and track the observables of a general solid state system driven by an incident laser field. The main result is a non-linear equation of motion for tracking an observable, together with a constraint on the size of expectations which may be reproduced via tracking. Among other applications, this model provides a potential route to the design of laser fields which cause photo-induced superconductivity in materials above their critical temperature. As a first test, the strategy is used to make the expectation value of the current conform to an arbitrary function under a range of model parameters. Additionally, using two reference spectra for materials in the conducting and insulating regimes respectively, the tracking algorithm is used to make each material mimic the optical spectrum of the other.

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
preprint2019arXiv

Controlling Arbitrary Observables in Correlated Many-body Systems

Here we present an expanded analysis of a model for the manipulation and control of observables in a strongly correlated, many-body system, which was first presented in [McCaul et al., eprint: arXiv:1911.05006]. A field-free, non-linear equation of motion for controlling the expectation value of an essentially arbitrary observable is derived, together with rigorous constraints that determine the limits of controllability. We show that these constraints arise from the physically reasonable assumptions that the system will undergo unitary time evolution, and has enough degrees of freedom for the electrons to be mobile. Furthermore, we give examples of multiple solutions to generating target observable trajectories when the constraints are violated. Ehrenfest theorems are used to further refine the model, and provide a check on the validity of numerical simulations. Finally, the experimental feasibility of implementing the control fields generated by this model is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Direct comparison of many-body methods for realistic electronic Hamiltonians

A large collaboration carefully benchmarks 20 first principles many-body electronic structure methods on a test set of 7 transition metal atoms, and their ions and monoxides. Good agreement is attained between the 3 systematically converged methods, resulting in experiment-free reference values. These reference values are used to assess the accuracy of modern emerging and scalable approaches to the many-electron problem. The most accurate methods obtain energies indistinguishable from experimental results, with the agreement mainly limited by the experimental uncertainties. Comparison between methods enables a unique perspective on calculations of many-body systems of electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Efficient and Stochastic Multireference Perturbation Theory for Large Active Spaces within a Full Configuration Interaction Quantum Monte Carlo Framework

Full Configuration Interaction Quantum Monte Carlo (FCIQMC) has been effectively applied to very large configuration interaction (CI) problems, and was recently adapted for use as an active space solver and combined with orbital optimisation. In this work, we detail an approach within FCIQMC to allow for efficient sampling of fully internally-contracted multireference perturbation theories within the same stochastic framework. Schemes are described to allow for the close control over the resolution of stochastic sampling of the effective higher-body intermediates within the active space. It is found that while CASPT2 seems less amenable to a stochastic reformulation, NEVPT2 is far more stable, requiring a similar number of walkers to converge the NEVPT2 expectation values as to converge the underlying CI problem. We demonstrate the application of the stochastic approach to the computation of NEVPT2 within a (24,24) active space in a biologically relevant system, and show that small numbers of walkers are sufficient for a faithful sampling of the NEVPT2 energy to chemical accuracy, despite the active space already exceeding the limits of practicality for traditional approaches. This raises prospects of an efficient stochastic solver for multireference chemical problems requiring large active spaces, with an accurate treatment of external orbitals.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Frequency-dependent and algebraic bath states for a Dynamical Mean-Field Theory with compact support

We demonstrate an algebraic construction of frequency-dependent bath orbitals which can be used in a robust and rigorously self-consistent DMFT-like embedding method, here called $ω-$DMFT, suitable for use with Hamiltonian-based impurity solvers. These bath orbitals are designed to exactly reproduce the hybridization of the impurity to its environment, while allowing for a systematic expansion of this bath space as impurity interactions couple frequency points. In this way, the difficult non-linear fit of bath parameters necessary for many Hamiltonian-formulation impurity solvers in DMFT is avoided, while the introduction of frequency dependence in this bath space is shown to allow for more compact bath sizes. This has significant potential use with a number of new, emerging Hamiltonian solvers which allow for the embedding of large impurity spaces within a DMFT framework. We present results of the $ω-$DMFT approach for the Hubbard model on the Bethe lattice, a 1D chain, and the 2D square lattice, which show excellent agreement with standard DMFT results, with fewer bath orbitals and more compact support for the hybridization representation in the key impurity model of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Assessment of Multireference Approaches to Explicitly Correlated Full Configuration Interaction Quantum Monte Carlo

The Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method has proved able to provide near-exact solutions to the electronic Schrödinger equation within a finite orbital basis set, without relying on an expansion about a reference state. However, a drawback to the approach is that being based on an expansion of Slater determinants, the FCIQMC method suffers from a basis set incompleteness error that decays very slowly with the size of the employed single particle basis. The FCIQMC results obtained in a small basis set can be improved significantly with explicitly correlated techniques. Here, we present a study that assesses and compares two contrasting `universal' explicitly correlated approaches that fit into the FCIQMC framework; the $[2]_{R12}$ method of Valeev {\em et al.}, and the explicitly correlated canonical transcorrelation approach of Yanai {\em et al}. The former is an {\em a posteriori} internally-contracted perturbative approach, while the latter transforms the Hamiltonian prior to the FCIQMC simulation. These comparisons are made across the 55 molecules of the G1 standard set. We found that both methods consistently reduce the basis set incompleteness, for accurate atomization energies in small basis sets, reducing the error from 28~mE$_{\text h}$ to 3-4~mE$_{\text h}$. While many of the conclusions hold in general for any combination of multireference approaches with these methodologies, we also consider FCIQMC-specific advantages of each approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From plane waves to local Gaussians for the simulation of correlated periodic systems

We present a simple, robust and black-box approach to the implementation and use of local, periodic, atom-centered Gaussian basis functions within a plane wave code, in a computationally efficient manner. The procedure outlined is based on the representation of the Gaussians within a finite bandwidth by their underlying plane wave coefficients. The core region is handled within the projected augment wave framework, by pseudizing the Gaussian functions within a cut-off radius around each nucleus, smoothing the functions so that they are faithfully represented by a plane wave basis with only moderate kinetic energy cutoff. To mitigate the effects of the basis set superposition error and incompleteness at the mean-field level introduced by the Gaussian basis, we also propose a hybrid approach, whereby the complete occupied space is first converged within a large plane wave basis, and the Gaussian basis used to construct a complementary virtual space for the application of correlated methods. We demonstrate that these pseudized Gaussians yield compact and systematically improvable spaces with an accuracy comparable to their non-pseudized Gaussian counterparts. A key advantage of the described method is its ability to efficiently capture and describe electronic correlation effects of weakly bound and low-dimensional systems, where plane waves are not sufficiently compact or able to be truncated without unphysical artefacts. We investigate the accuracy of the pseudized Gaussians for the water dimer interaction, neon solid and water adsorption on a LiH surface, at the level of second-order M\{o}ller--Plesset perturbation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An excited-state approach within full configuration interaction quantum Monte Carlo

We present a new approach to calculate excited states with the full configuration interaction quantum Monte Carlo (FCIQMC) method. The approach uses a Gram-Schmidt procedure, instantaneously applied to the stochastically evolving distributions of walkers, to orthogonalize higher energy states against lower energy ones. It can thus be used to study several of the lowest-energy states of a system within the same symmetry. This additional step is particularly simple and computationally inexpensive, requiring only a small change to the underlying FCIQMC algorithm. No trial wave functions or partitioning of the space is needed. The approach should allow excited states to be studied for systems similar to those accessible to the ground-state method, due to a comparable computational cost. As a first application we consider the carbon dimer in basis sets up to quadruple-zeta quality, and compare to existing results where available.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analytic nuclear forces and molecular properties from full configuration interaction quantum Monte Carlo

Unbiased stochastic sampling of the one- and two-body reduced density matrices is achieved in full configuration interaction quantum Monte Carlo with the introduction of a second, "replica" ensemble of walkers, whose population evolves in imaginary time independently from the first, and which entails only modest additional computational overheads. The matrices obtained from this approach are shown to be representative of full configuration-interaction quality, and hence provide a realistic opportunity to achieve high-quality results for a range of properties whose operators do not necessarily commute with the hamiltonian. A density-matrix formulated quasi-variational energy estimator having been already proposed and investigated, the present work extends the scope of the theory to take in studies of analytic nuclear forces, molecular dipole moments and polarisabilities, with extensive comparison to exact results where possible. These new results confirm the suitability of the sampling technique and, where sufficiently large basis sets are available, achieve close agreement with experimental values, expanding the scope of the method to new areas of investigation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Krylov-projected quantum Monte Carlo

We present an approach to the calculation of arbitrary spectral, thermal and excited state properties within the full configuration interaction quantum Monte Carlo framework. This is achieved via an unbiased projection of the Hamiltonian eigenvalue problem into a space of stochastically sampled Krylov vectors, thus enabling the calculation of real-frequency spectral and thermal properties and avoiding explicit analytic continuation. We use this approach to calculate temperature-dependent properties and one- and two-body spectral functions for various Hubbard models, as well as isolated excited states in ab initio systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Semi-stochastic full configuration interaction quantum Monte Carlo: developments and application

We expand upon the recent semi-stochastic adaptation to full configuration interaction quantum Monte Carlo (FCIQMC). We present an alternate method for generating the deterministic space without a priori knowledge of the wave function and present stochastic efficiencies for a variety of both molecular and lattice systems. The algorithmic details of an efficient semi-stochastic implementation are presented, with particular consideration given to the effect that the adaptation has on parallel performance in FCIQMC. We further demonstrate the benefit for calculation of reduced density matrices in FCIQMC through replica sampling, where the semi-stochastic adaptation seems to have even larger efficiency gains. We then combine these ideas to produce explicitly correlated corrected FCIQMC energies for the Beryllium dimer, for which stochastic errors on the order of wavenumber accuracy are achievable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spectral functions of strongly correlated extended systems via an exact quantum embedding

Density matrix embedding theory (DMET) [Phys. Rev. Lett., 109, 186404 (2012)], introduced a new approach to quantum cluster embedding methods, whereby the mapping of strongly correlated bulk problems to an impurity with finite set of bath states was rigorously formulated to exactly reproduce the entanglement of the ground state. The formalism provided similar physics to dynamical mean-field theory at a tiny fraction of the cost, but was inherently limited by the construction of a bath designed to reproduce ground state, static properties. Here, we generalize the concept of quantum embedding to dynamic properties and demonstrate accurate bulk spectral functions at similarly small computational cost. The proposed spectral DMET utilizes the Schmidt decomposition of a response vector, mapping the bulk dynamic correlation functions to that of a quantum impurity cluster coupled to a set of frequency dependent bath states. The resultant spectral functions are obtained on the real-frequency axis, without bath discretization error, and allows for the construction of arbitrary dynamic correlation functions. We demonstrate the method on the 1D and 2D Hubbard model, where we obtain zero temperature, thermodynamic limit spectral functions, and show the trivial extension to two-particle Green functions. This advance therefore extends the scope and applicability of DMET in condensed matter problems as a computationally tractable route to correlated spectral functions of extended systems, and provides a competitive alternative to dynamical mean-field theory for dynamic quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochastic multi-configurational self-consistent field theory

The multi-configurational self-consistent field theory is considered the standard starting point for almost all multireference approaches required for strongly-correlated molecular problems. The limitation of the approach is generally given by the number of strongly-correlated orbitals in the molecule, as its cost will grow exponentially with this number. We present a new multi-configurational self-consistent field approach, wherein linear determinant coefficients of a multi-configurational wavefunction are optimized via the stochastic full configuration interaction quantum Monte Carlo technique at greatly reduced computational cost, with non-linear orbital rotation parameters updated variationally based on this sampled wavefunction. This extends this approach to strongly-correlated systems with far larger active spaces than it is possible to treat by conventional means. By comparison with this traditional approach, we demonstrate that the introduction of stochastic noise in both the determinant amplitudes and the gradient and Hessian of the orbital rotations does not preclude robust and reliable convergence of the orbital optimization. It can even improve the ability to avoid convergence to local minima in the orbital space, and therefore aid in finding variationally lower-energy solutions. We consider the effect on the convergence of the orbitals as the number of walkers and the sampling time within the active space increases, as well as the effect on the final energy and error. The scope of the new protocol is demonstrated with a study of the increasingly strongly correlated electronic structure in a series of polycyclic aromatic hydrocarbons, up to the large coronene molecule in a complete active space of 24 pi electrons in 24 orbitals, requiring only modest computational resources.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectroscopic accuracy directly from quantum chemistry: application to ground and excited states of beryllium dimer

We combine explicit correlation via the canonical transcorrelation approach with the density matrix renormalization group and initiator full configuration interaction quantum Monte Carlo methods to compute a near-exact beryllium dimer curve, {\it without} the use of composite methods. In particular, our direct density matrix renormalization group calculations produce a well-depth of $D_e$=931.2 cm$^{-1}$ which agrees very well with recent experimentally derived estimates $D_e$=929.7$\pm 2$~cm$^{-1}$ [Science, 324, 1548 (2009)] and $D_e$=934.6~cm$^{-1}$ [Science, 326, 1382 (2009)]], as well the best composite theoretical estimates, $D_e$=938$\pm 15$~cm$^{-1}$ [J. Phys. Chem. A, 111, 12822 (2007)] and $D_e$=935.1$\pm 10$~cm$^{-1}$ [Phys. Chem. Chem. Phys., 13, 20311 (2011)]. Our results suggest possible inaccuracies in the functional form of the potential used at shorter bond lengths to fit the experimental data [Science, 324, 1548 (2009)]. With the density matrix renormalization group we also compute near-exact vertical excitation energies at the equilibrium geometry. These provide non-trivial benchmarks for quantum chemical methods for excited states, and illustrate the surprisingly large error that remains for 1$^1Σ^-_g$ state with approximate multi-reference configuration interaction and equation-of-motion coupled cluster methods. Overall, we demonstrate that explicitly correlated density matrix renormalization group and initiator full configuration interaction quantum Monte Carlo methods allow us to fully converge to the basis set and correlation limit of the non-relativistic Schrödinger equation in small molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The intermediate and spin-liquid phase of the half-filled honeycomb Hubbard model

We obtain the phase-diagram of the half-filled honeycomb Hubbard model with density matrix embedding theory, to address recent controversy at intermediate couplings. We use clusters from 2-12 sites and lattices at the thermodynamic limit. We identify a paramagnetic insulating state, with possible hexagonal cluster order, competitive with the antiferromagnetic phase at intermediate coupling. However, its stability is strongly cluster and lattice size dependent, explaining controver- sies in earlier work. Our results support the paramagnetic insulator as being a metastable, rather than a true, intermediate phase, in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unbiased Reduced Density Matrices and Electronic Properties from Full Configuration Interaction Quantum Monte Carlo

Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the Full Configuration Interaction Quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself), and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems.

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

Linear-scaling and parallelizable algorithms for stochastic quantum chemistry

For many decades, quantum chemical method development has been dominated by algorithms which involve increasingly complex series of tensor contractions over one-electron orbital spaces. Procedures for their derivation and implementation have evolved to require the minimum amount of logic and rely heavily on computationally efficient library-based matrix algebra and optimized paging schemes. In this regard, the recent development of exact stochastic quantum chemical algorithms to reduce computational scaling and memory overhead requires a contrasting algorithmic philosophy, but one which when implemented efficiently can often achieve higher accuracy/cost ratios with small random errors. Additionally, they can exploit the continuing trend for massive parallelization which hinders the progress of deterministic high-level quantum chemical algorithms. In the Quantum Monte Carlo community, stochastic algorithms are ubiquitous but the discrete Fock space of quantum chemical methods is often unfamiliar, and the methods introduce new concepts required for algorithmic efficiency. In this paper, we explore these concepts and detail an algorithm used for Full Configuration Interaction Quantum Monte Carlo (FCIQMC), which is implemented and available in MOLPRO and as a standalone code, and is designed for high-level parallelism and linear-scaling with walker number. Many of the algorithms are also in use in, or can be transferred to, other stochastic quantum chemical methods and implementations. We apply these algorithms to the strongly correlated Chromium dimer, to demonstrate their efficiency and parallelism.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An explicitly correlated approach to basis set incompleteness in Full Configuration Interaction Quantum Monte Carlo

By performing a stochastic dynamic in a space of Slater determinants, the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method has been able to obtain energies which are essentially free from systematic error to the basis set correlation energy, within small and systematically improvable errorbars. However, the weakly exponential scaling with basis size makes converging the energy with respect to basis set costly and in larger systems, impossible. To ameliorate these basis set issues, here we use perturbation theory to couple the FCIQMC wave function to an explicitly correlated strongly orthogonal basis of geminals, following the [2]_{\textrm{R12}} approach of Valeev {\em et al.}. The required one- and two-particle density matrices are computed on-the-fly during the FCIQMC dynamic, using a sampling procedure which incurs relatively little additional computation expense. The F12 energy corrections are shown to converge rapidly as a function of sampling, both in imaginary time, and number of walkers. Our pilot calculations on the binding curve for carbon dimer, which exhibits strong correlation effects as well as substantial basis set dependence, demonstrate that the accuracy of the FCIQMC-F12 method surpasses that of all previous FCIQMC calculations, and that the F12 correction improves accuracy equivalent to increasing the quality of the one-electron basis by two cardinal numbers.

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

Excited States, Dynamic Correlation Functions and Spectral Properties from Full Configuration Interaction Quantum Monte Carlo

In this communication, we propose a method for obtaining isolated excited states within the Full Configuration Interaction Quantum Monte Carlo framework. This method allows for stable sampling with respect to collapse to lower energy states and requires no uncontrolled approximations. In contrast with most previous methods to extract excited state information from Quantum Monte Carlo methods, this results from a modification to the underlying propagator, and does not require explicit orthogonalization, analytic continuation, transient estimators or restriction of the Hilbert space via a trial wavefunction. Furthermore, we show that the propagator can directly yield frequency-domain correlation functions and spectral functions such as the density of states which are difficult to obtain within a traditional Quantum Monte Carlo framework. We demonstrate this approach with pilot applications to the neon atom and beryllium dimer.

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