The Tin Pest Problem as a Test of Density Functionals Using High-Throughput Calculations

Michael J. Mehl,Mateo Ronquillo,David Hicks,Marco Esters,Corey Oses,Rico Friedrich,Andriy Smolyanyuk,Eric Gossett,Daniel Finkenstadt,Stefano Curtarolo
DOI: https://doi.org/10.1103/PhysRevMaterials.5.083608
2020-10-14
Abstract:At ambient pressure tin transforms from its ground-state semi-metal $\alpha$-Sn (diamond structure) phase to the compact metallic $\beta$-Sn phase at 13$^\circ$C (286K). There may be a further transition to the simple hexagonal $\gamma$-Sn above 450K. These relatively low transition temperatures are due to the small energy differences between the structures, $\approx 20$\,meV/atom between $\alpha$- and $\beta$-Sn. This makes tin an exceptionally sensitive test of the accuracy of density functionals and computational methods. Here we use the high-throughput Automatic-FLOW (AFLOW) method to study the energetics of tin in multiple structures using a variety of density functionals. We look at the successes and deficiencies of each functional. As no functional is completely satisfactory, we look Hubbard U corrections and show that the Coulomb interaction can be chosen to predict the correct phase transition temperature. We also discuss the necessity of testing high-throughput calculations for convergence for systems with small energy differences.
Materials Science
What problem does this paper attempt to address?
### The problems the paper attempts to solve This paper aims to solve the accuracy problem of density functional theory (DFT) in different phase - transition processes of tin. Specifically, the research focuses on: 1. **Phase transition from α - Sn to β - Sn**: Tin transforms from the semimetallic α - Sn (diamond structure) phase to the metallic β - Sn phase at normal pressure, and the transition temperature is 13°C (286K). This phase transition leads to the so - called "tin pest" problem, that is, significant volume expansion occurs when β - Sn is stored at low temperatures. 2. **Accuracy test of DFT in systems with extremely small energy differences**: Since the energy difference between α - Sn and β - Sn is very small (about 20 meV/atom), tin becomes an ideal system for testing the accuracy of DFT methods. Different DFT functionals may give different prediction results, so high - throughput calculations are required to evaluate the performance of various functionals. 3. **Application of Hubbard U correction**: In order to predict the phase - transition temperature more accurately, the author introduced Hubbard U correction and showed how to predict the correct phase - transition temperature by adjusting the Coulomb interaction. 4. **Convergence test of high - throughput calculations**: For systems with very small energy differences, the default settings of high - throughput calculations may not be sufficient to find the correct phase order. Therefore, the author emphasized the importance of testing convergence. ### Main content of the paper - **Introduction**: Introduced the phase - transition phenomenon of tin and its importance in practical applications, especially the "tin pest" problem in lead - free solders. - **Calculation details**: Described the first - principles calculations using the AFLOW platform and VASP, including the calculations of lattice energy, phonon spectra, and free energy. - **Crystal structures**: Described in detail the various crystal structures of tin involved in the research, including α - Sn, β - Sn, and γ - Sn. - **Selection and evaluation of DFT functionals**: Compared the prediction abilities of different functionals such as LDA, GGA, and meta - GGA for tin phase transitions, and especially emphasized the performance of the SCAN functional. - **Thermodynamic properties**: Analyzed the thermodynamic properties of α - Sn, β - Sn, and γ - Sn, including the prediction of phase - transition temperatures. - **Hubbard U correction**: Improved the prediction results of LDA and PBE functionals by introducing Hubbard U correction. - **Convergence discussion**: Discussed the convergence problem of high - throughput calculations in systems with very small energy differences. ### Conclusion This paper verified the superiority of the SCAN functional in predicting tin phase transitions through high - throughput calculations and comparison of multiple DFT functionals, but also pointed out its overestimation problem of energy differences. At the same time, the author emphasized the importance of ensuring calculation convergence when dealing with systems with extremely small energy differences.