LLG Denkleminin Boyutsuzlaştırması: Simülasyonun sayısal kararlılığını artırmak için fiziksel zaman, $t_0 = (1+\alpha^2)/(\gamma J)$ birimiyle $\tau = t / t_0$ şeklinde boyutsuzlaştırılmıştır. Enerji ölçeklemesi ($J=1$) ile de Hamiltoniyenden türetilen etkin alan ($\mathbf{h}_i^{\text{eff}} = -\partial \mathcal{H} / \partial \mathbf{S}_i$) mikromanyetik sabitlerden arınmış boyutsuz bir forma dönüştürülmüştür. Deterministik presesyon ve sönümleme terimlerini içeren denklem şu formdadır:
$$ \frac{d\mathbf{S}_i}{d\tau} = -\mathbf{S}_i \times \mathbf{h}_i^{\text{toplam}} - \alpha \, \mathbf{S}_i \times (\mathbf{S}_i \times \mathbf{h}_i^{\text{toplam}}) $$
Burada $\mathbf{h}_i^{\text{toplam}} = \mathbf{h}_i^{\text{eff}} + \mathbf{h}_i^{\text{th}}$, etkin alan ile termal alanın toplamıdır.
Termal Gürültü ve Langevin Dinamiği: Sonlu sıcaklık ($T > 0$), etkin alana stokastik bir termal alan ($\mathbf{h}_i^{\text{th}}$) eklenerek modellenir. Bu beyaz gürültü uzayda ve zamanda ilintisizdir ve şiddeti Fluktuasyon–Disipasyon teoremi ile belirlenir. Gilbert formundaki
$$\langle h^{\text{th}}_\mu(t)\, h^{\text{th}}_\nu(t')\rangle = (2\alpha k_B T/\gamma\mu_s)\,\delta_{\mu\nu}\,\delta(t-t')$$
ilintisi, $\tau$ birimine geçildiğinde $1/(1+\alpha^2)$ çarpanı alır. Bu nedenle her $d\tau$ adımında eklenen Gauss gürültüsünün standart sapması şudur:
$$ \sigma = \sqrt{\frac{2\alpha T}{(1+\alpha^2)\, d\tau}} $$
Gauss dağılımlı sayılar, düzgün dağılımlı rastgele sayılardan Box–Muller dönüşümü ile üretilir. (Kontrol: etkileşimsiz spinler için bu kalibrasyon, $\langle S_z \rangle$ değerinin Langevin fonksiyonu $\coth(B/T) - T/B$'ye uymasını sağlar.)
Stratonovich İntegrasyonu: Termal alan, denklemde spin vektörüyle çarpımsal olarak bulunur (çarpımsal gürültü). Bu durumda stokastik integralin yorumu önemlidir: LLG denklemi fiziksel olarak Stratonovich anlamında yorumlanmalıdır. Basit Euler şeması Itô yorumuna karşılık gelir ve ek bir sürüklenme düzeltmesi olmadan yanlış denge dağılımına yol açar. Bu simülasyon, termal alanı her iki aşamada da sabit tutan Heun (tahminci–düzeltici) yöntemini kullanır; Heun yöntemi Stratonovich çözümüne yakınsar. Her aşamadan sonra spinler yeniden normalize edilerek $|\mathbf{S}_i| = 1$ korunur. Sonlu $d\tau = 0.02$ adımı küçük bir sistematik hata getirir; Boltzmann dağılımı $d\tau \to 0$ limitinde tam olarak elde edilir.
$C_6$ Geometrisi ve DMI Vektörleri: Simülasyon, her spinin $\pi/3$ açılarıyla dizilmiş 6 komşuya sahip olduğu üçgen örgü üzerine kuruludur. DMI vektörleri ($\mathbf{D}_{ij}$), komşulara yönelen birim vektörler ($\hat{u}_{ij}$) üzerinden tanımlanır: Néel tipi için $\mathbf{D}_{ij} = D(\hat{u}_{ij} \times \hat{z})$, Bloch tipi için $\mathbf{D}_{ij} = D\hat{u}_{ij}$. Hamiltoniyen
$$\mathcal{H} = -J\sum_{\langle i,j\rangle}\mathbf{S}_i\cdot\mathbf{S}_j - \sum_{\langle i,j\rangle}\mathbf{D}_{ij}\cdot(\mathbf{S}_i\times\mathbf{S}_j) - \sum_i \mathbf{B}\cdot\mathbf{S}_i$$
olduğundan, her spine etkiyen etkin alan şöyledir:
$$ \mathbf{h}_i^{\text{eff}} = -\frac{\partial \mathcal{H}}{\partial \mathbf{S}_i} = \sum_{j} \left( J\,\mathbf{S}_j + \mathbf{S}_j \times \mathbf{D}_{ij} \right) + \mathbf{B} $$
Sistem Boyutu ve Periyodik Sınırlar: Simülasyon, eşkenar dörtgen biçiminde $L \times L$ ($L=30$) bir örgüde, periyodik sınır koşullarıyla (PBC) çalışır. Karakteristik sarmal periyodu (örgü sabiti biriminde) $\lambda \sim 2\pi J/|D|$ mertebesindedir. $\lambda$ sistem boyutuna uymadığında periyodik sınırlar dokuyu sıkıştırır ya da gerer; bu durum skyrmionların şeklini, sayısını ve düzenini etkiler. Belirgin skyrmionlar için $\lambda$'nın $L$'den yeterince küçük olması gerekir.
Topolojik Yük ($Q$): Topolojik yük, $C_6$ örgüsünü tamamen kaplayan üçgenler (örgü noktası başına bir yukarı ve bir aşağı üçgen) üzerinden Berg–Lüscher tanımıyla hesaplanır. Her üçgenin üç spininin birim küre üzerinde gerdiği katı açı şöyledir:
$$ \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ırlar altında bu tanım her konfigürasyon için tam sayı verir; $Q$ grafiği yalnızca bir skyrmion örgü ölçeğinde doğduğunda ya da yok olduğunda basamak yapar. $B_z > 0$ iken çekirdeği $-z$ olan her skyrmion $Q = -1$ katkı verir. (Not: Küçük açılarda basit üçlü çarpım $\mathbf{S}_1\cdot(\mathbf{S}_2\times\mathbf{S}_3)$, katı açının yaklaşık iki katıdır; bütün üçgenler üzerinden toplandığında $Q$'yu yaklaşık iki kat fazla verir.)
Zaman Adımı ve Termodinamik Hafıza: "İntegrasyon/Kare" parametresi, ekrana iki çizim arasında kaç adet $d\tau = 0.02$ büyüklüğünde Heun adımı yapılacağını belirler; grafiklerin yatay ekseni boyutsuz zaman $\tau$'yu gösterir. Parametre değiştirildiğinde sistemin anlık spin durumu korunur ve evrim kaldığı yerden devam eder (ani soğutma ya da tavlama deneyleri yapılabilir); yalnızca zaman serisi grafikleri sıfırlanır.
Önerilen Deney: Sayfa, rastgele bir durumdan $T = 0.1$, $B_z = 1$, $D = 1$ (Néel) parametreleriyle başlar ve birkaç on $\tau$ içinde ferromanyetik arka plan içinde $-z$ çekirdekli skyrmionlar oluşur ("Skyrmion örneği" düğmesi bu başlangıca geri döner). "Ferromanyetik (+z) başlat" ile düzgün durumdan başlandığında ise aynı parametrelerde skyrmion oluşmaz: düzgün durum yarı kararlıdır ve küçük termal dalgalanmalar topolojik bariyeri aşmaya yetmez. Bu karşılaştırma, skyrmionların topolojik korumasını doğrudan gösterir.
Görselleştirme ve Renk Kodlaması: Örgü projeksiyonu, 3B spin vektörlerini 2B düzleme yansıtır: