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}}) $$
Termal Gürültü ve Langevin Dinamiği: Sonlu sıcaklık ($T > 0$), Fluktuasyon-Disipasyon teoremi gereği etkin alana stokastik bir termal gürültü alanı ($\mathbf{h}_i^{\text{th}}$) eklenerek modellenir. Bu beyaz gürültü uzayda ve zamanda sıfır korelasyonludur. Zamanın kesiklileştirilmesi ile gerçekleştirilen integrasyonda, her adımda eklenen rastgele Gaussian gürültünün standart sapması $\sigma = \sqrt{2 \alpha T / d\tau}$ ile belirlenir. Tarayıcı tabanlı sözde-rastgele sayı üreteçlerinin zayıflığını aşmak için Box-Muller dönüşümü kullanılarak normal (Gauss) dağılımlı 3B vektörler üretilir.
Stratonovich İntegrasyonu: Termal alan, denklemde spin vektörüyle çarpımsal (multiplicative) olarak bulunur. Euler yöntemi (Itô limiti) kullanıldığında sistemin birim spin kısıtı ($|\mathbf{S}|=1$) bozulur. Fiziksel Boltzmann dağılımına yakınsaması için Stratonovich kalkülüsü gereklidir. Algoritma, Heun (Predictor-Corrector) integratörünü kullanır ve termal gürültü vektörünü ($\mathbf{h}_i^{\text{th}}$) her iki integrasyon aşamasında da sabit tutarak termodinamik dengeyi garanti eder.
$C_6$ Geometrisi ve DMI Vektörleri: Simülasyon, her spinin $\pi/3$ açılarıyla dizilmiş 6 komşuya sahip olduğu üçgensel (triangular) örgü üzerine inşa edilmiştir. DMI vektörleri ($\mathbf{D}_{ij}$) komşuların birim vektörleri ($\hat{u}_{ij}$) üzerinden analitik olarak izdüşümlenir. Néel tipi DMI $\mathbf{D}_{ij} = D(\hat{u}_{ij} \times \hat{z})$ ve Bloch tipi DMI $\mathbf{D}_{ij} = D\hat{u}_{ij}$ yönelimlerini takip eder.
Sistem Boyutu, PBC ve Topolojik Frustrasyon: Simülasyon, eşkenar dörtgen (rhombus) formunda kaplanan $L \times L$ boyutlarında ($L=30$) bir örgüde çalışır. Periyodik Sınır Koşulları (PBC) kullanılarak kenar etkileri elimine edilmiştir. Ancak, içsel skyrmion adım uzunluğu ($\lambda \approx 2\pi J / |D|$) sonlu örgü boyutuna sığmadığında, PBC skyrmionu yapay olarak keser ve topolojik frustrasyona neden olur. Net skyrmionik fazlar gözlemlemek için $|D|$ değeri $\lambda \le L$ koşulunu sağlamalıdır.
Topolojik Yük ($Q$) ve Sonlu Farklar Sapması: Topolojik invariant, $C_6$ örgüsünü tam kaplayan (tiling) üçgen alt-hücrelerin skaler üçlü çarpımları toplanarak kesikli formda hesaplanır:
$$ Q \approx \frac{1}{4\pi} \sum_{\triangle} \mathbf{S}_1 \cdot (\mathbf{S}_2 \times \mathbf{S}_3) $$
Sürekli uzayda $Q$ mutlak kuantize bir tamsayıdır. Ancak bu simülasyondaki skaler üçlü çarpım sonlu farklar yaklaşımı ile ele alınır. Yüksek sıcaklık dalgalanmalarında veya DMI'ın alan vektörleri arasındaki açıyı dramatik şekilde büktüğü keskin gradyan rejimlerinde bu yaklaşım yetersiz kalır ve $Q$ grafiğinde kesirli (fractional) sapmalar gözlemlenir.
Zaman Adımları (dt) ve Termodinamik Hafıza: İstemci taraflı işleme yükünü hafifletmek için fiziksel zaman ( $\tau$ ) ile görselleştirme kareleri (frames) birbirinden ayrılmıştır. "İntegrasyon/Kare" parametresi, ekrana çizim yapılmadan önce arka planda kaç adet $d\tau = 0.02$ büyüklüğünde Heun entegrasyonu yapılacağını belirler. Parametre değişimlerinde sistemin anlık spin mikro-durumu korunarak termodinamik evrim kaldığı yerden devam eder (simulated annealing), ancak grafiksel gösterimdeki sadeliği korımak için istatistik havuzları sıfırlanır, yani zaman serisi grafikleri parametreler değiştirildiğinde yeniden başlar.
Görselleştirme ve Renk Kodlaması: Örgü projeksiyonu, 3B spin vektörlerini 2B düzleme yansıtır: