Kuantum spin dalgaları · magnonlar

Magnon Laboratuvarı

Klasik spin dalgasının kuantum karşılığı magnondur: kristale yayılmış tek bir spin kuantumu. Bu sayfadaki simülasyonlar magnonu en basit yerden kurar. Önce ferromanyette tek magnon, sonra iki magnonun etkileşimi ve bağlı durumları, son olarak da Néel durumunun kuantum dalgalanmalarıyla eridiği ve magnonların yerini spinonlara bıraktığı antiferromanyet. Her simülasyon kesin sayısal sonucu yaklaşık teoriyle yan yana gösterir. Aşağıdaki seçiciden simülasyonu seçin.

\[ H = -J\sum_{n}\mathbf S_n\cdot\mathbf S_{n+1} \;-\; h\sum_n S^z_n ,\qquad S=\tfrac12 \]

\(J>0\) ferromanyetik değiş-tokuş · \(h\): \(\hat z\) yönünde alan · \(N\) siteli halka (periyodik sınır) · \(\hbar=a=1\), zaman \(\hbar/J\) biriminde. Simülasyon tek-magnon sektöründe Schrödinger denklemini kesin çözer (\(N\times N\) problem); hiçbir yaklaşıklık yoktur.

SimülasyonMagnon Laboratuvarı

1 · Taban durumu. Bütün spinler ↑. Bu durum Hamiltoniyenin kesin özdurumudur: hiçbir şey hareket etmez, magnon sayısı sıfırdır. Bir siteye tıklayarak o spini çevirin.

periyodik halka · N = 30

renk = genliğin fazı arg c_n · doygunluk = |c_n|

Şekil 1 · Magnon enerjisi ω_k ve durumun magnon içeriği

Şekil 1. Eğri \(\omega_k=h+J(1-\cos k)\). Gri noktalar halkada izinli \(k=2\pi m/N\) değerleridir. Renkli noktaların büyüklüğü, o anki durumun \(|k\rangle\) bileşenindeki ağırlıktır, \(|\tilde c_k|^2\). Tek spin çevirmede bütün \(k\)'lar eşit ağırlıklıdır; magnon tek bir noktadır; paket \(k_0\) çevresinde dar bir bölgedir. Bir noktaya tıklamak o magnonu seçer.

Şekil 2

Kısayollar: Boşluk oynat/duraklat · R sıfırla · +/− hız · F tam sayfa

Hamiltoniyen

Değiş-tokuş J

Alan h

h bütün magnon enerjilerini aynı miktarda kaydırır; olasılıkları değil yalnızca fazları değiştirir.

Durum

Site sayısı N

Dalga sayısı k = 2πm/N

Paket genişliği σ

Süperpozisyon

Tek magnon ağırlığı p

\(|\psi\rangle=\sqrt{1-p}\,|F\rangle+\sqrt p\,|k\rangle\)

Çizim büyütmesi

Okların enine bileşeni yalnızca çizimde büyütülür.

Koni açısı θ

Zaman

Hız

1×: gerçek zamanın her saniyesi 2 ħ/J.

Hesap: genlikler \(\tilde c_k(t)=\tilde c_k(0)\,e^{-i\omega_k t}\) ile kesin olarak ilerletilir ve \(c_n=N^{-1/2}\sum_k \tilde c_k e^{ikn}\) ile siteye dönülür. Sayısal hata yalnızca kayan nokta mertebesindedir.

1 · Uzak iki magnon. Halkanın karşıt noktalarında iki spin çevrildi. Her biri tek magnon gibi yayılır. Çubuklar sitedeki çevrik spin sayısı \(\langle n_i\rangle\)'dir (toplam 2). Halkalar iki etkileşimsiz magnonun öngörüsüdür. Sağdaki harita iki çevrik spinin ortak konum olasılığı \(P(n_1,n_2)\)'dir. Magnonlar karşılaşana kadar sonuç etkileşimsiz öngörüyle aynıdır; fark ancak karşılaştıklarında başlar. Zincire ya da haritaya tıklayarak başlangıç yerlerini değiştirin.

periyodik halka

köşegen: aynı site (yasak)

Şekil 1 · Göreli uzaklık dağılımı P(r)

Şekil 2

Kısayollar: Boşluk oynat/duraklat · R sıfırla · +/− hız · F tam sayfa

Hamiltoniyen

Değiş-tokuş J

Anizotropi Δ

Δ = 1 Heisenberg, Δ = 0 XY (yalnız sert-çekirdek), Δ ≫ 1 Ising.

Alan h

h iki-magnon enerjilerini 2h kaydırır; olasılıkları değiştirmez.

Durum

Site sayısı N

Toplam momentum K = 2πM/N

Durum (K bloğunda sıra)

0: en düşük enerji. Spektrumdaki bir noktaya tıklamak da seçer.

Paket genişliği σ

Zaman

Hız

1×: gerçek zamanın her saniyesi 2 ħ/J.

Hesap: her \(K\) için göreli koordinat \(r=n_2-n_1\) üzerindeki \(\sim N/2\) boyutlu gerçel simetrik matris Jacobi yöntemiyle köşegenleştirilir; zaman evrimi \(e^{-iE t}\) ile kesindir. Toplam boyut \(N(N-1)/2\).

