マイクロ磁気学のソースを表示

{{翻訳直後|[[:en:Special:Permalink/721248697|en:Micromagnetics]]|date=2016年10月}}
'''マイクロ磁気学'''（マイクロじきがく、'''マイクロマグネティクス''' {{Lang-en-short|Micromagnetics}} とも）とは[[物理学]]の一分野で、サブ[[マイクロメートル]]スケールの[[磁性|磁気]]的挙動の予測を取り扱う。ここで取り扱われる長さスケールは材料の[[原子]]的構造を無視できる（[[連続体力学|連続体近似]]）程度には大きいが、[[磁壁]]や[[磁気渦]]を取り扱える程度には小さい。

マイクロ磁気学は、磁気エネルギーを最安定化することで{{仮リンク|力学的平衡|en|Mechanical equilibrium}}を、また時間依存動力学方程式を解くことにより動的挙動を扱うことができる。

== 歴史 ==
物理学の一分野としてのマイクロ磁気学（すなわち、（[[強磁性|強]]）[[磁性体]]のサブマイクロメートル領域における振る舞いを扱う分野）は[[1963年]]、{{仮リンク|ウィリアム・フラー・ブラウン・ジュニア|en|William Fuller Brown, Jr.}}が反平行磁壁構造についての論文を発表したことに始まる。比較的最近まで、計算マイクロ磁気学には計算コスト的な困難が付き纏っていたが、現在では最新型のデスクトップ[[パーソナルコンピュータ|PC]]で小さめの問題なら解くことが可能になってきている。

== 静的マイクロ磁気学 ==
静的マイクロ磁気学の目的は、力学的平衡における[[磁化]] {{Mvar|'''M'''}} の空間的分布を解明することである。ほとんどの場合、扱われる材料の[[温度]]は[[キュリー温度]]よりもかなり低いため、磁化係数 {{Math|{{mabs|'''''M'''''}}}} はいたるところ{{仮リンク|飽和磁化|en|Saturation magnetization}} {{Math|''M''<sub>s</sub>}} と等しい。この場合、問題は、「磁化配向ベクトル」もしくは「規格化磁化」と呼ばれる {{Math|1='''''m''''' = '''''M'''''/''M''<sub>s</sub>}} により記述される、磁化の空間的配向を解明することに帰着する。

力学的平衡は、次のように表わされる磁気エネルギーを最安定化することで達成される。
: <math>E = E_\text{exch} + E_\text{anis} + E_\text{Z} + E_\text{demag}+E_\text{m-e}</math>
ここで、拘束条件として {{Math|1={{mabs|'''''M'''''}}=''M''<sub>s</sub>}} もしくは {{Math|1={{mabs|'''''m'''''}} = 1}} を課す。

このエネルギーへの寄与それぞれについて以下に説明する。

=== 交換エネルギー ===
交換エネルギーは、[[量子力学]]における[[交換相互作用]]を現象論的かつ連続体的に表わしたものであり、次のように書ける。
: <math>E_\text{exch} = A \int_V \left((\nabla m_x)^2 + (\nabla m_y)^2 + (\nabla m_z)^2\right) \mathrm{d}V</math>
ここで {{Mvar|A}} は「交換定数」、{{Mvar|m<sub>x</sub>, m<sub>y</sub>, m<sub>z</sub>}} はそれぞれ {{Mvar|'''m'''}} の方向成分である。積分範囲は試料の体積全体に渡るものとする。

交換エネルギーは試料中の磁化配向ができるだけゆっくり変化することを選好するように作用する。このエネルギーは磁化が完全に一様となる際に最小となる。

=== 異方性エネルギー ===
[[磁気異方性]]は[[結晶構造]]と[[スピン軌道相互作用]]の組み合わせから生じる。一般的に書くと以下のようになる。
: <math>E_\text{anis} = \int_V F_\text{anis}(\mathbf{m}) \mathrm{d}V</math>
ここで、{{Math|''F''<sub>anis</sub>}} は異方性エネルギー密度で、磁化配向の関数となる。{{Math|''F''<sub>anis</sub>}} が最低エネルギーとなるような配向が沿う軸は「磁化容易軸」と呼ばれる。

[[T対称性|時間反転対称性]]により、{{Math|''F''<sub>anis</sub>}} は {{Mvar|'''m'''}} の[[偶関数と奇関数|偶関数]]であることが保証される。これを満たす最も単純な関数は以下のような関数である。
: <math>F_\text{anis}(\boldsymbol{m}) = -K m_z^2</math>
ここで、{{Mvar|K}} は「異方性定数」と呼ばれる。この近似は「一軸磁気異方性」と呼ばれ、上式の場合の容易軸は {{Mvar|z}} 軸である。

