BeH2与HX(X=F,Cl, Br,I)形成的二氢键复合物的结构特征与本质

对BeH2与HX(X=F, Cl, Br, I)形成的二氢键复合物的结构特征及本质进行了探讨. 在MP2/6-311++G(3d,3p)水平优化、频率验证, 得到复合物的分子结构, 用分子间距离及电子密度拓扑理论确认BeH2与卤化氢已形成了二氢键型复合物. 在MP2/6-311++G(3d, 3p)水平下进行基函数重叠误差(BSSE)校正后的结合能在-14.468 kJ·mol-1到-5.464 kJ·mol-1之间.用对称匹配微扰理论(SAPT)对复合物的结合能进行分解, 结果表明, BeH2…HX二氢键复合物中静电能对总吸引能的贡献都是最主要的, 但交换排斥能、诱导能、色散能对总结合能的贡献也很重要. 从BeH2…HF到BeH2…HI, 诱导能对总吸引能的贡献从37.8%逐渐减小到24.0%. 而色散能对总吸引能的贡献从BeH2…HF体系中的16.0%逐渐增加到BeH2…HI体系中的33.8%.     

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.     

Estimated MP2 and CCSD(T) interaction energies of n-alkane dimers at the basis set limit: comparison of the methods of Helgaker et al. and Feller

The MP2 (the second-order M?ller-Plesset calculation) and CCSD(T) (coupled cluster calculation with single and double substitutions with noniterative triple excitations) interaction energies of all-trans n-alkane dimers were calculated using Dunning's [J. Chem. Phys. 90, 1007 (1989)] correlation consistent basis sets. The estimated MP2 interaction energies of methane, ethane, and propane dimers at the basis set limit [EMP2(limit)] by the method of Helgaker et al. [J. Chem. Phys. 106, 9639 (1997)] from the MP2/aug-cc-pVXZ (X=D and T) level interaction energies are very close to those estimated from the MP2/aug-cc-pVXZ (X=T and Q) level interaction energies. The estimated EMP2(limit) values of n-butane to n-heptane dimers from the MP2/cc-pVXZ (X=D and T) level interaction energies are very close to those from the MP2/aug-cc-pVXZ (X=D and T) ones. The EMP2(limit) values estimated by Feller's [J. Chem. Phys. 96, 6104 (1992)] method from the MP2/cc-pVXZ (X=D, T, and Q) level interaction energies are close to those estimated by the method of Helgaker et al. from the MP2/cc-pVXZ (X=T and Q) ones. The estimated EMP2(limit) values by the method of Helgaker et al. using the aug-cc-pVXZ (X=D and T) are close to these values. The estimated EMP2(limit) of the methane, ethane, propane, n-butane, n-pentane, n-hexane, n-heptane, n-octane, n-nonane, and n-decane dimers by the method of Helgaker et al. are -0.48, -1.35, -2.08, -2.97, -3.92, -4.91, -5.96, -6.68, -7.75, and -8.75 kcal/mol, respectively. Effects of electron correlation beyond MP2 are not large. The estimated CCSD(T) interaction energies of the methane, ethane, propane, and n-butane dimers at the basis set limit by the method of Helgaker et al. (-0.41, -1.22, -1.87, and -2.74 kcal/mol, respectively) from the CCSD(T)/cc-pVXZ (X=D and T) level interaction energies are close to the EMP2(limit) obtained using the same basis sets. The estimated EMP2(limit) values of the ten dimers were fitted to the form m0+m1X (X is 1 for methane, 2 for ethane, etc.). The obtained m0 and m1 (0.595 and -0.926 kcal/mol) show that the interactions between long n-alkane chains are significant. Analysis of basis set effects shows that cc-pVXZ (X=T, Q, or 5), aug-cc-pVXZ (X=D, T, Q, or 5) basis set, or 6-311G** basis set augmented with diffuse polarization function is necessary for quantitative evaluation of the interaction energies between n-alkane chains.     

Theoretical investigations on heteronuclear chalcogen-chalcogen interactions: on the nature of weak bonds between chalcogen centers

