NTO二聚体分子间相互作用的理论研究

在DFT-B3LYP/6-311＋＋G**水平上求得NTO二聚体势能面上六种优化构型和电子结构. 经基组叠加误差(BSSE)和零点能(ZPE)校正, 求得分子间最大相互作用能为－53.66 kJ/mol. 二子体系间的电荷转移很少. 由自然键轨道分析揭示了相互作用的本质. 对优化构型进行振动分析, 并基于统计热力学求得200.0～800.0 K温度范围从单体形成二聚体的热力学性质变化. 发现二聚主要由强氢键所贡献, 但结合能大小并不为氢键所完全决定. 二聚过程在较低温度或常温下能自发进行.     

2-氨基苯并噻唑的结构及激发态质子转移动力学

通过傅里叶变换红外光谱（FTIR）、傅里叶变换拉曼（FT-Raman）和488 nm拉曼光谱,结合密度泛函理论（DFT）计算研究了2-氨基苯并噻唑（ABT）在晶态和溶剂中的二聚体结构,并解释了质子性溶剂中ABT二聚体与溶剂分子间的氢键作用.电子光谱实验揭示了ABT二聚体的光物理和光化学反应;紫外吸收和荧光发射光谱结果表明,溶剂、激发波长和pH值对ABT二聚件激发态衰变具有调控作用;含时密度泛函理论（TD-DFT）解释了ABT二聚体双荧光现象,提出了高激发态的质子转移机理.     

(TCNQ)22-基氢键分子磁体表现出非常强的反铁磁偶合作用

合成并表征了一种新的离子对化合物[4-NH2-Py][TCNQ](其中4-NH2-Py+是4-氨基吡啶阳离子,TCNQ-为7,7,8,8-四氰基对苯二醌二甲烷自由基阴离子)。在该离子对化合物晶体中,2个TCNQ-离子形成了面对面堆积的二聚体;阴离子中的氰基分别和阳离子上的氨基、吡啶质子化氮原子之间存在非常强的分子间氢键。通过氢键作用,相邻的TCNQ-二聚体被阳离子连成三维氢键网络。变温磁化率测量表明,在2～400 K温度范围内,该离子对化合物表现为抗磁性。在密度泛函理论框架下,用对称性破损方法计算了化合物晶体中π二聚体内以及通过氢键连接的相邻的TCNQ-离子之间的磁交换常数,发现π二聚体内存在非常强的反铁磁交换作用,与之相比,通过氢键连接的TCNQ-离子之间的磁交换作用可以忽略。π二聚体内强反铁磁交换作用(J/kB≈1805 K)导致了该化合物基本表现为抗磁性。     

1H-ANTA二聚体分子间相互作用的密度泛函理论研究

在DFT-B3LYP/6-311＋＋G**水平上求得1H-3-硝基-5-氨基-1,2,4-三唑(1H-ANTA)二聚体势能面上5种优化构型和电子结构. 经基组叠加误差(BSSE)和零点能(ZPE)校正, 求得分子间最大结合能为70.63 kJ/mol. 二聚体的形成使电荷向三唑环转移. 由氢键强弱推断二聚体稳定性的顺序与结合能顺序相一致, 氢键是二聚体的主要作用形式. 对优化构型进行振动分析, 并基于统计热力学求得200.0～800.0 K温度范围内单体形成二聚体的热力学性质变化. 发现在该温度范围所有二聚过程均能自发进行.     

反相气相色谱法研究1,3,5-三氨基-2,4,6-三硝基苯的表面特性

采用反相气相色谱技术(IGC)研究了4种不同粒度的1,3,5-三氨基-2,4,6-三硝基苯(TATB)的表面性质。4种不同粒度的TATB表面自由能的色散分量（γds）随着温度的升高而增加;粒度越大的粒子,其色散自由能上升越快;在较高温度下,粗颗粒TATB显示了更强的色散作用(γds=193.2 mJ/m2,353 K),粒度最小的亚微米TATB显示了最弱的色散作用(γds=64.0 mJ/m2,353 K)。由于制备方法不同和粒子大小的差异,4种TATB的表面酸碱性质显示了明显的差别,细颗粒TATB表面有较强的亲电子特性;而其他3种TATB在极性探针分子的作用下的吸附均表现为吸热吸附,表现出在分子内和分子间具有强烈的相互作用,其Ka和Kb值均为负。     