1 · Néel durumu özdurum değildir. Simülasyon klasik taban durumu \(|\!\uparrow\downarrow\uparrow\downarrow\cdots\rangle\)'den başlar ve zinciri kesin olarak zamanda ilerletir. Çubuklar \(\langle S^z_i\rangle\), kesikli çizgiler Néel değerleri \(\pm S\)'dir. Ferromanyetteki \(|\!\uparrow\uparrow\cdots\rangle\)'nin aksine, değiş-tokuşun \(S^+_iS^-_{i+1}\) terimi komşu ↑↓ çiftlerini çevirir; bu yüzden düzen kendiliğinden erir. Şekil 2'de bağ enerjisinin \(zz\) kısmından \(xy\) kısmına aktığı görülür; toplamı sabittir.

S = 1/2 · N = 12

Şekil 1

Şekil 2

Kısayollar: Boşluk oynat/duraklat · R sıfırla · +/− hız · F tam sayfa

Zincir (tam köşegenleştirme)

Spin S

Site sayısı N

Lineer spin dalgası teorisi

Örgü

Spin S

Örgü boyutu (k toplamı)

\(\Delta S\) sonlu örgüde \(\mathbf k\) toplamıyla hesaplanır; \(\gamma_{\mathbf k}=\pm1\) modları (düzgün dönme) dışarıda bırakılır.

Spektrum

Seçilen q = 2πn/N

Genişletme η

Sonlu zincirin spektrumu kesiklidir; η yalnızca çizim için Lorentz genişletmesidir.

Zaman

Hız

1×: gerçek zamanın her saniyesi 2 ħ/J.

Okuma kılavuzu

Çubuklar, renkler ve oklar

Çubuk yüksekliği

1–4. adımlarda her sitedeki çubuk \(P_n=|c_n|^2=\tfrac12-\langle S^z_n\rangle\)'dir: o sitede çevrik spin bulunma olasılığı. Sağ eksen aynı büyüklüğü \(\langle S^z_n\rangle\) olarak okur. Dikey ölçek kendiliğinden ayarlanır; değerler sol eksendedir.

Renk

Renk, karmaşık genliğin fazı \(\arg c_n\)'dir; çark Spin Dalgası Laboratuvarı sayfasındakiyle aynıdır (0° kırmızı, 90° sarı-yeşil, 180° camgöbeği, 270° mor-mavi). Enerji taban durumundan ölçüldüğü için faz \(\omega_k\) hızıyla döner. 5. adımda aynı renk \(\arg\langle S^+_n\rangle\)'dir, yani Spin Dalgası Laboratuvarı sayfasındaki presesyon fazının birebir karşılığıdır.

Fazın anlamı

Tek başına bir durumun küresel fazı ölçülemez; yalnızca siteler arası faz farkları (\(k\)) ve iki durumun süperpozisyonundaki göreli faz anlamlıdır. 1–4. adımlardaki renk, \(|F\rangle\) ile süperpoze edildiğinde ortaya çıkacak fazdır. Bu süperpozisyon 5(a)'da gerçekten yapılır.

Oklar (yalnız 5. adım)

Ok \(\langle\mathbf S_n\rangle\) beklenen değeridir. Kuantum spinin klasik bir yönü yoktur; ok yalnızca ortalamadır. 5(a)'da okun boyu \(1/2\)'den kısadır. Bu, sitenin indirgenmiş durumunun karışık olduğunu, yani sitenin zincirin geri kalanıyla dolanık olduğunu gösterir.

Türetme

Taban durumundan magnona

Adım 1

Ferromanyetik taban durumu

\(\mathbf S_n\cdot\mathbf S_{n+1}=S^z_nS^z_{n+1}+\tfrac12(S^+_nS^-_{n+1}+S^-_nS^+_{n+1})\). \(|F\rangle=|\!\uparrow\uparrow\cdots\uparrow\rangle\) üzerinde \(S^+\) terimleri sıfır verir, \(S^zS^z\) ise \(\tfrac14\) verir:

\[ H|F\rangle=E_0|F\rangle,\qquad E_0=-\tfrac{NJ}{4}-\tfrac{Nh}{2}. \]

\(|F\rangle\) kesin özdurumdur ve içinde kuantum dalgalanması yoktur. Antiferromanyette Néel durumu böyle değildir: \(S^+_nS^-_{n+1}\) terimi ↓↑ çiftini ↑↓ yapar, bu yüzden AFM taban durumu Néel'den farklıdır.

Adım 2

Çevrik spin bir parçacıktır

\(|n\rangle=S^-_n|F\rangle\) durumunda iki bağ \(\tfrac14\)'ten \(-\tfrac14\)'e iner (\(+J\)), alan terimi \(+h\) ekler; değiş-tokuşun çevirme terimi çevrik spini komşuya taşır:

\[ \begin{aligned} H|n\rangle&=(E_0+J+h)\,|n\rangle\\ &\quad-\tfrac J2\big(|n+1\rangle+|n-1\rangle\big). \end{aligned} \]

Bu, sıkı-bağ modelindeki bir parçacığın denklemidir. \(S^z_{\rm top}\) korunduğu için çevrik spin sayısı tam olarak 1 kalır. \(j\) sitesinden başlayan sonsuz zincir çözümü

\[ c_n(t)=i^{\,n-j}J_{n-j}(Jt)\,e^{-i(J+h)t}, \] \[ \langle (n-j)^2\rangle=\tfrac{(Jt)^2}{2}. \]