To understand the intermolecular interactions between chalcogen centers (O, S, Se, Te), quantum chemical calculations on model systems were carried out. These model systems were pairs of monomers of the composition (CH3)2X1 (X1 = O, S, Se, Te) as the donors and CH3X2Z (with X2 = O, S, Se, Te and Z = Me, CN) as the acceptors. The variation of X1, X2, and Z leads to 32 pairs with 8 homonuclear cases (X1 = X2 = O, S, Se, Te) and 24 heteronuclear cases (X1 not equal X2). The MP2/SDB-cc-pVTZ, 6-311G* level of theory was used to derive the geometrical parameters and the interaction energies of the model systems. The pairs with Z = CN (17-32) show a considerably higher interaction energy than the pairs with CH3 groups only (1-16). Natural bond orbital (NBO) analysis revealed that the interaction of the dimers 1, 2, 5, 6, 9, 10, 13, 14, 17, 21, 25, and 29 is mainly due to weak hydrogen bonding between methyl groups and chalcogen centers. These systems all contain hard chalcogen atoms as acceptors. For all other systems, the chalcogen-chalcogen interaction dominates. The one-electron picture of an interaction between the lone pair of the donor chalcogen atom and the chalcogen-carbon antibonding sigma* orbital serves as a model to qualitatively rationalize trends found in many of these systems. However, it has to be applied with some amount of skepticism. A detailed analysis based on symmetry-adapted perturbation theory (SAPT) reveals that induction and dispersion forces dominate and contribute to the bonding in each case. Hydrogen-bonded compounds involve bonding electrostatic contributions. Compounds dominated by chalcogen-chalcogen interactions exhibit bonding due to electrostatic interactions only if one of the chalcogen atoms involved is sulfur or oxygen.     

Intermolecular Weak Interactions in HTeXH Dimers (X=O,S, Se,Te): Hydrogen Bonds,Chalcogen–Chalcogen Contacts and Chiral Discrimination

A theoretical study of the HTeXH (X=O, S, Se and Te) monomers and homodimers was carried out by means of second‐order Møller‐Plesset perturbation theory (MP2) computational methods. In the case of monomers, the isomerization energy from HTeXH to H₂Te=X and H₂X=Te (X=O, S, Se, and Te) and the rotational transition‐state barriers were obtained. Due to the chiral nature of these compounds, homo and heterochiral dimers were found. The electron density of the complexes was characterized with the atoms‐in‐molecules (AIM) methodology, finding a large variety of interactions. The charge transfer within the dimers was analyzed by means of natural bond orbitals (NBO). The density functional theory‐symmetry adapted perturbation theory (DFT‐SAPT) method was used to compute the components of the interaction energies. Hydrogen bonds and chalcogen–chalcogen interactions were characterized and their influence analyzed concerning the stability and chiral discrimination of the dimers.     

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.     

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.     

Ion pair formation in water. association constants of bolaform, bisquaternary ammonium, electrolytes by chemical trapping

The first and second association constants, K1 and K2, for ion pair formation in aqueous 0.02-3.5 M solutions of bis(trimethyl)-alpha,omega-alkanediammonium halides with variable spacer lengths, 1-n-1 2X (n = 2-4, X = Cl, Br) and bolaform salts and for tetramethylammonium halides (TMAX, X = Cl, Br), K(TMAX), were determined by the chemical trapping method. Values for K(TMAX) are small, K(TMABr) = 0.83 M(-1) and K(TMACl) = 0.29 M(-1), in agreement with literature values. For the bolaform salts, K1 depends on spacer length and counterion type, ranges from 0.4 to 17 M(-1), is 2-10 times larger than K2, is larger for Br- than Cl-, and decreases by a factor of approximately 3 for Cl- and approximately 10 for Br- as n increases from 2 to 4. K2, for the formation of bolaform dihalide pair, is essentially the same as that for ion pair formation in TMAX solutions, i.e., K2 approximately K(TMAX). Values of K1 and K(TMABr) obtained from changes in 79Br line widths are in good agreement with those obtained by chemical trapping. The results are consistent with a thermodynamic model in which the ion association depends on the balance of the ion specific hydration free energies of cations and anions and their ion specific and hydration interactions in ion pairs. Spacer length dependent ion pairing by bolaform electrolytes, which are analogues of the headgroups and counterions of gemini amphiphiles, suggests a new model for the spacer length dependent sphere-to-rod transitions of gemini micelles. Neutral, but polar, headgroup-counterion pairs have a lower demand for hydration that free headgroups and counterions, and headgroup-counterion pair formation releases interfacial water into the bulk aqueous phase, permitting tighter amphiphile packing in rodlike micelles.     

SiH4与HX(X=F,Cl,Br,I)形成的二氢键复合物的结构特征及本质

本文就SiH4与HX形成的二氢键复合物的结构特征及本质进行了探讨。在MP2/6-311++G(3d,3p)水平优化、频率验证得到复合物的分子结构,通过分子间距离及电子密度等值线图,我们确认SiH4与卤化氢已形成了二氢键复合物。MP2/6-311++G(3d,3p)水平下进行BSSE校正后的结合能为2.703－4.439 KJ/mol。用对称匹配微绕理论(SAPT)对结合能进行分解,分解结果显示,SiH4匟X(X=F,Cl,Br,I)二氢键复合物中静电能对总吸引能的贡献小于28％,并且相对稳定,这就是说SiH4匟X二氢键复合物的本质并非静电作用,而是静电能、诱导能、色散能、交换能对总结合能的贡献都非常重要。     

