Onde ''prog.c'' é o programa que deseja utilizar. É necessário que o programa e a biblioteca com funções para MC estejam na mesma pasta no momento da compilação. E execute da seguinte maneira
<source lang=sh>
<source lang=sh>
−
$ ./a.out TEMP SEED
+
$ ./a.out Q TEMP
</source>
</source>
−
onde o segundo termo é a temperatura do banho térmico e o terceiro é a semente utilizada pelo gerador de números pseudo-aleatórios. Os programas possuem uma diretiva de compilação para visualização do sistema ao decorrer da execução. Para utilizar é necessário ter o gnuplot <ref name=GNUPLOT>https://fiscomp.if.ufrgs.br/index.php/Gnuplot</ref> instalado e compilar da forma
+
onde o segundo termo é o estado Q do sistema e o terceiro é a temperatura do sistema. Os programas possuem uma diretiva de compilação para visualização do sistema ao decorrer da execução. Para utilizar é necessário ter o gnuplot <ref name=GNUPLOT>https://fiscomp.if.ufrgs.br/index.php/Gnuplot</ref> instalado e compilar da forma
<source lang=sh>
<source lang=sh>
−
$ gcc -DGNU prog.c -lm
+
$ gcc -DVIEW prog.c -lm
</source>
</source>
Linha 82: Linha 82:
<source lang=sh>
<source lang=sh>
−
$ ./a.out TEMP SEED | gnuplot
+
$ ./a.out Q TEMP | gnuplot
</source>
</source>
Linha 88: Linha 88:
+
[[Potts Metropolis]]
−
[[Biblioteca com funções em C para simulações de Monte Carlo]]
+
[[Potts Banho Térmico]]
+
[[Biblioteca com funções em C para simulações de Monte Carlo]]
+
[[Script em bash para gerar os dados]]
== Referências ==
== Referências ==
<references/>
<references/>
Edição das 18h36min de 16 de maio de 2021
Índice
1 Modelo de Potts
2 Método de Monte Carlo
2.1 Algorítmo de Metropolis-Hasting
2.2 Algoritmo de Banho Térmico
3 Programas Utilizados
4 Referências
Modelo de Potts
O "modelo de Potts de Q-estados" trata de um sistema de rede com N spins interagentes , onde um spin pode assumir valores discretos . Cada spin do sistema está limitado a interagir com outros spins em sua vizinhança e a energia da interação entre dois spins e é dada pelo potencial
onde é a função delta de Kronecker e é a constante de interação entre os spins. Dessa maneira, a interação entre dois spins vizinhos contabiliza um valor de energia ao sistema apenas se . A hamiltoniana do sistema é dada pela soma entre todas as interações entre spins vizinhos:
Este modelo é tido como uma generalização natural do Modelo de Ising e para o caso ambos modelos são equivalentes a menos de uma constante:
Nesse caso, a interação entre dois spins e assume a mesma dinâmica do modelo de Ising a contribuição para a energia do sistema será
Neste trabalho, o modelo de Potts foi estudado em uma rede quadrada 2D com vizinhança de von Neumann para primeiros vizinhos e condições de contorno periódicas. A quantidade de spins no modelo é com interações ferromagnéticas e , favorecendo vizinhanças de spins que compartilham o mesmo valor de para minimizar a energia do sistema.
Método de Monte Carlo
O método de Monte Carlo é aplicado ao modelo de Potts com o objetivo de gerar estados de equilíbrio para medir os observáveis do sistema. Os spins são iniciados com valores aleatórios de Q na rede e o método de Monte Carlo escolhe arbitrariamente um spin e gera um novo valor de para o spin. A partir disso, através de algum algoritmo específico, se escolhe como os estados serão gerados e quais serão aceitos ou não para o sistema transicionar, respeitando as condições de balanço detalhado e ergodicidade das cadeiras markovianas. Para este trabalho, foram estudados os algoritmos de Metropolis-Hasting e o algoritmo de banho térmico.
Algorítmo de Metropolis-Hasting
O primeiro algoritmo utilizado para gerar as configurações do sistema foi o algoritmo de Metropolis-Hasting. O algoritmo escolhe repetidamente um novo estado para o sistema e aceitando ou rejeitando ele de acordo com uma probabilidade de aceitação de transitar de um estado antigo para o novo estado . O algoritmo que iremos descrever utiliza a dinâmica de inversão única de spins.
Temos que a condição de balanceamento detalhado é dada por [1]:
onde é a diferença de energia entre o novo e o antigo estado.
Vamos supor que tenhamos os estados e e que temos a relação de energias: . Então, a maior das duas chances de aceitação é , portanto iremos igualar essa probabilidade a 1. Para que seja respeitada, iremos definir o valor de como . Temos, assim, o algoritmo de Metropolis-Hasting:
Dessa forma, sempre que tivermos um estado cuja energia seja menor do que a do estado atual, iremos aceitar a transição, mas se a energia for maior, teremos uma pequena probabilidade de trocarmos de estado.
Algoritmo de Banho Térmico
O algoritmo de Metropolis-Hasting para inversão única de spins é eficaz para o modelo de Potts em baixos valores de ou temperaturas acima da temperatura crítica, entretanto para valores altos de ou baixas temperaturas o algoritmo falha convergir o sistema rapidamente para a situação de equilíbrio.
Considerando um caso onde e um spin que possui 4 vizinhos, se todos os vizinhos do spin possuem valores diferentes uns do outro e do próprio spin, poderá levar em média passos de Monte Carlo para sortear um novo valor de que tem a transição aceita, e dessa forma o algoritmo irá demorar mais tempo para alcançar a configuração de equilíbrio do sistema. A dificuldade de aceitar transições é maior ainda para baixas temperaturas, onde a probabilidade de transicionar para um novo estado tem um peso maior para qualquer diferente dos spins da vizinhança, dessa maneira poderá demorar 96 passos para gerar um spin que seja igual a algum spin da vizinhança e realizar a transição. Para contornar este problema podemos utilizar o algoritmo de banho térmico.
O algoritmo de banho térmico, assim como o metropolis, é um algoritmo que mudaremos um spin por vez. O algoritmo segue as seguintes etapas: primeiro escolhemos um spin na rede ( ), e independente do seu valor atual, escolhemos um novo valor para . Esse novo valor será aceito, ou não, de acordo com a proporção dos pesos de Boltzmann. Temos que no algoritmo de Banho Térmico nós atribuímos um valor , entre 1 e , ao spin com uma probabilidade
Onde é a energia do sistema quando e o somatório é dado em todas energias possíveis. Como com esse algoritmo pode fazer que o spin assuma qualquer valor ele satisfaz a condição de ergodicidade. Temos que e , assim, pela descrição do algoritmo temos
Ou seja, o algoritmo de banho térmico respeita a condição de balanço detalhado.
Programas Utilizados
Programas na linguagem C. Para utilizar os programas, abra o terminal e compile da forma
$ gcc prog.c -lm
Onde prog.c é o programa que deseja utilizar. É necessário que o programa e a biblioteca com funções para MC estejam na mesma pasta no momento da compilação. E execute da seguinte maneira
$ ./a.out Q TEMP
onde o segundo termo é o estado Q do sistema e o terceiro é a temperatura do sistema. Os programas possuem uma diretiva de compilação para visualização do sistema ao decorrer da execução. Para utilizar é necessário ter o gnuplot [2] instalado e compilar da forma
$ gcc -DVIEW prog.c -lm
e então executar da maneira
$ ./a.out Q TEMP | gnuplot
Entretanto com a visualização no gnuplot o programa pode demorar mais para executar, então é recomendado diminuir os tempos de transiente (TRAN) e medidas (TMAX).
Potts Metropolis
Potts Banho Térmico
Biblioteca com funções em C para simulações de Monte Carlo
Script em bash para gerar os dados
Referências
↑ M. E. J. Newman, G. T. Barkema, "Monte Carlo Methods in Statistical Physics". Oxford University Press Inc., New York, 1999.