Yayılma balistiktir (\(\propto t\)), difüzif (\(\propto\sqrt t\)) değildir.

Adım 3

Magnon: momentum özdurumu

Öteleme simetrisi yüzünden özdurumlar düzlem dalgalardır:

\[ |k\rangle=\frac{1}{\sqrt N}\sum_n e^{ikn}|n\rangle, \] \[ \omega_k=E_k-E_0=h+J(1-\cos k), \]

\(k=2\pi m/N\). Bu, \(2JS(1-\cos k)\) ifadesinin \(S=\tfrac12\) halidir ve klasik lineer spin dalgası frekansıyla birebir aynıdır. FM'de tek magnon için kuantum ve klasik dağılım çakışır. Küçük \(k\)'da \(\omega_k\simeq h+\tfrac J2k^2\): magnon, etkin kütlesi \(1/J\) olan bir parçacıktır.

\(k=0\) özeldir: \(|0\rangle\propto S^-_{\rm top}|F\rangle\), yani döndürülmüş taban durumudur ve \(h=0\)'da enerjisi sıfırdır (Goldstone modu). \(k\neq0\) magnonların toplam spini \(N/2-1\)'dir. Tek spin çevirme bütün magnonların eşit süperpozisyonudur: \(|j\rangle=N^{-1/2}\sum_k e^{-ikj}|k\rangle\).

Adım 4

Dalga paketi ve grup hızı

Tek-magnon sektöründe konum operatörü \(X\), Heisenberg resminde tam olarak \(X(t)=X(0)+V\,t\) ile ilerler; hız operatörü \(V=\partial\omega_k/\partial k=J\sin k\)'dir. Gerçel ve simetrik bir zarf için

\[ \langle x\rangle(t)=x_0+\langle v\rangle t,\qquad \sigma^2(t)=\sigma_0^2+\mathrm{Var}(v)\,t^2, \]

ortalamalar \(|\tilde c_k|^2\) dağılımı üzerinden alınır. Dar bir paket için \(\langle v\rangle\simeq J\sin k_0\) ve \(\mathrm{Var}(v)\simeq J^2\cos^2k_0\,\sigma_k^2\). En büyük hız \(k_0=\pi/2\)'dedir; orada dağılımın eğriliği sıfır olduğu için paket en yavaş yayılır.

Adım 5

Klasik spin dalgası nereden çıkar?

\(S^+_n\) magnon sayısını bir azaltır. Bu yüzden magnon sayısı belli bir durumda \(\langle S^\pm_n\rangle=0\) olur. (a) \(|\psi\rangle=\sqrt{1-p}\,|F\rangle+\sqrt p\,|k\rangle\) için

\[ \langle S^+_n\rangle=\sqrt{\tfrac{p(1-p)}{N}}\;e^{i(kn-\omega_k t)},\qquad \langle S^z_n\rangle=\tfrac12-\tfrac pN . \]

Bu, klasik presesyon frekansıyla dönen bir spin dalgasıdır, ama genliği \(1/\sqrt N\) kadar küçüktür ve \(|\langle\mathbf S_n\rangle|<\tfrac12\) olur. (b) \(k=0\) koherent durumu \(\bigotimes_n\big(\cos\tfrac\theta2|\!\uparrow\rangle+e^{i\phi}\sin\tfrac\theta2|\!\downarrow\rangle\big)\) bir ürün durumudur. Magnon sayısı binom dağılımlıdır:

\[ \langle n\rangle=NS(1-\cos\theta),\qquad \frac{\Delta n}{\langle n\rangle}=\sqrt{\frac{1-s}{Ns}},\quad s=\sin^2\tfrac\theta2 . \]

\(N\to\infty\)'da göreli dalgalanma sıfıra gider: klasik dalga, fazı keskin, sayısı belirsiz çok sayıda magnondur. \(S=\tfrac12\)'de \(k\neq0\) için çok-magnonlu koherent durum kesin değildir, çünkü magnonlar etkileşir. Bu, İki magnon simülasyonunun konusudur.

Köprü

Spin Dalgası Laboratuvarı sayfasıyla karşılaştırma

Büyüklük	Klasik (LLG, lineer)	Kuantum (tek magnon)
Dağılım	\(\omega_k=h+2JS(1-\cos k)\)	Aynı (\(S=\tfrac12\)); FM'de tek-magnon için kesin
Tek spin uyarımı	Enine genlik \(\psi_n\propto i^nJ_n(2JSt)\)	Olasılık genliği \(c_n\propto i^nJ_n(Jt)\): aynı matematik, farklı anlam
Uyarım genliği	Sürekli: koni açısı \(\theta_0\)	Kesikli: magnon sayısı \(n\), \(\langle n\rangle=NS(1-\cos\theta)\)
Enine bileşen	Her zaman tanımlı	Sayı özdurumunda \(\langle S^\perp\rangle=0\); yalnız süperpozisyonda sıfırdan farklı
Spin uzunluğu	\(\|\mathbf S\|=S\)	\(\|\langle\mathbf S_n\rangle\|\le S\), \(\langle\mathbf S^2\rangle=S(S+1)\)
Termal dağılım	Eşbölüşüm \(k_BT/\omega_k\)	Bose–Einstein \(1/(e^{\omega_k/k_BT}-1)\) → Bloch \(T^{3/2}\)

Eleştirel notlar

