EQUAÇÃO GERAL DE GRACELI.[quantização de Graceli].

  G ψ = E ψ =  E [G+].... ..  =

G ψ = E ψ =  E [G+ψ ω /c] =   [/ ] /  /   = ħω [Ϡ ]  [ξ ] [,ς]   [ q G*]ψ μ / h/c ψ(x, t)  [x  t ]..

[ q G*] ==G ψ = E ψ =  E [G+].... .. 

SISTEMA GRACELI DE:

 TENSOR G+ GRACELI = SDCTIE GRACELI, DENSIDADE DE CARGA E DISTRIBUIÇÃO ELETRÔNICA, NÍVEIS DE ENERGIA, NÚMERO E ESTADO QUÂNTICO. + POTENCIAL DE SALTO QUÂNTICO RELATIVO AOS ELEMENTOS QUÍMICO COM O SEU RESPECTIVO  E ESPECÍFICO NÍVEL DE ENERGIA., POTENCIAL DE ENERGIA, POTENCIAL QUÍMICO,  SISTEMA GRACELI DO INFINITO DIMENSIONAL.

ONDE A CONFIGURAÇÃO ELETRÔNICA TAMBÉM PASSA A SER DIMENSÕES FÍSICO-QUÍMICA DE GRACELI. 

[ q G*] = energia quântica Graceli.

O modelo de Ising é muito empregado na física do estado sólido^[1] e consiste em um modelo teórico que descreve de forma clara as propriedades magnéticas, bem como, busca simular a estrutura de uma substância ferromagnética.^[2] Obtido através de uma simplificação do modelo de Heisenberg, este é tido como um modelo dinâmico de transição de fase que incorpora interações de curto alcance numa rede d-dimensional.^[3]

Mesmo sendo conhecido como um modelo simples, uma vez que considera apenas as interações entre os spins dos vizinhos mais próximos, o modelo de Ising demonstra resultados compensatórios,^[2] podendo descrever qualitativamente vários fenômenos relacionados às transições de fase. Sua usabilidade vem do fato de que as características essenciais dos fenômenos de transição (especialmente na vizinhança próxima do ponto crítico), não serão significativamente alteradas com a inclusão das interações dos spins mais distantes.^[4]

Do mesmo modo, a substituição de uma rede estrutural por outra, desde que sua dimensionalidade seja mantida, não irá alterar significativamente as características dos fenômenos de transição de fase. Por fim, essas soluções podem ser compartilhadas, com poucas modificações, por muitos outros sistemas físicos passando por tipos muito diferentes de transições (gás-líquido ao invés de paramagnético-ferromagnético, por exemplo).^[4] Devido a essa relação íntima entre as transições de fase e o modelo de Ising, serão inicialmente abordadas as transições de fase, tão importantes para toda a física, e posteriormente o modelo de Ising propriamente dito.

Transições de fase

Os fenômenos físicos aos quais o formalismo da mecânica estatística é aplicado podem, em geral, ser classificados em duas categorias. Na primeira categoria, os elementos microscópicos do sistema dado são considerados como não interagentes.^[4] Para essa situação, as funções termodinâmicas do sistema devem estar diretamente relacionadas ao conhecimento dos níveis de energia dos seus constituintes individuais.

A segunda categoria, por sua vez, é formada por sistemas nos quais a interação interatômica desempenha um importante papel físico. Particularmente em sólidos, desde que as posições dos átomos ao longo de uma grande faixa de temperaturas se mantenham próximas aos seus valores médios, é possível tratar o sólido escolhido como um “conjunto de osciladores harmônicos praticamente não interagentes”.^[4]

Nesta situação, os níveis de energia do sistema total não irão estar relacionados aos níveis de energia dos elementos individuais. Em vez disso, constata-se que, em circunstâncias favoráveis, um grande número de constituintes microscópicos do sistema podem exibir uma tendência de interagir uns com os outros de uma maneira bastante forte e cooperativa.^[4] Este comportamento, assume significado macroscópico a uma temperatura comumente chamada de temperatura crítica do sistema (${\displaystyle T_{c}}$).

Quando abordados os fenômenos pertencentes a essa segunda categoria é comum perceber-se uma mudança repentina do comportamento das propriedades termodinâmicas do sistema dado, ou seja, transições de fase. Assim como os fenômenos críticos, as transições de fases costumam ocorrer nos mais diversos sistemas: fluidos simples e misturass, materiais magnéticos e ferroelétricos, superfluidos, supercondutores, etc.^[3] Os sistemas físico-químicos que sofrem transições de fase podem ser representados em vários graus de precisão.

