Described in this dissertation is the use of quantum Monte Carlo methods to study two ideas in quantum many-body problems: superfluidity and entanglement. Density matrices are presented a central tool in the analysis, as are discussed in the review of path integral Monte Carlo (PIMC) and variational Monte Carlo (VMC) methods. PIMC is used to model a one-dimensional system of fermionic lithium atoms according to existing experiments, including a realistic temperature. New estimators of the pair momentum distribution are implemented, yielding in a clear in-situ signature of a pairing mechanism (dubbed FFLO after its first proposers) which implies a microscopic phase fluctuation in space between a normal fluid and a superfluid. VMC is used to model homonuclear diatomic molecules of period-2 elements. The degree of entanglement and the responsible electronic configurations in real space are quantified in terms of the entanglement spectra. Calculating the reduced denstity matrix as an intermediate step reveals a novel way of understanding chemical bonds, as exemplified by Be2 and C2 . Possible implications of these results in integrable many-body models and in quantum chemistry are discussed, as well as direction for future investigation.
Described in this dissertation is the use of quantum Monte Carlo methods to study two ideas in quantum many-body problems: superfluidity and entanglement. Density matrices are presented a central tool in the analysis, as are discussed in the review of path integral Monte Carlo (PIMC) and variational Monte Carlo (VMC) methods. PIMC is used to model a one-dimensional system of fermionic lithium atoms according to existing experiments, including a realistic temperature. New estimators of the pair momentum distribution are implemented, yielding in a clear in-situ signature of a pairing mechanism (dubbed FFLO after its first proposers) which implies a microscopic phase fluctuation in space between a normal fluid and a superfluid. VMC is used to model homonuclear diatomic molecules of period-2 elements. The degree of entanglement and the responsible electronic configurations in real space are quantified in terms of the entanglement spectra. Calculating the reduced denstity matrix as an intermediate step reveals a novel way of understanding chemical bonds, as exemplified by Be2 and C2 . Possible implications of these results in integrable many-body models and in quantum chemistry are discussed, as well as direction for future investigation.
This thesis is divided into two parts. In the first part we study quantum phases in bosonic systems in an optical lattice by using numerical unbiased Quantum Monte Carlo (QMC) method. Ultracold atoms experiments in optical lattices provide a possibility to simulate the materials of condensed matter in a clean and well-controled way. We focus on two components bosonic systems where exotic phases such as pair-superfluid and pair-supersolid appear. In the second part we provide numerical methods based on QMC simulations to study the entanglement properties of strong correlated systems. By employing the replica trick, we reconstruct the entanglement spectrum from the Renyi entropies. Furthermore, we study the trace of the power of the partial transposed reduced density matrix, which is related to the entanglement measurement negativity for mixed states. In the following we briefly summarize each chapter. In Chapter 1 we very briefly introduce the cold atom systems and the connection with QMC simulations. In Chapter 2 we review the QMC method based on the path integral (world line) representation. We focus on the directed worm algorithm which is an efficient algorithm with global updtae for bosonic and spin systems. All the following methods in this thesis – the method of two component bosonic systems, the reconstruction of entanglement spectrum and the measuring of partial transposed quantities, are based on directed worm algorithm introduced in this chapter. In Chapter 3 we discuss two component bosons in a square lattice. We show that the interspecies attraction and nearest-neighbor intraspecies repulsion result in the pair-supersolid phase, where a diagonal solid order coexists with an off-diagonal pair-superfluid order. The quantum and thermal transitions out of the pair-supersolid phase are characterized. It is found that there is a direct first-order transition from the pair-supersolid phase to the double-superfluid phase without an intermediate region. Furthermore, the melting of the pair-supersolid occurs in two steps. Upon heating, first the pair-superfluid is destroyed via a Kosterlitz-Thouless transition, then the solid order melts via an Ising transition. In Chapter 4 we represent a new method to reconstruct a subset of the entanglement spectrum of quantum many body systems by QMC, where the method can in principle be applied to two or higher dimension. The approach builds on the replica trick to evaluate particle number resolved traces of the first n of powers of a reduced density matrix. From this information we reconstruct first n entanglement spectrum levels using a polynomial root solver. We illustrate the power and limitations of the method by an application to the extended Bose-Hubbard model in one dimension where we are able to resolve the quasidegeneracy of the entanglement spectrum in the Haldane-insulator phase. In general, the method is able to reconstruct the largest few eigenvalues in each symmetry sector. In Chapter 5 we devise a Quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. In particular, we study several scale invariant combinations of the moments for the 1D hard-core boson model. For two adjacent intervals unusual finite size corrections are present, showing parity effects that oscillate with a filling dependent period. For large chains we find perfect agreement with conformal field theory (CFT) calculations. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. The finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderate large sizes.
Variational and diffusion Monte Carlo (VMC and DMC) are widely used techniques in the Monte Carlo study of many-body quantum systems. Through efficient use of ran- dom sampling to solve many-dimensional integrals, they can be used to reliably solve the Schrödinger equation. Within the wider ecosystems of electronic structure calculation, quantum chemistry, and quantum simulation, quantum Monte Carlo (QMC) methods are used in tandem with density functional theory (DFT) and other computational methods to deliver insight into a vast array of condensed matter systems. Of particular interest to this work are systems in the intermediate correlation regime, where methods such as simple band theory and the Hubbard model break down. This thesis proposes and tests an explicit three-body backflow term to improve the accuracy of VMC and DMC calculations by incorporating higher-order correlation in order to reduce fixed-node er- ror in DMC, finding statistically significant improvements in ground state energies and improved variance in a water molecule. Then, QMC is used to calculate accurate charge densities for weak and intermediately correlated materials, to test commonly used dens- ity functional approximations in a controlled environment, providing insight into where these approximations perform well, and where they fail. Finally, DMC is used to generate ground state energy data for a dipolar Fermi gas of dressed Rydberg atoms interacting in the intermediate correlation regime. This data is parametrised in terms of interaction strength, for use in generating exchange-correlation functionals, to enable the DFT study of inhomogeneous systems of dressed Rydberg atoms.
Several microscopic theories point out that Bose-Einstein condensation (BEC), i.e., a macroscopic occupation of the lowest energy single particle state in many-boson systems, may appear also in quantum fluids and solids and that it is at the origin of the phenomenon of superfluidity. Nevertheless, the connection between BEC and superfluidity is still matter of debate, since the experimental evidences indicating a non zero condensate fraction in superfluid helium are only indirect. In the theoretical study of BEC in quantum fluids and solids, perturbative approaches are useless because of the strong correlations between the atoms, arising both from the interatomic potential and from the quantum nature of the system. Microscopic Quantum Monte Carlo simulations provide a reliable description of these systems. In particular, the Path Integral Monte Carlo (PIMC) method is very suitable for this purpose. This method is able to provide exact results for the properties of the quantum system, both at zero and finite temperature, only with the definition of the Hamiltonian and of the symmetry properties of the system, giving an easy picture for superfluidity and BEC in many-boson systems. In this thesis, we apply PIMC methods to the study of several quantum fluids and solids. We describe in detail all the features of PIMC, from the sampling methods to the estimators of the physical properties. We present also the most recent techniques, such as the high-order approximations for the thermal density matrix and the worm algorithm, used in PIMC to provide reliable simulations. We study the liquid phase of condensed 4He, providing unbiased estimations of the one-body density matrix g1(r). We analyze the model for g1(r) used to fit the experimental data, highlighting its merits and its faults. In particular we see that, even if it presents some difficulties in the description of the overall behavior of g1(r), it can provide an accurate estimation of the kinetic energy K and of the condensate fraction n0 of the system. Furthermore, we show that our results for n0 as a function of the pressure are in a good agreement with the most recent experimental results. The study of the solid phase of 4He is the most significant part of this thesis. The recent observation of non classical rotational inertia (NCRI) effects in solid helium has generated big interest in the study of an eventual supersolid phase, characterized at the same time by crystalline order and superfluidity. Nevertheless, until now it has been impossible to give a theoretical model able to describe all the experimental evidences. In this work, we perform PIMC simulations of 4He at high densities, according to different microscopic configurations of the atoms. In commensurate crystals we see that BEC does not appear, our model being able to reproduce the momentum distribution obtained form neutron scattering experiments. In a crystal with vacancies, we have been able to see a transition to a superfluid phase at temperatures in agreement with experimental results if the vacancy concentration is low enough. In amorphous solids, superfluid effects are enhanced but appear at temperatures higher than the experimental estimation for the transition temperature. Finally, we study also metastable disordered configurations in molecular para-hydrogen at low temperature. The aim of this study is to investigate if a Bose liquid other than helium can display superfluidity. Choosing accurately a ¿quantum liquid¿ initial configuration and the dimensions of the simulation box, we have been able to frustrate the formation of the crystal and to calculate the temperature dependence of the superfluid density, showing a transition to a superfluid phase at temperatures close to 1 K.
Quantum Monte Carlo methods are a powerful suite of techniques for solving the quantum many-body problem. By using random numbers to stochastically sample quantum properties, QMC methods are capable of studying low-temperature quantum systems well beyond the reach of conventional deterministic techniques. QMC techniques have likewise been indispensible tools for augmenting our current knowledge of superfluidity and superconductivity. In this thesis, I present two new quantum Monte Carlo techniques, the Monte Carlo Power Method and Bose-Fermi Auxiliary-Field Quantum Monte Carlo, and apply previously developed Path Integral Monte Carlo methods to explore two new phases of quantum hard spheres and hydrogen. I lay the foundation for a subsequent description of my research by first reviewing the physics of quantum liquids in Chapter One and the mathematics behind Quantum Monte Carlo algorithms in Chapter Two. I then discuss the Monte Carlo Power Method, a stochastic way of computing the first several extremal eigenvalues of a matrix too memory-intensive to be stored and therefore diagonalized. As an illustration of the technique, I demonstrate how it can be used to determine the second eigenvalues of the transition matrices of several popular Monte Carlo algorithms. This information may be used to quantify how rapidly a Monte Carlo algorithm is converging to the equilibrium probability distribution it is sampling. I next present the Bose-Fermi Auxiliary-Field Quantum Monte Carlo algorithm. This algorithm generalizes the well-known Auxiliary-Field Quantum Monte Carlo algorithm for fermions to bosons and Bose-Fermi mixtures. Despite some shortcomings, the Bose-Fermi Auxiliary-Field Quantum Monte Carlo algorithm represents the first exact technique capable of studying Bose-Fermi mixtures of any size in any dimension. In Chapter Six, I describe a new Constant Stress Path Integral Monte Carlo algorithm for the study of quantum mechanical systems under high pressures. While the eventual hope is to apply this algorithm to the exploration of yet unidentified high-pressure, low-temperature phases of hydrogen, I employ this algorithm to determine whether or not quantum hard spheres can form a low-temperature bcc solid if exchange is not taken into account. In the final chapter of this thesis, I use Path Integral Monte Carlo once again to explore whether glassy para-hydrogen exhibits superfluidity. Physicists have long searched for ways to coax hydrogen into becoming a superfluid. I present evidence that, while glassy hydrogen does not crystallize at the temperatures at which hydrogen might become a superfluid, it nevertheless does not exhibit superfluidity. This is because the average binding energy per p-H2 molecule poses a severe barrier to exchange regardless of whether the system is crystalline. All in all, this work extends the reach of Quantum Monte Carlo methods to new systems and brings the power of existing methods to bear on new problems.
We study systems of particles with varying numbers of distinguishable species ranging from 4-species to the N-species limit, where N is the number of particles in the system itself. We examine these systems using various powerful simulation techniques that belong to the quantum Monte Carlo (QMC) family. The N-species limit is studied via Path Integral Monte Carlo (PIMC), which is a finite temperature, non-perturbative technique. At this limit, we explore the viability of two different thermal density matrices for systems that interact via hard-sphere, hard-cavity potentials. We then study the impact of finite-size effects on calculations of the energy, pressure and specific heat of the system when reaching the thermodynamic limit under periodic boundary conditions. We then move to studying a 4-species fermionic system using Variational Monte Carlo (VMC) and Diffusion Monte Carlo (DMC), which are ground state techniques. We apply these methods to the clustering problem of 4-particle and 8-particle systems. While the 4-particle, 4-species system is relatively simple to calculate as it is physically identical to a bosonic system of 4-particles, the 8-particle fermi system is more complicated. We detail the process of simulating a bound 8-particle state with respect to decay into two independent 4-particle clusters with DMC for the first time.
Several algorithms have been proposed to calculate the spatial entanglement spectrum from high order R\'enyi entropies and maximum entropy techniques. We present an alternative approach for computing the entanglement spectrum with quantum Monte Carlo for both continuum and lattice Hamiltonians. This method provides direct access to the matrix elements of the spatially reduced density matrix and we determine an estimator that can be used in variational Monte Carlo as well as other Monte Carlo methods. The algorithm is based on using a generalization of the swap operator, which can be extended to calculate a general class of density matrices that can include combinations of spin, space, particle, and even momentum degrees of freedom. We demonstrate the method by applying it to the H2 and N2 molecules and describe how the spatial entanglement spectrum encodes a covalent bond that includes all the many body correlations.