H-bonded complexes of aniline with HF/F- and anilide with HF in terms of symmetry-adapted perturbation, atoms in molecules, and natural bond orbitals theories

The hydrogen-bonded isoelectronic complexes of aniline with HF/F- and an ionic form of aniline with HF were investigated by use of computational methods: Symmetry-Adapted Perturbation Theory (SAPT), Atoms in Molecules (AIM), and Natural Bond Orbitals (NBO) approaches. All computations were based on structural models previously generated at the B3LYP/6-311+(d,p) level. The differences between neutral (Ph-NH2...HF)and anionic (Ph-NH2...F- and Ph-NH-...HF) complexes were clearly outlined. The discussed charged complexes serve as Lewis acids and base, HF and F-, respectively. It was found that electrostatic and induction energy terms, obtained as a result of the SAPT method, are most dependent on the type of H-bonding (i.e.,charged or neutral). The electrostatic term is the most distinctive between the neutral and charge-assisted hydrogen bonds in the investigated two-body systems, whereas the latter is more significant in the case of weaker interactions (larger H...B distances). Application of Principal Component Analysis (PCA) to energy components obtained from the SAPT procedure indicated that all of them are relatively well intercorrelated.The above-mentioned terms together with the exchange energy terms are the most important contributions ofthe main principal component, which describes 95% of the total variance. Comparison of AIM parameters in bond critical points for modeled H-bond systems shows a good agreement with those from equilibrium complexes, both experimental and calculated ones. It was found that charged H-bonded complexes exhibit larger fluctuation of electron density and its Laplacian in bond critical points, in line with SAPT analysis. NBO results confirmed the effect of the strength of interaction on property changes both in the region of H-bonding and outside of it. The latter, more distant consequences follow the Bent-Walsh rule for all studied complexes.     

HCN(HNC)与NH3, H2O和HF分子间相互作用的理论研究

在MP2/aug-cc-pVTZ水平上, 对HCN(HNC)与NH3, H2O和HF分子间可能存在的氢键型复合物进行了全自由度能量梯度优化, 通过在相同水平上的频率验证分析发现了稳定的分子间相互作用形式是HCN(HNC)作为质子供体或作为质子受体形成的复合物. 基组重叠误差对总相互作用能的影响均小于3.34 kJ/mol. 通过自然键轨道(NBO)分析, 研究了单体和复合物中的原子电荷和电荷转移对分子间相互作用的影响. 对称性匹配微扰理论(SAPT, Symmetry Adapted Perturbation Theory)能量分解结果表明, 在分子间相互作用中, 静电作用与诱导作用占主导地位, 而诱导作用与复合物的电荷转移之间具有良好的正相关性.     

传统氢键(X—H…Y)还是X—H…π相互作用？ ——无机苯与卤化氢相互作用本质的量子化学研究

在MP2/6-311＋＋G**水平对无机苯(B3N3H6)与卤化氢HX (X＝F, Cl, Br, I)相互作用体系进行了系统研究. 结果表明在B3N3H6-HX (X＝F, Cl, Br, I)体系的平衡几何结构中, HX的H原子倾向于指向B3N3H6环上的N原子, 且从HF到HI相互作用强度依次减弱. 与苯-卤化氢体系比较, 除与HF相互作用B3N3H6较C6H6强外, 其余体系B3N3H6均较C6H6弱(结合能数值相差4 kJ/mol左右). 对称匹配微扰理论(SAPT)能量分解结果说明静电、诱导和色散力对描述B3N3H6-卤化氢体系的相互作用都很重要, 从HF到HI静电能占总吸引作用能的百分比逐渐减少, 色散能占总吸引作用能的百分比逐渐增加, 这种变化趋势与苯-卤化氢体系比较类似, 表明B3N3H6与卤化氢的相互作用随着卤素原子序数的递增, 传统氢键作用趋势减弱, X—H…π相互作用趋势增强.     

Theoretical study of properties of H bonds and intermolecular interactions in linear cis-,trans-cyclotriazane clusters (n = 2-8)

