Pith. sign in

REVIEW 3 major objections 5 minor 62 references

Melting behavior of CaO at high temperature and pressure: a molecular dynamics study

T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read A molecular dynamics study computes the melting curve of calcium oxide up to 30 GPa without assuming a constant overheating ratio, and finds the ratio grows with pressure.

desk verdict Useful direct melting curve for CaO to 30 GPa, but internal number mismatches and a known short-range potential flaw need fixing before I'd trust the quantitative overheating ratios. read the letter →

arxiv 2412.06367 v1 pith:ENVFCCC5 submitted 2024-12-09 cond-mat.mtrl-sci physics.chem-ph

classification cond-mat.mtrl-sciphysics.chem-ph
keywords calciumoxidemeltingcurvemoleculardynamicstwo-phasecoexistenceoverheatingratioBorn-Meyer-Hugginspotentialhighpressureenthalpyoffusion
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Calcium oxide (CaO, lime) is a key component of planetary mantles and refractory ceramics, yet its melting point is experimentally uncertain by hundreds of kelvin because the material is reactive, volatile, and hard to contain at high temperature. This paper uses classical molecular dynamics to compute the melting temperature of CaO at ambient pressure with two independent techniques, obtaining 3066 ± 12 K by void-nucleated melting and 2940 ± 65 K by solid-liquid two-phase coexistence, and then extends the two-phase method to produce the first direct melting curve up to 30 GPa without assuming a constant Clapeyron slope or a constant overheating ratio. The central new result is that the thermal instability temperature Ts of the perfect crystal rises with pressure faster than the equilibrium melting temperature Tf, so the overheating ratio grows from about 19% at ambient pressure to about 32% at 30 GPa. That finding matters because previous high-pressure estimates of the CaO melting curve were obtained by scaling Ts by a fixed factor, an assumption the paper shows to be unjustified. If correct, the paper provides a direct computational melting curve for a geochemically central oxide and a cautionary result for similar scaling approaches.

What carries the argument

The argument is carried by classical molecular dynamics with the Born-Meyer-Huggins potential (ionic charges ±1.2e) combined with the two-phase solid-liquid coexistence technique. In that technique a supercell is half solid and half liquid, and after equilibration in the NPT ensemble the system is run in the isobaric-isenthalpic (NPH) ensemble: if the temperature drifts upward the liquid is recrystallizing, if it drifts downward the solid is melting, and the plateau temperature where neither happens is Tf. For the ambient-pressure cross-check the paper also uses void-nucleated melting, where a spherical cavity of increasing size is carved into the crystal until the apparent melting temperature plateaus. The high-pressure curve is obtained by repeating the coexistence simulation at 5 GPa intervals up to 30 GPa, with empirical fits $T_f(P) = 235.241 P^{0.558} + 2926.195$ K and $T_s(P) = 300.369 P^{0.618} + 3480.601$ K.

What would settle it

Measure the CaO melting temperature in a laser-heated diamond anvil cell at 10, 20, and 30 GPa; the paper predicts 3779 ± 45 K, 4224 ± 48 K, and 4444 ± 52 K. If the measured values fall outside roughly the quoted uncertainties, the potential-based melting curve is falsified; alternatively, a free-energy (thermodynamic integration) calculation with the same potential or an ab initio potential would show whether the slope is an artifact.

Watch

Extended reading notes

Core claim

The paper's central claim is that the melting curve of CaO can be computed directly, without empirical input about the Clapeyron slope or the overheating ratio, and that when this is done the melting curve deviates from the constant-scaling picture. At ambient pressure the two-phase coexistence method yields Tf = 2940 ± 65 K and the void-nucleated method yields 3066 ± 12 K, bracketing recent laser-heating measurements; the melting temperature then rises to 3432 ± 43 K at 5 GPa and 4444 ± 52 K at 30 GPa. The thermal instability temperature Ts, measured on defect-free crystals, rises from 3503 K to 5858 K over the same pressure range, so the overheating ratio η = Ts/Tf − 1 increases from 19.1% to 31.8%. This pressure dependence directly contradicts the assumption, used in the only previous high-pressure melting curve of CaO, that Ts can be scaled by a constant factor to obtain Tf. The paper therefore positions its two-phase coexistence calculation as the first direct determination of the CaO melting curve beyond ambient pressure.

Load-bearing premise

The load-bearing premise is that the Born-Meyer-Huggins potential with ionic charges ±1.2e stays accurate enough for computing melting at compressions up to 30 GPa, even though at very short Ca-O distances the potential becomes attractive and unphysical and can spoil the simulation; the paper itself treats 30 GPa as the upper limit for this reason.

Editorial extensions

