Theory of melting lines with a variable enthalpy of fusion

This paper presents a new analytical model for melting lines by modifying the Clausius-Clapeyron relation to account for a variable enthalpy of fusion driven by anharmonic solid-state effects, yielding parabolic solutions defined by fundamental thermophysical properties that corroborate recent universal models of liquid melting.

The Gist

Imagine you are trying to draw a map of when a solid turns into a liquid (like ice melting into water). In physics, this map is called a melting line. It shows how much heat you need to melt something, depending on how much you squeeze it (pressure).

For a long time, scientists used a simple rulebook (the Clausius-Clapeyron equation) to draw these maps. But there was a catch: they assumed that the "energy cost" to melt a substance (called latent heat) never changed, no matter how hot or how squeezed it got. This worked great for things turning into gas (like boiling water), but it was a terrible guess for solids turning into liquids. Solids and liquids are both dense and sticky, so the rules are much more complicated.

The New Idea: The "Stretchy" Energy Cost
This paper proposes a new way to draw that map. The author, Anthony Papathanassiou, suggests that the energy cost to melt a solid isn't a fixed number; it's more like a stretchy rubber band. As you heat the solid, the atoms inside start to jiggle wildly (anharmonicity), and the amount of energy needed to break them loose changes depending on how much space the atoms have (volume).

Think of it like this:

• Old View: Imagine trying to push a heavy box up a ramp. You assume the weight of the box stays exactly the same the whole time.
• New View: The box is actually made of a special material that gets lighter or heavier depending on how fast you are moving it and how much you squeeze it. To get the right answer, you have to account for that changing weight.

The "Volume" Connection
The paper uses a clever trick. It looks at how much a solid expands when it gets hot (thermal expansion) and how much heat it holds. It turns out that near the melting point, the "stretchy" part of the heat capacity is directly linked to the difference in size between the solid and the liquid.

By plugging this "stretchy" energy idea into the old rulebook, the author derives a new mathematical equation.

The Result: A Perfect Parabola
When the author solves this new equation, the shape of the melting line isn't a straight line or a weird squiggle. It turns out to be a parabola (the same U-shape you see when you toss a ball in the air).

• Why is this cool? It means that for many different materials (from helium to iron), the relationship between pressure and melting temperature follows this same simple, curved path.
• The "Double Confirmation": The author notes that another scientist (Trachenko) recently found the exact same parabolic shape, but they used a completely different theory based on how sound waves move through liquids. It's like two people climbing a mountain from opposite sides and meeting at the exact same peak. This suggests that the "parabolic melting line" is a fundamental truth of nature, not just a lucky guess.

What the Map Tells Us
The paper claims that if you know a few basic facts about a material—how squishy it is (bulk modulus), how much it expands when hot, and how much heat it holds—you can predict its entire melting curve without needing to run expensive experiments for every single point.

In Summary
This paper says: "Stop assuming the energy to melt things is constant. It changes based on how the atoms jiggle and expand. If you account for that change, the melting line for almost any material is a simple, predictable curve (a parabola), and we can calculate it using basic physics properties."

Technical Summary

Technical Summary: Theory of Melting Lines with a Variable Enthalpy of Fusion

Problem Statement
The thermodynamic description of solid-liquid phase boundaries (melting lines, MLs) has historically been challenging compared to sublimation and boiling curves. Conventional derivations using the Clausius-Clapeyron (CC) equation typically rely on the approximation of a constant latent heat of fusion ($\delta H$) to yield analytical solutions. While this assumption facilitates the derivation of universal equations for gas-phase transitions, it fails to capture the complex thermodynamics of the solid-liquid transition, where both phases are dense, strongly interacting condensed states. Existing models often focus exclusively on the structural instabilities of the solid (e.g., Lindemann criterion, Simon-Glatzel equation) or require a definitive consensus on liquid state theory, which remains elusive. Consequently, there is a lack of a universal analytical equation for MLs derived from a two-phase thermodynamic framework that accounts for the variability of latent heat.

