Expoente de Lyapunov para um Gás de Lennard–Jones - CBPFIndex

56 5.6. Estrutura do programa
2
1
0
−1
−2
2
1
0
−1
−2
2
1
0
−1
−2
U
4.2x10 5
K
K + U
U
4.2x10 5
K
K + U
U
4.2x10 5
4.4x10 5
4.4x10 5
4.4x10 5
4.6x10 5
4.6x10 5
4.6x10 5
Fase de Equilibração
4.8x10 5
t (passos)
Fase de Equilibração
4.8x10 5
t (passos)
Fase de Equilibração
4.8x10 5
t (passos)
5.0x10 5
5.0x10 5
5.0x10 5
Fase das Medições
5.2x10 5
Fase das Medições
5.2x10 5
Fase das Medições
5.2x10 5
ρ 0 = 0.01
5.4x10 5
ρ 0 = 0.05
5.4x10 5
ρ 0 = 0.50
5.4x10 5
Figura 5.4: Comportamento da energia durante as fases de equilibração e medições para três
valores de densidade. Os gráficos mostram a energia cinética K , energia potencial U
e a energia total por partícula, isto é, dividido por N = 108. A fase de equilibração
consistiu de 500 000 passos, sendo a velocidade ajustada a cada passo de acordo
com a equação (5.11). O ajuste da velocidade tem efeito sobre a energia cinética,
que permanece fixa durante este período e começa a flutuar quando o sistema entra
na fase das medições; o inverso ocorrendo com a energia total. Ver figura 5.9 mais
a frente onde é apresentado a continuação temporal destes três gráficos.