If this is right

  • The CaO melting curve is available from 0 to 30 GPa as a direct simulation output, so thermodynamic assessments of CaO-bearing planetary and ceramic systems no longer need a Clapeyron slope borrowed from low-pressure experiments.
  • The overheating ratio of CaO rises from about 19% at ambient pressure to about 32% at 30 GPa, so scaling the thermal instability temperature by a fixed factor underestimates the melting temperature by an amount that grows with pressure.
  • The ambient-pressure agreement between the two-phase (2940 ± 65 K) and void-nucleated (3066 ± 12 K) techniques indicates the computed melting temperature is not an artifact of a single nucleation protocol.
  • The classical-MD enthalpy of fusion, 80.37 kJ/mol, matches the assessed tabulated value of 79.5 kJ/mol, supporting the thermodynamic consistency of the potential used for the melting curve.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: The same constant-scaling shortcut that this paper invalidates for CaO is used elsewhere for other refractory oxides, so if the pressure dependence of the overheating ratio is generic, those melting curves may need revisiting.
  • Editorial inference: Because the paper's own potential is acknowledged to become unphysical at very short Ca-O distances, an immediate test would be to recompute the 30 GPa point with a potential that adds a steep repulsive wall and see whether the predicted flattening of the curve survives.
  • Editorial inference: The ambient-pressure bracket of 2940–3066 K could be sharpened with longer isobaric-isenthalpic runs, which would also separate the two recent laser-heating measurements more decisively.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This manuscript reports classical and ab initio molecular dynamics simulations of the melting behavior of CaO. At ambient pressure, the void-nucleated melting technique gives Tf = 3066 ± 12 K and the two-phase coexistence technique gives Tf = 2940 ± 65 K. Using classical MD with the two-phase method, the authors compute the high-pressure melting curve and the thermal instability curve up to 30 GPa, and they report that the overheating ratio η = Ts/Tf − 1 increases with pressure, from 19.1% at ambient pressure to 31.8% at 30 GPa in Table 3. This is used to argue that assuming a constant overheating ratio, as in earlier work by Sun et al., is not justified. The manuscript also reports an enthalpy of fusion of 80.37 kJ/mol from classical MD, in good agreement with tabulated values and with prior simulations.

Significance. If the reported results are correct, the paper provides a useful direct melting curve for CaO up to 30 GPa without imposing a constant Clapeyron slope or a constant overheating ratio, and it identifies a pressure-dependent overheating ratio that bears on extrapolation strategies for refractory oxides. The paper has notable strengths: it compares two independent melting techniques, validates the caloric curve against NIST-JANAF and earlier simulations, and combines classical and ab initio MD. The enthalpy of fusion result is robust and well aligned with literature values. However, the central quantitative claims currently rest on a numerical inconsistency between the text and Table 3, and on a classical potential whose short-range unphysical attraction is acknowledged by the authors; these issues need to be resolved before the main conclusions can be fully accepted.

major comments (3)
  1. [Sec. 3.3 and Table 3] The text states that the overheating ratio increases 'from about 16% at ambient pressure to more than 24% at P=30 GPa,' but Table 3 lists 19.1% at 0 GPa and 31.8% at 30 GPa for the same quantities. Since the pressure dependence of η is the main physical conclusion of the paper, this discrepancy must be reconciled. In addition, Eq. (2) does not reproduce the tabulated values at P = 0: the fit gives Tf = 2926.2 K and Ts = 3480.6 K, whereas Table 3 reports 2940 K and 3503 K. Please provide the residuals or goodness-of-fit for Eq. (2) and explain the anchoring choice.
  2. [Sec. 2.1 and Sec. 3.3] The BMH potential is described as 'attractive at very short distances for the cation-anion pairs,' with an 'unphysical' attractive well, and the paper states that the 30 GPa ceiling exists because of the 'lack of a short-range repulsive barrier.' Since both the melting curve Tf(P) and the thermal instability curve Ts(P) are computed with this potential, a short-range artifact that is sampled differently by the solid and liquid under compression could bias both quantities and hence the reported increase of η with pressure. The potential was validated for liquid structure and enthalpy at ambient pressure, but no high-pressure validation is provided. Please quantify the spurious close Ca–O contact populations at high pressure, or repeat key points (e.g., 15 and 30 GPa) with an added repulsive short-range term or with ab initio MD, and report whether the η(P) trend is preserved.
  3. [Sec. 3.1.1, Sec. 3.1.2, and Sec. 3.3] The two ambient-pressure melting temperatures, 3066 ± 12 K (void-nucleated) and 2940 ± 65 K (two-phase), differ by 126 K, which is roughly 1.9 combined standard errors. Calling these values 'consistent' is too strong without a quantitative reconciliation, especially because the paper chooses 2940 K as the zero-pressure anchor for the high-pressure melting curve while using 3066 K for the enthalpy-of-fusion calculation and for comparison with laser-heating experiments. An offset of this size propagates directly into the reported Tf(P) curve. Please justify the choice of the two-phase value as the anchor and discuss the systematic error this introduces.
