(57)【特許請求の範囲】
【請求項１】
コンピュータに、
有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する第１手順、及び
前記補間された磁化ベクトルに加わる力として作用する磁界を前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合磁界を算出する第２手順
を実行させることを特徴とする磁気交換結合エネルギー算出プログラム。
【請求項２】
前記第２手順に代えて、前記補間された磁化ベクトルに加わる力として作用する磁界を前記分割された直線上の領域毎に数値積分することにより、前記磁気交換結合磁界を算出する第３手順を備えることを特徴とする請求項１に記載の磁気交換結合エネルギー算出プログラム。
【請求項３】
前記算出された磁気交換結合磁界をLLG（Landau-Lifshitz-Gilbert）方程式の有効磁界に含まれる交換結合磁界として使用することにより、マイクロ磁化解析を行うことを特徴とする請求項１又は２に記載の磁気交換結合エネルギー算出プログラム。
【請求項４】
コンピュータに、
有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する第１手順、及び
前記補間された磁化ベクトルが占める体積のエネルギーを前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合エネルギーを算出する第２手順
を実行させることを特徴とする磁気交換結合エネルギー算出プログラム。
【請求項５】
前記第２手順に代えて、前記補間された磁化ベクトルが占める体積のエネルギーを前記分割された直線上の領域毎に数値積分することにより、前記磁気交換結合エネルギーを算出する第３手順を備えることを特徴とする請求項４に記載の磁気交換結合エネルギー算出プログラム。
【請求項６】
有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する第１手順、及び
前記補間された磁化ベクトルに加わる力として作用する磁界を前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合磁界を算出する第２手順
を備えることを特徴とする磁気交換結合エネルギー算出方法。
【請求項７】
前記第２手順に代えて、前記補間された磁化ベクトルに加わる力として作用する磁界を前記分割された直線上の領域毎に数値積分することにより、前記磁気交換結合磁界を算出する第３手順を備えることを特徴とする請求項６に記載の磁気交換結合エネルギー算出方法。
【請求項８】
前記算出された磁気交換結合磁界をLLG（Landau-Lifshitz-Gilbert）方程式の有効磁界に含まれる交換結合磁界として使用することにより、マイクロ磁化解析を行うことを特徴とする請求項６又は７に記載の磁気交換結合エネルギー算出方法。
【請求項９】
有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する第１手順、及び
前記補間された磁化ベクトルが占める体積のエネルギーを前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合エネルギーを算出する第２手順
を備えることを特徴とする磁気交換結合エネルギー算出方法。
【請求項１０】
前記第２手順に代えて、前記補間された磁化ベクトルが占める体積のエネルギーを前記分割された直線上の領域毎に数値積分することにより、前記磁気交換結合エネルギーを算出する第３手順を備えることを特徴とする請求項９に記載の磁気交換結合エネルギー算出方法。
【請求項１１】
有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する補間手段と、
前記補間された磁化ベクトルに加わる力として作用する磁界を前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合磁界を算出する算出手段とを備えることを特徴とする磁気交換結合エネルギー算出装置。
【請求項１２】
前記算出手段は、前記補間された磁化ベクトルが占める体積のエネルギーを前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合エネルギーを算出することを特徴とする請求項１１記載の磁気交換結合エネルギー算出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気交換結合エネルギー算出プログラム、算出方法及び算出装置に関する。
【背景技術】
【０００２】
従来より、磁性体を含む複数の物質により占められている等価要素内の平均磁束密度及び平均磁界を演算し、上記平均磁束密度及び平均磁界を用いて、上記等価要素よりも広い領域に生じる電磁場を解析する装置が知られている。また、従来より、電磁場解析対象領域を微小領域に分割し、微小領域における最小磁束密度と最大磁束密度との比αをパラメータとするＨ−Ｂ曲線及びＷ−Ｂ曲線をデータベースに格納し、電磁場の解析を実行するシステムが知られている。さらに、従来より、磁性体材料内において原子の磁気モーメントの向きが同じ向きに揃った領域である磁区の構造を段階的に変化させ、段階毎に磁性材料の磁気エネルギーを算出して、最も磁気エネルギーが小さくなる段階を求める装置が知られている。また、従来より、磁気異方性のような、従来の解析で考慮されていない磁気的性質を取り入れて、マイクロ磁化解析を実行するプログラムが知られている。
【０００３】
マイクロ磁化解析とは、図１に示すように、磁性体を小さな磁石の集合としてモデル化し、磁区の状態を数値シミュレーションする手法である。マイクロ磁化とは、個々の小さな磁石又は小さな磁性体要素のことである。マイクロ磁化解析では、計算コストの観点から実際の原子サイズオーダーのメッシュは使用せず、約１０ｎｍのメッシュを使用する。一般に用いられるメッシュサイズ（〜１０ｎｍ）では、隣接するメッシュの磁化ベクトルの成す角は５度以下と小さく、角度はほぼ連続とみなされる。 
【０００４】
マイクロ磁化の運動を支配する方程式（支配方程式）は下記の式(1)であり、LLG（Landau-Lifshitz-Gilbert）方程式と呼ばれている。
0005593934.tif 000002
ここでM、γ、α、Heffは、それぞれ磁化ベクトル、磁気回転比、摩擦係数、有効磁界である。
【０００５】
また、有効磁界Heffは、式(2)のように複数の磁界ベクトルの合成である。マイクロ磁化が受ける磁界は、外部磁界Hout、反磁界Hdemag、異方性磁界Han、及び磁気交換結合磁界Hexである。
0005593934.tif 000003
【０００６】
磁気交換結合磁界Hexは、本来隣接する原子間に作用する力である。原子間距離より大きいサイズでメッシュ分割された解析モデルを用いて、計算精度を維持したマイクロ磁化の解析を行うためには、隣接する磁化ベクトルの角度変化が１０度以下となるぐらい小さいサイズ（〜１０ｎｍ）でメッシュ分割された解析モデルが必要である。
【０００７】
マイクロ磁化解析は、主にミクロンオーダーの小さな磁性体を解析対象としてきたが、計算機の進歩によって近年では数十ミクロンオーダーの磁性体材料への適用が実現されている。今後は、モータや変圧器等のさらに大規模な磁性体領域を解析対象とする必要性が出てくると予想されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００５−４３３４０号公報
【０００９】
【特許文献２】特開２００４−３４７４８２号公報
【００１０】
【特許文献３】特開２００４−２１９１７８号公報
【００１１】
【特許文献４】特開２００５−１０００６７号公報
【発明の概要】
【発明が解決しようとする課題】
【００１２】
上述したように、マイクロ磁化解析を高精度で実現するために必要なメッシュサイズは、約１０ｎｍ以下であるため、解析対象領域を大きくするにしたがって取り扱わなければならないメッシュ数が増大するという問題がある。メッシュサイズを大きくすると、計算に必要な自由度が小さくなるため、計算時間の短縮化が可能となる。しかしながら、メッシュサイズを大きくすると、隣接する磁化ベクトルの成す角が１０度以上になることがあるため、安易にメッシュサイズを増加すると計算精度は保証されない。
【００１３】
また、従来のマイクロ磁化解析では、メッシュサイズが大きくなると（例えば１０ｎｍ以上）、磁化ベクトルの回転角度（隣接する磁化ベクトルの成す角）が増大するため、高精度なマイクロ磁化解析を実現できないという問題、即ち、磁気交換結合エネルギーや磁気交換結合磁界の計算精度が低下する問題がある。
【００１４】
上記課題に鑑み、明細書に開示された磁気交換結合エネルギー算出プログラム、算出方法及び算出装置は、従来よりも磁化ベクトルの回転角度が大きくなる場合でも、高精度な磁気交換結合磁界を算出できることを目的とする。
【課題を解決するための手段】
【００１５】
上記目的を達成するため、明細書に開示された磁気交換結合エネルギー算出プログラムは、コンピュータに、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する第１手順、及び前記補間された磁化ベクトルに加わる力として作用する磁界を前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合磁界を算出する第２手順を実行させる。
【００１６】
明細書に開示された磁気交換結合エネルギー算出プログラムは、コンピュータに、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する第１手順、及び前記補間された磁化ベクトルが占める体積のエネルギーを前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合エネルギーを算出する第２手順を実行させる。
【００１７】
明細書に開示された磁気交換結合エネルギー算出方法は、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する第１手順、及び前記補間された磁化ベクトルに加わる力として作用する磁界を前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合磁界を算出する第２手順を備える。
【００１８】
明細書に開示された磁気交換結合エネルギー算出方法は、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する第１手順、及び前記補間された磁化ベクトルが占める体積のエネルギーを前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合エネルギーを算出する第２手順を備える。
【００１９】
明細書に開示された磁気交換結合エネルギー算出装置は、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間を結ぶ直線を複数の領域に分割し、前記分割された領域毎の磁化ベクトルを、当該２つの磁化ベクトルに直交する回転軸を基準として、当該２つの磁化ベクトルに基づき補間する補間手段と、前記補間された磁化ベクトルに加わる力として作用する磁界を前記補間された領域毎の磁化ベクトルの回転角で積分することにより、磁気交換結合磁界を算出する算出手段とを備える。
【発明の効果】
【００２０】
明細書に開示された磁気交換結合エネルギー算出プログラム、算出方法及び算出装置は、従来よりも磁化ベクトルの回転角度が大きくなる場合でも、高精度な磁気交換結合磁界を算出できる。
【図面の簡単な説明】
【００２１】
【図１】マイクロ磁化解析で使用される磁性体モデルの一例を示す図である。
【図２】磁気交換結合エネルギー算出装置の構成を示すブロック図である。
【図３】２つの磁化ベクトルの成す角が９０度である場合の例を示す図である。
【図４】式(4)を用いて算出された磁気交換結合エネルギーと式(5)を用いて算出された磁気交換結合エネルギーの誤差と、磁化ベクトルの回転角度との関係を示すグラフである。
【図５】有限体積法において磁化ベクトルを配置した様子を示す図である。
【図６】隣接する２つの磁化ベクトルが占める体積を示す図である。
【図７】２つの磁化ベクトルに直交する回転軸uと回転角Θの例を示す図である。
【図８】２つの磁化ベクトルの成す角を補間し、２つの磁化ベクトルの間をN分割する状態を示す図である。
【図９】磁気交換結合エネルギーの計算手順を示すフローチャートである。
【図１０】磁気交換結合磁界の計算手順を示すフローチャートである。
【図１１】磁区状態の解答例を示す図である。
【図１２】（Ａ）〜（Ｃ）は、磁性体モデルを１０ｎｍのメッシュサイズで分割して、マイクロ磁化解析を実行した結果を示す図である。
【図１３】（Ａ）〜（Ｃ）は、磁性体モデルを３０ｎｍのメッシュサイズで分割して、マイクロ磁化解析を実行した結果を示す図である。
【発明を実施するための形態】
【００２２】
以下、図面を参照しながら本発明の実施の形態を説明する。
【００２３】
図２は、本実施の形態に係る磁気交換結合エネルギー算出装置の構成を示すブロック図である。
【００２４】
図２の磁気交換結合エネルギー算出装置は、コンピュータ１である。コンピュータ１は、コンピュータ１の全体動作を制御するＣＰＵ（Central Processing Unit）１１（補間手段、算出手段）、ワーキングエリアとして機能するメモリ１２、ＯＳ（Operating System）やシミュレーション用のプログラムを含むハードディスクドライブ（ＨＤＤ：Hard Disk Drive）１３、ネットワークインターフェイスカード（ＮＩＣ：Network Interface Card）１４、入力インターフェイス（Ｉ／Ｆ）１５及びビデオインターフェイス（Ｉ／Ｆ）１６を備えている。ＣＰＵ１１は、メモリ１２、ＨＤＤ１３、ＮＩＣ１４、入力Ｉ／Ｆ１５及びビデオＩ／Ｆ１６に接続されている。入力Ｉ／Ｆ１５にはキーボード・マウス１７が接続されている。ビデオＩ／Ｆ１６には、モニタ１８が接続されている。ＣＰＵ１１は、ＨＤＤ１３に格納されているシミュレーション用のプログラムを読み出し、実行することで、有限体積法を利用するマイクロ磁化解析を実行する。また、マイクロ磁化解析に使用される磁性体モデルのデータは、予めＨＤＤ１３に格納されているものとする。
【００２５】
以下、２つの磁化ベクトルを用いた簡単な計算によって、磁気交換結合エネルギーの計算精度の低下問題について説明する。
【００２６】
シミュレーション用のプログラムが上記式(1)の磁化ベクトルのＬＬＧ方程式を離散化すると、磁気交換結合エネルギーは下記式(3)となる。
0005593934.tif 000004
ここで、各パラメータは以下のとおりである。
ｅ：磁気交換結合エネルギー密度、Ａ：磁気交換結合係数、Ｍｓ：磁性材料の単位体積当たりの飽和磁化、Ｌ：隣接する磁化ベクトル同士の距離、Ｍｉ：ｉ番目の磁化ベクトル、dS：隣接するメッシュ間の断面積、Θ：隣接する磁化ベクトルのなす角
【００２７】
図３に示すように、２つの磁化ベクトルの成す角が９０度である場合、シミュレーション用のプログラムが、２つのメッシュ間で磁化ベクトルが不連続となる従来手法を用いて、上記式(1)の磁化ベクトルのＬＬＧ方程式を離散化すると、磁気交換結合エネルギーは下記式(4)になる。
0005593934.tif 000005
【００２８】
同様に２つの磁化ベクトルの成す角が９０度である場合、シミュレーション用のプログラムが、２つのメッシュ間で補間した磁化ベクトルを積分して磁気交換結合エネルギーを計算すると、磁気交換結合エネルギーは式(5)になる。この計算式は、上記２つの磁化ベクトルの回転角が線形かつ連続的に変化すると仮定し、２つの磁化ベクトルが占める領域を角度θで積分することで解析的に得られる。
0005593934.tif 000006
【００２９】
前者の磁気交換結合エネルギーは、後者の磁気交換結合エネルギーに比べて磁気交換結合エネルギーの絶対値が約２０％低下する。つまり、磁化ベクトルの回転角が例えば９０度である場合のように磁化ベクトルの回転角が大きい場合、従来の手法では、精度の良い磁気交換結合エネルギーが得られない。
【００３０】
図４は、式(4)を用いて算出された磁気交換結合エネルギーと式(5)を用いて算出された磁気交換結合エネルギーの誤差と、磁化ベクトルの回転角度との関係を示すグラフである。図４に示すように、誤差の大きさは、磁化ベクトルの回転角度の増加に従って増加する。つまり、メッシュサイズが大きくなると（例えば１０ｎｍ以上）、磁化ベクトルの回転角度が増大するため、磁気交換結合エネルギーの誤差も増大する。
【００３１】
本実施の形態では、従来のメッシュに比べて、大きいサイズ（例えば３０ｎｍ）のメッシュを含む磁性体モデルを作成した場合でも、シミュレーション用のプログラムは、隣接するメッシュ間の磁化ベクトルの回転角度が連続かつ線形的に変化すると仮定して、２つの磁化ベクトルが占める領域を積分（解析解）することで、精度の良い磁気交換結合エネルギーを算出する。
【００３２】
以下、精度の良い磁気交換結合エネルギーを算出する方法を具体的に説明する。
【００３３】
一般的なマイクロ磁化解析の離散化手法には、要素中心に磁化ベクトルを配置する有限体積法が用いられる。有限体積法には有限要素法のような補間関数が存在しないが、要素中心の値を線形補間することで、要素境界面の値を求めることができる。図５は、有限体積法において磁化ベクトルを配置した様子を示す図である。以下の説明では、隣接するメッシュの磁化ベクトルの回転角が連続的に変化すると仮定し、シミュレーション用のプログラムが、有限体積法に線形補間の考え方を適用して、磁気交換結合エネルギーを高精度に求める。
【００３４】
隣接する２つの磁化ベクトルが占める体積Vijは、図６の斜線の領域である。この領域の磁気交換結合エネルギーは、下記式(6)になる。
0005593934.tif 000007
ここでdSijは２つのメッシュ間の断面積であり、 Lijは２つの磁化ベクトル間の距離である。dVijは下記式(7)で示され、Lijは下記式(8)で示される。
0005593934.tif 000008
【００３５】
従来より一般的に用いられている２つの磁化ベクトルMiとMjの磁気交換結合エネルギーの式(9)は、隣接する磁化ベクトルの回転角度が小さい場合（例えば、１０度以下の場合）には使用できるが、隣接する磁化ベクトルの回転角度が大きくなると磁気交換結合エネルギーの誤差が増大する問題を持っている。
0005593934.tif 000009
【００３６】
磁気交換結合エネルギーを正確に計算するためには、シミュレーション用のプログラムが、２つの磁化ベクトルMiとMjの成す角を算出し、線形補間された回転角を使って、２つの磁化ベクトルMiとMjが占める体積のエネルギーを積分して、解析的に磁気交換結合エネルギーを求める必要がある。回転角は線形的に変化するため、２つの磁化ベクトルの成す角θがπ以下である場合に、図７に示す回転軸u及び回転角Θはそれぞれ下記式(10)、式(11)のように計算できる。ここで、Msは磁化ベクトルの大きさである。また、回転軸uは、２つの磁化ベクトルMi、Mj と直交するベクトルである。
0005593934.tif 000010
【００３７】
磁気交換結合エネルギーの１次元における距離の積分を回転角の積分に変数変換すると、式(6)は式(12)のように変形できる。
0005593934.tif 000011
【００３８】
ここで、シミュレーション用のプログラムが、右辺の角度θに関する積分を解析的に計算すると、正確な磁気交換結合エネルギーは下記式(13)となる。
0005593934.tif 000012
【００３９】
式(13)は、シミュレーション用のプログラムが、磁気交換結合エネルギーを解析的に計算した結果であるが、シミュレーション用のプログラムが、図８に示すように、磁化ベクトルMiとMj の成す角を補間し、式(6)の積分区間をN分割し、磁化ベクトルMiとMj の占める体積のエネルギーを数値積分しても、精度の良い磁気交換結合エネルギーを計算することができる。式(6)は数値積分の形式で表現されると式(14)となる。ここで、M1はMiに相当し、MNはMjに相当する。また、wはガウス積分の重み関数である。Ｎは積分区間の分割数であり、ユーザが指定することができる。
0005593934.tif 000013
【００４０】
図９は、シミュレーション用のプログラムが実行する磁気交換結合エネルギーの計算手順を示すフローチャートである。
【００４１】
まず、シミュレーション用のプログラムは、磁気交換結合エネルギーを初期化する、即ち磁気交換結合エネルギーEexを０に設定する（ステップＳ１）。次に、シミュレーション用のプログラムは、磁性体モデルに含まれる、計算対象の磁化ベクトルMiの「ｉ」の値を「１」から「全メッシュ数」まで順次インクリメントする（ステップＳ２−１〜ステップＳ２−２のループ）。「ｉ」の値をインクリメントするタイミングは、後述するステップＳ４及びＳ５の処理が終了するＳ２−２のタイミングである。シミュレーション用のプログラムが、磁化ベクトルMiの「ｉ」の値を「全メッシュ数」に設定し、ステップＳ４及びＳ５の処理が終了した場合には、本計算は終了する。
【００４２】
シミュレーション用のプログラムは、磁性体モデルに含まれる、計算対象の磁化ベクトルMiに隣接する磁化ベクトルMjの「ｊ」の値を「０」から「全メッシュ数−１」まで順次インクリメントする（ステップＳ３−１〜ステップＳ３−２のループ）。「ｊ」の値をインクリメントするタイミングは、後述するステップＳ４及びＳ５の処理が終了するＳ３−２のタイミングである。
【００４３】
シミュレーション用のプログラムは、ステップＳ３−１の後に、上記式(7)、(8)に従って、２つの磁化ベクトルMi，Mjが占める体積dVij及び２つの磁化ベクトルMi，Mj間の距離Lijを計算し、且つ上記式(11)に従って、回転角Θを計算する（ステップＳ４）。その後、シミュレーション用のプログラムは、計算された体積dVij、距離Lij、回転角Θ及び上記式(13)を使って、磁気交換結合エネルギーEijを計算し、計算された磁気交換結合エネルギーEijをそれまでの合計の磁気交換結合エネルギーEexに加算し、全体の磁気交換結合エネルギーEexを求める（ステップＳ５）。ステップＳ５では、シミュレーション用のプログラムは、磁化ベクトルMi，Mj の「ｉ」「ｊ」の値がインクリメントされる度に磁気交換結合エネルギーEijを計算する。尚、シミュレーション用のプログラムは、上記ステップＳ５において、式(14)を使って、磁気交換結合エネルギーEijを計算してもよい。
【００４４】
シミュレーション用のプログラムは、磁気交換結合エネルギーEijを磁化ベクトルで微分することで、磁気交換結合磁界を算出することができる（Hex,ij=−∂Eij／∂Mi）。
【００４５】
以上説明したように、シミュレーション用のプログラムは、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルMi，Mj の間の回転角を、当該２つの磁化ベクトルMi，Mjに直交する回転軸uを基準として、線形的に補間し、２つの磁化ベクトルMi，Mjが占める体積Vのエネルギーを線形的に補間された回転角で積分することにより、磁気交換結合エネルギーを算出する。よって、従来よりも精度の良い磁気交換結合エネルギーを算出することができる。
【００４６】
または、シミュレーション用のプログラムは、２つの磁化ベクトルMi，Mjの間の回転角を補間し、２つの磁化ベクトルMi，Mj間を複数の領域に分割して、２つの磁化ベクトルMi，Mjが占める体積Vのエネルギーを分割された領域毎に数値積分することによって、磁気交換結合エネルギーを算出するので、従来よりも精度の良い磁気交換結合エネルギーを算出することができる。
【００４７】
また、シミュレーション用のプログラムは、上記算出された磁気交換結合エネルギーを磁化ベクトルで微分することで得られる磁気交換結合磁界をLLG方程式の有効磁界に含まれる磁気交換結合磁界として用いることで、高精度なマイクロ磁化解析を実現できる。
【００４８】
次に、精度の良い磁気交換結合磁界を算出する方法を具体的に説明する。
【００４９】
マイクロ磁化のＬＬＧ方程式において、マイクロ磁化に加わる力として作用する磁気交換結合による磁界は式(15)になる。
0005593934.tif 000014
【００５０】
従来手法では、dlijをLと近似する計算を使って、式(16)のように磁気交換結合磁界を計算した。
0005593934.tif 000015
【００５１】
この近似計算は、磁化ベクトルの回転角が大きくなるにしたがって、計算精度が低下する問題を持っている。より高精度に磁気交換結合による磁界を計算するためには、シミュレーション用のプログラムが、解析的に積分を実行する必要がある。
【００５２】
式(15)の解析的な計算では、２つの磁化ベクトルMi，Mjの回転軸に対して回転する角度θで積分する。２つの磁化ベクトルMi，Mjに直交する回転軸uは、ベクトル積によって下記式(17)で表される。回転軸uの周りに磁化ベクトルMiを回転することによって、磁化ベクトルMjは回転角の関数として表現できるので、回転角による解析的な積分が可能となる。磁化ベクトルMiを回転軸uに対してθ回転する変換行列T（Mi→Mj）は下記式(18)となり、磁化ベクトルMj の変換式は下記式(19)となる。ここで、Iijは単位行列である。
0005593934.tif 000016
【００５３】
磁化ベクトルの変化量は変換行列Tを用いて表すと、下記式(20)になる。
0005593934.tif 000017
【００５４】
シミュレーション用のプログラムが、式(20)を用いて、式(15)の磁気交換結合磁界を解析的に積分すると、下記式(21)を得る。ここでは、シミュレーション用のプログラムは、式(15)の磁気交換結合磁界の１次元における距離の積分を回転角Θの積分に変数変換している。
0005593934.tif 000018
ここで、τijは下記式(22)である。
0005593934.tif 000019
【００５５】
シミュレーション用のプログラムは、式(21)、(22)を用いることで、従来手法の式(16)に比べて、より高精度に磁気交換結合磁界を計算することが可能となる。
【００５６】
式(21)、(22)は、解析的に磁気交換結合磁界を積分した結果であるが、シミュレーション用のプログラムは、磁化ベクトルMi，Mjの成す角を補間し、式(15)の積分区間をN分割し、マイクロ磁化に加わる力として作用する磁界を数値積分しても高精度の磁気交換結合磁界を計算することができる。式(15)を数値積分で表現すると、式(23)が得られる。ここで、M1はMiに相当し、MNはMjに相当する。また、wはガウス積分の重み関数である。尚、積分区間の分割数Ｎは、ユーザが指定することができる。
0005593934.tif 000020
【００５７】
シミュレーション用のプログラムは、LLG方程式の有効磁界である磁気交換結合磁界として、解析解である式(21)の値又は数値積分より得られる式(23)の値を用いることで、高精度なマイクロ磁化解析を実現できる。
【００５８】
図１０は、シミュレーション用のプログラムが実行する磁気交換結合磁界の計算手順を示すフローチャートである。
【００５９】
まず、シミュレーション用のプログラムは、磁気交換結合磁界を初期化する、即ち磁気交換結合磁界Hex,ijを０に設定する（ステップＳ１１）。次に、シミュレーション用のプログラムは、磁性体モデルに含まれる、計算対象の磁化ベクトルMiの「ｉ」の値を「１」から「全メッシュ数」まで順次インクリメントする（ステップＳ１２−１〜ステップＳ１２−２のループ）。「ｉ」の値をインクリメントするタイミングは、後述するステップＳ１４及びＳ１５の処理が終了するＳ１２−２のタイミングである。シミュレーション用のプログラムが、磁化ベクトルMiの「ｉ」の値を「全メッシュ数」に設定し、ステップＳ１４及びＳ１５の処理が終了した場合には、本計算は終了する。
【００６０】
シミュレーション用のプログラムは、磁性体モデルに含まれる、計算対象の磁化ベクトルMiに隣接する磁化ベクトルMj の「ｊ」の値を「０」から「全メッシュ数−１」まで順次インクリメントする（ステップＳ１３−１〜ステップＳ１３−２のループ）。「ｊ」の値をインクリメントするタイミングは、後述するステップＳ１４及びＳ１５の処理が終了するＳ１３−２のタイミングである。
【００６１】
シミュレーション用のプログラムは、ステップＳ１３−１の後に、上記式(7)、(8)に従って、２つの磁化ベクトルMi，Mjが占める体積ｄVij及び２つの磁化ベクトルMi，Mj間の距離Lijを計算し、上記(17)、(18)に従って、回転軸u及び変換行列Tを計算し、且つ上記式(11)に従って、回転角Θを計算する（ステップＳ１４）。その後、シミュレーション用のプログラムは、計算された体積dVij、距離Lij、回転軸u、変換行列T、回転角Θ及び上記式(21)、(22)を使って、磁気交換結合磁界Hex,ijを計算し、計算された磁気交換結合磁界Hex,ijをそれまでの合計の磁気交換結合磁界Hexに加算し、全体の磁気交換結合磁界Hexを求める（ステップＳ１５）。ステップＳ１５では、シミュレーション用のプログラムは、磁化ベクトルMi，Mjの「ｉ」「ｊ」の値がインクリメントされる度に磁気交換結合磁界Hex,ijを計算する。尚、シミュレーション用のプログラムは、上記ステップＳ１５において、式(23)を使って、磁気交換結合磁界Hex,ijを計算してもよい。
【００６２】
以上説明したように、シミュレーション用のプログラムは、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルMi，Mjの間の回転角を、当該２つの磁化ベクトルMi，Mjに直交する回転軸uを基準として、線形的に補間し、２つの磁化ベクトルMi，Mjに加わる力として作用する磁界を線形的に補間された回転角で積分することにより、磁気交換結合磁界を算出する。よって、従来よりも精度の良い磁気交換結合磁界を算出することができる。
【００６３】
または、シミュレーション用のプログラムは、２つの磁化ベクトルMi，Mjの間の回転角を補間し、２つの磁化ベクトルMi，Mj間を複数の領域に分割して、２つの磁化ベクトルMi，Mjに加わる力として作用する磁界を分割された領域毎に数値積分することによって、磁気交換結合磁界を算出するので、従来よりも精度の良い磁気交換結合磁界を算出することができる。
【００６４】
また、シミュレーション用のプログラムは、LLG方程式の有効磁界に含まれる磁気交換結合磁界として、上記算出された磁気交換結合磁界を用いることで、高精度なマイクロ磁化解析を実現できる。
【００６５】
本実施の形態で使用されるシミュレーション用のプログラムは、マイクロ磁化を扱う磁気シミュレーション用のアプリケーションソフトである。また、本実施の形態で実行されるマイクロ磁化解析は、HDD磁気ヘッド、MRAM（Magnetoresistive Random Access Memory）、超小型モータ等への適用が可能である。
【００６６】
（実施例）
本出願人は、上記方法で算出した磁気交換結合磁界をＬＬＧ方程式の有効磁界の交換結合磁界として使用し、マイクロ磁化解析を行った。このマイクロ磁化解析では、本出願人は、下記のモデル条件（磁性体サイズ、材料特性）を使用し、磁性体モデルのメッシュサイズを変えた場合でも、メッシュサイズの変更前後の磁区状態が同様であるかを検証した。図１１は、磁区状態の解答例である。
〔モデル条件〕
磁性体領域のサイズ（X×Y×Z）：500×500×10（nm）
磁化の大きさ    Ms：8E+5（A/m）
異方性磁界      Ku：5E+2（J/m3）
磁気交換結合係数A：1.3E-11（J/m）
減衰係数        α：1
【００６７】
図１２（Ａ）〜（Ｃ）は、磁性体モデルを１０ｎｍのメッシュサイズで分割して、マイクロ磁化解析を実行した結果を示す。図１２（Ａ）は、磁化ベクトルの状態を示す。図１２（Ｂ）は磁化ベクトルのＸ軸方向の成分を示す。図１２（Ｃ）は磁化ベクトルのＹ軸方向の成分を示す。また、図１２（Ｂ）、（Ｃ）においては、赤色がプラス方向の成分を示し、青色がマイナス方向の成分を示す。磁性体モデルを１０nｍのメッシュサイズで分割した場合には、図１１に示す解答例と同様の磁区状態が得られる。
【００６８】
図１３（Ａ）〜（Ｃ）は、磁性体モデルを３０ｎｍのメッシュサイズで分割して、マイクロ磁化解析を実行した結果を示す。図１３（Ａ）は、磁化ベクトルの状態を示す。図１３（Ｂ）は磁化ベクトルのＸ軸方向の成分を示す。図１３（Ｃ）は磁化ベクトルのＹ軸方向の成分を示す。図１３（Ｂ）、（Ｃ）においては、赤色がプラス方向の成分を示し、青色がマイナス方向の成分を示す。図１３（Ａ）に示すように、磁化ベクトルの回転の不連続性が見られるが、メッシュサイズ１０ｎｍの結果と同様に、図１１の解答例と同様の磁区状態が得られる。
【００６９】
これにより、メッシュサイズを粗くしても（１０ｎｍ→３０ｎｍ）、同様の磁区状態が得られるので、上述した磁気交換結合磁界を精度良く算出する方法が有効であることが分かる。また、メッシュサイズを粗くしても（１０ｎｍ→３０ｎｍ）、精度良くマクロ磁化解析を実施できる。さらに、メッシュサイズを３倍に増加すると、ＬＬＧ方程式の自由度を約１／２７に低減でき、計算に必要なメモリリソースの低減（約１／２７）に寄与できる。また、同時に計算時間の短縮化（約１／２７）にも寄与することが可能となる。
【００７０】
以上の実施の形態に関し、更に以下の付記を開示する。
【００７１】
（付記１）コンピュータに、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間の回転角を、当該２つの磁化ベクトルに直交する回転軸を基準として、補間する第１手順、及び前記２つの磁化ベクトルに加わる力として作用する磁界を前記補間された回転角で積分することにより、磁気交換結合磁界を算出する第２手順を実行させることを特徴とする磁気交換結合エネルギー算出プログラム。
【００７２】
（付記２）前記第２手順は、前記２つの磁化ベクトルに加わる力として作用する磁界を前記線形的に補間された回転角で積分することにより得られる式（i）、（ii）に従って、前記磁気交換結合磁界を算出することを特徴とする付記１に記載の磁気交換結合エネルギー算出プログラム。
0005593934.tif 000021
Hex,ij：磁気交換結合磁界、A：磁気交換結合係数、dSij：２つの要素間の断面積、Ms：単位体積あたりの飽和磁化、Lij：２つの磁化ベクトル間の距離、Θij：線形的に補間された回転角、Iij：単位行列、ui, uj：２つの磁化ベクトルに直交する回転軸、ux、uy、uz：回転軸のベクトル成分
【００７３】
（付記３）前記２つの磁化ベクトル間を複数の領域に分割して、前記２つの磁化ベクトルに加わる力として作用する磁界を前記分割された領域毎に数値積分することによって、前記磁気交換結合磁界を算出する第３手順を備えることを特徴とする付記１又は２に記載の磁気交換結合エネルギー算出プログラム。
【００７４】
（付記４）前記第３手順は、前記２つの磁化ベクトル間を複数の領域に分割して、前記２つの磁化ベクトルに加わる力として作用する磁界を前記分割された領域毎に数値積分する式（iii）に従って、前記磁気交換結合磁界を算出することを特徴とする請求項３に記載の磁気交換結合エネルギー算出プログラム。
0005593934.tif 000022
Hex,ij：磁気交換結合磁界、A：磁気交換結合係数、dSij：２つの要素間の断面積、Ms：単位体積あたりの飽和磁化、Lij：２つの磁化ベクトル間の距離、Ｎ：分割数、wi：ガウス積分の重み関数、M1、MN：隣接する２つの磁化ベクトル
【００７５】
（付記５）前記算出された磁気交換結合磁界をLLG（Landau-Lifshitz-Gilbert）方程式の有効磁界に含まれる交換結合磁界として使用することにより、マイクロ磁化解析を行うことを特徴とする付記１乃至４のいずれかに記載の磁気交換結合エネルギー算出プログラム。
【００７６】
（付記６）コンピュータに、有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間の回転角を、当該２つの磁化ベクトルに直交する回転軸を基準として、補間する第１手順、及び前記２つの磁化ベクトルが占める体積のエネルギーを前記補間された回転角で積分することにより、磁気交換結合エネルギーを算出する第２手順を実行させる磁気交換結合エネルギー算出プログラム。
【００７７】
（付記７）前記第２手順は、前記２つの磁化ベクトルが占める体積のエネルギーを前記補間された回転角で積分することにより得られる式（iv）に従って、前記磁気交換結合エネルギーを算出することを特徴とする付記６に記載の磁気交換結合エネルギー算出プログラム。
0005593934.tif 000023
Eij：磁気交換結合エネルギー、A：磁気交換結合係数、dSij：２つの要素間の断面積、Lij：２つの磁化ベクトル間の距離、Θ：線形的に補間された回転角
【００７８】
（付記８）前記２つの磁化ベクトル間を複数の領域に分割して、前記２つの磁化ベクトルが占める体積のエネルギーを前記分割された領域毎に数値積分することによって、前記磁気交換結合エネルギーを算出する第３手順を備えることを特徴とする付記６又は７に記載の磁気交換結合エネルギー算出プログラム。
【００７９】
（付記９）前記第３手順は、前記２つの磁化ベクトル間を複数の領域に分割して、前記２つの磁化ベクトルが占める体積のエネルギーを前記分割された領域毎に数値積分する式（v）に従って、前記磁気交換結合エネルギーを算出することを特徴とする付記８に記載の磁気交換結合エネルギー算出プログラム。
0005593934.tif 000024
Eij：磁気交換結合エネルギー、A：磁気交換結合係数、dSij：２つの要素間の断面積、Ms：単位体積あたりの飽和磁化、Lij：２つの磁化ベクトル間の距離、Ｎ：分割数、wi：ガウス積分の重み関数、M1、MN：隣接する２つの磁化ベクトル
【００８０】
（付記１０）有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間の回転角を、当該２つの磁化ベクトルに直交する回転軸を基準として、補間する第１手順、及び前記２つの磁化ベクトルに加わる力として作用する磁界を前記補間された回転角で積分することにより、磁気交換結合磁界を算出する第２手順を備えることを特徴とする磁気交換結合エネルギー算出方法。
【００８１】
（付記１１）前記２つの磁化ベクトル間を複数の領域に分割して、前記２つの磁化ベクトルに加わる力として作用する磁界を前記分割された領域毎に数値積分することによって、前記磁気交換結合磁界を算出する第３手順を備えることを特徴とする付記１０に記載の磁気交換結合エネルギー算出方法。
【００８２】
（付記１２）前記算出された磁気交換結合磁界をLLG（Landau-Lifshitz-Gilbert）方程式の有効磁界に含まれる交換結合磁界として使用することにより、マイクロ磁化解析を行うことを特徴とする付記１０又は１１に記載の磁気交換結合エネルギー算出方法。
【００８３】
（付記１３）有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間の回転角を、当該２つの磁化ベクトルに直交する回転軸を基準として、補間する第１手順、及び前記２つの磁化ベクトルが占める体積のエネルギーを前記補間された回転角で積分することにより、磁気交換結合エネルギーを算出する第２手順を備えることを特徴とする磁気交換結合エネルギー算出方法。
【００８４】
（付記１４）前記２つの磁化ベクトル間を複数の領域に分割して、前記２つの磁化ベクトルが占める体積のエネルギーを前記分割された領域毎に数値積分することによって、前記磁気交換結合エネルギーを算出する第３手順を備えることを特徴とする請求項１３に記載の磁気交換結合エネルギー算出方法。
【００８５】
（付記１５）有限体積法で使用される、隣接し合う要素の中心に配置される２つの磁化ベクトルの間の回転角を、当該２つの磁化ベクトルに直交する回転軸を基準として、補間する補間手段と、前記２つの磁化ベクトルに加わる力として作用する磁界を前記補間された回転角で積分することにより、磁気交換結合磁界を算出する算出手段とを備えることを特徴とする磁気交換結合エネルギー算出装置。
【００８６】
（付記１６）前記算出手段は、前記２つの磁化ベクトルが占める体積のエネルギーを前記補間された回転角で積分することにより、磁気交換結合エネルギーを算出することを特徴とする付記１５に記載の磁気交換結合エネルギー算出装置。
【図１】
0005593934.tif 000025
【図２】
0005593934.tif 000026
【図３】
0005593934.tif 000027
【図４】
0005593934.tif 000028
【図５】
0005593934.tif 000029
【図６】
0005593934.tif 000030
【図７】
0005593934.tif 000031
【図８】
0005593934.tif 000032
【図９】
0005593934.tif 000033
【図１０】
0005593934.tif 000034
【図１１】
0005593934.tif 000035
【図１２】
0005593934.tif 000036
【図１３】
0005593934.tif 000037