磁壁移動の粘弾性モデル

磁壁移動は、ピン止めポテンシャル中を動く磁壁の座標が、粘性抵抗と遅れを伴って外力に応答する現象として記述できる。粘弾性モデルは、時間応答の形（緩和・遅れ・非対称性）を方程式の構造として固定し、材料・形状・導電率・欠陥がどの項に現れるかを切り分ける枠組みである。

参考ドキュメント

1. B. Alessandro, C. Beatrice, G. Bertotti, and A. Montorsi, Domain-wall dynamics and Barkhausen effect in metallic ferromagnetic materials, Journal of Applied Physics 68, 2901 (1990).
2. G. Durin and S. Zapperi, The Barkhausen effect, arXiv:cond-mat/0404512 (2004).
3. 電磁鋼板における磁区と渦電流・鉄損に関する解説（日本語）, 特殊鋼倶楽部 会誌 PDF（2014年号内の記事を含む）.

1. 位置づけ

軟磁性材料の磁化過程では、磁壁移動と磁化回転が競合し、さらに導電体では渦電流が磁壁運動に反作用を与える。磁壁移動を粘弾性として扱う意義は次の2点にある。

• 時間領域の磁化変化を、粘性（散逸）と弾性（復元力）と遅れ（メモリ）に分解できること
• 同じ測定波形でも、熱活性化、渦電流遅れ、磁壁の内部自由度（角度自由度）を別々の署名として扱えること

2. 粘弾性という言葉を磁壁に写像する

粘弾性は本来、応力とひずみの関係が時間微分や畳み込み（メモリ）を含むことを指す。磁壁の場合、連続体のひずみの代わりに、磁壁中心位置や磁区幅などの集団座標を用いる。

磁壁位置を $q(t)$ とし、外部磁場や内部場により磁壁に働く駆動力を $F(t)$ とする。ピン止めを含むポテンシャルを $U(q)$ とすると、最も基本的な粘弾性対応の運動方程式は

$$m\ddot{q}(t)+\eta \dot{q}(t)+\frac{\partial U}{\partial q}(q)=F(t)+\xi (t)$$

である。ここで

• $m$ は慣性（磁壁質量、あるいは有効質量）に対応する係数
• $\eta$ は粘性抵抗（散逸）に対応する係数
• $\xi (t)$ は欠陥起源の確率的揺らぎや熱揺らぎを表す項

である。

ピン止めポテンシャルが局所的に放物線近似できる範囲では

$$U(q)\simeq \frac{1}{2}k(q-q_0{)}^2$$

となり、$k$ が弾性（復元力の強さ）を担う。このとき、外力が取り去られた後の緩和は、過減衰なら

$$q(t)-q_0\propto e^{-t/\tau },\tau =\frac{\eta }{k}$$

の指数緩和を与える。磁壁の時間定数を測って $\tau$ を得ることは、粘性と復元力の比を観測したことに等しい。

3. 磁壁に働く外力の表式

磁壁が ${180}^{\circ }$ の反転境界で、外部磁場 $H$ が磁化容易方向に沿うとき、ゼーマンエネルギー密度は $-{\mu }_0\mathbf{M}\cdot \mathbf{H}$ であり、両側磁区のエネルギー密度差は $2{\mu }_0{M}_{s}H$ である。したがって磁壁に働く圧力（単位面積当たりの駆動力）は

$$p_H=2{\mu }_0{M}_{s}H$$

となり、磁壁面積を $A$ とすれば駆動力は $F_H=p_{H}A$ である。形状や多磁壁系では、反磁場や相互作用が $U(q)$ と $F(t)$ の両方に入るため、単一磁壁の $p_H$ は局所近似として使うのが自然である。

4. 粘性項 $\eta$ の物理的起源

磁壁の粘性抵抗は単一機構ではなく、複数の散逸が同一の $\eta$ に畳み込まれる。支配機構を切り分けると、粘弾性モデルの解釈が安定する。

4.1 磁化ダイナミクス由来（ギルバート減衰）

連続体磁化 $\mathbf{M}(\mathbf{r},t)$ は LLG 方程式

$$\frac{d\mathbf{M}}{dt}=-\gamma \mathbf{M}\times {\mathbf{H}}_{\mathrm{eff}}+\frac{\alpha }{M_s}\mathbf{M}\times \frac{d\mathbf{M}}{dt}$$

で記述され、$\alpha$ が散逸の大きさを担う。磁壁を1次元集団座標で表すと、低磁場での定常運動（ウォーカー崩壊の手前）では、磁壁速度 $v=\dot{q}$ が

$$v\propto \frac{\gamma \mathrm{\Delta }}{\alpha }H$$

の形で磁場に比例し、$\alpha$ が粘性的な抵抗として現れる。

