Grupo - Ising 2D: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
Linha 43: Linha 43:
== Teoria do Campo Médio: Uma abordagem aproximada ==
== Teoria do Campo Médio: Uma abordagem aproximada ==


O método do campo médio pode ser utilizado para introduzir algumas propriedades de um sistema de spins, assim como uma primeira análise de transições de fase. Porém seus resultados não são quantitativamente exatos, sendo necessária uma abordagem diferente ao problema para fins de resultados melhores.\\
O método do campo médio pode ser utilizado para introduzir algumas propriedades de um sistema de spins, assim como uma primeira análise de transições de fase. Porém seus resultados não são quantitativamente exatos, sendo necessária uma abordagem diferente ao problema para fins de resultados melhores (ver seção sobre o método Monte Carlo).


A magnetização do sistema está relacionada ao alinhamento de spin médio $<s_i>$
A magnetização do sistema está relacionada ao alinhamento de spin médio <math><s_i></math>. A magnetização total a uma temperatura <math>T</math> para um sistema de <math>N</math> spins é dada por
 
<math>
M = \sum_i <s_i> = N <s_i>,
</math>
 
 
Se adicionarmos um campo magnético ao problema, a função de energia do sistema se torna
 
<math>
E = -J \sum_{<ij>} s_i s_j - \mu H \sum_i s_i,
</math>
 
onde <math>H</math> representa o campo magnético e <math>\mu</math> o momento magnético associado com cada spin. Este campo faz com que os spins tendam a se orientar paralelamente a <math>H</math>, visto que isso diminui a energia. Para obter a aproximação de campo médio, consideramos que o sistema é constituído de um único spin <math>s_i</math>, de tal forma que a única energia envolvida é a energia de campo. As probabilidades de encontrar o sistema de um único spin nos seus dois possíveis estados são dadas por
 
<math>
P_+ = C e^{+\mu H/k_b T},
P_- = C e^{-\mu H/k_b T},
</math>
 
onde <math>C</math> é um coeficiente que pode ser determinado fazendo a condição de que as duas probabilidades se somem a 1. Portanto,
 
<math>
C = \frac{1}{e^{+\mu H/k_b T} + e^{-\mu H/k_b T}},
</math>
 
A média de <math>s_i</math> pode ser calculada por
 
<math>
<s_i> = \sum_{s_i = \pm 1} s_i P_{\pm} = P_+ - P_- = \tanh(\mu H/k_b T).
</math>


== O Método de Monte Carlo ==
== O Método de Monte Carlo ==

Edição das 21h38min de 20 de janeiro de 2018

Grupo: Ânderson Rosa, Caetano Pires e Lucas Doria.

sepa falar algo aqui tb

Introdução(?)

-talvez falar sobre materiais ferromagnéticos;

-falar sobre os conceitos de mec estatística necessários?;

-falar sobre o sistema de spins (geometricamente)?;

O Modelo de Ising

O modelo de Ising é construído a partir da assunção de que os spins do sistema apontam apenas na direção ou . Assim, o -ésimo spin do sistema pode assumir dois valores, que por conveniência são assumidos Cada um desses "Ising spins" interage com outros spins do sistema.

Em um material magnético real, a interação é maior entre spins mais próximos. Com essa motivação, uma forma de representar a interação entre os spins do sistema é levar em conta a interação apenas entre um spin e seus vizinhos mais próximos da cadeia de spins. A energia de tal sistema pode ser expressa por[1]

onde a soma se dá sobre todos os pares de spins mais próximos entre si, e é a constante de correlação, que assumimos positiva.

Uma análise qualitativa da expressão para a energia do microestado acima já mostra, por exemplo, que se dois spins são paralelos entre si, a energia de interação entre eles é . Se os spins são antiparalelos, então o produto dentro da soma é negativo, de forma que Portanto, as interações favorecem um alinhamento paralelo entre spins vizinhos.

Embora a energia do sistema seja menor quando todos os spins são paralelos entre si, é preciso considerar o efeito da temperatura sobre o sistema. No modelo estudado em questão, é considerado que o sistema se encontra em equilíbrio com uma fonte de temperatura , de forma que o comportamento do sistema pode ser estudado a partir do ensemble canônico.[1]

Para um sistema que se encontra em equilíbrio com uma fonte em temperatura , a probabilidade de encontrar o sistema em um estado particular é proporcional ao fator de Boltzmann [1]

onde é a energia do estado correspondente e a constante de Boltzmann. Cada um desses estados é uma configuração particular do conjunto de spins, chamados microestados do sistema. Portanto, se temos uma cadeia com Ising spins, o sistema possui microestados possíveis. Essa interação do sistema com uma fonte à temperatura faz com que o sistema passe por transições de um microestado para outro, fazendo com que spins individuais alternem entre +1 e -1 enquanto ganham ou perdem energia devido a fonte.

Uma medida macroscópica do momento magnético total do sistema é chamada de magnetização, e é uma média dos diversos microestados que o sistema visita durante uma medida. O momento magnético de um microestado é a soma dos valores dos spins daquele estado em particular. Assim, a magnetização medida é dada por

Teoria do Campo Médio: Uma abordagem aproximada

O método do campo médio pode ser utilizado para introduzir algumas propriedades de um sistema de spins, assim como uma primeira análise de transições de fase. Porém seus resultados não são quantitativamente exatos, sendo necessária uma abordagem diferente ao problema para fins de resultados melhores (ver seção sobre o método Monte Carlo).

A magnetização do sistema está relacionada ao alinhamento de spin médio . A magnetização total a uma temperatura para um sistema de spins é dada por


Se adicionarmos um campo magnético ao problema, a função de energia do sistema se torna

onde representa o campo magnético e o momento magnético associado com cada spin. Este campo faz com que os spins tendam a se orientar paralelamente a , visto que isso diminui a energia. Para obter a aproximação de campo médio, consideramos que o sistema é constituído de um único spin , de tal forma que a única energia envolvida é a energia de campo. As probabilidades de encontrar o sistema de um único spin nos seus dois possíveis estados são dadas por

onde é um coeficiente que pode ser determinado fazendo a condição de que as duas probabilidades se somem a 1. Portanto,

A média de pode ser calculada por

O Método de Monte Carlo

Transições de fase(?)

Conclusões e Observações

O modelo de Ising estudado neste trabalho é um modelo de spin extremamente simples. Outros modelos podem ser estudados. Por exemplo, podemos considerar os spins como sendo vetores de comprimento constante mas que tenham movimento de rotação em um plano[2], ou até mesmo considerar vetores em três dimensões[3]. Além disso, o alcance das interações entre os spins do sistema pode ser incrementada para os segundos, terceiros ou até mais distantes vizinhos mais próximos de um spin. Todos esses modelos têm sido estudados extensivamente. Apesar disso, o modelo simples em 2D com spins +1 e -1 estudado ainda representa bem as propriedades do sistema, principalmente o fenômeno de transição de fase.

Referências

  1. 1,0 1,1 N. J. Giordano, "Computational Physics". Department of Physics, Purdue University. Upper Saddle River, New Jersey. Prentice-Hall, 1997.