Quantum Phase Transition with a Simple Variational Ansatz

We study the zero-temperature quantum phase transition between liquid and hcp solid helium-4. We use the variational method with a simple yet exchange-symmetric and fully explicit wavefunction. It is found that the optimized wavefunction undergoes spontaneous symmetry breaking and describes the quantum solidification of helium at 22 atm. The explicit form of the wavefunction allows to consider various contributions to the phase transition. We find that the employed wavefunction is an excellent candidate for describing both a first-order quantum phase transition and the ground state of a Bose solid.

A microscopic description of vacancies in solid 4He

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 vacancy 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 elastic constants of solid 4He under pressure: a diffusion Monte Carlo study

We study the elasticity of perfect 4He at zero-temperature using the diffusion Monte Carlo method and a realistic semi-empirical pairwise potential to describe the He-He interactions. Specifically, we calculate the value of the elastic constants of hcp helium C_{ij} as a function of pressure up to 110 bar. It is found that the pressure dependence of all five non-zero C_{ij} is linear and we provide accurate parametrization of each of them. Our elastic constants results are compared to previous variational calculations and low-temperature measurements and in general notably good agreement is found among them. Furthermore, we report T = 0 results for the Gruneisen parameters, sound velocities and Debye temperature of hcp 4He. This work represents the first of a series of computational studies aimed at thoroughly characterizing the response of solid helium to external stress-strain.

Ground-state properties and superfluidity of two- and quasi two-dimensional solid 4He

In a recent study we have reported a new type of trial wave function symmetric under the exchange of particles and which is able to describe a supersolid phase. In this work, we use the diffusion Monte Carlo method and this model wave function to study the properties of solid 4He in two- and quasi two-dimensional geometries. In the purely two-dimensional case, we obtain results for the total ground-state energy and freezing and melting densities which are in good agreement with previous exact Monte Carlo calculations performed with a slightly different interatomic potential model. We calculate the value of the zero-temperature superfluid fraction rho_{s} / rho of 2D solid 4He and find that it is negligible in all the considered cases, similarly to what is obtained in the perfect (free of defects) three-dimensional crystal using the same computational approach. Interestingly, by allowing the atoms to move locally in the perpendicular direction to the plane where they are confined to zero-point oscillations (quasi two-dimensional crystal) we observe the emergence of a finite superfluid density that coexists with the periodicity of the system.

The kinetics of homogeneous melting beyond the limit of superheating

Molecular dynamics simulation is used to study the time-scales involved in the homogeneous melting of a superheated crystal. The interaction model used is an embedded-atom model for Fe developed in previous work, and the melting process is simulated in the microcanonical $(N, V, E)$ ensemble. We study periodically repeated systems containing from 96 to 7776 atoms, and the initial system is always the perfect crystal without free surfaces or other defects. For each chosen total energy $E$ and number of atoms $N$, we perform several hundred statistically independent simulations, with each simulation lasting for between 500 ps and 10 ns, in order to gather statistics for the waiting time $tau_{rm w}$ before melting occurs. We find that the probability distribution of $tau_{rm w}$ is roughly exponential, and that the mean value $<tau_{rm w} >$ depends strongly on the excess of the initial steady temperature of the crystal above the superheating limit identified by other researchers. The mean $<tau_{rm w}>$ also depends strongly on system size in a way that we have quantified. For very small systems of $sim 100$ atoms, we observe a persistent alternation between the solid and liquid states, and we explain why this happens. Our results allow us to draw conclusions about the reliability of the recently proposed Z method for determining the melting properties of simulated materials, and to suggest ways of correcting for the errors of the method.

Supersolidity in quantum films adsorbed on graphene and graphite

Using quantum Monte Carlo we have studied the superfluid density of the first layer of $^4$He and H$_2$ adsorbed on graphene and graphite. Our main focus has been on the equilibrium ground state of the system, which corresponds to a registered $sqrt3 times sqrt3$ phase. The perfect solid phase of H$_2$ shows no superfluid signal whereas $^4$He has a finite but small superfluid fraction (0.67%). The introduction of vacancies in the crystal makes the superfluidity increase, showing values as large as 14% in $^4$He without destroying the spatial solid order.

Properties of Vacancy Formation in He-4 hcp Crystals at Zero Temperature and Fixed Pressure