4.2 導電体での渦電流由来（電磁場の反作用）

導電体中で磁壁が動くと、磁束変化により渦電流が誘起され、そのジュール損と自己誘導により磁壁の運動に反作用が生じる。重要なのは、渦電流の反作用が速度に比例する抵抗に見えるだけでなく、有限の電磁拡散時間を通して遅れ（メモリ）を伴う点である。

軟磁性電磁鋼板などでは、渦電流が磁壁近傍に集中し、磁区幅や磁壁配置が損失低減と直結することが工学的に整理されている。磁区幅を細分化すると磁壁がゆっくり動けばよくなり、渦電流が減って損失が下がるという説明は、磁壁速度が渦電流損に直結する事実を反映している。

5. 遅れ（メモリ）を含む粘弾性方程式

渦電流の反作用や、内部自由度の緩和が無視できない場合、抵抗は瞬時の速度だけでなく過去の速度履歴に依存する。一般化すると

$$m\ddot{q}(t)+{\eta }_0\dot{q}(t)+{\int }_0^t\mathrm{\Gamma }(t-t^{\prime })\dot{q}(t^{\prime })d{t}^{\prime }+\frac{\partial U}{\partial q}(q)=F(t)+\xi (t)$$

となる。$\mathrm{\Gamma }(t)$ はメモリ核であり、指数型

$$\mathrm{\Gamma }(t)=\frac{{\eta }_1}{\tau }{e}^{-t/\tau }$$

を仮定 distinguish すれば、これは標準線形固体（Zener 型）に相当する粘弾性応答となる。すなわち、短時間では ${\eta }_0+{\eta }_1$ の強い抵抗、長時間では ${\eta }_0$ の抵抗へと移る。

周波数領域では

$$[-m{\omega }^2+i\omega \eta (\omega )+k]q(\omega )=F(\omega )$$

と書け、複素粘性 $\eta (\omega )$ が損失（虚部）と遅れ（周波数分散）を担う。磁壁移動の位相遅れや、バルク磁化の複素透磁率の分散は、この形で統一的に解釈できる。

6. 慣性と有効質量

磁壁に慣性が現れる経路は少なくとも2つある。

• 磁壁内部の角度自由度（ウォール内部の磁化回転）が時間発展し、集団座標 $q$ に二階時間微分項が誘起されること
• 導電体での渦電流反作用が、メモリ核として現れる結果、短時間応答が慣性のように見えること

この文脈では、有効質量が負になる署名が議論されており、バークハウゼン信号のパルス形状の左右非対称性（左に偏った立ち上がり）を負の有効質量の効果として結びつける理論が提示されている。

7. ピン止め・熱活性化・クリープと粘弾性

欠陥・析出物・内部応力は、磁壁にランダムなピン止めポテンシャル $U(q)$ を与える。ここに温度が入ると、磁壁は障壁を熱活性化で越え、時間依存が強くなる。粘弾性モデルの $\xi (t)$ は熱揺らぎの効果も表現でき、平均速度 $v$ は熱活性化則を含む非線形応答になる。

弾性界面のランダム媒体中の運動として整理すると、駆動磁場 $H$ に対して

• クリープ領域（小 $H$）

$$v(H)\sim \exp [-\frac{U_c}{k_{B}T}{\left(\frac{H_c}{H}\right)}^{\mu }]$$

• デピニング近傍

$$v(H)\sim (H-H_c{)}^{\beta }$$

• フロー領域（大 $H$）

$$v(H)\propto H$$

のようなクロスオーバーが現れる。粘弾性モデルは、このうちフロー領域を線形粘性として最も単純に表し、クリープ領域を確率過程（活性化）として追加する形になる。

8. バルクハウゼン効果と粘弾性拡張（ABBM と遅れ）

バルク軟磁性体のバークハウゼンノイズは、磁壁がピン止め配置間を雪崩的に遷移する確率過程として理解されてきた。ABBM（Alessandro–Beatrice–Bertotti–Montorsi）モデルは、磁壁位置（または磁化）を1自由度で扱い、ランダムな有効力と粘性抵抗の下でのランジュバン記述から、雪崩統計を導く。

最も単純な AB BM 型の骨格は、過減衰近似で

$$\eta \dot{q}=F_{\mathrm{drive}}(t)-F_{\mathrm{pin}}(q)$$

のように書かれ、$F_{\mathrm{pin}}(q)$ をブラウン運動的なランダム力として与える。これに渦電流遅れを入れると、上で述べたメモリ核が追加され、

• 雪崩パルスの左右非対称性
• パルス後尾の指数的緩和
• 周波数スペクトルの形の変化

などが、粘弾性の時間定数として現れる。

