Theoretical investigation of the weakly dihydrogen bonded dimers H2 − nXnAlHHArF and H2 − nXnAlHHKrF (n = 0–2; X = F, Cl)

An ab initio computational study of the properties of the dihydrogen-bonded complexes of H_{2 − n}X_nAlH (n = 0–2; X = F, Cl) with the rare gas (Rg) compounds HArF and HKrF was carried out at the MP2(full)/6-311++G(2d,2p) level of theory. For all the studied complexes, we found a low zero-point corrected binding energies. Large red shifts of the H–Rg vibrational stretching frequency in both complexes were predicted. Electrostatic interactions between the individual monomers are also predicted.     

Intermolecular complexes of HArF and HKrF with CO

Stable linear weakly bound hydrogen-bonded complexes of HArF and HKrF with the CO molecule have been predicted by ab initio computations at the MP2/6- 311+ +G(2d,2p) level of theory. The complexes, having stabilities in the order, FArH...CO>FKrH...CO>FArH...OC>FKrH...OC are compared. They exhibit unusual blueshifts of the Ar-H (Kr-H) stretching frequency, as well as contraction of the Ar-H (Kr-H) bond, and these effects are rationalized mainly in terms of shifts in the electron density of HArF (HKrF) on complexation, caused mainly by a combination of the intermolecular electrostatic interaction, electron-electron (Pauli) and nuclear-nuclear repulsion and charge density transferred from the CO molecule to the rare-gas-containing molecule.     

Experimental and computational study of HXeY...HX complexes (X, Y = Cl and Br): an example of exceptionally large complexation effect

The complexes of xenon hydrides HXeY (Y = Cl and Br) with hydrogen halides HX (X = Cl and Br) have been studied both computationally and experimentally in a xenon matrix. The experiments revealed three new complexes: HXeBr...HBr, HXeBr...HCl, and HXeCl...HCl. The experimental assignments were done on the basis of the strong H-Xe stretching absorption of HXeY (Y = Cl and Br) molecules and supported by theoretical results. We experimentally obtained monomer-to-complex blue-shifts of this vibrational mode for all the studied systems (up to approximately 150 cm (-1)). The electronic structure calculations revealed three local structures for each HNgY...HX complexes and their computed interaction energies varied between -460 and -2800 cm (-1). The computational estimates of the vibrational shifts were in agreement with the experimental values. We also found possible experimental absorption belonging to HXeBr...(HBr) 2 trimer and its vibrational shift (+245 cm (-1)) is similar to the computational estimate of a cyclic ternary complex (+252 cm (-1)).     

How short can the H...H intermolecular contact be? New findings that reveal the covalent nature of extremely strong interactions

Ab initio calculations at the MP2/6-311++G(d,p) and MP2/aug-cc-pVDZ//MP2/aug-cc-pVTZ levels have been performed for the following complexes: H2OH+...HBeH, H2OH+...HBeBeH, H2OH+...HBeF, HClOH+...HBeH, Cl2OH+...HBeH, and Cl2OH+...HBeF. For all dimers considered, extremely short H...H intermolecular contacts (1.0-1.3 A) were obtained. These are the shortest intermolecular distances which have ever been reported, with binding energies within the range of 13.7-24.3 kcal/mol (MP2/aug-cc-pVDZ//MP2/aug-cc-pVTZ level). The interaction energies of the complexes analyzed were also extrapolated to the complete basis set (CBS) limit. To explain the nature of such strong interactions, the Bader theory was applied, and the characteristics of the bond critical points (BCPs) were analyzed. It was pointed out that for the major part of the H...H contacts considered here the Laplacian of the electron density at H...H BCP is negative indicating the partly covalent nature of such a connection. The term "covalent character of the hydrogen bond" used sometimes in recent studies is discussed. An analysis of the interaction energy components for dihydrogen bonded systems considered indicates that in contrast to conventional hydrogen bonded systems the attractive electrostatic term is outweighed by the repulsive exchange energy term and that the higher order delocalization energy term is the most important attractive term.     

