Inverse hydrogen bonds between SiH_4 and hydrides of Na,Mg and Be

The optimized geometries of the three complexes between MeHn (Me=Na,Mg,Be;n=1 or 2) and SiH4 have been calculated at the B3LYP/6-311++g**,MP2/6-311++g(3df,3pd) and MP2/aug-cc-pvtz levels,respectively.The red-shift inverse hydrogen bonds (IHBs) based on Si-H,an electron donor,were reported.The calculated binding energies with basis set super-position error (BSSE) correction of the three complexes are-5.98,-8.65 and-3.96 kJ mol-1 (MP2/6-311++g(3df,3pd)),respectively,which agree with the results obtained via M...     

PH~3...H~2O体系分子间相互作用的从头计算MO研究

在MP2/6-311++G(3d,3p)水平,对PH~3...H~2O体系可能存在的氢键复合物进行了全自由度能量梯度优化,发现PH~3...H~2O体系存在三个能量极小结构A,B和C。其中结构A和B以H~2O作为质子授体,结构C以PH~3作为质子授体,结构A较结构B和C分别稳定6.52kJ/mol和8.18kJ/mol。结构A具有C~s对称性,其结构中P原子和O原子间距离为354.78nm,键角OHP为171.35ⅲ,属于接近于直线的传统型氢键结构。进一步在高级电子相关校正的MP4SDTQ下,用6-311++G(3df,3pd)基加上键函数{3s3p2d1f},通过BSSE校正,精确计算了结构A的结合能为-10.84kJ/mol。     

Microwave spectrum of 3-butyne-1-thiol: evidence for intramolecular S-H...pi hydrogen bonding

The microwave spectrum of 3-butyne-1-thiol has been studied by means of Stark-modulation microwave spectroscopy and quantum-chemical calculations employing the B3LYP/6-311++G(3df,2pd), MP2/aug-cc-pVTZ, MP2/6-311++G(3df,2pd), and G3 methods. Rotational transitions attributable to two conformers of this molecule were assigned. One of these conformers possesses an antiperiplanar arrangement of the atoms S-C1-C2-C3, while the other is synclinal and stabilized by the formation of an intramolecular hydrogen bond between the H-atom of the thiol group and the pi-electrons of the C[triple bond]C triple bond. The energy difference between these conformers was estimated to be 1.7(4) kJ mol(-1) by relative intensity measurements, with the hydrogen-bonded conformer being lower in energy. The spectra of five vibrationally excited states of the synclinal conformer were observed, and an assignment of these states to particular vibrational modes was made with the aid of a density functional theory (DFT) calculation of the vibrational frequencies at the B3LYP/6-311++G(3df,2pd) level of theory.     

The molecular structure of fluoromalononitrile as determined by gas-phase electron diffraction and ab initio calculations

The molecular structure of fluoromalononitrile was studied by means of gas-phase electron diffraction and quantum mechanical methods using HF/6-31G(d), MP2/6-311++G(2df,2pd) and DFT/B3LYP/6-31G(d), B3PW91/6-31G(d), B3LYP/6-311++G(2df,2pd) and B3PW91/6-311++G(2df,2pd). The r(g) and angle(alpha) structural parameters we obtained from the present analysis are: CC=1.487(5) A, CN=1.157(3) A, CF=1.386(5) A, CH=1.096 A (ass.), angleCCC=106.7(1.0) degrees , angleCCF=108.0(0.7) degrees , angleCCN=177.6(2.0) degrees . Uncertainties in parenthesis are 3sigma.     

An ab initio correlated study of the potential energy surface for the HOBr.H2O complex

The potential energy surface (PES) for the HOBr.H(2)O complex has been investigated using second- and fourth-order M?ller-Plesset perturbation theory (MP2, MP4) and coupled cluster theory with single and doubles excitations (CCSD), and a perturbative approximation of triple excitations (CCSD-T), correlated ab initio levels of theory employing basis sets of triple zeta quality with polarization and diffuse functions up to the 6-311++G(3dp,3df ) standard Pople's basis set. Six stationary points being three minima, two first-order transition state (TS) structures and one second-order TS were located on the PES. The global minimum syn and the anti equilibrium structure are virtually degenerated [DeltaE(ele-nuc) approximately 0.3 kcal mol(-1), CCSD-T/6-311++G(3df,3pd) value], with the third minima being approximately 4 kcal mol(-1) away. IRC analysis was performed to confirm the correct connectivity of the two first-order TS structures. The CCSD-T/6-311++G(3df,3pd)//MP2/6-311G(d,p) barrier for the syn<-->anti interconversion is 0.3 kcal mol(-1), indicating that a mixture of the syn and anti forms of the HOBr.H(2)O complex is likely to exist.     

