HOME 教育状況公表 本日

化学式ワープロと分子軌道計算~量子力学もコンピュータで解こう~

情報処理概論 WebClass files syllabus 52210 C1 情報機器の準備 ウェブ講義の接続

工場では実際にモノを作る前に,CAEによる論理シミュレーションが行われる(工場のしくみ,p.166). 分子軌道計算は,材料の反応予測や材料物性の予測に利用される. プラントを建設する前に,論理シミュレータでプラントを設計する.

材料物性の予測・・・分子軌道計算。エネルギー最少。 ハードウェア込みで500万円超?、アカデミックディスカウント。 計算科学の歴史。1953年 MD、1973年、ガウシャンPDP-11、1985年80386DX。

特許と化学式

電子出願ソフトサポートサイト

特許 の明細書に分子式を記述するときの 画像ファイル の形式は?

検索

化学構造の表現・表記

SMILES記法

International Chemical Identifier(InChI)

SMILES記法チュートリアル

化学式ワープロ

ChemBioDraw/ChemDraw

BIOVIA Draw(バイオビアドロー)

ChemDraw

ChemDraw for iPad

コンピュータの進歩と計算化学

1970年 Gussian登場, PDP-11,非経験的分子軌道法,ab initio分子軌道法

1983年 MOPACが登場, 半経験的分子軌道法

分子軌道法

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

有機化学美術館

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

Fig クーロンポテンシャルとHOMOとLUMO
©Copyright Takaaki Shiraya all rights reserved.

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

Fig クーロンポテンシャルとHOMOのイメージ( マセマティカによる可視化)
PVDF (©赤間未行

一分子では,世の中にあるものは作れません.最低 1molは必要!

SCIGRESS(サイグレス)

Gaussian(ガウシャン)

Winmostar(ウィンモスター)

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

IT創薬

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

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

Reaxys (リアクシス)

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

化学工学系ソフトウェア

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

Aspen Plus

分子動力学

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

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

SCIGRESS MD-ME

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

結晶構造 可視化ソフト

VESTA

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

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

BIOVIA Draw

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

SCIGRESS

研究室のPCでの利用

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

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

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

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

同時使用ライセンス(フローティングライセンス)

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

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

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

シクロデキストリンの構造

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

エネルギー化学への応用

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

酸化還元電位の算出

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

Fe 2+ = Fe 3+ + e -
反応式- 1

HOMO軌道とLUMO軌道のエネルギー順位をSCIGRESS(PM6 geometry with PM6 wavefunction※1)を用いて計算すると下図のようになった. Fe3+のLUMO軌道のエネルギー準位は,-30.23eVとなり,Fe2+のエネルギー準位は-28.68eVとなった. さらに,水素分子(H2)の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
 

Fe3+のLUMO軌道とFe2+のHOMO軌道のエネルギー計算結果

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

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

新型コロナウイルスの影響による在宅課題

分子軌道計算ソフトを使ったα-シクロデキストリンのHOMO-LUMOエネルギーの計算

Winmostarのインストール
Winmostarのインストールをしてください.学生の方は学生版が利用できます.あとは,ウェブページの指示に従ってインストールしてください. 本学の電子メールの使い方がわからない方は,「学生向け電子メールの使い方」をご覧ください.
α-シクロデキストリンの3次元のmolファイルをダウンロード
Chemspiderを使ってα-シクロデキストリン(alpha-cyclodextrin)の3次元のmolファイルをダウンロードしてください.
Winmostarの起動
Winmostarを起動して,molファイルを開いてください.
構造最適化
構造最適化を行います.メニューバーの「半経験QM] -> [MOPAC] -> 「MOP7W70の実行」をクリックしてください.
MOPAC Setupのダイアログウィンドウが起動したら,OKをクリックしてください. ファイル名の保存ダイアログウィンドが立ち上がった[保存]をクリックしてください.なお,ダイアログウィンドウは,起動しないこともあります. メモ帳が軌道したら,最適化完了です.
エネルギー計算
エネルギー計算を行います.メニューバーの「半経験QM] -> [CNDO/S] -> 「実行」をクリックしてください.
CNDO/S Setupのダイアログウィンドウが起動したら,OKをクリックしてください. ファイル名の保存ダイアログウィンドが立ち上がった[保存]をクリックしてください.なお,ダイアログウィンドウは,起動しないこともあります. CNDO/S UV-Vis-Spectrumなど3つのダイアログウィンドウが起動したエネルギー計算完了です.
HOMO LUMOエネルギーの表示
Energy Level… ダイアログウィンドウを見てください.HOMO:の後の数字と同じUnitの値がHOMO軌道のエネルギーになります.その値を読みとりウェブクラスに提出してください. また,HOMO:のUnitの値に1を加えた値がLUMO軌道となります.LUMO軌道のエネルギーの値をウェブクラスに提出してください.

QRコード
https://edu.yz.yamagata-u.ac.jp/Public/52210/52210_06.asp

SSLの仕組み

このマークはこのページで 著作権が明示されない部分について付けられたものです。

山形大学 データベースアメニティ研究所
〒992-8510 山形県米沢市城南4丁目3-16
3号館(物質化学工学科棟) 3-3301
仁科・立花・伊藤研究室 准教授 伊藤智博
0238-26-3573
http://amenity.yz.yamagata-u.ac.jp/

Copyright ©1996- 2020 Databese Amenity Laboratory of Virtual Research Institute,  Yamagata University All Rights Reserved.