Modelo de Potts




Em mecânica estatística, o modelo de Potts, uma generalização do modelo Ising, é um modelo de spins em interação em um reticulado cristalino. Pelo estudo do modelo de Potts, pode-se ter uma visão do comportamento de ferromagnetos e outros fenômenos em física do estado sólido. A força do modelo de Potts está menos no fato de que modela bem estes sistemas físicos e mais no fato de que o caso unidimensional é exatamente solvável. O modelo de Potts tem uma rica formulação matemática, que tem sido extensivamente estudada.


O modelo recebe este nome em homenagem ao matemático australiano Renfrey Potts, que descreveu o modelo perto da conclusão de sua tese de doutorado em 1951. O modelo estava relacionado como o "modelo de Potts planar" ou "modelo do relógio", sugerido a Potts por seu orientador, o físico britânico Cyril Domb. O modelo de Potts planar de quatro estados é às vezes chamado de modelo Ashkin-Teller, em homenagem aos físicos Julius Ashkin e Edward Teller, que falaram sobre um modelo equivalente em 1943.[1]


O modelo de Potts está relacionado a e é generalizado por vários outros modelos, incluindo o modelo XY, o modelo de Heisenberg e o modelo n-vetor. O modelo de Potts de intervalo infinito é conhecido como modelo de Kac. Quando se assume que os spins interagem de maneira não-abeliana, o modelo está relacionado com o modelo do tubo de fluxo, usado para discutir confinamento em cromodinâmica quântica. Generalizações do modelo de Potts também têm sido usadas para modelar crescimento de grãos em metais e endurecimento de espumas. Uma generalização adicional destes métodos por James Glazier e François Graner, conhecida como modelo de Potts celular, tem sido usada para simular fenômenos estáticos e cinéticos em espuma e morfogênese biológica.[2]




Índice






  • 1 Descrição física


  • 2 Discussão


  • 3 Descrição teórica da medida


    • 3.1 Topologia do espaço de estados


    • 3.2 Energia de interação


    • 3.3 Função de partição e medida


    • 3.4 Solução do campo livre


    • 3.5 Modelo interativo




  • 4 O modelo de Potts no processamento de sinal e imagem


  • 5 Ver também


  • 6 Referências





Descrição física |


O modelo de Potts consiste em spins colocados em um reticulado. O reticulado é geralmente assumido como um reticulado euclidiano retangular bidimensional, mas é frequentemente generalizado a outras dimensões e outros reticulados. Domb sugeriu originalmente que o spin assume um de q{displaystyle q}q valores possíveis, distribuídos uniformemente pelo círculo, em ângulos


θn=2πnq,{displaystyle theta _{n}={frac {2pi n}{q}},}{displaystyle theta _{n}={frac {2pi n}{q}},}

em que n=1,...,q{displaystyle n=1,...,q}{displaystyle n=1,...,q} e o hamiltoniano da interação é dado por


Hc=Jc∑(i,j)cos⁡si−θsj){displaystyle H_{c}=J_{c}sum _{(i,j)}cos left(theta _{s_{i}}-theta _{s_{j}}right)}{displaystyle H_{c}=J_{c}sum _{(i,j)}cos left(theta _{s_{i}}-theta _{s_{j}}right)}

com a soma correndo pelos pares de vizinhos mais próximos (i,j){displaystyle (i,j)}(i,j) sobre todos os locais do reticulado. As cores do local si{displaystyle s_{i}}{displaystyle s_{i}} assumem valores em {1,...,q}{displaystyle {1,...,q}}{displaystyle {1,...,q}}. Aqui, Jc{displaystyle J_{c}}J_{c} é uma constante de acoplamento, que determina a força da interação. Este modelo é agora conhecido como modelo de Potts vetorial ou modelo do relógio. Potts forneceu a locação em duas dimensões da transição de fase para q=3{displaystyle q=3}{displaystyle q=3} e q=4{displaystyle q=4}{displaystyle q=4}.[3] No limite conforme q→{displaystyle qrightarrow infty }{displaystyle qrightarrow infty }, este se torna o modelo XY.


O que é agora conhecido como o modelo de Potts padrão foi sugerido por Potts na disciplina sobre seu estudo acima e usa um hamiltoniano mais simples, dado por:


