遷移状態のCBS-QB3計算ができません
光化学反応を研究している大学院生です．Gaussian16で遷移状態の高精度エネルギー計算を行おうと思っているのですが，うまくいきません．
DFT b3LYP 6-31++G(d,p)を用い，QST2法で反応物と生成物から遷移状態の構造を計算しました．その構造を用いて，CBS-QB3法でエネルギー計算を行いました．すると，遷移状態の構造が構造最適化されてしまい，生成物の構造が出力されました．また，エネルギーも生成物のものが出力されてしまいました．もちろんキーワードとしてoptを入れているわけではありません．
構造最適化を済ませたchkファイルを読み込ませれば良いとの情報を得て試したのですが，同様の結果になりました．
遷移状態構造のCBS-QB3法によるエネルギー一点計算はどのように行うのでしょうか？ある論文では，遷移状態のエネルギーもCBS-QB3法によって求めており，原理的には不可能ではないと思っています．
また，そもそもエネルギー計算しかできない計算方法なのに，なぜインプットした構造が変わるのか疑問に思いました．
お忙しいところ恐れ入りますが，ご存知でしたらご教示お願い致します．
最後になりましたが，このようなサイトを作っていただき，大変感謝しております．いつも参考にさせていただいています．
そのキーワードだけだと chk ファイルの読み込みが出来ていないと思います。
freq=readfc とか guess=read とか geom=check とか chk ファイルを読み込むキーワードを色々調べてみて下さい。
あと、キーワードセクションと電荷・多重度の間にタイトルセクションを入れて下さい。
参考として貼って頂いた blog をみますと、”If you want to obtain a G3 energy of a transition state, or a minimum calculated using any basis set other than 6-31G(d), you will need to start the G3 job from a checkpoint file from a frequency calculation containing a geometry and force constant matrix (Hessian) computed previously at your preferred level of theory.” と書いてあります。
一応確認ですが、今回使用している chk ファイルは、qst2 計算の後に freq 計算した時の chk ファイルですか？
また、geom=connectivity は不要だと思います。（ONIOM などで MM を使う場合は、geom=connectivity を使った方が良いそうですが。。。）chk ファイルから座標を読み込むときは、geom=check を使い、座標情報を input に書く必要はありません。参考に貼って頂いた blog では geom=allcheck を使用しています。これだと、座標以外にも電荷、多重度、タイトルセクションなども読み込まれるため、input ファイルをさらに省略することができます。
ご参考になれば幸いです。
