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
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
3. Método Estocástico 31<br />
forma, a média da hessiana se escreve:<br />
V =<br />
Ou, ainda:<br />
⎛<br />
⎜<br />
⎝<br />
Vq1q1<br />
Vq2q1<br />
.<br />
VqNq1<br />
⎞<br />
Vq1q2 · · · Vq1qN<br />
⎟<br />
Vq2q2 · · · ⎟ Vq2qN ⎟<br />
.. ⎟<br />
. . . ⎠<br />
<br />
VqNq2 · · · VqNqN<br />
=<br />
= 1<br />
3 Tr Vq1q1<br />
⎛ ⎞<br />
13 O · · · O<br />
⎜ ⎟<br />
⎜O<br />
13 · · · O⎟<br />
⎜<br />
⎝<br />
..<br />
⎟<br />
. . . . ⎠<br />
O O · · · 13<br />
+ 1<br />
3 Tr Vq1q2<br />
1<br />
V =<br />
3 Tr <br />
Vq1q1 13N + 1<br />
3 Tr <br />
Vq1q2 (Y3N − 13N) =<br />
= 1<br />
3 Tr <br />
Vq1q1 13N + 1<br />
3 Tr <br />
Vq1q2<br />
Y3N<br />
on<strong>de</strong> as matrizes Y3N e Y3N são <strong>de</strong>finidas como a seguir:<br />
Y3N =<br />
⎛ ⎞<br />
13 13 · · · 13<br />
⎜ ⎟<br />
⎜13<br />
13 · · · 13⎟<br />
⎜<br />
⎝<br />
..<br />
⎟<br />
. . . . ⎠<br />
13 13 · · · 13<br />
Notar que (Y3N) 2 = N Y3N , como conseqüência:<br />
e<br />
Y3N =<br />
<br />
Y3N<br />
2<br />
= (Y3N − 13N) 2 = (Y3N) 2 − 2Y3N + 13N =<br />
= (N − 2) Y3N + (N − 1)13N<br />
⎛ ⎞<br />
O 13 · · · 13<br />
⎜ ⎟<br />
⎜13<br />
O · · · 13⎟<br />
⎜<br />
⎝<br />
..<br />
⎟<br />
. . . . ⎠<br />
13 13 · · · O<br />
⎛ ⎞<br />
O 13 · · · 13<br />
⎜<br />
⎟<br />
⎜13<br />
O · · · 13⎟<br />
⎜<br />
⎝<br />
..<br />
⎟<br />
. . . . ⎠<br />
13 13 · · · O<br />
(3.41)<br />
(3.42)<br />
resultado que será usado mais adiante. Um outro resultado importante, que também será<br />
empregado em cálculos subseqüentes, é obtido da multiplicação da hessiana V com a ma-<br />
triz Y3N . Ao utilizarmos a invariância translacional <strong>de</strong> U, expressa através da equação (3.31),<br />
esta multiplicação resulta:<br />
VY3N = O<br />
Conseqüentemente:<br />
V Y3N = V (Y3N − 13N) = −V<br />
δV Y3N = V − V Y3N<br />
= V Y3N − V <br />
Y3N = −δV<br />
(3.43)<br />
on<strong>de</strong> a segunda equação acima advém <strong>de</strong> Y3N ser <strong>um</strong>a matriz constante, não sendo afetada