目次

第一原理計算実習

未整理計算結果

第一原理計算は強力な1)道具ではあるが、決して万能ではない。常に、次のようなことに留意せよ。

第一原理計算の注意点

Lesson 1: 全エネルギー計算

最も簡単な計算の例として、一酸化炭素(CO)の一点計算を行う→実習 CO in a box

Lesson 2: 断熱ポテンシャル面

- 一酸化炭素(CO)の断熱ポテンシャル 時刻<latex>{t}</latex>における炭素の位置を${x(t)}$、酸素の位置を${y(t)}$、酸素と炭素がバネでつながっていると考えてそのバネ定数を${K}$とする。ただし、酸素と炭素の位置はバネが伸びも縮みもしていない状態すなわち力が働いていない状態からのずれで表わすことにする。このとき、全エネルギーは $$ E = \frac{1}{2}M\left( {\frac{{{\rm{d}}x}}{{{\rm{d}}t}}} \right)^2 + \frac{1}{2}m\left( {\frac{{{\rm{d}}y}}{{{\rm{d}}t}}} \right)^2 + \frac{1}{2}K\left( {x\left( t \right) - y\left( t \right)} \right)^2 $$ で与えられる。 $$ M\displaystyle\frac{{{\rm{d}}^2 }}{{{\rm{d}}t^2 }}x\left( t \right) = - K\left( {x\left( t \right) - y\left( t \right)} \right) \\ $$

$$ m\displaystyle\frac{{{\rm{d}}^2 }}{{{\rm{d}}t^2 }}y\left( t \right) = - K\left( {y\left( t \right) - x\left( t \right)} \right) \\ $$

$$ x\left( t \right) = x_0 {\rm{e}}^{ - i\omega t} $$

$$ y\left( t \right) = y_0 {\rm{e}}^{ - i\omega t} $$

$$ - M\omega ^2 x_0 {\rm{e}}^{ - i\omega t} = - K\left( {x_0 {\rm{e}}^{ - i\omega t} - y_0 {\rm{e}}^{ - i\omega t} } \right) $$

$$ - m\omega ^2 y_0 {\rm{e}}^{ - i\omega t} = - K\left( {y_0 {\rm{e}}^{ - i\omega t} - x_0 {\rm{e}}^{ - i\omega t} } \right) $$

演習問題

  1. CO2の伸縮振動の振動数を求めよ。→演習
  2. COとCO2でバネ定数に相当するものを比較せよ。→演習

Lesson 3: 構造最適化

 - 一酸化炭素(CO)の構造最適化

演習問題

  1. H2OのO-H間距離とH-O-Hの角度を求めよ。→演習
  2. SiH4のSi-H間距離を求めよ。→演習

Lesson 4: 結晶

演習問題

  1. シリコン結晶の安定な格子常数を求めよ。→演習
  2. 面心立方構造の鉄結晶について磁気的性質を調べよ。→演習

Lesson 5: 結晶表面

  1. Al(001)面
  2. Al(111)面
  3. Al(110)面

Lesson 6: 表面緩和

Lesson 7: 表面吸着

Lesson 8: 電子の分散関係(バンド計算)

1)
強力すぎるほど強力な
2)
吸着エネルギーや断熱ポテンシャルなど
3)
不純物・吸着子・欠陥などに由来する局在状態など
4)
バルクのエネルギーギャップやエネルギーポケット中に存在する表面状態の分散関係など