Equation of state of He-4 hcp crystals with vacancies is determined at zero temperature using the diffusion Monte Carlo technique, an exact ground state zero-temperature method. This allows us to extract the formation enthalpy and isobaric formation energy of a single vacancy in otherwise perfect helium solid. Results were obtained for pressures up to 160 bar. The isobaric formation energy is found to reach a minimum near 57 bar where it is equal to $10.5pm 1.2$ K. At the same pressure, the vacancy formation volume exhibits a maximum and reaches the volume of the unit cell. This pressure coincides with the pressure interval over which a peak in the supersolid fraction of He-4 was observed in a recent experiment.

Melting properties of a simple tight-binding model of transition metals: I.The region of half-filled d-band

We present calculations of the free energy, and hence the melting properties, of a simple tight-binding model for transition metals in the region of d-band filling near the middle of a d-series, the parameters of the model being designed to mimic molybdenum. The melting properties are calculated for pressures ranging from ambient to several Mbar. The model is intended to be the simplest possible tight-binding representation of the two basic parts of the energy: first, the pairwise repulsion due to Fermi exclusion; and second, the d-band bonding energy described in terms of an electronic density of states that depends on structure. In addition to the number of d-electrons, the model contains four parameters, which are adjusted to fit the pressure dependent d-band width and the zero-temperature pressure-volume relation of Mo. We show that the resulting model reproduces well the phonon dispersion relations of Mo in the body-centred-cubic structure, as well as the radial distribution function of the high-temperature solid and liquid given by earlier first-principles simulations. Our free-energy calculations start from the free energy of the liquid and solid phases of the purely repulsive pair-potential model, without d-band bonding. The free energy of the full tight-binding model is obtained from this by thermodynamic integration. The resulting melting properties of the model are quite close to those given by earlier first-principles work on Mo. An interpretation of these melting properties is provided by showing how they are related to those of the purely repulsive model.

Comment on Molybdenum at High Pressure and Temperature: Melting from Another Solid Phase

There has been a major controversy over the past seven years about the high-pressure melting curves of transition metals. Static compression (diamond-anvil cell: DAC) experiments up to the Mbar region give very low melting slopes dT_m/dP, but shock-wave (SW) data reveal transitions indicating much larger dT_m/dP values. Ab initio calculations support the correctness of the shock data. In a very recent letter, Belonoshko et al. propose a simple and elegant resolution of this conflict for molybdenum. Using ab initio calculations based on density functional theory (DFT), they show that the high-P/high-T phase diagram of Mo must be more complex than was hitherto thought. Their calculations give convincing evidence that there is a transition boundary between the normal bcc structure of Mo and a high-T phase, which they suggest could be fcc. They propose that this transition was misinterpreted as melting in DAC experiments. In confirmation, they note that their boundary also explains a transition seen in the SW data. We regard Belonoshko et al.s Letter as extremely important, but we note that it raises some puzzling questions, and we believe that their proposed phase diagram cannot be completely correct. We have calculated the Helmholtz and Gibbs free energies of the bcc, fcc and hcp phases of Mo, using essentially the same quasiharmonic methods as used by Belonoshko et al.; we find that at high-P and T Mo in the hcp structure is more stable than in bcc or fcc.

Two-dimensional molecular para-hydrogen and ortho-deuterium at zero temperature

We study molecular para-hydrogen (p-${rm H_{2}}$) and ortho-deuterium (o-${rm D_{2}}$) in two dimensions and in the limit of zero temperature by means of the diffusion Monte Carlo method. We report energetic and structural properties of both systems like the total and kinetic energy per particle, radial pair distribution function, and Lindemanns ratio in the low pressure regime. By comparing the total energy per particle as a function of the density in liquid and solid p-${rm H_{2}}$, we show that molecular para-hydrogen, and also ortho-deuterium, remain solid at zero temperature. Interestingly, we assess the quality of three different symmetrized trial wave functions, based on the Nosanow-Jastrow model, in the p-${rm H_{2}}$ solid film at the variational level. In particular, we analyze a new type of symmetrized trial wave function which has been used very recently to describe solid $^{4}$He and found that also characterizes hydrogen satisfactorily. With this wave function, we show that the one-body density matrix $varrho_{1} (r)$ of solid p-${rm H_{2}}$ possesses off-diagonal long range order, with a condensate fraction that increases sizably in the negative pressure regime.