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

R. J. Needs

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

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.quant-gascond-mat.mes-hallcond-mat.str-elphysics.comp-phphysics.chem-phastro-ph.SRcond-mat.stat-mechcond-mat.supr-con

Catalog footprint

What is connected

44works
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

44 published item(s)

preprint2016arXiv

Pseudopotential for the 2D contact interaction

We propose a smooth pseudopotential for the contact interaction acting between ultracold atoms confined to two dimensions. The pseudopotential reproduces the scattering properties of the repulsive contact interaction up to 200 times more accurately than a hard disk potential, and in the attractive branch gives a 10-fold improvement in accuracy over the square well potential. Furthermore, the new potential enables diffusion Monte Carlo simulations of the ultracold gas to be run 15 times quicker than was previously possible.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Trail-Needs pseudopotentials in quantum Monte Carlo calculations with plane-wave/blip basis sets

We report a systematic analysis of the performance of a widely used set of Dirac-Fock pseudopotentials for quantum Monte Carlo (QMC) calculations. We study each atom in the periodic table from hydrogen (Z=1) to mercury (Z=80), with the exception of the 4f elements (57 <= Z <= 70). We demonstrate that ghost states are a potentially serious problem when plane-wave basis sets are used in density functional theory (DFT) orbital-generation calculations, but that this problem can be almost entirely eliminated by choosing the s channel to be local in the DFT calculation; the d channel can then be chosen to be local in subsequent QMC calculations, which generally leads to more accurate results. We investigate the achievable energy variance per electron with different levels of trial wave function and we determine appropriate plane-wave cutoff energies for DFT calculations for each pseudopotential. We demonstrate that the so-called "T-move" scheme in diffusion Monte Carlo is essential for many elements. We investigate the optimal choice of spherical integration rule for pseudopotential projectors in QMC calculations. The information reported here will prove crucial in the planning and execution of QMC projects involving beyond-first-row elements.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Trion formation in a two-dimensional hole-doped electron gas

The interaction between a single hole and a two-dimensional, paramagnetic, homogeneous electron gas is studied using diffusion quantum Monte Carlo simulations. Calculations of the electron-hole correlation energy, pair-correlation function, and the electron-hole center-of-mass momentum density are reported for a range of electron--hole mass ratios and electron densities. We find numerical evidence of a crossover from a collective Mahan exciton to a trion-dominated state in a density range in agreement with that found in recent experiments on quantum well heterostructures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computational searches for iron oxides at high pressures

We have used density-functional-theory methods and the ab initio random structure searching (AIRSS) approach to predict stable structures and stoichiometries of mixtures of iron and oxygen at high pressures. Searching was performed for 12 different stoichiometries at pressures of 100, 350 and 500 GPa, which involved relaxing more than 32,000 structures. We find that Fe$_2$O$_3$ and FeO$_2$ are the only phases stable to decomposition at 100 GPa, while at 350 and 500 GPa several stoichiometries are found to be stable or very nearly stable. We report a new structure of Fe$_2$O$_3$ with $P2_12_12_1$ symmetry which is found to be more stable than the known Rh$_2$O$_3$(II) phase at pressures above $\sim$233 GPa. We also report two new structures of FeO, with $Pnma$ and $R\bar{3}m$ symmetries, which are found to be stable within the ranges 195-285 GPa and 285-500 GPa, respectively, and two new structures of Fe$_3$O$_4$ with $Pca2_1$ and $P2_1/c$ symmetries, which are found to be stable within the ranges 100-340 GPa and 340-500 GPa, respectively. Finally, we report two new structures of Fe$_4$O$_5$ with $P4_2/n$ and $P\bar{3}m1$ symmetries, which are found to be stable within the ranges 100-231 GPa and 231-500 GPa, respectively. Our new structures of Fe$_3$O$_4$ and Fe$_4$O$_5$ are found to have lower enthalpies than their known structures within their respective stable pressure ranges.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pseudopotential for the electron-electron interaction

