よくわかる計算工学最前線 「原子間ポテンシャルの 基礎から応用まで 」

Slides:



Advertisements
Similar presentations
『わかりやすいパターン認 識』 第 5 章 特徴の評価とベイズ誤り確率 5.4 ベイズ誤り確率と最近傍決定則 発表日: 5 月 23 日(金) 発表者:時田 陽一.
Advertisements

相対論的場の理論における 散逸モードの微視的同定 斎藤陽平( KEK ) 共同研究者:藤井宏次、板倉数記、森松治.
相の安定性と相転移 ◎ 相図の特徴を熱力学的考察から説明 ◎ 以下の考察
◎ 本章  化学ポテンシャルという概念の導入   ・部分モル量という種類の性質の一つ   ・混合物の物性を記述するために,化学ポテンシャルがどのように使われるか   基本原理        平衡では,ある化学種の化学ポテンシャルはどの相でも同じ ◎ 化学  互いに反応できるものも含めて,混合物を扱う.
環境表面科学講義 村松淳司 村松淳司.
自己重力多体系の 1次元シミュレーション 物理学科4年 宇宙物理学研究室  丸山典宏.
平成18年度 構造有機化学 講義スライド 復習: 混成軌道 軌道のs性とその応用 奥野 恒久.
W e l c o m ! いい天気♪ W e l c o m ! 腹減った・・・ 暑い~ 夏だね Hey~!! 暇だ。 急げ~!!
基盤科学への招待 クラスターの不思議 2005年6月3日  横浜市立大学 国際総合科学部  基盤科学コース 野々瀬真司.
小笠原智博A*、宮永崇史A、岡崎禎子A、 匂坂康男A、永松伸一B、藤川高志B 弘前大学理工学部A 千葉大大学院自然B
スパッタ製膜における 膜厚分布の圧力依存性
学年 名列 名前 福井工業大学 工学部 環境生命化学科 原 道寛
前回の内容 結晶工学特論 第4回目 格子欠陥 ミラー指数 3次元成長 積層欠陥 転位(刃状転位、らせん転位、バーガーズベクトル)
シリカガラスの分子動力学シミュレーション
電子物性第1 第6回 ー原子の結合と結晶ー 電子物性第1スライド6-1 目次 2 はじめに 3 原子の結合と分子 4 イオン結合
TTF骨格を配位子に用いた 分子性磁性体の開発 分子科学研究所 西條 純一.
空孔の生成 反対の電荷を持つイオンとの安定な結合を切る必要がある 欠陥の生成はエンタルピーを増大させる
Ⅰ 孤立イオンの磁気的性質 1.電子の磁気モーメント 2.イオン(原子)の磁気モーメント 反磁性磁化率、Hund結合、スピン・軌道相互作用
Ⅲ 結晶中の磁性イオン 1.結晶場によるエネルギー準位の分裂 2.スピン・ハミルトニアン
Ⅳ 交換相互作用 1.モット絶縁体、ハバード・モデル 2.交換相互作用 3.共有結合性(covalency)
微粒子合成化学・講義 村松淳司
セラミックス 第4回目 5月 7日(水)  担当教員:永山 勝久.
基礎無機化学 期末試験の説明と重要点リスト
◎ 本章  化学ポテンシャルという概念の導入   ・部分モル量という種類の性質の一つ   ・混合物の物性を記述するために,化学ポテンシャルがどのように使われるか   基本原理        平衡では,ある化学種の化学ポテンシャルはどの相でも同じ ◎ 化学  互いに反応できるものも含めて,混合物を扱う.
原子核物理学 第4講 原子核の液滴模型.
分子動力学冬季講習 ポテンシャル3 共有結合性材料
非エルミート 量子力学と局在現象 羽田野 直道 D.R. Nelson (Harvard)
Dissociative Recombination of HeH+ at Large Center-of-Mass Energies
原子核物理学 第8講 核力.
シリカガラスのバルク,表面及び融着界面の構造
平成18年度 構造有機化学 講義スライド テーマ:炭素陽イオン 奥野 恒久.
今後の予定 4日目 10月22日(木) 班編成の確認 講義(2章の続き,3章) 5日目 10月29日(木) 小テスト 4日目までの内容
前回の講義で水素原子からのスペクトルは飛び飛びの「線スペクトル」
原子で書いた文字「PEACE ’91 HCRL」.白い丸はMoS2結晶上の硫黄原子.走査型トンネル顕微鏡写真.
研究課題名 研究背景・目的 有機エレクトロニクス材料物質の基礎電子物性の理解 2. 理論 3. 計算方法、プログラムの現状
原子核物理学 第2講 原子核の電荷密度分布.
Introduction to Soft Computing (第11回目)
プラズマ発光分光による銅スパッタプロセス中の原子密度評価
量子系における 確率推論の平均場理論 田中和之 東北大学大学院情報科学研究科
量子力学の復習(水素原子の波動関数) 光の吸収と放出(ラビ振動)
分子軌道理論(Molecular Orbital theory, MO理論)
(昨年度のオープンコースウェア) 10/17 組み合わせと確率 10/24 確率変数と確率分布 10/31 代表的な確率分布
最小 6.1.The [SiO4] tetrahedron
自己集積的手法による 原子膜複合体の形成と評価 法政大学 生命科学部 マイクロ・ナノテクノロジー研究センター 高井 和之
相の安定性と相転移 ◎ 相図の特徴を熱力学的考察から説明 ◎ 以下の考察
G Knebel et al J. Phys,: Condens. Matter 16 (2004)
Bruciteの(001)面における真の接触面での摩擦特性
学年   名列    名前 物理化学 第2章 2-1、2-2 Ver. 2.1 福井工業大学  原 道寛 HARA2005.
チャネル結合AMDによる sd殻Ξハイパー核の研究
第一原理計算でひも解く合金が示す長周期積層欠陥構造の形成メカニズム
原子核物理学 第5講 原子核の振動と回転.
ベイジアンネットワーク概説 Loopy Belief Propagation 茨城大学工学部 佐々木稔
タンパク質-リガンド複合体への共溶媒効果の系統的解析
経営学研究科 M1年 学籍番号 speedster
サポートベクターマシン Support Vector Machine SVM
1:Weak lensing 2:shear 3:高次展開 4:利点 5:問題点
水素の室温大量貯蔵・輸送を実現する多孔性材料の分子ダイナミクスに基づく解明と先導的デザイン
格子ゲージ理論によるダークマターの研究 ダークマター(DM)とは ダークマターの正体を探れ!
今後の予定 7日目 11月12日 レポート押印 1回目口頭報告についての説明 講義(4章~5章),班で討論
α decay of nucleus and Gamow penetration factor ~原子核のα崩壊とGamowの透過因子~
原子核物理学 第7講 殻模型.
課題研究 P4 原子核とハドロンの物理 (理論)延與 佳子 原子核理論研究室 5号館514号室(x3857)
原子核物理学 第6講 原子核の殻構造.
相の安定性と相転移 ◎ 相図の特徴を熱力学的考察から説明 ◎ 以下の考察
2008年 電気学会 全国大会 平成20年3月19日 福岡工業大学 放電基礎(1)
弱電離気体プラズマの解析(LXXVI) スプラインとHigher Order Samplingを用いた 電子エネルギー分布のサンプリング
北大MMCセミナー 第23回 Date:2014年3月6日(木) 16:30~18:00 ※通常と曜日が異なります
複合アニオンに起因した多軌道性と低次元性からうまれる 強相関電子物性の研究
学年   名列    名前 物理化学 第2章 2-1、2-2 Ver. 2.0 福井工業大学  原 道寛 HARA2005.
軽い原子核の3粒子状態 N = 11 核 一粒子エネルギー と モノポール a大阪電気通信大学 b東京工業大学
Presentation transcript:

