Path Integral Ground State Research Articles

Path-integral ground state and optical response of the charged Bose polaron

Path-integral ground state and optical response of the charged Bose polaron

Effect of neighbouring molecules on ground-state properties of many-body polar linear rotor systems

A path integral ground state approach has been used to estimate the ground-state energy and structural properties of hydrogen fluoride molecules pinned to a one-dimensional lattice. In the simulations, the molecules are assumed to be rigid, and only the continuous rotational degrees of freedom are considered. The constituents of a many-body system interact through the dipole-dipole interaction because the molecules have a permanent dipole moment. The workability of our approach has been demonstrated by estimating the ground-state energy, order parameter and nearest neighbour correlation for the systems of 2 and 3 HF molecules using quantum Monte Carlo simulations based on the path integral ground state methodology. The results agree satisfactorily with those obtained from exact Hamiltonian matrix diagonalisation. In addition, the effect of neighbours on the ground state properties has been investigated for larger systems. The converged ground-state energy per neighbour as a function of inter-nuclear separation is considered to be the equation of state for the system. The chemical potential of the rotor chains also supports the convergence test.

Properties of fermionic systems with the path-integral ground state method

We investigate strongly correlated many-body systems composed of bosons and fermions with a fully quantum treatment using the path-integral ground state method, PIGS. To account for the Fermi-Dirac statistics, we implement the fixed-node approximation into PIGS, which we then call FN-PIGS. In great detail, we discuss the pair density matrices we use to construct the full density operator in coordinate representation, a vital ingredient of the method. We consider the harmonic oscillator as a proof-of-concept and, as a platform representing quantum many-body systems, we explore helium atoms. Pure ^44He systems demonstrate most of the features of the method. Complementarily, for pure ^33He, the fixed-node approximation resolves the ubiquitous sign problem stemming from anti-symmetric wave functions. Finally, we investigate ^33He-^44He mixtures, demonstrating the method’s robustness. One of the main features of FN-PIGS is its ability to estimate any property at temperature T = 0T=0 without any additional bias apart from the FN approximation; biases from long simulations are also excluded. In particular, we calculate the correlation function of pairs of equal and opposite spins and precise values of the ^33He kinetic energy in the mixture.

A path integral ground state Monte Carlo algorithm for entanglement of lattice bosons

A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to L=256L=256 sites at unit-filling without any restrictions on site occupancy, far beyond the reach of exact diagonalization. The favorable scaling of the algorithm is demonstrated through a further measurement of the Rényi entanglement entropy at the two dimensional superfluid-insulator critical point for large system sizes, confirming the existence of the expected entanglement boundary law in the ground state. The Rényi estimator is extended to measure the symmetry resolved entanglement that is operationally accessible as a resource for experimentally relevant lattice gases with fixed total particle number.

Absence of Off-Diagonal Long-Range Order in hcp ^{4}He Dislocation Cores.

The mass transport properties along dislocation cores in hcp ^{4}He are revisited by considering two types of edge dislocations as well as a screw dislocation, using a fully correlated quantum simulation approach. Specifically, we employ the zero-temperature path-integral ground state (PIGS) method together with ergodic sampling of the permutation space to investigate the fundamental dislocation core structures and their off-diagonal long-range order properties. It is found that the Bose-Einstein condensate fraction of such defective ^{4}He systems is practically null (≤10^{-6}), just as in the bulk defect-free crystal. These results provide compelling evidence for the absence of intrinsic superfluidity in dislocation cores in hcp ^{4}He and challenge the superfluid dislocation-network interpretation of the mass-flux-experiment observations, calling for further experimental investigation.

Path-integral quantum Monte Carlo calculations of light nuclei

We describe a path-integral ground-state quantum Monte Carlo method for light nuclei in continuous space. We show how to efficiently update and sample the paths with spin-isospin dependent and spin-orbit interactions. We apply the method to the triton and alpha particle using both local chiral interactions with next-to-next-to-leading-order %(N$^2$LO) and the Argonne interactions. For operators, like the total energy, that commute with the Hamiltonian, our results agree with Green's function Monte Carlo and auxiliary field diffusion Monte Carlo calculations. For operators that do not commute with the Hamiltonian and for Euclidean response functions, the path-integral formulation allows straightforward calculation without forward walking or the increased variance typical of diffusion methods. We demonstrate this by calculating density distributions, root mean square radii, and Euclidean response functions for single-nucleon couplings.

On the accuracy and efficiency of different methods to calculate Raman vibrational shifts of parahydrogen clusters.

The Raman vibrational frequency shifts of pure parahydrogen and orthodeuterium clusters of sizes N = 4-9 are calculated using the Langevin equation path integral ground state method. The shifts are calculated using three different methods; the results obtained from each are compared to experiment and variance properties are assessed. The first method requires the direct calculation of energies from two simulations: one when the cluster is in the v = 0 vibrational state and one when the cluster has v = 1 total quantum of vibration. The shift is directly calculated from the difference in those two energies. The second method requires only a v = 0 simulation to be performed. The ground state energy is calculated as usual and the excited state energy is calculated by using the distribution of the v = 0 simulation and the ratio of the density matrices between the v = 1 state and the v = 0 state. The shift is calculated from the difference in those two energies. These first two are both exact methods. The final method is based on perturbation theory where the shift is calculated by averaging the pairwise difference potential over the pair distribution function. However, this is an approximate approach. It is found that for large enough system sizes, despite the approximations, the perturbation theory method has the strongest balance between accuracy and precision when weighing against computational cost.

Ground state chemical potential of parahydrogen clusters of size N = 21-40.

We report the ground state chemical potential of parahydrogen clusters between N = 21-40 calculated using the Langevin equation Path Integral Ground State method. There has been much debate in the past whether the chemical potential size evolution in this region is jagged (indicating magic number cluster sizes) or if it is smooth (indicating some quantum melting below 1K). We compare to previous diffusion Monte Carlo and Path Integral Ground State (PIGS) results, including very recent Variational Path Integral Molecular Dynamics (VPIMD) calculations [S. Miura, J. Chem. Phys. 148, 102333 (2018)]. We find that the ground state chemical potential is not a smooth curve and that magic number clusters are present, consistent with VPIMD and PIGS Monte Carlo results.

Path integral simulations of confined parahydrogen molecules within clathrate hydrates: Merging low temperature dynamics with the zero-temperature limit.

Clathrate hydrates, or cages comprised solely of water molecules, have long been investigated as a clean storage facility for hydrogen molecules. A breakthrough occurred when hydrogen molecules were experimentally placed within a structure-II clathrate hydrate, which sparked much interest to determine their feasibility for energy storage [Mao et al., Science 297, 2247-2249 (2002)]. We use Path Integral Molecular Dynamics (PIMD) and Langevin equation Path Integral Ground State (LePIGS) for finite temperature and zero-temperature studies, respectively, to determine parahydrogen occupancy properties in the small dodecahedral (512) and large hexakaidecahedral (51264) sized cages that comprise the structure-II unit cell. We look at energetic and structural properties of small clusters of hydrogen, treated as point-like particles, confined within each of the different sized clathrates, and treated as rigid, to determine energetic and structural properties in the zero-temperature limit. Our predicted hydrogen occupancy within these two cage sizes is consistent with previous literature values. We then calculate the energies as a function of temperature and merge the low temperature results calculated using finite temperature PIMD with the zero-temperature results using LePIGS, demonstrating that the two methods are compatible.

A path integral ground state approach for asymmetric top rotors with nuclear spin symmetry: Application to water chains

A path integral ground state (PIGS) approach for the simulation of asymmetric top rotors is presented. The method is based on Monte Carlo sampling of angular degrees of freedom. A symmetry-adapted rotational density matrix is used to account for nuclear spin statistics. To illustrate the method, ground-state properties of collections of para-water molecules confined to a one-dimensional lattice are computed. Those include energetic and structural observables. An advantage of the PIGS method is that expectation values can be obtained directly since the square of the wavefunction is sampled during a simulation. To benchmark the method, ground state energies and orientational distributions are computed using exact diagonalization for a single para-water molecule in an external field using a finite basis of symmetric top eigenfunctions. Benchmark results are also provided for N = 2 para-water molecules pinned to lattice sites at various distances to sample the crossover from hydrogen bonding to the dipole-dipole interaction regime. Excellent agreement between the PIGS results and the finite basis set calculations is observed. A thorough analysis of the convergence in terms of the imaginary time propagation length and systematic Trotter error is performed. The PIGS approach is then applied to a chain of N = 11 water molecules, and an equation of state is constructed in terms of the intermolecular separation. Ordering effects are also studied, and a transition between hydrogen bonding to dipole-dipole alignment is observed. The method is scalable and can also be applied in higher dimensions.

A path integral ground state replica trick approach for the computation of entanglement entropy of dipolar linear rotors.

