No Arabic abstract
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.
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.
Several research groups have recently reported {em ab initio} calculations of the melting properties of metals based on density functional theory, but there have been unexpectedly large disagreements between results obtained by different approaches. We analyze the relations between the two main approaches, based on calculation of the free energies of solid and liquid and on direct simulation of the two coexisting phases. Although both approaches rely on the use of classical reference systems consisting of parameterized empirical interaction models, we point out that in the free energy approach the final results are independent of the reference system, whereas in the current form of the coexistence approach they depend on it. We present a scheme for correcting the predictions of the coexistence approach for differences between the reference and {em ab initio} systems. To illustrate the practical operation of the scheme, we present calculations of the high-pressure melting properties of iron using the corrected coexistence approach, which agree closely with earlier results from the free energy approach. A quantitative assessment is also given of finite-size errors, which we show can be reduced to a negligible size.
Understanding the behavior of molecular systems under pressure is a fundamental problem in condensed matter physics. In the case of nitrogen, the determination of the phase diagram and in particular of the melting line, are largely open problems. Two independent experiments have reported the presence of a maximum in the nitrogen melting curve, below 90 GPa, however the position and the interpretation of the origin of such maximum differ. By means of ab initio molecular dynamics simulations based on density functional theory and thermodynamic integration techniques, we have determined the phase diagram of nitrogen in the range between 20 and 100 GPa. We find a maximum in the melting line, related to a transformation in the liquid, from molecular N_2 to polymeric nitrogen accompanied by an insulator-to-metal transition.
First-principles calculation predict that olivine Li4MnFeCoNiP4O16 has ferrotoroidic characteristic and ferrimagnetic configuration with magnetic moment of 1.56 muB per formula unit. The ferrotoroidicity of this material makes it a potential candidate for magnetoelectric materials . Based on the orbital-resolved density of states for the transtion-metal ions in Li4MnFeCoNiP4O16, the spin configuration for Mn2+,Fe3+,Co2+, and Ni2+ is t2g3eg2, t2g3eg2,t2g1t2g3eg1eg2, and t2g2t2g3eg1eg2, respectively. Density functional theory plus U (DFT+U) shows a indirect band gap of 1.25 eV in this predicted material, which is not simply related to the electronic conductivity in terms of being used as cathode material in rechargeable Li-ion batteries.
We present a novel ab initio non-equilibrium approach to calculate the current across a molecular junction. The method rests on a wave function based full ab initio description of the central region of the junction combined with a tight binding approximation for the electrodes in the frame of the Keldysh Greens function formalism. Our procedure is demonstrated for a dithiolethine molecule between silver electrodes. The main conducting channel is identified and the full current-voltage characteristic is calculated.