Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Physics Condensed Matter

A quantum Monte Carlo method for non-parabolic electron bands in semiconductor heterostructures

John B. Shumway

Published Apr 15, 2005• 4 citations
Publisher

Abstract

We present a many-body computational technique for simulating interacting electrons in non-parabolic semiconductor bands. The technique uses an imaginary time propagator for a non-parabolic electron band that is described by an energy dependent effective mass. This derivation exploits a mathematical analogy between the kinetic energy with effective mass corrected to first order in energy and the relativistic kinetic energy. The propagator can be used in ground state and finite temperature quantum Monte Carlo (QMC) algorithms. We give a demonstration of this path integral QMC technique applied to interacting electrons in a self-assembled InGaAs quantum dot.

DOI
10.1088/0953-8984/17/17/004

Related papers

Quantum Monte Carlo simulations of electrons and holes

John B. Shumway

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.

•Jan 1, 1999•1 citation

Quantum Monte Carlo using a Stochastic Poisson Solver

Debashish Das, Richard M. Martin, Malvin H. Kalos

Quantum Monte Carlo (QMC) is an extremely powerful method to treat many-body systems. Usually quantum Monte Carlo has been applied in cases where the interaction potential has a simple analytic form, like the 1/r Coulomb potential. However, in a complicated environment as in a semiconductor heterostructure, the evaluation of the interaction itself becomes a non-trivial problem. Obtaining the potential from any grid-based finite-difference method, for every walker and every step is unfeasible. We demonstrate an alternative approach of solving the Poisson equation by a classical Monte Carlo within the overall quantum Monte Carlo scheme. We have developed a modified ''Walk On Spheres'' algorithm using Green's function techniques, which can efficiently account for the interaction energy of walker configurations, typical of quantum Monte Carlo algorithms. This stochastically obtained potential can be easily incorporated within popular quantum Monte Carlo techniques like variational Monte Carlo (VMC) or diffusion Monte Carlo (DMC). We demonstrate the validity of this method by studying a simple problem, the polarization of a helium atom in the electric field of an infinite capacitor.

Physical Review E•May 6, 2005•0 citations

Quantum Monte Carlo method using a stochastic Poisson solver

Dyutiman Das, Richard M. Martin, Malvin H. Kalos

The quantum Monte Carlo (QMC) technique is an extremely powerful method to treat many-body systems. Usually the quantum Monte Carlo method has been applied in cases where the interaction potential has a simple analytic form, like the 1/r Coulomb potential. However, in a complicated environment as in a semiconductor heterostructure, the evaluation of the interaction itself becomes a nontrivial problem. Obtaining the potential from any grid-based finite-difference method for every walker and every step is infeasible. We demonstrate an alternative approach of solving the Poisson equation by a classical Monte Carlo calculation within the overall quantum Monte Carlo scheme. We have developed a modified "walk on spheres" algorithm using Green's function techniques, which can efficiently account for the interaction energy of walker configurations, typical of quantum Monte Carlo algorithms. This stochastically obtained potential can be easily incorporated with variational, diffusion, and other Monte Carlo techniques. We demonstrate the validity of this method by studying a simple problem, the polarization of a helium atom in the electric field of an infinite capacitor.

Physical Review E•Apr 18, 2006•0 citations

Green’s function Monte Carlo method with exact imaginary-time propagation

Kevin E. Schmidt, Parhat Niyaz, Amanda L. Vaught +1 more

We present a general formulation of the Green's function Monte Carlo method in imaginary-time quantum Monte Carlo which employs exact propagators. This algorithm has no time-step errors and is obtained by minimal modifications of the time-independent Green's function Monte Carlo method. We describe how the method can be applied to the many-body Schrödinger equation, lattice Hamiltonians, and simple field theories. Our modification of the Green's function Monte Carlo algorithm is applied to the ground state of liquid 4He. We calculate the zero-temperature imaginary-time diffusion constant and relate that to the effective mass of a mass-four "impurity" atom in liquid 4He.

Physical Review E•Jan 1, 2005•44 citations

Continuous-Time Quantum Monte Carlo Algorithm for the Lattice Polaron

P. E. Kornilovitch

An efficient continuous-time path-integral quantum Monte Carlo algorithm for the lattice polaron is presented. It is based on Feynman's integration of phonons and subsequent simulation of the resulting single-particle self-interacting system. The method is free from the finite-size and finite-time-step errors and works in any dimensionality and for any range of electron-phonon interaction. The ground-state energy and effective mass of the polaron are calculated for several models. The polaron spectrum can be measured directly by Monte Carlo, which is of general interest.

arXiv:cond-mat/9808155•Physical Review Letters•Dec 14, 1998•139 citations

