Термостат Нозе — Гувера

Термостат Нозе — Гувера — это алгоритм молекулярной динамики, предназначенный для поддержания заданной температуры системы путём термостатирования, то есть искусственного контроля кинетической энергии частиц. Он относится к классы детерминированных термостатов, в отличие от стохастических методов (например, термостата Андерсена), и основан на введении дополнительной степени свободы, которая моделирует взаимодействие системы с тепловым резервуаром. Метод был разработан в 1984 году голландским физиком Шуити Нозе и усовершенствован американским физиком Уильямом Гувером. Термостат Нозе — Гувера широко применяется в компьютерном моделировании физических, химических и биологических систем, обеспечивая корректное воспроизведение канонического ансамбля (NVT) в равновесной статистической механике.

История

Предпосылки создания

До появления термостата Нозе — Гувера в молекулярной динамике использовались простые методы поддержания температуры, такие как масштабирование скоростей (алгоритм Берендсена), который не обеспечивал строгого соответствия каноническому ансамблю. В 1980 году Андерсен предложил стохастический термостат, где скорости частиц случайным образом переприсваивались из распределения Максвелла — Больцмана, что нарушало детерминированность траекторий. Это ограничивало применение методов для изучения динамических свойств, таких как транспортные коэффициенты.

Разработка Нозе

В 1984 году Шуити Нозе (Токийский университет) опубликовал работу, в которой предложил расширенный лагранжиан, включающий дополнительную координату \( s \), связанную с тепловым резервуаром. Эта координата вводила масштабирование времени, что позволяло системе естественным образом флуктуировать вокруг заданной температуры. Нозе показал, что при усреднении по траекториям такой системы можно получить распределение канонического ансамбля. Однако его исходный формализм приводил к нефизическим временным масштабам.

Усовершенствование Гувера

В 1985 году Уильям Гувер (Калифорнийский университет в Дэвисе) модифицировал подход Нозе, заменив масштабирование времени на введение дополнительной переменной — «трения» \( \zeta \), которая управляет скоростями частиц. Это позволило получить уравнения движения, не требующие переопределения времени. Так появился термостат Нозе — Гувера в современной форме. Впоследствии метод был расширен для работы с несколькими термостатами (термостат Нозе — Гувера — Мартыны) и для моделирования неравновесных процессов.

Принцип работы

Основная идея

Термостат Нозе — Гувера вводит в систему дополнительную степень свободы, которая взаимодействует с частицами через фиктивную силу трения. Эта сила регулируется так, чтобы средняя кинетическая энергия системы соответствовала заданной температуре \( T \). Уравнения движения имеют вид:

\[ \dot{\mathbf{r}}_i = \frac{\mathbf{p}_i}{m_i}, \quad \dot{\mathbf{p}}_i = \mathbf{F}_i - \zeta \mathbf{p}_i, \quad \dot{\zeta} = \frac{1}{Q} \left( \sum_{i=1}^N \frac{\mathbf{p}_i^2}{m_i} - N_f k_B T \right), \]

где:

• \( \mathbf{r}_i \) и \( \mathbf{p}_i \) — координаты и импульсы частиц,
• \( m_i \) — массы частиц,
• \( \mathbf{F}_i \) — силы, действующие на частицы,
• \( \zeta \) — термостатная переменная (коэффициент трения),
• \( Q \) — эффективная масса термостата (параметр, определяющий инерцию),
• \( N_f \) — число степеней свободы системы,
• \( k_B \) — постоянная Больцмана.

Роль параметра \( Q \)

Параметр \( Q \) (или его обратная величина \( \tau \)) определяет, насколько быстро термостат реагирует на отклонения температуры. При малом \( Q \) термостат быстро корректирует температуру, но может вызвать нефизические осцилляции. При большом \( Q \) система медленно приходит к равновесию. Оптимальное значение \( Q \) выбирается эмпирически, часто на основе времени релаксации \( \tau \approx 0.1-1 \) пс.

Отличие от других термостатов

• Термостат Берендсена: масштабирует скорости, но не даёт правильного канонического распределения (не сохраняет флуктуации энергии).
• Термостат Андерсена: стохастический, нарушает динамику.
• Термостат Нозе — Гувера: детерминированный, сохраняет траектории в фазовом пространстве, но может быть неэргодичным для малых систем (например, для гармонического осциллятора), что потребовало создания модификаций — термостата Нозе — Гувера — Мартыны.

Применение

Молекулярная динамика

Термостат Нозе — Гувера является стандартным методом в программных пакетах молекулярного моделирования, таких как LAMMPS, GROMACS, NAMD, AMBER. Он используется для:

• Моделирования жидкостей и твёрдых тел при постоянной температуре (например, воды, белков, полимеров).
• Изучения фазовых переходов и критических явлений.
• Расчёта термодинамических свойств (теплоёмкость, энтропия) через флуктуации.

Статистическая механика

Метод позволяет генерировать траектории, соответствующие каноническому ансамблю (NVT), что необходимо для вычисления средних по ансамблю величин, таких как свободная энергия, с помощью методов термодинамического интегрирования или зонтичной выборки.

Неравновесные процессы

В модификациях (например, с несколькими термостатами) метод применяется для моделирования теплопроводности, вязкости и других транспортных свойств, где требуется поддержание градиента температуры.

Ограничения и критика

Эргодичность

Для простых систем с малым числом степеней свободы (например, одномерный гармонический осциллятор) термостат Нозе — Гувера может быть неэргодичным, то есть не охватывать все возможные состояния фазового пространства. Это приводит к неверным распределениям. Для решения этой проблемы были разработаны:

• Термостат Нозе — Гувера — Мартыны (1988): использует цепочку из нескольких термостатов, соединённых последовательно.
• Термостат с кинетической энергией (например, метод Брауна — Кларка).

Выбор параметров

Неправильный выбор массы термостата \( Q \) может привести к осцилляциям температуры или медленной сходимости. Для сложных систем (например, биомолекул) часто требуется калибровка.

Вычислительные затраты

Введение дополнительной переменной увеличивает число уравнений, что незначительно повышает вычислительную сложность, но для больших систем (миллионы частиц) это может быть заметно.

Интересные факты

• Шуити Нозе первоначально работал над теорией фазовых переходов, и его термостат был побочным результатом исследований.
• Уильям Гувер известен также работами в области гидродинамики и теории хаоса.
• Термостат Нозе — Гувера является одним из немногих детерминированных методов, строго воспроизводящих каноническое распределение.
• В 2010-х годах были предложены квантовые обобщения метода для моделирования систем при низких температурах.

Источники

• Nosé, S. (1984). «A unified formulation of the constant temperature molecular dynamics methods». Journal of Chemical Physics, 81(1), 511–519.
• Hoover, W. G. (1985). «Canonical dynamics: Equilibrium phase-space distributions». Physical Review A, 31(3), 1695–1697.
• Martyna, G. J., Klein, M. L., & Tuckerman, M. (1992). «Nosé–Hoover chains: The canonical ensemble via continuous dynamics». Journal of Chemical Physics, 97(4), 2635–2643.
• Frenkel, D., & Smit, B. (2002). Understanding Molecular Simulation: From Algorithms to Applications. Academic Press.
• Allen, M. P., & Tildesley, D. J. (2017). Computer Simulation of Liquids. Oxford University Press.