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:

• Néel Tipi (Arayüzey): Kiral Néel skyrmionlarını veya spin sikloidlerini kararlı kılar. $\mathbf{D}_{ij} = D (\hat{u}_{ij} \times \hat{z})$ olarak tanımlanır.
• Bloch Tipi (Hacimsel): Bloch skyrmionlarını veya spin sarmallarını (helix) kararlı kılar. $\mathbf{D}_{ij} = D \hat{u}_{ij}$ olarak tanımlanır.

Topolojik Yük ($Q$) ve Küsuratlı Değerleri: Topolojik invaryant $Q$, spin konfigürasyonunun birim küreyi kaç kez sardığını sayar. Ayrık örgü üzerinde komşu spinlerin üçlü skaler çarpımı kullanılarak sayısal olarak şu şekilde yaklaşık olarak hesaplanabilir:

$$ Q \approx \frac{1}{4\pi} \sum_{i,j} \mathbf{S}_{i,j} \cdot (\mathbf{S}_{i+1,j} \times \mathbf{S}_{i,j+1}) $$

Mutlak sıfırda ve sürekli termodinamik limitte, $Q$ kesinlikle kuantize olmuş bir tamsayı ($0, \pm 1, \pm 2, \dots$) olmalıdır. Ancak bu simülasyonda, topolojik yük zaman serisi sıklıkla tamsayı olmayan değerler gösterebilir. Bu sorun iki temel nedenden kaynaklanır:

1. Uyumluluk ve periyodik sınır koşulları Kesilmesi: Skyrmion büyüklüğü ($\lambda \approx 2\pi J / |D|$), sonlu örgü boyutuyla ($L=30$) uyumlu (commensurate) değilse, Periyodik Sınır Koşulları (PBC) skyrmionu sınırlarda yapay olarak keser. Tamamlanmamış, engellenmiş bir topolojik doku matematiksel olarak kesirli bir $Q$ değeri üretir. Net skyrmionik dokuları gözlemleyebilmek için $|D|$ değeri $\lambda \le L$ koşulunu sağlayacak büyüklükte seçilmelidir.
2. Sonlu Farklar Yaklaşımı: Üçlü skaler çarpım, sürekli uzaysal integralin bir yaklaşımıdır. DMI güçlü olduğunda veya termal dalgalanmalar ($T > 0$) yüksek olduğunda, komşu spinler arasındaki açı büyür. Keskin spin gradyanlarının olduğu bu rejimlerde bu yaklaşım doğru sonuçlar üretemez. Kesin tamsayılar elde etmek için hesaplama açısından çok daha maliyetli geometrik yöntemler (örn. Berg-Lüscher) gereklidir.

Monte Carlo Algoritması ve Marsaglia Örneklemesi: Simülasyon Metropolis-Hastings algoritmasını kullanır. Birim küre üzerinde düzgün rastgele deneme spinleri üretmek için Marsaglia'nın reddetme (rejection) yöntemi (1972) kullanılır. İki değişken $u, v \in [-1, 1]$, $s = u^2 + v^2 < 1$ olana kadar üretilir. Yeni spin bileşenleri $S_x = 2u\sqrt{1-s}, S_y = 2v\sqrt{1-s}, S_z = 1-2s$ şeklinde elde edilir.

Güncelleme Sıklığı (MCS/kare): Zaman, Monte Carlo Adımı (MCS) cinsinden ölçülür. Tarayıcının görsel çizim darboğazını aşmak için, bu parametre fiziksel Markov zincirini görsel kare (frame) hızından ayırır. Daha yüksek bir değer, sistemi faz uzayında hızla ilerletir, görünür kareler arasındaki otokorelasyonu en aza indirir ve termodinamik yakınsamayı hızlandırır.

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.

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.