よくわかる計算工学最前線 「原子間ポテンシャルの 基礎から応用まで 」 よくわかる計算工学最前線 「原子間ポテンシャルの 基礎から応用まで 」 2003-01-17 東京大学工学系研究科機械工学専攻 泉 聡志

日本材料学会分子動力学部門委員会 ポテンシャルWG 局所構造を反映した原子間ポテンシャル作成手法の調査 Tersoff型、 MEAM(modified EAM)型、 Tight-binding法 データベース化 分子動力学冬季講座 fcc, bcc系のポテンシャル HCP系のポテンシャル 共有結合系のポテンシャル

概要 ポテンシャルの基礎(思想) ポテンシャルの応用 二体ポテンシャル イオンポテンシャル 多体ポテンシャル(共有結合系) クラスターポテンシャル ボンドオーダーポテンシャル Tight-Binding法 ポテンシャルの応用 シリコンクラスター用Tersoffポテンシャルの作成 SiH4表面反応用ポテンシャルの作成

原子間ポテンシャルの階層 Density Functional Theory ( GGA, LDA ) CASTEP, Wien97 C.Z.Wang, D. G. Pettifor … (Order-N) Tight-Binding Tersoff, Finnis… Empirical Potential (Bond-order) Vashishita… Charge-Variable Potential Lennard-Jones, Morse Empirical 2-body Potential 近似 速度