Bu simülasyonun kapsamı ve sınırları

Neden kesin?

FM taban durumu ve tek-magnon durumları Heisenberg modelinin kesin özdurumlarıdır. Bu yüzden burada Holstein–Primakoff ya da lineer spin dalgası yaklaşıklığı gerekmez. Bu, FM'ye özgü bir kolaylıktır. AFM'de taban durumu bile kesin bilinmez ve magnonlar ancak yaklaşık olarak (LSWT) tanımlanır.

Magnon bozon mudur?

Tek magnon düzeyinde \(b_k^\dagger=N^{-1/2}\sum_n e^{ikn}S^-_n\) bir bozon gibi davranır. \(S=\tfrac12\)'de bir siteye iki çevirme sığmaz (\((S^-_n)^2=0\)); ayrıca iki magnon \(S^zS^z\) terimi üzerinden etkileşir ve 1D'de bağlı durumlar oluşturur (Bethe, 1931). Bozon resmi seyrek magnon limitinde (\(\langle n\rangle/N\ll 2S\)) geçerlidir.

Sonlu halka

\(k\) kuantizedir, spektrum kesiklidir. Tek spin uyarımının cephesi \(t_s=N/2J\)'de halkanın karşı noktasına ulaşır ve kendisiyle girişim yapar. Bundan sonra sonsuz zincir formülleri (Bessel, \(\sigma(t)\)) geçerli değildir, ama simülasyon halka için kesin kalır. Uzun zamanda kısmi canlanmalar görülür.

Kuantum mu, klasik mi?

Tek magnonun dinamiği klasik lineer spin dalgasıyla matematiksel olarak aynıdır. Kuantum doğası, uyarımın kesikliliğinde (\(\Delta S^z_{\rm top}=-1\)), \(\langle S^\perp\rangle=0\) olmasında ve dolanıklıkta görünür. Bu nedenle tek bir magnonun “dalgası” klasik bir ok dalgası değil, bir olasılık genliği dalgasıdır.

Önerilen alıştırmalar

Balistik yayılma. Adım 2, \(N=60\). Şekil 2'de \(\sigma(t)\)'nin \(Jt/\sqrt2\) doğrusunu izlediğini gösterin. \(J\)'yi iki katına çıkarın: cephe hızı ne olur? Difüzyonla (\(\sigma\propto\sqrt t\)) farkı nedir?
Alan neyi değiştirir? Adım 2'de \(h\)'yi 0'dan 1'e çıkarın. \(P_n\) neden hiç değişmez? (İpucu: \([H_{\rm Zeeman},H_{\rm ex}]=0\) ve sektörde \(S^z_{\rm top}\) sabittir.)
Goldstone modu. Adım 3, \(h=0\), \(m=0\). Renk neden dönmez? \(|k=0\rangle\)'ın \(S^-_{\rm top}|F\rangle\) ile orantılı olduğunu gösterin ve \(k\neq0\) magnonun toplam spinini bulun.
Grup hızı. Adım 4, \(N=60\), \(\sigma=4\). \(k_0=\pi/4,\pi/2,3\pi/4\) için ölçülen hızı \(J\sin k_0\) ile karşılaştırın. Paket hangi \(k_0\)'da en yavaş yayılır? Neden?
Ok neden kısa? Adım 5(a), \(p=0.5\). \(N=20\) ve \(N=60\) için enine genliği okuyun ve \(\propto1/\sqrt N\) olduğunu doğrulayın. \(|\langle\mathbf S_n\rangle|<\tfrac12\) olmasını tek sitenin indirgenmiş yoğunluk matrisi üzerinden açıklayın.
Klasik dalgada kaç magnon var? Spin Dalgası Laboratuvarı sayfasındaki \(N=80\), \(\theta_0=5^\circ\) düzlem dalgası, \(S=\tfrac12\)'de ortalama kaç magnona karşılık gelir? Sonucu (\(\langle n\rangle\approx0.15\)) Adım 5(b) ile tartışın: “küçük genlikli klasik dalga” kuantumda neye karşılık gelir?

Okuma kılavuzu

Çubuklar, harita ve spektrum

Çubuklar

\(\langle n_i\rangle=\tfrac12-\langle S^z_i\rangle\): \(i\) sitesinde çevrik spin bulunma olasılığı. Toplam her an tam 2'dir. Halkalar, iki magnonun hiç etkileşmediği varsayımıyla hesaplanan yoğunluktur: \(|\varphi_a(i,t)|^2+|\varphi_b(i,t)|^2\), burada \(\varphi\) Tek magnon simülasyonundaki çözümdür.

İki parçacık haritası

Yatay eksen birinci, dikey eksen ikinci çevrik spinin yeri. Harita simetriktir, çünkü magnonlar özdeştir. Köşegen (\(n_1=n_2\)) gridir: \(S=\tfrac12\)'de bir siteye iki çevirme sığmaz. Köşegenin hemen yanındaki hücreler komşu çiftlerdir. Bağlı durum köşegene yapışık bir şerit, serbest magnonlar ise köşegenden uzaklaşan bir desen verir. Periyodik sınır yüzünden sol üst ve sağ alt köşeler de köşegene komşudur.

Göreli uzaklık

\(P(r)\), iki çevrik spin arasındaki halka uzaklığının (\(r=1,\dots,N/2\)) olasılığıdır. Bağlı durumda \(P(r)\) üstel olarak azalır; serbest magnonlarda uzaklık zamanla büyür.

