Expoente de Lyapunov para um Gás de Lennard–Jones - CBPFIndex
Expoente de Lyapunov para um Gás de Lennard–Jones - CBPFIndex
Expoente de Lyapunov para um Gás de Lennard–Jones - CBPFIndex
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
6. Dinâmica Molecular: Aplicação 81<br />
repetindo os resultados simulacionais <strong>para</strong> os parâmetros τ (k)<br />
c . Realizamos este procedimento<br />
pois, como discutimos na seção anterior, toda informação que dispomos referente aos tempos<br />
característicos τ (k)<br />
c é proveniente do ajuste gaussiano da função <strong>de</strong> correlação fc (τ) obtida<br />
a partir dos dados da simulação. Além da utilização dos mesmos valores <strong>para</strong> os tempos<br />
característicos, <strong>um</strong> outro fator com<strong>um</strong> entre os resultados simulação e teoria é que ambos<br />
foram obtidos através da solução do polinômio característico completo (3.62), diferentemente<br />
do resultado teoria aproximada.<br />
Figura 6.8: <strong>Expoente</strong> <strong>de</strong> <strong>Lyapunov</strong> em função da <strong>de</strong>nsida<strong>de</strong>. Simulação: todos os parâmetros<br />
obtidos através da simulação. Teoria: parâmetros µ e σ 2<br />
λ obtidos dos cálculos<br />
analíticos e os tempos τ (k)<br />
c , através da simulação. Ambos os gráficos construídos<br />
com a solução do polinômio característico (3.62) completo. Teoria aproximada:<br />
calculado com a aproximação do polinômio característico (3.64). Linha cheia correspon<strong>de</strong><br />
a λ ∝ ρ 1/3<br />
0 .<br />
Denominamos teoria aproximada os pontos obtidos <strong>de</strong> acordo com a aproximação <strong>para</strong> λ<br />
discutida no final do capítulo 3 cujo resultado repetimos na equação (6.17) mais acima.<br />
Através <strong>de</strong>sta aproximação, alternativa à solução do polinômio característico completo, o<br />
expoente <strong>de</strong> <strong>Lyapunov</strong> é expresso em função apenas dos parâmetros σ 2<br />
λ<br />
e τ(1)<br />
c . Os pon-<br />
tos correspon<strong>de</strong>ntes ao resultado teoria aproximada presentes na figura 6.8 foram calculados<br />
utilizando os resultados analíticos <strong>para</strong> o parâmetro σ 2<br />
λ<br />
e os resultados do ajuste gaussi-