Visualisation du Modèle d'Ising

Transitions de Phase et Phénomènes Critiques en Mécanique Statistique - Simulation Monte Carlo Interactive

Wilhelm Lenz (1920) · Ernst Ising (1925) · Lars Onsager (1944)

Température (T) 2.27
T/Tc 1.00
Énergie E -1.50
Aimantation |M| 0.85
Pas Monte Carlo 0

Contrôles de Simulation

Basse 2.27 Haute
Température Critique Tc ≈ 2.269
-2.0 0.0 +2.0
Antiferro -1.0 1.0 Ferro +1.0
20 50 100
1 20 100

Évolution de l'Énergie

Évolution de l'Aimantation

Théorie du Modèle d'Ising

Le modèle d'Ising est l'un des modèles les plus emblématiques en mécanique statistique, décrivant le comportement d'interaction des spins sur un réseau.

Hamiltonien

H = -J Σ<ij> sᵢsⱼ - h Σᵢ sᵢ
  • sᵢ = ±1 - Direction du spin (haut/bas)
  • J - Constante de couplage (J>0 ferromagnétique, J<0 antiferromagnétique)
  • h - Intensité du champ magnétique externe
  • Σ<ij> - Somme sur les spins les plus proches voisins

Jalons

  • 1920 - Wilhelm Lenz propose le modèle
  • 1925 - Ernst Ising résout le cas 1D (pas de transition de phase)
  • 1944 - Lars Onsager résout exactement le cas 2D (découvre la transition de phase)

Phénomènes de Transition de Phase

Température Critique

Tc = 2/ln(1+√2) ≈ 2.269

Près de la température critique, le système subit une transition de l'état ordonné à désordonné.

Trois Régimes

Basse Température T < Tc

Phase ferromagnétique ordonnée. Rupture spontanée de symétrie, la plupart des spins pointent dans la même direction, aimantation |M| > 0.

Point Critique T ≈ Tc

Fluctuations critiques. Apparition d'amas à grande échelle, ralentissement critique, susceptibilité magnétique divergente.

Haute Température T > Tc

Phase paramagnétique désordonnée. Spins orientés aléatoirement, aimantation moyenne M = 0.

Algorithme de Metropolis-Hastings

Utilisant les méthodes Monte Carlo pour simuler le comportement thermodynamique du système.

Étapes de l'Algorithme

  1. Sélectionner aléatoirement un spin sᵢ
  2. Calculer le changement d'énergie ΔE si retourné
  3. Si ΔE ≤ 0, accepter le retournement
  4. Si ΔE > 0, accepter avec probabilité exp(-ΔE/kT)
  5. Répéter N×N fois pour un pas Monte Carlo

Probabilité d'Acceptation

P(accept) = min(1, e^(-ΔE/kT))
Note : Près du point critique, le système présente le phénomène de 'ralentissement critique', avec une convergence significativement plus lente. L'algorithme de Wolff (retournement d'amas) peut être utilisé pour accélérer.

Guide d'Observation

Ferromagnétisme Basse T (T < 2.0)

Régler T ≈ 1.5, observer de grandes régions de même couleur. C'est l'état ordonné ferromagnétique avec rupture spontanée de symétrie.

Fluctuations Critiques (T ≈ 2.27)

Régler T = 2.27, observer la formation et la mort d'amas à grande échelle. C'est la région la plus intéressante!

Paramagnétisme Haute T (T > 3.0)

Régler T ≈ 4.0, observer les retournements aléatoires de spins. C'est l'état paramagnétique désordonné.

Effet de Champ Externe

Ajuster le champ externe h, observer le biais dans la direction des spins. h > 0 favorise le haut, h < 0 favorise le bas.

Phase Antiferromagnétique (J < 0)

Régler J = -1.0, un état ordonné antiferromagnétique en forme de rayures se forme à basse température.

Compréhension Interactive des Formules

H

Énergie Totale

L'Hamiltonien du système, représentant l'énergie totale. Le système tend vers l'état d'énergie la plus basse.

-J Σ sᵢsⱼ

Terme d'Interaction

Énergie d'interaction des spins les plus proches voisins. J > 0 : même direction a une énergie plus basse (ferromagnétique); J < 0 : direction opposée a une énergie plus basse (antiferromagnétique).

-h Σ sᵢ

Terme de Champ Externe

Énergie du champ externe agissant sur les spins. h > 0 : les spins vers le haut ont une énergie plus basse; h < 0 : les spins vers le bas ont une énergie plus basse.