Энергия Гельмгольца воды

[alenov]

Добрый день! Впервые пишу в этой ветке, пожалуйста, не судите строго, если что-то не в формате, то сразу же исправлю. Мучаюсь с расчетом энергии Гельмгольца для воды. Использую алгоритм Hill, 1990. В общем старый и хорошо известный алгоритм. С его помощью я могу расчитать значения безразмерной функции Гельмгольца и ее производные. У Хилла (стр. 1235) написано, что справедливо следующее соотношение: $\psi_f=A/(T R), где $\psi_f$ - функция Гельмгольца, А - энергия Гельмгольца, Т - текущая температура и R - газовая постоянная. Из энергии Гельмгольца мне нужно расчитать энергию Гиббса. Для этого я использую уравнение 52 на странице 26 (Oelkers et al., 2009). Из этого уравнения при стандартных условиях Т $ = 298.15$°К и Р $ = 1 $ бар выражение $A+\rho (\frac{dA}{d\rho$$})_T = - 398 cal mol^{-1}$. Исследовав эту функцию, я пришел к выводу, что большая часть этого значения должна приходиться именно на А. Проблема в том, что при данных условиях $\psi_f = -0.0339$ и я уверен, что результат верный (проверял с помощью стороннего алгоритма). Если просто умножить на 298.15°К и 1.9872041 cal/°K/mol, то получается -20 cal/mol - примерно в 20 раз меньше требуемого значения. (Значение производной положительное и на 2 порядка меньше значения энергии Гельмгольца.) Может быть кто-нибудь может указать на мою ошибку?

Ссылки на цитированные статьи:
https://drive.google.com/file/d/0B08REa ... sp=sharing
https://drive.google.com/file/d/0B08REa ... sp=sharing

[alenov]

Я разобрался. Уравнение в Oelkers et al., 2009 неправильное. Нужно вводить поправки на значения энтальпии, энтропии и энергии Гиббса в тройной точке. Очень помогли следующие статьи:
https://drive.google.com/file/d/0B08REa ... sp=sharing
http://www.ajsonline.org/content/291/6/541.citation
Последнюю мне удалось посмотреть только в читальном зале - отсканировать не дали, к сожалению.