異方性エネルギーは容易軸に沿った磁化配向が選好されるように作用する。

=== ゼーマンエネルギー ===
ゼーマンエネルギーは磁化と外部[[磁場]]との相互作用エネルギーである。これは次のように書ける。
: <math>E_\text{Z} = -\mu_0 \int_V \boldsymbol{M}\cdot\boldsymbol{H}_\text{a} \boldsymbol{d}V</math>
ここで、{{Math|'''''H'''''<sub>a</sub>}} は外部磁場、{{Math|''µ''<sub>0</sub>}} は{{仮リンク|真空透磁率|en|Vacuum permeability}}である。

ゼーマンエネルギーは外部磁場と平衡な磁化配向が選好されるように作用する。

=== 反磁場エネルギー ===
[[ファイル:Micromagnetics_by_Zureks.svg|サムネイル|マイクロ磁気配向の例。 一様状態と比べて、磁束閉鎖構造により反磁場エネルギーは低下するが交換エネルギーは上昇する。]]
反磁場とは磁性試料が自身に及ぼす磁場である。反磁場エネルギーは以下のように書ける。
: <math>E_\text{demag} = -\frac{\mu_0}{2} \int_V \boldsymbol{M}\cdot\boldsymbol{H}_\text{d} \mathrm{d}V</math>
ここで、{{Math|'''''H'''''<sub>d</sub>}} は{{仮リンク|反磁場|en|Demagnetizing field}}である。この場は磁化配向そのものに依存し、次の方程式系を解くことにより得られる。
: <math>\nabla\cdot\boldsymbol{H}_\text{d} = -\nabla\cdot\boldsymbol{M}</math>
: <math>\nabla\times\boldsymbol{H}_\text{d} = 0</math>
ここで、 {{Math|&minus;&nabla; &middot; '''''M'''''}} は「磁荷密度」と呼ばれることがある。これらを解くと、以下のようになる（{{仮リンク|静磁場|en|Magnetostatics|preserve=yes}}の項も参照）。
: <math>\boldsymbol{H}_\text{d} = -\frac{1}{4\pi} \int_V \nabla\cdot\boldsymbol{M} \frac{\boldsymbol{r}}{r^3} \mathrm{d}V</math>
ここで、 {{Math|'''''r'''''}} は積分点から計算する {{Math|'''''H'''''<sub>d</sub>}} の場所へ向うベクトルである。

試料の端では {{Math|'''''M'''''}} が試料内における有限値から試料の外におけるゼロへと[[連続 (数学)|不連続]]に変化することから、磁荷密度が[[無限|無限大]]になりうることは特筆に値する。この現象は、適切な[[境界値問題|境界条件]]を用いることによって取り扱うことが通常は可能である。

反磁場エネルギーは磁荷密度が最小となるような磁荷配向が選好されるように作用する。特に、試料の端では磁荷は表面と平行となる傾向がある。ほとんどの場合、他のエネルギーを最低にしつつ同時にこのエネルギーを最低にすることは不可能である。この結果、力学的平衡は個々の項を最小化するのではなく、総磁気エネルギーを最小化するためにそれぞれの間で妥協して成り立つことになる。

=== 磁気弾性エネルギー  ===
磁気[[弾性エネルギー]]とは、[[結晶構造|結晶格子]]の[[弾性|弾性変形]]により蓄えられるエネルギーを指す。磁気弾性に関わる効果が無視できる場合は無視できる。磁荷配向 {{Mvar|'''m'''}} に付随する、結晶性固体の局所[[ひずみ|歪み]]には選好がある。単純なモデルとして、この歪みを等積的で、横方向には完全に[[等方的と異方的|等方的]]であるとすると、次の偏差モデル{{訳語疑問点|date=2016年10月|}} (deviatoric ansatz) を得る。

<math> \mathbf{\varepsilon}_0(\boldsymbol{m}) = \frac{3}{2} E\, [\boldsymbol{m}\otimes \boldsymbol{m} - \frac{1}{3}\mathbf{1}]</math>

ここで、材料のパラメータ {{Math|''E'' > 0}} は磁気歪み定数と呼ばれる。明らかに、{{Mvar|E}} は磁化により {{Mvar|'''m'''}} の方向に誘起される歪みである。このモデルを用いれば、弾性エネルギー密度は弾性的で応力を伴う歪み{{訳語疑問点|date=2016年10月|}} (elastic, stress-producing strains) <math> \mathbf{\varepsilon}_e := \mathbf{\varepsilon} -\mathbf{\varepsilon}_0</math> の関数と考えることができる。[[二次形式]]の磁気弾性エネルギーは以下のようになる。

