Molecular para-hydrogen has been proposed theoretically as a possible candidate for superfluidity, but the eventual superfluid transition is hindered by its crystallization. In this work, we study a metastable non crystalline phase of bulk p-H2 by me
ans of the Path Integral Monte Carlo method in order to investigate at which temperature this system can support superfluidity. By choosing accurately the initial configuration and using a non commensurate simulation box, we have been able to frustrate the formation of the crystal in the simulated system and to calculate the temperature dependence of the one-body density matrix and of the superfluid fraction. We observe a transition to a superfluid phase at temperatures around 1 K. The limit of zero temperature is also studied using the diffusion Monte Carlo method. Results for the energy, condensate fraction, and structure of the metastable liquid phase at T=0 are reported and compared with the ones obtained for the stable solid phase.
The changes that vacancies produce in the properties of hcp solid 4He are studied by means of quantum Monte Carlo methods. Our results show that the introduction of vacancies produces significant changes in the behavior of solid 4He, even when the va
cancy concentration is very small. We show that there is an onset temperature where the properties of incommensurate 4He change significantly. Below this temperature, we observe the emergence of off-diagonal long range order and a complete spatial delocalization of the vacancies. This temperature is quite close to the temperature where non-classical rotational inertia has been experimentally observed. Finally, we report results on the influence of vacancies in the elastic properties of hcp 4He at zero temperature.
The feasibility of path integral Monte Carlo ground state calculations with very few beads using a high-order short-time Greens function expansion is discussed. An explicit expression of the evolution operator which provides dramatic enhancements in
the quality of ground-state wave-functions is examined. The efficiency of the method makes possible to remove the trial wave function and thus obtain completely model-independent results still with a very small number of beads. If a single iteration of the method is used to improve a given model wave function, the result is invariably a shadow-type wave function, whose precise content is provided by the high-order algorithm employed.
We have investigated the ground state properties of solid $^4$He with the Shadow Path Integral Ground State method. This exact T=0 K projector method allows to describes quantum solids without introducing any a priori equilibrium position. We have fo
und that the efficiency in computing off-diagonal properties in the solid phase sensibly improves when the direct sampling of permutations, in principle not required, is introduced. We have computed the exact one-body density matrix (obdm) in large commensurate 4He crystal finding a decreasing condensate fraction with increasing imaginary time of projection, making our result not conclusive on the presence of Bose-Einstein condensation in bulk solid 4He. We can only give an upper bound of 2.5 times 10^-8 on the condensate fraction. We have exploited the SPIGS method to study also 4He crystal containing grain boundaries by computing the related surface energy and the obdm along these defects. We have found that also highly symmetrical grain boundaries have a finite condensate fraction. We have also derived a route for the estimation of the true equilibrium concentration of vacancies x_v in bulk T=0 K solid 4He, which is shown to be finite, x_v=0.0014(1) at the melting density, when computed with the variational shadow wave function technique.