Hp=−Jp∑(i,j)δ(si,sj),{displaystyle H_{p}=-J_{p}sum _{(i,j)}delta (s_{i},s_{j}),}{displaystyle H_{p}=-J_{p}sum _{(i,j)}delta (s_{i},s_{j}),}

em que δ(si,sj){displaystyle delta (s_{i},s_{j})}{displaystyle delta (s_{i},s_{j})} é o delta de Kronecker, que é igual a 1 sempre que si=sj{displaystyle s_{i}=s_{j}}{displaystyle s_{i}=s_{j}} e 0 de outro modo.


O modelo de Potts padrão q=2{displaystyle q=2}{displaystyle q=2} é equivalente ao modelo Ising e ao modelo de Potts vetorial de dois estados com Jp=−2Jc{displaystyle J_{p}=-2J_{c}}{displaystyle J_{p}=-2J_{c}}. O modelo de Potts padrão q=3{displaystyle q=3}{displaystyle q=3} é equivalente ao modelo de Potts vetorial de três estados com Jp=−(3/2)Jc{displaystyle J_{p}=-(3/2)J_{c}}{displaystyle J_{p}=-(3/2)J_{c}}.


Uma generalização comum consiste em introduzir um termo de "campo magnético" externo h{displaystyle h}h, movendo os parâmetros no interior das somas e permitindo que variem através do modelo


βHg=−β(i,j)Jijδ(si,sj)−ihisi,{displaystyle beta H_{g}=-beta sum _{(i,j)}J_{ij}delta (s_{i},s_{j})-sum _{i}h_{i}s_{i},}{displaystyle beta H_{g}=-beta sum _{(i,j)}J_{ij}delta (s_{i},s_{j})-sum _{i}h_{i}s_{i},}

em que β=1/kT{displaystyle beta =1/kT}{displaystyle beta =1/kT} é a temperatura inversa, k{displaystyle k}k é a constante de Boltzmann e T{displaystyle T}T é a temperatura. A soma pode correr por vizinhos mais distantes no reticulado ou pode, na verdade, ser uma força de intervalo infinito.


Textos diferentes podem adotar convenções ligeiramente diferentes, o que pode alterar H{displaystyle H}H e a função de partição associada pelas constantes aditivas ou multiplicativas.



Discussão |


Apesar de sua simplicidade como um modelo de sistema físico, o modelo de Potts é útil como um sistema de modelo para o estudo de transições de fase. Por exemplo, reticulados bidimensionais com J>0{displaystyle J>0}{displaystyle J>0} exibem uma transição de primeira ordem se q>4{displaystyle q>4}{displaystyle q>4}. Quando q≤4{displaystyle qleq 4}{displaystyle qleq 4}, uma transição contínua é observada, assim como no modelo Ising em que q=2{displaystyle q=2}{displaystyle q=2}. Outro uso é encontrado na relação do modelo com problemas de percolação, além polinômios de Tutte e polinômios cromáticos encontrados em combinatória.[4]


O modelo tem uma relação íntima com o modelo de grupo aleatório de Fortuin–Kasteleyn, outro modelo em mecânica estatística. A compreensão desta relação tem ajudado a desenvolver eficientes métodos de Monte Carlo de cadeia de Markov para exploração numérica do modelo em q{displaystyle q}q pequeno.


Para valores inteiros de q{displaystyle q}q, q≥3{displaystyle qgeq 3}{displaystyle qgeq 3}, o modelo exibe o fenômeno da adsorção interfacial com propriedades de molhabilidade críticas intrigantes quando se fixam fronteiras opostas em dois estados diferentes.[5]



Descrição teórica da medida |


O modelo de Potts unidimensional pode ser expresso em termos de um subdeslocamento de tipo finito, ganhando assim acesso a todas as técnicas matemáticas associadas com este formalismo. Em particular, pode ser resolvido exatamente pelo uso de técnicas de operadores de transferência. Entretanto, Ernst Ising usava métodos combinatórios para resolver o modelo Ising, o "ancestral" do modelo de Potts, em sua tese de doutorado de 1924. Esta seção desenvolve o formalismo matemático, baseado na teoria da medida, por trás desta solução.


