Coordination Bond

遷移状態計算って、やってみると意外となんとかできるもんですが、それでも基底状態と比較して、やや面倒なことも多いと思います。 私はラジカル系の計算をすることが多いのですが、ラジカル的な解を求めようとしているのに、勝手にイオン的な（閉殻的な）遷移状態がでてきてしまったり、難しいことが多いです。 なかなか、ノウハウを共有している資料も多くないので、解決に向けたノウハウを少しずつ貯めて行きたいと思います。（Gaussian09を用いたDFT法です。） （もしかしたら、 こちら も役立つかも？） [閉殻分子同士のラジカル性の反応を見たいのに、勝手に求核攻撃のような解に収束してしまう] 一重項分子の結合のホモリシスを伴いながら進行する系で、どう考えても遷移状態ではラジカル性がないとおかしいのに、求まったTSの軌道を見ると完全に閉殻、つまりイオニックな反応になってしまったことがありました。 三重項として計算をすすめて構造を見つけた後に、その構造を出発点に、もう一度一重項として計算を行う方法があるようです。とんでもないことをやっているようにも聞こえますが、αスピン、βスピンの距離が離れた状態だと、 βスピンがひっくり返って 三重項となった状態と、あまりエネルギーは変わらないようです。これを利用すると、閉殻的な、ニセの一重項にトラップされることなく、望みの電子配置に極めて近い構造、軌道、から計算を始められるようです。 （新たに求まった、開核性の一重項の電子配置にある遷移状態が、閉殻性の一重項よりも不安定だったら、、どうしたらいいんでしょうね。最適な汎関数を探す旅に出なければならないかもしれません） Gaussian先生の構造最適化は、他の計算ソフトと比較して、とても優秀だそうです。ですが、その能力を過信してはいけません。（どんなアルゴリズムで最適化を進めているか、理解してませんよね？私もですが、、気になるあなたはBernyアルゴリズム、について検索すると幸せになれるかもしれません。） 他に、より安定な解があるのに、全然違う準安定な解を拾ってくる可能性については、常に気を配らなければなりません。 結合が安定にできる原子間距離にあるときの計算は、結構精度よく再現できているそうですが、原子間距離が離れてくると、誤差が大きくなってきます。 DFT法ではなく、