Interactions in diatomic dimers involving closed-shell metals

Interaction energies of dimers containing alkaline earth (Be, Mg, and Ca) metals have been investigated using symmetry-adapted perturbation theory (SAPT) and supermolecular (SM) methods. Also, to enable broader comparisons, some calculations have been performed on the Zn dimer and on the He-Mg dimer. Although all of the investigated metallic atoms have closed electronic shells, the quasidegeneracy of the ground states of these atoms with the lowest-lying excited states leads to convergence problems in theories based on a single-determinant reference state. The main goal of the present work was to establish how the quality of the interaction energies computed using various electronic-structure methods changes across the range of atoms. We show that although the convergence problems become somewhat less severe with the increase of the atomic number, single-determinant-based methods do not provide reliable interaction energies for any of the investigated metallic dimers even at the level of the coupled-cluster method with single, double, and noniterative triple excitations [CCSD(T)]. However, interaction energies accurate to within a few percent can be obtained if CCSD(T) calculations in large basis sets are extrapolated to the complete basis set limit and followed by full configuration interaction (FCI) calculations with a frozen-core (FC) approximation. Since the systems considered contain only two valence electrons, FCI/FC calculations have been feasible for all of them except for Zn2, providing the best theoretical estimates of the binding energies to date. We found that a large part of the error of the SAPT results originates from limiting some exchange components to terms proportional to the squares of the intermonomer orbital overlap integrals. When the neglected terms were approximately accounted for, the accuracy improved significantly and became comparable to that of CCSD(T), allowing us to obtain for the first time a physical interpretation of the interaction energies in metallic dimers.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Physical nature of silane⋯carbene dimers revealed by state-of-the-art ab initio calculations

Using the SAPT2 + 3(CCD)δMP2 method in complete basis set (CBS) limit, it is shown that the interactions in the recently studied silane⋯carbene dimers are mainly dispersive in nature. Consequently, slow convergence of dispersion energy also forces slow convergence of the interaction energy. Therefore, obtaining very accurate values requires extrapolation of the correlation part to the CBS limit. The most accurate values obtained at the CCSD(T)/CBS level of theory show that the studied silane⋯carbene dimers are rather weakly bound, with interaction energies ranging from about −1.9 to −1.3 kcal/mol. Comparing to CCSD(T)/CBS, it will be shown that SCS-MP2 and MP2C methods clearly underestimate and methods based on SAPT2+ and having some third-order corrections, as well as the MP2 method, overestimate values of interaction energies. Popular SAPT(DFT) method performs better than SCS-MP2 and MP2C; nevertheless, underestimation is still considerable. The underestimation is slightly quenched if third-order dispersion energy and its exchange counterpart is added to the SAPT(DFT). The closest value of CCSD(T)/CBS has been given by the SAPT2 + (3)(CCD)δMP2 method in quadruple-ζ basis set. © 2019 Wiley Periodicals, Inc.     

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…     

Symmetry-adapted perturbation theory analysis of the N...HX hydrogen bonds

The main aim of the study was the detailed investigation of the interaction energy decomposition in dimers and trimers containing N...HX bonds of different types. The study of angular dependence of interaction energy terms partitioned according to the symmetry-adapted perturbation theory (SAPT) was performed for the dimers containing N...HX bonds as mentioned above: ammonia-HX (X = F, Cl, Br) and pyridine-HF complexes. It was found that the electrostatic and induction terms exhibit strong angular dependence, while the exchange contributions are less affected. The dispersion terms are virtually nondirectional. In addition, the three-body SAPT interaction energy analysis for the mixed acid-base NH3...(HF)2 and (NH3)2...HF trimers revealed strong differences between interactions of similar strength but different types (i.e., hydrogen bond and general electrostatic interaction). The importance of the induction terms for the nonadditivity of the interaction energy in strongly polar systems was confirmed.     

Spectroscopic determination of dipole-dipole interaction in solutions of trioctylammonium halides

Dimerization energies of trioctylammonium halides in carbon tetrachloride solution, associated as a result of dipole-dipole interaction, were calculated from the infrared absorption spectra. The values obtained were 5, 4, and 1.5 kcal/mole for the chloride, bromide, and iodide, respectively. The calculation showed that the dimers of these salts are contact dimers.     

Inter-Alkali-Metal Dative Bond in the MMN3− (M=Alkali Metal) Cluster

Alkali metals are generally Lewis acids. On the contrary, Lewis basic character of alkali metals forming donor – acceptor complexes is a very rare phenomenon. In this contribution, I have theoretically designed an anionic cluster MMN₃⁻ (M=alkali metals) on the basis of experimentally known reagent, alkali salt of azide ion MN₃, which shows unprecedented M:⁻→M donor-acceptor interaction. To the best of author's knowledge, the characterization of such donor-acceptor interaction among alkali metals is unprecedented. Formation of the 2c–2e donor-acceptor bonds have been confirmed by quantum theory of atoms in molecules and electron localization function analyses. The calculated bond dissociation energies are significant suggesting their possible spectroscopic identification.     