Ainda que o exemplo abaixo seja desenvolvido para o caso unidimensional, muitos dos argumentos e quase toda a notação podem ser generalizados facilmente a qualquer número de dimensões. Algo do formalismo também é amplo o bastante para lidar com modelos relacionados, como o modelo XY, o modelo de Heisenberg e o modelo n-vetor.



Topologia do espaço de estados |


Considere Q={1,...,q}{displaystyle Q={1,...,q}}{displaystyle Q={1,...,q}} um conjunto finito de símbolos e considere


QZ={s=(…,s−1,s0,s1,…):sk∈Q∀k∈Z}{displaystyle Q^{mathbf {Z} }={s=(ldots ,s_{-1},s_{0},s_{1},ldots ):s_{k}in Q;forall kin mathbf {Z} }}{displaystyle Q^{mathbf {Z} }={s=(ldots ,s_{-1},s_{0},s_{1},ldots ):s_{k}in Q;forall kin mathbf {Z} }}

o conjunto de todas as cordas de valores bi-infinitas a partir do conjunto Q{displaystyle Q}Q. Este conjunto é chamado de deslocamento total. Para a definição do modelo de Potts, tanto este espaço inteiro, como um certo subconjunto do mesmo, um subdeslocamento de tipo finito, podem ser usados. Os deslocamentos recebem este nome porque há um operador natural neste espaço, o operador de deslocamento τ:QZ→QZ{displaystyle tau :Q^{mathbf {Z} }rightarrow Q^{mathbf {Z} }}{displaystyle tau :Q^{mathbf {Z} }rightarrow Q^{mathbf {Z} }}, agindo como


τ(s)k=sk+1.{displaystyle tau (s)_{k}=s_{k+1}.}{displaystyle tau (s)_{k}=s_{k+1}.}

Este conjunto tem uma topologia produto natural. A base para este topologia é composta pelos conjuntos cilindros


Cm[ξ0,…k]={s∈QZ:sm=ξ0,…,sm+k=ξk},{displaystyle C_{m}[xi _{0},ldots ,xi _{k}]={sin Q^{mathbf {Z} }:s_{m}=xi _{0},ldots ,s_{m+k}=xi _{k}},}{displaystyle C_{m}[xi _{0},ldots ,xi _{k}]={sin Q^{mathbf {Z} }:s_{m}=xi _{0},ldots ,s_{m+k}=xi _{k}},}

isto é, o conjunto de todas as cordas possíveis em que k+1{displaystyle k+1}k+1 spins correspondem exatamente a um dado conjunto específico de valores ξ0,...,ξk{displaystyle xi _{0},...,xi _{k}}{displaystyle xi _{0},...,xi _{k}}. Representações explícitas para os conjuntos cilindros podem ser apreendidas ao notar que a corda de valores corresponde ao número q-ádico e assim, intuitivamente, a topologia produto lembra aquela da linha de números reais.[6]



Energia de interação |


A interação entre os spins é então dada por uma função contínua V:QZ→R{displaystyle V:Q^{mathbf {Z} }rightarrow mathbf {R} }{displaystyle V:Q^{mathbf {Z} }rightarrow mathbf {R} } nesta topologia. Qualquer função fará. Por exemplo,


V(s)=−(s0,s1){displaystyle V(s)=-Jdelta (s_{0},s_{1})}{displaystyle V(s)=-Jdelta (s_{0},s_{1})}

será visto para descrever a interação entre os vizinhos mais próximos. É claro que funções diferentes dão interações diferentes. Então, uma função de s0{displaystyle s_{0}}s_{0}, s1{displaystyle s_{1}}{displaystyle s_{1}} e s2{displaystyle s_{2}}{displaystyle s_{2}} descreverá uma interação entre os vizinhos mais próximos seguintes. Uma função V{displaystyle V}V dá a energia de interação entre um conjunto de spins, que não é o hamiltoniano, mas é usada em sua construção. O argumento para a função V{displaystyle V}V é um elemento s∈QZ{displaystyle sin Q^{mathbf {Z} }}{displaystyle sin Q^{mathbf {Z} }}, isto é, uma corda infinita de spins. No exemplo acima, a função V{displaystyle V}V apenas selecionou dois spins da corda infinita: os valores s0{displaystyle s_{0}}s_{0} e s1{displaystyle s_{1}}{displaystyle s_{1}}. Em geral, a função V{displaystyle V}V pode depender de alguns ou todos os spins. Atualmente, apenas aqueles que dependem de um número finito são exatamente solváveis.


