Abstract and subjects
The unified analytic melt-shear model that we introduced a decade ago is generalized to multi-phase materials. A new scheme for calculating the values of the model parameters for both the cold (T = 0) shear modulus (G) and the melting temperature at all densities (rho) is developed. The generalized melt-shear model is applied to molybdenum, a multi-phase material with a body-centered cubic (bcc) structure at low rho which loses its dynamical stability with increasing pressure (P) and is therefore replaced by another (dynamically stable) solid structure at high rho. One of the candidates for the high-rho structure of Mo is face-centered cubic (fcc). The model is compared to (i) our ab initio results on the cold shear modulus of both bcc-Mo and fcc-Mo as a function of rho, and (ii) the available theoretical results on the melting of bcc-Mo and our own quantum molecular dynamics (QMD) simulations of one melting point of fcc-Mo. Our generalized model of G(rho, T) is used to calculate the shear modulus of bcc-Mo along its principal Hugoniot. It predicts that G of bcc-Mo increases with P up to similar to 240 GPa and then decreases at higher P. This behavior is intrinsic to bcc-Mo and does not require the introduction of another solid phase such as Phase II suggested by Errandonea et al. Generalized melt-shear models for Ta and W also predict an increase in G followed by a decrease along the principal Hugoniot, hence this behavior may be typical for transition metals with ambient bcc structure that dynamically destabilize at high P. Thus, we concur with the conclusion reached in several recent papers (Nguyen et al., Zhang et al., Wang et al.) that no solid-solid phase transition can be definitively inferred on the basis of sound velocity data from shock experiments on Mo. Finally, our QMD simulations support the validity of the phase diagram of Mo suggested by Zeng et al.