Our calculations based upon Becke's three-parameter functional of density-functional theory (DFT) with the correlation of Lee, Yang, and Parr (B3LYP), natural bond orbital, and atoms in molecule indicate that in drastic contrast to most H-bonded systems, the anticooperative and cooperative effects coexist in the linear H-bonded cis-,trans (c,t)-cyclotriazane clusters (n = 2-8). As cluster size increases, the properties along the H-bonded chains at trans-positions take on the unexpectedly anticooperative changes which are reflected in elongation of the N...H hydrogen bonds, frequency blueshift in the N-H stretching vibrations, decay in the n(N)-->sigma*(N-H) charge transfers, and weakening of strengths of the N...H bonds. And the cooperative changes in the corresponding properties for the cis- H-bonded chains are observed to be concurrent with the anticooperativities. The rise and fall in the n(N)-->sigma*(N-H) interactions cause increment and decrement in capacities of the clusters to concentrate electrons at the bond critical points of the N...H bonds, and thereby leading to the cooperative and the anticooperative changes especially in the N...H lengths and the N-H stretching frequencies. In terms of three-body symmetry-adapted perturbation theory (three-body SAPT), the first exchange nonadditivity plays a more important role in stabilizing trimer than the nonadditive induction. However, the dominance of the first exchange nonadditivity in three-body interaction unexpectedly triggers the anticooperative effect that counteracts the concurrent cooperative effect. According to the SAPT(DFT), which is a combination of SAPT with asymptotically corrected DFT, DFT/B3LYP is able to succeed in describing the electrostatic, exchange, and induction components, but fails to yield satisfactory interaction energies due to the fact that about 40% of short-range dispersion energy is neglected by the DFT, which is different from many H-bonded described well by the DFT. A quantum cluster equilibrium model illustrates that the c,t-cyclotriazane liquid phase exhibits a weak cooperative effect.     

The properties of weak and strong dihydrogen-bonded D-H...H-A complexes.

The properties of six dihydrogen-bonded (DHB) dimers with the BeH2 molecule as a proton acceptor were calculated by MP2, CCSD(T) and B3LYP methods. The structural, energetic and spectroscopic parameters are presented and analyzed in terms of their possible correlation with the interaction energy and the intermolecular H...H separation. The symmetry-adapted perturbation theory (SAPT) calculations were performed to gain more insight into the nature of the H...H interactions. The studied complexes are divided into three groups based on the calculated intermolecular distances and the interaction energies which range from approximately -1 to -42 kJ mol(-1). The analysis of the interaction energy components indicates that, in contrast to conventional hydrogen bonds, the induction energy is the most important term in the BeH2NH4+ complex. On the other hand, there is no sharp boundary between the DHB complexes classified as hydrogen bonded and van der Waals systems. The complexation-induced changes in vibrational frequencies and in proton shielding constants show a relationship with the interaction energy. The values of the 2hJXH and 3hJBeX coupling constants correlate well with the interaction energy and with the intermolecular distance.     

Fluorine as a Lewis acid: A symmetry-adapted perturbation theory based on density functional theory and interacting quantum atoms study of noncovalent interactions in the NCF⋯NH3 complex

Different computational methods are used to investigate the nature of interaction in the NCF⋯NH₃ model complex, in which the fluorine atom acts as a Lewis acid and forms a noncovalent bond with the ammonia (Lewis base). Symmetry-adapted perturbation theory based on density functional theory (SAPT(DFT)) indicates that the noncovalent interaction in the NCF⋯NH₃ complex is mainly electrostatics. However, dispersion and induction terms also play important roles. Although fluorine noncovalent interactions are typically classified as halogen bonds, they are somewhat different from the well-known halogen bonds of iodine, bromine, and chlorine. The halogen bonds of NCCl⋯NH₃ and NCBr⋯NH₃ are directional and the C  X  N (X = Cl or Br) angle tends to be linear. In contrast, the fluorine interaction in NCF⋯NH₃ is not directional; the interaction energy shows no sensitivity to the angular (C  F  N ) distortions, and the energy profile is flat over a wide angular range (from 180° to about 140° ). However, for the angles less than 130° , the energy curve shows a clear angular dependence and the interaction between NCF and NH₃ becomes stronger as the C  F  N angle decreases. It seems that at the tighter angles, a tetrel-bonded NCF⋯NH₃ complex is preferred. Moreover, interacting quantum atoms (IQA) analysis shows that the competition between different intra-atomic and interatomic interactions determines the geometry of NCF⋯NH₃ complex.     

