====== 一酸化炭素の計算 ======
===== 実習1 =====
  * 結合エネルギーをエレクトロンボルト単位で求める。
    - [[seminar:remote#coteditを使ってmac側にcopyを作成|前回の実習]]を参照し、分子のエネルギーをJacapoで計算する。
    - （結合エネルギー）＝（分子のエネルギー）ー（原子のエネルギーの和）\\ となることを用いて結合エネルギーをもとめる。

^炭素原子のエネルギー| $-$10.746116094 Ry|$-$146.20827010491 eV |
^酸素原子のエネルギー| $-$31.634296294 Ry|$-$430.4062692673 eV |
 (1 Ry(リュドベリー) = 13.605684958731 eV（エレクトロンボルト）) 

==== 分子のエネルギー ====
<code text>
$ grep TOT: out.txt
 TOT:  CPU time                  Total energy
 TOT:                   LDA            GGA            LDA
 TOT:             PerdewZunger   PerdewWang91      VosWilNus
 TOT:  CPU time   non-selfcons       selfcons   non-selfcons
 TOT:   seconds          eV             eV             eV
 TOT:       8.8      -588.7320      -591.9119      -588.7573
 TOT:       9.8      -588.2469      -591.4679      -588.3085
 TOT:      10.3      -587.8200      -591.1384      -587.8709
 TOT:      10.8      -587.7559      -591.1019      -587.7985
 TOT:      11.4      -587.7044      -591.1017      -587.7405
 TOT:      11.9      -587.7184      -591.1058      -587.7516
 TOT:      12.4      -587.7194      -591.1045      -587.7524
 TOT:      13.2      -587.7166      -591.1047      -587.7497
 TOT:      14.1      -587.7172      -591.1048      -587.7503
 TOT:      15.0      -587.7173      -591.1048      -587.7504
 TOT:      15.9      -587.7171      -591.1048      -587.7502
 TOT:      16.8      -587.7173      -591.1048      -587.7504
 TOT:      17.6      -587.7173      -591.1048      -587.7504
 TOT:      stopping    0.000000
$ 
                                        ↑これ
</code>
''selfcons''と書かれている欄の一番下の数字が全エネルギー（単位はeV）である。
===== 実習2 =====
  * C-O間距離を$0.8$Åから$1.8$Åまで変化させながらエネルギーを計算して、横軸C-O間距離、縦軸結合エネルギーのグラフを描く（余裕があれば他の分子についても計算をしてみよ）

  * [[モースポテンシャル]]
  * [[CO分子の振動状態]]