化学式ワープロと分子軌道計算～量子力学もコンピュータで解こう～ / 情報処理概論 / 化学・バイオ工学科 / 工学部 / 山形大学

原子と原子は静電気力で結びつく。化学結合にはイオン結合、共有結合、金属結合などがある。固体には、塩化ナトリウムなどのイオン結晶、ダイヤモンドやケイ素などの共有結合からなる共有結晶、銀などの金属結合からなる金属結晶、供給結合からなる分子が弱い力でくっつくドライアイスなどの分子結晶などがある。液体は水のように分子からなる液体、塩化ナトリウムを高温でナトリウムイオンと塩化物イオンがイオン液体（熔融塩）、銀が自由に動き回る溶融金属、などがある。ドライアイスやヨウ素などの分子結晶は昇華しやすく、液体にはなりにくい。

有機化学美術館

分子全体に広がる一電子空間軌道関数である分子軌道によって、分子を構成する個々の電子の状態が記述される。この分子軌道を計算して、分子の電子状態を求める方法が分子軌道法である。コンピュータで計算する手法である．なお，分子軌道法については，物理化学IIIの授業で詳しく触れます．

Fig クーロンポテンシャルとＨＯＭＯとＬＵＭＯ
©Copyright Takaaki Shiraya all rights reserved.

上図を３Ｄで可視化したものが下図です。クーロンポテンシャルは原子核のまわりにできるへこみみたいなものです。電子はそこにたまる水みたいなものです。

一分子では，世の中にあるものは作れません．最低 1molは必要！

SCIGRESS（サイグレス）

Gaussian(ガウシャン)

Winmostar(ウィンモスター)

Gaussian < パッケージプログラム | 計算科学研究センター

IT創薬

オンライン反応経路および化学構造の探索

クラウド上のＡＩコンピューティングサービス。これから台頭してくると思われる。

Reaxys (リアクシス)

SciFinder，化学構造による特許の検索可

化学工学系ソフトウェア

化学工学系ソフトウェア | モルシス

Aspen Plus

分子動力学

分子動力学シミュレーション | 京

リチウムイオン電池　～分子の宇宙から未来の電池へ～

SCIGRESS MD-ME

生命科学系分子科学計算ソフトウェア

結晶構造可視化ソフト

VESTA

CrystalMaker (クリスタルメーカー)

学術情報基盤センターが提供するソフト

BIOVIA Draw

分子の構造，分子式を書くときに使います．

SCIGRESS

研究室のPCでの利用

山形大学工学部学術情報基盤センタ－が提供するソフトウェアリスト

学術情報基盤センターでの利用

実習室のWindows PCにログオン後
[スタート]　→ [SCIGRESS 2.8.1] → [Workspace]

計算例）フェリシアニド([Fe(CN)₆]^3-)錯体の分子軌道計算

同時使用ライセンス（フローティングライセンス）

80プログラムの同時起動が可能

分子軌道計算では初期配座が重要

ChemSpider:http://www.chemspider.com/
alpha-Cyclodextrinで検索する。3D表示する。mol形式でダウンロードする。サイグレスで開く。

構造最適化を行う。MOに関してはライセンスの上限があるので先着順です！構造最適化が終わったらＨＯＭＯとＬＵＭＯの波動関数を計算する。

エネルギー化学への応用

量子力学計算をエネルギー化学に応用できます。波動関数のエネルギー準位と酸化還元電位の関係を調べることができます。

酸化還元電位の算出

Fe3+/Fe2+の半反応式は次式で表されている．この半反応式の酸化還元電位は標準水素電極(NHE)に対し+0.771 Vである．

{Fe}^{2+} = {Fe}^{3+} + e^{-}

（反応式- 1）

HOMO軌道とLUMO軌道のエネルギー順位をSCIGRESS(PM6 geometry with PM6 wavefunction※1)を用いて計算すると下図のようになった． Fe³⁺のLUMO軌道のエネルギー準位は，-30.23eVとなり，Fe²⁺のエネルギー準位は-28.68eVとなった．さらに，水素分子(H₂)のHOMO軌道のエネルギー準位は，-15.04eVとなった．水素イオン(H⁺)のLUMO軌道のエネルギー準位をSCIGRESSで計算すると，「No surface was found for molecular orbital 1 since it does not contain a value as large as 0.060.」とエラーが表示されLUMO軌道のエネルギー準位を計算できない．

Waiting for resources to become available...
Keywords used in this job are as follows:
 PM6 NOMM NOXYZ XYZ NODIIS GRAPH T=10D CHARGE=1  

MO-GManager 1.1.1.0 Apr  2 2018 21:50:58

MO-GComputeEngine 2.9.0 Apr  2 2018 21:50:15

Final heat of formation =     311.4630 kcal/mol
Calculation done.
TabulatorManager 2.9.0 Apr  2 2018 21:45:33

TabulatorComputeEngine 2.9.0 Apr  2 2018 21:44:57

3	3	3	(529169.25,529169.25,529169.25) Angstrom	(-529169.25,-529169.25,-529169.25) Angstrom	Molecular orbital values:  0.060, -0.060	Using molecular orbitals 1 - 1
No surface was found for molecular orbital 1 since it
 does not contain a value as large as  0.060.
Computation time: 00:00:00.388

ギブスの自由エネルギーの計算はどうやるの？

Gaussianでギブスの自由エネルギーを求める方法は下記のキーワード検索すると見つかることがあります．検索キーワード：gibbs free energy gaussian Google
文献などのExampleを探しましょう．基本は，物理化学