A microwave spectroscopic and quantum chemical study of 3-butyne-1-selenol (HSeCH2CH2C[triple bond]CH)

The microwave spectrum of 3-butyne-1-selenol has been studied by means of Stark-modulation microwave spectroscopy and quantum chemical calculations employing the B3LYP/aug-cc-pVTZ and MP2/6-311++G(3df,3pd) methods. Rotational transitions attributable to the H80SeCH2CH2C[triple bond]CH and H78SeCH2CH2C[triple bond]CH isotopologues of two conformers of this molecule were assigned. One of these conformers possesses an antiperiplanar arrangement for the atoms Se-C-C-C, while the other is synclinal and seems to be stabilized by the formation of a weak intramolecular hydrogen bond between the hydrogen atom of the selenol group and the pi electrons of the CC triple bond. The energy difference between these conformers was determined to be 0.2(5) kJ/mol by relative intensity measurements, and the hydrogen-bonded form was slightly lower in energy.     

Properties of the halogen-hydride interaction: an ab initio and "atoms in molecules" analysis

X-Cl...H-Y interactions are analyzed by applying ab initio methods as well as the Bader theory. All calculations were performed using Pople's basis sets (6-311++G(2df,2pd) and 6-311++G(3df,3pd)) as well as the Dunning-type bases (aug-cc-pVDZ and aug-cc-pVTZ) within the MP2 method. For the complexes analyzed here, X-Cl and H-Y may be treated as a Lewis acid and a Lewis base, respectively. The Cl...H interactions are rather weak or at most moderate since, for the strongest interaction of the F3...HLi complex, the binding energy calculated at the MP2/6-311++G(3df,3pd) level of approximation amounts to -3.4 kcal/mol, and the H...Cl distance is equal to 2.65 A, less than the corresponding sum of van der Waals radii. These interactions may be classified as halogen-hydride interactions. However, some of the complexes analyzed, especially F3SiCl...HBeF and F3SiCl...HBeF, are very weakly bound, probably by typical van der Waals interactions.     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

Theoretical Comparison of the Linear and Bent Structures for the Weakly Bound CO2—HF Complex

Density‐functional theory method (DFT) B3LYP/6–311++G(3df,2pd) and Moller‐Plesset perturbation method (MP2) MP2/6–311++G(3df,2pd) of Gaussian 03 were selected for the theoretical study of weakly bound CO₂—HF complex. In addition to the well‐known linear structure, the various bent structure complexes were also found in this work. The self‐consistent energy differences were only around 0.02 kJ/mol between the bent structure and linear structure by comparison. From the results of H‐bonding distance, dHF elongation and red shift of VHF vibration frequency, all the evidence shows that the H‐bonding effect in the bent structure is stronger than the linear structure. However, if one compares the Gibbs energy of the complex formation by temperature variation, it is very easily found that the linear form is favored under the thermal conditions of most temperatures whenever T ≥ 40 K. Such a fact is consistent with the former spectroscopic observed result of Klemperer et al.     

Computational investigation of the weakly bound dimers HOX...SO(3) (X = F, Cl, Br)

The electronic structure and thermochemical stability of the HOX-SO(3) (X = F, Cl, Br) complexes is studied using second-order M?ller-Plesset perturbation theory (MP2). The calculated dissociation energies of the HOF-SO(3), HOCl-SO(3), and HOBr-SO(3) complexes are 5.43, 6.02, and 5.98 kcal mol(-1) at MP2/6-311++G(3df,3pd) level, respectively. Anharmonic OH stretching frequencies of the HOX (X = F, Cl, Br) moieties along with the frequency shifts upon complex formation are calculated at the MP2/6-311++G(2df,2p) level. AIM and NBO analyses were also performed. Theoretical data strongly encourage performing of matrix-isolation studies of the title complexes and their spectroscopic identification.     

高级量子化学从头计算法研究N2和H2O分子间相互作用

在MP2/6-311++G(3d,3p)电子相关校正水平上,对N₂和H₂O分子间可能存在的氢键复合物进行全自由度能量梯度优化,发现了一个接近于直线的弱氢键总能量极小结构(1),进一步在高级电子相关校正的MP4SDTQ和CCSD(T)水平,用6-311++G(3d,3p)基组加上(3s3p2d1f)键函数,用MP4和CCSD(T)计算的结构1的结合能分别为-5.061kJ/mol和-4.715kJ/mol.     

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

