Search

We present a combined theoretical and experimental study of the high-pressure behavior of thallium. X-ray diffraction experiments have been carried out at room temperature up to 125 GPa using diamond-anvil cells, nearly doubling the pressure range of previous experiments. We have confirmed the hcp-fcc transition at 3.5 GPa and determined that the fcc structure remains stable up to the highest pressure attained in the experiments. In addition, HP-HT experiments have been performed up to 8 GPa and 700 K by using a combination of x-ray diffraction and a resistively heated diamond-anvil cell. Information on the phase boundaries is obtained, as well as crystallographic information on the HT bcc phase. The equation of state for different phases is reported. Ab initio calculations have also been carried out considering several potential high-pressure structures. They are consistent with the experimental results and predict that, among the structures considered in the calculations, the fcc structure of thallium is stable up to 4.3 TPa. Calculations also predict the post-fcc phase to have a close-packed orthorhombic structure above 4.3 TPa., Comment: 29 pages, 14 figures

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.

We determined the zero-temperature phase diagram of D$_2$ physisorbed on graphane using the diffusion Monte Carlo method. The substrate used was C-graphane, an allotropic form of the compound that has been experimentally obtained through hydrogenation of graphene. We found that the ground state is the $\delta$ phase, a commensurate structure observed experimentally when D$_2$ is adsorbed on graphite, and not the registered $\sqrt 3 \times \sqrt 3$ structure characteristic of H$_2$ on the same substrate., Comment: 6 figures

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., Comment: 13 pages, accepted for publication in J. Low Temp. Phys., Special Issue on Supersolids

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., Comment: 21 pages, 7 figures

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., Comment: 16 pages, 8 figures

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., Comment: 19 pages, 8 figures

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., Comment: 5 pages, accepted for publication in PRB

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.5\pm 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.

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: 34 pages, 12 figures. Accepted for publication in Journal of Chemical Physics

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., Comment: 1 page, 1 figure. submitted to Phys. Rev. Lett

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 Lindemann's 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., Comment: 11 pages, 9 figures

The ground state of solid $^4$He is studied using the diffusion Monte Carlo method and a new trial wave function able to describe the supersolid. The new wave function is symmetric under the exchange of particles and reproduces the experimental equation of state. Results for the one-body density matrix show the existence of off-diagonal long-range order with a very small condensate fraction $\sim 10^{-4}$. The superfluid density of the commensurate system is below our resolution threshold, $\rho_s/\rho < 10^{-5}$. With a 1% concentration of vacancies the superfluid density is manifestly larger, $\rho_s/\rho=3.2(1) \cdot 10^{-3}$., Comment: 5 pages, 4 figures

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., Comment: 28 pages and 10 figures. Submitted to Phys. Rev. B

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, Lindemann's 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., Comment: 20 pages, 7 figures. To be published in Physical Review B

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., Comment: 29 pages, 9 figures. To be published in Journal of Physics: Condensed Matter

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., Comment: 8 pages, 6 figures

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., Comment: 27 pages, 10 figures. to be published in Journal of Chemical Physics

We present a zero-temperature quantum Monte Carlo calculation of liquid $^4$He immersed in an array of confining potentials. These external potentials are centered in the lattice sites of a fcc solid geometry and, by modifying their well depth and range, the system evolves from a liquid phase towards a progressively localized system which mimics a solid phase. The superfluid density decreases with increasing order, reaching a value $ \rho_{\rm s}/\rho = 0.079(16)$ when the Lindemann's ratio of the model equals the experimental value for solid $^4$He., Comment: 5 pages,5 figures

A diffusion Monte Carlo simulation of superfluid $^4$He at zero temperature and pressures up to 275 bar is presented. Increasing the pressure beyond freezing ($\sim$ 25 bar), the liquid enters the overpressurized phase in a metastable state. In this regime, we report results of the equation of state and the pressure dependence of the static structure factor, the condensate fraction, and the excited-state energy corresponding to the roton. Along this large pressure range, both the condensate fraction and the roton energy decrease but do not become zero. The roton energies obtained are compared with recent experimental data in the overpressurized regime., Comment: 5 pages, accepted for publication in Phys. Rev. Lett

Hydrogen single-particle dynamics in solid LiH at T=20 K has been studied through the incoherent inelastic neutron scattering technique. A careful analysis of the scattering data has allowed for the determination of a reliable hydrogen-projected density of phonon states and, from this, of three relevant physical quantities: mean squared displacement, mean kinetic energy, and Einstein frequency. In order to interpret these experimental findings, a fully-quantum microscopic calculation has been carried out using the variational Monte Carlo method. The agreement achieved between neutron scattering data and Monte Carlo estimates is good. In addition, a purely harmonic calculation has been also performed via the same Monte Carlo code, but anharmonic effects in H dynamics were not found relevant. The possible limitations of the present semi-empirical potentials are finally discussed., Comment: RevTex, 19 pages, accepted for publication in PRB