Algoritmo de Wang Landau: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
Sem resumo de edição
Sem resumo de edição
 
Linha 1: Linha 1:
== Wang Landau ==
== Modelo de Ising ==


Lorem ipsum dolor sit amet, consectetur adipiscing elit. Etiam porttitor augue vel blandit rhoncus. Vestibulum augue ipsum, interdum a mauris ac, finibus accumsan lorem. Vestibulum molestie efficitur nunc, a pellentesque nisi facilisis vitae. Donec pulvinar, mi in dapibus egestas, lacus orci rhoncus massa, hendrerit porttitor mauris nunc eu diam. Vestibulum scelerisque risus tellus, vel rutrum est consequat non. Nullam tincidunt vel leo vel aliquet. Donec efficitur sapien ac sem venenatis posuere. Suspendisse vulputate eget metus vitae faucibus. Etiam euismod vestibulum arcu eget fringilla. Proin vitae nibh nunc. Nulla facilisi. Proin molestie ligula et justo accumsan pellentesque. Donec at vehicula quam. Donec ut cursus mi, in pharetra nisi. Fusce finibus, ipsum at facilisis blandit, est quam mattis mi, vel suscipit ex quam et lacus. Suspendisse sit amet enim vitae erat iaculis tincidunt ut at mi.
Uma rede 2D que conssite de uma variável discreta em cada sítio que pode ser usada para representar o momento de dipolo magnético de um átomo
Cada sítio pode ter o valor de spin +1 ou -1.


Sed venenatis egestas eros. Pellentesque habitant morbi tristique senectus et netus et malesuada fames ac turpis egestas. Nunc volutpat id tellus a viverra. Nam fermentum lacinia tempor. Maecenas vel nulla ante. In quis varius velit, ut egestas tortor. Aliquam placerat interdum lacus quis convallis. Nam tempor id leo a congue. Integer vitae mattis est. Praesent sed egestas est.
O hamiltoniano pode ser calculado por
 
<math> H = -J \sum_{\langle ij \rangle}\sigma_i \sigma_j </math>
 
A soma ocorre sobre todos sítios vizinhos
 
Se considerarmos J>0 a interação é ferromagnética
 
J<0 é antiferromagnética
 
Uma maneira de estudar esse sistema é pelo algoritmo de Metropolis
 
1. Escolhe uma configuração inicial
 
2. Escolhe aleatoriamente um sitio
 
3. Flipa o spin desse sitio, recalcula a energia e a variação de energia
 
4. Gera um numero aleatorio 0 < r < 1
 
5. Se <math> r < e^{-\Delta E/k_BT} </math>, mantem o spin flipado, se não, volta
 
6. Volta para o passo 2
 
O modelo de ising também pode ser estudado pelo algoritmo de wang-landau
 
== Wang-Landau ==
 
Consideramos somente os 4 vizinhos mais proximos
 
Método para obter densidade de estados de um sistema.
 
O algoritmo de wang-landau faz uma random walk no espaço de energias
 
Ele obtem a densidade de estados como uma função da energia, g(E). Essa função é incrivelmente útil...
 
Vamos manter o g(E) e também o histograma de visitas para cada energia
 
1. setamos g(E) = 1 e um fator de modificação f=e
 
2. Aleatoriamente flipa um spin com probabilidade: <math> p(E_1 \to E_2) = min(g(E_1)/g(E_2), 1) </math>
 
3. Modifica a densidade de estados <math> g(E) \to g(E) \times f </math> e atualizamos o histograma
 
4. Continuamos até o histograma estiver reto, diminui o valor de f e reseta o histograma
 
5. Repito 2-4 até <math> \ln f \approx 1 </math>
 
=== Funções ===
 
Função de partição:
 
<math> Z = \sum g(E) e^{-E/k_BT} </math>
 
Energia interna:
 
<math> U(T) = \frac{\sum_E Eg(E) e^{-E/k_BT}}{\sum_E g(E) e^{-E/k_BT}} = \langle E \rangle </math>
 
Calor específico:
 
<math> C(T) = \frac{\partial U(T)}{\partial T} = \frac{\langle E^2 \rangle - \langle E \rangle ^2}{k_BT^2} </math>
 
Energia livre de Helmoltz:
 
<math> F(T) = -k_BT\ln(Z) = -k_BT\ln\left( \sum_E g(E) e^{-E/k_BT} \right) </math>
 
Entropia:
 
<math> S(T) = \frac{U(T) - F(T)}{T} </math>





Edição atual tal como às 21h01min de 21 de novembro de 2021

Modelo de Ising

Uma rede 2D que conssite de uma variável discreta em cada sítio que pode ser usada para representar o momento de dipolo magnético de um átomo Cada sítio pode ter o valor de spin +1 ou -1.

O hamiltoniano pode ser calculado por

A soma ocorre sobre todos sítios vizinhos

Se considerarmos J>0 a interação é ferromagnética

J<0 é antiferromagnética

Uma maneira de estudar esse sistema é pelo algoritmo de Metropolis

1. Escolhe uma configuração inicial

2. Escolhe aleatoriamente um sitio

3. Flipa o spin desse sitio, recalcula a energia e a variação de energia

4. Gera um numero aleatorio 0 < r < 1

5. Se , mantem o spin flipado, se não, volta

6. Volta para o passo 2

O modelo de ising também pode ser estudado pelo algoritmo de wang-landau

Wang-Landau

Consideramos somente os 4 vizinhos mais proximos

Método para obter densidade de estados de um sistema.

O algoritmo de wang-landau faz uma random walk no espaço de energias

Ele obtem a densidade de estados como uma função da energia, g(E). Essa função é incrivelmente útil...

Vamos manter o g(E) e também o histograma de visitas para cada energia

1. setamos g(E) = 1 e um fator de modificação f=e

2. Aleatoriamente flipa um spin com probabilidade:

3. Modifica a densidade de estados e atualizamos o histograma

4. Continuamos até o histograma estiver reto, diminui o valor de f e reseta o histograma

5. Repito 2-4 até

Funções

Função de partição:

Energia interna:

Calor específico:

Energia livre de Helmoltz:

Entropia:


Método

Subtítulo