研究紹介
計算理論・各種資料
こちらのページをご覧ください.
材料組織シミュレーション
磁気記録材料の組織形成
Fe-45at%Pt合金の923KにおけるFePtナノグラニュラー組織形成の計算結果です．島状の部分がFePt相で，黒色と黄色がそれぞれ不規則状態と規則状態を表しています．FePt相のサイズの増大に伴い規則化が進行する様子が再現されています．
T. Koyama, H. Onodera, Materials Transactions, 44, 1523-1528 (2003).
強誘電体の分極ドメイン組織形成
強誘電体BaTiO3の298Kにおける構造相転移の2次元計算結果です．黒い部分が立方晶，色のついた部分が正方晶です．正方晶のc軸方向が分極モーメントの方向に対応します．構造相転移における弾性拘束や分極モーメント間の双極子－双極子相互作用の影響を受けて分極ドメイン組織が形成する様子が再現されています．
T. Koyama, H. Onodera, Materials Transactions, 50, 970-976 (2009).
原子力燃料棒の界面反応相の安定性
1次元計算により，原子力燃料棒のUO2 (Fluorite)/β-Zr (Bcc)界面における反応相の安定性を解析した結果です．温度は1873Kとし，反応相として液相(L)とα-Zr (Hcp)相を想定しています．解析には平衡状態図の熱力学データベースや各元素の拡散係数を用いています．現在，この手法を昇温過程に適用することで，燃料棒の液相化温度の推定を試みています．
Y. Nishida, Y. Tsukada, T. Koyama, M. Kurata, Journal of Nuclear Materials, 466, 551-559 (2015).
Mg合金の長周期積層構造
Mg-Y-Zn系の573Kにおける組織形成の2次元計算結果です．hcp相の二相分離組織形成（図a）やhcp相中の積層欠陥にYやZnが偏析した場合の組織形成（図b）を解析し，この系で観察されている長周期積層構造の形成メカニズムの解明を試みています．
K. Narita, T. Koyama, Y. Tsukada, Materials Transactions, 54, 661-667 (2013).
T. Koyama, Y. Tsukada, K. Narita, Materials Transactions, 56, 937-942 (2015).
材料外部領域を熱力学的矛盾なく扱えるCALPHAD連携フェーズフィールドモデルの開発
焼結など、材料の相とその外部領域が共存する多相・多成分材料プロセスにおけるフェーズフィールド（PF）モデリングには，熱力学的に一貫した定式化が欠けています．この問題に対処するため，本研究では材料外部領域（PED）を備えた新しい多相フィールド（MPF）モデルを開発しました．本MPF-PEDモデルは，外部領域の取り扱いにおける従来モデルの限界を克服し，複雑な材料プロセスの現実的なシミュレーションを可能にする，堅牢かつ物理的に一貫性のある枠組みを提供する手法です．
A. Ishii, Y. Matsuoka, T. Koyama, STAM: Methods. 6 [1] (2026) 2653418.
MgO上にスパッタ成膜されたFePt-Cナノ粒子膜組織形成のフェーズフィールドモデリング
本研究では，スパッタリング法を用いてMgO基板上に成膜したFePt-Cナノ粒子膜の微細組織形成を、フェーズフィールド法を用いてモデル化し，ハードディスクドライブに求められる高記録密度に不可欠な「島状」微細構造の形成を促進する主要因を特定することに成功しました．
Y. Matsuoka, M. Ode, T. Abe, T. Koyama, Y. K. Takahashi, Materials & Design. 261 (2026) 115314.
材料特性計算
組織形態情報を用いた磁気ヒステリシス計算
フェーズフィールド法により得たFe-45at%Pt合金のFePtナノグラニュラー組織をベースに，外部磁場を80kOeから-80kOeまで変化させたときの磁気モーメント場の変化を計算した結果です（左図）．磁気モーメントが右方向（赤色）から左方向（青色）へ反転する過程が再現されています．また，外部磁場方向の磁気モーメント成分を積分することにより，磁化過程における磁気ヒステリシスを得ることができます（右図）．
小山敏幸, まぐね, 3, 564-571 (2008).
分極ドメイン組織変化に対応した分極ヒステリシス計算
強誘電体BaTiO3に交番電場を作用させたときの分極ドメイン組織変化の2次元計算結果です（左図）．正方晶のc軸方向（分極モーメントの方向）が異なる領域を色で区別しています．外部電場と分極ドメイン間に作用する双極子－双極子相互作用の影響で，赤色のドメインが縮小する様子が分かります．また，外部電場方向の分極モーメント成分を積分することにより，組織変化に対応した分極ヒステリシスを得ることができます（右図）．
T. Koyama, H. Onodera, Materials Transactions, 50, 970-976 (2009).
改良セカント法に基づく応力－ひずみ曲線計算
単相の応力－ひずみ曲線が既知である前提のもとに，二相組織の応力－ひずみ曲線を算出した結果です．まず，フェーズフィールド法を用いて形態が異なる組織を用意し，その組織をベースに応力－ひずみ曲線を計算しています．左右の計算結果を比較すると，析出相の形態が応力－ひずみ曲線に影響を及ぼしていることが分かります．
T. Koyama, ISIJ International, 52, 723-728 (2012).
AIエージェント
自由エネルギー計算およびフェーズフィールドシミュレーションにおけるAIエージェントの簡易実装例
ギブズエネルギーの計算や拡散シミュレーションを例に，シンプルなAIエージェントの実装方法とその有用性について解説した研究で，サンプルのPythonコードも補足資料として配布しています．またジャルジンスキー等式の拡散対への適用は世界初の試みです．
T. Koyama, Y. Matsuoka, A. Ishii, STAM: Methods. 6 [1] (2026) 2601436.