minor comments (5)
  1. [Abstract] The word 'Claperyon' appears in the abstract and should be 'Clapeyron.'
  2. [Sec. 3.1.2] The phrase 'Contrary to what is reported in Ref.,13' contains a stray comma before the reference number; please correct the citation format.
  3. [Sec. 3.3] The sentence 'This can be explained by empirical laws described in which have the same form of the Equations 2 and 3 above' is incomplete; a reference or derivation appears to be missing.
  4. [Sec. 3.2 and Table 2] Table 2 lists the AIMD values at 3000 K and 3200 K under 'enthalpy of fusion,' while the text calls these values 'solution enthalpy' at temperatures other than the melting temperature; please clarify the distinction and relabel the table entries accordingly.
  5. [Sec. 3.1.2] Please specify how the uncertainty of ±65 K for the two-phase Tf was obtained; the description of the NPH runs does not indicate whether this is a standard deviation of time averages, a fit uncertainty, or an estimate of finite-size effects.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the melting temperatures and high-pressure curves come from direct coexistence and void-nucleated MD simulations, not from fitted inputs; Eq. (2) is a descriptive fit to the authors' own data, and the self-citations are peripheral.

full rationale

The paper's central results, Tf at ambient pressure and the Tf(P) and Ts(P) curves, are produced by direct molecular dynamics simulations using the two-phase coexistence and void-nucleated techniques. Neither method fits a parameter to the quantity being predicted: the interatomic potential is taken from prior work by Alvares et al. and is validated against structural, dynamic, and thermodynamic data, not against the CaO melting point. The ambient-pressure Tf values are cross-checked against experimental laser-heating measurements and against independent MD results, and the enthalpy of fusion is computed from the caloric curves at the simulated Tf. The high-pressure melting curve is obtained by running the same coexistence protocol at successive pressures, and the thermal instability curve Ts is obtained by heating a defect-free crystal. The conclusion that Ts/Tf increases with pressure rests on the directly simulated values in Table 3, not on any fitted relation. Equation (2) is a post-hoc empirical fit to those computed points and is presented only for compact representation, so it does not constitute a prediction from fitted inputs. The self-citations (Refs. 2, 59-61) appear in contextual or uncertainty discussions and are not load-bearing for the main derivation. The unphysical short-range attraction of the BMH potential, acknowledged by the authors, is a potential source of systematic error in the high-pressure results, but that is an accuracy concern, not a circularity of the derivation chain.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The central results rest primarily on the BMH potential parameters from prior work and on the assumptions about MD protocols. The two-phase coexistence melting curve uses the BMH potential, which is fitted to Ca-O interactions in earlier studies and has a known unphysical short-range attraction. The thermal instability curve uses rapid heating of defect-free crystals, a kinetic procedure. No new physical entities are introduced. The empirical fits of Eq. (2) are post-hoc parametrizations of the computed data points.

free parameters (2)
  • BMH potential parameters (Aij, rho_ij, sigma_ij, Cij, Dij) for Ca-Ca, Ca-O, O-O; ionic charges q(Ca)=+1.2e, q(O)=-1.2e = Listed in Table 1 (Aij in kcal/mol, rho_ij and sigma_ij in Angstrom, Cij in kcal/mol Angstrom^6, Dij=0)
    These parameters come from prior work (Refs 12 and 35). They are fitted quantities that entirely determine the computed melting temperatures and the 30 GPa melting curve; they are not re-derived in this paper.
  • Empirical fit coefficients in Eq. (2) = Tf: 235.241, 0.558, 2926.195; Ts: 300.369, 0.618, 3480.601
    Coefficients are least-squares fits to the authors' computed data points (Table 3). They are used to present the curves and to quantify the slope comparison, but they carry no error bars and are not predictions.
