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

N. D. Drummond

N. D. Drummond contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-scicond-mat.quant-gascond-mat.mes-hallcond-mat.str-elphysics.comp-phcond-mat.otherquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
7topics
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

14 published item(s)

preprint2022arXiv

Low-density Phase Diagram of the Three-Dimensional Electron Gas

Variational and diffusion quantum Monte Carlo methods are employed to investigate the zero-temperature phase diagram of the three-dimensional homogeneous electron gas at very low density. Fermi fluid and body-centered cubic Wigner crystal ground state energies are determined using Slater-Jastrow-backflow and Slater-Jastrow many-body wave functions at different densities and spin polarizations in finite simulation cells. Finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Unlike previous studies, our results show that the electron gas undergoes a first-order quantum phase transition directly from a paramagnetic fluid to a body-centered cubic crystal at density parameter $r_\text{s} = 86.6(7)$, with no region of stability for an itinerant ferromagnetic fluid. However there is a possible magnetic phase transition from an antiferromagnetic crystal to a ferromagnetic crystal at $r_\text{s}=93(3)$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Point defect formation energies in graphene from diffusion quantum Monte Carlo and density functional theory

Density functional theory (DFT) is widely used to study defects in monolayer graphene with a view to applications ranging from water filtration to electronics to investigation of radiation damage in graphite moderators. To assess the accuracy of DFT in such applications, we report diffusion quantum Monte Carlo (DMC) calculations of the formation energies of some common and important point defects in monolayer graphene: monovacancies, Stone-Wales defects, and silicon substitutions. We find that standard DFT methods underestimate monovacancy formation energies by around 1 eV. The disagreement between DFT and DMC is somewhat smaller for Stone-Wales defects and silicon substitutions. We examine vibrational contributions to the free energies of formation for these defects, finding that vibrational effects are non-negligible. Finally, we compare the DMC atomization energies of monolayer graphene, monolayer silicene, and bulk silicon, finding that bulk silicon is significantly more stable than monolayer silicene by 0.7522(5) eV per atom.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Diffusion quantum Monte Carlo and GW study of the electronic properties of monolayer and bulk hexagonal boron nitride

We report diffusion quantum Monte Carlo (DMC) and many-body $GW$ calculations of the electronic band gaps of monolayer and bulk hexagonal boron nitride (hBN). We find the monolayer band gap to be indirect. $GW$ predicts much smaller quasiparticle gaps at both the single-shot $G_0W_0$ and the partially self-consistent $GW_0$ levels. In contrast, solving the Bethe-Salpeter equation on top of the $GW_0$ calculation yields an exciton binding energy for the direct exciton at the $K$ point in close agreement with the DMC value. Vibrational renormalization of the electronic band gap is found to be significant in both the monolayer and the bulk. Taking vibrational effects into account, DMC overestimates the band gap of bulk hBN, while $GW$ theory underestimates it.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ground state properties of the one-dimensional electron liquid

We present calculations of the energy, pair correlation function (PCF), static structure factor (SSF), and momentum density (MD) for the one-dimensional electron gas using the quantum Monte Carlo method. We are able to resolve peaks in the SSF at even-integer-multiples of the Fermi wave vector, which grow as the coupling is increased. Our MD results show an increase in the effective Fermi wave vector as the interaction strength is raised in the paramagnetic harmonic wire; this appears to be a result of the vanishing difference between the wave functions of the paramagnetic and ferromagnetic systems. We have extracted the Luttinger liquid exponent from our MDs by fitting to data around the Fermi wave vector, finding good agreement between the exponent of the ferromagnetic infinitely-thin wire and the ferromagnetic harmonic wire.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum Monte Carlo calculation of the zero-temperature phase diagram of the two-component fermionic hard-core gas in two dimensions

Motivated by potential realizations in cold-atom or cold-molecule systems, we have performed quantum Monte Carlo simulations of two-component gases of fermions in two dimensions with hard-core interactions. We have determined the gross features of the zero-temperature phase diagram, by investigating the relative stabilities of paramagnetic and ferromagnetic fluids and crystals. We have also examined the effect of including a pairwise, long-range r^(-3) potential between the particles. Our most important conclusion is that there is no region of stability for a ferromagnetic fluid phase, even if the long-range interaction is present. We also present results for the pair-correlation function, static structure factor, and momentum density of two-dimensional hard-core fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Bulk and surface energetics of lithium hydride crystal: benchmarks from quantum Monte Carlo and quantum chemistry

We show how accurate benchmark values of the surface formation energy of crystalline lithium hydride can be computed by the complementary techniques of quantum Monte Carlo (QMC) and wavefunction-based molecular quantum chemistry. To demonstrate the high accuracy of the QMC techniques, we present a detailed study of the energetics of the bulk LiH crystal, using both pseudopotential and all-electron approaches. We show that the equilibrium lattice parameter agrees with experiment to within 0.03 %, which is around the experimental uncertainty, and the cohesive energy agrees to within around 10 meV per formula unit. QMC in periodic slab geometry is used to compute the formation energy of the LiH (001) surface, and we show that the value can be accurately converged with respect to slab thickness and other technical parameters. The quantum chemistry calculations build on the recently developed hierarchical scheme for computing the correlation energy of a crystal to high precision. We show that the hierarchical scheme allows the accurate calculation of the surface formation energy, and we present results that are well converged with respect to basis set and with respect to the level of correlation treatment. The QMC and hierarchical results for the surface formation energy agree to within about 1 %.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Continuum variational and diffusion quantum Monte Carlo calculations

