Algoritmo de Wang-Landau

De Física Computacional
Revisão de 19h31min de 28 de novembro de 2021 por Rafabel (discussão | contribs) (→‎Introdução)
Ir para navegação Ir para pesquisar

Introdução

Simulações computacionais, como o método de Monte Carlo, são vastamente utilizadas para estudar transições de fase e fenômenos críticos. O método padrão para simulações de Monte Carlo é o algoritmo de Metropolis, entretanto, algoritmos novos e mais eficientes são usados em simulações modernas, como o algoritmo de Wang-Landau. Ao contrário dos métodos convencionais de Monte Carlo, que geram diretamente uma distribuição canônica a uma dada temperatura , a abordagem de Wang-Landau estima a densidade de estados diretamente por meio de um passeio aleatório, que produz um histograma plano no espaço de energia . -----------REF WANG LANDAU-------- Mesmo para modelos que podem ser resolvidos analiticamente, a densidade de estados não pode ser determinada para sistemas maiores [1]


Amostragem de Wang-Landau

No início da simulação, é desconhecido e fazemos uma estimativa inicial para ele. A abordagem mais simples é definir para todas as energias possíveis . A configuração de spin inicial para toda a rede pode ser escolhida arbitrariamente. Então, uma caminhada aleatória no espaço de energia é iniciada pela formação de estados de teste, cada um dos quais é produzido escolhendo aleatoriamente um spin e alterando aleatoriamente seu estado.

Cada vez que uma energia é visitada, a entrada correspondente em é incrementada em 1. A estimativa de é então modificada por um fator multiplicativo , e o valor atualizado realiza um passeio aleatório adicional no espaço de .

Se e são energias antes e depois de um valor de spin ser alterado, a probabilidade de transição da energia para é

A razão das probabilidades de transição de para e de a podem ser calculados como

Logo, o algoritmo de passeio aleatório satisfaz o equilíbrio detalhado:

onde é a probabilidade na energia e é a probabilidade de transição de para .

Algoritmo

Resumindo, o passo a passo do algoritmo pode ser escrito como:

1. Seto e um fator de modificação ;

2. Aleatoriamente, flipo um spin com probabilidade ;

3. Modifico a densidade de estados e atualizo o histograma ;

4. Continuo até o histograma estar reto, então diminuo o valor de e reseto o valor de ;

5. Repito 2-4 até .

Aplicação ao Modelo de Ising 2D

Modelo de Ising

Uma rede 2D que consiste 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 ou .

O hamiltoniano pode ser calculado por

A soma ocorre sobre todos sítios vizinhos e a interação é ferromagnética para e antiferromagnética para .

Implementação

Resultados

Fig1.png Fig2.png Fig3.png

  1. Exact Distribution of Energies in the Two-Dimensional Ising Model, P. D. Beale, Phys. Rev. Lett. 76,78, Jan. 1996. https://doi.org/10.1103/PhysRevLett.76.78