Spektrum

3. adımda her nokta bir iki-magnon özdurumudur. Gölgeli bant, iki serbest magnonun mümkün enerji toplamlarıdır. Noktaların çoğu bu bandın içindedir. Bandın dışına düşen her nokta etkileşimin doğrudan izidir.

Türetme

İki magnon problemi

Adım 1

Etkileşimin iki kaynağı

Tek çevrik spin iki bağı “bozar” ve enerjiyi \(J\Delta\) artırır. Birbirinden uzak iki çevrik spin dört bağ bozar: \(2J\Delta\). Yan yana iki çevrik spin ise yalnızca iki bağ bozar:

\[ \underbrace{\uparrow\downarrow\uparrow\cdots\uparrow\downarrow\uparrow}_{2J\Delta}\qquad \underbrace{\uparrow\downarrow\downarrow\uparrow}_{J\Delta}. \]

Komşu olmak enerjiyi \(J\Delta\) kadar düşürür: çekici bir etkileşim. İkinci kaynak kinematiktir: \((S^-_i)^2=0\), yani bir siteye iki çevirme sığmaz (sert çekirdek). Hamiltoniyen yerel olduğu için uzaktaki iki magnon kesin olarak bağımsız hareket eder.

Adım 2

Kütle merkezi ve göreli hareket

Toplam momentum korunur. Bu yüzden

\[ \psi(n_1,n_2)=e^{iK(n_1+n_2)/2}f(r),\qquad r=n_2-n_1 \]

yazılınca problem göreli koordinatta tek boyutlu bir sıkı-bağ problemine dönüşür:

\[ \begin{aligned} E f(r)=&\big(2h+2J\Delta-J\Delta\,\delta_{r,1}\big)f(r)\\ &-J\cos\tfrac K2\,\big[f(r+1)+f(r-1)\big],\quad f(0)=0. \end{aligned} \]

Göreli parçacık \(r=0\)'da sert bir duvar, \(r=1\)'de \(-J\Delta\) derinliğinde bir kuyu görür. Atlama genliği \(J\cos(K/2)\) toplam momentuma bağlıdır. Halkada \(f(r)=e^{iKN/2}f(N-r)\) koşulu her \(K\) bloğunu \(\sim N/2\) boyuta indirir.

Adım 3

Süreklilik ve bağlı durum

Saçılma durumları \(f(r)\sim\sin(qr+\delta_q)\) biçimindedir. Enerjileri iki serbest magnonun enerjileri toplamıdır:

\[ E=\omega_{K/2+q}+\omega_{K/2-q}=2h+2J\Delta-2J\cos\tfrac K2\cos q . \]

Bağlı durum \(f(r)\propto x^{\,r-1}\) ile aranır. \(r=1\) denkleminden

\[ x=\frac{\cos(K/2)}{\Delta},\qquad E_b(K)=2h+J\Delta-\frac{J\cos^2(K/2)}{\Delta}, \]

ve bağlı durum \(|x|<1\), yani \(|\cos(K/2)|<\Delta\) olduğunda vardır. Heisenberg zincirinde \(E_b=2h+J\sin^2(K/2)\) olur. Bu, Bethe'nin 1931 sonucudur. Bağlı durum \(K\neq0\) için her zaman bandın altındadır. Ising limitinde bağlanma enerjisi \(J\Delta\)'ya, yani “bir bağ kazancına” gider.

Adım 4

Bağlı durum bir parçacıktır

Bağlı durumun kendi dağılım bağıntısı vardır, dolayısıyla kendi grup hızı ve kütlesi de vardır:

\[ v_b=\frac{dE_b}{dK}=\frac{J}{2\Delta}\sin K,\qquad E_b\simeq 2h+\tfrac J4K^2\ (\Delta=1). \]

Heisenberg zincirinde bağlı çiftin kütlesi \(2/J\), yani tek magnonun kütlesinin (\(1/J\)) tam iki katıdır. Yan yana iki çevrik spinden başlanırsa bağlı durumlara giden toplam pay

\[ W_b=\frac1{2\pi}\int dK\,\big(1-x_K^2\big)=1-\frac{1}{2\Delta^2}\quad(\Delta\ge1) \]

olur. \(\Delta=1\)'de \(W_b=1/2\)'dir: çiftin yarısı birlikte kalır, yarısı dağılır. \(\Delta<1\)'de integral yalnız bağlı durumun var olduğu \(|K|>2\arccos\Delta\) bölgesi üzerinden alınır.

Eleştirel notlar

Magnonlar ne ölçüde bozondur?

Lineer spin dalgası teorisi neyi kaçırır?

Lineer teoride magnonlar serbest bozonlardır ve iki-magnon spektrumu yalnızca süreklilikten oluşur. Bağlı durum, Holstein–Primakoff açılımındaki dördüncü dereceden (magnon–magnon) terimlerden gelir. Büyük \(S\)'de bu terimler \(1/S\) mertebesinde küçülür, ama bağlı durum ortadan kalkmaz (Wortis, 1963).

Δ = 0: etkileşim nerede?