As teorias clássicas das mudanças de fase em sólidos e líquidos envolvem problemas matemáticos admiráveis,^[5] que vêm sendo estudados desde 1873, quando van der Waals publicou sua tese de doutorado acerca do comportamento crítico em sistemas fluidos. Pierre Curie e Pierre Weiss, por sua vez, se dedicaram ao estudo da transição para o ferromagnetismo. Teorias dessa natureza, envolvem experiências nas vizinhanças dos pontos críticos e têm sido usadas para descrever os aspectos qualitativos de vários tipos de transições de fases.^[3]

Na região crítica (que envolve o ponto crítico), grandezas termodinâmicas como os calores específicos, a compressibilidade e a susceptibilidade magnética, apresentam um comportamento peculiar caracterizado por divergências assintóticas. Esse comportamento, por sua vez, é identificado através dos expoentes críticos.

A partir da década de 60, quando o avanço das técnicas e tecnologias permitiu que as teorias clássicas passassem por um processo mais rigoroso de análise, percebeu-se através das experiências nas vizinhanças da criticalidade que os expoentes críticos assumem valores universais, bem definidos, mas que não coincidem com as previsões das “teorias clássicas”.  Assim sendo, enquanto as teorias clássicas para sistemas distintos afirmavam que estes expoentes possuíam um caráter universal (um conjunto clássico de valores), tanto os dados experimentais quanto diversos resultados teóricos apontavam para a existência de classes de universalidade (definidas por parâmetros como a dimensionalidade dos sistemas sob consideração).^[3]

Atualmente se reconhece que os expoentes críticos são determinados por pouquíssimos fatores, entre os quais estão:^[3]

1. A dimensionalidade dos sistemas físicos;
2. A dimensionalidade do parâmetro de ordem; e
3. O alcance das interações microscópicas (nos sistemas de interesse físico as interações microscópicas serão, em geral, de curto alcance).

Tendo em vista essas questões, será discutido agora um dos mais simples e mais estudados modelos da mecânica estatística:^[6]: o modelo de Ising.

Definição

Um dos fenômenos mais interessantes da física do estado sólido é o ferromagnetismo.^[1] Do mesmo modo, a magnetização espontânea dos materiais ferromagnéticos constituiu por muito tempo um problema intransponível para a física como ciência.

As primeiras teorias clássicas que tentaram explicar este problema surgiram em meados do século XIX,^[7] no entanto seu desenvolvimento completo se deu apenas com o surgimento da mecânica quântica e do conceito de spin. A partir deste momento, para explicar o ferromagnetismo, tornou-se interessante olhar para as características microscópicas do sistema e modelar o comportamento de alguns dos seus elementos.

O modelo de Ising, primeiramente proposto por Wilhelm Lenz em 1920, com o objetivo de estudar o magnetismo em materiais, foi resolvido para uma dimensão pelo seu aluno de doutorado, o estudante Ernst Ising, em 1925.^[2] Embora Ernst tenha obtido com o modelo de Ising resultados interessantes acerca das propriedades termodinâmicas, ele pode perceber que o caso em uma dimensão não possuía transição de fase ou magnetização espontânea.^[8]

A demonstração da transição de fase ocorreu, por sua vez, em 1944, quando Lars Onsager resolveu por aproximação o modelo bidimensional de Ising na ausência de um campo magnético.^[8] O cálculo por aproximação, feito para funções termodinâmicas com descontinuidades é bastante complicado visto que, a convergência da aproximação por funções analíticas em tais casos é lenta.^[5] Sem nenhuma aproximação, a solução analítica tanto para duas como para três dimensões seria impossível.^[8]

O modelo de Ising bidimensional cede a um tratamento exato na mecânica estatística. Esse é o único exemplo não trivial de uma transição de fase que pode ser trabalhada com rigor matemático.^[3] O nome modelo de Ising foi popularizado após o artigo de R. Peierls, publicado em 1936 com o título “On Ising’s Model of Ferromagnetism”.^[8]

A proposta de explicação que o modelo de Ising usa para a magnetização espontânea é baseada no conceito de interação de troca.^[9] Para modelarmos essa situação teremos que considerar então os termos referentes à energia no estado fundamental, àqueles relativos à interação coulombiana (interação entre partículas eletricamente carregadas) além, é claro, do termo associado à interação de troca.

A hamiltoniana de Heisenberg, assim como segue, é a responsável maior por definir a interação de troca entre os elétrons dando origem ao já citado modelo isotrópico de Heisenberg.

