Temperature dependence of the microscopic structure and density anomaly of the SPC/E and TIP4P-Ew water models. Molecular dynamics simulation results

Abstract

We have investigated temperature trends of the microscopic structure of the SPC/E and TIP4P-Ew water models in terms of the pair distribution functions, coordination numbers, the average number of hydrogen bonds, the distribution of bonding states of a single molecule as well as the angular distribution of molecules by using the constant pressure molecular dynamics simulations. The evolution of the structure is put in correspondence with the dependence of water density on high temperatures down to the region of temperatures where the system becomes supercooled. It is shown that the fraction of molecules with three and four bonds determine the maximum density for both models. Moreover, the temperature dependence of the dielectric constant is obtained and analyzed.

Ми дослiдили температурнi залежностi мiкроскопiчної структури моделей води SPC/E i TIP4P-Ew в термiнах парних функцiй розподiлу, координацiйних чисел, середнього числа водневих зв’язкiв, розподiлу зв’язаних станiв однiєї молекули, використовуючи метод молекулярної динамiки при постiйногму тиску. Еволюцiя структури поставлена у вiдповiднiсть до залежностi густини води вiд температури в областi, що знаходиться мiж високими температурами i температурами, де система стає переохолодженою. Показано, що частка молекул з трьома i чотирма зв’язками визначає максимальну густину для обох моделей. Бiльш того, отримано i проаналiзовано температурну залежнiсть дiелектричної сталої.