Termostato de Nosé-Hoover: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
Rafabel (discussão | contribs)
Sem resumo de edição
Rafabel (discussão | contribs)
Sem resumo de edição
Linha 10: Linha 10:
== Método ==
== Método ==


Para entender o termostado de Nóse-Hoover, primeiramente será mostrado o termostato de Nosé<ref name=nose> A molecular dynamics method for simulations in the canonical ensemble, Molecular Physics, 1984, Vol. 52, No. 2, 255-268 </ref>.  
Para entender o termostado de Nóse-Hoover, primeiramente será mostrado o termostato de Nosé<ref name=nose> NOSÉ, Shuichi, '''A molecular dynamics method for simulations in the canonical ensemble''', Molecular Physics, 1984, Vol. 52, No. 2, 255-268 </ref>.  


Este termostato atribui coordenadas generalizados adicionais <math> s </math> e o seu momento conjugado <math> p_s </math> ao banho térmico. O fator <math> s </math> é definido como um fator de escala das velocidades, onde:
Este termostato atribui coordenadas generalizados adicionais <math> s </math> e o seu momento conjugado <math> p_s </math> ao banho térmico. O fator <math> s </math> é definido como um fator de escala das velocidades, onde:
Linha 22: Linha 22:
<math> \mathcal K_s = \frac{1}{2}Q\dot s^2 = \frac{p_s^2}{2Q} </math>
<math> \mathcal K_s = \frac{1}{2}Q\dot s^2 = \frac{p_s^2}{2Q} </math>


onde <math> Q </math> é entendido como a "inércia térmica".
onde <math> Q </math> é 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:
 
<math> \mathcal L = \mathcal K + \mathcal K_s - \mathcal U - \mathcal U_s = \sum_i \frac{\bold p_i^2}{2m_is^2} + \frac{p_s^2}{2Q} - \mathcal U(\bold r) - (N_f + 1)k_BTln(s)</math>


== Resultados ==
== Resultados ==

Edição das 01h01min de 25 de maio de 2021

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[1].

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é[2].

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

𝐯˙=s𝐫˙=s𝐩/m

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

𝒰s=(Nf+1)kBTln(s) e 𝒦s=12Qs˙2=ps22Q

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:

=𝒦+𝒦s𝒰𝒰s=i𝐩i22mis2+ps22Q𝒰(𝐫)(Nf+1)kBTln(s)

Resultados

Programas Utilizados

Referências

  1. https://www2.ph.ed.ac.uk/~dmarendu/MVP/MVP03.pdf
  2. NOSÉ, Shuichi, A molecular dynamics method for simulations in the canonical ensemble, Molecular Physics, 1984, Vol. 52, No. 2, 255-268