BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

J. S. Spencer

J. S. Spencer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elphysics.chem-phphysics.comp-phcond-mat.mtrl-scicond-mat.stat-mechSoftware Engineering

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

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

Understanding and Improving the Efficiency of Full Configuration Interaction Quantum Monte Carlo

Within Full Configuration Interaction Quantum Monte Carlo, we investigate how the statistical error behaves as a function of the parameters which control the stochastic sampling. We define the inefficiency as a measure of the statistical error per particle sampling the space and per timestep and show there is a sizeable parameter regime where this is minimised. We find that this inefficiency increases sublinearly with Hilbert space size and can be reduced by localising the canonical Hartree--Fock molecular orbitals, suggesting that the choice of basis impacts the method beyond that of the sign problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hubbard-like Hamiltonians for interacting electrons in s, p and d orbitals

Hubbard-like Hamiltonians are widely used to describe on-site Coulomb interactions in magnetic and strongly-correlated solids, but there is much confusion in the literature about the form these Hamiltonians should take for shells of p and d orbitals. This paper derives the most general s, p and d orbital Hubbard-like Hamiltonians consistent with the relevant symmetries, and presents them in ways convenient for practical calculations. We use the full configuration interaction method to study p and d orbital dimers and compare results obtained using the correct Hamiltonian and the collinear and vector Stoner Hamiltonians. The Stoner Hamiltonians can fail to describe properly the nature of the ground state, the time evolution of excited states, and the electronic heat capacity.

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

Minimising biases in Full Configuration Interaction Quantum Monte Carlo

We show that Full Configuration Interaction Quantum Monte Carlo (FCIQMC) is a Markov Chain in its present form. We construct the Markov matrix of FCIQMC for a two determinant system and hence compute the stationary distribution. These solutions are used to quantify the dependence of the population dynamics on the parameters defining the Markov chain. Despite the simplicity of a system with only two determinants, it still reveals a population control bias inherent to the FCIQMC algorithm. We investigate the effect of simulation parameters on the population control bias for the neon atom and suggest simulation setups to in general minimise the bias. We show a reweighting scheme to remove the bias caused by population control commonly used in Diffusion Monte Carlo [J. Chem. Phys. 99, 2865 (1993)] is effective and recommend its use as a post processing step.

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

Density matrix quantum Monte Carlo

We present a quantum Monte Carlo method capable of sampling the full density matrix of a many-particle system at finite temperature. This allows arbitrary reduced density matrix elements and expectation values of complicated non-local observables to be evaluated easily. The method resembles full configuration interaction quantum Monte Carlo but works in the space of many-particle operators instead of the space of many-particle wave functions. One simulation provides the density matrix at all temperatures simultaneously, from $T=\infty$ to $T=0$, allowing the temperature dependence of expectation values to be studied. The direct sampling of the density matrix also allows the calculation of some previously inaccessible entanglement measures. We explain the theory underlying the method, describe the algorithm, and introduce an importance-sampling procedure to improve the stochastic efficiency. To demonstrate the potential of our approach, the energy and staggered magnetization of the isotropic antiferromagnetic Heisenberg model on small lattices, the concurrence of one-dimensional spin rings, and the Renyi $S_2$ entanglement entropy of various sublattices of the $6 \times 6$ Heisenberg model are calculated. The nature of the sign problem in the method is also investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The sign problem and population dynamics in the full configuration interaction quantum Monte Carlo method

The recently proposed full configuration interaction quantum Monte Carlo method allows access to essentially exact ground-state energies of systems of interacting fermions substantially larger than previously tractable without knowledge of the nodal structure of the ground-state wave function. We investigate the nature of the sign problem in this method and how its severity depends on the system studied. We explain how cancelation of the positive and negative particles sampling the wave function ensures convergence to a stochastic representation of the many-fermion ground state and accounts for the characteristic population dynamics observed in simulations.

0 citations0 reviewsNo code linkedOpen access