方程式に対する差分スキームについて1Landau-Lifshitz方程式について.pdfVIP

数理解析研究所講究録 第 1543 巻 2007 年 245-253 245 Landau-Lifshitz 方程式に対する差分スキームについて 武蔵野大学 秋山高宏 (Takahiro Akiyama)l Musashino university 岐阜大学 石渡哲哉 (Tetsuya Ishiwata)2 Gifu university みずほ情報総研 不破敦 (Atsushi Fuwa) Mizuho InformationResearch Institute 早稲田大学 堤正義 (Masayoshi Tsutsumi) Waseda university 概要 強磁性体中のス???ンの運動を記述するランダウ. リフシッツ方和式に対する差分 スキームを提案する。 この方程式は (1) 解の長さを保存する、 (2) エネルギー保存系あるいは散逸系である、 という性質をもつ。本講演では、提案する差分スキームがこの両方の性質を元の方程 式から継承していることを示す。 またそのスキームは陰的非線形であるので、 差分解 の一意可解性等についても言及する。 最後に、 厳密解に対する数値計算結果を紹介 する。 1Landau-Lifshitz 方程式について 強磁性体の中の電子スピンの動きを考える。歳差運動による制動を考慮した電子スピン の運動モデルとして、 Landau-Lifshitz 方程式 $\frac{\partial u}{\partial t}=u\cross\Delta u-\mu ux(u\cross\Delta u)$ (1) が知られている。 ここで, $\Omega\subseteq \mathbb{R}^{N}$ は強磁性体領域を表し、 $u=u(x, t)=(u_{1}(x, t)$ , $u_{2}(x, t),u_{3}(x, t))$ : $\Omega x(0, \infty)arrow \mathbb{R}^{3}$ は, 電子スピンの向きを表すベクトルとする。 ま た. $\Delta=\sum_{:=1}^{N}\frac{\partial^{2}}{\partial x_{1}^{2}}$ であり、 $\cross$ は外積を表すものとす乱 この方程式の右辺第二項は Gilbert 制動項と呼ばれ、係数 $\mu$ は非負とする。 $\mu=0$ の場合は歳差運動から生じる制動を考慮し ていないモデルであり、 Heiscnberg 方程式と呼ばれる。 以下、 Landau-Lifshitz 方程式の解がもつ重要な性質として、 長さとエネルギーに関す る性質について述べる。初期値