Infrared laser spectroscopy of CH3...HF in helium nanodroplets: The exit-channel complex of the F + CH4 reaction

High-resolution infrared laser spectroscopy is used to study the CH3...HF and CD3...HF radical complexes, corresponding to the exit-channel complex in the F + CH4 --> HF + CH3 reaction. The complexes are formed in helium nanodroplets by sequential pickup of a methyl radical and a HF molecule. The rotationally resolved spectra presented here correspond to the fundamental v = 1 <-- 0 H-F vibrational band, the analysis of which reveals a complex with C(3v) symmetry. The vibrational band origin for the CH3...HF complex (3797.00 cm(-1)) is significantly redshifted from that of the HF monomer (3959.19 cm(-1)), consistent with the hydrogen-bonded structure predicted by theory [E. Ya. Misochko et al., J. Am. Chem. Soc. 117, 11997 (1995)] and suggested by previous matrix isolation experiments [M. E. Jacox, Chem. Phys. 42, 133 (1979)]. The permanent electric dipole moment of this complex is experimentally determined by Stark spectroscopy to be 2.4+/-0.3 D. The wide amplitude zero-point bending motion of this complex is revealed by the vibrational dependence of the A rotational constant. A sixfold reduction in the line broadening associated with the H-F vibrational mode is observed in going from CH3...HF to CD3...HF. The results suggest that fast relaxation in the former case results from near-resonant intermolecular vibration-vibration (V-V) energy transfer. Ab initio calculations are also reported (at the MP2 level) for the various stationary points on the F + CH4 surface, including geometry optimizations and vibrational frequency calculations for CH3...HF.     

A semi-empirical and ab initio combined approach for the full conformational searches of gaseous lysine and lysine–H2O complex

A full structural search of the canonical, zwitterionic, protonated and deprotonated lysine conformers in gas phase is presented. A total of 17,496 canonical, 972 zwitterionic, 11,664 protonated and 1458 trial deprotonated structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were initially optimized at the AM1 level, and the resulting structures were determined at the B3LYP/6-311G^* level. A total of 927 canonical, 730 protonated and 193 deprotonated conformers were found, but there were no stable zwitterionic structures in the gas phase. The most stable conformers of the canonical, protonated and deprotonated lysine were further optimized at the B3LYP/6-311++G^** level. The energies of the most stable structures were determined at the MP2/6-311G(2df,p) level and the vibrational frequencies were calculated at the B3LYP/6-311++G^** level. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, vertical ionization energies, enthalpies, Gibbs free energies and conformational distributions of gaseous lysine were presented. Numerous new structures are found and the lowest-energy lysine conformer is more stable than the existing one by 1.1 kcal/mol. Hydrogen bonds are classified and may cause significant red-shifts to the associated vibrational frequencies. The calculated proton affinity/dissociation energy and gas-phase basicity/acidity are in good agreement with the experiments. Calculations are also presented for the canonical lysine–H₂O and zwitterionic lysine–H₂O clusters. Interaction between lysine and H₂O significantly affects the relative conformational stabilities. Only one water molecule is sufficient to produce the stable zwitterionic structures in gas phase. The lowest-energy structure is found to be zwitterions when applying the conductor-like polarized continuum solvent model (CPCM) to the lysine–H₂O complexes.     

New scale factors for harmonic vibrational frequencies using the B3LYP density functional method with the triple-zeta basis set 6-311+G(d,p)

