6.4. Magneetilised korrastused

6.4.1. Ferromagnetism keskmise välja lähenduses

Eksisteerivad erinevad magneetilised korrastused. Pildil tähistavad nooled magneetiku kristallvõre sõlmedel lokaliseeritud magnetmomente.

../_images/magnetism.png
  • Ferromagneetikus on need magnetmomendid allpool teatud temperatuuri orienteeritud kõik ühes suunas

  • Kui ferromagneetikus on kaks alamvõret, siis nende alamvõrede sõlmedel lokaliseeritud magnetmomendid orienteeruvad jälle ühes suunas, kuid nende magnetmomentide moodulid on erinevad. Kahel esimesel juhul tekkib nullist erinev makroskoopiline magnetmoment ehk magneetumus.

  • Antiferromagneetikus on lokaliseeritud magnetmomendid orienteeritud antiparalleelselt. See tähendab, et makroskoopilist nullist erinevat magnetmomenti ei eksisteeri, aga mikroskoopilisel tasandil on siiski teatud magneetiline korrastus olemas.

  • Ferrimagneetikus on antiparalleelselt orienteeritud lokaalsed magnetmomentide moodulid erinevad

  • Saavad olla ka magneetiliste korrastuste keerukamad juhud

Olgu \(\hat{\mathbf{S}}_i\) on võre sõlmes nr \(i\) paikneva iooni osaliselt täidetud elektronkihi summaarse spinni operaator. Kasutame magneetilise korrastuse kirjeldamiseks hästi tuntud Heisenbergi mudeli, mis kasutab välises magnetväljas oleva spinnsüsteemi kirjeldamiseks hamiltoniaani