We calculate the second Rényi entanglement entropy for systems of interacting linear rotors in their ground state as a measure of entanglement for continuous rotational degrees of freedom. The entropy is defined in relation to the purity of a subsystem in a bipartite quantum system, and to compute it, we compare two sampling ensembles based on the path integral ground state (PIGS) formalism. This scheme centers on the replica trick and is aided by the ratio trick, both developed in this context by Hastings et al. [Phys. Rev. Lett. 104, 157201 (2010)]. We study a system composed of linear quantum rotors on a lattice in one dimension, interacting via an anisotropic dipole-dipole potential. The ground state second Rényi entropies estimated by PIGS are benchmarked against those from the density matrix renormalization group for various interaction strengths and system sizes. We find that the entropy grows with an increase in interaction strength, and for large enough systems, it appears to plateau near log(2). We posit that the limiting case of many strongly interacting rotors behaves akin to a lattice of two-level particles in a cat state, in which one naturally finds an entanglement entropy of log(2).

Path integral molecular dynamic simulation of flexible molecular systems in their ground state: Application to the water dimer.

We extend the Langevin equation Path Integral Ground State (LePIGS), a ground state quantum molecular dynamics method, to simulate flexible molecular systems and calculate both energetic and structural properties. We test the approach with the H2O and D2O monomers and dimers. We systematically optimize all simulation parameters and use a unity trial wavefunction. We report ground state energies, dissociation energies, and structural properties using three different water models, two of which are empirically based, q-TIP4P/F and q-SPC/Fw, and one which is ab initio, MB-pol. We demonstrate that our energies calculated from LePIGS can be merged seamlessly with low temperature path integral molecular dynamics calculations and note the similarities between the two methods. We also benchmark our energies against previous diffusion Monte Carlo calculations using the same potentials and compare to experimental results. We further demonstrate that accurate vibrational energies of the H2O and D2O monomer can be calculated from imaginary time correlation functions generated from the LePIGS simulations using solely the unity trial wavefunction.

Quantum phases of dipolar rotors on two-dimensional lattices.

The quantum phase transitions of dipoles confined to the vertices of two-dimensional lattices of square and triangular geometry is studied using path integral ground state quantum Monte Carlo. We analyze the phase diagram as a function of the strength of both the dipolar interaction and a transverse electric field. The study reveals the existence of a class of orientational phases of quantum dipolar rotors whose properties are determined by the ratios between the strength of the anisotropic dipole-dipole interaction, the strength of the applied transverse field, and the rotational constant. For the triangular lattice, the generic orientationally disordered phase found at zero and weak values of both dipolar interaction strength and applied field is found to show a transition to a phase characterized by net polarization in the lattice plane as the strength of the dipole-dipole interaction is increased, independent of the strength of the applied transverse field, in addition to the expected transition to a transverse polarized phase as the electric field strength increases. The square lattice is also found to exhibit a transition from a disordered phase to an ordered phase as the dipole-dipole interaction strength is increased, as well as the expected transition to a transverse polarized phase as the electric field strength increases. In contrast to the situation with a triangular lattice, on square lattices, the ordered phase at high dipole-dipole interaction strength possesses a striped ordering. The properties of these quantum dipolar rotor phases are dominated by the anisotropy of the interaction and provide useful models for developing quantum phases beyond the well-known paradigms of spin Hamiltonian models, implementing in particular a novel physical realization of a quantum rotor-like Hamiltonian that possesses an anisotropic long range interaction.

Dipolar Bose Supersolid Stripes.

We study the superfluid properties of a system of fully polarized dipolar bosons moving in the XY plane. We focus on the general case where the polarization field forms an arbitrary angle α with respect to the Z axis, while the system is still stable. We use the diffusion MonteCarlo and the path integral ground state methods to evaluate the one-body density matrix and the superfluid fractions in the region of the phase diagram where the system forms stripes. Despite its oscillatory behavior, the presence of a finite large-distance asymptotic value in the s-wave component of the one-body density matrix indicates the existence of a Bose condensate. The superfluid fraction along the stripes direction is always close to 1, while in the Y direction decreases to a small value that is nevertheless different from zero. These two facts confirm that the stripe phase of the dipolar Bose system is a clear candidate for an intrinsic supersolid without the presence of defects as described by the Andreev-Lifshitz mechanism.

Number-squeezed and fragmented states of strongly interacting bosons in a double well

We present path integral ground state (PIGS) quantum Monte Carlo calculations for the ground state ($T = 0$) properties of repulsively interacting bosons in a three-dimensional external double well potential over a range of interaction strengths and potential parameters. We focus our calculations on ground state number statistics and the one-body density matrix in order to understand the level of squeezing and fragmentation that the system exhibits as a function of interaction strength. We compare our PIGS results to both a two-mode model and a recently-proposed eight-mode model. For weak interactions, the various models agree with the numerically exact PIGS simulations. However, the models fail to correctly predict the amount of squeezing and fragmentation exhibited by the PIGS simulations for strong interactions. One novel and somewhat surprising result from our simulations involves the relationship between squeezing and interaction strength: rather than a monotonic relationship between these quantities, we find that for certain barrier heights the squeezing increases as a function of interaction strength until it reaches a maximum, after which it decreases again. We also see a similar relationship between fragmentation and interaction strength. We discuss the physical mechanisms that account for this behavior and the implications for the design of atom interferometers, which can use squeezed states to reduce measurement uncertainty.