We have calculated optimal frequency scaling factors for the B3LYP/ 6-311+G(d,p) method for fundamental vibrational frequencies on the basis of a set of 125 molecules. Using the new scaling factor, the vibrational frequencies calculated with the triple-zeta basis set 6-311+G(d,p) give significantly better accuracy than those calculated with the double-zeta 6-31G(d) basis set. Scale factors were also determined for low-frequency vibrations using the molecular set of 125 molecules and for zero-point energies using a smaller set of 40 molecules. We have studied the effect on the calculated vibrational frequencies for various combinations of diffuse and polarization functions added to the triple-zeta 6-311G basis set. The 6-311+G(d,p) basis set is found to give almost converged frequencies for most molecules, and we conclude that our optimum scaling factors are valid for the basis sets 6-311G(d,p) to 6-311++G(3df,3pd). The new scale factors are 0.9679 for vibrational frequencies, 1.0100 for low-frequency vibrations, and 0.9877 for zero-point vibrational energies.     

Probes of the metal-to-ligand charge-transfer excited states in ruthenium-Am(m)ine-bipyridine complexes: the effects of NH/ND and CH/CD isotopic substitution on the 77 K luminescence

The effects of ligand perdeuteration on the metal-to-ligand charge-transfer (MLCT) excited-state emission properties at 77 K are described for several [Ru(L)(4)bpy](2+) complexes in which the emission process is nominally [uIII,bpy-] --> [RuII,bpy]. The perdeuteration of the 2,2'-bipyridine (bpy) ligand is found to increase the zero-point energy differences between the ground states and MLCT excited states by amounts that vary from 0 +/- 10 to 70 +/- 10 cm(-1) depending on the ligands L. This indicates that there are some vibrational modes with smaller force constants in the excited states than in the ground states for most of these complexes. These blue shifts increase approximately as the energy difference between the excited and ground states decreases, but they are otherwise not strongly correlated with the number of bipyridine ligands in the complex. Careful comparisons of the [Ru(L)(4)(d(8)-bpy)](2+) and [Ru(L)(4)(h(8)-bpy](2+) emission spectra are used to resolve the very weak vibronic contributions of the C-H stretching modes as the composite contributions of the corresponding vibrational reorganizational energies. The largest of these, 25 +/- 10 cm(-1), is found for the complexes with L = py or bpy/2 and smaller when L = NH(3). Perdeuteration of the am(m)ine ligands (NH(3), en, or [14]aneN(4)) has no significant effect on the zero-point energy difference, and the contributions of the NH stretching vibrational modes to the emission band shape are too weak to resolve. Ligand perdeuteration does increase the excited-state lifetimes by a factor that is roughly proportional to the excited-state-ground-state energy difference, even though the CH and NH vibrational reorganizational energies are too small for nuclear tunneling involving these modes to dominate the relaxation process. It is proposed that metal-ligand skeletal vibrational modes and configurational mixing between metal-centered, bpy-ligand-centered, and MLCT excited states are important in determining the zero-point energy differences, while a large number of different combinations of relatively low-frequency vibrational modes must contribute to the nonradiative relaxation of the MLCT excited states.     

On the contribution of vibrational anharmonicity to the binding energies of water clusters

The second-order vibrational perturbation theory method has been used together with the B3LYP and MP2 electronic structure methods to investigate the effects of anharmonicity on the vibrational zero-point energy (ZPE) contributions to the binding energies of (H2O)n, n = 2-6, clusters. For the low-lying isomers of (H2O)6, the anharmonicity correction to the binding energy is calculated to range from -248 to -355 cm(-1). It is also demonstrated that although high-order electron correlation effects are important for the individual vibrational frequencies, they are relatively unimportant for the net ZPE contributions to the binding energies of water clusters.     

低温基质隔离红外光谱和从头算研究CF~2O...IF

在10KAr和Xe基质中测得CF~2O...IF的两种分子间配合物A和B的红外光谱,观察到A和B之间的异构化反应。采用量子化学从头计算,在MP2/LANL2DZ水平上计算得到A和B的平衡构型的结构参数和振动频率。IF与CF~2O的相互作用引起CF~2O的C=O双键伸长、C-F键缩短,导致v~C~=~O的红移和v~C~F~2的蓝移。由势能曲线可看出,A的势能比B高23kJ/mol,表明B比A稳定,升温能够发生A向B的转化。     