Defina a função Hn:QZ→R{displaystyle H_{n}:Q^{mathbf {Z} }rightarrow mathbf {R} }{displaystyle H_{n}:Q^{mathbf {Z} }rightarrow mathbf {R} } como


Hn(s)=∑k=0nV(τks).{displaystyle H_{n}(s)=sum _{k=0}^{n}V(tau ^{k}s).}{displaystyle H_{n}(s)=sum _{k=0}^{n}V(tau ^{k}s).}

Pode-se ver que esta função consiste em duas partes: a auto-energia de uma configuração [s0,s1,...,sn]{displaystyle [s_{0},s_{1},...,s_{n}]}{displaystyle [s_{0},s_{1},...,s_{n}]} de spins, além da energia de interação deste conjunto e todos os outros spins no reticulado. O limite n→{displaystyle nrightarrow infty }{displaystyle nrightarrow infty } desta função é o hamiltoniano do sistema. Para n{displaystyle n}n finito, estes são às vezes chamados de hamiltonianos de estado finito.[6]



Função de partição e medida |


A função de partição de estado finito correspondente é dada por


Zn(V)=∑s0,…,sn∈Qexp⁡(−βHn(C0[s0,s1,…,sn])){displaystyle Z_{n}(V)=sum _{s_{0},ldots ,s_{n}in Q}exp(-beta H_{n}(C_{0}[s_{0},s_{1},ldots ,s_{n}]))}{displaystyle Z_{n}(V)=sum _{s_{0},ldots ,s_{n}in Q}exp(-beta H_{n}(C_{0}[s_{0},s_{1},ldots ,s_{n}]))}

com C0{displaystyle C_{0}}{displaystyle C_{0}} sendo os conjuntos cilindros definidos acima. Aqui, β=1/kT{displaystyle beta =1/kT}{displaystyle beta =1/kT}, em que k{displaystyle k}k é a constante de Boltzmann e T{displaystyle T}T é a temperatura. É muito comum em tratamentos matemáticas configurar β=1{displaystyle beta =1}{displaystyle beta =1}, já que é facilmente recuperado ao reescalar a energia de interação. A função de partição é escrita como uma função da interação V{displaystyle V}V para enfatizar que isto é apenas uma função da interação e não de qualquer configuração específica de spins. A função de partição e o hamiltoniano são usados para definir uma medida sobre a sigma-álgebra de Borel da seguinte maneira: a medida de um conjunto cilindro, isto é, um elemento da base, é dada por


μ(Ck[s0,s1,…,sn])=1Zn(V)exp⁡(−βHn(Ck[s0,s1,…,sn])).{displaystyle mu (C_{k}[s_{0},s_{1},ldots ,s_{n}])={frac {1}{Z_{n}(V)}}exp(-beta H_{n}(C_{k}[s_{0},s_{1},ldots ,s_{n}])).}{displaystyle mu (C_{k}[s_{0},s_{1},ldots ,s_{n}])={frac {1}{Z_{n}(V)}}exp(-beta H_{n}(C_{k}[s_{0},s_{1},ldots ,s_{n}])).}

Pode-se então estender por aditividade contável à sigma-álgebra total. Esta medida é uma medida de probabilidade, que dá a probabilidade de que uma dada configuração ocorra no espaço de configuração QZ{displaystyle Q^{mathbf {Z} }}{displaystyle Q^{mathbf {Z} }}. Ao dotar o espaço de configuração com uma medida de probabilidade construída a partir de um hamiltoniano desta forma, o espaço de configuração se transforma em um conjunto canônico.


A maioria das propriedades termodinâmicas pode ser expressa diretamente em termos de uma função de partição. Assim, por exemplo, a energia livre de Helmholtz é dada por


An(V)=−kTlog⁡Zn(V).{displaystyle A_{n}(V)=-kTlog Z_{n}(V).}{displaystyle A_{n}(V)=-kTlog Z_{n}(V).}

Outra importante quantidade relacionada é a pressão topológica, definida como