Intramolecular homolytic substitution of sulfinates and sulfinamides--a computational study

Ab initio and density functional theory (DFT) calculations predict that intramolecular homolytic substitution by alkyl radicals at the sulfur atom in sulfinates proceeds through a smooth transition state in which the attacking and leaving radicals adopt a near collinear arrangement. When forming a five-membered ring and the leaving radical is methyl, G3(MP2)-RAD//ROBHandHLYP/6-311++G(d,p) calculations predict that this reaction proceeds with an activation energy (ΔE(1)(?)) of 43.2 kJ mol(-1). ROBHandHLYP/6-311++G(d,p) calculations suggest that the formation of five-membered rings through intramolecular homolytic substitution by aryl radicals at the sulfur atom in sulfinates and sulfinamides, with expulsion of phenyl radicals, proceeds with the involvement of hypervalent intermediates. These intermediates further dissociate to the observed products, with overall energy barriers of 45-68 kJ mol(-1), depending on the system of interest. In each case, homolytic addition to the phenyl group competes with substitution, with calculated barriers of 51-78 kJ mol(-1). This computational study complements and provides insight into previous experimental observations.     

FC(O)O自由基与NO2反应的量子化学研究

用量子化学DFT, MP2, G3和G3MP2方法对FC(O)O自由基与NO₂的反应机理进行了理论研究. 优化了反应势能面上各驻点的几何结构, 通过内禀反应坐标(IRC)计算和振动分析, 确认了反应中的过渡态, 并用过渡态理论(TST)计算了相关反应的速率常数.     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

A theoretical study of nonlinearity in heterocyclic hydrogen-bonded complexes

Ab initio calculations are performed at the MP2/6-311++G(d,p) and DFT/B3LYP/6-311++G(d,p) theoretical levels to obtain geometries, H-bond energies and harmonic infrared vibrational properties for the Cs symmetry structures of heterocyclic hydrogen-bonded complexes, CnHmY-HX. The H-bond lengths in DFT/B3LYP calculation level are in better agreement with the experimental values than the MP2 results. The geometry optimization are interpreted in terms of hydrogen bond nonlinearity represented by theta; and phi angles, once the hydrogen bond is formed among n-electrons pairs of the heteroatom in heterocyclic and the hydrogen atom in HX. The hydrogen bond energy after of the zero-point vibrational energy (ZPE) and basis set superposition error (BSSE) corrections are overestimated at DFT/B3LYP, whereas the MP2 BSSE corrections are very large than corresponding DFT/B3LYP. For example, the BSSE corrections for the C2H4S-HNC complex are 7.60 and 0.09 kJ mol(-1) in MP2 and DFT/B3LYP calculations levels, respectively. The new vibrational modes in infrared harmonic spectrum arising from complexation show several interesting features, especially the intermolecular stretching mode.     

Dihydrogen trioxide clusters, (HOOOH)n (n = 2-4), and the hydrogen-bonded complexes of HOOOH with acetone and dimethyl ether: implications for the decomposition of HOOOH

Hydrogen-bonded gas-phase molecular clusters of dihydrogen trioxide (HOOOH) have been investigated using DFT (B3LYP/6-311++G(3df,3pd)) and MP2/6-311++G(3df,3pd) methods. The binding energies, vibrational frequencies, and dipole moments for the various dimer, trimer, and tetramer structures, in which HOOOH acts as a proton donor as well as an acceptor, are reported. The stronger binding interaction in the HOOOH dimer, as compared to that in the analogous cyclic structure of the HOOH dimer, indicates that dihydrogen trioxide is a stronger acid than hydrogen peroxide. A new decomposition pathway for HOOOH was explored. Decomposition occurs via an eight-membered ring transition state for the intermolecular (slightly asynchronous) transfer of two protons between the HOOOH molecules, which form a cyclic dimer, to produce water and singlet oxygen (Delta (1)O 2). This autocatalytic decomposition appears to explain a relatively fast decomposition (Delta H a(298K) = 19.9 kcal/mol, B3LYP/6-311+G(d,p)) of HOOOH in nonpolar (inert) solvents, which might even compete with the water-assisted decomposition of this simplest of polyoxides (Delta H a(298K) = 18.8 kcal/mol for (H 2O) 2-assisted decomposition) in more polar solvents. The formation of relatively strongly hydrogen-bonded complexes between HOOOH and organic oxygen bases, HOOOH-B (B = acetone and dimethyl ether), strongly retards the decomposition in these bases as solvents, most likely by preventing such a proton transfer.     

