No Arabic abstract
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.
Hydrogen is the most abundant element in the universe, and its properties under conditions of high temperature and pressure are crucial to understand the interior of of large gaseous planets and other astrophysical bodies. At ultra high pressures solid hydrogen has been predicted to transform into a quantum fluid, because of its high zero point motion. Here we report first principles two phase coexistence and Z method determinations of the melting line of solid hydrogen in a pressure range spanning from 30 to 600 GPa. Our results suggest that the melting line of solid hydrogen, as derived from classical molecular dynamics simulations, reaches a minimum of 367 K at about 430 GPa, at higher pressures the melting line of the atomics Cs IV phase regain a positive slope. In view of the possible importance of quantum effects in hydrogen at such low temperatures, we also determined the melting temperature of the atomic CsIV phase at pressures of 400, 500, 600 GPa, employing Feynman path integral simulations. These result in a downward shift of the classical melting line by about 100 K, and hint at a possible secondary maximum in the melting line in the region between 500 and 600 GPa, testifying to the importance of quantum effects in this system. Combined, our results imply that the stability field of the zero temperature quantum liquid phase, if it exists at all, would only occur at higher pressures than previously thought.
This Comment points out a number of errors in the recent paper by Zarechnaya, Dubrovinskaia, Dubrovinsky, et al. (Phys. Rev. Lett. 102, 185501 (2009)). Results and conclusions presented by Zarechnaya et al. (2009) are either incorrect or have been presented before.
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.
The phase diagram of Zn has been explored up to 140 GPa and 6000 K, by combining optical observations, x-ray diffraction, and ab-initio calculations. In the pressure range covered by this study, Zn is found to retain a hexagonal close-packed crystal symmetry up to the melting temperature. The known decrease of the axial ratio of the hcp phase of Zn under compression is observed in x-ray diffraction experiments from 300 K up to the melting temperature. The pressure at which the axial ratio reaches the square root of 3 value, around 10 GPa, is slightly affected by temperature. When this axial ratio is reached, we observed that single crystals of Zn, formed at high temperature, break into multiple polycrystals. In addition, a noticeable change in the pressure dependence of the axial ratio takes place at the same pressure. Both phenomena could be caused by an isomorphic second-order phase transition induced by pressure in Zn. The reported melt curve extends previous results from 24 to 135 GPa. The pressure dependence obtained for the melting temperature is accurately described up to 135 GPa by using a Simon-Glatzel equation. The determined melt curve agrees with previous low-pressure studies and with shock-wave experiments, with a melting temperature of 5060 K at 135 GPa. Finally, a thermal equation of state is reported, which at room-temperature agrees with the literature.
This paper reports an investigation on the phase diagram and compressibility of wolframite-type tungstates by means of x-ray powder diffraction and absorption in a diamond-anvil cell and ab initio calculations. The diffraction experiments show that monoclinic wolframite-type MgWO4 suffers at least two phase transitions, the first one being to a triclinic polymorph with a structure similar to that of CuWO4 and FeMoO4-II. The onset of each transition is detected at 17.1 and 31 GPa. In ZnWO4 the onset of the monoclinic-triclinic transition has been also found at 15.1 GPa. These findings are supported by density-functional theory calculations, which predict the occurrence of additional transitions upon further compression. Calculations have been also performed for wolframite-type MnWO4, which is found to have an antiferromagnetic configuration. In addition, x-ray absorption and diffraction experiments as well as calculations reveal details of the local-atomic compression in the studied compounds. In particular, below the transition pressure the ZnO6 and equivalent polyhedra tend to become more regular, whereas the WO6 octahedra remain almost unchanged. Fitting the pressure-volume data we obtained the equation of state for the low-pressure phase of MgWO4 and ZnWO4. These and previous results on MnWO4 and CdWO4 are compared with the calculations, being the compressibility of wolframite-type tungstates systematically discussed. Finally Raman spectroscopy measurements and lattice dynamics calculations are presented for MgWO4.