Electron and hole excitations in semiconductors may be approximated as particles with effective masses which interact via Coulomb potentials. We study systems of electrons and holes with quantum Monte Carlo (QMC) techniques, covering three related areas: (1) elastic scattering of excitons, (2) thermodynamics of electron-hole plasmas, and (3) electrons confined in a quantum dot. Excitons are bound states of an electron and a hole, and obey Bose statistics. A low density exciton gas is an experimentally realizable dilute Bose gas. The scattering length a8 of a dilute Bose gas determines its properties, but is difficult to calculate. We present an essentially exact QMC treatment of exciton-exciton scattering, and find scattering lengths for different spin orientations of the excitons. At some mass ratios mh/me the scattering lengths diverge in conjunction with the appearance of biexciton vibrational states, an effect not found by earlier perturbative treatments. Path integral Monte Carlo (PIMC) is used to model the thermodynamics of the electron-hole plasma. Our primary interest is the Bose condensation of an excitonic gas. At low density and low temperature the spin-unpolarized system forms biexcitons. Since we are interested in Bose condensation, we study a spin-polarized system, which has no biexcitons. Restricted paths are used to handle the Fermion sign problem. With an appropriate choice of paired nodes for the restricted path approximation we find an excitonic Bose condensate. The energy of the low temperature, low density exciton gas determined from PIMC agrees well with the theory of dilute Bose gases, in which our previously calculated scattering length is used to model the exciton-exciton interactions. At higher densities the excitons are less well defined and the transition changes character. Finally, we study electrons in the inhomogeneous environment of a self-assembled InAs-GaAs quantum dot. We combine our ground state QMC treatment with another method, density functional theory (DFT) within the local spin density approximation (LSDA). Our comparison shows that LSDA is acceptable for treating interactions in the case considered, but recommend further tests for application of LSDA to larger dots or coupled dot systems.
Electron and hole excitations in semiconductors may be approximated as particles with effective masses which interact via Coulomb potentials. We study systems of electrons and holes with quantum Monte Carlo (QMC) techniques, covering three related areas: (1) elastic scattering of excitons, (2) thermodynamics of electron-hole plasmas, and (3) electrons confined in a quantum dot. Excitons are bound states of an electron and a hole, and obey Bose statistics. A low density exciton gas is an experimentally realizable dilute Bose gas. The scattering length a8 of a dilute Bose gas determines its properties, but is difficult to calculate. We present an essentially exact QMC treatment of exciton-exciton scattering, and find scattering lengths for different spin orientations of the excitons. At some mass ratios mh/me the scattering lengths diverge in conjunction with the appearance of biexciton vibrational states, an effect not found by earlier perturbative treatments. Path integral Monte Carlo (PIMC) is used to model the thermodynamics of the electron-hole plasma. Our primary interest is the Bose condensation of an excitonic gas. At low density and low temperature the spin-unpolarized system forms biexcitons. Since we are interested in Bose condensation, we study a spin-polarized system, which has no biexcitons. Restricted paths are used to handle the Fermion sign problem. With an appropriate choice of paired nodes for the restricted path approximation we find an excitonic Bose condensate. The energy of the low temperature, low density exciton gas determined from PIMC agrees well with the theory of dilute Bose gases, in which our previously calculated scattering length is used to model the exciton-exciton interactions. At higher densities the excitons are less well defined and the transition changes character. Finally, we study electrons in the inhomogeneous environment of a self-assembled InAs-GaAs quantum dot. We combine our ground state QMC treatment with another method, density functional theory (DFT) within the local spin density approximation (LSDA). Our comparison shows that LSDA is acceptable for treating interactions in the case considered, but recommend further tests for application of LSDA to larger dots or coupled dot systems.
The work in this thesis is concerned with the application and development of quantum Monte Carlo (QMC) methods. We begin by proposing a technique to maximise the efficiency of the extrapolation of DMC results to zero time step, finding that a relative time step ratio of 1:4 is optimal. We discuss the post-processing of QMC data and the calculation of accurate error bars by reblocking, setting out criteria for the choice of block length. We then quantify the effects of uncertainty in the correlation length on estimated error bars, finding that the frequency of outliers is significantly increased for short runs. We then report QMC calculations of biexciton binding energies in bilayer systems. We have also calculated exciton-exciton interaction potentials, and radial distribution functions for electrons and holes in bound biexcitons. We find a larger region of biexciton stability than other recent work [C. Schindler and R. Zimmermann, Phys. Rev. B 78, 045313 (2008)]. We also find that individual excitons retain their identity in bound biexcitons for large layer separations. Finally, we give details of a QMC study of the one-dimensional homogeneous electron gas (1D HEG). We present calculations of the energy, pair correlation function, static structure factor (SSF), and momentum density (MD) for the 1D HEG. We observe 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 exponents of the ferromagnetic infinitely-thin and harmonic wires.
The bilayer Hubbard model with electron-hole doping is an ideal platform to study excitonic orders due to suppressed recombination via spatial separation of electrons and holes. However, suffering from the sign problem, previous quantum Monte Carlo studies could not arrive at an unequivocal conclusion regarding the presence of phases with clear signatures of excitonic condensation in bilayer Hubbard models. Here, we develop a determinant quantum Monte Carlo algorithm for the bilayer Hubbard model that is sign-problem-free for equal and opposite doping in the two layers and study excitonic order and charge and spin density modulations as a function of chemical potential difference between the two layers, on-site Coulomb repulsion, and interlayer interaction. In the intermediate coupling regime and in proximity to the SU(4)-symmetric point, we find a biexcitonic condensate phase at finite electron-hole doping, as well as a competing (π,π) charge density wave state. We extract the Berezinskii-Kosterlitz-Thouless transition temperature from superfluid density and a finite-size scaling analysis of the correlation functions and explain our results in terms of an effective biexcitonic hard-core boson model.
We review the path integral method wherein quantum systems are mapped with Feynman's path integrals onto a classical system of "ring-polymers" and then simulated with the Monte Carlo technique. Bose or Fermi statistics correspond to possible "cross-linking" of polymers. As proposed by Feynman, superfluidity and Bose condensation result from macroscopic exchange of bosons. To map fermions onto a positive probability distribution, one must restrict the paths to lie in regions where the fermion density matrix is positive. We discuss a recent application to the two-component electron-hole plasma. At low temperature excitons and bi-excitons form. We have used nodal surfaces incorporating paired fermions and see evidence of a Bose condensation in the energy, specific heat and superfluid density. In the restricted path integral picture, pairing appears as intertwined electron-hole paths. Bose condensation occurs when these intertwined paths wind around the periodic boundaries.
The properties of the macroscopic world around us, and of which we are a part, are largely determined by the low energy, collective behaviour of many interacting particles, including the nuclei and, especially, the electrons present. Although the fundamental laws governing the behaviour of these many-body systems are believed to be known in principle, the practical solution of the equations of quantum mechanics remains a challenging area of research. This thesis is concerned with the application of quantum Monte Carlo methods to two model systems: the spin-polarised homogeneous electron gas, and a hole-doped electron gas. Electronic structure theory is briefly reviewed before discussing in more detail the quantum Monte Carlo methods used in this thesis. A study of the three-dimensional spin-polarised homogeneous electron gas (HEG) is then reported, where the relatively new technique of twist averaging is investigated in detail and accurate energies and pair correlation functions are obtained over densities rs=0.5–20 a.u. and the full range of spin-polarisation, allowing comparison with the Perdew-Zunger interpolation scheme used in local spin density approximation exchange-correlation functionals. Following this, an impurity is added to the electron gas in the form of a positively charged hole, and the interaction is studied. Relaxation energies, pair correlation functions and momentum densities are reported. Trion formation is observed over a range of carrier densities and electron-hole mass ratios in agreement with experiment. Isolated trions are also studied, where the diffusion Monte Carlo method is exact. Methodological innovations developed while carrying out this work are discussed, including a variance reduction technique for twist-averaged calculations and a new trial wave function for impurity-in-HEG calculations.
"Quantum Monte Carlo (QMC) is one of the most powerful methods for solving many body problems in Physics and Chemistry. MC sampling allows one to measure physical properties systems with very complicated wave functions that incorporate correlation between the electrons. Additionally, QMC has O(N3) scaling verses exponential scaling for methods such as Conguration Interaction and “exact diagonalization"", and it’s perfect parallelism allows one to speed up computation with little effort on large parallel computers. Our work has focused on two systems. The first is a 1D ionic Hubbard model with a complex phase diagram that we find to exhibit a ferroeletric \bond ordered"" phase. A new many-body berry's phase formulation of the electronic polarization and localization is used in conjunction with QMC to analyze whether a bond order phase exists between the customary band and Mott insulating phases. The bond order and bond order correlation functions are measured to illustrate the stability of this intermediate state, and phase swapping between two degenerate states of bond order (§B) is observed. We propose that the Mott state it is unstable to bond ordering for any degree of ionicity. The second system we have studied is a set of idealized (harmonic) and realistic quantum dots. Two idealized dots, spherical and squashed, quasi-2D structures, are analyzed. In the first QMC is contrasted with a large set of commonly used electronic structure methods including MSFT, KLI, LSDA, and CCSD, and in the second we scrutinize whether SDW states truly exist as predicted by LSDA. In conjunction with the Computational Electronics Group at the Beckman Institute, we have performed QMC calculations on realistically modeled quantum dots in the effective mass approximation. The potentials of these dots are determined using the FEM method to solve Poisson's equation and the effective mass is allowed to spatially vary in our QMC simulations. We find that in arrays of QD's, LSDA may be insufficiently accurate to determine the proper spin states of the electrons. Additionally, we illustrate that QMC has a bright future ahead of it in modeling such devices using arbitrary functions represented on grids using splines."
We evaluate the zero temperature properties and Bose-Einstein condensation (BEC) of hard sphere bosons in a harmonic trap using Quantum Monte Carlo methods. The chief purpose of this work is to go beyond the dilute limit, to test the limits of the Gross-Pitaevskii equation and related mean field approximations and to explore the role of interactions in determining the zero temperature properties of the trapped Bose gas. The density is increased by adjusting both the number of trapped bosons, N , and the s-wave scattering length, a, to study systems from the highly dilute region corresponding to early experiments in trapped BEC [2, 3] up to liquid He densities and beyond. Rather than relying on the perturbative Bogoliubov approximation (which assumes a large condensate fraction) to describe the condensate, we use the one-body density matrix (OBDM) formulation of BEC so that the properties of the condensate for systems with arbitrarily large interactions may be studied. In this formulation of BEC, condensate properties are obtained by diagonalizing the OBDM and obtaining the corresponding single particle “natural orbitals” and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single particle orbital(s) with macroscopic occupation (N0 >> 1). Within this framework, we calculate the effects of interactions and increased density on the ground state energy, the density profile, the momentum distribution, the condensate fraction and condensate “wave-function” and several other properties.