${\displaystyle H=-J\textstyle \sum _{i,j}\displaystyle (\sigma _{i}^{x}\sigma _{j}^{x}+\sigma _{i}^{y}\sigma _{j}^{y}+\sigma _{i}^{z}\sigma _{j}^{z})}$ /

/ G ψ = E ψ =  E [G+].... .. 

Na equação acima J é a interação de troca entre os spins ${\displaystyle \sigma _{i}}$ e ${\displaystyle \sigma _{j}}$, bem como x, y e z são as direções das suas componentes.^[9] A hamiltoniana de Ising, por sua vez, surge de uma particularização da hamiltoniana de Heisenberg, ou seja, para descrever um sistema cujos spins apontam para uma única direção (direção z, por exemplo). Veja a seguir a hamiltoniana de Ising.

${\displaystyle H=-J\textstyle \sum _{i,j}\displaystyle \sigma _{i}^{z}\sigma _{j}^{z}}$ / G ψ = E ψ =  E [G+].... .. 

Por ter uma posição favorecida, o modelo de Ising é dito anisotrópico, ou seja, suas propriedades físicas devem variar com a direção.

Aproximação de ordem zero

A fim de descrever a teoria microscópica das transições de fase foram desenvolvidos alguns modelos matemáticos. O modelo de Ising, como citado anteriormente, busca descrever tais transições, para tanto considera as interações de curto alcance entre os spins dos átomos imediatamente mais próximos (chamados primeiros vizinhos) de um átomo específico em uma rede d-dimensional.

Algumas técnicas foram desenvolvidas com o intuito de encontrar a solução para o modelo de Ising, essas técnicas formam o pequeno conjunto de ferramentas disponíveis para para descrever sistemas complexos. Duas técnicas podem ser destacadas, são elas a aproximação de ordem zero e a primeira aproximação.^[3]

A partir das tentativas de desenvolver estudos estatísticos de transição de fase do tipo ordem-desordem, Gorsky em 1928 supôs que o trabalho gasto para mudar um átomo de uma posição ordenada para uma desordenada era diretamente proporcional ao grau de ordem do sistema.^[4] A ideia foi desenvolvida por Bragg e Williams em 1934 e 1935 respectivamente, que pela primeira vez introduziram a ideia de parâmetro de longo alcance.

Nessa aproximação, a ideia central consiste em isolar um spin da rede e considerar a existência de um valor fixo para a magnetização dos demais spins, de modo que é possível supor a existência de um campo magnético interno médio, que interage com o spin isolado, permitindo assim desprezar as interações entre os pares de spins vizinhos e as flutuações causadas na função de correlação, que se estendem fora da célula fixada da rede. Desta forma a aproximação de ordem zero considera somente as flutuações que ocorrem na célula fixada, ou seja, envolve somente as flutuações que ocorrem com um único spin fixo.^[3][10]

Ao fixar um único spin da rede e considerar a interação deste com os demais, a partir de um campo médio de magnetização, o problema de mecânica estatística de muitos corpos é reduzido para um problema de um único corpo. Essa aproximação, por sua vez, carrega certo grau de inexatidão, especialmente quando são envolvidos os sistemas próximos do ponto crítico, em que as flutuações se estendem a grandes distâncias. Mas, apesar dessa aproximação não fornecer resultados exatos, ela pode ser muito útil, como por exemplo, na descrição de propriedades termodinâmicas de sistemas em regiões de altas temperaturas, onde não é esperado a existência de correlação entre o spin e seus primeiros vizinhos.^[10] Existem ainda algumas versões mais sofisticadas do modelo da aproximação de ordem zero, que consistem em fixar mais de um spin da rede, a exemplo de um par de spins vizinhos, reduzindo o problema de vários corpos a um problema de dois corpos e fornecendo resultados mais precisos.^[10]

Em um modelo de spin ½ , para uma rede com N spins, N representa o número total de spins dado pela soma de todos os spins Up da rede (${\displaystyle N_{+}}$) com todos spins Down (${\displaystyle N_{-}}$).^[4] Defini-se o parâmetro de ordem de longo alcance (${\displaystyle L}$) nesse sistema na forma,

${\displaystyle L=\left({\frac {1}{N}}\right)\sum _{i}\sigma _{i}}$ / G ψ = E ψ =  E [G+].... .. 

note que ${\displaystyle L}$ é um parâmetro médio e representa a média da correlação de spin (${\displaystyle {\overline {\sigma _{i}}}}$). A equação acima pode ser escrita como:

${\displaystyle L=\left({\frac {N_{+}-N_{-}}{N}}\right).}$/ G ψ = E ψ =  E [G+].... .. 

