Fiziksel Model ve Hamiltoniyen: 3B dış manyetik alanlar ve kiral etkileşimleri içeren sürekli spin vektörleri ($|\mathbf{S}_i| = 1$) için Hamiltoniyen şu şekilde tanımlanır:
$$ \mathcal{H} = -J \sum_{\langle i,j \rangle} \mathbf{S}_i \cdot \mathbf{S}_j - \sum_i \mathbf{B} \cdot \mathbf{S}_i - \sum_{\langle i,j \rangle} \mathbf{D}_{ij} \cdot (\mathbf{S}_i \times \mathbf{S}_j) $$
Dzyaloshinskii-Moriya Etkileşimi (DMI): Antisimetrik değiş-tokuş etkileşimi topolojik dokuları kararlı kılar. $\mathbf{D}_{ij}$ vektörünün yönelimi yapısal tersinme asimetrisine (inversion asymmetry) bağlıdır:
Topolojik Yük ($Q$): Topolojik invaryant $Q$, spin konfigürasyonunun birim küreyi kaç kez sardığını sayar. Sürekli limitte $Q = \frac{1}{4\pi}\int \mathbf{S}\cdot(\partial_x \mathbf{S} \times \partial_y \mathbf{S})\,dx\,dy$ ile verilir. Bu simülasyon ayrık örgüde Berg–Lüscher tanımını kullanır: her kare plaket iki üçgene bölünür ve her üçgenin üç spininin birim küre üzerinde gerdiği katı açı hesaplanır:
$$ \Omega_{\triangle} = 2\arctan\frac{\mathbf{S}_1\cdot(\mathbf{S}_2\times\mathbf{S}_3)}{1+\mathbf{S}_1\cdot\mathbf{S}_2+\mathbf{S}_2\cdot\mathbf{S}_3+\mathbf{S}_3\cdot\mathbf{S}_1}, \qquad Q = \frac{1}{4\pi}\sum_{\triangle}\Omega_{\triangle} $$
Periyodik sınır koşulları altında örgü bir tordur ve bu tanım, sıcaklıktan ve gradyanların keskinliğinden bağımsız olarak her konfigürasyon için tam sayı verir. Bu yüzden grafikteki $Q$ basamaklar halinde değişir. Her basamak, bir skyrmionun örgü üzerinde doğduğu ya da yok olduğu bir olaydır. Sürekli kuramda topolojik olarak yasak olan bu süreç, örgü ölçeğinde (komşu spinlerin neredeyse zıt yöneldiği tek bir hücrede) mümkündür. Basit üçlü çarpım yaklaşımı $\frac{1}{4\pi}\sum \mathbf{S}_i\cdot(\mathbf{S}_{i+\hat x}\times\mathbf{S}_{i+\hat y})$ ise yalnızca yavaş değişen dokular için doğrudur; birkaç örgü sabiti büyüklüğündeki skyrmionların sayısını belirgin biçimde eksik verir.
İşaret ve Uzunluk Ölçeği: $B_z > 0$ iken arka plan $+z$ yönündedir ve çekirdeği $-z$ olan her skyrmion $Q = -1$ katkı verir; yani $|Q|$ skyrmion sayısıdır. Karakteristik sarmal (helix) periyodu küçük $D$ için $\lambda \approx 2\pi J/|D|$, kare örgüde ise tam olarak $\tan(2\pi a/\lambda) = |D|/J$ ile verilir ($D = J$ için $\lambda = 8a$). $\lambda$ sistem boyutuna ($L=30$) uymadığında periyodik sınırlar dokuyu sıkıştırır veya gerer. Bu bozulma skyrmionların şeklini ve sayısını etkiler, ancak $Q$'yu kesirli yapmaz.
Monte Carlo Algoritması: Simülasyon Metropolis–Hastings algoritmasını kullanır. Deneme spini, mevcut spin etrafındaki bir "koni" içinde $\mathbf{S}' = (\mathbf{S} + \delta\,\mathbf{r})/|\mathbf{S} + \delta\,\mathbf{r}|$ şeklinde üretilir. Burada $\mathbf{r}$, birim küre üzerinde düzgün dağılmış bir vektördür ve Marsaglia yöntemi (1972) ile elde edilir: $u, v \in [-1, 1]$, $s = u^2 + v^2 < 1$ olana kadar çekilir ve $\mathbf{r} = (2u\sqrt{1-s},\, 2v\sqrt{1-s},\, 1-2s)$ alınır. Bu öneri simetriktir, dolayısıyla ayrıntılı denge korunur. Koni genişliği $\delta$, kabul oranı yaklaşık %50 olacak şekilde her süpürmeden sonra yavaşça ayarlanır. Tamamen rastgele deneme spinleri düşük sıcaklıkta neredeyse hiç kabul edilmez (örneğin $T = 0.15$'te yaklaşık %3); koni hamleleri sistemin yapay olarak donmasını önler.
Güncelleme Sıklığı (MCS/kare): Zaman, Monte Carlo süpürmesi (MCS) cinsinden ölçülür. Bu parametre, ekrana iki çizim arasında kaç MCS yapılacağını belirler. Daha yüksek bir değer dengeye MCS cinsinden daha hızlı ulaştırmaz; ancak gerçek (duvar saati) zamanında daha hızlı ilerlemeyi ve görünür kareler arasında daha az ilinti olmasını sağlar.
Termodinamik Hafıza: Simülasyon sırasında sıcaklık veya manyetik alan parametrelerinin değiştirilmesi, örgüyü rastgele bir duruma sıfırlamaz. Sistem mevcut mikro-durumunu korur ve yeni koşullar altında termodinamik evrimine devam eder. Bu, ani soğutma (thermal quench) veya tavlama (simulated annealing) süreçlerinin doğrudan simüle edilmesini sağlar. Ancak, görselleştirmedeki sadeliği korumak adına, parametre değişimlerinde manyetizasyon dizileri, topolojik yük dizileri ve grafikler otomatik olarak temizlenir.
Önerilen Deneyler: Sayfa, rastgele bir durumdan $T = 0.15$, $B_z = 0.5$, $D = 1$ (Néel) parametreleriyle başlar. Birkaç yüz MCS içinde $-z$ çekirdekli skyrmionlardan oluşan bir doku ve $|Q| \approx 15$–$20$ civarında bir topolojik yük oluşur; "Skyrmion örneği" düğmesi bu başlangıca geri döner. $B_z = 0$ ayarlayıp "Rastgele başlat" düğmesine basıldığında skyrmion yerine labirent benzeri sarmal (spiral) alanlar oluşur ve $Q$ sıfıra yakın kalır. Buna karşılık, var olan skyrmionlar alan sıfırlansa bile uzun süre yaşamaya devam eder; bu, topolojik korumanın doğrudan bir göstergesidir. $B_z$ artırıldıkça skyrmionlar küçülür ve sayıları azalır; $B_z \gtrsim 1$ olduğunda sistem ferromanyetik duruma geçer. Néel tipinde düzlem içi okların skyrmion merkezine göre radyal, Bloch tipinde ise teğetsel yöneldiğine dikkat edin.
Görselleştirme ve Renk Kodlaması: 3B spin vektörleri 2B örgü düzlemine yansıtılır. Bir hücrenin arka plan rengi $S_z$ bileşenini temsil eder (Pozitif $\to$ Koyu Mavi, Negatif $\to$ Kırmızı). İçerisindeki siyah ok ise $(S_x, S_y)$ izdüşümünü gösterir.