Fiziksel Model ve Hamiltoniyen: Bu simülasyon, “$C_6$ Simetrili Örgüde Stokastik Landau-Lifshitz-Gilbert (LLG) Dinamiği” sayfasındaki modelin aynısını kullanır: $C_6$ simetrili üçgen örgüde birim uzunluklu klasik spinler ($|\mathbf{S}_i| = 1$), en yakın komşu değiş-tokuş etkileşimi, Dzyaloshinskii–Moriya etkileşimi (DMI) ve $z$ yönünde dış manyetik alan. Tek fark, alanın artık uzayda düzgün olmamasıdır: alanın şiddeti örgü satırına ($y$) bağlı olarak sinüzoidal biçimde değişir.
$$ \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 B_z(y_i)\, S_i^z, \qquad B_z(y) = \pm\left[B_0 + A \sin\!\left(\frac{2\pi y}{L}\right)\right] $$
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}$ alınır. Örgü satırları ekranda yatay olduğundan alan yalnızca dikey konuma bağlıdır: gradyan dikey, eş-alan çizgileri yataydır. Alanın en zayıf ve en güçlü olduğu satırlar ($y = 3L/4$ ve $y = L/4$) örgü üzerinde kesikli çizgilerle gösterilir; sağdaki grafik de $|B_z|$ profilini gösterir.
LLG Denklemi: Spinlerin zaman evrimi, “$C_6$ Simetrili Örgüde Stokastik Landau-Lifshitz-Gilbert (LLG) Dinamiği” sayfasındaki boyutsuz Landau–Lifshitz–Gilbert denklemiyle hesaplanır ($\tau = t/t_0$, $t_0 = (1+\alpha^2)/(\gamma J)$):
$$ \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}}), \qquad \mathbf{h}_i^{\text{eff}} = \sum_{j} \left( J\,\mathbf{S}_j + \mathbf{S}_j \times \mathbf{D}_{ij} \right) + B_z(y_i)\,\hat{z} $$
Burada $\mathbf{h}_i^{\text{toplam}} = \mathbf{h}_i^{\text{eff}} + \mathbf{h}_i^{\text{th}}$'dir. Termal alan “$C_6$ Simetrili Örgüde Stokastik Landau-Lifshitz-Gilbert (LLG) Dinamiği” sayfasındakiyle aynı biçimde tanımlanır; ancak bu sayfa varsayılan olarak $T = 0$'da çalışır.
Skyrmion Neden Hareket Eder? Skyrmionun çekirdeği alana zıt yöndedir. Bu nedenle skyrmionun Zeeman enerjisi, bulunduğu yerdeki alan şiddeti $|B|$ ile artar. Alan uzayda değiştiğinde skyrmionun enerjisi konumuna bağlı hale gelir ve skyrmion, bir parçacık gibi, daha zayıf alana doğru bir kuvvet hisseder:
$$ \mathbf{F} = -\nabla U(\mathbf{R}) $$
Burada $\mathbf{R}$ skyrmionun konumu, $U$ ise enerjisidir. Örgü üzerindeki mavi ok bu kuvvetin yönünü gösterir.
Thiele Denklemi ve Skyrmion Hall Etkisi: Şeklini koruyarak hareket eden bir skyrmionun hızı $\mathbf{v}$, Thiele denklemiyle ifade edilir:
$$ \mathbf{G} \times \mathbf{v} + \alpha\,\mathcal{D}\,\mathbf{v} = \mathbf{F}, \qquad \mathbf{G} = 4\pi Q\,\hat{z}, \qquad \mathcal{D} = \frac{1}{2}\int \left[ (\partial_x \mathbf{S})^2 + (\partial_y \mathbf{S})^2 \right] dx\,dy $$
İlk terim (jiroskopik terim) topolojik yük $Q$ ile orantılıdır ve kuvvete dik bir hız bileşeni üretir. İkinci terim sönümden gelir. Sonuç olarak skyrmion kuvvet yönünde değil, kuvvetle şu açıyı yapan bir doğrultuda hareket eder:
$$ \tan\theta_{\text{SkH}} = \frac{4\pi |Q|}{\alpha\,\mathcal{D}} $$
Bu, skyrmion Hall etkisidir. Sönüm küçükken ($\alpha \ll 1$) açı $90^\circ$'ye yaklaşır: skyrmion daha zayıf alana doğru "yuvarlanmak" yerine neredeyse eş-alan çizgileri boyunca, yani yatay olarak kayar. Sapmanın yönü yalnızca $Q$'nun işaretine bağlıdır. Bu nedenle Néel ve Bloch skyrmionları aynı yörüngeyi izler; alan ters çevrildiğinde ($Q \to -Q$) ise sapma yönü tersine döner. Hareketi skyrmionun kiralitesi değil, topolojisi belirler.
Ölçüm: Skyrmion çekirdeği her karede $S_z$'nin yerel minimumu çevresindeki ağırlık merkezi olarak belirlenir. Hız, son $\Delta\tau = 10$ süresindeki yer değiştirmeden hesaplanır (yeşil ok). Kontrol panelinde ölçülen Hall açısı, Thiele denkleminin tahminiyle yan yana gösterilir. Tahmin için $Q$ (Berg–Lüscher tanımı) ve $\mathcal{D}$, o anki spin konfigürasyonundan sonlu farklarla hesaplanır. Varsayılan parametrelerde ölçülen açı ile tahmin birkaç derece içinde uyuşur (örneğin $\alpha = 0.1$ için $83^\circ$ ve $84^\circ$).
Periyodik Alan Profili ve Manyetik Vadi: Sabit (doğrusal) bir gradyan periyodik sınır koşullarıyla uyumlu değildir; bu yüzden alan sinüzoidal seçilmiştir. Skyrmion, gradyanın en büyük olduğu orta satırdan ($y = L/2$) başlar ve yatay yönde hareket ederek periyodik sınırdan dolanır. Hızın küçük bir bileşeni kuvvet yönündedir; bu yüzden skyrmion her turda biraz daha zayıf alana doğru kayar. Sonunda alanın en zayıf olduğu satıra ("manyetik vadi", $y = 3L/4$) ulaşır; orada gradyan sıfır olduğu için durur. Bu yavaş iniş, sağdaki profil grafiğinde sarı noktanın hareketiyle izlenebilir.
Başlangıç Durumu ve Etkileşim: "Tek skyrmion" düğmesi, ferromanyetik arka plan içinde tek bir skyrmion oluşturur ve onu gradyan uygulanmadan önce düzgün $B_0$ alanında güçlü sönümle ($\alpha = 0.5$, $\Delta\tau = 30$) dengeye durultur. Örgüye tıklayarak istediğiniz yere yeni skyrmionlar ekleyebilirsiniz. Yeni skyrmion, mevcut bir skyrmiona yaklaşık 4.5 örgü sabitinden daha yakın olamaz; daha yakın bir tıklama iki skyrmionu üst üste bindirirdi. "Alanı ters çevir" düğmesi alanı ve tüm spinleri birlikte ters çevirir; bu, $Q$'nun işareti değişmiş bir eş sistem verir.
Skyrmionlar Arası Etkileşim: Simülasyon tüm spinleri çözdüğü için skyrmionlar arasındaki etkileşim ayrıca eklenmemiştir; kendiliğinden ortaya çıkar. İki skyrmion birbirine yaklaştığında kuyrukları üst üste biner ve aralarında itici bir kuvvet oluşur. Bu kuvvet birkaç örgü sabiti içinde üstel olarak azalır. Thiele denklemine göre skyrmionlar bu kuvvete de doğrudan uzaklaşarak değil, büyük ölçüde ona dik hareket ederek tepki verir. Sonuç olarak düzgün bir alanda ($A = 0$) birbirine yakın iki skyrmion ortak merkezleri etrafında dönmeye başlar ve bu sırada yavaşça birbirinden uzaklaşır; dönme yönü $Q$'nun işaretine bağlıdır. Gradyanlı alanda ise skyrmionlar dikey yönde de birbirini iter: alt satırdaki skyrmion, tek başına olsaydı yukarı doğru vadiye kayacakken, komşusu tarafından aşağı itilebilir. Birden fazla skyrmion varken kontrol panelinde ölçülen Hall açısı yalnızca gradyan kuvvetini değil, etkileşimi de yansıtır.
Sayısal Yöntem: Integrasyon, “$C_6$ Simetrili Örgüde Stokastik Landau-Lifshitz-Gilbert (LLG) Dinamiği” sayfasındaki Heun (tahminci–düzeltici) yöntemiyle $d\tau = 0.02$ adımıyla yapılır. "İntegrasyon/Kare" parametresinin varsayılan değeri 20'dir; bu değerle skyrmion saniyede yaklaşık iki örgü sabiti ilerler.
Görselleştirme: Renk kodlaması “$C_6$ Simetrili Örgüde Stokastik Landau-Lifshitz-Gilbert (LLG) Dinamiği” sayfasıyla aynıdır:
Soldaki grafik, izlenen skyrmionun başlangıca göre yatay ve dikey yer değiştirmesini gösterir (periyodik sınırdan geçişler açılmış olarak). Sağdaki grafik $|B_z|$ profilini ve skyrmionun bulunduğu satırı gösterir.
Önerilen Deneyler:
Sınırlamalar: Bu bir tanıtım modelidir. $D = J$ seçimiyle skyrmionlar yalnızca birkaç örgü sabiti büyüklüğündedir ve ayrık örgü etkileri içerir; Thiele denklemi skyrmionun şeklini koruduğunu varsayar. Alanın en zayıf olduğu bölgede alan çok düşerse ($B_0 - A \lesssim 0.7$) skyrmion uzayıp şeklini kaybedebilir. Deneylerde skyrmionlar genellikle alan gradyanıyla değil, spin-polarize akımla sürülür; skyrmion Hall etkisi bu tür deneylerde gözlenmiştir (Jiang ve ark., Nature Physics 13, 162 (2017); Litzius ve ark., Nature Physics 13, 170 (2017)). Akımla sürüş de aynı Thiele çerçevesine oturur; yalnızca kuvvet terimi değişir.