Melting curves and entropy of melting of iron under Earth's core conditions

WJ Zhang and ZY Liu and ZL Liu and LC Cai, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 244, 69-77 (2015).

DOI: 10.1016/j.pepi.2014.10.011

The melting curves of iron are determined up to 365 GPa via molecular dynamic (MD) simulations combining with the embedded atom model (EAM) potential developed by Ackland et al. We simulated the melting with three approaches, the hysteresis, two-phase and recently modified Z methods. All three techniques can produce satisfying results, consistent well with most of static compression measurements and shock experiments. Hence, we recommend that these three techniques and this EAM potential are reliable techniques and potential for simulating melting properties of iron. Fitting the well-known Simon equation to our two-phase data we yield the analytical melting curve for iron: 1825(1 + P/57.723)(0.654), which gives a melting point at the inner core boundary of 6345 K, very close to the recent diamond anvil cell (DAC) extrapolated value and other ab initio calculations. Furthermore, the analyses of our entropy of melting and solid-liquid interfacial energy gamma(sl) indicate that at high pressure, the entropy of fusion shows weak pressure effect. The gamma(sl) increases monotonically with pressure, and can be described as a second-order polynomial relation. (C) 2014 Elsevier B.V. All rights reserved.

Return to Publications page