Clusterização: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
Mgteus (discussão | contribs)
Sem resumo de edição
Mgteus (discussão | contribs)
Linha 18: Linha 18:


=== Dinâmica do Algoritmo ===
=== Dinâmica do Algoritmo ===
[[File:wiki.png|thumb|Wikipedia logo]]

Edição das 14h00min de 28 de maio de 2021

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(xy) 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: 1Padd; 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 (1Padd)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 (1Padd)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μ=2mJ2nJ=2J(mn)

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: (1Padd)mA(μν)(1Padd)nA(νμ)=eβ(EνEμ)=(1Padd)mnA(μν)A(νμ), tal que, A(μν)A(νμ)=[e2βJ(1Padd)]nmPadd=1e2βJ(1Padd)mA(μν)(1Padd)nA(νμ)=1

Algoritmo de Wolf

Dinâmica do Algoritmo

Arquivo:Wiki.png
Wikipedia logo