Modelo de Ising 2D
Simulação Estocástica e Transições de Fase

01. Fundamentação Teórica

O Modelo de Ising Bidimensional é um dos paradigmas mais fundamentais da Mecânica Estatística moderna. Ele foi formulado para descrever a natureza de materiais ferromagnéticos em função da temperatura. Em sua essência, o sistema consiste em uma rede quadrada bidimensional de dimensão $n \times n$.

Cada sítio da malha abriga uma variável aleatória chamada spin ($\sigma_i \in \{-1, +1\}$), que pode estar orientada "para cima" ou "para baixo". A energia total de uma configuração do sistema é ditada pela sua função Hamiltoniana:

Onde $J$ dita a constante de acoplamento (se positiva, os spins vizinhos tendem a se alinhar, caracterizando ferromagnetismo) e $h$ é o campo magnético externo. Ao cruzar a temperatura crítica (Ponto de Curie do sistema), a simulação revela o desabrochar espontâneo da magnetização — a clássica transição de fase de segunda ordem.

02. Dinâmica de Metropolis e Performance JIT

Dado o crescimento exponencial dos microestados $2^N$, calcular a Função de Partição de forma analítica exata é inviável para redes amplas. Foi adotado o rigor estocástico dos Métodos de Monte Carlo. O coração da atualização computacional baseia-se no Algoritmo de Metropolis-Hastings, que assegura o alcance da distribuição de equilíbrio de Gibbs por meio de cadeias de Markov regulares.

Diferenciais Técnicos Alcançados:

• Aceleração JIT extrema: Uso cirúrgico do framework `numba` para compilar as funções críticas em código de máquina via LLVM, eliminando o overhead do interpretador Python (CPython) nos laços de atualização de spins.
• Multilinguagem: Estudo arquitetado e testado usando Python para a interatividade e a linguagem R como pilar de validação analítica da amostragem base.
• UX Científico Moderno: Criação de front-end em Streamlit puramente via código, fornecendo a estudantes e curiosos um laboratório virtual em tempo real (painel lateral de interatividade).

03. Análise Físico-Estatística Avançada

O relatório técnico incorporado apresenta o aprofundamento rigoroso das propriedades do modelo e da engenharia de simulação implementada. Dentre as descobertas e formulações fundamentais descritas, destacam-se:

• Balanço Detalhado e Ergodicidade: Demonstração matemática de como as probabilidades de aceitação no critério de Metropolis satisfazem a condição de balanço detalhado, garantindo a convergência assintótica do microestado da rede para a distribuição de equilíbrio de Gibbs.
• Grandezas Observáveis Termodinâmicas: Dedução física das flutuações microscópicas de energia e magnetização para estimativa de variáveis macroscópicas, como a Magnetização Média, Calor Específico ($C$) e Suscetibilidade Magnética ($\chi$).
• Desaceleramento Crítico (Critical Slowing Down): Análise da divergência do tempo de autocorrelação integrado ($\tau_{\text{int}} \sim |T - T_c|^{-z\nu}$) próximo ao ponto de transição de fase ($T_c \approx 2,269 \,J/k_B$), ilustrando as limitações computacionais de algoritmos de atualização local nas proximidades de fenômenos críticos.
• Invariância de Escala e Dimensão Fractal: Estudo geométrico da fronteira dos domínios magnéticos na temperatura crítica, evidenciando visualmente a emergência de comportamento fractal, cuja dimensão teórica exata é conhecida como $D_f = 1,375$ (Evolução de Schramm-Loewner).

03. Simulação Interativa (Live App)

Sinta a termodinâmica rodando na nuvem. Altere o parâmetro $\beta$ (inverso da Temperatura) e veja a rede de spins ferromagnéticos quebrar sua simetria visualmente (em baixas temperaturas formam-se domínios magnéticos densos; em altas temperaturas prevalece o caos e a desordem térmica).