<math> E_\text{m-e} = \frac{1}{2} [\mathbf{\varepsilon} -\mathbf{\varepsilon}_0(\boldsymbol{m})] : \mathbb{C} :
[\mathbf{\varepsilon} -\mathbf{\varepsilon}_0(\boldsymbol{m})] </math>

ここで、<math>\mathbb{C} :=\lambda \mathbf{1}\otimes \mathbf{1} + 2\mu \mathbb{I}</math> は四次の弾性[[テンソル]]である。さらに、弾性応答は（二つの[[ラメ定数]] λ と μに基いて）等方的であることを仮定すると、{{Mvar|'''m'''}} の長さが一定であることを考慮に入れれば、次の不変量に基く表現が可能である。

<math>
E_\text{m-e}= 
\frac{\lambda}{2} \mbox{tr}^2[\mathbf{\varepsilon}] 
+ \mu \, \mbox{tr}[\mathbf{\varepsilon}^2] 
- 3\mu E \big\{ \mbox{tr}[\mathbf{\varepsilon}(\boldsymbol{m}\otimes\boldsymbol{m})] 
- \frac{1}{3}\mbox{tr}[\mathbf{\varepsilon}] \big\}
</math>

このエネルギー項が磁気歪みに寄与する。

== 動的マイクロ磁気学 ==
動的マイクロ磁気学の目的は、例えば[[パルス]]磁場や[[交流]]磁場などの[[非定常状態|非定常]]な条件下での試料の磁気配向の[[時間発展]]を予測することである。このためには、磁化の時間発展を試料に作用する局所「有効磁場」を用いて記述する[[偏微分方程式]]である[[ランダウ＝リフシッツ＝ギルバート方程式|ランダウ・リフシッツ・ギルバート方程式]]を解くことになる。

=== 有効磁場 ===
'''有効磁場'''とは、磁化が「感じる」局所的磁場を指す。これは、厳密ではないが、次のように磁気エネルギーの磁化配向に関する[[微分|導関数]]を用いて表わすことができる。
: <math>\boldsymbol{H}_\mathrm{eff} = - \frac{1}{\mu_0 M_s} \frac{\mathrm{d}^2 E}{\mathrm{d}\boldsymbol{m}\mathrm{d}V}</math>
ここで、{{Math|d''E''/d''V''}} はエネルギー密度である。[[変分法]]を用いると、磁化の変化 {{Math|d'''''m'''''}} に対応する磁気エネルギー変化 {{Math|d''E''}} の関係は以下のように表わせる。
: <math>\mathrm{d}E = -\mu_0 M_s \int_V (\mathrm{d}\boldsymbol{m})\cdot\boldsymbol{H}_\text{eff}\,\mathrm{d}V</math>
{{Mvar|'''m'''}} は[[単位ベクトル]]であるから、{{Math|d'''''m'''''}} は常に {{Mvar|'''m'''}} に[[直交]]する。そのため、上記の定義は {{Math|'''''H'''''<sub>eff</sub>}} の {{Mvar|'''m'''}} に[[平行]]な成分を含んでいない。この成分は通常、磁化の[[動力学]]に影響しないため問題にはならないことが多い。

磁気エネルギーに対するそれぞれの寄与項の表式を代入すると、有効磁場は次のように表わされる。
: <math>\boldsymbol{H}_\mathrm{eff} = \frac{2A}{\mu_0 M_s} \nabla^2 \boldsymbol{m} - \frac{1}{\mu_0 M_s} \frac{\partial
F_\text{anis}}{\partial \boldsymbol{m}} + \boldsymbol{H}_\text{a} + \boldsymbol{H}_\text{d}</math>