密度汎関数法(DFT) 様々な近似(擬ポテンシャル・交換相関エネルギ)に基づき、電子状態をセルフコンシストに算出 原子個数 50~100個程度(状態の緩和+エネルギ計算) 汎用ソフトが整備されつつある VASP CASTEP

密度汎関数法(DFT)

VASPによるc-Si,SiO2(β-c)電荷密度の計算 密度汎関数法によるKS方程式の自己無撞着場計算 平面波‐擬ポテンシャル法 ウルトラソフト型擬ポテンシャル 相関交換項の計算 : GGA , LDA 固有関数の計算 : 共役勾配法 , RMM-DIIS 分子動力学 : 共役勾配法, Verlet法 c-Si(64個)の電荷密度の様子[100]面 c-Si(64個)の電荷密度の様子[110]面 SiO2(24個)の電荷密度の様子[100]面 SiO2(24個)の電荷密度の様子[110]面

分子動力学 電子状態の計算を経験的ポテンシャル関数で近似し、原子系のダイナミクス・統計集合を解く方法 ポテンシャルの精度が課題

分子動力学ポテンシャル 二体ポテンシャル イオンポテンシャル 多体ポテンシャル(共有結合) Tight-Binding法 固定電荷型 電荷移動型 多体ポテンシャル(共有結合) クラスターポテンシャル ボンドオーダーポテンシャル Tight-Binding法

1.二体ポテンシャル1 Lennard-Jonesポテンシャル 原子の電荷の偏りによる双極子間の引力相互作用(van der Waals力) 希ガス(Ar, Ne)などに適用

1.二体ポテンシャル2 Morseポテンシャル クラスタの結合のために提案 金属等に応用 Expの関数系は後にEAMに使われる(第一原理計算と良く合う~Erolessi(1994)) でφ=-B (結合エネルギ) 引力項∝h , 斥力項∝|h|2 となるWolfsberg-Helmholtz近似を満たしている。

1、二体ポテンシャルの問題点 弾性定数の記述が貧弱 空孔生成エネルギや表面などback-bondの記述ができない 必ずCauchyの関係式C12=C44が成り立つ 空孔生成エネルギや表面などback-bondの記述ができない 環境に依存したポテンシャルが必要!!

分子動力学ポテンシャル 二体ポテンシャル イオンポテンシャル 多体ポテンシャル(共有結合) Tight-Binding法 固定電荷型 電荷移動型 多体ポテンシャル(共有結合) クラスターポテンシャル ボンドオーダーポテンシャル Tight-Binding法

2. イオンポテンシャル(固定電荷) BHSポテンシャル(SiO2など) 長距離ポテンシャルのため、カットオフ(スクリーニング)もしくはEwald和が必要 Qは実効電荷、σはイオン半径、bはイオンの堅さ 長距離力の計算により、計算時間は10倍程度になる

2. Charge-Variable MD SiO2, Si3N4の電荷の移動 バルクに対しては固定電荷で良いが、界面・表面に対しては、局所的な電荷の移動を考慮する必要がある。 電荷の移動+原子間力の二つを解く必要がある。 Vashishita, 保川等 計算時間100倍程度

2, Charge-Variable MD

分子動力学ポテンシャル 二体ポテンシャル イオンポテンシャル 多体ポテンシャル(共有結合) Tight-Binding法 固定電荷型 電荷移動型 多体ポテンシャル(共有結合) クラスターポテンシャル ボンドオーダーポテンシャル Tight-Binding法

共有結合性材料の特徴 sp2,sp3結合等、方向性を持った結合 C, Siの最外殻はs, px, py, pzの4軌道 sp3混成軌道

3.Si(C)系の多体ポテンシャル Stillinger-Waberポテンシャル Tersoffポテンシャル “多体クラスターポテンシャル” ダイヤモンド構造と液体構造にフィッティング Tersoffポテンシャル “ボンドオーダーポテンシャル” 多形態(ダイマ・グラファイト・ダイヤモンド・単純立方・体心立方など)のシリコンの構造にフィッティング 発展形としてBrenner(C-H)、Marty(Si-H)、大平(Si-H)ポテンシャル

3-1. 多体クラスターポテンシャル ある特定の構造について、結合距離(二体項)、結合角(三体項)、二面角(四体項)を合わせこむ。 合わせ込んだ構造に以外に遷移する場合(化学反応・相転移等)に対応できない!

SWポテンシャルの形 式(1)~(6) 二体項はrαβのみの関数 三体項はrαβ、 rαγ、 cosθの関数 hは1/3。 θ=109.47°でh=cosθ (ダイヤモンド構造が安定) α γ β θ

SWポテンシャルの合わせこみ ダイヤモンド構造の凝集エネルギ・格子定数・融点 液体構造

SWポテンシャルの物性値1 テキスト 表5 Diamond struct. Exp. SW T3 Lattice constants 5.429[Å] 5.431 5.432 Cohesive energy 4.63[eV] 4.63 Elastic constants C11 167.0[GPa] 161.6 142.5 C12 65.0[GPa] 81.6 75.4 C44 81.0[GPa] 60.3 69.0 Vacancy energy (3-4) [eV] 2.82 3.70 Interstitial energy Td (5-6) [eV] 5.25 3.45 H (4-5) [eV] 6.95 4.61 Bond-centered 5.99 5.86 Surface energy (100) 1.57[/1×1 cell] 1.416 1.367 Surface energy (111) 1.39[/1×1 cell] 1.158 0.957

SWポテンシャルの物性値2 クラスターの形状は合わない テキスト 図3

SWポテンシャルの物性値 バルク(ダイヤモンド構造)の性質は概ね合っている ダイヤモンド構造以外の他の形態は保証なし クラスターの形状は合わない Si系のポテンシャルの評価はH. Balamaneらによって詳細に行われている。Phys. Rev. B46(1992), 2250 文献[7]

SWポテンシャルの応用 格子間シリコンの平衡濃度・拡散挙動(Brown, Maroudas 1994) 転位の移動度(Bulatov, Yip 1996) Epitaxy結晶成長(Gilmer 1992) インプランテーション(Rubia, Gilmer 1995) SiO2系への拡張(Jiang, Brown 1995)

3-2. ボンドオーダーポテンシャル 見かけは二体ポテンシャルbαβに多体項が含まれる。 配位数が大きくなれば、結合を作るための十分な価電子がなくなり、電子が非局在化して結合間で共鳴し、結合を弱める効果を表現出来る。 局所状態密度の二次モーメント近似より導かれる。