\[\hat{H}=-\sum_{i,j}'J_{ij}\hat{\mathbf{S}}_{i}\hat{\mathbf{S}}_{j} -\mu_\mathrm{B}g\mathbf{H}\sum_{i}\hat{\mathbf{S}}_{i}.\]

Selle operaatori esimene panus kujutab endast kristallvõre sõlmedes nr \(i\) ja \(j\) paiknevate spinnide interaktsiooni. See interaktsioon põhjustabki magneetilist korrastust. Teine liige arvestab välise magnetvälja mõjuga spinnidele. Lisaks

  • Summa \(\sum_{i,j}'\) tähistab, et \(i\neq j\) summeerimisel

  • \(J_{ij}=J_{ji}\) on spinn-spinn interaktsiooni konstant

  • \(\mu_\mathrm{B}\) on Bohri magneton

  • \(\mathbf{H}\) on välise magnetvälja tugevus

  • \(g<0\) on konstant, mis arvestab, millise spinniga on tegemist

Heisenbergi mudel ei ole täpselt lahendatav. See tähendab, et ei saa näiteks täpselt leida makroskoopilist magneetumust, mis on tekitatud spinn-spinn interaktsiooni poolt. Kõige levinum ligikaudne meetod antud hamiltoniaaniga kirjeldatud spinnsüsteemi omaduste uurimiseks on keskmise välja lähendus (tuntud ka kui molekulaarvälja lähendus). Selles skeemis asendatakse spinnidevaheline interaktsioon spinnile mõjuva teiste spinnide poolt tekitatud efektiivse (keskmise) väljaga

\[\hat{H}\rightarrow \hat{H}_{\mathrm{mf}}=-\sum_{i}\left(\mu_\mathrm{B}g\mathbf{H} +2\sum_{j}'J_{ij}\langle\hat{\mathbf{S}}_{j}\rangle\right)\hat{\mathbf{S}}_{i} =-\mu_\mathrm{B}g\sum_{i}\left(\mathbf{H}+\mathbf{H}_{\mathrm{eff}}^{(i)}\right)\hat{\mathbf{S}}_{i} ,\]

kus

\[\mathbf{H}_{\mathrm{eff}}^{(i)}=\frac{2}{\mu_\mathrm{B}g}\sum_{j}'J_{ij}\langle\hat{\mathbf{S}}_{j}\rangle\]

on efektiivne magnetväli. Statistiline keskväärtus on defineeritud standartselt statistilise operaatori kaudu

\[\langle\hat{\mathbf{S}}_{j}\rangle =\mathrm{Tr}\{\hat{\mathbf{S}}_{j}\hat\rho(\hat{H}_\mathrm{mf})\},\qquad\hat\rho(\hat{H}_\mathrm{mf})= Z^{-1}\exp\left(-\frac{\hat{H}_{\mathrm{mf}}}{k_\mathrm{B}T}\right),\qquad Z=\mathrm{Tr} \exp\left(-\frac{\hat{H}_{\mathrm{mf}}}{k_\mathrm{B}T}\right).\]

Teeme veel järgmised lihtsustused

  • Me võtame arvesse interaktsiooni ainult antud võre sõlme lähimate naabritega: \(J_{ij}=J\neq 0\), kui \(i,j\) on lähimad naabrid, ja \(J_{ij}=0\), kui on teisiti.

  • Keskväärtuse \(\langle\hat{\mathbf{S}}_{j}\rangle\) loeme sõltumatuks võre sõlmest, \(\langle\hat{\mathbf{S}}_{j}\rangle\to \langle\hat{\mathbf{S}}\rangle\).

Sellisel juhul

\[\mathbf{H}_{\mathrm{eff}}=\frac{2n_0J}{\mu_\mathrm{B}g}\langle\hat{\mathbf{S}}\rangle,\qquad \hat{H}_{\mathrm{mf}}=-\mu_\mathrm{B}g\left(\mathbf{H}+\mathbf{H}_{\mathrm{eff}}\right)\sum_{i=1}^N\hat{\mathbf{S}}_{i}\]

kus \(n_0\) on lähimate naabrite arv ja \(N\) on võre sõlmede arv. Siin on kõik magnetvälja suunad võrdväärsed ja saab valida \(\mathbf{H}=(0,0,H)\) ja \(\mathbf{H}_\mathrm{eff}=(0,0,H_\mathrm{eff})\). Sellisel juhul võrdub keskmise välja hamiltoniaan

(6.8)\[\hat{H}_{\mathrm{mf}}=-\mu_\mathrm{B}g\left(H+H_{\mathrm{eff}}\right)\sum_{i=1}^N\hat{S}_{iz}.\]

Leitud Hamiltoni operaatori omavektoriteks on \(|S_{1z},S_{2z}\ldots S_{Nz}\rangle=\Pi_{i=1}^N|S_{iz}\rangle\), kus \(|S_{iz}\rangle\) on operaatori \(\hat{S}_{iz}\) ortonormeeritud omavektorid, see tähendab \(\hat{S}_{iz}|S_{iz}\rangle=S_{iz}|S_{iz}\rangle\), kus \(S_{iz}\) on vastav omaväärtus. Tänu ortonormeerimisele kehtib \(\langle S_{iz}|S_{iz}'\rangle=\delta_{S_{iz},S_{iz}'}\).

Teeme ülesande konkreetsemaks. Olgu meil tegemist spinniga \(S=1/2\), see tähendab, et võre sõlmel lokaliseeritud spinn võrdub ühe elektroni spinniga. Sellisel juhul on operaatori \(\hat{S}_z\) omaväärtuseks \(S_z=\pm1/2\) ning \(g=-2.002\). Üldjuhul, kui on tegemist spinniga \(1/2\), siis spinnoperaatoril on kuju \(\hat{\mathbf{S}}=(\frac{1}{2}\hat{\sigma}_x,\frac{1}{2}\hat{\sigma}_y,\frac{1}{2}\hat{\sigma}_z)\), kus \(\hat{\sigma}_{x,y,z}\) on Pauli maatriksid

\[\begin{split}\hat{\sigma}_x=\left(\begin{array}{cc}0&1\\1&0\end{array}\right),\quad \hat{\sigma}_y=\left(\begin{array}{cc}0&-i\\i&0\end{array}\right),\quad \hat{\sigma}_z=\left(\begin{array}{cc}1&0\\0&-1\end{array}\right).\end{split}\]

Operaatori \(\hat{S}_z\) omavektoriteks, mis vastavad omaväärtustele \(S_z=\pm1/2\), on

\[\begin{split}|+1/2\rangle=\left(\begin{array}{c}1\\0\end{array}\right),\qquad |-1/2\rangle=\left(\begin{array}{c}0\\1\end{array}\right),\end{split}\]

nii, et \(\hat{S}_z|\pm1/2\rangle=\frac{1}{2}\hat{\sigma}_z|\pm1/2\rangle=\pm\frac{1}{2}|\pm1/2\rangle\).

Valitud baasil võrdub tihedusmaatriks

\[\begin{split}&\langle S_{1z},S_{2z}\ldots S_{Nz}|\hat{\rho}(\hat{H}_\mathrm{mf})|S_{1z}',S_{2z}'\ldots S_{Nz}'\rangle=\\ &Z^{-1}\langle S_{1z},S_{2z}\ldots S_{Nz}|\exp\left(-\frac{\hat{H}_{\mathrm{mf}}}{k_\mathrm{B}T}\right)|S_{1z}',S_{2z}'\ldots S_{Nz}'\rangle=\\ &Z^{-1}\prod\limits_{i=1}^N \exp\left(\frac{\mu_\mathrm{B}g(H+H_{\mathrm{eff}})S_{iz}'}{k_\mathrm{B}T}\right)\prod\limits_{j=1}^N\delta_{S_{jz},S_{jz}'}.\end{split}\]

Selle tulemuse saamiseks oli kasutatud (6.8), eksponendi \(\exp\left(-\frac{\hat{H}_{\mathrm{mf}}}{k_\mathrm{B}T}\right)\) rittaarendamine, omavektorite ortonormeerimine \(\langle S_{iz}|S_{iz}'\rangle=\delta_{S_{iz},S_{iz}'}\) ning saadud rea eksponendiks tagasi võtmine.

Meil on veel tarvis leida statistiline summa

\[\begin{split}Z=&\mathrm{Tr}\langle S_{1z},S_{2z}\ldots S_{Nz}|\exp\left(-\frac{\hat{H}_{\mathrm{mf}}}{k_\mathrm{B}T}\right)|S_{1z}',S_{2z}'\ldots S_{Nz}'\rangle=\sum_{S_{1z}}\ldots\sum_{S_{Nz}} \prod\limits_{i=1}^N \exp\left(\frac{\mu_\mathrm{B}g(H+H_{\mathrm{eff}})S_{iz}}{k_\mathrm{B}T}\right)=\\ &\prod\limits_{i=1}^N\sum_{S_{iz}} \exp\left(\frac{\mu_\mathrm{B}g(H+H_{\mathrm{eff}})S_{iz}}{k_\mathrm{B}T}\right)= \prod\limits_{i=1}^N\left[ \exp\left(\frac{\mu_\mathrm{B}g(H+H_{\mathrm{eff}})}{2k_\mathrm{B}T}\right)+\exp\left(-\frac{\mu_\mathrm{B}g(H+H_{\mathrm{eff}})}{2k_\mathrm{B}T}\right) \right]=\\ &\prod\limits_{i=1}^N 2\cosh\left(\frac{\mu_\mathrm{B}g(H+H_{\mathrm{eff}})}{2k_\mathrm{B}T}\right)=\left[2\cosh\left(\frac{\mu_\mathrm{B}g(H+H_{\mathrm{eff}})}{2k_\mathrm{B}T}\right)\right]^N\end{split}\]

Kasutades neid tulemusi, leiame spinni keskväärtus

(6.9)\[\langle \hat{S}_{iz}\rangle=\mathrm{Tr}\left\{\hat{S}_{iz}\hat{\rho}(\hat{H}_\mathrm{mf})\right\}=\frac{1}{2}\tanh\left(\frac{\mu_\mathrm{B}g(H+H_{\mathrm{eff}})}{2k_\mathrm{B}T}\right)=\langle \hat{S}_{z}\rangle.\]

Tuletuskäik on siin eelnevaga analoogne.

Magneetumus ehk makroskoopilise magnetmomendi tihedus on defineeritud, kui \(M=\frac{N}{V}\mu_\mathrm{B}g\langle \hat{S}_{z}\rangle\). Keskmise välja lähenduses on magneetumus võrdeline efektiivse magnetväljaga

\[H_\mathrm{eff}=\frac{2n_0JV}{(\mu_\mathrm{B}g)^2N}M.\]

Avaldis (6.9) defineerib meil ka magneetumuse, kuna \(\langle \hat{S}_{z}\rangle\sim M\). Samal ajal seisab avaldises (6.9) magneetumus eksponendis, sest \(H_\mathrm{eff}\sim M\). Seega tekkib meil võrrand magneetumuse määramiseks

\[M=\lambda\tanh\left(\frac{\mu_\mathrm{B}|g|H+\Lambda M)}{2k_\mathrm{B}T}\right),\qquad \lambda=\frac{1}{2}\frac{N}{V}\mu_{\mathrm{B}}|g|,\quad \Lambda = \frac{2n_{0}JV}{\mu_{\mathrm{B}}|g|N}.\]

Juhul, kui väline magnetväli puudub \(H=0\), saadakse magneetumuse jaoks võrrand

(6.10)\[M=\lambda \tanh \left(\frac{\Lambda M}{2k_{\mathrm{B}}T}\right).\]

Sellel võrrandil on olemas nullist erinev reaalarvuline lahend \(M\neq 0\), kui \(J\sim\Lambda>0\) ja \(T<T_c\), kus faasisiirde temperatuuri

\[T_c=\frac{\lambda\Lambda}{2k_\mathrm{B}}=\frac{Jn_{0}}{2k_{\mathrm{B}}}\]

saab leida piirprotsessis \(M\to0\). II liiki ferromagneetilise faasisiirde temperatuuri \(T_c\) nimetatakse ka Curie-Weissi temperatuuriks.

Sõltuvust \(M=M(T)\) saab leida ainult numbriliselt ja see on tüüpiline korrastusparameetri jaoks II liiki faasisiirde puhul

  • kui \(T>T_c\), siis oleme paramagneetilises faasis, kus \(M=0\)

  • kui \(T<T_c\), siis realiseerub ferromagneetiline faas ja \(M\neq 0\)

  • faasisiirde punkti lähedal \(T\to T_c-\) kehtib \(M\sim\sqrt{T_c-T}\)

6.4.2. Antiferromagnetism keskmise välja lähenduses

Olgu meil jälle olukord, kui \(S=1/2\), ja puudub väline magnetväli. Ferromagneetiku magneetumus \(M\) rahuldab sel juhul keskmise välja lähenduses võrrandit (6.10). Üldistame seda võrrandit kahe alamvõrega süsteemi jaoks. Olgu \(M_{a}\) ja \(M_{b}\) alamvõrede \(a\) ja \(b\) magneetumused. Arvestame ainult erinevatel alamvõredel lokaliseeritud spinnide vahelist interaktsiooni ja seda omakorda ainult lähimate naabrite vahel. Alamvõredevaheliseks spinn-spinn interaktsioonikonstandiks on \(J_{ab}=J_{ba}\). Alamvõrede magneetumused rahuldavad siis võrrandisüsteemi

(6.11)\[M_{a}=\lambda \tanh \left(\frac{\Lambda_{ab} M_{b}}{2k_{\mathrm{B}}T}\right),\qquad M_{b}=\lambda \tanh \left(\frac{\Lambda_{ba} M_{a}}{2k_{\mathrm{B}}T}\right),\]

kus

\[\lambda = \frac{1}{2}\frac{N}{V}\mu_{\mathrm{B}}|g|,\qquad \Lambda_{ab} = J_{ab} \frac{2n_{0}V}{\mu_{\mathrm{B}}|g|N}\]

ja \(N\) on sõlmede arv ühel alamvõrel: \(N_{a}=N_{b}\equiv N\).

Võrranditest (6.11) järeldub, et

  • alamvõrede magneetumused on samamärgilised \(\mathrm{sgn}(M_a)=\mathrm{sgn}(M_b)\), kui \(J_{ab}>0\). Tegemist on ferromagnetismiga.

  • kui \(J_{ab}<0\), siis alamvõrede magneetumused on vastasmärgilised \(\mathrm{sgn}(M_a)\neq\mathrm{sgn}(M_b)\). See uus lahend vastab antiferromagnetismile.

Tehtud eelduste korral saame võrranditest (6.11), et \(|M_a|=|M_b|\). Antiferromagnetilisele korrastusele vastav lahend \(M_a=-M_b\) (\(J_{ab}<0\)) realiseerub, kui \(T<T_{\mathrm{N}}\), kus Neeli temperatuuri

\[T_{\mathrm{N}}=-\frac{J_{ab}n_{0}}{2k_{\mathrm{B}}}=\frac{|J_{ab}|n_{0}}{2k_{\mathrm{B}}}\]

saab leida, kui minna piirile \(M_{a,b}\to0\). Temperatuur \(T_\mathrm{N}\) on antiferromagneetilise faasisiirde temperatuur. Seega, kui \(T<T_{\mathrm{N}}\), siis realiseerub antiferromagneetiline faas.