We propose a pseudopotential for the electron-electron Coulomb interaction to improve the efficiency of many-body electronic structure calculations. The pseudopotential accurately replicates the scattering properties of the Coulomb interaction, and recovers the analytical solution for two electrons in a parabolic trap. A case study for the homogeneous electron gas using the diffusion Monte Carlo and configuration interaction methods recovers highly accurate values for the ground state energy, and the smoother potential reduces the computational cost by a factor of ~30. Finally, we demonstrate the use of the pseudopotential to study isolated lithium and beryllium atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures

Establishing the phase diagram of hydrogen is a major challenge for experimental and theoretical physics. Experiment alone cannot establish the atomic structure of solid hydrogen at high pressure, because hydrogen scatters X-rays only weakly. Instead our understanding of the atomic structure is largely based on density functional theory (DFT). By comparing Raman spectra for low-energy structures found in DFT searches with experimental spectra, candidate atomic structures have been identified for each experimentally observed phase. Unfortunately, DFT predicts a metallic structure to be energetically favoured at a broad range of pressures up to 400 GPa, where it is known experimentally that hydrogen is nonmetallic. Here we show that more advanced theoretical methods (diffusion quantum Monte Carlo calculations) find the metallic structure to be uncompetitive, and predict a phase diagram in reasonable agreement with experiment. This greatly strengthens the claim that the candidate atomic structures accurately model the experimentally observed phases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparing electron-phonon coupling strength in diamond, silicon and silicon carbide: First-principles study

Renormalization of the electronic band gap due to electron-phonon coupling in the tetrahedral semiconductors diamond, silicon and cubic silicon carbide is studied from first principles. There is a marked difference between the coupling of the vibrational state to the valence band maximum and to the conduction band minimum. The strength of phonon coupling to the valence band maximum is similar between the three systems and is dominated by vibrations that change the bond length. The coupling strength to the conduction band minimum differs significantly in diamond, silicon carbide and silicon. In diamond, the coupling is dominated by six small pockets of vibrational states in the phonon Brillouin zone, that are ultimately responsible for the stronger electron-phonon coupling in this material. Our results represent a first step towards the development of an a priori understanding of electron-phonon coupling in semiconductors and insulators, that should aid the design of materials with tailored electron-phonon coupling properties.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dissociation of high-pressure solid molecular hydrogen: Quantum Monte Carlo and anharmonic vibrational study

A theoretical study is reported of the molecular-to-atomic transition in solid hydrogen at high pressure. We use the diffusion quantum Monte Carlo method to calculate the static lattice energies of the competing phases and a density-functional-theory-based vibrational self-consistent field method to calculate anharmonic vibrational properties. We find a small but significant contribution to the vibrational energy from anharmonicity. A transition from the molecular Cmca-12 direct to the atomic I4_1/amd phase is found at 374 GPa. The vibrational contribution lowers the transition pressure by 91 GPa. The dissociation pressure is not very sensitive to the isotopic composition. Our results suggest that quantum melting occurs at finite temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Finite-size errors in continuum quantum Monte Carlo calculations

We analyze the problem of eliminating finite-size errors from quantum Monte Carlo (QMC) energy data. We demonstrate that both (i) adding a recently proposed [S. Chiesa et al., Phys. Rev. Lett. 97, 076404 (2006)] finite-size correction to the Ewald energy and (ii) using the model periodic Coulomb (MPC) interaction [L. M. Fraser et al., Phys. Rev. B 53, 1814 (1996); P. R. C. Kent et al., Phys. Rev. B 59, 1917 (1999); A. J. Williamson et al., Phys. Rev. B 55, 4851 (1997)] are good solutions to the problem of removing finite-size effects from the interaction energy in cubic systems, provided the exchange-correlation (XC) hole has converged with respect to system size. However, we find that the MPC interaction distorts the XC hole in finite systems, implying that the Ewald interaction should be used to generate the configuration distribution. The finite-size correction of Chiesa et al. is shown to be incomplete in systems of low symmetry. Beyond-leading-order corrections to the kinetic energy are found to be necessary at intermediate and high densities, and we investigate the effect of adding such corrections to QMC data for the homogeneous electron gas. We analyze finite-size errors in two-dimensional systems and show that the leading-order behavior differs from that which has hitherto been supposed. We compare the efficiency of different twist-averaging methods for reducing single-particle finite-size errors and we examine the performance of various finite-size extrapolation formulas. Finally, we investigate the system-size scaling of biases in diffusion QMC.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-fidelity pseudopotentials for the contact interaction