Microwave spectrum, conformational composition, and intramolecular hydrogen bonding of (2-chloroethyl)amine (ClCH2CH2NH2)

The microwave spectrum of (2-chloroethyl)amine, ClCH(2)CH(2)NH(2), has been investigated in the 22-120 GHz region. Five rotameric forms are possible for this compound. In two of these conformers, denoted I and II, the Cl-C-C-N chain of atoms is antiperiplanar, with different orientations of the amino group. The link of the said atoms is synclinal in the three remaining forms, III-V, which differ with respect to the orientation of the amino group. The microwave spectra of four of these conformers, I-IV, have been assigned. In two of these rotamers, III and IV, the amino group is oriented in such a manner that rare and weak five-membered N-H···Cl intramolecular hydrogen bonds are formed. The geometries of conformers I and II preclude a stabilization by this interaction. The energy differences between the conformers were obtained from relative intensity measurements of spectral lines. The hydrogen-bonded conformer IV represents the global energy minimum. This rotamer is 0.3(7) kJ/mol more stable than the other hydrogen-bonded conformer III, 4.1(11) kJ/mol more stable than II, and 5.5(15) kJ/mol more stable than I. The spectroscopic work has been augmented by quantum chemical calculations at the CCSD/cc-pVTZ and MP2/6-311++G(3df,3pd) levels of theory. The CCSD rotational constants and energy differences are in good agreement with their experimental counterparts.     

Substituent effects on the acidity of weak acids. 3. Phenols and benzyl alcohols

The gas-phase acidities of meta- and para-substituted phenols have been calculated at the B3LYP/6-311+G, MP2/6-311+G, MP2/6-311++G, and MP2/6-311+G(2df,2pd) theoretical levels. The larger basis sets give the more satisfactory DeltaHacid values that are correlated with the observed acidities with a slope close to unity. They are systematically about 2 kcal/mol too small. The acidities of most substituted phenols are linearly related to those of the corresponding substituted benzoic acids and benzyl alcohols. Here, the substituent effect is a Coulombic interaction (field effect) of the distributed charge in the benzene ring with the negative charge of the anionic center. The exceptions are strong para-substituted pi-acceptors such as NO(2) and CHO, and to a smaller extent, CN and CF(3). Here there is direct charge transfer from the phenoxide oxygen to the substituent.     

Electronic properties of multifurcated bent hydrogen bonds CH3...Y and CH2...Y

H-bonding angle angleYHX has an important effect on the electronic properties of the H-bond Y...HX, such as intra- and intermolecular hyperconjugations and rehybridization, and topological properties of electron density. We studied the multifurcated bent H-bonds of the proton donors H3CZ (Z = F, Cl, Br), H2CO and H2CF2 with the proton acceptors Cl(-) and Br(-) at the four high levels of theory: MP2/6-311++G(d,p), MP2/6-311++G(2df,2p), MP2/6-311++G(3df,3pd) and QCISD/6-311++G(d,p), and found that they are all blue-shifted. These complexes have large interaction energies, 7-12 kcal mol(-1), and large blue shifts, delta r(HC) = -0.0025 --0.006 A and delta v(HC) = 30-90 cm(-1). The natural bond orbital analysis shows that the blue shifts of these H-bonds Y...HnCZ are mainly caused by three factors: rehybridization; indirect intermolecular hyperconjugation n(Y) -->sigma*(CZ), in that the electron density from n(Y) of the proton acceptor is transferred not to sigma*(CH), but to sigma*(CZ) of the donor; intramolecular hyperconjugation n(Z) -->sigma*(CH), in that the electron density in sigma*(CH) comes back to n(Z) of the donor such that the occupancy in sigma*(CH) decreases. The topological properties of the electron density of the bifurcated H-bonds Y...H2CZ are similar to those of the usual linear H-bonds, there is a bond critical point between Y and each hydrogen, and a ring critical point inside the tetragon YHCH. However, the topological properties of electron density of the trifurcated H-bonds Y...H3CZ are essentially different from those of linear H-bonds, in that the intermolecular bond critical point, which represents a closed-shell interaction, is not between Y and hydrogen, but between Y and carbon.