あるケミストの独り言(winchemwinの日記)

ケミスト(化学者)の視点で、面白そうな情報(シミュレーション関係など)を発信

分子軌道情報(HOMO/LUMO)自動算出アプリ(CustumTkinter) その7

 これまでのシリーズでは軌道情報の自動計算アプリのコードの紹介をしてきました。今回はアプリの利用方法について簡単に紹介したいと思います。
 以下はこのシリーズの最初の回にも紹介したアプリの起動画面になります。

分子軌道情報自動算出アプリ 起動画面

 画面構成ですが、ファイルの選択、計算方法の設定、算出する軌道レベルの設定からなっています。
 一番上はSMILES記載ファイルの選択部で、右端の「select」ボタンを押すことで、ファイル選択が画面が現れますので、必要なファイルを選択してOKを押すことで読み込みファイルの設定が行なえます。ファイル形式は基本的にはテキスト形式ですが拡張子として「.txt」と「.smi」が読み込み可能です。
 例えば以下のようなSMILESリストを読み込んで連続的に計算させることが可能です。
 

 量子化学計算はMethod(計算手法:HF, DFTなど)、Function (汎関数B3LYP, PBEなど)、基底関数(3-21G, 6-31G(d)など)をプルダウンメニューから選択することができます。選択できる汎関数や基底関数を増やしたい場合はコード内のリストに付け加えることで実施できるようになるかと思います。
 Options設定として、計算機側の設定(Thred数、メモリ)と分子の状態(電荷、スピン多重度)の設定も行えるようにしています。こちらは直接数値を入力する形です。

 あと、算出する軌道レベルの数を右側の「HOMO/LUMO Level」の箇所で設定します。デフォルトでは5となっており、HOMO-5、LUMO+5まで計算する設定になっています。

 すべての設定が終えたらアプリ下部の「Quantum Calculation Start」ボタンを押すことで計算がスタートします。アプリを終了したい場合は右下の「Quit」から終了できます。

 計算結果はそれぞれ「Calcd_HOMO_List.csv」、「Calcd_LUNO_List.csv」に保存されます。
 実際の結果は以下のような形で保存されています。

分子軌道リスト計算結果 1
<>
分子軌道リスト結果 2

 それぞれのSMILESに対して、各軌道のエネルギー情報のリストが自動的に作成されてきています。また右端には振動計算(FreqNega)の結果も追加されており、「FALSE」は負の振動数がない、すなわち一応安定構造で計算は実施されているであろうという目安の結果も示しています。

 機械学習等で多くの量子化学計算データを必要とする際にある程度便利に使っていただけるのではないかと思いますので、参考に紹介させていただきました。
 
 分子軌道情報の算出アプリについては今回で終了としたいと思います。
 次回からは屈折率(リスト)の自動算出アプリについて紹介してゆければと思っています。