Modelo de Ising

Autores: André Guimarães, Filssen Schereiber, João Roth e Lucas Oliveira

Este trabalho tem como objetivo estudar o Modelo de Ising...

Introdução

O modelo de Ising surgiu na década de 1920 como uma simplificação do modelo de Heisenberg. Inicialmente desenvolvido para compreender o ferromagnetismo, passou, ao longo dos tempo, a ser aplicado em diversas áreas interdisciplinares, tais como biologia, neurociência, dinâmica social e dinâmica, devido a sua simplicidade.

Este é, provavelmente, o modelo mais comumente estudado na mecânica estatística. É possível dizer que o modelo de Ising é, para a mecânica estatística, o que a mosca-da-fruta é para a genética. O modelo possibilita ilustrar o conceito essencial de como o equilíbrio entre energia e entropia conduz a uma mudança de fase. (falta a referência)

Foram realizados diversos estudos sobre matemática, física e as potenciais aplicações desse modelo. Neste trabalho, pretendemos explorar alguns aspectos fundamentais e os resultados possíveis de serem obtidos por meio de simulações de Monte Carlo em sistemas magnéticos com diversas dimensionalidades.

Problema com o modelo de Heisenberg

Em 1924, Wilhelm Lenz delegou a seu aluno de graduação Ernst Ising a tarefa de compreender e resolver o modelo de Heisenberg para o ferromagnetismo. Heisenberg propõe um modelo de rede que busca explicar o ferromagnetismo em sólidos a partir da energia de interação entre os spins. Por exemplo, podemos considerar uma rede quadrada com contornos periódicos onde os elétrons do sólido são distribuídos uniformemente. Dessa forma, podemos descrever a energia resultante da interação entre quaisquer dois elétrons como:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{ij} = -J_{ij}(\vec{S}_{i} \cdot \vec{S}_{j}) }

(1)

Na equação acima, Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle J_{ij}} representa uma constante de troca que denota a intensidade da interação $J$ entre dois spins vizinhos. Caso os spins não sejam vizinhos, esse coeficiente de interação é nulo. Além disso, os vetores Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{S}_{i}} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{S}_{j}} representam os momentos magnéticos. Portanto, através do produto escalar, a energia depende da projeção de um momento magnético em relação ao outro: será Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle +1} se estiverem alinhados e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle -1} se estiverem em sentidos contrários.

Vale ressaltar que o coeficiente J é positivo, o que implica em contribuições negativas na energia quando os vetores momento magnético são paralelos, devido ao sinal negativo da expressão. Por conseguinte, para spins antiparalelos, observa-se um aumento na energia. Como esperado, o sistema tende a permanecer no estado de menor energia. Assim, os spins do sistema tendem a se alinhar, como é característico de uma fase ferromagnética.

Esse sistema tem um Hamiltoniano da forma:

${\mathcal {H}}=-J\sum _{\langle ij\rangle }({\vec {S}}_{i}\cdot {\vec {S}}_{j})$

(2)

Na expressão, os símbolos de bra-ket são usados para representar a soma sobre os primeiros vizinhos. Contudo, desse sistema, temos uma função de partição da forma:

${\mathcal {Z}}=\sum _{E_{i}}e^{\frac {E_{i}}{K_{B}T}}$

(3)

Essa expressão representa o somatório de todas as configurações da exponencial da energia de cada configuração dividida pela constante de Boltzmann e pela temperatura.

Nesse momento, Ising se deparou com um problema, pois, além de possuir poucos recursos matemáticos na época, o número de configurações é infinito, uma vez que é possível observar qualquer direção entre os vetores momento magnético de forma contínua.

Modelo de Ising

Ising resolveu o problema das infinitas configurações assumindo que, ao invés de lidar com spins que rotacionam 360 graus nas três coordenadas espaciais, ele selecionaria uma coordenada e consideraria a projeção dos spins nela. Na prática, ele conseguiu reduzir o problema de infinitas configurações para apenas duas: ou o spin aponta para o sentido positivo do eixo, ou para o sentido negativo.

Dessa forma, o uso do produto escalar para os vetores de momento magnético não é mais necessário, uma vez que esses agora assumem valores entre Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle -1} e $+1$. Assim, podemos escrever o Hamiltoniano como:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{H} = -J \sum_{\langle ij \rangle} S_i S_j }

(4)

Com essas considerações, torna-se possível descrever a transição de fase entre ferromagnetismo e paramagnetismo em sistemas ao menos bidimensionais.

Além disso, podemos introduzir um termo no Hamiltoniano responsável por descrever o sistema na presença de um campo magnético externo:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{H} = -J \sum_{\langle ij \rangle} S_i S_j - H \sum_{i} S_i }

(5)

onde H é um coeficiente relacionado ao campo magnético externo. O acréscimo desse termo faz com que o estado de menor energia seja aquele favorável à direção desse campo.

Para o estudo da transição de fase magnética definimos uma grandeza chamada de magnetização. Essa grandeza é responsável para computar a média de todos os spins do sistema e, nesse caso, será o parâmetro de ordem. Seja Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle N} o número de spins do sistema, escrevemos a magnetização como:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle M = \frac{1}{N} \sum_{i} S_i }

(6)

ou ainda, é possível expressarmos a magnetização em função do número de spins para cima Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle N_{\uparrow}} ou para baixo Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle N_{\downarrow}} :

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle M = \frac{N_{\uparrow} - N_{\downarrow}}{N} }

(7)

Com essa notação, é possível observar um comportamento característico da magnetização. A magnetização será próxima de zero quando os números de spins para cima e para baixo forem próximos, ou seja, quando o sistema está em um estado de alta desordem. Por outro lado, para valores do módulo da magnetização próximos da unidade, temos um sistema bem ordenado, com praticamente todos os spins paralelos entre si.

Esse comportamento, conforme esperado do parâmetro de ordem, auxilia na distinção clara da transição de fase de um estado ordenado para um estado desordenado. Tal visualização pode ser realizada fixando alguns parâmetros físicos do sistema e analisando o comportamento da magnetização em relação à temperatura ou ao campo externo, por exemplo.