assumptions (5)
  • domain assumption The BMH potential, including its short-range behavior, accurately describes CaO energetics at pressures up to 30 GPa and temperatures up to roughly 6000 K.
    Invoked in Sec. 2.1 and Sec. 3.3; the paper itself notes the potential has an unphysical attractive well at short Ca-O distances and no steep repulsive barrier, which limits simulations to 30 GPa and may bias the high-pressure melting curve.
  • domain assumption The two-phase coexistence NPT/NPH protocol yields the equilibrium melting temperature within the quoted uncertainty.
    Invoked in Sec. 3.1.2; the method assumes interface motion and finite-size effects are negligible and that the Berendsen barostat with fixed in-plane lattice parameter is appropriate.
  • domain assumption The thermal instability temperature Ts measured by heating a defect-free crystal at ~10^12 K/s corresponds to the true mechanical melting limit.
    Invoked in Sec. 3.3; Ts depends on heating rate and simulation time, so its use for the overheating ratio requires this kinetic quantity to be a meaningful thermodynamic limit.
  • domain assumption PBEsol is sufficiently accurate for the enthalpies of solid and liquid CaO.
    Invoked in Sec. 2.2 for AIMD; the paper notes AIMD underestimates CMD enthalpy by up to 15% near melting, so this assumption affects only the AIMD enthalpy values, not the central melting-curve claim.
  • standard math Statistical mechanics of the NPT/NPH ensembles, Ewald summation, and Nosé-Hoover thermostats provide unbiased averages over 100-300 ps trajectories.
    Standard MD background assumed throughout Sec. 2; accepted practice, though the short trajectories (100-250 ps) and high heating rates are not independently validated for convergence.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Melting behavior of CaO at high temperature and pressure: a molecular dynamics study." pith.science (2026). https://pith.science/paper/ENVFCCC5

@misc{pith2026241206367,
  author       = {Pith},
  title        = {Pith review of: Melting behavior of CaO at high temperature and pressure: a molecular dynamics study},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ENVFCCC5}},
  note         = {Machine review of arXiv:2412.06367}
}
abstract

The thermodynamic behavior of calcium oxide (\ce{CaO}) under high temperature and pressure conditions is critical for understanding the physics of planetary interiors. This study employs molecular dynamics (MD) simulations, including both classical and ab-initio approaches, to investigate the melting behavior of CaO. We calculate the melting temperature of \ce{CaO} by the void-nucleated melting and two-phase coexistence techniques, aiming to resolve discrepancies in experimental data on the melting point, which range from 2843~K to 3223~K in different studies due to the high reactivity and vapor pressure of the substance. The obtained results are $T_f = 3066\pm12$~K and $T_f = 2940\pm65$~K using the void-nucleated melting and the two-phase coexistence method, respectively. Additionally, we calculate the enthalpy of fusion and the high-pressure melting curve, for the first time without making any assumption on the Clapeyron slope. This is extremely important since in experiments the Claperyon slope of the melting curve is estimated from low pressure measurements and the overheating ratio (i.e. $\eta=\frac{T_s}{T_f}-1$, where $T_s$ represents the thermal instability limit corresponding to the homogeneous melting temperature of the solid) is often assumed to be constant in simulations. Our MD results show that $T_s$ increases more rapidly with pressure than $T_f$ and thus that the overheating ratio sensibly depends upon pressure. These findings contribute to the accurate modeling of the CaO phase diagram, which is essential for geochemistry, cosmochemistry, and materials science.

Figures

Figures reproduced from arXiv: 2412.06367 by the authors.