Tersoffポテンシャル ボンドオーダーに配位数依存(η,δで調整)とともに、角度依存性g(θ)を入れ、シリコンの物性値にフィッティングした。

Abellポテンシャル 現在提案されてる、EAMポテンシャル・FSポテンシャル・Tersoffポテンシャルは、すべて上の式で表される二次モーメント近似によるポテンシャルである。 →二次モーメント近似の導出(資料参照)

角度依存性・配位数依存性 120°付近でg(θ)が極小→ボンドオーダー増加 配位数が増大→ボンドオーダー減少 Bond-order g(θ) 配位数増加 g(θ) Bond-order 126° ボンドオーダの角度・配位数依存性 角度依存項

Tersoffポテンシャルの合わせこみ 異なる結晶構造(結合状態)のエネルギ・格子間距離 ダイマ→グラファイト→ダイヤモンド→単純立方→BCC→FCC 異なる二つのバージョンが存在 表面構造合わせこみ(T2)→あまり使われていない 弾性定数合わせこみ(T3)

Tersoffポテンシャルの物性値1 テキスト 表5 Diamond struct. Exp. SW T3 Lattice constants 5.429[Å] 5.431 5.432 Cohesive energy 4.63[eV] 4.63 Elastic constants C11 167.0[GPa] 161.6 142.5 C12 65.0[GPa] 81.6 75.4 C44 81.0[GPa] 60.3 69.0 Vacancy energy (3-4) [eV] 2.82 3.70 Interstitial energy Td (5-6) [eV] 5.25 3.45 H (4-5) [eV] 6.95 4.61 Bond-centered 5.99 5.86 Surface energy (100) 1.57[/1×1 cell] 1.416 1.367 Surface energy (111) 1.39[/1×1 cell] 1.158 0.957

Tersoffポテンシャルの物性値2 クラスターの形状は合わない テキスト 図3

Tersoffポテンシャルの物性値 バルク(ダイヤモンド構造)の性質は概ね合っている ダイヤモンド構造以外についてもエネルギは保証→相転移にも対応できる クラスターの形状は合わない(T2はクラスターの形状は合うが、エネルギは×)。 Si系のポテンシャルの評価はH. Balamaneらによって詳細に行われている。Phys. Rev. B46(1992), 2250 文献[7]

Tersoffポテンシャルの応用 照射欠陥(Kitabatake 1993) アモルファス化・固相成長(Motooka 1993, 2000) 圧力相変態(Mizushima, Yip 1994) インプランテーション(Nordlund 1998) アモルファス化(Bond-defect)(Marques 2001) 水素系への応用 Si-H Marty(1995), 大平(1994), 泉(2002) C-H Brenner(1990), (2002)

シリコン系へのTersoffポテンシャルの応用へ

FSとTersoffポテンシャルの類似性

金属系のボンドオーダーポテンシャル Tersoffポテンシャル FSポテンシャル EAMポテンシャル ※金属結合は不飽和な共有結合 結合角・距離(弱い)と配位数でボンドオーダーを決定 FSポテンシャル 結合距離と配位数でボンドオーダーを決定 EAMポテンシャル ただし、埋め込み関数(ボンドオーダーの効果に相当する項)を多項式やlog形にして汎用性・柔軟性を持たせる → FSと傾向は同じ  

分子動力学ポテンシャル 二体ポテンシャル イオンポテンシャル 多体ポテンシャル(共有結合) Tight-Binding法 固定電荷型 電荷移動型 多体ポテンシャル(共有結合) クラスターポテンシャル ボンドオーダーポテンシャル Tight-Binding法

4. Tight-Binding分子動力学 経験的ポテンシャルでは、バルク・クラスターなど多種多様な物性を表現するのが難しい。扱えてもパラメータが多くなる。 量子効果を最もシンプルに取り込むTightーBinding法(TB法)が用いられている TB法はマトリクスの対角化が必要になるため、計算時間はO(N2)。

