17.04.2013 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!