First-principles modeling of temperature and concentration dependent solubility in the phase separating Fe$_x$Cu$_{1-x}$ alloy system

D. Reith,M. Stöhr,R. Podloucky,T. C. Kerscher,S. Müller
DOI: https://doi.org/10.1103/PhysRevB.86.020201
2012-04-16
Abstract:We present a novel cluster-expansion (CE) approach for the first-principles modeling of temperature and concentration dependent alloy properties. While the standard CE method includes temperature effects only via the configurational entropy in Monte Carlo simulations, our strategy also covers the first-principles free energies of lattice vibrations. To this end, the effective cluster interactions of the CE have been rendered genuinely temperature dependent, so that they can include the vibrational free energies of the input structures. As a model system we use the phase-separating alloy Fe-Cu with our focus on the Fe-rich side. There, the solubility is derived from Monte Carlo simulations, whose precision had to be increased by averaging multiple CEs. We show that including the vibrational free energy is absolutely vital for the correct first-principles prediction of Cu solubility in the bcc Fe matrix: The solubility tremendously increases and is now in quantitative agreement with experimental findings.
Materials Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **How to accurately predict the temperature - and concentration - dependent alloy solubility through first - principles modeling, especially the solubility of Cu in the bcc Fe matrix in the Fe - Cu alloy system**. Specifically, the traditional Cluster Expansion (CE) method only considers the temperature effect through the configurational entropy in Monte Carlo (MC) simulations, while ignoring the important contribution of the vibrational free energy. This has led to a significant deviation between the predicted solubility of Cu in the Fe matrix and the experimental results. Therefore, this study aims to develop an improved CE method that incorporates the temperature - dependent vibrational free energy into the model, thereby predicting the solubility of Cu in the Fe matrix more accurately. ### Main problems: 1. **Limitations of the traditional CE method**: The traditional CE method only considers the temperature effect through configurational entropy, ignoring other temperature - dependent terms such as the vibrational free energy. 2. **Inaccurate prediction of the solubility of Cu in the Fe matrix**: Due to the ignorance of the vibrational free energy, the Cu solubility predicted by the traditional method is much lower than the experimental value. 3. **Need for higher prediction accuracy**: Especially in the case of low - concentration Cu, a high - precision CE model is required to accurately predict the solubility. ### Solutions: - **Introduce temperature - dependent effective cluster interaction energies (ECIs)**: Make ECIs a function of temperature \( J_f(T) \), thereby incorporating the vibrational free energy \( F_{\text{vib},\text{DFT}}(\sigma, T) \) into the model. - **Combine with MC simulations**: Use MC simulations to handle the phase stability and related properties of a large number of interacting atoms and consider the temperature effect. - **Averaging strategy**: In order to improve the prediction accuracy, use multiple CE models for averaging to reduce the fluctuations in a single simulation. Through these improvements, the researchers have demonstrated the importance of including the vibrational free energy and proved its significant impact on the Cu solubility, making the prediction results quantitatively consistent with the experimental data. ### Formula summary: 1. **Cluster expansion formula**: \[ E_{\text{CE}}(\sigma)=\sum_{f \in F} D_f J_f \Pi_f(\sigma) \] where \( \Pi_f(\sigma) \) is the geometric correlation, \( D_f \) is the symmetry degeneracy, and \( J_f \) is the effective cluster interaction energy. 2. **Temperature - dependent energy expression**: \[ E_{\text{DFT}}(\sigma, T)=E_0^{\text{DFT}}(\sigma)+F_{\text{vib},\text{DFT}}(\sigma, T) \] where \( E_0^{\text{DFT}}(\sigma) \) is the total ground - state energy obtained from standard DFT calculations, and \( F_{\text{vib},\text{DFT}}(\sigma, T) \) is the vibrational free energy. 3. **Improved CE formula**: \[ E_{\text{CE}}(\sigma, T)=\sum_{f \in F(T)} D_f J_f(T) \Pi_f(\sigma) \] where \( F(T) \) and \( J_f(T) \) are functions of temperature. 4. **Averaged CE formula**: \[ \bar{E}_{\text{CE}}(\sigma, T)=\frac{1}{n} \sum_{i = 1}^{n} E_{\text{CE},i}(\sigma, T)=\sum_{f \in \bar{F}(T)} D_f \bar{J}_f(T) \Pi_f(\sigma) \] Through these methods, the researchers have successfully improved the prediction accuracy of the solubility of Cu in the Fe matrix and verified the important influence of the vibrational free energy on the phase stability and solubility.