Luttinger parameter of quasi-one-dimensional para−H2

We have studied the ground-state properties of para-hydrogen in one dimension and in quasi-one-dimensional configurations using the path integral ground state Monte Carlo method. This method produces zero-temperature exact results for a given interaction and geometry. The quasi-one-dimensional setup has been implemented in two forms: the inner channel inside a carbon nanotube coated with H$_2$ and a harmonic confinement of variable strength. Our main result is the dependence of the Luttinger parameter on the density within the stable regime. Going from one dimension to quasi-one dimension, keeping the linear density constant, produces a systematic increase of the Luttinger parameter. This increase is however not enough to reach the superfluid regime and the system always remain in the quasi-crystal regime, according to Luttinger liquid theory.

Droplets of Trapped Quantum Dipolar Bosons.

Strongly interacting systems of dipolar bosons in three dimensions confined by harmonic traps are analyzed using the exact path integral ground-state MonteCarlo method. By adding a repulsive two-body potential, we find a narrow window of interaction parameters leading to stable ground-state configurations of droplets in a crystalline arrangement. We find that this effect is entirely due to the interaction present in the Hamiltonian without resorting to additional stabilizing mechanisms or specific three-body forces. We analyze the number of droplets formed in terms of the Hamiltonian parameters, relate them to the corresponding s-wave scattering length, and discuss a simple scaling model for the density profiles. Our results are in qualitative agreement with recent experiments showing a quantum Rosensweig instability in trapped Dy atoms.

Roton Excitations and the Fluid–Solid Phase Transition in Superfluid 2D Yukawa Bosons

We compute several ground state properties and the dynamical structure factor of a 0-temperature system of Bosons interacting with the 2D screened Coulomb (2D-SC) potential. We resort to the exact shadow path integral ground state (SPIGS) quantum Monte Carlo method to compute the imaginary-time correlation function of the model, and to the genetic algorithm via falsification of theories (GIFT) to retrieve the dynamical structure factor. We provide a detailed comparison of ground state properties and collective excitations of 2D-SC and 4He atoms. The roton energy of the 2D-SC system is an increasing function of density, and not a decreasing one as in 4He. This result is in contrast with the view that the roton is the soft mode of the fluid-solid transition. We uncover a remarkable quasi-universality of backflow and of other properties when expressed in terms of the amount of short range order as quantified by the height of the first peak of the static structure factor.

Raman Vibrational Shifts of Small Clusters of Hydrogen Isotopologues.

Raman vibrational shifts of small parahydrogen (pH2), orthodeuterium (oD2), and paratritium (pT2) clusters with respect to the free molecules are calculated by combining a first order perturbation theory approach with Langevin equation Path Integral Ground State (LePIGS) simulations [ J. Phys. Chem. A 2013 , 117 , 7461 ]. Our theoretical predictions are compared to existing cryogenic free jet expansion results for pure (pH2)N clusters [ Phys. Rev. Lett. 2004 , 92 , 223401 ] and to new measurements for (oD2)N clusters reported here. This method has been successfully used before to predict the Raman vibrational shifts of (pH2)N clusters [ J. Chem. Phys. 2014 , 141 , 014310 ]. The 6-D interaction potential of Hinde [ J. Chem. Phys. 2008 , 128 , 154308 ] is reduced to 1-D using the Adiabatic Hindered Rotor approximation to yield effective pair potentials for both molecules being in the ground vibrational state, and for one of them carrying one quantum of vibrational excitation. These reduced 1-D potentials are fitted to a Morse Long Range analytic form for later convenience. Good agreement between experiment and theory is found for the smaller clusters, but significant deviations remain for the larger ones.

Phase diagram of dipolar bosons in two dimensions with tilted polarization

We analyze the ground state of a system of dipolar bosons moving in the $XY$ plane and such that their dipolar moments are all aligned in a fixed direction in space. We focus on the general case where the polarization field forms a generic angle $\alpha$ with respect to the $Z$ axis. We use the Path Integral Ground State method to analyze the static properties of the system as both $\alpha$ and the density $n$ vary over a wide range were the system is stable. We use the maximum of the static structure function as an order parameter to characterize the different phases and the transition lines among them. We find that aside of a superfluid gas and a solid phase, the system reaches a stripe phase at large tilting angles that is entirely induced by the anisotropic character of the interaction. We also show that the quantum phase transition from the gas to the stripe phase is of second order, and report approximate values for the critical exponents.