別件で調べ物をしていたら10年くらい前に検討したエントリを見つけたので、これをpsi4のSAPTでやったらどういう結果になるかな？と。
今回psi4コマンドで流すのでinputはこんな感じです(いずれはpsikitでサクッと計算できるようにしたい)。EmacsやめちゃったからVimでgamessのインプットいじってつくった。繰り返しの技を覚えたけど、こういうの手になじまないとめんどくさいですね。
molecule cation_pi { 1 1 N 26.56861 22.88472 72.47206 C 25.05591 23.05763 72.38527 C 27.18276 23.04159 71.08205 C 27.14794 23.94295 73.40648 C 26.89124 21.49601 73.01604 H 24.65761 22.29324 71.71973 H 24.84300 24.05067 71.99229 H 26.95089 24.04166 70.71972 H 26.73892 22.29197 70.42917 H 28.26387 22.89551 71.15590 H 26.89847 24.92550 73.00913 H 28.22884 23.81270 73.44995 H 26.70726 23.81001 74.39332 H 26.46855 20.75492 72.33939 H 26.44825 21.40230 74.00640 H 27.97398 21.38820 73.06968 H 24.63622 22.94695 73.38408 -- 0 1 C 30.61324 20.49644 72.40614 C 30.47183 21.03406 71.19420 C 30.79531 21.59051 73.36029 C 30.75216 22.77934 72.62888 C 31.00908 21.63948 74.74382 N 30.43745 22.45766 71.27884 C 30.92671 24.02761 73.23103 C 31.18435 22.86580 75.34238 C 31.14962 24.05129 74.58925 H 30.33993 20.55500 70.23268 H 31.04481 20.72825 75.32842 H 30.86231 23.00634 70.52231 H 30.90815 24.93911 72.64846 H 31.36065 22.92555 76.40870 H 31.30378 24.99933 75.08996 H 30.61446 19.44515 72.64676 units angstrom no_reorient symmetry c1 } set basis 6-31g_d_p_ energy('sapt0')
10年も経つと6-31G**でもサクッと終わるのね。subaracy!
結果はこのあたりを見ればいいんですかね。ESも効いてるけどInduction(CT+mixみたいな項?)とDispersionも効いてるから、次はカチオンの有無でトリプトファンのそれぞれの原子のRESP chargeがどのくらい変化するかとか双極子どうなっているかとか確認すればいいんかな。
Special recipe for scaled SAPT0 (see Manual): Electrostatics sSAPT0 -13.53632077 [kcal/mol] Exchange sSAPT0 10.77916633 [kcal/mol] Induction sSAPT0 -5.12922298 [kcal/mol] Dispersion sSAPT0 -5.69289464 [kcal/mol] Total sSAPT0 -13.57927206 [kcal/mol] -------------------------------------------------------