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.

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.

Zero-temperature generalized phase diagram of the 4d transition metals under pressure

We use an accurate implementation of density functional theory (DFT) to calculate the zero-temperature generalized phase diagram of the 4$d$ series of transition metals from Y to Pd as a function of pressure $P$ and atomic number $Z$. The implementation used is full-potential linearized augmented plane waves (FP-LAPW), and we employ the exchange-correlation functional recently developed by Wu and Cohen. For each element, we obtain the ground-state energy for several crystal structures over a range of volumes, the energy being converged with respect to all technical parameters to within $sim 1$ meV/atom. The calculated transition pressures for all the elements and all transitions we have found are compared with experiment wherever possible, and we discuss the origin of the significant discrepancies. Agreement with experiment for the zero-temperature equation of state is generally excellent. The generalized phase diagram of the 4$d$ series shows that the major boundaries slope towards lower $Z$ with increasing $P$ for the early elements, as expected from the pressure induced transfer of electrons from $sp$ states to $d$ states, but are almost independent of $P$ for the later elements. Our results for Mo indicate a transition from bcc to fcc, rather than the bcc-hcp transition expected from $sp$-$d$ transfer.

Atomic kinetic energy, momentum distribution and structure of solid neon at zero-temperature

We report on the calculation of the ground-state atomic kinetic energy, $E_{k}$, and momentum distribution of solid Ne by means of the diffusion Monte Carlo method and Aziz HFD-B pair potential. This approach is shown to perform notably for this crystal since we obtain very good agreement with respect to experimental thermodynamic data. Additionally, we study the structural properties of solid Ne at densities near the equilibrium by estimating the radial pair-distribution function, Lindemanns ratio and atomic density profile around the positions of the perfect crystalline lattice. Our value for $E_{k}$ at the equilibrium density is $41.51(6)$ K, which agrees perfectly with the recent prediction made by Timms {it et al.}, $41(2)$ K, based on their deep-inelastic neutron scattering experiments carried out over the temperature range $4 - 20$ K, and also with previous path integral Monte Carlo results obtained with the Lennard-Jones and Aziz HFD-C2 atomic pairwise interactions. The one-body density function of solid Ne is calculated accurately and found to fit perfectly, within statistical uncertainty, to a Gaussian curve. Furthermore, we analyze the degree of anharmonicity of solid Ne by calculating some of its microscopic ground-state properties within traditional harmonic approaches. We provide insightful comparison to solid $^4$He in terms of the Debye model, in order to size the relevance of anharmonic effects in Ne.

Zero-temperature equation of state of solid 4He at low and high pressures

We study the zero-temperature equation of state (EOS) of solid 4He in the hexagonal closed packet (hcp) phase over the 0-57 GPa pressure range by means of the Diffusion Monte Carlo (DMC) method and the semi-empirical Aziz pair potential HFD-B(HE). In the low pressure regime (P ~ 0-1 GPa) we assess excellent agreement with experiments and we give an accurate description of the atomic kinetic energy, Lindemann ratio and Debye temperature over a wide range of molar volumes (22-6 cm^{3}/mol). However, on moving to higher pressures our calculated P-V curve presents an increasingly steeper slope which ultimately provides differences within ~40 % with respect to measurements. In order to account for many-body interactions arising in the crystal with compression which are not reproduced by our model, we perform additional electronic density-functional theory (DFT) calculations for correcting the computed DMC energies in a perturbative way. We explore both generalized gradient and local density approximations (GGA and LDA, respectively) for the electronic exchange-correlation potential. By proceeding in this manner, we show that discrepancies with respect to high pressure data are reduced to 5-10 % with few computational extra cost. Further comparison between our calculated EOSs and ab initio curves deduced for the perfect crystal and corrected for the zero-point motion of the atoms enforces the reliability of our approach.

Melting curve and Hugoniot of molybdenum up to 400 GPa by ab initio simulations

We report ab initio calculations of the melting curve and Hugoniot of molybdenum for the pressure range 0-400 GPa, using density functional theory (DFT) in the projector augmented wave (PAW) implementation. We use the ``reference coexistence technique to overcome uncertainties inherent in earlier DFT calculations of the melting curve of Mo. Our calculated melting curve agrees well with experiment at ambient pressure and is consistent with shock data at high pressure, but does not agree with the high pressure melting curve from static compression experiments. Our calculated P(V) and T(P) Hugoniot relations agree well with shock measurements. We use calculations of phonon dispersion relations as a function of pressure to eliminate some possible interpretations of the solid-solid phase transition observed in shock experiments on Mo.

Ab initio melting curve of molybdenum by the phase coexistence method

We report ab initio calculations of the melting curve of molybdenum for the pressure range 0-400 GPa. The calculations employ density functional theory (DFT) with the Perdew-Burke-Ernzerhof exchange-correlation functional in the projector augmented wave (PAW) implementation. We present tests showing that these techniques accurately reproduce experimental data on low-temperature b.c.c. Mo, and that PAW agrees closely with results from the full-potential linearized augmented plane-wave implementation. The work attempts to overcome the uncertainties inherent in earlier DFT calculations of the melting curve of Mo, by using the ``reference coexistence technique to determine the melting curve. In this technique, an empirical reference model (here, the embedded-atom model) is accurately fitted to DFT molecular dynamics data on the liquid and the high-temperature solid, the melting curve of the reference model is determined by simulations of coexisting solid and liquid, and the ab initio melting curve is obtained by applying free-energy corrections. Our calculated melting curve agrees well with experiment at ambient pressure and is consistent with shock data at high pressure, but does not agree with the high pressure melting curve deduced from static compression experiments. Calculated results for the radial distribution function show that the short-range atomic order of the liquid is very similar to that of the high-T solid, with a slight decrease of coordination number on passing from solid to liquid. The electronic densities of states in the two phases show only small differences. The results do not support a recent theory according to which very low dTm/dP values are expected for b.c.c. transition metals because of electron redistribution between s-p and d states.