0分子間の相互作用は、計算化学が扱う主要な問題の一つです。単純に考えると、例えばAという分子とBという分子が相互作用してABという錯体を形成する場合、A,Bが単独で存在したときのエネルギーと錯体ABのエネルギーを比較して、その相互作用がどれぐらい強いかを議論するということになります。ここには思いがけない落とし穴が存在します。Error:基底関数重なり誤差)です。以下、GAMESSでのBSSEの計算方法とモデル化学による差について見てみたいと思います。A-B]における計算をもう少し詳しく見ると、Aを単独で計算するのと、A-Bという錯体で計算するのに違いがあることがわかります。錯体A-Bの中のAは、Bの基底関数を取り込むことができます。普通はこの「取り込み」はイコール相互作用ということになりますが、小規模の基底関数(例えばSTO-3Gや3-21Gなど)では、Bの基底関数をAの基底関数の分極関数やDiffuse関数のように使われてしまい、結果、相互作用のエネルギーだけでなく基底関数の改良による安定化エネルギーが取り込まれた形で数値が出ます。この「基底関数の改良による安定化エネルギー」こそがBSSEです。どんなに大規模な基底関数でもこの寄与が(大幅に縮小するとはいえ)完全に取り除けないことは、容易に想像できると思います。BSSEを求めることは、相互作用を扱う場合必須であると言えます。(※これから紹介するBSSEの計算(counterpoise補償法)は、BSSEと完全にイコールではありません。よって、”BSSEを補正した”と言っても完全に除けたわけではなく、その比率が目的とする相互作用のエネルギーに比べて十分小さくするぐらいの意味合いで考えるべきでしょう。時には補正しすぎることもあるようです。)BSSEの推定メチルイソシアニド(CH3-NC)と三フッ化ホウ素(BF3)の錯体を例に、実際にBSSEを計算してみましょう。本来、2分子間の相互作用はGibbs自由エネルギーで取り扱うべきですが、BSSEは全電子エネルギーに影響を与えるということと、計算時間の都合で振動計算は省略します(手抜きじゃないですよ?)。原系も生成系も対称性が高い構造をしていますので、メチルイソシアニドと錯体はC3vで、三フッ化ホウ素はD3hで計算を行うことにします。それぞれの入力構造($DATAグループ)は以下のとおりです。非常に単純な分子なので、モデリングソフトを使うよりも手計算で組んだほうが早いでしょう。$ENDこれらについて、理論モデルとしてRHFとRPBE0(hybrid-DFT)を、基底関数として3-21Gと6-31G(d)を使い(つまり4種類の異なるモデル化学で)構造最適化を行ってみました。OPTTOL=1.0E-05に設定しています($STATPT内)。それぞれ見積もられた相互作用エネルギーΔEは以下の通りです。-286652.6213ΔE-14.27117145モデル化学によって、相互作用エネルギー(ΔE)の値は大きく異なります。3-21Gを使ったときに絶対値が大きくなる傾向にありますが、これはBSSEが多く含まれているためです。問題のBSSEですが、上記計算だけでは求められません。BoysのCounterpoise補償法では、あと4つエネルギー計算をする必要があります。Counterpoise補償法では、「ゴースト原子」と呼ばれる「原子核がない、空の軌道だけがある原子」を利用してBSSEを見積もります(下図)。最適化された錯体から一方のみをそっくり取り出し、単独でエネルギー計算(RUNTYP=ENERGY)を行います(E1)。次に、錯体の中で相手の方をゴースト化し、エネルギー計算を行い(E2)、両者の差を取ります(E1-E2)。構造は全く同じで、相互作用する相手分子の軌道がある時とない時のエネルギー差が「BSSE」であるとするわけです。complexの最適化構造を抜き出します。ATOMS」の方を使います。-2.0710495435これをCH3-NCとBF3とに分けます。これら2つのエネルギー計算の和がE1の計算に相当します。-2.0710495435次に、CH3-NCとBF3とをそれぞれゴースト化します。陽子数(2列目)を負の値にすればOKです。これら2つのエネルギー計算の和がE2の計算に相当します。-2.0710495435あとは、E1-E2でBSSE(counterpoise)が算出されます。前出の4つのモデル化学についてBSSEと補正ΔEを示します。5.530825435corr.ΔE-8.7403460153-21G基底を使った場合、BSSEが非常に大きくなっていることが一目瞭然です。6-31G(d)基底ではBSSEが大幅に減少しています。RPBE0ではBSSE補正がきちっと効いているのか、補正後のエネルギーが3-21Gと6-31G(d)でほとんど変わりません。Exchangeから入手)でも同様に計算した結果を示します。2.460337491corr.ΔE-12.40591446RPBE0/6-311+G(2d)の結果は流石に優れています。kcal/mol以内。DZVPの結果も次いで良く、しかも6-311+G(2d)よりも3倍以上速く計算することができます。ゴースト化錯体のエネルギー計算で比較)のですが、BSSEはDZVPが半分、最終的な補正エネルギーも6-311+G(2d)に近いということで、本系に対して非常に有効な基底関数のようです。諸熊分割によるRHF法でのBSSE計算RHF法に限り、上記のような複数の計算を行わなくてもCounterpoise補償法によるBSSEを計算することができます。IATM(1)=n」と記述します。SCFはできないので$SCFグループでDIRSCF=.F.にしています(デフォルトが.F.ですので、何も書かなくても構いません)。$ENDIATM(1)=nの「n」は、モノマーの区切りを示すもので、上記のケースでは6番目までがメチルイソシアニドですので「6」としています。3成分以上の錯体だった場合は、各モノマーの最後の原子をカンマ区切りの番号で指定します(最後の1つは必要ありません)。出力は一番最後で、以下のような内容です(抜粋)。CORRECTIONが目的の値です。先の手計算とも合致しています。コメントをキャンセルコメントを投稿するにはログインしてください。