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

N. D. Drummond

N. D. Drummond 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.mes-hallcond-mat.quant-gascond-mat.str-elphysics.comp-phcond-mat.otherastro-ph.SRcond-mat.stat-mechphysics.atm-clusquant-ph

Catalog footprint

What is connected

29works
10topics
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

29 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
preprint2016arXiv

Quasiparticle and excitonic gaps of one-dimensional carbon chains

We report diffusion quantum Monte Carlo (DMC) calculations of the quasiparticle and excitonic gaps of hydrogen-terminated oligoynes and polyyne. The electronic gaps are found to be very sensitive to the atomic structure in these systems. We have therefore optimised the geometry of polyyne by directly minimising the DMC energy with respect to the lattice constant and the Peierls-induced carbon-carbon bond-length alternation. We find the bond-length alternation of polyyne to be 0.136(2) Å and the excitonic and quasiparticle gaps to be 3.30(7) and 3.4(1) eV, respectively. The DMC zone-centre longitudinal optical phonon frequency of polyyne is 2084(5) cm$^{-1}$, which is consistent with Raman spectroscopic measurements for large oligoynes.

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

Quantum Monte Carlo Calculation of the Binding Energy of Bilayer Graphene

We report diffusion quantum Monte Carlo calculations of the interlayer binding energy of bilayer graphene. We find the binding energies of the AA- and AB-stacked structures at the equilibrium separation to be 11.5(9) and 17.7(9) meV/atom, respectively. The out-of-plane zone-center optical phonon frequency predicted by our binding-energy curve is consistent with available experimental results. As well as assisting the modeling of interactions between graphene layers, our results will facilitate the development of van der Waals exchange-correlation functionals for density functional theory calculations.

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

Electrons and phonons in single layers of hexagonal indium chalcogenides from ab initio calculations

We use density functional theory to calculate the electronic band structures, cohesive energies, phonon dispersions, and optical absorption spectra of two-dimensional In$_2$X$_2$ crystals, where X is S, Se, or Te. We identify two crystalline phases (alpha and beta) of monolayers of hexagonal In$_2$X$_2$, and show that they are characterized by different sets of Raman-active phonon modes. We find that these materials are indirect-band-gap semiconductors with a sombrero-shaped dispersion of holes near the valence-band edge. The latter feature results in a Lifshitz transition (a change in the Fermi-surface topology of hole-doped In$_2$X$_2$) at hole concentrations $n_{\rm S}=6.86\times 10^{13}$ cm$^{-2}$, $n_{\rm Se}=6.20\times 10^{13}$ cm$^{-2}$, and $n_{\rm Te}=2.86\times 10^{13}$ cm$^{-2}$ for X=S, Se, and Te, respectively, for alpha-In$_2$X$_2$ and $n_{\rm S}=8.32\times 10^{13}$ cm$^{-2}$, $n_{\rm Se}=6.00\times 10^{13}$ cm$^{-2}$, and $n_{\rm Te}=8.14\times 10^{13}$ cm$^{-2}$ for beta-In$_2$X$_2$.

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

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

Electrically Tunable Band Gap in Silicene

We report calculations of the electronic structure of silicene and the stability of its weakly buckled honeycomb lattice in an external electric field oriented perpendicular to the monolayer of Si atoms. We find that the electric field produces a tunable band gap in the Dirac-type electronic spectrum, the gap being suppressed by a factor of about eight by the high polarizability of the system. At low electric fields, the interplay between this tunable band gap, which is specific to electrons on a honeycomb lattice, and the Kane-Mele spin-orbit coupling induces a transition from a topological to a band insulator, whereas at much higher electric fields silicene becomes a semimetal.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Strategies for improving the efficiency of quantum Monte Carlo calculations

We describe a number of strategies for minimizing and calculating accurately the statistical uncertainty in quantum Monte Carlo calculations. We investigate the impact of the sampling algorithm on the efficiency of the variational Monte Carlo method. Our finding that a relative time-step ratio of 1:4 is optimal in DMC is in agreement with the result of [J. Vrbik and S.M.Rothstein, Intern. J. Quantum Chem. 29, 461-468 (1986)]. Finally, we discuss the removal of serial correlation from data sets by reblocking, setting out criteria for the choice of block length and quantifying the effects of the uncertainty in the estimated correlation length.

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
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, time-dependent density functional theory, and density functional theory calculations of diamondoid excitation energies and Stokes shifts

We have computed the absorption and emission energies and hence Stokes shifts of small diamondoids as a function of size using different theoretical approaches, including density functional theory and quantum Monte Carlo (QMC) calculations. The absorption spectra of these molecules were also investigated by time-dependent density functional theory (TD-DFT) and compared with experiment. We have analyzed the structural distortion and formation of a self-trapped exciton in the excited state, and we have studied the effects of these on the Stokes shift as a function of size. Compared to recent experiments, QMC overestimates the excitation energies by about 0.8(1) eV on average. Benefiting from a cancellation of errors, the optical gaps obtained in DFT calculations with the B3LYP functional are in better agreement with experiment. It is also shown that TD-B3LYP calculations can reproduce most of the features found in the experimental spectra. According to our calculations, the structures of diamondoids in the excited state show a distortion which is hardly noticeable compared to that found for methane. As the number of diamond cages is increased, the distortion mechanism abruptly changes character. We have shown that the Stokes shift is size-dependent and decreases with the number of diamond cages. The rate of decrease in the Stokes shift is on average 0.1 eV per cage for small diamondoids.

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