A magnetização por sua vez, é expressa como o produto da diferença entre a quantidade de spins Up e Down e pelo momento magnético (${\displaystyle \mu }$), matematicamente tem-se que:

${\displaystyle M=(N_{+}-N_{-})\mu ,}$ / G ψ = E ψ =  E [G+].... .. 

que ainda, pode ser reescrita na forma:

${\displaystyle M=N.L.\mu }$ / G ψ = E ψ =  E [G+].... .. 

evidenciando a dependência que a magnetização possui em relação ao parâmetro de ordem de longo alcance.

O hamiltoniano deste sistema, que possui a forma dada por:

${\displaystyle H\{\sigma _{i}\}=-J\textstyle \sum _{nm}\displaystyle \sigma _{i}\sigma _{j}-\mu B\textstyle \sum _{i}\displaystyle \sigma _{i}}$ / G ψ = E ψ =  E [G+].... .. 

pode ser reescrito, ficando da seguinte maneira:

${\displaystyle H\{\sigma _{i}\}=-{\frac {jq{\overline {L}}}{2}}\textstyle \sum _{i}\displaystyle \sigma _{i}-\mu B\textstyle \sum _{i}\displaystyle \sigma _{i}.}$ / G ψ = E ψ =  E [G+].... .. 

Nesta equação, q representa o número de coordenação, isto é, a quantidade de vizinhos imediatamente  mais próximos do spin fixado (em uma rede unidimensional q = 2 e em uma rede bidimensional q = 4) e B representa o campo magnético externo.

O hamiltoniano pode ser resolvido, de forma a denotar o valor da energia do sistema,  que vale:

${\displaystyle E=-{\frac {NqJ{\overline {L}}}{2}}L-\mu BNL.}$/ G ψ = E ψ =  E [G+].... .. 

De modo que a energia média para esse sistema é dada de forma semelhante:

${\displaystyle E_{(}med)=U=-{\frac {NqJ{\overline {L}}^{2}}{2}}-\mu BN{\overline {L}}.}$/ G ψ = E ψ =  E [G+].... .. 

A diferença de energia (${\displaystyle \Delta \varepsilon }$) para trocar um spin Up por um Down, por sua vez, é descrita como:

${\displaystyle \Delta \varepsilon \approx 2\mu {\Biggl (}{\frac {qJ\sigma }{\mu }}+B{\Biggr )}.}$/ G ψ = E ψ =  E [G+].... .. 

Fazendo ${\displaystyle B\prime ={\frac {qJ\sigma }{\mu }}}$, a equação passa a ter a seguinte forma:

${\displaystyle \Delta \varepsilon \approx 2\mu {\Biggl (}B\prime +B{\Biggr )}}$/ G ψ = E ψ =  E [G+].... .. 

onde ${\displaystyle B\prime }$ representa o campo magnético interno do sistema.

As condições para mudança de fase são determinadas a partir da relação a seguir,

${\displaystyle {\frac {\overline {N_{-}}}{\overline {N_{+}}}}=\exp({\frac {-\Delta \varepsilon }{kT}}).}$/ G ψ = E ψ =  E [G+].... .. 

Quando realizados os tratamentos matemáticos adequados é obtida a seguinte relação,

${\displaystyle {\overline {L_{o}}}=\tanh {\Biggl (}{\frac {qJ{\overline {L}}+\mu B}{kT}}{\Biggr )}}$/ G ψ = E ψ =  E [G+].... .. 

que, na ausência do campo magnético externo se reduz a

${\displaystyle {\overline {L_{o}}}=\tanh {\Biggl (}{\frac {qJ{\overline {L_{o}}}}{kT}}{\Biggr )}.}$/ G ψ = E ψ =  E [G+].... .. 

A equação acima pode ser resolvida graficamente, de modo que sua solução indica uma temperatura crítica (${\displaystyle T_{c}}$), a partir da qual é possível definir o início da transição de fase ferromagnética. Essa temperatura recebe o nome de temperatura de Curie, que pode ser expressa matematicamente por,

${\displaystyle T_{c}={\frac {qJ}{k}}.}$/ G ψ = E ψ =  E [G+].... .. 

Quando o sistema sai de um estado com temperaturas maiores que ${\displaystyle T_{c}}$ e passa para um estado com temperatura menores, ocorre uma transição de fase. Durante este fenômeno os spins da rede passam por um processo chamado de acoplamento, se alinhando na mesma direção e sentido, de modo que o sistema passa a apresentar magnetização espontânea.^[4]

O calor específico para este sistema, por sua vez, pode ser expresso da seguinte maneira,

