Researcher profile

R. J. Needs

R. J. Needs contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
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

12 published item(s)

preprint2010arXiv

Benchmark all-electron ab initio quantum Monte Carlo calculations for small molecules

We study the efficiency, precision and accuracy of all-electron variational and diffusion quantum Monte Carlo calculations using Slater basis sets. Starting from wave functions generated by Hartree-Fock and density functional theory, we describe an algorithm to enforce the electron-nucleus cusp condition by linear projection. For the 55 molecules in the G2 set, the diffusion quantum Monte Carlo calculations recovers an average of 95% of the correlation energy and reproduces bond energies to a mean absolute deviation of 3.2 kcal/mol. Comparing the individual total energies with essentially exact values, we investigate the error cancellation in atomization and chemical reaction path energies, giving additional insight into the sizes of nodal surface errors.

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.

preprint2010arXiv

Methods for calculating forces within quantum Monte Carlo simulations

Atomic force calculations within the variational and diffusion quantum Monte Carlo (VMC and DMC) methods are described. The advantages of calculating DMC forces with the "pure" rather than the "mixed" probability distribution are discussed. An accurate and practical method for calculating forces using the pure distribution is presented and tested for the SiH molecule. The statistics of force estimators are explored and violations of the Central Limit Theorem are found in some cases.

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

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.

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.

preprint2009arXiv

Hydrogen/nitrogen/oxygen defect complexes in silicon from computational searches

Point defect complexes in crystalline silicon composed of hydrogen, nitrogen, and oxygen atoms are studied within density-functional theory (DFT). Ab initio Random Structure Searching (AIRSS) is used to find low-energy defect structures. We find new lowest-energy structures for several defects: the triple-oxygen defect, {3O}, triple oxygen with a nitrogen atom, {N, 3O}, triple nitrogen with an oxygen atom, {3N,O}, double hydrogen and an oxygen atom, {2H,O}, double hydrogen and oxygen atoms, {2H,2O} and four hydrogen/nitrogen/oxygen complexes, {H,N,O}, {2H,N,O}, {H,2N,O} and {H,N,2O}. We find that some defects form analogous structures when an oxygen atom is replaced by a NH group, for example, {H,N,2O} and {3O}, and {H,N} and {O}. We compare defect formation energies obtained using different oxygen chemical potentials and investigate the relative abundances of the defects.

preprint2009arXiv

Post-Cotunnite phase of TeO2

We have used first-principles density-functional-theory methods with a random-structure-searching technique to determine the structure of the previously unidentified post-cotunnite phase of TeO2. Our calculations indicate a transition from the cotunnite to post-cotunnite phase at 130 GPa. The predicted post-cotunnite structure has P2_1/m space group symmetry and its calculated x-ray diffraction pattern is in good agreement with the available experimental data. We find that the cotunnite phase re-enters at about 260 GPa.

preprint2009arXiv

Ultracold atoms at unitarity within quantum Monte Carlo

Variational and diffusion quantum Monte Carlo (VMC and DMC) calculations of the properties of the zero-temperature fermionic gas at unitarity are reported. The ratio of the energy of the interacting to the non-interacting gas for a system of 128 particles is calculated to be 0.4517(3) in VMC and 0.4339(1) in the more accurate DMC method. The spherically-averaged pair-correlation functions, momentum densities, and one-body density matrices are very similar in VMC and DMC, but the two-body density matrices and condensate fractions show some differences. Our best estimate of the condensate fraction of 0.51 is a little smaller than values from other quantum Monte Carlo calculations.

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.

preprint2008arXiv

Hydrogen/silicon complexes in silicon from computational searches

Defects in crystalline silicon consisting of a silicon self-interstitial atom and one, two, three, or four hydrogen atoms are studied within density-functional theory (DFT). We search for low-energy defects by starting from an ensemble of structures in which the atomic positions in the defect region have been randomized. We then relax each structure to a minimum in the energy. We find a new defect consisting of a self-interstitial and one hydrogen atom (denoted by {I,H}) which has a higher symmetry and a lower energy than previously reported structures. We recover the {I,H_2} defect found in previous studies and confirm that it is the most stable such defect. Our best {I,H_3} defect has a slightly different structure and lower energy than the one previously reported, and our lowest energy {I,H_4} defect is different to those of previous studies.

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.