The contact interaction is often used in modeling ultracold atomic gases, although it leads to pathological behavior arising from the divergence of the many-body wavefunction when two particles coalesce. This makes it difficult to use this model interaction in quantum Monte Carlo and other popular numerical methods. Researchers therefore model the contact interaction with pseudopotentials, such as the square well potential, whose scattering properties deviate markedly from those of the contact potential. In this article, we propose a family of pseudopotentials that reproduce the scattering phase shifts of the contact interaction up to a hundred times more accurately than the square well potential. Moreover, the pseudopotentials are smooth, resulting in significant improvements in efficiency when used in numerical calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ionic ammonia ice

We report experimental and theoretical evidence that solid molecular ammonia becomes unstable at room temperature and high pressures and transforms into an ionic crystalline form. This material has been characterised in both hydrogenated NH3 and deuterated ND3 ammonia samples up to about 180 and 200 GPa, respectively, by infrared absorption, Raman spectroscopy and x-ray diffraction. The presence of a new, strong IR absorption band centered at 2500 cm-1 in NH3 (1900 cm-1 in ND3) signals the ionization of ammonia molecules into NH2- and NH4+ ions, in line with previous theoretical predictions. We find experimental evidence for a coexistence of two crystalline ionic forms, which our ab-initio structure searches predict to be the most stable at the relevant pressures. The ionic crystalline form of ammonia is stable at low temperatures, which contrasts with the behaviour of water in which no equivalent crystalline ionic phase has been found.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anharmonic vibrational properties in periodic systems: energy, electron-phonon coupling, and stress

A unified approach is used to study vibrational properties of periodic systems with first-principles methods and including anharmonic effects. Our approach provides a theoretical basis for the determination of phonon-dependent quantities at finite temperatures. The low-energy portion of the Born-Oppenheimer energy surface is mapped and used to calculate the total vibrational energy including anharmonic effects, electron-phonon coupling, and the vibrational contribution to the stress tensor. We report results for the temperature dependence of the electronic band gap and the linear coefficient of thermal expansion of diamond, lithium hydride, and lithium deuteride.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electron-Phonon Coupling and the Metalization of Solid Helium at Terapascal Pressures

Solid He is studied in the pressure and temperature ranges 1-40 TPa and 0-10,000 K using first-principles methods. Anharmonic vibrational properties are calculated within a self-consistent field framework, including the internal and free energies, density-pressure relation, stress tensor, thermal expansion, and the electron-phonon coupling renormalization of the electronic band gap. We find that an accurate description of electron-phonon coupling requires us to use a non-perturbative approach. The metalization pressure of 32.9 TPa at 0 K is larger than found previously. The vibrational effects are large; for example at P=30 TPa the band gap is increased by 2.8 eV by electron-phonon coupling and a further 0.1 eV by thermal expansion compared to the static value. The implications of the calculated metalization pressure for the cooling of white dwarfs are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lithiation of silicon via lithium Zintl-defect complexes

An extensive search for low-energy lithium defects in crystalline silicon using density-functional-theory methods and the ab initio random structure searching (AIRSS) method shows that the four-lithium-atom substitutional point defect is exceptionally stable. This defect consists of four lithium atoms with strong ionic bonds to the four under-coordinated atoms of a silicon vacancy defect, similar to the bonding of metal ions in Zintl phases. This complex is stable over a range of silicon environments, indicating that it may aid amorphization of crystalline silicon and form upon delithiation of the silicon anode of a Li-ion rechargeable battery.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pseudopotentials for correlated electron systems