Theoretical studies on vibrational spectra of some halides of Group IIB elements

The vibrational spectra of Group IIB elements halides MX2 and their dimers M2X4 (M=Zn(II), Cd(II) and Hg(II); X=F, Cl, Br and I) have been systematically investigated by ab initio RHF and B3LYP methods with LanL2MB, LanL2DZ and SDD basis sets. The optimized geometries, calculated vibrational frequencies are evaluated via comparison with the experimental data. The vibrational frequencies, calculated by these methods with different basis sets, are compared to each other too. The best results can be obtained by RHF/SDD method, with this method, the deviations for MX2 and Hg2X4 are <7%. Some vibrational frequencies of M2X4 that have not been experimentally reported are also predicted.     

Effects of various halogen anions and cations of alkali metals on energetics of excess charge recombination in stilbene donor-acceptor capped DNA hairpins

DNA hairpin conjugates with a stilbenedicarboxamide (Sa) hole donor and a stilbenediether (Sd) hole acceptor are considered as model systems for studying charge recombination (CR) of excess charges in DNA. Using the method of thermodynamic integration, we estimated the relative free energies of this process in hairpins with three adenine:thymine pairs between Sa and Sd surrounded by 1 M aqueous solutions of ionic compounds M(+)Cl(-) (M = Li, Na, K) and Na(+)X(-) (X = F, Cl, Br, I). The values of this quantity were calculated with respect to the free energy for the same hairpin in the 1 M NaCl aqueous solution. Based on the results obtained, we conclude that halogen anions have no significant influence on the rate of the CR reaction. By contrast, cations of other alkali metals can considerably change the potential barrier of the process, thus affecting the reaction rate. Different results obtained for cations and anions were attributed to the fundamental distinction in the electrostatic interactions of M(+) and X(-) species with negatively charged phosphate groups of the hairpin. In addition, our results show that the relative free energy of CR is larger for cations that are able to be closer to Sd and Sa structural units. The latter correlation suggests that the replacement of Na(+) by cations of other alkali metals enables one to change the CR rate modifying it in either direction.     

Hydration energies in the gas phase of select (MX)mM+ ions, where M+ = Na+, K+, Rb+, Cs+, NH4+ and X- = F-, Cl-, Br-, I-, NO2-, NO3-. Observed magic numbers of (MX)mM+ ions and their possible significance

The sequential hydration energies and entropies with up to four water molecules were obtained for MXM(+) = NaFNa(+), NaClNa(+), NaBrNa(+), NaINa(+), NaNO(2)Na(+), NaNO(3)Na(+), KFK(+), KBrK(+), KIK(+), RbIRb(+), CsICs(+), NH(4)BrNH(4)(+), and NH(4)INH(4)(+) from the hydration equilibria in the gas phase with a reaction chamber attached to a mass spectrometer. The MXM(+) ions as well as (MX)(m)M(+) and higher charged ions such as (MX)(m)M(2)(2+) were obtained with electrospray. The observed trends of the hydration energies of MXM(+) with changing positive ion M(+) or the negative ion X(-) could be rationalized on the basis of simple electrostatics. The most important contribution to the (MXM-OH(2))(+) bond is the interaction of the permanent and induced dipole of water with the positive charge of the nearest-neighbor M(+) ion. The repulsion due to the water dipole and the more distant X(-) has a much smaller effect. Therefore, the bonding in (MXM-OH(2))(+) for constant M and different X ions changes very little. Similarly, for constant X and different M, the bonding follows the hydration energy trends observed for the naked M(+) ions. The sequential hydration bond energies for MXM(H(2)O)(n)(+) decrease with n in pairs, where for n = 1 and n = 2 the values are almost equal, followed by a drop in the values for n = 3 and n = 4, that again are almost equal. The hydration energies of (MX)(m)M(+) decrease with m. The mass spectra with NaCl, obtained with electrospray and observed in the absence of water vapor, show peaks of unusually high intensities (magic numbers) at m = 4, 13, and 22. Experiments with variable electrical potentials in the mass spectrometer interface showed that some but not all of the ion intensity differentiation leading to magic numbers is due to collision-induced decomposition of higher mass M(MX)(m)(+) and M(2)(MX)(m)(2+) ions in the interface. However, considerable magic character is retained in the absence of excitation. This result indicates that the magic ions are present also in the saturated solution of the droplets produced by electrospray and are thus representative of particularly stable nanocrystals in the saturated solution. Hydration equilibrium determinations in the gas phase demonstrated weaker hydration of the magic ion (NaCl)(4)Na(+).