\(\Delta=0\)'da yalnız sert-çekirdek kalır. 1D'de sert-çekirdekli bozonlar Jordan–Wigner dönüşümüyle serbest fermiyonlara eşlenir. Bu yüzden 1. adımda \(\Delta=0\) iken, magnonlar karşılaşıp birbirinin içinden “geçtikten” sonra bile site yoğunluğu etkileşimsiz öngörüyle örtüşür. Sert-çekirdeğin izi iki-parçacık korelasyonlarında (\(P(r)\), harita) görülür. Halkada ince bir istisna var: Jordan–Wigner dönüşümü iki fermiyon için antiperiyodik sınır koşulu üretir. Bu yüzden dalga cepheleri halkayı dolaştığında (\(t\approx N/2J\) civarı) yoğunlukta da küçük bir fark belirir. Bu bir topoloji etkisidir; sonsuz zincirde fark hiç oluşmaz.

Boyut ve anizotropi

1D Heisenberg zincirinde bağlı durum her \(K\neq0\) için vardır. 2D ve 3D'de izotropik durumda yalnızca Brillouin bölgesinin sınırına yakın bir kısımda bulunur. Kolay-eksen anizotropisi (\(\Delta>1\)) bağlanmayı güçlendirir; kolay-düzlem (\(\Delta<1\)) zayıflatır.

Δ < 1 ve taban durumu