P(V)=limn→1nlog⁡Zn(V),{displaystyle P(V)=lim _{nto infty }{frac {1}{n}}log Z_{n}(V),}{displaystyle P(V)=lim _{nto infty }{frac {1}{n}}log Z_{n}(V),}

que aparecerá como o logaritmo do autovalor condutor do operador de transferência da solução.[6]



Solução do campo livre |


O modelo mais simples é o modelo em que não há qualquer interação, então, V=c{displaystyle V=c}{displaystyle V=c} e Hn=c{displaystyle H_{n}=c}{displaystyle H_{n}=c} (com c{displaystyle c}c constante e independente de qualquer configuração de spin). A função de partição se torna


Zn(c)=e−s0,…,sn∈Q1.{displaystyle Z_{n}(c)=e^{-cbeta }sum _{s_{0},ldots ,s_{n}in Q}1.}{displaystyle Z_{n}(c)=e^{-cbeta }sum _{s_{0},ldots ,s_{n}in Q}1.}

Se todos os estados forem permitidos, isto é, o conjunto subjacente de estados for dado por um deslocamento total, então, a soma pode ser trivialmente avaliada como


Zn(c)=e−qn+1.{displaystyle Z_{n}(c)=e^{-cbeta }q^{n+1}.}{displaystyle Z_{n}(c)=e^{-cbeta }q^{n+1}.}

Se spins avizinhantes forem permitidos apenas em certas configurações específicas, então, o espaço de estados é dado por um subdeslocamento de tipo finito. A função de partição pode então ser escrita como


Zn(c)=e−|Fixτn|=e−TrAn,{displaystyle Z_{n}(c)=e^{-cbeta }|{mbox{Fix}},tau ^{n}|=e^{-cbeta }{mbox{Tr}}A^{n},}{displaystyle Z_{n}(c)=e^{-cbeta }|{mbox{Fix}},tau ^{n}|=e^{-cbeta }{mbox{Tr}}A^{n},}

em que c{displaystyle c}c é a cardinalidade ou contagem de um conjunto e Fix{displaystyle {mbox{Fix}}}{displaystyle {mbox{Fix}}} é o conjunto de pontos fixos da função iterada de deslocamento:


Fixτn={s∈QZ:τns=s}.{displaystyle {mbox{Fix}},tau ^{n}={sin Q^{mathbf {Z} }:tau ^{n}s=s}.}{displaystyle {mbox{Fix}},tau ^{n}={sin Q^{mathbf {Z} }:tau ^{n}s=s}.}

A matriz q x q A{displaystyle A}A é a matriz de adjacência que especifica quais valores de spins avizinhantes são permitidos.[6]



Modelo interativo |


O caso mais simples de modelo interativo é o modelo Ising, em que o spin pode assumir apenas um de dois valores, sn∈{−1,1}{displaystyle s_{n}in {-1,1}}{displaystyle s_{n}in {-1,1}}, e apenas spins dos vizinhos mais próximos interagem. O potencial de interação é dado por


V(σ)=−Jps0s1.{displaystyle V(sigma )=-J_{p}s_{0}s_{1}.}{displaystyle V(sigma )=-J_{p}s_{0}s_{1}.}

Este potencial pode ser capturado em uma matriz 2 x 2 com elementos de matriz


σ′=exp⁡Jpσσ′){displaystyle M_{sigma sigma '}=exp left(beta J_{p}sigma sigma 'right)}{displaystyle M_{sigma sigma '}=exp left(beta J_{p}sigma sigma 'right)}

com o índice σ′∈{−1,1}{displaystyle sigma ,sigma 'in {-1,1}}{displaystyle sigma ,sigma 'in {-1,1}}. A função de partição é então dada por


Zn(V)=TrMn.{displaystyle Z_{n}(V)={mbox{Tr}},M^{n}.}{displaystyle Z_{n}(V)={mbox{Tr}},M^{n}.}

A solução geral para um número arbitrário de spins e uma interação arbitrária de intervalo finito é dada pela mesma forma geral. Neste caso, a expressão precisa para a matriz M{displaystyle M}M é um pouco mais complexa.


O objetivo ao resolver um modelo tal como o modelo de Potts é dar uma expressão de forma fechada para a função de partição e uma expressão para os estados de Gibbs ou estados de equilíbrio no limite de n→{displaystyle nrightarrow infty }{displaystyle nrightarrow infty }, o limite termodinâmico.[6]



