Shock Induced Order-disorder Transformation in Ni3Al

The Hugoniot of Ni3Al with L12 structure is calculated with an equation of state (EOS) based on a cluster expansion and variation method from first principles. It is found that an order-disorder transition occurs at a shock pressure of 205GPa, corresponding to 3750K in temperature. On the other hand, an unexpected high melting temperature about 6955K is obtained at the same pressure, which is completely different from the case at ambient pressure where the melting point is slightly lower than the order-disorder transition temperature, implying the high pressure phase diagram has its own characteristics. The present work also demonstrates the configurational contribution is more important than electronic excitations in alloys and mineral crystals within a large range of temperature, and an EOS model based on CVM is necessary for high pressure metallurgy and theoretical Earth model.

💡 Research Summary

This paper presents a first‑principles based equation of state (EOS) for the intermetallic compound Ni₃Al in its L1₂ ordered structure, incorporating configurational effects through a combined cluster expansion (CE) and Cluster Variation Method (CVM) framework. Density‑functional theory (DFT) calculations were performed for a series of atomic configurations ranging from perfectly ordered L1₂ to completely disordered solid solutions. The resulting total energies were fitted to a CE, yielding effective cluster interaction (ECI) parameters that capture pair, triplet, and quadruplet interactions.

The CE Hamiltonian was then fed into a CVM calculation using the tetrahedron (four‑point) approximation, which provides a statistical‑mechanical description of the configurational entropy and internal energy as functions of temperature and pressure. By minimizing the CVM free‑energy functional, equilibrium order parameters and configurational contributions to the thermodynamic potentials were obtained across a wide temperature range.

These CVM results were incorporated into an EOS that also includes the usual electronic and vibrational (phonon) contributions. The combined EOS was used to compute the Hugoniot curve—i.e., the locus of states satisfying the Rankine‑Hugoniot jump conditions—under shock compression up to 300 GPa.

The calculations reveal a sharp order‑disorder transition at a shock pressure of approximately 205 GPa, corresponding to a temperature of about 3750 K. At this point the L1₂ ordering collapses, the material expands abruptly, and the Hugoniot exhibits a distinct kink. Remarkably, the same pressure yields a melting temperature of roughly 6955 K, which is more than twice the order‑disorder temperature. This high‑pressure melting point is opposite to the ambient‑pressure behavior where the melting temperature lies just below the order‑disorder transition. The authors attribute the elevated melting point to the strong suppression of lattice vibrations under compression, which raises the thermal energy required for disordering the crystal lattice into a liquid.

A separate analysis of the individual contributions shows that, above ~1000 K, the configurational entropy associated with ordering dominates over both electronic excitations and phonon entropy. Consequently, conventional EOS models that neglect configurational degrees of freedom (e.g., Mie‑Grüneisen or Debye models) would significantly underestimate the thermodynamic response of Ni₃Al under extreme conditions.

The paper also discusses broader implications. In high‑pressure metallurgy, the ability to predict order‑disorder transitions and melting points accurately is crucial for designing shock‑resistant alloys and for interpreting dynamic compression experiments. In geophysics, the findings suggest that intermetallic phases similar to Ni₃Al could remain solid at temperatures far higher than previously thought in the deep Earth, potentially affecting models of mantle and core dynamics.

In summary, the study makes three major contributions: (1) it introduces a robust CE‑CVM‑based EOS that captures configurational effects in a quantitative manner; (2) it predicts a pressure‑induced order‑disorder transition at 205 GPa and an unexpectedly high melting temperature at the same pressure for Ni₃Al; and (3) it demonstrates that configurational contributions outweigh electronic and vibrational terms over a wide temperature range, underscoring the necessity of CVM‑type EOS for accurate high‑pressure material modeling.