Tight-Bindingの解法の基礎1 n番目の固有状態ψ(n)の波動関数を原子軌道φiα(原子i, 軌道α)の線形結合で表す 波動関数に代入して永年方程式を得る Tight-Binding法ではHiαjβの計算にSlater-Kosterの式を使う。よって、積分不要。

Tight-Bindingの解法の基礎2 ハミルトニアンマトリックスHiα,jβの対角化(固有値問題)を行い、固有値(エネルギ)を求める。Order N3 小さいエネルギ準位から電子を二つずつ詰めて行くと、バンドエネルギが得られる。

Order-N Tb法 対角化をなくしてエネルギを求める Bond Order Potential法 Density Matrix法 Fermi Operator法 Global Density of State法 BOP法は、Abell-Tersoffポテンシャルの基盤となる考え方。

計算をはじめるにあたって Tersoffポテンシャル http://www.fml.t.u-tokyo.ac.jp/~izumi/MD/  にてソースコード+マニュアル配布 Tight-Bindingポテンシャル L.Colombo, Comp. Mater. Sci, 12(1998),278 ElsevierのHP上にソースコードあり。 ポテンシャルWG(日本材料学会分子動力学部門委員会WG「局所構造を反映した原子間ポテンシャル作成手法の調査」 ) http://www.fml.t.u-tokyo.ac.jp/wg-cmd/

ポテンシャルの開発の試み

角度/距離依存性・配位数依存性でどこまで高精度なポテンシャルが出来るか? シリコンクラスタ用Tersoffポテンシャル T2,T3→T4の開発 SiH4のシリコン表面反応ポテンシャル ボンドオーダーポテンシャルに基づいた形式 ポテンシャルの合わせこみデータの収集 Saddle-Pointの合わせこみ

Tersoffポテンシャルのクラスタへの応用 T2, T3ともにクラスターの構造が記述できない T3は角度依存性が大きく、表面に向かない T2は金属的結合で角度依存性が小さい Tersoff3の結晶成長計算 表面原子が109°の角度(ダイヤモンド構造の位置)で付着

クラスタ用ポテンシャルの作成 クラスターの構造間エネルギ差 Si3~Si6までのクラスタの形状間エネルギ差をエネルギ構造差定理を使うことによりfitting。 バルク(グラファイト・ダイヤモンド・単純立方・BCC)の格子間距離・エネルギを最低条件 → ボンドオーダーのみを調節して、T3より角度依存性・配井数依存性を弱くして、T2の性質に近づける。 どの構造がより安定化?を重視

結果(クラスターエネルギ) Tersoffと比較して、エネルギの順序の再現が大幅に改善!

考察(Global Min. ) ●S1b Si4, Si5の最安定構造はAb-initioと一致する。Si6は三角柱が最安定となり、Ab-initioの八面体とは異なる。Si10もよりOpenな五角柱となり、Ab-initioと傾向が異なる。 ●S1c Si6, Si10の最安定構造の傾向があうポテンシャルを作成した。しかしながら、バルクの性質との共存は難しく。バルクの性質を犠牲にせざるを得なかった。具体的には、ダイヤモンド構造が最安定にならず、BCCが最安定となってしまった。

結論(クラスタ用ポテンシャル) ボンドオーダーを調整することによって、クラスタについては、T2・T3より優れたポテンシャルが作成できた。 しかしながら、バルクとクラスタの性質の双方を満足するポテンシャルは角度と配位数の調節だけでは不可能であった。 シリコンクラスタはダングリングボンドの数が多く、古典MDでは難しい。

実用的なCVD表面反応用Si-Hポテンシャルの開発 ポテンシャル作成における課題 どのような形状(安定or鞍点?)をあわせ込むのか? フィッティングデータをどうやって集めるのか? どのような関数形状を選ぶのか? フィッティングアルゴリズムをどうするのか?