阳离子-π体系相互作用的理论研究——Ⅰ.铵离子-苯复合物构型及相互作用的密度泛函研究

运用密度泛函B3LYP/6 31G 方法对铵离子 苯复合物的 6种可能构型进行了计算研究 ,发现铵离子的 2个氢指向苯环且不具对称性的结构为能量优势结构 .根据计算得到的复合物结构特点 ,分子间作用方式及热力学参数 ,得出铵离子和苯之间的作用是氢键相互作用的结论 .     

TATB基PBX结合能的分子动力学模拟

用分子动力学（MD）方法, 模拟计算了四种氟聚合物(聚偏二氟乙烯（PVDF）、聚三氟氯乙烯（PCTFE）、氟橡胶（F₂₃₁₁)、氟树脂(F₂₃₁₄))与TATB(1,3,5- 三氨基- 2,4,6- 三硝基苯)晶体的相互作用. 结果发现, 四种氟聚物与TATB的结合能大小排序为PVDF>F₂₃₁₁>F₂₃₁₄>PCTFE, 各氟聚物在TATB不同晶面上的结合能大小排序为(001)>(010)>(100), 结合能主要由分子间氢键决定.     

3,6-二氨基-1,2,4,5-四嗪二聚体分子间相互作用的理论研究

在DFT-B3LYP/6-31G(d)水平下,求得3,6-二氨基-1,2,4,5-四嗪二聚体势能面上3种优化几何构型和电子结构。经基组叠加误差(BSSE)和零点能(ZPE)校正,求得分子间最大相互作用能为-38.88kJ/mol。电荷分布与转移分析表明,二子体系间的电荷转移很少,但接触点上氮原子和氢原子电荷变化比较大。由自然键轨道(NBO)分析揭示了分子间相互作用的本质。对优化构型进行振动分析,并基于统计热力学求得200.0~800.0 K温度范围从单体形成二聚体的热力学性质变化,发现二聚主要由强氢键所贡献,二聚过程在较低温度或常温下能自发进行。     

6-烷基鸟嘌呤与DNA碱基间氢键作用的理论研究

用密度泛函B3LYP方法在6-311＋G**基组水平上对鸟嘌呤及顺(cis-)、反式(anti-)-6-烷基鸟嘌呤(O⁶-AlkylG)与DNA碱基(胸腺嘧啶T、胞嘧啶C、腺嘌呤A、鸟嘌呤G)的氢键二聚体结构进行了优化. 在MP2/cc-pVXZ(X＝D,T)// B3LYP/6-311＋G**水平上, 采用完全基组外推方法校正了氢键二聚体的相互作用能, 并用完全均衡校正法(CP)校正了基组重叠误差(BSSE). 在B3LYP/6-311＋G**水平上计算了各氢键碱基对的全电子波函数, 并用分子中的原子理论(AIM)分析了碱基间的弱相互作用. 计算结果显示, 鸟嘌呤6-O烷基化改变了碱基间的氢键作用模式, 使碱基对发生了明显的螺旋桨式扭转和不同程度的位移, 碱基间的电子密度分布和氢键作用能明显减小. O⁶-AlkylG对DNA碱基间的氢键作用是去稳定化的, 去稳定化影响的顺序为GC＞GG＞GA≈GT. 计算结果与文献给出的实验结论基本一致.     

6-甲基-4-羟基嘧啶单体及二聚体质子转移过程的理论研究

应用密度泛函理论的B3LYP/6-311+G(d)方法研究了6-甲基-4-羟基嘧啶单体及二聚体质子转移的异构化反应.对反应势能面的研究发现,该化含物可能存在9种单体异构体,对其最稳定的单体构型进行分析.各单体间异构化反应的过渡态共有9种,反应的活化能最小为22.06 kJ/mol,最大为356.55 kJ/mol,最可能的反应路径在室温下即可进行. 研究了2种二聚体及其异构化反应的过渡态,发现二聚体均比其对应的单体稳定,而且质子转移所需要的活化能仅为20.13 kJ/mol,比单体低很多. 氢键在这种变化中起了主要作用,由单体和二聚体的总能量计算了氢键的键能.     

