Abstract
We develop an algorithm for calculating the normal modes of vibration of mechanical systems with constraints, particularly of molecules with rigid bonds and models of rigid molecules, and use it to obtain the harmonic free energy of a crystal. The anharmonic correction is then calculated by the conventional thermodynamic integration over temperature in the NVT ensemble. Attention is paid to finite-size errors, tail corrections, thermostat choice, ergodicity, and other sources of inaccuracies. The calculated free energy of ice XIV modeled by the TIP4P/2005 potential agrees with the previously reported value and is by one order more accurate.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.