Theoretical investigation of the weakly dihydrogen bonded complexes FArCCH...HBeX (X = H, F, Cl, Br)

An ab initio computational study of the properties of four linear dihydrogen-bonded complexes formed between the first compound with an Ar-C chemical bond (FArCCH) and HBeX (X = H, F, Cl, and Br) molecules was undertaken at the MP2/6-311++G(2d,2p) level of theory. The calculated complexation energy at MP2 and G2(MP2) levels decreases in the order HBeH...HCCArF > BrBeH...HCCArF > ClBeH...HCCArF > FBeH...HCCArF. The intermolecular stretching frequency, and shifts within the monomers, are compared with the energetic strength of complexation.     

Zero-point corrections and temperature dependence of HD spin-spin coupling constants of heavy metal hydride and dihydrogen complexes calculated by vibrational averaging

Vibrational corrections (zero-point and temperature dependent) of the H-D spin-spin coupling constant J(HD) for six transition metal hydride and dihydrogen complexes have been computed from a vibrational average of J(HD) as a function of temperature. Effective (vibrationally averaged) H-D distances have also been determined. The very strong temperature dependence of J(HD) for one of the complexes, [Ir(dmpm)Cp*H2]2 + (dmpm = bis(dimethylphosphino)methane) can be modeled simply by the Boltzmann average of the zero-point vibrationally averaged JHD of two isomers. For this complex and four others, the vibrational corrections to JHD are shown to be highly significant and lead to improved agreement between theory and experiment in most cases. The zero-point vibrational correction is important for all complexes. Depending on the shape of the potential energy and J-coupling surfaces, for some of the complexes higher vibrationally excited states can also contribute to the vibrational corrections at temperatures above 0 K and lead to a temperature dependence. We identify different classes of complexes where a significant temperature dependence of J(HD) may or may not occur for different reasons. A method is outlined by which the temperature dependence of the HD spin-spin coupling constant can be determined with standard quantum chemistry software. Comparisons are made with experimental data and previously calculated values where applicable. We also discuss an example where a low-order expansion around the minimum of a complicated potential energy surface appears not to be sufficient for reproducing the experimentally observed temperature dependence.     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

(N_2…CO)~+体系相互作用的密度泛函理论研究

在2种密度泛函方法和适宜基组水平上,对(N2…CO)+体系可能存在的相互作用复合物进行了全自由度能量梯度优化,发现势能面上存在2个能量极小点,均为共平面型。 比较了它们之间的相对稳定性,并对其进行了轨道成键分析,同时探讨了最稳定结构A的正则振动模式。 通过消除基函数引起的重叠误差(BSSE)和零点振动能(ZPVE)的校正,精确求算出复合物结构A、B的相互作用能DE分别为125.0和61.0 kJ/mol, 同等电子体(CO…CO)+相比,二者存在较大的差异。     

F原子与瞬态自由基CH_2SH反应的理论研究

用量子化学从头算和密度泛函理论（DFT）对F原子与自由基CH_2SH在势能面上 的反应进行了研究。在B3LYP/6-311G水平上计算出了各物种的优化构型、振动频率 和零点振动能(ZPVE);各物种的总能量由B3LYP/6-311 + G(2df, pd)//B3LYP/6- 311G计算,另外对反应物和产物还计算了其G3能量。结果表明：首先F通过与C或S 结合的两种途径与CH_2SH相配位,再通过H（4）原子转移形成甲基,然后甲基再旋 转,甲基中H（4）原子最终与F结合,反应产物为HF和CH_2S。反应为放热反应,分 别为ΔH_r = -370.7 kJ/mol (DFT)和-396.94 kJ/mol (G3)。此外依据计算出的反 应热,可得自由基·CH_2SH的生成热Δ_fH°_(298.15) = 146.44 kJ/mol (DFT), 而Δ_fH°_0 = 167.36 kJ/mol (G3)。它们与以前的实验和理论值是一致的。     

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.     

氧负离子自由基与苯的反应机理研究