Towards a force field based on density fitting

Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n = 16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed.     

Model chemistry calculations of thiophene dimer interactions: origin of pi-stacking

The intermolecular interaction energies of thiophene dimers have been calculated by using an aromatic intermolecular interaction (AIMI) model (a model chemistry for the evaluation of intermolecular interactions between aromatic molecules). The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium-size basis set. The calculated interaction energies of the parallel and perpendicular thiophene dimers are -1.71 and -3.12 kcal/mol, respectively. The substantial attractive interaction in the thiophene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge transfer but rather long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases the attraction significantly. The dispersion interaction is found to be the major source of attraction in the thiophene dimer. The calculated total interaction energy of the thiophene dimer is highly orientation dependent. Although electrostatic interaction is substantially weaker than dispersion interaction, it is highly orientation dependent, and therefore electrostatic interaction play an important role in the orientation dependence of the total interaction energy. The large attractive interaction in the perpendicular dimer is the cause of the preference for the herringbone structure in the crystals of nonsubstituted oligothiophenes (alpha-terthienyls), and the steric repulsion between the beta-substituents is the cause of the pi-stacked structure in the crystals of some beta-substituted oligothiophenes.     

Theoretical study on mesogenic core structures of nematic liquid crystalline compounds

The structures and intermolecular interaction energies of 10 dimers, included in the mesogenic core structures of typical liquid crystalline (LC) compounds, are obtained at the MP2/6-31G(d) level of theory. It is proved that the dispersion energy significantly contributes to the total interaction energy of these dimers. Even when bulky substituents are introduced into the core part, the interaction energy is still large. It is also revealed that when a long intermolecular distance is provided by a high steric repulsion originating from the linkage of two phenyl groups, the dispersion energy is significantly small. However, in this range of intermolecular distances, the electrostatic energy caused by a strong quadrupole-quadrupole attractive interaction plays a dominant role, and as a result, a rather stable dimer is formed. In all 10 dimers, the dispersion, electrostatic, and exchange-repulsion energies strongly depend on the geometrical orientation of the molecules. The calculated interaction energies of these dimers are also compared with the corresponding experimentally measured viscosities. The results suggest an explicit linear relationship between the interaction energies and viscosities.     

Critical examination of the supermolecule density functional theory calculations of intermolecular interactions

The results of calculations employing twelve different combinations of exchange and correlation functionals are compared with results of ab initio calculations for two different configurations of the water dimer and three different configurations of the thymine-adenine complex. None of the density functional theory (DFT) treatments could properly reproduce the results of coupled-cluster calculations for all configurations examined. The DFT approaches perform well when the interaction energy is dominated by the electrostatic component and the dispersion energy is less important. Two mechanisms that compensate for the missing dispersion component were identified. The first one is the decrease of the magnitude of the intermolecular exchange-repulsion and the second one is the increase of the magnitude of the attractive deformation energy. For some functionals both effects are observed together, but for some other ones only the second effect occurs. The three correlation functionals that were examined were found to make only very small contributions to the deformation energy. The examination of angular and distance dependence of the interactions shows that the currently available DFT approaches are not suitable for developing intermolecular potential energy surfaces. They could however be used to find global minima on potential energy surfaces governed by intermolecular electrostatic interactions. Additional single point ab initio calculations are recommended as the means of validating optimized structures.     

Differences in structure, energy, and spectrum between neutral, protonated, and deprotonated phenol dimers: comparison of various density functionals with ab initio theory