Assessing the Accuracy of SAPT(DFT) Interaction Energies by Comparison with Experimentally Derived Noble Gas Potentials and Molecular Crystal Lattice Energies

The density functional version of symmetry‐adapted perturbation theory, SAPT(DFT), is a computationally efficient method for calculating intermolecular interaction energies. We evaluate its accuracy by comparison with experimentally determined noble gas interaction potentials and sublimation enthalpies, most of which have not been previously calculated using this method. In order to compare the results with wavefunction methods, we also calculate these quantities using MP2 and, for noble gas dimers, using CCSD(T). For the crystal lattice energy calculations, we include corrections to the dispersion, electrostatic, and induction energies that account for the finite interaction distance cutoff and higher‐order induction contributions. Overall, the energy values extrapolated to the complete basis set limit show that SAPT(DFT) achieves significantly better agreement with experiment than MP2.     

Substituent effects on CL···N, S···N, and P···N noncovalent bonds

Cl, S, and P atoms have previously been shown as capable of engaging in a noncovalent bond with the N atom on another molecule. The effects of substituents B on the former atoms on the strength of this bond are examined, and it is found that the binding energy climbs in the order B = CH(3) < NH(2) < CF(3) < OH < Cl < NO(2) < F. However, there is some variability in this pattern, particularly for the NO(2) group. The A···N bonds (A = Cl, S, P) can be quite strong, amounting to as much as 10 kcal/mol. The binding energy arises from approximately equal contributions from its induction and electrostatic components, although the former becomes more dominant for the stronger bonds. The induction energy is due in large measure to the transfer of charge from the N lone pair to a B-A σ* antibonding orbital of the electron-acceptor molecule containing Cl, S, or P. These A···N bonds typically represent the lowest-energy structure on each potential energy surface, stronger than H-bonds such as NH···F, CH···N, or SH···N.     

Structure,properties, and nature of the BrF‐HX complexes: An Ab initio study

Structure and properties of complexes (energies and charge transfer) of complexes BrF‐HX (X = F, Cl, Br, I) have been investigated at the MP2/aug‐cc‐pVDZ (aug‐cc‐pVDZ‐pp basis sets for I) level. Two types of geometries (hydrogen‐bonded and halogen‐bonded) are observed. The calculated interaction energies show that the halogen bonded structures are more stable than the corresponding hydrogen‐bonded structures. To study the nature of the intermolecular interactions, symmetry‐adapted perturbation theory (SAPT) energy decomposition analysis reveals that the BrF‐HX complexes are dominantly electrostatic in nature. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Orbital relaxation and the third-order induction energy in symmetry-adapted perturbation theory

Theoretical investigations of the induction interaction between closed-shell molecules which fully account for the orbital relaxation effects are presented. Explicit expressions for the third-order induction energy in terms of molecular integrals and orbital energies are given and implemented within the sapt2008 program for symmetry-adapted perturbation theory (SAPT) calculations. Numerical investigations for the He–He, He–LiH, Ar–Ar, H₂–CO, H₂O–H₂O, and H₂O–NH₃ model dimers show that the orbital relaxation increases the third-order induction interaction by 15 to 50% at near-equilibrium geometries, with the largest effect observed for complexes involving highly polar monomers. At large intermonomer separations, the relaxed third-order induction energy perfectly recovers the difference $\delta E^{\rm HF}_{\rm int}$ between the Hartree–Fock interaction energy and the sum of the uncorrelated SAPT contributions through second order in the intermolecular interaction operator. At the near-equilibrium geometries, the sum of the relaxed third-order induction and exchange-induction energies reproduces, however, only a small fraction (6 to 15%) of $\delta E^{\rm HF}_{\rm int}$ for the nonpolar systems and about 40 to 60% for the polar ones. A comparison of the complete SAPT calculations with the coupled-cluster treatment with single, double, and noniterative triple excitations [CCSD(T)] suggests that the pure SAPT approach with all the available third-order corrections is more accurate for nonpolar systems while for the polar ones the hybrid approach including $\delta E^{\rm HF}_{\rm int}$ gives better results.     