O modelo de Potts no processamento de sinal e imagem |


O modelo de Potts tem aplicações em reconstrução de sinal. Assuma que há uma observação com ruído de um sinal constante por partes g{displaystyle g}g em Rn{displaystyle mathbf {R} ^{n}}{mathbf  {R}}^{n}. Para recuperar g{displaystyle g}g a partir do vetor f{displaystyle f}f da observação com ruído em Rn{displaystyle mathbf {R} ^{n}}{mathbf  {R}}^{n}, procura-se um minimizador do problema inverso correspondente, a funcional de Lp-Potts (u){displaystyle P_{gamma }(u)}{displaystyle P_{gamma }(u)}, que é definida por


(u)=γu‖0+‖u−f‖pp=γ#{i:ui≠ui+1}+∑i=1n|ui−fi|p.{displaystyle P_{gamma }(u)=gamma |nabla u|_{0}+|u-f|_{p}^{p}=gamma #{i:u_{i}neq u_{i+1}}+sum _{i=1}^{n}|u_{i}-f_{i}|^{p}.}{displaystyle P_{gamma }(u)=gamma |nabla u|_{0}+|u-f|_{p}^{p}=gamma #{i:u_{i}neq u_{i+1}}+sum _{i=1}^{n}|u_{i}-f_{i}|^{p}.}

A penalidade do salto u‖0{displaystyle |nabla u|_{0}}{displaystyle |nabla u|_{0}} força soluções constantes por partes e o termo de dado u−f‖pp{displaystyle |u-f|_{p}^{p}}{displaystyle |u-f|_{p}^{p}} acopla o candidato minimizador u{displaystyle u}u ao dado f{displaystyle f}f. O parâmetro γ>0{displaystyle gamma >0}{displaystyle gamma >0} controla o tradeoff entre regularidade e fidelidade dos dados. Há rápidos algoritmos para a exata minimização da funcional de L1 e L2-Potts.[7]


Em processamento de imagem, a funcional de Potts está relacionada com o problema de segmentação. Entretanto, em duas dimensões, o problema é NP-difícil.[8]



Ver também |


  • Modelo Ising


Referências |





  1. Ashkin, J. (1943). «Statistics of Two-Dimensional Lattices with Four Components». Physical Review. 64 (5-6): 178–184. doi:10.1103/physrev.64.178 


  2. Graner, François (1992). «Simulation of biological cell sorting using a two-dimensional extended Potts model». Physical Review Letters. 69 (13): 2013–2016. doi:10.1103/physrevlett.69.2013 


  3. Potts, R. B. (janeiro de 1952). «Some generalized order-disorder transformations». Mathematical Proceedings of the Cambridge Philosophical Society. 48 (1): 106–109. ISSN 1469-8064. doi:10.1017/s0305004100027419 


  4. «Code for Computing Tutte Polynomials». homepages.mcs.vuw.ac.nz (em inglês). Consultado em 14 de julho de 2017 


  5. Selke, Walter; Huse, David A. (1 de junho de 1983). «Interfacial adsorption in planar potts models». Zeitschrift für Physik B Condensed Matter (em inglês). 50 (2): 113–116. ISSN 0722-3277. doi:10.1007/bf01304093 


  6. abcde Wu, F. Y. (1982). «The Potts model». Reviews of Modern Physics. 54 (1): 235–268. doi:10.1103/revmodphys.54.235 


  7. Friedrich, F.; Kempe, A.; Liebscher, V.; Winkler, G. (2008). «Complexity Penalized M-Estimation: Fast Computation». Journal of Computational and Graphical Statistics. 17 (1): 201–224 


  8. Boykov, Y.; Veksler, O.; Zabih, R. (novembro de 2001). «Fast approximate energy minimization via graph cuts». IEEE Transactions on Pattern Analysis and Machine Intelligence. 23 (11): 1222–1239. ISSN 0162-8828. doi:10.1109/34.969114 


































































Popular posts from this blog

404 Error Contact Form 7 ajax form submitting

How to know if a Active Directory user can login interactively

TypeError: fit_transform() missing 1 required positional argument: 'X'