遷移状態のCBS-QB3計算ができません
光化学反応を研究している大学院生です.Gaussian16で遷移状態の高精度エネルギー計算を行おうと思っているのですが,うまくいきません.
DFT b3LYP 6-31++G(d,p)を用い,QST2法で反応物と生成物から遷移状態の構造を計算しました.その構造を用いて,CBS-QB3法でエネルギー計算を行いました.すると,遷移状態の構造が構造最適化されてしまい,生成物の構造が出力されました.また,エネルギーも生成物のものが出力されてしまいました.もちろんキーワードとしてoptを入れているわけではありません.
構造最適化を済ませたchkファイルを読み込ませれば良いとの情報を得て試したのですが,同様の結果になりました.
遷移状態構造のCBS-QB3法によるエネルギー一点計算はどのように行うのでしょうか?ある論文では,遷移状態のエネルギーもCBS-QB3法によって求めており,原理的には不可能ではないと思っています.
また,そもそもエネルギー計算しかできない計算方法なのに,なぜインプットした構造が変わるのか疑問に思いました.
お忙しいところ恐れ入りますが,ご存知でしたらご教示お願い致します.
最後になりましたが,このようなサイトを作っていただき,大変感謝しております.いつも参考にさせていただいています.
管理人様
ご返信ありがとうございました.
前回貼り付けたinputファイルは,遷移状態の構造を入力して計算したものを貼り付けました.また,タイトルセクションはこの掲示板に貼る際に,省略してしまいました.紛らわしくなってしまい申し訳ありません.
調べる中で,CBS-QB3計算は,構造最適化計算が一連の計算の中に入っていることが分かりました.
ご指摘の通り,chkファイルを読み込めるようにし,様々なキーワードを試してみましたが,うまくいったものはありませんでした.
以下のように数撃ちゃ当たる精神で(間違っているものもあると思います…)色々試してみましたが,すべて構造最適化されるもしくは計算が走りませんでした.
chkファイル読み込み
・# cbs-qb3 geom=allcheck
・# cbs-qb3 noopt geom=allcheck
・# cbs-qb3 SP geom=allcheck
・# cbs-qb3 opt=maxcyc=0 geom=allcheck
・# freq b3lyp/6-31++g(d,p)//cbs-qb3 geom=allcheck
・# cbs-qb3 guess=read geom=allcheck
TS構造をGauss Viewで読み込み,新たにchkファイルを作る
・# cbs-qb3 guess=read geom=connectivity
・# cbs-qb3 guess=only geom=connectivity
・# cbs-qb3 opt=TS geom=connectivity
・# freq b3lyp/6-31++g(d,p)//cbs-qb3 geom=connectivity
QST2計算にCBS-QB3も入れてしまう
・# opt=qst2 freq b3lyp/6-31++g(d,p)//cbs-qb3 geom=connectivity
まだ何か検討すべき余地,もしくは間違いなどがありましたら,ご指摘いただけると嬉しいです.
Gaussianのサポートに直接問い合わせてみようと思います.回答がいただけたら,こちらに投稿いたします.
お時間を割いてくださり,大変ありがとうございました.