Correlated intermolecular interaction components from asymptotically corrected Kohn-Sham orbitals

In these years there was considerable interest inunderstanding of intermolecular forces in energetic(explosive) systems[1—3]. The supermolecular approach(SM) is widely adopted for calculating ab initio in-termolecular interactions. Nevertheless, it is unable toprovide physically meaningful interaction contribu-tions such as electrostatic, induction, repulsion anddispersion energies. In contrast, the symmetry-adaptedperturbation theory (SAPT)[4—8] has the ability to de-rive these correlated…     

Cooperative effects of noncovalent bonds to the Br atom of halogen-bonded H(3)N[ellipsis (horizontal)]BrZ and HCN[ellipsis (horizontal)]BrZ (Z = F, Br) complexes

A series of complexes formed between halogen-bonded H(3)N∕HCN[ellipsis (horizontal)]BrZ (Z = Br, F) dimers and H(3)N∕HCN[ellipsis (horizontal)]BrZ[ellipsis (horizontal)]XY (XY = HF, ClF, BeH(2), LiF) trimers were investigated at the MP2 and B3LYP levels of theory using a 6-31++G(d,p) basis set. Optimized structures, interaction energies, and other properties of interest were obtained. The addition of XY to the H(3)N∕HCN[ellipsis (horizontal)]BrZ dyad leads to enhanced intermolecular binding with respect to the isolated monomers. This enhanced binding receives contributions from the electrostatic and inductive forces between the constituent pairs, with, in some instances, substantial three-body non-additive contributions to the binding energy. It was found that the XY = LiF interaction causes the greatest distortion of the H(3)N∕HCN[ellipsis (horizontal)]BrZ halogen bond from the preferred linear orientation and also provides the strongest binding energy via the nonadditive energy.     

Application of symmetry-adapted perturbation theory to small ionic systems

The application of symmetry-adapted perturbation theory (SAPT) to small ionic systems was investigated in the context of the accuracy of calculated interaction energies for alkali halides. Two forms of alkali halides were considered: ion pairs M(+)X(-) (M = Li, Na, K, Rb, and X = F, Cl, Br, I) and dimers (MX)(2). The influence of the order of energy correction terms included in SAPT and the effect of the so-called hybrid approach to SAPT on the accuracy of the calculated energies (such as the interaction energies in the ion pairs and the binding energies in the dimers with respect to two free monomers) were studied. The effects of the size of basis sets, combined with SAPT, on the accuracy were also established.     

Monte Carlo simulation of liquid hydroxylamine using ab initio intermolecular potential functions

A potential model for intermolecular interactions between hydroxylamine (NH₂OH) molecules based on ab initio quantum mechanical calculations is reviewed and critically assessed by analyzing results from a Monte Carlo simulation of liquid hydroxylamine. The liquid structure is studied in detail using radial, energy, and angular distribution functions, coordination numbers, and their distribution. Results indicate a large first solvation shell (5.3 Å), which contains 13 molecules, out of which only 4 are truly bonded by nonlinear, low-energy hydrogen bonds. These are of either the OH…O or the OH…N type, as NH…O and NH…N linear bonds are considerably suppressed, and no cyclic dimers are found. The dependence of the structural and physical properties on the simulation characteristics has also been investigated.     

Intermolecular potentials based on symmetry-adapted perturbation theory with dispersion energies from time-dependent density-functional calculations

Recently, three of us have proposed a method [Phys. Rev. Lett. 91, 33201 (2003)] for an accurate calculation of the dispersion energy utilizing frequency-dependent density susceptibilities of monomers obtained from time-dependent density-functional theory (DFT). In the present paper, we report numerical calculations for the helium, neon, water, and carbon dioxide dimers and show that for a wide range of intermonomer separations, including the van der Waals and short-range repulsion regions, the method provides dispersion energies with accuracies comparable to those that can be achieved using the current most sophisticated wave-function methods. If the dispersion energy is combined with (i) the electrostatic and first-order exchange interaction energies as defined in symmetry-adapted perturbation theory (SAPT) but computed using monomer Kohn-Sham (KS) determinants, and (ii) the induction energy computed using the coupled KS static response theory, (iii) the exchange-induction and exchange-dispersion energies computed using KS orbitals and orbital energies, the resulting method, denoted by SAPT(DFT), produces very accurate total interaction potentials. For the helium dimer, the only system with nearly exact benchmark values, SAPT(DFT) reproduces the interaction energy to within about 2% at the minimum and to a similar accuracy for all other distances ranging from the strongly repulsive to the asymptotic region. For the remaining systems investigated by us, the quality of the SAPT(DFT) interaction energies is so high that these energies may actually be more accurate than the best available results obtained with wave-function techniques. At the same time, SAPT(DFT) is much more computationally efficient than any method previously used for calculating the dispersion and other interaction energy components at this level of accuracy.     