1、どのような形状(安定or鞍点?)をあわせ込むのか? クラスタ→バルクへの成長を表現するため、様々な大きさのクラスタを合わせこみ対象に選定 汎用性を高めるために、同一の大きさでも、異なる形状のクラスタを合わせこみ対象に選定 化学反応過程を考え、表面反応の素過程に重要な現象のエネルギ障壁を合わせる

合わせ込んだクラスタ

合わせ込んだ反応過程 平衡値だけの合わせこみで良いのか? 結晶成長の表面反応をシミュレーションするためには、Saddle-Pointのエネルギを合わせる必要がある。 ※結晶成長に重要な素過程 SiH4の 解離吸着 Hの脱離

2、フィッティングデータをどうやって集めるのか? 実験値の収集は大変・実験によりバラツキ 特に0Kでの結合エネルギの値の収集は困難 実験値は用いず、すべてGAUSSIAN98のG2セットの計算結果を用いた(実験値との平均誤差1.5%) 1 HF/6-31G(d) Opt Freq 予備計算 2 MP2(full)/6-31G(d) Opt 形状最適化 3 MP4/6-311G(d,p) 以降エネルギ計算(補正) 4 MP4/6-311+G(d,p) 5 MP4/6-311G(2df,p) 6 QCISD(T)/6-311G(d,p) 7 MP2/6-311+G(3df,2p) G2の計算レシピ

3、どのような関数形状を選ぶのか? Tersoff3をベースにしたEBOPを関数形状に選定

4、フィッティングアルゴリズムをどうするのか? SiH4のエネルギ・結合距離・振動数を合わせこみ ( A, B, λA, λB, c, d ) SiH,SiH2,SiH3のエネルギ・結合距離を合わせこみ( F1Si-H, F2Si-H で補正) Si2H6, Si2H4, Si2H2, HSiSiH3のエネルギと結合間距離を合わせこみ ( GijSi-Siで補正) Si3Hx~Si7Hxのエネルギがよく一致するように(特にエネルギの順序)合わせこみ( GijSi-Siで補正) 4まで出来た後に、水素脱離とSiH4解離吸着をH-H-Si, H-Si-H, H-Si-Hの三体項を調整することによって合わせこむ 1、はcrude MCによって、最適値を見つける。2,3,4は1の中からtry and errorで合わせこむ。5は4で設定していない三体項を使って合わせこむ。

誤差は、ほとんどが2%以下、平均で1.1%(Martyは3.3%) 環状のクラスタは苦手(Tersoff3の特性) Fitting結果(誤差) 誤差は、ほとんどが2%以下、平均で1.1%(Martyは3.3%) 環状のクラスタは苦手(Tersoff3の特性)

Fitting結果 (Hの解離1) Eb=2.37eV (DFT:2.49) d1 d2 H解離のポテンシャルマップ Si-Hのカットオフ距離と H-H-Si, H-Si-Hの三体項

Fitting結果 (SiH4の解離吸着) -0.70eV Eb=+0.56eV (ab-inito 0.5~0.6eV) Saddle-point! Eb Si-Siのカットオフ距離と H-Si-Hの三体項

結論 どのような形状(安定or鞍点?)をあわせ込むのか? →クラスタ形状と表面反応の活性化エネルギ SiH4の表面反応に使用する分子動力学ポテンシャルを作成 どのような形状(安定or鞍点?)をあわせ込むのか? →クラスタ形状と表面反応の活性化エネルギ フィッティングデータをどうやって集めるのか? →GAUSSIAN98による計算データ どのような関数形状を選ぶのか? →EBOPに補正項を入れる。パラメータが多ければあわせこみの可能性も増える フィッティングアルゴリズムをどうするのか? →Crude MC+try and error

ご清聴ありがとうございました