This topical review describes the methodology of continuum variational and diffusion quantum Monte Carlo calculations. These stochastic methods are based on many-body wave functions and are capable of achieving very high accuracy. The algorithms are intrinsically parallel and well-suited to petascale computers, and the computational cost scales as a polynomial of the number of particles. A guide to the systems and topics which have been investigated using these methods is given. The bulk of the article is devoted to an overview of the basic quantum Monte Carlo methods, the forms and optimisation of wave functions, performing calculations within periodic boundary conditions, using pseudopotentials, excited-state calculations, sources of calculational inaccuracy, and calculating energy differences and forces.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Diamond -> beta-tin phase transition in Si within diffusion quantum Monte Carlo

We have studied the diamond -> beta-tin phase transition in Si using diffusion quantum Monte Carlo (DMC) methods. Slater-Jastrow-backflow trial wave functions give lower DMC energies than Slater-Jastrow ones, and backflow slightly favors the beta-tin phase with respect to the diamond phase. We have investigated the changes in the equation of state that result from the use of different pseudopotentials, the inclusion of either zero-point motion or finite-temperature vibrations, and the application of corrections for finite-size effects. Our tests indicate that the choice of pseudopotential can significantly affect the equation of state. Using a Dirac-Fock pseudopotential leads to an overestimation of the transition pressure but an empirical pseudopotential designed for use in correlated calculations gives a transition pressure in quite good agreement with experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

First-Principles Study of a Positron Immersed in an Electron Gas

Calculations of the relaxation energy, contact pair-correlation function, and annihilating-pair momentum density for a single positron immersed in a homogeneous electron gas are presented. We achieve an accurate description of the electron-positron correlation effects by working in the reference frame in which the positron is stationary and using a mean-field approach based on single-component density functional theory. Our positron relaxation energies and annihilation rates are similar to those from the best existing many-body calculations. Our annihilating-pair momentum densities are significantly different from previous data, and include a "tail" beyond the Fermi edge.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase Diagram of the Low-Density Two-Dimensional Homogeneous Electron Gas

We have used quantum Monte Carlo methods to calculate the zero-temperature phase diagram of the two-dimensional homogeneous electron gas. We find a transition from a paramagnetic fluid to an antiferromagnetic triangular Wigner crystal at density parameter r_s=31(1) a.u. and a transition to a ferromagnetic crystal at r_s=38(5) a.u. The fully spin-polarized fluid is never stable. We searched for, but did not find, the ferromagnetic hybrid phase proposed by H. Falakshahi and X. Waintal [Phys. Rev. Lett. 94, 046801 (2005)].

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum Monte Carlo calculation of the energy band and quasiparticle effective mass of the two-dimensional Fermi fluid

We have used the diffusion quantum Monte Carlo method to calculate the energy band of the two-dimensional homogeneous electron gas (HEG), and hence we have obtained the quasiparticle effective mass and the occupied bandwidth. We find that the effective mass in the paramagnetic HEG increases significantly when the density is lowered, whereas it decreases in the fully ferromagnetic HEG. Our calculations therefore support the conclusions of recent experimental studies [Y.-W. Tan et al., Phys. Rev. Lett. 94, 016405 (2005); M. Padmanabhan et al., Phys. Rev. Lett. 101, 026402 (2008); T. Gokmen et al., Phys. Rev. B 79, 195311 (2009)]. We compare our calculated effective masses with other theoretical results and experimental measurements in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid

We have used the variational and diffusion quantum Monte Carlo methods to calculate the energy, pair correlation function, static structure factor, and momentum density of the ground state of the two-dimensional homogeneous electron gas. We have used highly accurate Slater-Jastrow-backflow trial wave functions and twist averaging to reduce finite-size effects where applicable. We compare our results with others in the literature and construct a local-density-approximation exchange-correlation functional for 2D systems.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Exciton-exciton interaction and biexciton formation in bilayer systems

We report quantum Monte Carlo calculations of biexciton binding energies in ideal two-dimensional bilayer systems with isotropic electron and hole masses. We have also calculated exciton-exciton interaction potentials, and pair distribution functions for electrons and holes in bound biexcitons. Comparing our data with results obtained in a recent study using a model exciton-exciton potential [C. Schindler and R. Zimmermann, Phys. Rev. B \textbf{78}, 045313 (2008)], we find a somewhat larger range of layer separations at which biexcitons are stable. We find that individual excitons retain their identity in bound biexcitons for large layer separations.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A variance-minimization scheme for optimizing Jastrow factors

We describe a new scheme for optimizing many-electron trial wave functions by minimizing the unreweighted variance of the energy using stochastic integration and correlated-sampling techniques. The scheme is restricted to parameters that are linear in the exponent of a Jastrow correlation factor, which are the most important parameters in the wave functions we use. The scheme is highly efficient and allows us to investigate the parameter space more closely than has been possible before. We search for multiple minima of the variance in the parameter space and compare the wave functions obtained using reweighted and unreweighted variance minimization.

0 citations0 reviewsNo code linkedOpen access