Accurate Treatments of Electronic Correlation: Phase Transitions in an Idealized One-Dimensional Ferroelectric and Modelling Experimental Quantum Dots

Timothy James Wilkens

Published Dec 1, 2001• 1 citation

Abstract

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 modelled 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 modelling such devices using arbitrary functions represented on grids using splines.

Related papers

QMC with a Stochastic Poisson Solver: An application to realistic models of quantum dots

Dyutiman Das, Jeongnim Kim, Richard M. Martin

Bulletin of the American Physical SocietyMar 21, 20050 citations

Accurate Treatments of Electronic Correlation: Phase Transitions in an Idealized 1D Ferroelectric and Modelling Experimental Quantum Dots

Timothy James Wilkens

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

Jan 1, 19990 citations

Quantum Monte Carlo study of the phase diagram of strongly correlated two-dimensional electron systems

Hideo Aoki, 青木秀夫

Jan 1, 19950 citations

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, 19991 citation

Numerical studies on critical states in two dimensions and quantum phase transitions of correlated electrons in low dimensions

Yasuhiro Hatsugai

Jan 1, 19980 citations

First principles quantum Monte Carlo study of correlated electronic systems

Huihuo Zheng

The many-body correlation between electrons is the origin of many fascinating phenomena in condensed matter systems, such as high temperature superconductivity, superfluidity, fractional quantum Hall effect, and Mott insulator. Strongly correlated systems have been an important subject of condensed matter physics for several decades, especially after the discovery of high temperature cuprate superconductors. In this thesis, we apply first principles quantum Monte Carlo (QMC) method to several representative systems to study the electron correlations in transition metal oxides (vanadium dioxide) and low dimensional electronic systems (graphene and graphene-like two dimensional systems). Vanadium dioxide (VO2) is a paradigmatic example of a strongly correlated system that undergoes a metal-insulator transition at a structural phase transition. To date, this transition has necessitated significant post-hoc adjustments to theory in order to be described properly. We apply first principles quantum Monte Carlo (QMC) to study the structural dependence of the properties of VO2. Using this technique, we simulate the interactions between electrons explicitly, which allows for the metal-insulator transition to naturally emerge, importantly without ad-hoc adjustments. The QMC calculations show that the structural transition directly causes the metal-insulator transition and a change in the coupling of vanadium spins. This change in the spin coupling results in a prediction of a momentum-independent magnetic excitation in the insulating state. While two-body correlations are important to set the stage for this transition, they do not change significantly when VO2 becomes an insulator. These results show that it is now possible to account for electron correlations in a quantitatively accurate way that is also specific to materials. Electron correlation in graphene is unique because of the interplay of the Dirac cone dispersion of pi electrons with long range Coulomb interaction. The random phase approximation predicts no metallic screening at long distances and low energies because of the zero density of states at Fermi level. It is thus interesting to see how screening takes place in graphene at different length scales. We addressed this problem by computing the structure factor S(q) and S(q, ω) of freestanding graphene using ab initio fixed-node diffusion Monte Carlo and the random phase approximation. The X-ray measured structure factor is reproduced very accurately using both techniques, provided that sigma-bonding electrons are included in the simulations. Strong dielectric screening from sigma electrons are observed, which redshifts the pi plasmons resonance frequency at long distance and reduces the effective interactions between pi electrons at short distance. The short distance screening makes suspended graphene a weakly correlated semimetal which otherwise would be an insulator. The third piece of works is dedicated to studying the low energy excitation of many-body systems using extended Koopmans' theorem (EKT). The EKT provides a straight forward way to compute charge excitation spectra, such as ionization potentials, electron affinities from any level of theory. We implemented the EKT within the QMC framework, and performed systematic benchmark studies of ionization potentials of the second- and third-row atoms, and closed- and open-shell molecules. We also applied it to compute the quasiparticle band structure of solids (graphene). For complex correlated systems, identifying relevant low energy physics degrees of freedom is extremely important to understanding the system's collective behavior at different length scales. In this sense, bridging the realistic systems to lower energy effective lattice models that involve fewer but important degrees of freedom is significant to understanding correlated systems. We have formulated three ab initio density matrix based downfolding (AIDMD) methods to downfold the ab initio systems into effective lattice models. We have demonstrated the successfulness of these methods by applying them to molecules (H2) and periodic systems (hydrogen chain and graphene).

Jul 12, 20161 citation