注意!(Cutoff関数の問題点) x 新しい結合が出来ると、急激にエネルギ低下=z局所的に強い引力が働く!! H Si 4.0Å Si-Hのcutoff 2.4Å Si H 新しい結合が出来ると、急激にエネルギ低下=z局所的に強い引力が働く!! 結合を切るときは逆に局所的に大きな障壁 Cutoffにより、短距離の間に大きなエネルギ変化が生じる!!

Saddle-point energy計算 双方向Gradual Ascent法(GA法) Bulatovらが開発したGA法(初期状態と最終状態をバネで結び、剛性を除々に強くすることによって、saddle-pointを見つける)を双方向に拡張 ② ④ 初期状態 GA 1 GA 3 GA 4 ① GA 2 ③ ⑤ GAバネ 最終状態 初期状態・最終状態ともに、少しずつGAで近づけてsaddle-pointを探す

Fitting手順 (Hの解離2) Eb=2.9eV (DFT:2.9) d1 d2 H解離のポテンシャルマップ Si-Hのカットオフ距離と H-H-Si, H-Si-Hの三体項

Tight-Bindingの解法 n番目の固有状態ψ(n)の波動関数を原子軌道φiα(原子i,軌道α)の線形結合で表す 波動関数に代入して永年方程式を得る ハミルトニアンマトリックスHiα,jβの対角化(固有値問題)を行い、固有値(エネルギ)を求める。

状態密度と結合エネルギー Energy 原子 固体 エネルギの状態密度分布 ntotal(E) 結合は、原子状態より低いエネルギの分子軌道が作られることによって生じる。 → 状態密度の分布形状が結合エネルギを決める。

状態密度のモーメント E n(E) 二次モーメント近似(軌道数1) E n(E) W/2 -W/2 W/2 -W/2 1/W しかし、これではLOCALな結合の記述が出来ない

局所状態密度の定義 個々の分子軌道からの全体の結合エネルギへの寄与

局所状態密度と結合形態の関係 モーメント定理 局所状態密度のモーメントは、該当原子から出発し、戻ってくるすべてのパスのN回の“HOP”(ハミルトニアンマトリックスの積)の和によって決まる。 4-Hops 四次モーメント 2-Hops 二次モーメント

二次モーメント近似 局所状態密度の二次モーメントは近接原子数によってのみ決まる。よって、結合エネルギは以下の式のように配位数の1/2乗となる。→二次モーメント近似 Z Ebind

局所構造を反映した 原子間ポテンシャル作成手法の調査       主旨   材料を構成する原子構造体の外力に対する応答として,き裂先端場から射出される転位挙動や結晶粒界近傍に生じる種々の力学現象などが,分子動力学法により正確に描写できるようになった.しかしながら,非可逆な変形応答やその極限状態 としての破壊事象は,本質的に不均質な構造を対象とした現象であり,これまでの平均的で均質系を仮定したマクロ系の材料特性のみを表現しうる従来の原子間相互作用の記述では,本質的に捉えられない. そこで,本WGでは従来からあるポテンシャルの問題点などを整理するとともに,不均 質な局所構造の影響を反映した原子間ポテンシャル作成のための調査を目的とする.

活動内容 現状の経験的,半経験的ポテンシャルの問題点の整理     Tersoff型、 MEAM(modified EAM)型、 Tight-binding法 非平衡現象に対応するポテンシャルのフィッテング手法の開発 汎用第一原理MDの利用技術の確立 ポテンシャルの第一原理MDへの繰み込み手法 ポテンシャルのデータベース構築の検討 データベースの作成 大学/企業で必要と考えられるポテンシャルの調査

活動方法 1. 年4回程度のミーティングを実施し,あらかじめ決められた報告者が,事前にE-mailで議事内容を連絡するとともに当日会議の進行を行う. 2. 報告の内容は,設定された調査の経過報告や文献紹介など,特に制約は設けない. 3. 学会のオーガナイズドセッションやワークショップ,フォーラムなどを利用した 活動結果の公表を定期的に行う.