2016 年 15 巻 3 号 p. 60-62
We have developed a method for quantifying the intermolecular forces in supramolecular systems by modifying our original method of coarse-graining intermolecular vibrations. We evaluated the true and apparent intermolecular stiffness constants for 21 dimers composed of formic acid-, acetic acid-, trichloroacetic acid-, formamide-, formamidine- or urea-monomer. In this method, the atomic displacement vectors of a dimer were projected onto a subspace spanned by bases corresponding to 12 relative translational and rotational motions and several intramolecular vibrations that are coupled with intermolecular vibrations. The intermolecular stiffness constants showed moderate linearity with the corresponding dimerization energies. The apparent stiffness constant can be explained by a mechanical model using inter- and intramolecular stiffness constants of the constituent monomers.
水素結合やハロゲン結合といった弱い分子間相互作用は分子間振動の起源でもあり,その定量化はTHz分光学の発展に伴い注目されてきた [1,2].最近我々は,分子軌道法に基づく基準振動解析から得られる振動数と変位ベクトルを使って分子間結合の剛性定数を定量化する方法を開発した [3].この方法では,N個の原子からなるホモ二量体の基準振動を表す3N次元変位ベクトルを12種の基本変位ベクトルが張る部分空間(粗視化空間)に射影することで振動の自由度を12次元に縮約する.しかし柔軟な分子では分子間振動に分子内振動が強く結合するため [4],12次元粗視化空間上では振動を精度よく表せない.この場合,各構成分子の最低波数の振動モードを基底に追加することにより粗視化空間の質が劇的に向上する [5].
粗視化空間の次元の増加に伴い剛性定数は増加する挙動を示す.3N次元で算出した値は"真の”剛性定数
ヘテロ二量体の場合,分子間振動に対する分子内振動の結合の強さは構成分子の種類によって異なる.各構成分子の最低波数の分子内振動を基底に追加しても適切な粗視化空間が得られるとは限らないため,基底系を構築する際の基準が必要である.
本論文では,上述の解析方法をヘテロ二量体にも適用できるように理論を一般化したことについて述べ,これを数種の水素結合性ホモ/ヘテロ二量体の解析に用いた結果を報告する.
粗視化空間の座標系における慣性行列
二量体の基準振動解析によって3N次元変位ベクトルが得られる.これを12個の基本変位ベクトルで張られる粗視化空間に射影したときの加重変位行列は近似的に
分子軌道計算にはGaussian03またはGaussian09を使い,各構成分子の単量体および二量体の構造最適化,基準振動解析をHF/6-311G (d,p)レベルで行った.基準振動の粗視化解析と分子内/分子間剛性の算出には,当研究室で開発したプログラム [5]を改良して用いた.
ギ酸,酢酸,トリクロロ酢酸,ホルムアミド,ホルムアミジン,尿素の単量体,これらの分子からなる6組のホモ二量体,および15組のヘテロ二量体について基準振動解析を行った.剛性行列の算出は,水素結合の方向がx軸となるような座標系(Figure 1)のもとで行った.分子内振動を基底に追加する際には,構成分子が単量体の時と二量体を形成した時とで配向を一致させる必要がある.そこで構成分子の慣性主軸を一致させるように回転することで,座標系を一意的に決定できるようにした.
Coordinate system selected for formic acid + urea dimer. Thick arrows indicate the principal axes of inertia for each molecule.
この座標系のもとで,x軸上の反対称並進を表す基本変位に対応する剛性定数(
The stiffness constant (
擬二原子近似のもとで算出した分子間剛性
Correlation between
分子間振動の粗視化にあたっては分子内の柔軟性を考慮したみかけの剛性定数を簡便に算出することが重要となる.連成振動子をモデルとした数値計算により,ホモ/ヘテロ二量体の
Comparison between