負の有効質量がパルス非対称性を与えるという議論は、遅れの存在を「慣性に見える成分」として抽出した例と位置づけられる。

9. 何を測ると何が分かるか

粘弾性モデルのパラメータは、測定量に対応づけると設計指針として使える。

9.1 速度–磁場特性 $v(H)$

低磁場フローで $v\propto H$ が成立する範囲の傾きは、実効移動度（mobility）に対応し、$\eta$ の逆数として働く。材料や温度、応力で傾きが変わる場合、ギルバート由来の散逸か、渦電流由来か、ピン止め由来かを分けて考える必要がある。

9.2 時間波形（MBN パルスの立ち上がりと減衰）

単発パルスの立ち上がりと後尾に指数関数が見える場合、$e^{-t/\tau }$ の時間定数は $\tau =\eta /k$（過減衰）あるいはメモリ核の $\tau$ を反映する。非対称性が一貫して観測される場合、遅れ（メモリ）を含む運動方程式が必要になる。

9.3 厚さ・導電率依存

渦電流に起因する遅れは、導電率と形状（特に厚さ）に強く依存する。厚さや抵抗率を系統的に変えることで、${\eta }_0$（磁化ダイナミクス由来）とメモリ核（渦電流由来）を分離しやすくなる。

9.4 応力依存と磁気粘性

応力誘起異方性が変わると、磁壁の弾性（$k$）とピン止め障壁（$U_c$）が変わり、クリープ則や雪崩統計の指数が変化し得る。応力下バークハウゼンをデピニング像で解釈する研究もあり、粘弾性モデルはその時間領域表現として整合する。

10. ばね・ダンパ・メモリ核として見た整理

力学要素	方程式の形	磁壁での対応	観測されやすい署名
ばね	$F=-k(q-q_0)$	局所ピン止め剛性、反磁場による復元力	指数緩和、共鳴条件の移動
ダンパ（瞬時）	$F=-{\eta }_0\dot{q}$	ギルバート散逸、微視的散逸の集約	$v\propto H$ の線形域
メモリ核	$F=-\int \mathrm{\Gamma }(t-t^{\prime })\dot{q}(t^{\prime })d{t}^{\prime }$	渦電流反作用、内部自由度の緩和	パルス非対称、位相遅れ、周波数分散
質量	$F=-m\ddot{q}$	磁壁内部自由度起因の慣性、遅れの短時間極限	高速駆動での立ち上がりの曲率

11. 鉄損・動的損失との接続

バルク導電体では、磁壁近傍で渦電流が発生し、磁壁の運動が遅くなるほど損失は小さくなる方向に働く一方、磁化過程を同じ時間内に完了するためには駆動を強める必要が出る。したがって、損失最小化は

• 磁区構造（磁区幅、磁壁密度）による渦電流経路の制御
• ピン止め分布（欠陥密度、内部応力）による雪崩性の制御
• 散逸機構（導電率、磁気減衰、磁歪・磁気弾性結合）の相対重みの設計

の同時最適化になる。磁壁移動の粘弾性モデルは、この同時最適化を「どの項を変えればよいか」という形に翻訳する役割を持つ。

12. まとめと展望

磁壁移動の粘弾性モデルは、磁壁座標に対する運動方程式を、弾性（復元力）、粘性（散逸）、慣性（有効質量）、遅れ（メモリ核）として分解して記述する枠組みである。バルク導電体では渦電流がメモリ核として現れ、バークハウゼン波形の非対称性や緩和時間として観測され得る。

今後は、時間分解の高いバークハウゼン計測や局所磁区観察と組み合わせ、メモリ核の形（単一緩和か多重緩和か）を材料・厚さ・応力で系統的に同定することが重要になる。さらに、粘弾性項を第一原理計算やマルチフィジックス電磁場計算と接続し、散逸の起源を電子論・形状・組織へ分配することで、低損失材料設計へ直接つながる記述へ発展すると期待される。

その他参考にしたsources

• S. Zapperi, P. Cizeau, G. Durin, and H. E. Stanley, Avalanches, depinning transition, and the Barkhausen effect, Physical Review B 58, 6353 (1998).
• S. Zapperi et al., Signature of effective mass in crackling noise asymmetry, arXiv:cond-mat/0507607.
• C. Beatrice and G. Bertotti, Influence of magnetic viscosity on domain wall dynamics and Barkhausen effect in metallic ferromagnetic systems, Journal of Magnetism and Magnetic Materials 92, 817 (1992).
• A. Barman and A. Haldar, Time-domain study of magnetization dynamics in ferromagnetic nanostructures（書籍章PDF内の磁壁移動度と粘性の記述を参照）.
• P. Le Doussal and K. J. Wiese ほか, disordered elastic systems/depinning に関する総説 PDF.