We have carried out extensive calculations for neutral, cationic protonated, anionic deprotonated phenol dimers. The structures and energetics of this system are determined by the delicate competition between H-bonding, H-π interaction and π-π interaction. Thus, the structures, binding energies and frequencies of the dimers are studied by using a variety of functionals of density functional theory (DFT) and M?ller-Plesset second order perturbation theory (MP2) with medium and extended basis sets. The binding energies are compared with those of highly reliable coupled cluster theory with single, double, and perturbative triple excitations (CCSD(T)) at the complete basis set (CBS) limit. The neutral phenol dimer is unique in the sense that its experimental rotational constants have been measured. The geometry of the neutral phenol dimer is governed by the hydrogen bond formed by two hydroxyl groups and the H-π interaction between two aromatic rings, while the structure of the protonated/deprotonated phenol dimers is additionally governed by the electrostatic and induction effects due to the short strong hydrogen bond (SSHB) and the charges populated in the aromatic rings in the ionic systems. Our salient finding is the substantial differences in structure between neutral, protonated, and deprotonated phenol dimers. This is because the neutral dimer involves in both H(π)···O and H(π)···π interactions, the protonated dimer involves in H(π)···π interactions, and the deprotonated dimer involves in a strong H(π)···O interaction. It is important to compare the reliability of diverse computational approaches employed in quantum chemistry on the basis of the calculational results of this system. MP2 calculations using a small cc-pVDZ basis set give reasonable structures, but those using extended basis sets predict wrong π-stacked structures due to the overestimation of the dispersion energies of the π-π interactions. A few new DFT functionals with the empirical dispersion give reliable results consistent with the CCSD(T)/CBS results. The binding energies of the neutral, cationic protonated, and anionic deprotonated phenol dimers are estimated to be more than 28.5, 118.2, and 118.3 kJ mol(-1), respectively. The energy components of the intermolecular interactions for the neutral, protonated and deprotonated dimers are analyzed.     

Is Electrostatics Sufficient to Describe Hydrogen‐Bonding Interactions?

The stability and geometry of a hydrogen‐bonded dimer is traditionally attributed mainly to the central moiety A?H???B, and is often discussed only in terms of electrostatic interactions. The influence of substituents and of interactions other than electrostatic ones on the stability and geometry of hydrogen‐bonded complexes has seldom been addressed. An analysis of the interaction energy in the water dimer and several alcohol dimers—performed in the present work by using symmetry‐adapted perturbation theory—shows that the size and shape of substituents strongly influence the stabilization of hydrogen‐bonded complexes. The larger and bulkier the substituents are, the more important the attractive dispersion interaction is, which eventually becomes of the same magnitude as the total stabilization energy. Electrostatics alone are a poor predictor of the hydrogen‐bond stability trends in the sequence of dimers investigated, and in fact, dispersion interactions predict these trends better.     

Role of long-range intermolecular forces in the formation of inorganic nanoparticle clusters

An understanding of the role played by intermolecular forces in terms of the electron density distribution is fundamental to the understanding of the self-assembly of molecules in the formation of a molecular crystal. Using ab initio methods capable of describing both short-range intramolecular interactions and long-range London dispersion interactions arising from electron correlation, analyses of inorganic dimers of As(4)S(4) and As(4)O(6) molecules cut from the structures of realgar and arsenolite, respectively, reveal that the molecules adopt a configuration that closely matches that observed for the crystal. Decomposition of the interaction energies using symmetry-adapted perturbation theory reveals that both model dimers feature significant stabilization from electrostatic forces as anticipated by a Lewis acid/Lewis base picture of the interaction. London dispersion forces also contribute significantly to the interaction, although they play a greater role in the realgar structure near equilibrium than in arsenolite.     

An intermolecular perturbation approach to hydrogen bondings in systems in the ground and excited electronic states

The intermolecular interaction energy for binary systems in the ground and excited electronic states was partitioned into the Coulomb, exchange-repulsion, induction, dispersion and charge-transfer interaction terms by the perturbation expansion method. The various interaction terms were evaluated for the hydrogen bondings in (HF)₂, (H₂O)₂, (CH₃OH)₂, (RCOOH)₂, and HF·H₂O in various geometrical configurations. It has been found that the Coulombic interaction plays a dominant role in the stability of these hydrogen bonded systems. The method was further applied to the HCOOH·H₂O codimer in both the ground and excited singlet electronic states. The results were in accord with the well-known water solvent effects on the shifts of absorption spectral bands.