=== ランダウ・リフシッツ・ギルバート方程式 ===
[[ファイル:Damped_Magnetization_Precession.jpg|サムネイル|ランダウ・リフシッツ・ギルバート方程式に表われる項：歳差（赤）、減衰（青）。点線で表わした磁化の[[トラジェクトリ]]は有効磁場 {{Math|'''''H'''''<sub>eff</sub>}} を一定とする単純化された仮定のもとでは[[螺旋]]を描く。]]
ランダウ・リフシッツ・ギルバート方程式は磁化の運動方程式である。有効磁場周りの磁化の[[ラーモア歳差運動]]と、磁性系と環境との間の相互作用により生じる[[減衰振動|減衰]]を組み合わせた運動を記述する。この方程式のいわゆる「ギルバート形式」（もしくは陰形式）は次のように書き下される。
: <math>\frac{\partial \boldsymbol m}{\partial t} = - |\gamma| \boldsymbol{m} \times \boldsymbol{H}_\mathrm{eff} + \alpha \boldsymbol{m}\times\frac{\partial \boldsymbol{m}} {\partial t}</math>
ここで、{{Mvar|γ}} は電子の[[磁気回転比]]、{{Mvar|α}} はギルバート減衰定数である。

上記の形式は次の「ランダウ・リフシッツ形式」（もしくは陽形式）と等価であることを数学的に示すことができる。
: <math>\frac{\partial\boldsymbol m}{\partial t} = - \frac{|\gamma|}{1+\alpha^2} \boldsymbol{m} \times \boldsymbol{H}_\mathrm{eff} - \frac{\alpha|\gamma|}{1+\alpha^2} \boldsymbol{m}\times(\boldsymbol{m}\times\boldsymbol{H}_\text{eff})</math>

== 応用 ==
マイクロ磁気学と力学の組み合わせにより、[[超音波]]スピーカーや高周波磁気歪み[[トランスデューサー]]などの磁気音響共鳴の関わる産業分野への応用が可能であり興味を集めている。この分野では上述した磁気歪み効果を取り込んだマイクロ磁気学モデルに基づく[[有限要素法|FEM]]シミュレーションが重要である<ref name=Miehe2011/>。

従来の磁区と磁壁の他にも、磁気[[渦]]や磁気反渦状態<ref name=Komineas2007/>や、例えば磁化が全ての点で原点と逆もしくは原点に向いている配向やトポロジカルに等価な配向である三次元ブロッホ点<ref name=Thiaville2003/><ref name="Döring1968"/>等のトポロジカルな線状・点状配向にもマイクロ磁気学は適用されている。このため、空間的にも、時間的にもナノ（さらにはピコ）スケールまで利用されている。

対応する[[トポロジカル量子数]]<ref name="Döring1968"/>を情報担体として応用することが、最新の、既に進行中の[[情報技術]]的計画において構想されている。

== 関連項目 ==
* [[磁性]]
* [[磁気スキルミオン]]

== 出典 ==
{{reflist|refs=
<ref name=Miehe2011> {{cite journal | last         = Miehe | first        = Christian | last2        = Ethiraj | first2       = Gautam | date         = 2011-10-15 | title        = A geometrically consistent incremental variational formulation for phase field models in micromagnetics | url          = http://www.sciencedirect.com/science/article/pii/S0045782512000977 | journal      = Computer Methods in Applied Mechanics and Engineering | publisher    = Elsevier | volume       = 245–246 | issue        = | pages        = 331–347 | doi          = 10.1016/j.cma.2012.03.021 | accessdate   = |bibcode = 2012CMAME.245..331M }}</ref>
<ref name=Komineas2007>{{cite arXiv |last=Komineas |first=Stavros|author2=Papanicolaou, Nikos |authorlink= |eprint=0712.3684v1 |title=Dynamics of vortex-antivortex pairs in ferromagnets |class=cond-mat.mtrl-sci |year=2007}}</ref>
<ref name=Thiaville2003>{{cite journal|last=Thiaville|first=André|author2=García, José |author3=Dittrich, Rok |author4=Miltat, Jacques |author5= Schrefl, Thomas |title=Micromagnetic study of Bloch-point-mediated vortex core reversal|journal=Physical Review B|date=March 2003|volume=67|issue=9|doi=10.1103/PhysRevB.67.094410|bibcode = 2003PhRvB..67i4410T }}</ref>
<ref name="Döring1968">{{cite journal|last=Döring|first=W.|title=Point Singularities in Micromagnetism|journal=Journal of Applied Physics|year=1968|volume=39|issue=2|pages=1006|doi=10.1063/1.1656144|bibcode = 1968JAP....39.1006D }}</ref>
}}

== 関連文献 ==

=== 外部リンク ===
* [http://www.ctcms.nist.gov/mumag/mumag.org.html µMAG -- Micromagnetic Modeling Activity Group].
* [http://math.nist.gov/oommf OOMMF -- Micromagnetic Modeling Tool].
* [http://mumax.github.io MuMax -- GPU-accelerated Micromagnetic Modeling Tool].

{{デフォルトソート:まいくろしきかく}}
[[Category:力学系]]