${\displaystyle C_{o}(T)\approx 4Nk{\Biggl (}{\frac {T_{c}}{T}}{\Biggr )}^{2}\exp {\Biggl (}{\frac {-2T_{c}}{T}}{\Biggr )}.}$/ G ψ = E ψ =  E [G+].... .. 

Por intermédio dessa equação é possível notar que quando ${\displaystyle T\rightarrow T_{c}}$ o calor específico do sistema desaparece, isto ocorre pois para temperaturas elevadas não há magnetização espontânea.

Aproximação de primeira ordem

Aglomerado com spin central e seus spins vizinhos.

As aproximações estabelecidas para o modelo de Ising de ordem zero são bastante genéricas, podendo ser aprimoradas quando considera-se a correlação de um dado spin entre os vizinhos próximos. Essas ponderações definem a chamada aproximação de Bethe, desenvolvida por Bethe e Rushbrooke,^[4] também conhecida como aproximação de primeira ordem do modelo de Ising, que considera que um dado spin, tem q vizinhos próximos, e cada ponto da rede é também associado aos seus respectivos spins (${\displaystyle \sigma _{0},...,\sigma _{j},...,}$)^[4], conforme ao lado, levando em conta as trocas de energia (${\displaystyle J_{0i},...,J_{ij},...,}$) entre seus pares, através de suas interações.

Para esse modelo considera-se ${\displaystyle J_{iJ}}$ desprezível,^[4] exceto para seus pares mais próximos, já que essa interação fica menos evidente conforme aumentam-se as distâncias entre o spin central (${\displaystyle \sigma _{0}}$) e os demais (${\displaystyle \sigma _{j}}$). Essa observação permite que essa troca de energia, no caso de interações isotrópicas, seja definida simplesmente pelo símbolo J, tendo em vista que para J>0 tem-se o caso do ferromagnetismo (caso de interesse) e para J<0 o antiferromagnetismo.^[3]

O hamiltoniano deste sistema, onde o spin central interage a curto alcance com seus primeiros vizinhos é dado por,