Figure 1
Figure 1. Apparent melting temperature of crystal CaO as a function of the defect volume (in [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Temperature as a function of simulation time in a 6 [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Snapshot of the simulation box with coexisting solid and liquid phases, visualized [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Enthalpy change (relative to the standard-state enthalpy of the solid, i.e. [PITH_FULL_IMAGE:figures/full_fig_p017_4.png]
Figure 5
Figure 5. Figure 5: Relative volume thermal expansion as a function of temperature. The lower branch [PITH_FULL_IMAGE:figures/full_fig_p019_5.png]
Figure 6
Figure 6. Figure 6: High-pressure melting curve of CaO (Tf , black squares) and thermal instability temperature curve (Ts, black dots) as calculated in this work. The melting curve estimated by Sun et al.15 is shown for comparison (in blue). Empty symbols represent different values of Tf …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

62 extracted references · 42 canonical work pages

  1. [1]

    Anderson, D. L. Theory of the Earth; Blackwell Scientific Publications, 1989

  2. [2]

    Belmonte, D.; Ottonello, G.; Zuccolini, M. V. Ab initio-assisted assessment of the CaO - SiO2 system under pressure. Calphad 2017, 59, 12–30, DOI: doi:10.1016/j.calphad.2017.07.009

  3. [3]

    Cement Chemistry; Emerald Publishing Limited, 1997

    Taylor, H. Cement Chemistry; Emerald Publishing Limited, 1997

  4. [4]

    Manara, D.; Böhler, R.; Capriotti, L.; Quaini, A.; Bao, Z.; Boboridis, K.; Luzzi, L.; Janssen, A.; Pöml, P.; Eloirdi, R. et al. On the melting behaviour of calcium monoxide under different atmospheres: A laser heating study. Journal of the European Ceramic Society 2014, 34, 1623–1636, DOI: doi:10.1016/j.jeurceramsoc.2013.12.018

  5. [5]

    Complete thermodynamic description of the Mg-Ca-O phase diagram including the Ca-O , Mg-O and CaO-MgO subsystems

    Liang, S.-M.; Schmid-Fetzer, R. Complete thermodynamic description of the Mg-Ca-O phase diagram including the Ca-O , Mg-O and CaO-MgO subsystems. Journal of the European Ceramic Society 2018, 38, 4768–4785, DOI: doi:10.1016/j.jeurceramsoc.2018.06.015

  6. [6]

    Measurement of the solidification points of several refractory oxides by means of a solar furnace

    Foex, M. Measurement of the solidification points of several refractory oxides by means of a solar furnace. Solar Energy 1965, 9, 61–67, DOI: doi:10.1016/0038-092x(65)90162-3

  7. [7]

    Reinvestigation of the Solidification Point of CaO by Digital Pyrometry

    Yamada, T.; Yoshimura, M.; Sōmiya, S. Reinvestigation of the Solidification Point of CaO by Digital Pyrometry. Journal of the American Ceramic Society 1986, 69, C--243--C--245, DOI: doi:10.1111/j.1151-2916.1986.tb07350.x

  8. [8]

    Laser‐pulse melting of calcium oxide and some peculiarities of its high‐temperature behavior

    Bgasheva, T.; Falyakhov, T.; Petukhov, S.; Sheindlin, M.; Vasin, A.; Vervikishko, P. Laser‐pulse melting of calcium oxide and some peculiarities of its high‐temperature behavior. Journal of the American Ceramic Society 2021, 104, 3461–3477, DOI: doi:10.1111/jace.17676

Show all 62 references
  1. [9]

    Molecular Dynamics Simulation of the Thermodynamic and Structural Properties for the CaO-SiO2 System

    Seo, W.-G.; Tsukihashi, F. Molecular Dynamics Simulation of the Thermodynamic and Structural Properties for the CaO-SiO2 System. ISIJ International 2004, 44, 1817–1825, DOI: doi:10.2355/isijinternational.44.1817

  2. [10]

    B.; Skorodumova, N

    Belonoshko, A. B.; Skorodumova, N. V.; Rosengren, A.; Johansson, B. Melting and critical superheating. Physical Review B 2006, 73, 012201, DOI: doi:10.1103/physrevb.73.012201

  3. [11]

    Theory of Melting

    Kuhlmann-Wilsdorf, D. Theory of Melting. Phys. Rev. 1965, 140, A1599--A1610, DOI: doi:10.1103/PhysRev.140.A1599

  4. [12]

    Alvares, C. M. S.; Deffrennes, G.; Pisch, A.; Jakse, N. Thermodynamics and structural properties of CaO : A molecular dynamics simulation study. The Journal of Chemical Physics 2020, 152, 084503, DOI: doi:10.1063/1.5141841

  5. [13]

    Prediction of the melting curve and phase diagram for CaO using newly developed interatomic potentials

    Wang, X.-W.; Sun, X.-W.; Song, T.; Tian, J.-H.; Liu, Z.-J. Prediction of the melting curve and phase diagram for CaO using newly developed interatomic potentials. Vacuum 2023, 209, 111717, DOI: doi:10.1016/j.vacuum.2022.111717

  6. [14]

    Ab initio construction of full phase diagram of MgO - CaO eutectic system using neural network interatomic potentials

    Lee, K.; Park, Y.; Han, S. Ab initio construction of full phase diagram of MgO - CaO eutectic system using neural network interatomic potentials. Physical Review Materials 2022, 6, 113802, DOI: doi:10.1103/physrevmaterials.6.113802

  7. [15]

    The high-pressure melting curve of CaO

    Sun, X.; Song, T.; Chu, Y.; Liu, Z.; Zhang, Z.; Chen, Q. The high-pressure melting curve of CaO . Solid State Communications 2010, 150, 1785–1788, DOI: doi:10.1016/j.ssc.2010.07.015

  8. [16]

    NIST-JANAF Thermochemical Tables, 4th Edition; American Institute of Physics, 1998

    M, C. NIST-JANAF Thermochemical Tables, 4th Edition; American Institute of Physics, 1998

  9. [17]

    Sun, J.; Ruzsinszky, A.; Perdew, J. P. Strongly Constrained and Appropriately Normed Semilocal Density Functional. Physical Review Letters 2015, 115, 036402, DOI: doi:10.1103/physrevlett.115.036402

  10. [18]

    Thermodynamic and phase diagram data for the CaO-SiO2 system

    Taylor, J.; Dinsdale, A. Thermodynamic and phase diagram data for the CaO-SiO2 system. Calphad 1990, 14, 71–88, DOI: doi:10.1016/0364-5916(90)90041-w

  11. [19]

    An assessment of the CaO - SiO2 system

    Hillert, M.; Sundman, B.; Wang, X. An assessment of the CaO - SiO2 system. Metallurgical Transactions B 1990, 21, 303–312, DOI: doi:10.1007/bf02664198

  12. [20]

    Eriksson, G.; Pelton, A. D. Critical evaluation and optimization of the thermodynamic properties and phase diagrams of the CaO - Al2O3 , Al2O3 - SiO2 , and CaO - Al2O3 - SiO2 systems. Metallurgical Transactions B 1993, 24, 807–816, DOI: doi:10.1007/bf02663141

  13. [21]

    Hudon, P.; Jung, I.-H.; Baker, D. R. Experimental Investigation and Optimization of Thermodynamic Properties and Phase Diagrams in the Systems CaO - SiO2 , MgO - SiO2 , CaMgSi2O6 - SiO2 and CaMgSi2O6 - Mg2SiO4 to 10 GPa . Journal of Petrology 2005, 46, 1859–1880, DOI: doi:10.1...

  14. [22]

    Y.; Reinhardt, A

    Chew, P. Y.; Reinhardt, A. Phase diagrams - Why they matter and how to predict them. The Journal of Chemical Physics 2023, 158, 030902, DOI: doi:10.1063/5.0131028

  15. [23]

    Accurate Melting Temperatures for Neon and Argon from Ab Initio Monte Carlo Simulations

    Pahl, E.; Calvo, F.; Koči, L.; Schwerdtfeger, P. Accurate Melting Temperatures for Neon and Argon from Ab Initio Monte Carlo Simulations. Angewandte Chemie International Edition 2008, 47, 8207--8210, DOI: doi:https://doi.org/10.1002/anie.200802743

  16. [24]

    Alfè, D.; Cazorla, C.; Gillan, M. J. The kinetics of homogeneous melting beyond the limit of superheating. The Journal of Chemical Physics 2011, 135, 024102, DOI: doi:10.1063/1.3605601

  17. [25]

    Modified Z method to calculate melting curve by molecular dynamics

    Wang, S.; Zhang, G.; Liu, H.; Song, H. Modified Z method to calculate melting curve by molecular dynamics. The Journal of Chemical Physics 2013, 138, 094114, DOI: doi:10.1063/1.4798225

  18. [26]

    R.; Hummel, F.; Kresse, G.; Kahl, G.; Dellago, C

    Pedersen, U. R.; Hummel, F.; Kresse, G.; Kahl, G.; Dellago, C. Computing Gibbs free energy differences by interface pinning. Phys. Rev. B 2013, 88, 094101, DOI: doi:10.1103/PhysRevB.88.094101

  19. [27]

    Kirkwood, J. G. Statistical Mechanics of Fluid Mixtures. The Journal of Chemical Physics 1935, 3, 300–313, DOI: doi:10.1063/1.1749657

  20. [28]

    Understanding Molecular Simulation: From Algorithms to Applications, 2nd ed.; Computational Science Series; Academic Press: San Diego, 2002; Vol

    Frenkel, D.; Smit, B. Understanding Molecular Simulation: From Algorithms to Applications, 2nd ed.; Computational Science Series; Academic Press: San Diego, 2002; Vol. 1

  21. [29]

    Alavi, S.; Thompson, D. L. Molecular Dynamics Simulations of the Melting of Aluminum Nanoparticles. The Journal of Physical Chemistry A 2005, 110, 1518–1523, DOI: doi:10.1021/jp053318s

  22. [30]

    Zhang, Y.; Maginn, E. J. A comparison of methods for melting point calculation using molecular dynamics simulations. The Journal of Chemical Physics 2012, 136, 144116, DOI: doi:10.1063/1.3702587

  23. [31]

    B.; Lukinov, T.; Burakovsky, L.; Preston, D

    Belonoshko, A. B.; Lukinov, T.; Burakovsky, L.; Preston, D. L.; Rosengren, A. Melting of a polycrystalline material: Melting of real materials. The European Physical Journal Special Topics 2013, 216, 199–204, DOI: doi:10.1140/epjst/e2013-01743-1

  24. [32]

    Melting temperature prediction via first principles and deep learning

    Hong, Q.-J. Melting temperature prediction via first principles and deep learning. Computational Materials Science 2022, 214, 111684, DOI: doi:10.1016/j.commatsci.2022.111684

  25. [33]

    Brodholt, J

    Di Paola, C.; P. Brodholt, J. Modeling the melting of multicomponent systems: the case of MgSiO3 perovskite under lower mantle conditions. Scientific Reports 2016, 6, 29830, DOI: doi:10.1038/srep29830

  26. [34]

    L.; Mayer, J

    Huggins, M. L.; Mayer, J. E. Interatomic Distances in Crystals of the Alkali Halides. The Journal of Chemical Physics 1933, 1, 643–646, DOI: doi:10.1063/1.1749344

  27. [35]

    Structural and dynamic properties of calcium aluminosilicate melts: A molecular dynamics study

    Bouhadja, M.; Jakse, N.; Pasturel, A. Structural and dynamic properties of calcium aluminosilicate melts: A molecular dynamics study. The Journal of Chemical Physics 2013, 138, 224510, DOI: doi:10.1063/1.4809523

  28. [36]

    Fast Parallel Algorithms for Short-Range Molecular Dynamics

    Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. Journal of Computational Physics 1995, 117, 1–19, DOI: doi:10.1006/jcph.1995.1039

  29. [37]

    P.; Aktulga, H

    Thompson, A. P.; Aktulga, H. M.; Berger, R.; Bolintineanu, D. S.; Brown, W. M.; Crozier, P. S.; in ’t Veld, P. J.; Kohlmeyer, A.; Moore, S. G.; Nguyen, T. D. et al. LAMMPS -- a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum sc...

  30. [38]

    Hoover, W. G. Canonical dynamics: Equilibrium phase-space distributions. Physical Review A 1985, 31, 1695–1697, DOI: doi:10.1103/physreva.31.1695

  31. [39]

    L.; Cococcioni, M.; Dabo, I

    Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I. et al. Quantum ESPRESSO : a modular and open-source software project for quantum simulations of materials. Journal of Physics: Condensed Matter...

  32. [40]

    Giannozzi, P.; Andreussi, O.; Brumme, T.; Bunau, O.; Buongiorno Nardelli, M.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Cococcioni, M. et al. Advanced capabilities for materials modelling with Quantum ESPRESSO . Journal of Physics: Condensed Matter 2017, 29, 465901, ...

  33. [41]

    Giannozzi, P.; Baseggio, O.; Bonfà, P.; Brunato, D.; Car, R.; Carnimeo, I.; Cavazzoni, C.; de Gironcoli, S.; Delugas, P.; Ferrari Ruffino, F. et al. Quantum ESPRESSO toward the exascale. The Journal of Chemical Physics 2020, 152, 154105, DOI: doi:10.1063/5.0005082

  34. [42]

    Troullier, N.; Martins, J. L. Efficient pseudopotentials for plane-wave calculations. Physical Review B 1991, 43, 1993–2006, DOI: doi:10.1103/physrevb.43.1993

  35. [43]

    P.; Ruzsinszky, A.; Csonka, G

    Perdew, J. P.; Ruzsinszky, A.; Csonka, G. I.; Vydrov, O. A.; Scuseria, G. E.; Constantin, L. A.; Zhou, X.; Burke, K. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. Physical Review Letters 2008, 100, 136406, DOI: doi:10.1103/physrevlett.100.136406

  36. [44]

    I.; Perdew, J

    Csonka, G. I.; Perdew, J. P.; Ruzsinszky, A.; Philipsen, P. H. T.; Lebègue, S.; Paier, J.; Vydrov, O. A.; Ángyán, J. G. Assessing the performance of recent density functionals for bulk solids. Physical Review B 2009, 79, 155107, DOI: doi:10.1103/physrevb.79.155107

  37. [45]

    H.; Gumbsch, P.; Lu, K.; Ma, E

    Jin, Z. H.; Gumbsch, P.; Lu, K.; Ma, E. Melting Mechanisms at the Limit of Superheating. Physical Review Letters 2001, 87, 055703, DOI: doi:10.1103/physrevlett.87.055703

  38. [46]

    F.; Wolf, D.; Phillpot, S

    Lutsko, J. F.; Wolf, D.; Phillpot, S. R.; Yip, S. Molecular-dynamics study of lattice-defect-nucleated melting in metals using an embedded-atom-method potential. Physical Review B 1989, 40, 2841–2855, DOI: doi:10.1103/physrevb.40.2841

  39. [47]

    J.; Steinebrunner, G.; Kirchner, B.; Huber, H

    Solca, J.; Dyson, A. J.; Steinebrunner, G.; Kirchner, B.; Huber, H. Melting curve for argon calculated from pure theory. Chemical Physics 1997, 224, 253–261, DOI: doi:10.1016/s0301-0104(97)00317-0

  40. [48]

    J.; Steinebrunner, G.; Kirchner, B.; Huber, H

    Solca, J.; Dyson, A. J.; Steinebrunner, G.; Kirchner, B.; Huber, H. Melting curves for neon calculated from pure theory. The Journal of Chemical Physics 1998, 108, 4107–4111, DOI: doi:10.1063/1.475808

  41. [49]

    M.; Rice, B

    Agrawal, P. M.; Rice, B. M.; Thompson, D. L. Molecular dynamics study of the melting of nitromethane. The Journal of Chemical Physics 2003, 119, 9617–9627, DOI: doi:10.1063/1.1612915

  42. [50]

    Molecular-dynamics study of liquid nickel above and below the melting point

    Jakse, N.; Pasturel, A. Molecular-dynamics study of liquid nickel above and below the melting point. The Journal of Chemical Physics 2005, 123, 244512, DOI: doi:10.1063/1.2145759

  43. [51]

    R.; Wang, C

    Morris, J. R.; Wang, C. Z.; Ho, K. M.; Chan, C. T. Melting line of aluminum from simulations of coexisting phases. Physical Review B 1994, 49, 3109–3115, DOI: doi:10.1103/physrevb.49.3109

  44. [52]

    Solid-liquid coexistence in small systems: A statistical method to calculate melting temperatures

    Hong, Q.-J.; van de Walle, A. Solid-liquid coexistence in small systems: A statistical method to calculate melting temperatures. The Journal of Chemical Physics 2013, 139, DOI: doi:10.1063/1.4819792

  45. [53]

    Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. Molecular dynamics with coupling to an external bath. The Journal of Chemical Physics 1984, 81, 3684–3690, DOI: doi:10.1063/1.448118

  46. [54]

    Physics of solid and liquid alkali halide surfaces near the melting point

    Zykova-Timan, T.; Ceresoli, D.; Tartaglino, U.; Tosatti, E. Physics of solid and liquid alkali halide surfaces near the melting point. The Journal of Chemical Physics 2005, 123, 164701, DOI: doi:10.1063/1.2035096

  47. [55]

    Why Are Alkali Halide Surfaces Not Wetted by Their Own Melt? Physical Review Letters 2005, 94, 176105, DOI: doi:10.1103/physrevlett.94.176105

    Zykova-Timan, T.; Ceresoli, D.; Tartaglino, U.; Tosatti, E. Why Are Alkali Halide Surfaces Not Wetted by Their Own Melt? Physical Review Letters 2005, 94, 176105, DOI: doi:10.1103/physrevlett.94.176105

  48. [56]

    VMD : Visual molecular dynamics

    Humphrey, W.; Dalke, A.; Schulten, K. VMD : Visual molecular dynamics. Journal of Molecular Graphics 1996, 14, 33–38, DOI: doi:10.1016/0263-7855(96)00018-5

  49. [57]

    a rmen von CaO , Al2O3 , CaO.Al2O3 , 3CaO.Al2O3 , 2CaO.SiO2 , 3CaO.SiO2 , 2CaO.Al2O3.SiO2 von 20° bis 1500° C . Zeitschrift f \

    von Gronow, H. E.; Schwiete, H. Die spezifischen W \"a rmen von CaO , Al2O3 , CaO.Al2O3 , 3CaO.Al2O3 , 2CaO.SiO2 , 3CaO.SiO2 , 2CaO.Al2O3.SiO2 von 20° bis 1500° C . Zeitschrift f \"u r anorganische und allgemeine Chemie 1933, 216, 185--195

  50. [59]

    V.; Natali, M

    Ottonello, G.; Attene, M.; Ameglio, D.; Belmonte, D.; Zuccolini, M. V.; Natali, M. Thermodynamic investigation of the CaO-Al2O3-SiO2 system at high P and T through polymer chemistry and convex-hull techniques. Chemical Geology 2013, 346, 81--92, DOI: doi:https://doi.org/10.101...

  51. [60]

    Belmonte, D.; Ottonello, G.; Zuccolini, M. V. Melting of - Al2O3 and vitrification of the undercooled alumina liquid: Ab initio vibrational calculations and their thermodynamic implications . The Journal of Chemical Physics 2013, 138, 064507, DOI: doi:10.1063/1.4790612

  52. [61]

    V.; Attene, M

    Belmonte, D.; Ottonello, G.; Zuccolini, M. V.; Attene, M. The system MgO - Al2O3 - SiO2 under pressure: A computational study of melting relations and phase diagrams. Chemical Geology 2017, 461, 54–64, DOI: doi:10.1016/j.chemgeo.2016.11.011

  53. [62]

    High-temperature thermal expansion of lime, periclase, corundum and spinel

    Fiquet, G.; Richet, P.; Montagnac, G. High-temperature thermal expansion of lime, periclase, corundum and spinel. Physics and Chemistry of Minerals 1999, 27, 103–111, DOI: doi:10.1007/s002690050246

  54. [63]

    Melting of - Al2O3 and vitrification of the undercooled alumina liquid: Ab initio vibrational calculations and their thermodynamic implications

    Gómez, L.; Gazza, C.; Dacharry, H.; Peñaranda, L.; Dobry, A. Pressure dependence of the melting mechanism at the limit of overheating in Lennard-Jones crystals. Physical Review B 2005, 71, 134106, DOI: doi:10.1103/physrevb.71.134106 mcitethebibliography main.out000066400000000...

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.