Kinetic analysis of phase transformations during continuous heating: Crystallization of glass-forming liquids

Basic Information

• Paper ID: 2504.14029
• Title: Kinetic analysis of phase transformations during continuous heating: Crystallization of glass-forming liquids
• Author: O.S. Houghton (MIT Materials Research Laboratory)
• Classification: cond-mat.mtrl-sci (Materials Science)
• Publishing Institution: Massachusetts Institute of Technology
• Paper Link: https://arxiv.org/abs/2504.14029

Abstract

This paper presents a breakthrough theoretical framework for analyzing phase transformation kinetics during continuous heating. While phase transformation kinetics under isothermal conditions have been well-described by the Johnson-Mehl-Avrami-Kolmogorov (JMAK) rate equation, phase transformation analysis under continuous heating conditions has previously been achievable only through numerical methods. This study derives for the first time the analytical solution to the JMAK rate equation for phase transformations during continuous heating and applies it to analyze crystallization behavior in glass-forming liquids. The model's validity is verified using experimental data from Fe₈₀B₂₀ metallic glass.

Research Background and Motivation

Importance of the Problem

1. Technical Requirements: Glass materials have important applications in metallic, organic, inorganic, and hybrid material fields. Understanding crystallization kinetics is crucial for developing high-performance glass materials.
2. Process Requirements: Continuous heating experiments offer advantages over isothermal experiments, including faster processing, wider temperature ranges, and ability to study rapid transformations, making them more representative of actual production processes.
3. Theoretical Deficiency: Despite widespread application of continuous heating experiments, a corresponding analytical theoretical framework is lacking, requiring reliance on numerical simulations.

Limitations of Existing Methods

1. Numerical Method Dependence: Current analysis requires complex numerical simulations to fit the shape and position of crystallization exothermic peaks.
2. Low Computational Efficiency: Numerical methods are computationally expensive, hindering parametric modeling and material design.
3. Limited Physical Insight: The absence of analytical expressions restricts deeper understanding of physical behavior.

Core Contributions

1. Theoretical Breakthrough: First derivation of the analytical solution to the JMAK rate equation for phase transformations during continuous heating.
2. Model Verification: Validation of model accuracy through experimental data and numerical simulations of Fe₈₀B₂₀ metallic glass.
3. Theoretical Foundation for Kissinger Method: Provides rigorous theoretical basis for the widely-used Kissinger kinetic analysis method.
4. Non-Arrhenius Behavior Analysis: Demonstrates that the model can better determine the non-Arrhenius temperature dependence of crystal growth rates in glass-forming liquids.
5. Broad Applicability: Provides analytical tools for design and development of glass-forming systems and can be extended to other phase transformation processes.

Detailed Methodology

Theoretical Foundation

This research is based on four fundamental assumptions of JMAK theory:

1. Spatially random nucleation
2. Growth rate of the new phase depends only on temperature, not on time
3. Transformation rate dX/dT is independent of thermal history, depending only on transformation fraction X and temperature T
4. Growth of the new phase is not limited by interfacial heat transfer

Core Equation Derivation

Isothermal JMAK equation:

X_iso(t) = 1 - exp(-(Kt)^n)

Integral equation under continuous heating conditions:

∫[0 to T] dT'/[K(T')]^(1/n) = n/Φ ∫[0 to T] K(T')dT'

where Φ is the constant heating rate.

Kinetic coefficient expression (fixed nucleation density, 3D spherical growth, n=3):

K(T) = (4π/3)N^(4/3)U(T)

Key Innovation: Analytical Solution

Analytical solution under VFT viscosity model:

X(T,Φ) = 1 - exp[-(αΦ)^(1/n)]

where α is a complex temperature-dependent parameter containing exponential integral functions:

α = complex expression involving D*, T₀ and other VFT parameters

Simplified solution under Arrhenius approximation:

α = simplified exponential integral expression involving activation energy E_η

Technical Innovations

1. Mathematical Breakthrough: Successfully solved complex integrals containing temperature-dependent exponential functions.
2. Physical Modeling: Organically combined the VFT viscosity model with crystal growth kinetics.
3. Parameter Decoupling: Achieved independent analysis of nucleation density N and growth kinetics.
4. Fragility Effect: Introduced viscosity decoupling factor ε to describe behavior in highly fragile systems.

Experimental Setup

Verification Material

• Primary Material: Fe₈₀B₂₀ metallic glass (atomic percentage)
• Selection Rationale: The crystallization behavior of this material strictly satisfies model assumptions, and detailed experimental and numerical studies by Greer provide a reliable comparison benchmark.

Experimental Conditions

• Heating Rates: 10-80 K/min
• Testing Method: Differential Scanning Calorimetry (DSC)
• Sample Form: Melt-spun ribbon

Model Parameters

• Arrhenius Model: E_η = 241 kJ/mol
• VFT Model: D* = 684, T₀ = 40 K
• Fragility Decoupling Model: D* = 94, T₀ = 190 K

Experimental Results

Main Results

1. Model Verification Accuracy:
   • Peak temperature deviation between analytical model and experimental data is only 4-6 K
   • Successfully reproduced peak shape and position at different heating rates
   • VFT model provides better fitting than Arrhenius model
2. Theoretical Verification of Kissinger Equation:

   ln(Φ/T_x²) = (n/3)ln(4π/3) + ln(fa₀³/kη_∞) + ln α - ln(η)
   

   
   Demonstrates the relationship between Kissinger plot slope and effective activation energy.
3. Fragility Effect Analysis:
   • Crystal growth in highly fragile liquids (m > 80) decouples from viscosity
   • Decoupling factor ε = 1.1 - 0.005m significantly affects crystallization behavior prediction

Parameter Sensitivity Analysis

1. Effect of Activation Energy E_η:
   • Higher E_η leads to earlier crystallization onset and sharper exothermic peaks
   • Requires reduced nucleation density N to maintain the same T_x
2. Effect of Nucleation Density N:
   • Increased N leads to lower T_x but minimal peak shape change
   • Effect on peak shape is smaller than fragility parameter m
3. Effect of Avrami Index n:
   • Decreased n leads to slight reduction in T_x
   • Peak becomes broader with reduced intensity

Experimental Findings

1. Non-Arrhenius Behavior Detection: By fitting peak shape, position, and height, non-Arrhenius behavior can be detected within the width range of individual crystallization exothermic peaks.
2. Extended Temperature Range: Peak fitting method provides crystal growth rate data over a wider temperature range compared to Kissinger method.
3. Fragility-Crystallization Relationship: Establishes quantitative relationship between liquid fragility and crystallization resistance.

Related Work

Historical Development

1. JMAK Theory: Theoretical foundation for isothermal phase transformation kinetics established over 80 years ago.
2. Kissinger Method: Continuous heating kinetic analysis method proposed in 1957.
3. Numerical Simulation: Detailed numerical studies of Fe₈₀B₂₀ by Greer and colleagues.

Advantages of This Work

1. First Analytical Solution: Fills the gap in continuous heating phase transformation theory.
2. Computational Efficiency: Significantly improves computational efficiency compared to numerical methods.
3. Physical Insight: Provides direct relationships between parameters and physical properties.

Conclusions and Discussion

Main Conclusions

1. Theoretical Achievement: Successfully derives the analytical solution to the JMAK equation for phase transformations during continuous heating.
2. Practical Value: Provides powerful theoretical tools for glass material design and process optimization.
3. Method Verification: Experimental verification using Fe₈₀B₂₀ demonstrates method accuracy and practicality.
4. Broad Applicability: The method can be extended to other phase transformation processes satisfying JMAK assumptions.

Limitations

1. Assumption Requirements: Must satisfy the four fundamental assumptions of JMAK theory.
2. Temperature Range: Thermodynamic corrections may be needed at high heating rates.
3. Nucleation Mode: Primarily applicable to site-saturation conditions; continuous nucleation requires numerical treatment.
4. Boundary Conditions: Assumptions about integral boundary conditions may affect accuracy.

Future Directions

1. Extended Applications: Generalization to other material systems and phase transformation types.
2. Process Optimization: Application to non-isothermal processes such as additive manufacturing and rapid sintering.
3. Multi-step Phase Transformations: Development of theoretical frameworks for complex multi-step transformations.
4. Real-time Prediction: Integration with machine learning for real-time material property prediction.

In-Depth Evaluation

Strengths

1. Theoretical Breakthrough: First solution to the analytical description of phase transformation kinetics during continuous heating, with significant theoretical value.
2. Mathematical Rigor: Derivation is rigorous, employing advanced mathematical tools such as exponential integral functions.
3. Sufficient Experimental Verification: Verified through experimental data at multiple heating rates with convincing results.
4. High Practical Value: Provides directly applicable tools for material design and process optimization.
5. Deep Physical Insight: Reveals quantitative effects of parameters such as fragility and nucleation density on crystallization behavior.

Limitations

1. Restricted Applicability: Strictly depends on JMAK assumptions, potentially inapplicable to certain complex systems.
2. Missing High-Temperature Corrections: Lacks detailed treatment of thermodynamic corrections at high heating rates.
3. Difficult Parameter Determination: Accurate determination of VFT model parameters requires extensive experimental data.
4. Multi-phase Competition: Insufficient consideration of complex multi-phase competitive crystallization.

Impact

1. Academic Value: Expected to become an important theoretical foundation in phase transformation kinetics.
2. Industrial Application: Direct application value for glass industry, metal material processing, and related fields.
3. Methodological Contribution: Provides a model for analytical modeling of other complex physical processes.
4. Educational Significance: Provides a complete theoretical framework for phase transformation kinetics teaching.

Applicable Scenarios

1. Material Design: Composition optimization and performance prediction of novel glass materials.
2. Process Control: Parameter optimization and quality control of heat treatment processes.
3. Rapid Screening: Fast performance evaluation of numerous candidate materials.
4. Theoretical Research: Deeper understanding of phase transformation kinetics mechanisms.

References

The paper cites 36 important references covering the historical development of JMAK theory, glass science fundamentals, and numerical simulation methods. Particularly important works include:

• Avrami (1939, 1940): Foundational work on JMAK theory
• Greer (1982): Classical study of Fe₈₀B₂₀ crystallization kinetics
• Kissinger (1957): Continuous heating kinetic analysis method
• Henderson (1979): Early theoretical work on non-isothermal crystallization kinetics

Overall Assessment: This is an excellent paper with significant theoretical value and practical importance. It solves for the first time the problem of analytical description of phase transformation kinetics during continuous heating, providing powerful theoretical tools for the materials science field. The theoretical derivation is rigorous, experimental verification is comprehensive, and the work is expected to have important impacts on phase transformation kinetics research and glass material development.