Clusterização

De Física Computacional
Revisão de 13h55min de 28 de maio de 2021 por Mgteus (discussão | contribs)
Ir para navegação Ir para pesquisar

PÁGINA EM CONSTRUÇÃO

Clusterização do Modelo de Ising

Clusterização

Balanço Detalhado

Para respeitarmos o Balanço Detalhado, precisamos que a mudança da rede de um estado $ν para um estado μ ocorra com a mesma probabilidade da mudança de um estado μ para ν, denotamos essa mudança por: A(μ→ν)=A(ν→μ), com A(x→y) sendo a razão de aceitação da mudança de um estado x para um estado y.

Supondo que estamos mudando de um estado ν para outro estado μ, temos que a diferença de energia entre esses dois é resultado da quebra das ligações entre pares de spins orientados na mesma direção que não foram adicionados ao cluster, já que, não há uma garantia que a ida de ν→μ quebre a mesma quantidade de ligações que a volta de μ→ν. A probabilidade de não adicionarmos um spin vizinho ao cluster é dada por: 1−Padd; uma vez que Padd é a probabilidade de incluir esse spin no cluster.

Supondo que existam m ligações a serem quebradas na ida de ν→μ, a probabilidade desse evento é dada por (1−Padd)m. Porém, o mesmo pode não valer para a volta de μ→ν, em razão disso, precisamos analisar o caso em que não há o mesmo número de ligações a serem quebradas na volta e então a probabilidade será dada por (1−Padd)n com n sendo o número de ligações a serem quebradas de μ→ν.

Consideramos agora que Eν e Eμ sejam as energias associadas aos estados ν e μ, respectivamente, temos que: a cada m ligações que são quebradas de μ→ν a energia aumenta com +2J e para cada n novas ligações geradas de μ→ν a energia diminui com −2J. Pode-se escrever então que a diferença de energia entre μ e ν é dada por: Eν−Eμ=2mJ−2nJ=2J(m−n)

Seguindo a definição do Balanço Detalhado e impondo que o Processo Markoviano dessas mudanças de estados descritas acima respeite a Distribuição de Boltzmann, precisamos que: (1−Padd)mA(μ→ν)(1−Padd)nA(ν→μ)=e−β(Eν−Eμ)=(1−Padd)m−nA(μ→ν)A(ν→μ), tal que, A(μ→ν)A(ν→μ)=[e2βJ(1−Padd)]n−m⟹Padd=1−e−2βJ⟹(1−Padd)mA(μ→ν)(1−Padd)nA(ν→μ)=1

Algoritmo de Wolf

Dinâmica do Algoritmo