Molecular para-hydrogen (p-H 2 ) 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 noncrystalline phase of bulk p-H 2 by means 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 noncommensurate 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 Ewald summation technique is generalised to power-law 1/|r| k potentials in three-, twoand one-dimensional geometries with explicit formulae for all the components of the sums. The cases of short-range, long-range and "marginal" interactions are treated separately. The jellium model, as a particular case of a charge-neutral system, is discussed and the explicit forms of the Ewald sums for such system are presented. A generalised form of the Ewald sums for a noncubic (nonsquare) simulation cell for three-(two-) dimensional geometry is obtained and its possible field of application is discussed. A procedure for the optimisation of the involved parameters in actual simulations is developed and an example of its application is presented.

We report a quantum Monte Carlo calculation of the phase diagram of bosons interacting with a repulsive inverse sixth power pair potential, a model for assemblies of Rydberg atoms in the local van der Waals blockade regime. The model can be parametrized in terms of just two parameters, the reduced density and temperature. Solidification happens to the fcc phase. At zero temperature the transition density is found with the diffusion Monte Carlo method at density $\rho = 3.9 (\hbar^2/m C_6)^{3/4} $, where $C_6$ is the strength of the interaction. The solidification curve at non-zero temperature is studied with the path integral Monte Carlo approach and is compared with transitions in corresponding harmonic and classical crystals. Relaxation mechanisms are considered in relation to present experiments, especially pertaining to hopping of the Rydberg excitation

We study the zero-temperature phase diagram of bosons interacting via screened Coulomb (Yukawa) potential by means of the diffusion Monte Carlo method. The Yukawa potential is used as a model interaction in the neutron matter, dusty plasmas, and charged colloids. As shown by Petrov et al. [Phys. Rev. Lett. 99, 130407 (2007)], interactions between weakly bound molecules of heavy and light fermionic atoms are described by an effective Yukawa potential with a strength related to the heavy-light mass ratio M/m, which might lead to crystallization in a two-dimensional geometry if the mass ratio of heavy-light fermions exceeds a certain critical value. In the present work we do a thorough study of the quantum three-dimensional Yukawa system. For strong interactions (equivalently, large mass ratios) the system experiences several phase transitions as the density is increased, passing from gas to solid and to gas phase again. Weakly interacting Yukawa particles do not crystallize at any density. We find the minimal interaction strength at which the crystallization happens. In terms of the two-component fermionic system, this strength corresponds to a heavy-light mass ratio of M/m ∼ 180, so that it is impossible to realize the gas-crystal transition in a conventional bulk system. For the Yukawa model of fermionic mixtures we also analyze the possibility of building molecular systems with very large effective mass ratios by confining the heavy component to a sufficiently deep optical lattice. We show how the effective mass of the heavy component can be made arbitrarily large by increasing the lattice depth, thus leading to a tunable effective mass ratio that can be used to realize a molecular superlattice.