472216 First-Principles Phase Diagrams: Iron at Earth’s Inner Core Conditions with Full Inclusion of Anharmonic and Finite-Size Effects
In this work , we utilize a novel machinery developed by our group (called “harmonically-mapped averaging (HMA)” method [6,7]) to tackle this long standing problem using DFT method as implemented in VASP. The speed up gained by using HMA method to measure thermodynamic properties is at least an order of magnitude faster than “direct” methods. This enabled measuring of the full (harmonic and anharmonic) Gibbs free energy of Iron at several states. In addition, the finite-size effects in free energy were estimated based on a harmonic-based technique described here . Applying these techniques to the three known iron candidate structures (hcp, fcc, and bcc), we were able, to the first time, to build a completely ab initio-based phase diagram of iron at Earth’s IC conditions in the thermodynamic limit.
 F. Birch, J. Geophys. Res. 57, 227 (1952).
 P. Soderlind, J. Moriarty, and J. Wills, Phys. Rev. B 53, 14063 (1996).
 L. Stixrude, Phys. Rev. Lett. 108, 055505 (2012).
 W. Luo, B. Johansson, O. Eriksson, S. Arapan, P. Souvatzis, M. I. Katsnelson, and R. Ahuja, Proc. Natl. Acad. Sci. 107,9962 (2010).
 S. G. Moustafa, A. J. Schultz, E. Zurek, and D. A. Kofke, In preparation, (2016).
 S. G. Moustafa, A. J. Schultz, D. A. Kofke, Phys. Rev. E 2015, 92, 043303.
 A. J. Schultz, S. G. Moustafa, L. Weisong, S. J. Weinstein, and D. A. Kofke, J. Chem. Theory Comput., 2016, 12 (4), pp 1491-1498.
 T. B. Tan, A. J. Schultz, and D. A. Kofke, J. Chem. Phys. 133, 134104 (2010).