Abstract:The calculation of free energies from first principles in materials is a formidable task which enables the prediction of phase stability with high accuracy; these calculations are complicated in magnetic materials by the interplay of electronic, magnetic, and vibrational degrees of freedom. In this work, we show the feasibility and accuracy of the calculation of phase stability in magnetic systems with ab initio methods and thermodynamic integration by sampling the magnetic and vibrational phase space with coupled atomistic spin dynamics-ab initio molecular dynamics (ASD-AIMD) simulations [Stockem et al., PRL 121, 125902 (2018)], where energies and interatomic forces are calculated with density functional theory (DFT). We employ the method to calculate the phase stability of Fe at ambient pressure from 800 K up to 1800 K. The Gibbs free energy difference between fcc and bcc Fe at zero pressure as a function of temperature is calculated carrying out thermodynamic integration over temperature on the energies at the DFT level from ASD-AIMD, using a reference free energy difference calculated in the paramagnetic state at temperatures much higher than the magnetic transition temperatures with thermodynamic integration over stress-strain variables with disordered local moment (DLM)-AIMD simulations. We show the importance of the magnetic ordering temperature of bcc Fe on the $\alpha$ to $\gamma$ structural transition temperature, whereas the $\gamma$ to $\delta$ transition is well reproduced independently of the exchange interactions. The Gibbs free energy difference between the two structures is within 5 meV/atom from the CALPHAD estimate, and both transition temperatures are reproduced within 150 K. The present work paves the way to free energy calculations in magnetic materials from first principles with accuracy in the order of 1 meV/atom.

What problem does this paper attempt to address?

The problem this paper attempts to solve is: from first principles, to calculate the phase stability of iron (Fe) at different temperatures, particularly the phase stability under the interaction of magnetic and vibrational degrees of freedom. Specifically, the authors aim to calculate the phase stability of iron at ambient pressure from 800 K to 1800 K, especially the Gibbs free energy difference between the face-centered cubic (fcc) and body-centered cubic (bcc) phases, by combining atomic spin dynamics (ASD) and ab initio molecular dynamics (AIMD) simulations with thermodynamic integration methods. ### Main Research Content: 1. **Calculation Methods**: - Use density functional theory (DFT) to calculate energy and interatomic forces. - Combine ASD-AIMD simulations to sample magnetic and vibrational phase space. - Calculate the free energy differences at different temperatures using thermodynamic integration methods. 2. **Research Background**: - At low temperatures, iron is in the body-centered cubic (bcc) ferromagnetic phase (α phase). At 1043 K, it undergoes a magnetic phase transition to a paramagnetic bcc phase, followed by a structural phase transition to a face-centered cubic (fcc) phase (γ phase) at 1185 K, and finally transitions back to a bcc phase (δ phase) at 1667 K. - Experimentally estimated Gibbs free energy difference between the fcc and bcc phases within the stability range of the γ phase is about 1 meV/atom, requiring extremely high precision in the calculation methods. 3. **Research Objectives**: - Calculate the phase stability of iron at different temperatures using ASD-AIMD and DLM-AIMD simulations combined with thermodynamic integration methods. - Pay special attention to the impact of magnetic exchange interactions on the phase transition temperature and the coupling effect of magnetic and vibrational degrees of freedom. ### Research Methods: 1. **ASD-AIMD Simulations**: - At each AIMD time step, use DFT to calculate interatomic forces and consider the direction of magnetic moments. - The ASD part simulates the evolution of magnetic moments, considering the exchange interactions between neighboring magnetic moments. 2. **DLM-AIMD Simulations**: - Simulate in the ideal paramagnetic state, ignoring the correlations between magnetic moments. - Used to calculate the reference free energy differences at high temperatures. 3. **Thermodynamic Integration Path**: - Start from the DLM state at 1800 K to calculate the reference free energy differences. - Calculate the free energy differences at different temperatures through temperature integration at constant volume. - Finally, add the contributions of thermal expansion/contraction. ### Research Results: - The calculated Gibbs free energy difference between the fcc and bcc phases differs from the CALPHAD estimated value by about 5 meV/atom. - The phase transition temperature differs from the experimental value by about 150 K. - Magnetic exchange interactions significantly affect the α to γ structural phase transition temperature, while the effect on the γ to δ transition is smaller. ### Conclusion: - This study demonstrates the feasibility and accuracy of combining ASD-AIMD and thermodynamic integration methods in calculating the phase stability of magnetic materials. - It provides a new approach for calculating the free energy of magnetic materials from first principles in the future.

Phase stability of Fe from first-principles: atomistic spin dynamics coupled with ab initio molecular dynamics simulations and thermodynamic integration

The magnetic phase transition in Mn2-xFexB alloys: First-principles calculations

The Mechanism of Structure Phase Transition from Α Fe to Ε Fe under Uniaxial Strain: First-principles Calculations

Ab Initio Study Of The Anharmonic Lattice Dynamics Of Iron At The Gamma-Delta Phase Transition

A first-principles theory of ferromagnetic phase transitions in metals

Thermodynamic Integration for Dynamically Unstable Systems Using Interatomic Force Constants without Molecular Dynamics

Longitudinal Spin Fluctuations and Magnetic Ordering Temperature in Metals: First-Principle Modeling and Phase Space Integration Measure

Ab initio studies on structural and thermodynamic properties of magnetic Fe

Electronic correlations at the alpha-gamma structural phase transition in paramagnetic iron

Effective Ab Initio Exchange Integrals And Magnetic Phase Transition In Fcc Fe And Mn Antiferromagnets

Crystal structure, thermodynamics, magnetics and disorder properties of Be-Fe-Al intermetallics

Phase stability of ternary fcc and bcc Fe-Cr-Ni alloys

Ab Initio Study of the β-Fe2O3 Phase

Structural dynamics of first-order phase transition in giant magnetocaloric La(Fe,Si) 13 : The free energy landscape

Thermodynamic Calculation Of Phase Diagram And Phase Stability With Nano-Size Particles

Analysis of the determination of magnetic phase transition temperatures by the first-principles calculations

First-principles predictions of structural and magnetic phase stability in irradiated α-Fe

Ab initio phase stabilities of Ce-based hard magnetic materials and comparison with experimental phase diagrams

Energetics and thermodynamics of α-iron from first-principles and machine-learning potentials

Subphase exploration for SmFe12-based permanent magnets by Gibbs energies obtained with first-principles cluster-expansion method

Understanding the phase stability in multi-principal-component AlCuFeMn alloy