${\displaystyle H_{q+1}=-\mu B\sigma _{0}-\mu (B+B')\sum _{j=1}^{q}\sigma _{j}-J\sum _{j=1}^{q}\sigma _{0}\sigma _{j}}$/ G ψ = E ψ =  E [G+].... .. 

onde, ${\displaystyle B}$ representa o campo magnético externo atuando sobre a rede de spins e ${\displaystyle B'}$ o campo magnético interno, que obedece ao princípio da auto consistência.^[3] Esse princípio está relacionado a afirmação de que ambas expressões devem ser assintoticamente equivalentes impondo que o valor médio de spin central da rede  ${\displaystyle {\overline {\sigma _{0}}}}$, seja igual ao valor médio de spin ${\displaystyle {\overline {\sigma _{J}}}}$, de qualquer um dos ${\displaystyle q}$ vizinhos^[4]. Assim a função partição (${\displaystyle Z}$) do sistema fica:

${\displaystyle Z=\sum _{\sigma _{0},\sigma _{j}=\pm 1}exp{\bigg [}{\frac {1}{kT}}{\bigg \{}\mu B\sigma _{0}+\mu (B+B')\sum _{j=1}^{q}\sigma _{j}+J\sum _{j=1}^{q}\sigma _{0}\sigma _{j}{\bigg \}}{\bigg ]}}$/ G ψ = E ψ =  E [G+].... .. 

de modo que considerando as seguintes variáveis

${\displaystyle \alpha ={\frac {\mu B}{kT}}}$, ${\displaystyle \alpha '={\frac {\mu B'}{kT}}}$ e ${\displaystyle \gamma ={\frac {J}{kT}}}$/ G ψ = E ψ =  E [G+].... .. 

obtém-se

${\displaystyle Z=\sum _{\sigma _{0},\sigma _{j}=\pm 1}exp{\bigg [}\alpha \sigma _{0}+(\alpha +\alpha ')\sum _{j=1}^{q}\sigma _{j}+\gamma \sum _{j=1}^{q}\sigma _{0}\sigma _{j}{\bigg ]}}$/ G ψ = E ψ =  E [G+].... .. 

A função partição canônica, pode ser escrita como o somatório de dois termos, pertencentes a ${\displaystyle \sigma _{0}=1}$ (spin Up) e ${\displaystyle \sigma _{0}=-1}$ (spin Down),

${\displaystyle Z=Z_{+}+Z_{-}}$

o que resulta em

${\displaystyle Z_{\pm }=\sum _{\sigma _{j}=\pm 1}exp{\bigg [}\pm \alpha +(\alpha +\alpha '\pm \gamma )\sum _{j=1}^{q}\sigma _{j}{\bigg ]}}$,/ G ψ = E ψ =  E [G+].... .. 

${\displaystyle Z_{\pm }=e^{\pm \alpha }[2cosh(\alpha +\alpha '\pm \gamma )]^{q}}$.

O valor médio do spin central é então dado por

${\displaystyle {\overline {\sigma _{0}}}={\frac {Z_{+}-Z_{-}}{Z}}}$/ G ψ = E ψ =  E [G+].... .. 

e, conforme segue, os spins dos ${\displaystyle q}$ vizinhos ficam

${\displaystyle {\overline {\sigma _{0}}}={\frac {1}{Z}}{\big \{}Z_{+}tanh(\alpha +\alpha '+\gamma )+Z_{-}tanh(\alpha +\alpha '+\gamma ){\big \}}}$/ G ψ = E ψ =  E [G+].... .. 

que, considerando a condição de autoconsistência do campo B’^[3] será

${\displaystyle {\overline {\sigma _{0}}}={\overline {\sigma _{j}}}={\overline {\sigma }}}$/ G ψ = E ψ =  E [G+].... .. 

Substituindo ${\displaystyle Z_{+}}$ e ${\displaystyle Z_{-}}$ nas expressões obtidas anteriormente para a função partição, na expressão acima para ${\displaystyle \sigma }$,

${\displaystyle Z_{+}=e^{+\alpha }[2cosh(\alpha +\alpha '+\gamma )]^{q}}$

${\displaystyle Z_{-}=e^{-\alpha }[2cosh(\alpha +\alpha '-\gamma )]^{q}}$

Temos a seguinte relação

${\displaystyle e^{2\alpha '}={\bigg \{}{\frac {cosh(\alpha +\alpha '+\gamma )}{cosh(\alpha +\alpha '-\gamma )}}{\bigg \}}^{q-1}}$/ G ψ = E ψ =  E [G+].... .. 

Esta expressão leva a equação de estado desta aproximação e representa o comportamento magnético da rede de spins. Pode-se obter algumas informações quando levado em conta as variáveis ${\displaystyle \alpha }$, ${\displaystyle \alpha '}$ e ${\displaystyle \gamma }$, que estão diretamente relacionadas à temperatura ${\displaystyle T}$, campo externo ${\displaystyle B}$, campo interno ${\displaystyle B'}$ e a energia de interação J.^[3]

Considerando o campo magnético externo ${\displaystyle B=0}$, tem-se que ${\displaystyle \alpha =0}$ ^[4] na expressão acima, que resultará na condição de magnetização espontânea. A equação de estado se torna então:

${\displaystyle \alpha '={\frac {q-1}{2}}ln{\bigg \{}{\frac {cosh(\alpha '+\gamma )}{cosh(\ \alpha '-\gamma )}}{\bigg \}}}$. / G ψ = E ψ =  E [G+].... .. 

Se o sistema não apresenta interações, se tem${\displaystyle \gamma =0}$, e ${\displaystyle \alpha '=0}$. Por sua vez, quando as interações são levadas em conta, com ${\displaystyle \gamma \neq 0}$, a expressão para ${\displaystyle \alpha '}$ deverá ser escrita em termos de uma expansão em série de Taylor, o que resulta na relação:

${\displaystyle \alpha '=(q-1)tanh\gamma {\Big \{}\alpha '-sech^{2}\gamma {\frac {\alpha '^{3}}{3}}+...{\Big \}}}$/ G ψ = E ψ =  E [G+].... .. 

que levará após um tratamento matemático elaborado, a relação:

${\displaystyle T<T_{c}={\frac {\frac {2J}{k}}{ln({\frac {q}{q-2}})}}}$/ G ψ = E ψ =  E [G+].... .. 

que determina a temperatura de Curie da rede.^[4]

Pode-se perceber que o parâmetro q está relacionado com a configuração dimensional da rede, onde em uma dimensão (isto é, para ${\displaystyle q=2}$), não ocorre transição de fase, ${\displaystyle T_{c}=0}$. Já para ${\displaystyle q=4}$ (o que corresponde a uma rede quadrada), temos ${\displaystyle {\frac {kT_{c}}{j}}={\frac {2}{ln2}}=2,885...}$ , sendo inferior ao valor previsto por Bragg-Williams ${\displaystyle {\frac {kT_{c}}{J}}=4}$ e ainda acima dos previsto por Onsager ${\displaystyle {\frac {kT_{c}}{J}}={\frac {2}{ln(1+{\sqrt {2}})}}=2,269...}$ .^[3]/ G ψ = E ψ =  E [G+].... .. 

Tomando a afirmação para a aproximação de Bragg-Williams,^[11]  que afirma que “não há nenhuma ordem de curto alcance além daquela estabelecida na ordem de longo alcance”, é possível expressar o parâmetro de ordem de longo alcance ${\displaystyle {\overline {L}}}$ igual ao parâmetro ordem de curto alcance ${\displaystyle {\overline {\sigma }}}$ na rede, logo:^[4]

${\displaystyle {\overline {L}}={\frac {(Z_{+}+Z_{-})-1}{(Z_{+}+Z_{-})+1}}={\frac {sinh(2\alpha +2\alpha ')}{cosh(2\alpha +2\alpha ')+exp(-2\gamma )}}}$./ G ψ = E ψ =  E [G+].... .. 

E no limite ${\displaystyle B\rightarrow 0}$ considerando temperaturas menores mais próximas da temperatura de Curie, chega-se a:

${\displaystyle {\overline {L_{0}}}\simeq {\bigg [}3{\frac {q}{q-1}}{\bigg \{}{\frac {q}{2}}ln{\bigg (}{\frac {q}{q-2}}{\bigg )}{\bigg \}}{\bigg (}1-{\frac {T}{T_{c}}}{\bigg )}{\bigg ]}^{1/2}}$/ G ψ = E ψ =  E [G+].... .. 

Uma correspondência com a aproximação de ordem zero pode ser obtida se for considerado ${\displaystyle q\gg 1}$ na equação acima. Também nota-se que a curva de magnetização espontânea, apresenta comportamento similar ao de ordem zero , ainda que dependendo explicitamente do valor ${\displaystyle q}$.^[4]

Se consideradas as correlações existentes entre os spins vizinhos na rede, pode-se avaliar os números ${\displaystyle N_{++}}$, ${\displaystyle N_{--}}$ e ${\displaystyle N_{+-}}$ (onde cada um dos termos representa o número total de pares de spins conforme indicado nos índices), em função dos parâmetros ${\displaystyle \alpha }$, ${\displaystyle \alpha '}$, e ${\displaystyle \gamma }$. Modificando os somatórios na expressão da função partição canônica do sistema, é possível escrever para ${\displaystyle \sigma _{0}=\sigma _{1}=+1}$ , ${\displaystyle \sigma _{0}=-\sigma _{1}=-1}$  e ${\displaystyle \sigma _{0}=-\sigma _{1}=\pm 1}$ ^[4]:

${\displaystyle Z=Z_{++}+Z_{--}+Z_{+-}}$

Como a função partição está diretamente ligada aos números ${\displaystyle N_{++}}$, ${\displaystyle N_{--}}$ e ${\displaystyle N_{+-}}$, pode-se reescrever, usando a equação de estado da aproximação, considerando que são diretamente proporcionais as funções partição ${\displaystyle Z_{++}}$, ${\displaystyle Z_{--}}$ e ${\displaystyle Z_{+-}}$ com a devida normalização^[4],

${\displaystyle ({\overline {N_{++}}},{\overline {N_{--}}},{\overline {N_{+-}}})={\frac {1}{2}}qN{\frac {(e^{2\alpha +2\alpha '+\gamma },e^{-2\alpha -2\alpha '+\gamma },2e^{-\gamma })}{2{\big \{}e^{\gamma }cosh(2\alpha +2\alpha '){\big \}}}}}$/ G ψ = E ψ =  E [G+].... .. 

onde

${\displaystyle {\frac {{\overline {N_{++}}}{\overline {N_{--}}}}{\overline {N_{}+-}}}={\frac {1}{4}}e^{4\gamma }={\frac {1}{4}}^{4J/kT}}$/ G ψ = E ψ =  E [G+].... .. 

A presença do fator ${\displaystyle exp(4J/kT)}$ para o caso do ferromagnetismo, favorece a formação de pares de spin paralelos ${\displaystyle \uparrow \uparrow }$ e ${\displaystyle \downarrow \downarrow }$, ao invés de pares de spins antiparalelos ${\displaystyle \uparrow \downarrow }$ e ${\displaystyle \downarrow \uparrow }$. Em termos de correlação isso significa dizer: positivo para spins vizinhos ${\displaystyle \uparrow \uparrow }$ ou ${\displaystyle \downarrow \downarrow }$ e negativo para spins vizinhos ${\displaystyle \uparrow }$ e ${\displaystyle \downarrow }$.

Comparações entre as aproximações

A aproximação de ordem zero (chamada também de aproximação de campo médio ou de Bragg-Williams) supõe que a energia de um átomo em um determinado sistema é definida pelo grau de ordem que domina o mesmo e não pelas configurações (flutuações) com os átomos vizinhos.^[4] Já a aproximação de primeira ordem (ou aproximação de Bethe), considera a interação de um spin central e seus vizinhos mais próximos, de modo que o hamiltoniano do grupo é dado pela interação do spin central com seus vizinhos somado à interação desses vizinhos com o resto da rede por um campo molecular médio, como considerado na aproximação zero.^[4]

A diferença na forma de construção das aproximações de ordem zero e primeira ordem, conduz a um resultado para o calor específico e energia interna do sistema, que torna o modelo de primeira ordem mais próximo dos sistemas físicos reais. Pode-se obter facilmente a equação para energia interna de ordem zero a partir da energia média do sistema, quando ${\displaystyle B\rightarrow 0}$,

${\displaystyle C_{o}={\frac {Nk{\overline {L_{o}}}^{2}}{{\frac {{\Biggl (}{\frac {T}{T_{c}}}{\Biggr )}^{2}}{1-{\overline {L_{o}}}^{2}}}-{\frac {T}{T_{c}}}}}}$/ G ψ = E ψ =  E [G+].... .. 

Quando a temperatura se aproxima da temperatura crítica a partir de temperaturas menores que ela, o calor específico acima apresenta uma descontinuidade no ponto de transição, de acordo com a equação abaixo.

${\displaystyle C_{o}={\frac {3}{2}}Nk}$/ G ψ = E ψ =  E [G+].... .. 

e quando ${\displaystyle T\rightarrow 0}$, o calor específico desaparece pois:

${\displaystyle C_{o}(T)\approx 4Nk{\Biggl (}{\frac {T_{c}}{T}}{\Biggr )}^{2}\exp {\Biggl (}{\frac {-2T_{c}}{T}}{\Biggr )}}$/ G ψ = E ψ =  E [G+].... .. 

Pode-se obter a energia interna e o calor específico da aproximação de primeira ordem, através da energia de configuração média na ausência de um campo externo, dada por:

${\displaystyle U_{o}(T)=-J{\Biggl (}{\frac {1}{2}}qN-2q{\overline {N}}_{+}+4{\overline {N}}++{\Biggr )}}$/ G ψ = E ψ =  E [G+].... .. 

A partir da  medida do parâmetro de longo alcance no sistema (${\displaystyle {\overline {L}}}$), obtemos os ${\displaystyle {\overline {N}}_{+}}$ e ${\displaystyle {\overline {N}}++}$, que aplicadas na expressão de energia e desenvolvidas conduzem ao resultado:

${\displaystyle U_{o}=-{\frac {1}{2}}qJN{\frac {\cosh(2\alpha ^{\prime })-\exp(-2\gamma )}{\cosh(2\alpha ^{\prime })+\exp(-2\gamma )}}.}$/ G ψ = E ψ =  E [G+].... .. 

Para temperatura ${\displaystyle T>T_{c}}$ e ${\displaystyle \alpha \prime =0}$, a energia interna e o calor específico serão dados como segue:

${\displaystyle U_{o}=-{\frac {1}{2}}qJN\tanh(\gamma ),}$/ G ψ = E ψ =  E [G+].... .. 

${\displaystyle {\frac {C_{o}}{Nk}}={\frac {1}{2}}q\gamma ^{2}\operatorname {sech} ^{2}(\gamma ).}$/ G ψ = E ψ =  E [G+].... .. 

Na equação do calor específico para primeira ordem e ordem zero, quando  ${\displaystyle T\rightarrow \infty }$, o calor específico tende a zero. Entretanto para ${\displaystyle T>T_{c}}$, o calor específico obtido pela aproximação de Bethe conduz a um valor não nulo, devido a interação de curto alcance do spin central com seus spins vizinhos.^[3][10] Em contraste, a aproximação de Bragg-Williams, apresenta uma descontinuidade no calor específico em ${\displaystyle T=T_{c}}$, se anulando em valores acima da temperatura crítica. Para todo q a primeira aproximação produz uma temperatura de transição mais próxima do valor correto da temperatura crítica, do que a aproximação de ordem zero.

Assim como o calor específico, é possível analisar também a energia interna em ambas aproximações. A energia interna é contínua na aproximação de Bethe, acima da temperatura crítica, como era de se esperar, dado que o calor específico se mantém nessa aproximação, devido a interação de curto alcance entre os spins. Já aproximação de Bragg-Williams, a energia interna é nula, para toda temperatura acima da temperatura crítica, resultado esse consistente com o gráfico de calor específico.