Electronic Energy Functionals: Levy-Lieb principle within the Ground State Path Integral Quantum Monte Carlo

Luigi Delle Site, Luca M. Ghiringhelli, David M. Ceperley

We propose a theoretical/computational protocol based on the use of the Ground State (GS) Path Integral (PI) Quantum Monte Carlo (QMC) for the calculation of the kinetic and Coulomb energy density for a system of $N$ interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the $N-1$-conditional probability density within the framework of the Levy-Lieb constrained search principle. The consequences for the development of energy functionals within the context of Density Functional Theory (DFT) are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the $N-1$-conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an {\it on-the-fly} determination of the energy functional for any system.

arXiv:1205.6503•May 29, 2012•0 citations
We present a many-body computational technique for simulating interacting electrons in non-parabolic semiconductor bands. The technique uses an imaginary time propagator for a non-parabolic electron band that is described by an energy dependent effective mass. This derivation exploits a mathematical analogy between the kinetic energy with effective mass corrected to first order in energy and the relativistic kinetic energy. The propagator can be used in ground state and finite temperature quantum Monte Carlo (QMC) algorithms. We give a demonstration of this path integral QMC technique applied to interacting electrons in a self-assembled InGaAs quantum dot.
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.
Quantum Monte Carlo (QMC) is an extremely powerful method to treat many-body systems. Usually quantum Monte Carlo has been applied in cases where the interaction potential has a simple analytic form, like the 1/r Coulomb potential. However, in a complicated environment as in a semiconductor heterostructure, the evaluation of the interaction itself becomes a non-trivial problem. Obtaining the potential from any grid-based finite-difference method, for every walker and every step is unfeasible. We demonstrate an alternative approach of solving the Poisson equation by a classical Monte Carlo within the overall quantum Monte Carlo scheme. We have developed a modified ''Walk On Spheres'' algorithm using Green's function techniques, which can efficiently account for the interaction energy of walker configurations, typical of quantum Monte Carlo algorithms. This stochastically obtained potential can be easily incorporated within popular quantum Monte Carlo techniques like variational Monte Carlo (VMC) or diffusion Monte Carlo (DMC). We demonstrate the validity of this method by studying a simple problem, the polarization of a helium atom in the electric field of an infinite capacitor.
The quantum Monte Carlo (QMC) technique is an extremely powerful method to treat many-body systems. Usually the quantum Monte Carlo method has been applied in cases where the interaction potential has a simple analytic form, like the 1/r Coulomb potential. However, in a complicated environment as in a semiconductor heterostructure, the evaluation of the interaction itself becomes a nontrivial problem. Obtaining the potential from any grid-based finite-difference method for every walker and every step is infeasible. We demonstrate an alternative approach of solving the Poisson equation by a classical Monte Carlo calculation within the overall quantum Monte Carlo scheme. We have developed a modified "walk on spheres" algorithm using Green's function techniques, which can efficiently account for the interaction energy of walker configurations, typical of quantum Monte Carlo algorithms. This stochastically obtained potential can be easily incorporated with variational, diffusion, and other Monte Carlo techniques. We demonstrate the validity of this method by studying a simple problem, the polarization of a helium atom in the electric field of an infinite capacitor.
We present a general formulation of the Green's function Monte Carlo method in imaginary-time quantum Monte Carlo which employs exact propagators. This algorithm has no time-step errors and is obtained by minimal modifications of the time-independent Green's function Monte Carlo method. We describe how the method can be applied to the many-body Schrödinger equation, lattice Hamiltonians, and simple field theories. Our modification of the Green's function Monte Carlo algorithm is applied to the ground state of liquid 4He. We calculate the zero-temperature imaginary-time diffusion constant and relate that to the effective mass of a mass-four "impurity" atom in liquid 4He.
An efficient continuous-time path-integral quantum Monte Carlo algorithm for the lattice polaron is presented. It is based on Feynman's integration of phonons and subsequent simulation of the resulting single-particle self-interacting system. The method is free from the finite-size and finite-time-step errors and works in any dimensionality and for any range of electron-phonon interaction. The ground-state energy and effective mass of the polaron are calculated for several models. The polaron spectrum can be measured directly by Monte Carlo, which is of general interest.
We propose a theoretical/computational protocol based on the use of the Ground State (GS) Path Integral (PI) Quantum Monte Carlo (QMC) for the calculation of the kinetic and Coulomb energy density for a system of NNN interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the N−1N-1N−1-conditional probability density within the framework of the Levy-Lieb constrained search principle. The consequences for the development of energy functionals within the context of Density Functional Theory (DFT) are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the N−1N-1N−1-conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an {\it on-the-fly} determination of the energy functional for any system.