Methodology
This work develops a two-phase analytical model by modifying the CC equation to incorporate a variable enthalpy of fusion ($\delta H$) along the melting line. The methodology is rooted in recent theoretical progress regarding the anharmonicity of crystalline solids at elevated temperatures. Specifically, the approach utilizes the finding that the isobaric heat capacity ($C_P$) of a solid comprises a temperature-dependent vibrational component and a volume-dependent (dilatational) anharmonic component.

The derivation proceeds as follows:

1. Variable Latent Heat: The enthalpy of fusion is expressed as a function of the specific volumes of the coexisting liquid ($v_L$) and solid ($v_S$) phases, derived from the volume-dependent component of the solid's heat capacity ($C_{TE} \propto \varepsilon \beta$). This yields the relation $\delta H = \varepsilon \ln(v_L/v_S)$, where $\varepsilon$ is a solid-state parameter related to the ratio of heat capacity to thermal expansion.
2. Modified CC Equation: This variable $\delta H$ is substituted into the CC equation ($dP/dT_m = \delta H / \delta v$), creating a differential equation where the slope depends on the specific volumes of both phases.
3. Second-Order Differential Equation: By differentiating the modified CC equation with respect to temperature and applying thermodynamic relations for thermal expansion ($\beta$) and isothermal bulk modulus ($B$), a second-order linear non-homogeneous ordinary differential equation governing the melting line is derived.
4. Boundary Conditions and Approximations: The solution is constrained by physical boundary conditions: the zero-pressure melting point ($T_0, 0$) and the triple point ($T_c, P_c$). The analysis assumes that the parameter $\varepsilon$ is weakly dependent on pressure and that the term $d$ (related to the difference in thermal expansion and bulk moduli) is significantly larger than the term $b$ (related to the logarithmic derivative of the volume difference).

Key Results
The solution to the derived differential equation yields an approximate parabolic function for the melting pressure as a function of temperature, $P(T_m)$.

• Parabolic Scaling: Under the condition that $d \gg b$, the second-order differential equation simplifies to $d^2P/dT_m^2 \approx d > 0$. The general solution is a quadratic function: $P(T_m) \approx \frac{1}{2}d T_m^2 + c_3 T_m + c_4$.
• Parameter Definition: The coefficients of this parabolic function are defined exclusively by fundamental thermophysical properties, including the bulk moduli ($B_L, B_S$), thermal expansion coefficients ($\beta_L, \beta_S$), specific volumes ($v_L, v_S$), and the isobaric heat capacity of the solid ($C_{P,S}$). No empirical fitting parameters are introduced.
• Slope Prediction: The initial slope of the melting line is derived as $dP/dT_m \approx \beta_L B_L$. This result indicates that the slope is determined by the product of the liquid's thermal expansion coefficient and its isothermal bulk modulus.
• Convexity and Monotonicity: The derived function maintains a strictly monotonically increasing and convex profile for all $P \geq 0$ and $T \geq T_0$, consistent with physical expectations for melting curves.

Significance and Claims
The paper claims that this derivation provides a universal analytical framework for melting lines that overcomes the limitations of the constant-latent-heat approximation. The primary significance lies in the convergence of two distinct theoretical foundations:

1. Theoretical Convergence: The parabolic scaling law and the specific prediction for the slope ($dP/dT_m \approx \beta_L B_L$) derived here from solid-state anharmonicity and heat capacity analysis align strikingly with a recent universal model derived by Trachenko from the Phonon Theory of Liquids (PTL).
2. Two-Phase Nature: Unlike earlier one-phase models that focused solely on solid instability, this approach is a genuine two-phase theory that explicitly incorporates the thermodynamics of the liquid phase through the specific volume and bulk modulus differences.
3. Experimental Consistency: The authors note that the predicted parabolic nature and slope magnitudes are consistent with experimental data for various materials, including Ar, He, H$_2$, H$_2$O, In, and Fe, as previously discussed in the context of Trachenko's work.

The study concludes that the universal parabolic nature of melting curves is a robust feature supported by independent theoretical routes, offering a unified perspective on the thermodynamics of the solid-liquid transition without relying on arbitrary empirical fitting.