A scheme is developed for creating pseudopotentials for use in correlated-electron calculations. Pseudopotentials for the light elements H, Li, Be, B, C, N, O, and F, are reported, based on data from high-level quantum chemical calculations. Results obtained with these correlated electron pseudopotentials (CEPPs) are compared with data for atomic energy levels and the dissociation energies, molecular geometries and zero-point vibrational energies of small molecules obtained from coupled cluster single double triple (CCSD(T)) calculations with large basis sets. The CEPPs give better results in correlated-electron calculations than Hartree-Fock-based pseudopotentials available in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum Monte Carlo calculation of the Fermi liquid parameters of the two-dimensional homogeneous electron gas

Fermi liquid theory is the basic paradigm within which we understand the normal behavior of interacting electron systems, but quantitative values for the parameters that occur in this theory are currently unknown in many important cases. One such case is the two-dimensional homogeneous electron gas (2D HEG), which is realized in a wide variety of semiconductor devices. We have used quantum Monte Carlo (QMC) methods to calculate the Landau interaction functions between pairs of quasiparticles. We use these to study the Fermi liquid parameters, finding that finite-size effects represent a serious obstacle to the direct determination of Fermi liquid parameters in QMC calculations. We have used QMC data in the literature for other properties of the 2D HEG to assemble a set of "best available" values for the Fermi liquid parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Theory of temperature dependent phonon-renormalized properties

We present a general harmonic theory for the temperature dependence of phonon-renormalized properties of solids. Firstly, we formulate a perturbation theory in phonon-phonon interactions to calculate the phonon renormalization of physical quantities. Secondly, we propose two new schemes for extrapolating phonon zero-point corrections from temperature dependent data that improve the accuracy by an order of magnitude compared to previous approaches. Finally, we consider the low-temperature limit of the class of observables that includes the electronic band gap, obtaining a $T^4$ dependence in three dimensions, $T^2$ in two dimensions, and $T^{3/2}$ in one dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diffusion quantum Monte Carlo calculation of the quasiparticle effective mass of the two-dimensional homogeneous electron gas

The quasiparticle effective mass is a key quantity in the physics of electron gases, describing the renormalization of the electron mass due to electron-electron interactions. Two-dimensional electron gases are of fundamental importance in semiconductor physics, and there have been numerous experimental and theoretical attempts to determine the quasiparticle effective mass in these systems. In this work we report quantum Monte Carlo results for the quasiparticle effective mass of a two-dimensional homogeneous electron gas. Our calculations differ from previous quantum Monte Carlo work in that much smaller statistical error bars have been achieved, allowing for an improved treatment of finite-size effects. In some cases we have also been able to use larger system sizes than previous calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Energetics of hydrogen/lithium complexes in silicon analyzed using the Maxwell construction

We have studied hydrogen/lithium complexes in crystalline silicon using density-functional-theory methods and the ab initio random structure searching (AIRSS) method for predicting structures. A method based on the Maxwell construction and convex hull diagrams is introduced which gives a graphical representation of the relative stabilities of point defects in a crystal and enables visualization of the changes in stability when the chemical potentials are altered. We have used this approach to study lithium and hydrogen impurities in silicon, which models aspects of the anode material in the recently-suggested lithium-ion batteries. We show that hydrogen may play a role in these anodes, finding that hydrogen atoms bind to three-atom lithium clusters in silicon, forming stable {H,3Li} and {2H,3Li} complexes, while the {H,2Li} complex is almost stable.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ab initio Random Structure Searching