Symmetry-adapted perturbation-theory calculations of intermolecular forces employing density-functional description of monomers

A symmetry-adapted perturbation theory based on Kohn-Sham determinants [SAPT(KS)] and utilizing asymptotically corrected exchange-correlation potentials has been applied to the He2, Ne2, (H2O)2, and (CO2)2 dimers. It is shown that SAPT(KS) is able to recover the electrostatic, first-order exchange, second-order induction, and exchange-induction energies with an accuracy approaching and occasionally surpassing that of regular SAPT at the currently programmed theory level. The use of the asymptotic corrections is critical to achieve this accuracy. The SAPT(KS) results can be obtained at a small fraction of the time needed for regular SAPT calculations. The robustness of the SAPT(KS) method with respect to the basis set size is also demonstrated. A theoretical justification for high accuracy of SAPT(KS) predictions for the electrostatic, first-order exchange, and second-order induction energies has been provided.     

Can Dispersion Forces Govern Aromatic Stacking in an Organic Solvent?

Experimental support for the dominance of van der Waals dispersion forces in aromatic stacking interactions occurring in organic solution is surprisingly limited. The size‐dependence of aromatic stacking in an organic solvent was examined. The interaction energy was found to vary by about 7.5 kJ mol⁻¹ on going from a phenyl–phenyl to an anthracene–pyrene stack. Strikingly, the experimental data were highly correlated with dispersion energies determined using symmetry‐adapted perturbation theory (SAPT), while the induction, exchange, electrostatic, and solvation energy components correlated poorly. Both the experimental data and the SAPT‐dispersion energies gave high‐quality correlations with the change in solvent accessible area upon complexation. Thus, the size‐dependence of aromatic stacking interactions is consistent with the dominance of van der Waals dispersion forces even in the presence of a competing polarizable solvent.     

Can Dispersion Forces Govern Aromatic Stacking in an Organic Solvent?

Experimental support for the dominance of van der Waals dispersion forces in aromatic stacking interactions occurring in organic solution is surprisingly limited. The size‐dependence of aromatic stacking in an organic solvent was examined. The interaction energy was found to vary by about 7.5 kJ mol^?1 on going from a phenyl–phenyl to an anthracene–pyrene stack. Strikingly, the experimental data were highly correlated with dispersion energies determined using symmetry‐adapted perturbation theory (SAPT), while the induction, exchange, electrostatic, and solvation energy components correlated poorly. Both the experimental data and the SAPT‐dispersion energies gave high‐quality correlations with the change in solvent accessible area upon complexation. Thus, the size‐dependence of aromatic stacking interactions is consistent with the dominance of van der Waals dispersion forces even in the presence of a competing polarizable solvent.     

Two-photon absorption of hydrogen-bonded octupolar molecule clusters

Charge-transfer octupolar molecules can form clusters in solution through intermolecular hydrogen bonds. In the present work we explore the role of such clustering on two-photon absorption (TPA) spectra assuming 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as a model system. Using density functional quadratic response theory we examine different cluster structures of TATB dimers, trimers, and tetramers taken from snapshots of molecular dynamics simulations. In comparison with the TPA spectrum of a monomer, significant red shifts of charge-transfer states are predicted for all chosen clusters, which mainly is the result of the distortion of the structures induced by the aggregation. The TPA spectra for dimers and trimers show strong conformation dependence, whereas they turn out to be more stable for tetramers. Enhancements of TPA absorption have also been found for clusters containing less distorted molecules connected by hydrogen bonds.