Mo高压熔化的分子动力学模拟

 采用第一原理方法计算了钼在零温下的结构,表明钼在500 GPa以下一直保持bcc结构（常温）,与实验一致。在零压附近计算了E-V关系,利用Murnaghan物态方程拟合得到了零压体积及其模量,与实验结果符合得很好。采用第一原理分子动力学模拟了钼的高压熔化性质。采用NVT系综计算了128个原子的系统,初始构形为bcc结构,体积分别为0.015 48、0.012 19、0.010 98、0.009 84、0.009 10 nm³/atom,计算了几个温度点,拟合得到了熔化曲线,熔化温度明显高于金刚石压砧（DAC）实验结果;将初始构形改变为fcc结构,模拟其熔化特性,得到的熔化温度明显下降,与激光加载DAC实验结果一致,认为可能的原因是钼熔化后形成的液体结构类似于fcc结构,而不是常态时的bcc结构。     

Molybdenum at high pressure and temperature: melting from another solid phase

The Gibbs free energies of bcc and fcc Mo are calculated from first principles in the quasiharmonic approximation in the pressure range from 350 to 850 GPa at room temperatures up to 7500 K. It is found that Mo, stable in the bcc phase at low temperatures, has lower free energy in the fcc structure than in the bcc phase at elevated temperatures. Our density-functional-theory-based molecular dynamics simulations demonstrate that fcc melts at higher than bcc temperatures above 1.5 Mbar. Our calculated melting temperatures and bcc-fcc boundary are consistent with the Mo Hugoniot sound speed measurements. We find that melting occurs at temperatures significantly above the bcc-fcc boundary. This suggests an explanation of the recent diamond anvil cell experiments, which find a phase boundary in the vicinity of our extrapolated bcc-fcc boundary.     

Quasi-Ab initio molecular dynamic study of Fe melting

We have investigated the melting of hcp Fe at high pressure by employing molecular dynamics simulations in conjunction with the full potential linear muffin tin orbital method. Apart from being of fundamental value, the melting of iron at high pressure is also important for our understanding of the Earth. The subject of iron melting at high pressures is controversial. The experimental data for the iron melting temperature can be separated into two regions, "low" and "high." Here we present an ab initio simulated iron melting curve which is in agreement with the low temperatures at lower pressures, but is in excellent agreement with the high-mostly shockwave-temperatures at high pressures. A comparison with available data lends support to the presented iron melting curve.     

An ab initio molecular dynamics study of iron phases at high pressure and temperature

The crystal structure of iron, the major component of the Earth's inner core (IC), is unknown for the IC high pressure (P; 3.3-3.6 Mbar) and temperature (T; 5000-7000 K). There is mounting evidence that the hexagonal close-packed (hcp) phase of iron, stable at the high P of the IC and a low T, might be unstable under the IC conditions due to the impact of high T and impurities. Experiments at the IC P and T are difficult and do not provide a conclusive answer as regards the iron stability at the pressure of the IC and temperatures close to the iron melting curve. Recent theory provides contradictory results regarding the nature of the stable Fe phase. We investigated the possibility of body-centered cubic (bcc) phase stabilization at the P and T in the vicinity of the Fe melting curve by using ab initio molecular dynamics. Thermodynamic calculations, relying on the model of uncorrelated harmonic oscillators, provide nearly identical free energies within the error bars of our calculations. However, direct simulation of iron crystallization demonstrates that liquid iron freezes in the bcc structure at the P of the IC and T = 6000 K. All attempts to grow the hcp phase from the liquid failed. The mechanism of bcc stabilization is explained. This resolves most of the earlier confusion.     

Melting of dense sodium

High-pressure high-temperature synchrotron diffraction measurements reveal a maximum on the melting curve of Na in the bcc phase at approximately 31 GPa and 1000 K and a steep decrease in melting temperature in its fcc phase. The results extend the melting curve by an order of magnitude up to 130 GPa. Above 103 GPa, Na crystallizes in a sequence of phases with complex structures with unusually low melting temperatures, reaching 300 K at 118 GPa, and an increased melting temperature is observed with further increases in pressure.     

A thermodynamic lattice theory on melting curve and eutectic point of binary alloys. Application to fcc and bcc structure

A thermodynamic lattice theory has been developed for determination of the melting curves and eutectic points of binary alloys. Analytical expressions for the melting curves of binary alloys composed of constituent elements with the same structure have been derived from expressions for the ratio of root mean square fluctuation in atomic positions on the equilibrium lattice positions and the nearest neighbor distance. This melting curve provides information on Lindemann’s melting temperatures of binary alloys with respect to any proportion of constituent elements, as well as on their eutectic points. The theory has been applied to fcc and bcc structure. Numerical results for some binary alloys provide a good correspondence between the calculated and experimental phase diagrams, where the calculated results for Cu_1−x Ni _x agree well with the measured ones, and those for the other alloys are found to be in a reasonable agreement with experiment.     

Thermal equation of state, and melting and thermoelastic properties of bcc tantalum from molecular dynamics

We performed molecular dynamics simulations with the extended Finnis-Sinclair (EFS) potential to investigate thermal equation of state (EOS), and melting and thermoelastic properties of tantalum. The agreement of the obtained thermal EOS with experiments at ambient conditions is reasonably good. The EFS potential with the two-phase method also reproduced very satisfyingly the high-pressure melting curve, excellently consistent with both the experiments of melting temperature at ambient pressure and shock melting at high pressure. From molecular dynamics simulations, we also obtained the thermoelastic properties of Ta for temperatures up to 3000 K at ambient pressure. Fully including anharmonic effects in molecular dynamics, our calculated elastic constants are in excellent agreement with experimental data. Shear modulus G decreases quickly with increasing temperature.     

Properties of MgO at high pressures: Shell-model molecular dynamics simulation

Shell-model molecular dynamics (MD) simulation has been performed to investigate the melting of the major Earth-forming mineral: periclase (MgO), at elevated temperatures and high pressures, based on the thermal instability analysis. The interatomic potential is taken to be the sum of pair-wise additive Coulomb, van der Waals attraction, and repulsive interactions. The MD simulation with selected Lewis–Catlow (LC) potential parameters is found to be very successful in describing the melting behavior for MgO, by taking account of the overheating of a crystalline solid at ambient pressure. The thermodynamic melting curve is estimated on the basis of the thermal instability MD simulations and compared with the available experimental data and other theoretical results in the pressure ranges 0–150 GPa. Our simulated melting curve of MgO is consistent with results obtained from Lindemann melting equation and two-phase simulated data at constant pressure by Belonoshko and Dubrovinsky, in the pressure below 20 GPa. The extrapolated melting temperatures in the lower mantle are in good agreement with the results obtained from Wang's empirical model up to 100 GPa. Compared with experimental measurements, our results are substantially higher than that determined by Zerr and Boehler, and the discrepancy between DAC and MD melting temperatures may be well explained with different melting mechanisms. Meanwhile, the radial distribution functions (RDFs) of Mg–Mg, O–Mg, and O–O ion pairs near the melting temperature have been investigated.     

Melting and Grüneisen parameters of NaCl at high pressure

The Buckingham potential has been employed to simulate the melting and thermodynamic parameters of sodium chloride (NaCl) using the molecular dynamics (MD) method. The constant-volume heat capacity and Grüneisen parameters have been obtained in a wide range of temperatures. The calculated thermodynamic parameters are found to be in good agreement with the available experimental data. The NaCl melting simulations appear to validate the interpretation of superheating of the solid in the one-phase MD simulations. The melting curve of NaCl is compared with the experiments and other calculations at pressure 0－30GPa range.     

Lattice dynamics and melting features of Li and Na

The high-pressure melting of Li and Na has been studied using ab initio calculations of the lattice dynamics. It has been shown that the recently discovered anomalous melting of Na is adequately explained by the phonon spectrum behavior and, accordingly, the thermal vibration amplitudes under compression. In a simple approach using the Lindemann criterion, the nonmonotonic behavior of the melting curve T _m(p) of Na has been quantitatively described within very wide pressure and temperature ranges, and, in particular, the melting temperature drop at p ～ 1 Mbar down to values lower than those at normal pressure. This approach leads to a nonphysical discontinuity of the melting curve T _m(p) of Li near the bcc-fcc-liquid triple point. This is due to the “softness” of the phonon spectrum of the bcc phase of Li that is the necessary condition for the existence of the high-temperature bcc phase. The melting of Na and Li is used as an example to determine why the Lindemann criterion is efficient in some situations and is inapplicable in the other cases.     

Atomic mechanism of homogeneous melting of bcc Fe at the limit of superheating

Atomic mechanism of homogeneous melting of bcc Fe is studied via monitoring spatiotemporal arrangements of the liquid-like atoms, which are detected by the Lindemann criterion of melting, during the heating process. Calculations are performed by molecular dynamics (MD) simulations. Calculations show that liquid-like atoms occur randomly in the crystalline matrix at temperature far below the melting point due to local instability of the crystalline lattice. Number of liquid-like atoms increases with increasing temperature and they have a tendency to form clusters. Subsequently, a single percolated liquid-like cluster is formed in the crystalline model and at the melting point 99% atoms in the model become liquid-like to form a liquid phase. Melting is also accompanied by the sudden changes in various static and thermodynamic quantities. However, total melting is reached just at the point above the melting one. Three characteristic temperatures of the homogeneous melting of bcc Fe are determined.     

Melting Curve of Iron: The Never-Ending Story?

For the better general understanding of melting behavior at high pressure, we investigated the influence of both crystallographic and electronic structure, and compressibility on melting temperatures for a large class of materials. In particular, we have established a large data base for melting of transition metals to megabar pressures. In general, bcc metals ( e.g. W, Ta, Mo, V, Cr) have very flat melting curves, and the initially steeper melting curves of fcc metals (Fe, Co, Ni) flatten significantly at high pressure. We also observed this trend for the more complicated alkaline earth, and rare earth metals. The flattening of the melting curves is due to the similarity of the solid and liquid structures and volumes. For iron there may exist an additional complication which may explain the reported results on both melting temperatures and structures. Due to the similarity in the free energies of its high pressure structures, these may coexist over a large pressure range. This phenomenon has been recently documented for noble gases with similar structures.     

高压下MgO熔化特性的分子动力学模拟

利用分子动力学方法,模拟研究了高压下MgO的熔化特性．通过晶体的现代熔化理论,对MgO的分子动力学模拟熔化温度进行了修正,得到了高温高压下MgO的熔化温度．计算得到的MgO熔化曲线和已有的实验及其它理论结果在0-135 GPa进行了比较,发现修正得到的MgO熔化温度和由Lindemann熔化方程及两相方法得到的结果在压力低于15 GPa时符合很好．同时,MgO熔化模拟有效解释了一阶相变分子动力学过程中出现的过热熔化现象．     

Melting and isothermal bulk modulus of the rocksalt phase of ZnO with molecular dynamics simulation

The Buckingham potential has been employed to simulate the melting of the rocksalt phase of ZnO over a wide pressure–temperature range and isothermal bulk modulus with shell model molecular dynamics method. The melting curve from molecular dynamics simulations is in good agreement with the results obtained from Lindemann melting equation in the pressure below 7 GPa. The calculated null compression bulk modulus is compared with the available experimental data and other theoretical results. At extended temperature ranges, the isothermal bulk modulus has also been predicted up to 3000 K.     

Melting temperature of lead and sodium at high pressures

We have redetermined the melting temperature of lead and sodium as a function of pressure in a new pressure cell made up of low strength materials. Many of the anomalies reported in a prior determination of the melting curve of lead have largely disappeared. The slope of our new melting curve for lead is in close agreement with the slope computed from thermochemical data.One new curve for sodium differs only slightly from the prior published curves.     

Self‐consistent theory of dynamic melting of a vortex lattice

The dynamic melting of vortex lattices in type II superconductors is considered. A field‐theoretic formulation of the pinning problem allows the average over the quenched disorder to be performed exactly. A self‐consistent theory is constructed using a functional method for the effective action, allowing a determination of the pinning force and the vortex fluctuations. The phase diagram for the dynamic melting transition is determined numerically. In contrast to perturbation theory, the self‐consistent theory is in quantitative agreement with the prediction of a recent phenomenological theory and simulations and experimental data.     

Molecular dynamics simulations of He bubble nucleation at grain boundaries

The nucleation behavior of He bubbles in single-crystal (sc) and nano-grain body-centered-cubic (bcc) Mo is simulated using molecular dynamics (MD) simulations, focusing on the effects of the grain boundary (GB) structure. In sc Mo, the nucleation behavior of He bubbles depends on irradiation conditions. He bubbles nucleate by either clustering of He atoms with pre-existing vacancies or self-interstitial-atom (SIA) punching without initial vacancies. In nano-grain Mo, strong precipitation of He at the GBs is observed, and the density, size and spatial distribution of He bubbles vary with the GB structure. The corresponding He bubble density is higher in nano-grain Mo than that in sc Mo and the average bubble size is smaller. In the GB plane, He bubbles distribute along the dislocation cores for GBs consisting of GB dislocations and randomly for those without distinguishable dislocation structures. The simulation results in nano-grain Mo are in agreement with previous experiments in metal nano-layers, and they are further explained by the effect of excess volume associated with the GBs.     

Phonon dispersion for the bcc metals Mo and Cr

In this paper the phonon dispersion for the bcc metals Mo and Cr is calculated based on the pair potentials obtained from cohesive energies and the Slater-Kirkwood- type three body interaction. In the calculation of the pair potentials the Möbius transform in the number theory is used and the cohesive energy is evaluated by the LMTO method. The results show a good agreement with inelastic neutron scattering data and indicate that the three-body interaction is necessary to account for the phonon dispersion.     

单轴应变条件下Fe从α到ε结构相变机制的第一性原理计算

采用基于密度泛函理论的平面波赝势方法,研究了沿［001］方向单轴应变条件下Fe从体心立方结构(bcc,α相)到六角密排结构(hcp,ε相)相变的临界压力、相变路径、相变势垒以及相变过程中原子磁性的变化.结果发现:单轴应变条件下Fe从α到ε结构的相变路径与以前理论计算模拟给出的静水压力条件下的相变路径明显不同;原子磁矩沿着相变路径突然降低,同时伴随着能量和体积的突然变化,是典型的一阶磁性相转变,表明原子磁性的丧失导致了bcc结构不稳定而向hcp结构转变.对单轴应变下吉布斯自由能的计算表明,相变势垒随着单轴应 关键词： 相变 单轴应变 第一性原理 铁     

The melting curve of CaSiO3 perovskite under lower mantle pressures

The high pressure melting curve of CaSiO₃ perovskite is simulated by using the constant temperature and pressure molecular dynamics method combined with effective pair potentials for the first time. The simulated results for the partial radial distribution function all compare well with experiment. The calculated equation of state is very successful in accurately reproducing the recent experimental data over a wide pressure range. The predicted high pressure melting curve is in good agreement with the experimental ones, and the melting curve up to the core–mantle boundary pressure, being very steep at lower pressures, rapidly flattens on increasing pressure. The present results also suggest the validity of the experimental data of Zerr and Boehler.