Termostato de Nosé-Hoover

Grupo: Gabriel Azevedo, Rafael Abel e Thierre F. Conceição.

Termostato de Nosé-Hoover

O termostato de Nosé-Hoover é um algoritmo utilizado para simulação de dinâmica molecular. Este algoritmo utiliza um ensemble NVT, onde o número de partículas (N), o volume (V) e a temperatura (T) são mantidas constantes. Esse ensemble é relevante quando o sistema em estudo está em contato com um banho térmico.

A maneira que o algoritmo de Nosé-Hoover mantém a temperatura constante é a partir da adição de uma variável dinâmica fictícia (um "agente" externo), que atua sobre as velocidades das partículas no sistema, as acelerando ou desacelerando até que estas atinjam a temperatura desejada.

Método

Para entender o termostado de Nóse-Hoover, primeiramente será mostrado o termostato de Nosé.

Este termostato atribui coordenadas generalizados adicionais $s$ e o seu momento conjugado $p_{s}$ ao banho térmico. O fator $s$ é definido como um fator de escala das velocidades, onde:

${\dot {\mathbf {v}}}=s{\dot {\mathbf {r}}}=s{\mathbf {p}}/m$

E também são definidas as energia potenciais e cinética associadas a $s$ como:

${\mathcal {U}}_{s}=(N_{f}+1)k_{B}Tln(s)$ e ${\mathcal {K}}_{s}={\frac {1}{2}}Q{\dot {s}}^{2}={\frac {p_{s}^{2}}{2Q}}$

onde $Q$ é entendido como a "inércia térmica", ele determina a escala do tempo da flutuação de temperatura.

O Lagrangiano do sistema extendida (consistente das partículas e do banho térmico) então é postulado como:

${\mathcal {L}}={\mathcal {K}}+{\mathcal {K}}_{s}-{\mathcal {U}}-{\mathcal {U}}_{s}=\sum _{i}{\frac {{\mathbf {p}}_{i}^{2}}{2m_{i}s^{2}}}+{\frac {p_{s}^{2}}{2Q}}-{\mathcal {U}}({\mathbf {r}})-(N_{f}+1)k_{B}Tln(s)$

Como não é explicitamente dependente do tempo:

${\mathcal {H}}_{N}={\mathcal {K}}+{\mathcal {K}}_{s}+{\mathcal {U}}+{\mathcal {U}}_{s}=\sum _{i}{\frac {{\mathbf {p}}_{i}^{2}}{2m_{i}s^{2}}}+{\frac {p_{s}^{2}}{2Q}}+{\mathcal {U}}({\mathbf {r}})+(N_{f}+1)k_{B}Tln(s)$

Como ${\mathcal {H}}_{N}$ se conserva, esse sistema é numericamente estável