采用密度泛函理论研究了氧负离子自由基与苯分子的反应. 针对可能存在的抽氢反应、置换氢反应、生成水反应及羟基化反应通道, 在B3LYP/6-31+G(d, p)的水平上得到了反应物、产物、中间体复合物和过渡态的优化构型、振动频率以及能量, 证实了这个反应的多产物通道性质, 较完整地描述了反应的机理和产物生成途径. 并且利用G2MP2能量计算结果, 从反应途径的能垒高度定性说明了实验观测结果的合理性, 同时也解释了以往实验的矛盾结果.     

Validation of the density-functional based tight-binding approximation method for the calculation of reaction energies and other data

We investigated the performance of the approximative density functional method DFTB versus BLYP and G2 with respect to zero-point corrected reaction energies, vibrational frequencies, and geometry parameters for a set of 28 reactions and 22 representative molecules containing C, H, N, and O (DFTB--density-functional based tight-binding approximation). The DFTB reaction energies show a mean absolute deviation versus the G2 reference of 4.3 kcalmol only. The corresponding value for the vibrational frequencies amounts to 75 cm(-1) versus BLYP/cc-pVTZ. With very few exceptions bond lengths and angles are in excellent agreement with the results of higher-level methods.     

The structure, properties, and nature of HArF-benzene complex: redshift and blueshift of Ar-H stretch frequency and rare gas atomic number dependence of hydrogen bonds

Ab initio calculations have been performed for the complexes of benzene with HArF, HKrF, and HXeF. The computed results indicate that the complexes of benzene-HArF exist in different conformations and among them those with π-hydrogen bonds are the more stable than those with C-H···F hydrogen bonds. Interestingly, the Ar-H stretching frequency is redshifted in the more stable isomer and blueshifted in the less stable form. The Ng (Ng=Ar, Kr, and Xe) atomic number dependence of the Ng-H···π and C-H···F hydrogen bonds has been explored. The result indicates that the strength of Ng-H···π and C-H···F hydrogen bonds is weakened with the increase of Ng atomic number. Natural bond orbital analysis has been performed to understand the interaction nature, frequency shift of H-Ng stretch, and dependence of Ng-H···π and C-H···F hydrogen bonds on the Ng atomic number.     

Noncovalent interactions between cytosine and SWCNT: curvature dependence of complexes via pi...pi stacking and cooperative CH...pi/NH...pi

Fragments of C24H12, adapted from a variety of armchair [(n,n), (n = 5, 7, and 8)] and zigzag [(m,0) (m = 8, 10, and 12)] single-walled carbon nanotube (SWCNT), are used to model corresponding SWCNTs with different diameters and electronic structures. The parallel binding mainly through pi...pi stacking interaction, as well as the perpendicular binding via cooperative NH...pi and CH...pi between cytosine and the fragments of SWCNT have been extensively investigated with a GGA type of DFT, PW91LYP/6-311++G(d,p). The eclipsed tangential (ET) conformation with respect to the six-membered ring of cytosine and the central ring of SWCNT fragments is less stable than the slipped tangential (ST) conformation for the given fragment; perpendicular conformations with NH2 and CH ends have higher negative binding energy than those with NH and CH ends. At PW91LYP/6-311++G(d,p) level, two tangential complexes are less bound than perpendicular complexes. However, as electron correlation is treated with MP2/6-311G(d,p) for PW91LYP/6-311++G(d,p) optimized complexes, it turns out there is an opposite trend that two tangential complexes become more stable than three perpendicular complexes. This result implies that electron correlation, a primary source to dispersion energy, has more significant contributions to the pi...pi stacking complexes than to the complexes via cooperative NH...pi and CH...pi interactions. In addition, it was found for the first time that binding energies for two tangential complexes become more negative with increasing nanotube diameter, while those for three perpendicular complexes have a weaker dependence on the curvature; i.e., binding energies are slightly less and less negative. The performance of a novel hybrid DFT, MPWB1K, was also discussed.