Calculating the radial distribution functions of supercritical methanol by means of Car-Parrinello and classical molecular dynamics

Radial distribution functions and the average number of hydrogen bonds per methanol molecule under standard, subcritical, and supercritical conditions are obtained via classical molecular dynamics and Car-Parrinello nonempirical molecular dynamics. It is shown that independent methods of modeling yield close results. It is noted that the calculated radial distribution functions agree well with the experimental data only at T = 298 K and P = 0.1 MPa, while at high temperatures and pressures, considerable divergence from the experimental functions known from the literature is observed. It is concluded that both modeling methods reproduce the degree of hydrogen bonding in methanol and its variations depending on the state parameters and correspond closely to the experimental results.