\(\Delta<1\) ve \(h

Sonlu halka

Spektrum kesiklidir ve bağlı durum sonlu boyuta üstel küçüklükte duyarlıdır. \(x\to1\) (\(K\to0\) ya da \(\Delta\to\cos(K/2)\)) yakınında bağlı durumun boyu \(N\)'yi aşar ve süreklilikten ayırt edilemez. Komşu çiftin uzaklaşan kısmı halkanın öbür yanından \(t_s\approx N/2J\) civarında geri döner.

Bethe ansatz

Bu simülasyondaki problem, Bethe'nin 1931'de Heisenberg zincirini tam çözerken attığı ilk adımdır. İki-magnon saçılma faz kayması \(\delta_q\) ve bağlı durumlar (karmaşık momentumlu “string” çözümleri), çok-magnonlu tam çözümün yapı taşlarıdır.

Önerilen alıştırmalar

Ne zaman etkileşirler? Adım 1, \(\Delta=1\), \(N=48\). Şekil 2'deki sapma \(D(t)\) hangi anda büyümeye başlıyor? Bunu iki cephenin karşılaşma zamanı \(t\approx N/4J\) ile karşılaştırın.
Jordan–Wigner. Aynı deneyi \(\Delta=0\) ile tekrarlayın. Karşılaşmadan sonra site yoğunluğu neden etkileşimsiz öngörüden sapmıyor? Haritaya ve \(P(r)\)'ye bakın: fark nerede? Daha uzun beklediğinizde ortaya çıkan küçük sapmayı halkadaki antiperiyodik sınır koşuluyla açıklayın.
Yarısı kalır. Adım 2, \(\Delta=1\). \(P(r\le3)\)'ün asimptota yaklaştığını ve bağlı durum ağırlığının \(W_b\simeq1/2\) olduğunu gösterin. \(\Delta=2\) için \(W_b=1-1/(2\Delta^2)=0.875\)'i doğrulayın.
Bethe bağlı durumu. Adım 3, \(\Delta=1\), \(h=0\). \(K=\pi\)'de bağlı durum enerjisinin \(J\), dalga fonksiyonunun \(f(r)=\delta_{r,1}\) olduğunu gösterin. Neden \(K=\pi\)'de göreli atlama sıfırdır?
Bağlı durumun kaybolması. \(\Delta=0.5\)'te bağlı durum yalnız \(|K|>2\arccos\Delta=2\pi/3\) için vardır. Spektrumda bunu bulun ve sınırda dalga fonksiyonunun nasıl yayıldığına bakın.
Ağır parçacık. Adım 4, \(K_0=\pi/2\). Bağlı paketin hızını \(J\sin K_0\) (tek magnon) ile karşılaştırın. \(\Delta=2\)'de hız nasıl değişir? Ising limitinde neden bağlı çift neredeyse hareketsiz kalır?

Okuma kılavuzu

Antiferromanyet simülasyonu nasıl okunur?

Çubuklar (1. adım)

\(\langle S^z_i(t)\rangle\): yukarı (camgöbeği) ve aşağı (turuncu). Kesikli çizgiler Néel değerleri \(\pm S\)'dir. Düzen parametresi \(m_s=\frac{1}{NS}\sum_i(-1)^i\langle S^z_i\rangle\), Néel'de 1, taban durumunda 0'dır.

Konfigürasyon listesi (2. adım)

Her satır \(S^z\) tabanından bir durumdur (↑, ↓; \(S=1\)'de ○ \(m=0\)). Sağdaki çubuk, kesin taban durumunda o konfigürasyonun olasılığı \(|\langle c|\psi_0\rangle|^2\)'dir. Ötelemeyle birbirine dönüşen konfigürasyonlar aynı olasılığı taşır.

Bogoliubov eğrileri (3. adım)

Üstte lineer teorinin magnon enerjisi \(\varepsilon_{\mathbf k}\), altta taban durumundaki bozon sayısı \(v_{\mathbf k}^2\). \(v_{\mathbf k}^2\), \(\varepsilon_{\mathbf k}\to0\) olan noktalarda (\(\Gamma\) ve \(M\)/\(k=\pi\)) ıraksar. Bu ıraksamanın integrale katkısı boyuta bağlıdır.

Spektrum haritası (4. adım)

Yatay eksen momentum transferi \(q\), dikey eksen enerji \(\omega\). Renk \(S^{zz}(q,\omega)\) şiddetidir (karekök ölçeği). Sonlu zincirde her sütun kesikli kutuplardan oluşur. Çizgiler: lineer teori (kesikli), \(S=\tfrac12\) için iki-spinon sürekliliğinin alt ve üst sınırları.

Türetme

Néel durumundan spinonlara

Adım 1

Néel durumu neden erir?

\(\mathbf S_i\cdot\mathbf S_{i+1}=S^z_iS^z_{i+1}+\tfrac12(S^+_iS^-_{i+1}+S^-_iS^+_{i+1})\). FM durumunda çevirme terimleri sıfır verir. Néel durumunda ise her ↑↓ bağı çevrilebilir:

\[ H|{\rm N}\rangle=-NJS^2|{\rm N}\rangle+\tfrac J2\sum_i\big(S^+_iS^-_{i+1}+S^-_iS^+_{i+1}\big)|{\rm N}\rangle . \]

İkinci terim sıfır olmadığı için \(|{\rm N}\rangle\) özdurum değildir ve zamanla değişir. Enerjisi korunur (\(-JS^2\) bağ başına), ama \(zz\) korelasyonundan \(xy\) korelasyonlarına aktarılır. Néel durumu ayrıca toplam spinin özdurumu da değildir. Kesin taban durumu ise bir singlettir (\(S_{\rm top}=0\), Lieb–Mattis).

Adım 2

Kesin taban durumu

Néel durumu bir varyasyonel deneme durumudur ve gerçek enerjinin üst sınırını verir: \(E_{\rm N}/N=-JS^2\). Kesin sonuç \(S=\tfrac12\) için Bethe ansatz (Hulthén, 1938) ile bilinir:

\[ \frac{E_0}{NJ}=\tfrac14-\ln2\simeq-0.4431\quad\text{(Néel: }-0.25\text{)} . \]

Kuantum dalgalanmaları enerjiyi neredeyse iki kat düşürür. Korelasyonlar \(S=\tfrac12\)'de \(\langle\mathbf S_0\cdot\mathbf S_r\rangle\sim(-1)^r\sqrt{\ln r}/r\) ile, \(S=1\)'de üstel olarak (\(\xi\approx6\)) azalır. 1D'de uzun menzilli Néel düzeni yoktur. Néel durumları en olası konfigürasyonlar olsa da ağırlıkları \(N\) ile üstel olarak küçülür.

Adım 3

Holstein–Primakoff ve Bogoliubov

A alt örgüsünde \(S^z=S-a^\dagger a\), B'de \(S^z=-S+b^\dagger b\). İkinci derecede

\[ H_2=JzS\sum_{\mathbf k}\Big[a^\dagger_{\mathbf k}a_{\mathbf k}+b^\dagger_{\mathbf k}b_{\mathbf k}+\gamma_{\mathbf k}\big(a_{\mathbf k}b_{-\mathbf k}+{\rm h.c.}\big)\Big], \]

\(\gamma_{\mathbf k}=\frac1z\sum_\delta e^{i\mathbf k\cdot\delta}\). \(ab\) terimi bozon sayısını korumaz. Bu yüzden \(\alpha_{\mathbf k}=u_{\mathbf k}a_{\mathbf k}-v_{\mathbf k}b^\dagger_{-\mathbf k}\) dönüşümü gerekir:

\[ \varepsilon_{\mathbf k}=JzS\sqrt{1-\gamma_{\mathbf k}^2},\qquad v_{\mathbf k}^2=\tfrac12\Big(\tfrac{1}{\sqrt{1-\gamma_{\mathbf k}^2}}-1\Big). \]

Yeni vakumda \(\langle a^\dagger a\rangle=\Delta S=\overline{v_{\mathbf k}^2}\neq0\). 1D'de \(\Delta S\) logaritmik ıraksar, 2D'de \(0.197\), 3D'de \(0.078\) olur. \(S=\tfrac12\) kare örgüde \(m=0.303\) çıkar (QMC: \(0.307\)). Dağılım klasik lineer teoriyle aynıdır; fark sıfır-nokta dalgalanmalarıdır.

Adım 4

Spinonlar ve Haldane boşluğu

Nötron saçılmasının ölçtüğü büyüklük

\[ S^{zz}(q,\omega)=\sum_n|\langle n|S^z_q|0\rangle|^2\,\delta(\omega-E_n+E_0) \]

dinamik yapı faktörüdür. LSWT tek bir keskin dal verir: \(\omega=2JS|\sin q|\). \(S=\tfrac12\) zincirde ise bir spin çevirmesi (\(\Delta S=1\)) iki \(S=\tfrac12\) spinona ayrışır ve ağırlık bir sürekliliğe dağılır:

\[ \tfrac\pi2 J|\sin q|\;\le\;\omega\;\le\;\pi J|\sin\tfrac q2| . \]

Alt sınır des Cloizeaux–Pearson sonucudur. LSWT bu sınırı (\(J|\sin q|\)) \(\pi/2\) kadar küçük tahmin eder. \(S=1\)'de spektrum boşlukludur (\(\Delta\simeq0.4105J\), \(q=\pi\)) ve \(q=\pi\) civarında iyi tanımlı bir üçlü mod vardır: “magnon” geri döner, ama kütlelidir (Haldane, 1983). 2D kare örgüde magnonlar yaşar; LSWT hızı yalnızca yaklaşık %18 düzeltilir.

Eleştirel notlar

Neyi gösteriyor, neyi göstermiyor?

Sonlu boyut

Tam köşegenleştirme \(N\le16\) (\(S=\tfrac12\)) ve \(N\le10\) (\(S=1\)) ile sınırlıdır. Spektrum kesiklidir ve \(S=\tfrac12\)'de \(q=\pi\)'deki en düşük uyarılma \(\sim1/N\) ile sıfıra gider. Spinon sürekliliğini sonlu zincirde tek tek kutupların sürekliliğin sınırları arasına yayılması olarak görürüz. Haldane boşluğu \(N=10\)'da henüz tam olarak yakınsamaz (değeri biraz büyük çıkar).

Zaman evrimi

1. adım bir “kuantum söndürme” deneyidir: Néel durumu taban durumu değildir, bu yüzden sistem dengeye gitmez, enerjisi sabit kalır. Sonlu halkada yarı-periyodik canlanmalar görülür. Uzun zaman davranışı (termalleşme, integrabilite) bu simülasyonun kapsamı dışındadır.

LSWT ne zaman güvenilir?

Lineer teori \(1/S\) açılımının ilk terimidir ve \(\Delta S\ll S\) olduğunda tutarlıdır. 2D ve 3D'de \(S=\tfrac12\) için bile niceliksel olarak iyidir. 1D'de ise \(\Delta S\) ıraksadığı için niteliksel olarak yanlıştır: uzun menzilli düzen ve keskin magnonlar yoktur. Yarım tam sayılı spinlerde zincir boşluksuz kalır (Lieb–Schultz–Mattis), tam sayılılarda boşluklu olur (Haldane).

Deneysel bağlantı

Yarı-1D \(S=\tfrac12\) antiferromanyetlerde (ör. KCuF\(_3\)) nötron saçılması iki-spinon sürekliliğini doğrudan göstermiştir. \(S=1\) Ni-zincir bileşiklerinde Haldane boşluğu ölçülmüştür. Kare örgü \(S=\tfrac12\) bileşiklerinde ise (ör. La\(_2\)CuO\(_4\)) spin dalgası teorisiyle uyumlu magnonlar gözlenir.

Önerilen alıştırmalar

FM ile karşılaştırma. Adım 1. Néel durumunun neden zamanla değiştiğini, FM durumunun neden değişmediğini \(H\)'nin iki durumdaki etkisiyle gösterin. \(zz\) ve \(xy\) bağ enerjilerinin toplamının neden sabit kaldığını açıklayın.
Varyasyonel sınır. Adım 2, \(S=\tfrac12\). \(N=8,12,16\) için \(E_0/N\)'yi okuyun ve \(\tfrac14-\ln2\)'ye yakınsamasını Néel (\(-0.25\)) ve LSWT (\(-0.4317\)) değerleriyle karşılaştırın.
Néel ağırlığı. İki Néel durumunun taban durumundaki toplam ağırlığını \(N\)'nin fonksiyonu olarak not edin. Neden üstel olarak azalır?
Boyut ve sıfır-nokta. Adım 3. 1D'de örgü boyutunu artırın ve \(\Delta S\)'nin \(\ln N\) ile büyüdüğünü gösterin. 2D'de \(\Delta S\)'nin yakınsadığını gösterin. \(S=\tfrac12\) kare örgü için \(m=S-\Delta S\)'yi QMC değeri \(0.307\) ile karşılaştırın.
Spinon sürekliliği. Adım 4, \(S=\tfrac12\), \(N=16\). \(q=\pi/2\)'de ağırlığın tek bir kutupta toplanmadığını ve en düşük kutbun LSWT'den çok des Cloizeaux–Pearson eğrisine yakın olduğunu gösterin.
Haldane boşluğu. \(S=1\)'e geçin. \(q=\pi\)'de en düşük uyarılmanın sıfırdan uzak kaldığını ve ağırlığın büyük kısmının tek bir kutupta toplandığını gösterin. \(N=6,8,10\) için boşluğu \(0.4105J\) ile karşılaştırın.

Yol haritası

Simülasyonlar

Simülasyon	İçerik	Yöntem	Durum
Tek magnon	FM taban durumu, çevrik spin, \(\|k\rangle\), dalga paketi, klasik limit	Tek-magnon sektörü, kesin	hazır
İki magnon	Etkileşimin kaynağı, saçılma sürekliliği, Bethe bağlı durumu, XXZ	İki-magnon sektörü, kesin	hazır
Antiferromanyet	Néel durumunun erimesi, kesin taban durumu, Holstein–Primakoff/Bogoliubov, spinonlar ve Haldane boşluğu	Tam köşegenleştirme (Lanczos, Krylov) ve LSWT	hazır
Termal magnonlar	Bose istatistiği, Bloch \(T^{3/2}\), klasik eşbölüşümle karşılaştırma	Analitik ve sayısal toplam	planlı

Magnon Laboratuvarı · Tüm hesaplar tarayıcıda yapılır · Adres satırına ?sim=tek, ?sim=iki ya da ?sim=afm eklenerek simülasyon, #simulasyon eklenerek tam sayfa görünüm doğrudan açılır. Klasik karşılık: Spin Dalgası Laboratuvarı.