Path integral Monte Carlo and the electron gas

Brown, Ethan

Permalink

https://hdl.handle.net/2142/50560 Copy

Description

Title

Path integral Monte Carlo and the electron gas

Author(s)

Brown, Ethan

Issue Date

2014-09-16

Director of Research (if dissertation) or Advisor (if thesis)

Ceperley, David M.

Doctoral Committee Chair(s)

Hughes, Taylor L.

Committee Member(s)

• Ceperley, David M.
• Cooper, S. Lance
• DuBois, Jonathan
• Schleife, Andre

Department of Study

Physics

Discipline

Physics

Degree Granting Institution

University of Illinois at Urbana-Champaign

Degree Name

Ph.D.

Degree Level

Dissertation

Keyword(s)

• Path integral Monte Carlo
• electron gas
• quantum Monte Carlo
• electronic structure

Abstract

Path integral Monte Carlo is a proven method for accurately simulating quantum mechanical systems at finite-temperature. By stochastically sampling Feynman's path integral representation of the quantum many-body density matrix, path integral Monte Carlo includes non-perturbative effects like thermal fluctuations and particle correlations in a natural way. Over the past 30 years, path integral Monte Carlo has been successfully employed to study the low density electron gas, high-pressure hydrogen, and superfluid helium. For systems where the role of Fermi statistics is important, however, traditional path integral Monte Carlo simulations have an exponentially decreasing efficiency with decreased temperature and increased system size. In this thesis, we work towards improving this efficiency, both through approximate and exact methods, as specifically applied to the homogeneous electron gas. We begin with a brief overview of the current state of atomic simulations at finite-temperature before we delve into a pedagogical review of the path integral Monte Carlo method. We then spend some time discussing the one major issue preventing exact simulation of Fermi systems, the sign problem. Afterwards, we introduce a way to circumvent the sign problem in PIMC simulations through a fixed-node constraint. We then apply this method to the homogeneous electron gas at a large swatch of densities and temperatures in order to map out the warm-dense matter regime. The electron gas can be a representative model for a host of real systems, from simple medals to stellar interiors. However, its most common use is as input into density functional theory. To this end, we aim to build an accurate representation of the electron gas from the ground state to the classical limit and examine its use in finite-temperature density functional formulations. The latter half of this thesis focuses on possible routes beyond the fixed-node approximation. As a first step, we utilize the variational principle inherent in the path integral Monte Carlo method to optimize the nodal surface. By using a ansatz resembling a free particle density matrix, we make a unique connection between a nodal effective mass and the traditional effective mass of many-body quantum theory. We then propose and test several alternate nodal ansatzes and apply them to single atomic systems. Finally, we propose a method to tackle the sign problem head on, by leveraging the relatively simple structure of permutation space. Using this method, we find we can perform exact simulations this of the electron gas and $^3$He that were previously impossible.

Graduation Semester

2014-08

Permalink

http://hdl.handle.net/2142/50560

Copyright and License Information

Copyright 2014 Ethan Brown

Owning Collections