relax_box calculation style
Lucas M. Hale, lucas.hale@nist.gov, Materials Science and Engineering Division, NIST.
Introduction
The free_energy_liquid calculation style uses thermodynamic integration to evaluate the absolute Helmholtz and Gibbs free energies of a liquid phase by integrating from a reference Uhlenbeck-Ford model potential.
Version notes
2022-09-20: Calculation first added to iprPy
Additional dependencies
Disclaimers
The simulations used by this method are all fixed volume and therefore will not relax the system to a given pressure. Instead, it is expected that the given system is already relaxed to the target pressure of interest. Specifically, the final dump files generated by the relax_liquid calculation serve as ideal input configurations for this calculation.
Method and Theory
The calculation performs two different simulations: an nvt simulation to estimate Einstein solid spring constants for each atom type, and thermodynamic integrations from the selected interatomic potential to the Einstein solid model and back. The method follows what is described in Freitas, Asta, de Koning, Computational Materials Science 112 (2016) 333–341.
The Einstein solid spring constants, \(k_i\), are evaluated using an nvt simulation run and measuring the mean squared displacements, \(\left<\left( \Delta r_i \right)^2\right>\), averaged for each atom type \(i\) and over time
For the integration, TBD…