Cation-π相互作用をSAPT0してみた。

Sun, 03 Mar 2019 16:20:38 +0919

別件で調べ物をしていたら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]     
-------------------------------------------------------

Psi4 install error

Sat, 08 Sep 2018 07:49:51 +0919

I got a this error. Conda always irritates me ;-)

>>> import psi4
Traceback (most recent call last):
  File "/Users/kzfm/anaconda3/lib/python3.6/site-packages/psi4/__init__.py", line 55, in <module>
    from . import core
ImportError: dlopen(/Users/kzfm/anaconda3/lib/python3.6/site-packages/psi4/core.so, 2): Library not loaded: @rpath/libiomp5.dylib
  Referenced from: /Users/kzfm/anaconda3/lib/python3.6/site-packages/psi4/core.so
  Reason: image not found

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/Users/kzfm/anaconda3/lib/python3.6/site-packages/psi4/__init__.py", line 60, in <module>
    raise ImportError("{0}".format(err))
ImportError: dlopen(/Users/kzfm/anaconda3/lib/python3.6/site-packages/psi4/core.so, 2): Library not loaded: @rpath/libiomp5.dylib
  Referenced from: /Users/kzfm/anaconda3/lib/python3.6/site-packages/psi4/core.so
  Reason: image not found

@rpathという変数がMacでは何らかの理由により正しく設定されないらしい

朝からイラツキすぎてまともに解決する気力もないので、当面放置するかも。これから買い物いかないといけないし。それにしてもCondaよくわからんエラー多すぎるわ。RDKitもそうだったしさー。あんまよくないんじゃないの？職場では使ってないしー

forum見たけど報告されてないっぽいし、あとで気力が戻ったら投げておきます（未定）。

追記 0800くらい

エントリ書いて朝ごはん食べたら冷製になったのでなおした。

$ install_name_tool -change @rpath/libiomp5.dylib /Users/kzfm/anaconda3/pkgs/intel-openmp-2018.0.0-h8158457_8/lib/libiomp5.dylib /Users/kzfm/anaconda3/lib/python3.6/site-packages/psi4/core.so

psi4 / Drkcore

Cation-π相互作用をSAPT0してみた。

Psi4 install error

追記 0800くらい