スピントロニクス理論の基礎/5-6 の変更点

更新 
[[[前の章へ]>スピントロニクス理論の基礎/5-5]] <<<<　
[[スピントロニクス理論の基礎]]（目次）　>>>>
[[[次の章へ]>スピントロニクス理論の基礎/5-7]]

#contents

* 5-6 磁壁粒子 [#jb28ce63]

教科書だとこれまでの章との関係が分かりづらいのでちょっと整理してみる。

まず、この章では磁壁があたかも粒子のように振る舞う条件を考える。~
ただし、「粒子のように振る舞う」とは、質量を持っていて慣性が見られるような運動をすることを指す。

- 外場を加えたら徐々に速度が変化する
- 運動エネルギーを定義できて、ポテンシャルとのやりとりで調和振動が生じる

などが「粒子のような振る舞い」の代表例。

これまで、スピンダイナミクスに取り入れる相互作用として

+ スピン・スピン相互作用 &math(J=\lambda K);
+ 容易軸相互作用 &math(K);
+ 外部磁場 &math(B);
+ 困難軸相互作用 &math(K_\perp);
+ ピン止め相互作用 &math(V_0);

の４つを考え、4-5, 5-3〜5-4 において４つのケースについて見てきた。

** 4-5 スピン／スピン相互作用＋容易軸 [#p1cc55f1]

１次元磁壁解が存在する。

もちろんその他の解も多数存在する。

磁壁解の境界条件は、系の両端で容易軸方向かつ逆向きの磁化を持つこと。

** 5-3 困難軸相互作用のない場合 (5.33) [#k333c792]

1 + 2 + 3 だけで 4, 5 が重要でない場合の結果が (5.33) で、その性質は

- 「磁場に比例」した速度で等速並進運動し、回転も等速になる
- 摩擦が大きくなればなるほど並進速度は上がり、回転速度は下がる

であった。

磁場を掛けた瞬間に等速運動が開始するような運動で、
「粒子的」とはほど遠いと言える。

** 5-5 困難軸があり、磁場が弱いとき (5.42) [#ybd7fca8]

ここではピン止めは考えなかった。

この運動の性質は、

- 回転は止まっている
- 磁場に比例した速度で等速並進運動する解がある
- 摩擦が大きくなるほど並進速度は低下する

であり、(5.33) との大きな違いは、
「等速並進運動」は解の１つであり、それ以外の解も存在することである。

特にこのケースは、&math(B\ll K_\perp); の条件では後に &math(X); 
粒子と呼ばれるように磁壁があたかも粒子として運動し、
(5.42) で見た等速並進運動は外部磁場による加速運動の極限として理解されることになる。

** 5-5 困難軸があり、磁場が強いとき (5.47) [#qe18e8bf]

ここでもピン止めは考えなかった。

この運動の性質は、

- &math(\phi_0); は振動的
- &math(X); は振動しつつ並進する
- 並進運動、回転運動とも、磁場が強くなると遅くなる
- 磁場の強い極限では、困難軸異方性は無視できて (5.33) に帰着する
- すなわち並進運動の速度は、
磁場がそれほど強くなければ摩擦が大きくなると遅くなるが、
磁場が強いところでは摩擦が大きい方が速くなる

であった。

この場合も、これ以外の解は存在せず、
「粒子的」とはほど遠い運動と言える。

* 「粒子」の運動方程式 [#r275c74c]

質量を持った「粒子」の運動方程式は、ニュートンの方程式は

&math(
F=M\ddot X
);

であるが、特に力 &math(F); が保存力であり、ポテンシャル &math(U); が定義できる時は、

&math(F=-\frac{\PD U}{\PD X});

となる。

ハミルトニアンは、運動エネルギーとポテンシャルエネルギーの和で書けるが、

&math(H=\frac{1}{2M}P^2+U(X));

その正準方程式は、

&math(
&\dot X=\frac{\PD H}{\PD P}=\frac{P}{M}\\
&\dot P=-\frac{\PD H}{\PD X}=-\frac{\PD U}{\PD X}
);

であり、この第１式を微分した上で &math(\dot P); に第２式を代入すれば、

&math(
\ddot X=\frac{\dot P}{M} =-\frac{1}{M}\frac{\PD U}{\PD X}=\frac{M}{F}
);

のように、上記運動方程式が得られることを確認できる。

** 緩和項がある場合 [#v915c277]

緩和項がある場合には、

&math(
F-\frac{M}{\tau}\dot X=M\ddot X
);

で表される。ここに現れた &math(\tau); は緩和時定数で、
外力 &math(F=0); の時に

&math(
M\ddot X=-\frac{M}{\tau}\dot X\\
\ddot X=-\frac{1}{\tau}\dot X\\
\dot X=\dot X_0e^{-\frac{t}{\tau}}
);

のようにして、速度の減衰時間を表す定数である。

** 時間についての２次微分方程式 [#va89131e]

５章に現れた「粒子的でない場合」の運動方程式と比べると、粒子的な方程式は
&math(\ddot X); として時間に対する２次微分成分が現れているところが大きく異なる。

この２次微分の項は、&math(X); と &math(P); の絡み合う連立１次微分方程式から、
&math(\dot P); を消去するために片方の方程式を微分した際に現れたことに注意せよ。

* (5.40), (5.41) が粒子的になる２つの条件 [#sbf1ca05]

(5.40)

&math(
\dot \phi_0+\alpha\frac{\dot X}{\lambda}=\gamma B-\frac{2\lambda V_0}{\hbar S\xi^2}X\theta(\xi-|X|)
);

(5.41)

&math(
\dot X-\alpha\lambda\dot \phi_0=\frac{K_\perp\lambda}{2\hbar}S\sin 2\phi_0
);

を見ると、&math(\theta(\xi-|X|)); や &math(\sin 2\phi_0); があることが、
&math(\dot X); や &math(\dot \phi_0); を消去する妨げになっている。

そこで、これらを無視できる条件、すなわち、

- &math(\phi_0\ll 1); により &math(\sin 2\phi_0\sim 2\phi_0); と見なせる条件 = &math(X); 粒子
- &math(|X|\ll\xi); により &math(\theta(\xi-|X|)\sim 1); と見なせる条件 = &math(\phi_0); 粒子

では、それぞれ &math(\dot\phi_0);、&math(\dot X); を消去して、
「粒子的な」運動方程式を得ることができる。

** X 粒子 (φ__0__<< 1) [#ucdda2c6]

(5.40) を元の &math(H_B); および &math(V_\mathrm{pin}); を使って書き、
(5.41) に代入して &math(\dot \phi_0); を消去することを試みる。

&math(
&\dot X-\alpha\lambda\dot\phi_0=\frac{K_\perp\lambda}{2\hbar}S\sin 2\phi_0\\
&\dot \phi_0=-\alpha\frac{\dot X}{\lambda}-\frac{\lambda}{N_w\hbar S}\frac{\PD}{\PD X}(H_B+V_\mathrm{pin})\\
&\dot X-\alpha\lambda\left\{
-\alpha\frac{\dot X}{\lambda}-\frac{\lambda}{N_w\hbar S}\frac{\PD}{\PD X}(H_B+V_\mathrm{pin})
\right\}=\frac{K_\perp\lambda}{2\hbar}S\sin 2\phi_0\\
&(1+\alpha^2)\dot X+\frac{\alpha\lambda^2}{N_w\hbar S}\frac{\PD}{\PD X}(H_B+V_\mathrm{pin})=\frac{K_\perp\lambda}{2\hbar}S\sin 2\phi_0\\
);

&math(\sin 2\phi_0); があるために &math(X); のみの微分方程式を得ることができない。

困難軸異方性があり、外力がなければ &math(\phi_0); は困難軸方向と直交する方向 (&math(\phi_0=0);) に固定されるが、外力があっても、困難軸異方性に比べて小さいときは、&math(\phi_0); の &math(0); からのずれは小さいことを仮定できる。

このとき、&math(\phi_0\ll 1); より

&math(\sin 2\phi_0\sim 2\phi_0);

と近似できるため、

&math(
&(1+\alpha^2)\dot X+\frac{\alpha\lambda^2}{N_w\hbar S}\frac{\PD}{\PD X}(H_B+V_\mathrm{pin})=\frac{K_\perp\lambda}{\hbar}S\phi_0\\
&(1+\alpha^2)\ddot X+\frac{\alpha\lambda^2}{N_w\hbar S}\frac{d}{dt}\frac{\PD}{\PD X}(H_B+V_\mathrm{pin})=\frac{K_\perp\lambda}{\hbar}S\dot\phi_0=\frac{K_\perp\lambda}{\hbar}S\left\{-\alpha\frac{\dot X}{\lambda}-\frac{\lambda}{N_w\hbar S}\frac{\PD}{\PD X}(H_B+V_\mathrm{pin})\right\}\\
&(1+\alpha^2)\frac{\hbar^2N_w}{K_\perp\lambda^2}\left[\ddot X
+\frac{\alpha}{1+\alpha^2}\frac{1}{\hbar S}\left\{
\frac{\lambda^2}{N_w}\frac{d}{dt}\frac{\PD V_\mathrm{pin}}{\PD X}
+K_\perp S^2\dot X
\right\}
+\frac{1}{1+\alpha^2}\frac{\lambda^2K_\perp}{N_w\hbar^2}\frac{\PD V_\mathrm{pin}}{\PD X}
\right]=-\frac{\PD H_B}{\PD X}-\frac{\alpha\hbar}{K_\perp S}\frac{d}{dt}\frac{\PD H_B}{\PD X}\\
&\red{(1+\alpha^2)}\frac{\hbar^2N_w}{K_\perp\lambda^2}\left[\ddot X
+\frac{\alpha}{\red{1+\alpha^2}}\frac{1}{\hbar\red S}\left\{
2V_0\frac{\lambda^2}{\xi^2}+K_\perp \red {S^2}\right\}\dot X
+\red{\frac{1}{1+\alpha^2}}\frac{2V_0K_\perp\lambda^2}{\hbar^2\xi^2}
\right]=\frac{N_w \hbar S\gamma B}{\lambda }
\,\red{+\frac{\alpha\hbar^2 N_w}{K_\perp\lambda}\gamma \dot B}\\
);

※ &math(\frac{\PD V_\mathrm{pin}}{\PD X}); を &math(\theta(\xi-|X|)); 
を無視して上記のように近似するのは &math(X); 
の大きなところでむちゃくちゃなことになるため注意が必要。

上記より、

(5.51)

&math(
&M_w\left[\ddot X+\frac{1}{\tau_w}\dot X+\Omega_\mathrm{pin}^2X\right]=F
);

(5.52)

&math(
M_w=\red{(1+\alpha^2)}\frac{\hbar^2N_w}{K_\perp\lambda^2}
);

&math(
\Omega_\mathrm{pin}=\sqrt{\red{\frac{1}{1+\alpha^2}}\frac{2V_0K_\perp\lambda^2}{\hbar^2\xi^2}}
);

&math(
&\frac{1}{\tau_w}=
\frac{\alpha}{\red{1+\alpha^2}}\frac{1}{\hbar \red S}\left\{
2V_0\frac{\lambda^2}{\xi^2}
+K_\perp \red {S^2}
\right\}
);

&math(
F=\frac{N_w}{\lambda}\left\{\hbar S\gamma B\,\red{+\alpha\frac{\hbar S}{K_\perp S^2}\hbar S\gamma \dot B}\right\}
);

となりそうに思える。

教科書とは大きく異なってしまったが、
- (5.40), (5.41) からすれば &math(\ddot X); の係数に &math(1+\alpha^2); 
が現れるのは間違いないような？？？~
(&math(\alpha\ll 1); として &math(\alpha^2); を無視しているのかもしれない？)
- &math(S); は必ず &math(\hbar S); 
あるいは &math(K_\perp S^2); の形、あるいはその比 &math(K_\perp/\hbar^2); 
の形で現れるはずで、&math(1/\tau_w); の &math(S); の付き方も上記が正しいのではないか
- 後に磁場の強さを変化させる場合について考察するなら、その時間微分成分を落とすべきではない~
(&math(K_\perp); が大きいとして無視しているのかもしれない？)

*** 摩擦のある時の粒子の運動 [#m96dc9f2]

ピン止めポテンシャルを無視すれば、

&math(
&\red{(1+\alpha^2)}\frac{\hbar^2N_w}{K_\perp\lambda^2}\left[\ddot X
+\frac{\alpha}{\red{1+\alpha^2}}\frac{K_\perp \red S}{\hbar}\dot X
\right]=\frac{N_w \hbar S\gamma B}{\lambda }\,\red{+\frac{\alpha\hbar^2 N_w}{K_\perp\lambda}\gamma \dot B}\\
);

となって、これは

&math(M_w\left(\ddot X+\frac{1}{\tau_w}\dot X\right)=F(t));

の形で摩擦のある時の粒子の方程式そのものである（ように見える）。

これをそのまま信じれば、教科書にあるとおり初速を持った磁壁は
外力を加えない限り慣性により等速運動をすることになる。
摩擦を無視できるとき慣性により等速運動をすることになる。

*** 等速運動をするときのスピンの動き [#ac0b69f5]

上の条件で &math(X); が等速運動をするとき、個々のスピンの動きはどうなるか？

- 磁壁が近づいてくるまで、スピンは容易軸方向（例えば上向き）に静止している
- 磁壁が通過するに伴い &math(\phi_0); をほぼゼロに保ちながら徐々に下を向く
- 磁壁が通過した後は下を向いて静止する

となるはずだ。

個々のスピンにそのような動きをさせる要素が存在するというのはにわかには信じられない。
このような複雑な動きを説明するには
「歳差運動を始めたスピンは回り続ける」というような簡単な話では片付かない。

「歳差運動を始めたスピンは回り続ける」というような簡単な話では片付かないような？

*** 「等速運動を続ける」は限定的にしか成立しない [#xfaf1080]

(5.40) より

&math(\dot \phi_0=-\alpha\frac{1}{\lambda}\dot X+\gamma B);

であるから、外力ゼロで &math(\dot X); が一定なら &math(\dot \phi_0); も一定値となる。

すなわち時間が経つと上で用いた

&math(\phi_0\ll 1);

の条件を満たさなくなってしまう。

「外力を切ってしまえば慣性が働いて &math(\dot X); が（ゼロでない）一定値に保たれる」
というのは ''動き始めの短い時間しか成り立たない'' ことに注意が必要だ。

ピン止めポテンシャルがある場合、
上式の右辺がゼロを中心として振動し、
結果的に &math(\phi_0\ll 1); の条件が崩れないことが、
長時間経過後にも粒子的な運動方程式が成り立つための条件になっている。

*** 短い間だけでも慣性が働く理由？ [#hd423bec]

教科書は磁壁が慣性質量を持つ理由として、

- 磁化（局在スピン）が角運動量でありトルクの運動方程式によって記述されること

としているが、それだと (5.33) や (5.47) で慣性の見られない解が得られていることと矛盾しないか？

教科書の記述をよく見てみる。

- 静止した磁壁内で、局在スピンは磁化困難軸に垂直な平面内を向く~
（向いている方向は静止磁壁内での、その位置における有効磁場方向と一致する）
- 容易軸方向の外部磁場で磁壁が並進運動をする際には、(なぜ条件を限定？)
-- 局在スピンはほぼその面内で回転するが(&math(\phi_0\sim 0);の条件)
-- これはトルクの運動であるために必ず面に垂直な方向にまず立ち上がる~
(外部磁場の印加で有効磁場が容易軸方向へ倒れ、その周りで歳差運動を始める。
容易軸方向へ倒れた磁場は同じ平面内にあるため、
倒れた有効磁場と磁化との外積は平面に垂直になる)
- 垂直方向の局在スピンの立ち上がりと、それによる磁壁の並進運動を組み合わせて考える~
(どういうことか？)
-- いったん運動が始まれば外部磁場を切った後も局在スピンの垂直成分が残り、~
(局在スピンの方向が平面外へ出た後、外部磁場が切れたらどうなるか？
単なる磁壁の並進運動であれば、並進後も局在スピンは平面内にある。
至る所局在スピンが平面外へ出ているときに、各所の有効磁場がどうなるかがすぐには分からない。)
-- 摩擦が無ければ等速度で磁壁が運動することがすぐに分かる~
(分からない)

この部分は良く分からない。

実際にこのモデルを使ってシミュレーションで確かめてみた結果はこちら → 
[[スピントロニクス理論の基礎/X-6]]

** φ__0__粒子 [#z4982fc3]

ピン止め力が &math(X); に比例するようなポテンシャル中の振動を考える。
すなわち、(5.40) で &math(X\theta(\xi-|X|)); となっていた部分を
単なる &math(X); で置き換えられる場合、ということになる。

どういう条件でそのような近似が使えるのかはまた後で考えるが、
教科書では「ピン止めが強い領域」とされている。

また、&math(X); 粒子の方程式として紹介された (5.51) を得る手順に於いても、
実は暗黙のうちに上記の近似が使われているので、
(5.51) は (5.53) にさらに近似を加えたことになる？

(5.40) を

&math(\dot \phi_0+\alpha\frac{1}{\lambda}\dot X=\gamma B-\frac{2\lambda V_0}{\hbar S\xi^2}X);

と書き直し、(5.41) を代入すれば、

(5.53)

&math(
&\dot \phi_0+\alpha\frac{1}{\lambda}\dot X=\gamma B-\frac{2\lambda V_0}{\hbar S\xi^2}X\\
&\dot \phi_0+\alpha\frac{1}{\lambda}\left\{
\alpha\lambda\dot\phi_0+\frac{K_\perp\lambda}{2\hbar}S\sin 2\phi_0
\right\}=\gamma B-\frac{2\lambda V_0}{\hbar S\xi^2}X\\
&(1+\alpha^2)\dot \phi_0+
\frac{\alpha K_\perp S}{2\hbar}\sin 2\phi_0
=\gamma B-\frac{2\lambda V_0}{\hbar S\xi^2}X\\
&(1+\alpha^2)\ddot \phi_0+
\frac{\alpha K_\perp S}{2\hbar}2\dot\phi_0\cos 2\phi_0
=-\frac{2\lambda V_0}{\hbar S\xi^2}\dot X\\
&=-\frac{2\lambda V_0}{\hbar S\xi^2}\left\{
\alpha\lambda\dot\phi_0+\frac{K_\perp\lambda}{2\hbar}S\sin 2\phi_0
\right\}\\
&(1+\alpha^2)N_w\frac{\hbar^2 S^2\xi^2}{2\lambda^2 V_0}\ddot \phi_0
+\alpha N_w\hbar S\frac{\xi^2}{\lambda^2}\frac{\frac{1}{2}K_\perp S^2}{V_0}\dot\phi_0\cos 2\phi_0\\
&=-N_w\frac{\hbar S}{\lambda}\left\{
\alpha\lambda\dot\phi_0+\frac{K_\perp\lambda}{2\hbar}S\sin 2\phi_0
\right\}\\
&=-\alpha\hbar SN_w\dot\phi_0
-N_w\frac{K_\perp S^2}{2}\sin 2\phi_0\\
&(1+\alpha^2)N_w\frac{\hbar^2 S^2\xi^2}{2\lambda^2 V_0}\ddot \phi_0
+\alpha\hbar SN_w\dot\phi_0
+N_w\frac{K_\perp S^2}{2}\sin 2\phi_0
{\color{magenta}{+\alpha N_w\hbar S\frac{\xi^2}{\lambda^2}\frac{\frac{1}{2}K_\perp S^2}{V_0}\dot\phi_0\cos 2\phi_0}}=0\\
);

最後の式で色つきの部分を &math(K_\perp/V_0\ll 1); の条件により無視すれば、
ほぼ教科書の式となるが、正しくは

(5.54)

&math(
M_\phi=\red{(1+\alpha^2)}N_w\frac{\hbar^2S^2\xi^2}{2V_0\lambda^2}
);

である。

(5.52) でも &math(1+\alpha^2); のファクターが抜けていたところを見ると、
&math(\alpha\ll 1); を仮定して &math(\alpha^2); を無視しているのかもしれない？

* 質問・コメント [#c31d3b3c]

#article_kcaptcha
Counter: 6085 (from 2010/06/03), today: 1, yesterday: 4