It is essential to know the arrangement of the atoms in a material in order to compute and understand its properties. Searching for stable structures of materials using first-principles electronic structure methods, such as density functional theory (DFT), is a rapidly growing field. Here we describe our simple, elegant and powerful approach to searching for structures with DFT which we call ab initio random structure searching (AIRSS). Applications to discovering structures of solids, point defects, surfaces, and clusters are reviewed. New results for iron clusters on graphene, silicon clusters, polymeric nitrogen, hydrogen-rich lithium hydrides, and boron are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Predicted formation of superconducting platinum-hydride crystals under pressure in the presence of molecular hydrogen

Noble metals adopt close-packed structures at ambient pressure and rarely undergo structural transformation at high pressures. Platinum (Pt), in particular, is normally considered to be unreactive and is therefore not expected to form hydrides under pressure. We predict that platinum hydride (PtH) has a lower enthalpy than its constituents solid Pt and molecular hydrogen at pressures above 21.5 GPa. We have calculated structural phase transitions from tetragonal to hexagonal close-packed or face-centered cubic (fcc) PtH between 70 and 80 GPa. Linear response calculations indicate that PtH is a superconductor at these pressures with a critical temperature of about 10--25 K. These findings help to shed light on recent observations of pressure-induced metallization and superconductivity in hydrogen-rich materials. We show that formation of fcc metal hydrides under pressure is common among noble metal hydrides and examine the possibility of superconductivity in these materials.

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
preprint2011arXiv

Quantum Monte Carlo Study of a Positron in an Electron Gas

Quantum Monte Carlo calculations of the relaxation energy, pair-correlation function, and annihilating-pair momentum density are presented for a positron immersed in a homogeneous electron gas. We find smaller relaxation energies and contact pair-correlation functions in the important low-density regime than predicted by earlier studies. Our annihilating-pair momentum densities have almost zero weight above the Fermi momentum due to the cancellation of electron-electron and electron-positron correlation effects.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum Monte Carlo study of the first-row atoms and ions

Quantum Monte Carlo calculations of the first-row atoms Li-Ne and their singly-positively-charged ions are reported. Multi-determinant-Jastrow-backflow trial wave functions are used which recover more than 98% of the correlation energy at the Variational Monte Carlo (VMC) level and more than 99% of the correlation energy at the Diffusion Monte Carlo (DMC) level for both the atoms and ions. We obtain the first ionization potentials to chemical accuracy. We also report scalar relativistic corrections to the energies, mass-polarization terms, and one- and two-electron expectation values.

0 citations0 reviewsNo code linkedOpen access
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.

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

Doping driven magnetic instabilities and quantum criticality of NbFe$_{2}$

Using density functional theory we investigate the evolution of the magnetic ground state of NbFe$_{2}$ due to doping by Nb-excess and Fe-excess. We find that non-rigid-band effects, due to the contribution of Fe-\textit{d} states to the density of states at the Fermi level are crucial to the evolution of the magnetic phase diagram. Furthermore, the influence of disorder is important to the development of ferromagnetism upon Nb doping. These findings give a framework in which to understand the evolution of the magnetic ground state in the temperature-doping phase diagram. We investigate the magnetic instabilities in NbFe$_{2}$. We find that explicit calculation of the Lindhard function, $χ_{0}(\mathbf{q})$, indicates that the primary instability is to finite $\mathbf{q}$ antiferromagnetism driven by Fermi surface nesting. Total energy calculations indicate that $\mathbf{q}=0$ antiferromagnetism is the ground state. We discuss the influence of competing $\mathbf{q}=0$ and finite $\mathbf{q}$ instabilities on the presence of the non-Fermi liquid behavior in this material.

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

Fundamental high pressure calibration from all-electron quantum Monte Carlo calculations

