We present extensive ab initio path-integral Monte Carlo (PIMC) simulations of two-dimensional quantum dipole few-body systems ($2\ensuremath{\le}N\ensuremath{\le}7$) in a harmonic confinement, taking into account both Bose- and Fermi-statistics. This allows us to study the nonclassical rotational inertia, which can lead to a negative superfluid fraction in the case of fermions [Phys. Rev. Lett. 112, 235301 (2014)]. Moreover, we study in detail the structural characteristics of such systems and are able to clearly resolve the impact of quantum statistics on density profiles and the respective shell structure. Further, we present results for a more advanced center-two-particle correlation function [Phys. Rev. E 91, 043104 (2015)], which allows detection of differences between Fermi and Bose systems even when they are almost absent in other observables like the density. Overall, we find that bosonic systems sensitively react to even small values of the dipole-dipole coupling strength, whereas such a weak interaction is effectively masked for fermions by the Pauli exclusion principle. In addition, the abnormal superfluid fraction for fermions is not reflected by the structural properties of the system, which are equal to the bosonic case even though the moments of inertia diverge from each other. Lastly, we have explored the possibility of fermionic PIMC simulations of quantum dipole systems despite the notorious fermion sign problem, which can be further extended in future investigations in this field.
We present extensive ab initio path-integral Monte Carlo (PIMC) simulations of two-dimensional quantum dipole few-body systems (2\ensuremath≤N\ensuremath≤7) in a harmonic confinement, taking into account both Bose- and Fermi-statistics. This allows us to study the nonclassical rotational inertia, which can lead to a negative superfluid fraction in the case of fermions [Phys. Rev. Lett. 112, 235301 (2014)]. Moreover, we study in detail the structural characteristics of such systems and are able to clearly resolve the impact of quantum statistics on density profiles and the respective shell structure. Further, we present results for a more advanced center-two-particle correlation function [Phys. Rev. E 91, 043104 (2015)], which allows detection of differences between Fermi and Bose systems even when they are almost absent in other observables like the density. Overall, we find that bosonic systems sensitively react to even small values of the dipole-dipole coupling strength, whereas such a weak interaction is effectively masked for fermions by the Pauli exclusion principle. In addition, the abnormal superfluid fraction for fermions is not reflected by the structural properties of the system, which are equal to the bosonic case even though the moments of inertia diverge from each other. Lastly, we have explored the possibility of fermionic PIMC simulations of quantum dipole systems despite the notorious fermion sign problem, which can be further extended in future investigations in this field.
The ab initio path integral Monte Carlo (PIMC) method is one of the most successful methods in statistical physics, quantum chemistry and related fields, but its application to quantum degenerate Fermi systems is severely hampered by an exponential computational bottleneck: the notorious fermion sign problem. Very recently, Xiong and Xiong [J. Chem. Phys. 157, 094112 (2022)] have suggested to partially circumvent the sign problem by carrying out simulations of fictitious systems guided by an interpolating continuous variable ξ ∈ [-1, 1], with the physical Fermi- and Bose-statistics corresponding to ξ = -1 and ξ = 1. It has been proposed that information about the fermionic limit might be obtained by calculations within the bosonic sector ξ > 0 combined with an extrapolation throughout the fermionic sector ξ < 0, essentially bypassing the sign problem. Here, we show how the inclusion of the artificial parameter ξ can be interpreted as an effective penalty on the formation of permutation cycles in the PIMC simulation. We demonstrate that the proposed extrapolation method breaks down for moderate to high quantum degeneracy. Instead, the method constitutes a valuable tool for the description of large Fermi-systems of weak quantum degeneracy. This is demonstrated for electrons in a 2D harmonic trap and for the uniform electron gas (UEG), where we find excellent agreement (∼0.5%) with exact configuration PIMC results in the high-density regime while attaining a speed-up exceeding 11 orders of magnitude. Finally, we extend the idea beyond the energy and analyze the radial density distribution (2D trap), as well as the static structure factor and imaginary-time density-density correlation function (UEG).
We present a Quantum Monte Carlo study of the ground state properties of a two dimensional system of Bose particle with dipole moment in a harmonic trap. The direction of the dipoles is assumed to be fixed by an external field. We study how the system behaves when the direction of the dipole moments is changed. Our analysis is made using Path Integral Monte Carlo simulations in the Grand Canonical ensemble, using the Worm Algorithm. We study systems of increasing size, focusing on the spatial distribution of the atoms inside the trap and on the superfluid fraction of the sample. We find that vertical dipoles form a crystal, characterized by a shell structure, while tilting the dipoles has the effect of inducing a striped structure. We find that similar behaviour can be found in the analogue classical system. Vertical dipoles display nonclassical inertia while a strong tilting decreases superfluidity.
The ab initio thermodynamic simulation of correlated Fermi systems is of central importance for many applications, such as warm dense matter, electrons in quantum dots, and ultracold atoms. Unfortunately, path integral Monte Carlo (PIMC) simulations of fermions are severely restricted by the notorious fermion sign problem (FSP). In this paper, we present a hands-on discussion of the FSP and investigate in detail its manifestation with respect to temperature, system size, interaction-strength and -type, and the dimensionality of the system. Moreover, we analyze the probability distribution of fermionic expectation values, which can be non-Gaussian and fat-tailed when the FSP is severe. As a practical application, we consider electrons and dipolar atoms in a harmonic confinement, and the uniform electron gas in the warm dense matter regime. In addition, we provide extensive PIMC data, which can be used as a reference for the development of new methods and as a benchmark for approximations.
Bose-Einstein condensation has been experimentally found to take place in finite trapped systems when one of the confining frequencies is increased until the gas becomes effectively two-dimensional (2D). We confirm the plausibility of this result by performing path-integral Monte Carlo (PIMC) simulations of trapped Bose gases of increasing anisotropy and comparing them to the predictions of finite-temperature many-body theory. PIMC simulations provide an essentially exact description of these systems; they yield the density profile directly and provide two different estimates for the condensate fraction. For the ideal gas, we find that the PIMC column density of the squeezed gas corresponds quite accurately to that of the exact analytic solution and, moreover, is well mimicked by the density of a 2D gas at the same temperature; the two estimates for the condensate fraction bracket the exact result. For the interacting case, we find 2D Hartree-Fock solutions whose density profiles coincide quite well with the PIMC column densities and whose predictions for the condensate fraction are again bracketed by the PIMC estimates.
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.
Path integral Monte Carlo approach to the structural properties and collective excitations of liquid 3He without fixed nodes
Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic [Formula: see text] exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic [Formula: see text] is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid [Formula: see text] without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor [Formula: see text], the momentum distribution [Formula: see text], and the static density response function [Formula: see text]. In addition, the dynamic structure factor [Formula: see text] is rigorously reconstructed from imaginary-time PIMC data. From simulations of [Formula: see text], we derived the familiar phonon-maxon-roton dispersion function that is well-known for [Formula: see text] and has been reported previously for two-dimensional [Formula: see text] films (Nature 483:576-579 (2012)). The comparison of our new results for both [Formula: see text] and [Formula: see text] with neutron scattering measurements reveals an excellent agreement between theory and experiment.