We develop an all-electron quantum Monte Carlo (QMC) method for solids that does not rely on pseudopotentials, and use it to construct a primary ultra-high pressure calibration based the equation of state of cubic boron nitride(c-BN). We compute the static contribution to the free energy with QMC, and obtain the phonon contribution from density functional theory, yielding a high-accuracy calibration up to 900 GPa usable directly in experiment. Furthermore, we compute the anharmonic Raman frequency shift with QMC as a function of pressure and temperature, allowing optical pressure calibration in table-top experiments. In contrast to present experimental approaches, small systematic errors in the theoretical EOS do not increase with pressure, and no extrapolation is needed. This all-electron methodology is generally applicable to first-row solids, and can be used to provide a new reference for ab initio calculations of solids and to benchmark pseudopotential accuracy.

0 citations0 reviewsNo code linkedOpen access
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.

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 for minerals at high pressure: Phase stability, equations of state, and elasticity of silica

Silica is an abundant component of the Earth whose crystalline polymorphs play key roles in its structure and dynamics. As the simplest silicates, understanding pure silica is a prerequisite to understanding the rocky part of the Earth, its majority. First principle density functional theory (DFT) methods have often been used to accurately predict properties of silicates. Here, we study silica with quantum Monte Carlo (QMC), which until now was not computationally possible for such complex materials, and find that QMC overcomes the failures of DFT. QMC is a benchmark method that does not rely on density functionals, but rather explicitly treats the electrons and their interactions via a stochastic solution of Schrodinger's equation. Using ground state QMC plus phonons within the quasiharmonic approximation from density functional perturbation theory, we obtain the thermal pressure and equations of state of silica phases up to Earth's core-mantle boundary. Our results provide the most well-constrained equations of state and phase boundaries available for silica. QMC indicates a transition to the most dense alpha-PbO2 structure above the core-insulating D" layer, suggesting the absence of significant free silica in the bulk lower mantle, which has been assumed but never proven. We also find an accurate shear elastic constant and its geophysically important softening with pressure.

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

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

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
preprint2005arXiv

Quantum Monte Carlo calculations of the structural properties and the B1-B2 phase transition of MgO

We report diffusion Monte Carlo (DMC) calculations on MgO in the rock-salt and CsCl structures. The calculations are based on Hartree-Fock pseudopotentials, with the single-particle orbitals entering the correlated wave function being represented by a systematically convergeable cubic-spline basis. Systematic tests are presented on system-size errors using periodically repeating cells of up to over 600 atoms. The equilibrium lattice parameter of the rock-salt structure obtained within DMC is almost identical to the Hartree-Fock result, which is close to the experimental value. The DMC result for the bulk modulus is also in good agreement with the experimental value. The B1-B2 transition pressure (between the rock-salt and CsCl structures) is predicted to be just below 600 GPa, which is beyond the experimentally accessible range, in accord with other predictions based on Hartree-Fock and density functional theories.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Diamond and $β$-tin structures of Si studied with quantum Monte Carlo calculations

We have used diffusion quantum Monte Carlo (DMC) calculations to study the pressure-induced phase transition from the diamond to $β$-tin structure in silicon. The calculations employ the pseudopotential technique and systematically improvable B-spline basis sets. We show that in order to achieve a precision of 1 GPa in the transition pressure the non-cancelling errors in the energies of the two structures must be reduced to 30 meV/atom. Extensive tests on system size errors, non-local pseudopotential errors, basis-set incompleteness errors, and other sources of error, performed on periodically repeated systems of up to 432 atoms, show that all these errors together can be reduced to well below 30 meV/atom. The calculated DMC transition pressure is about 3-4 GPa higher than the accepted experimental range of values, and we argue that the discrepancy may be due to the fixed-node error inherent in DMC techniques.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

A Quantum Monte Carlo Approach to the Adiabatic Connection Method

We present a new method for realizing the adiabatic connection approach in density functional theory, which is based on combining accurate variational quantum Monte Carlo calculations with a constrained optimization of the ground state many-body wavefunction for different values of the Coulomb coupling constant. We use the method to study an electron gas in the presence of a cosine-wave potential. For this system we present results for the exchange-correlation hole and exchange-correlation energy density, and compare our findings with those from the local density approximation and generalized gradient approximation.

0 citations0 reviewsNo code linkedOpen access