Interaction of anions with perfluoro aromatic compounds

The complexes formed by a variety of anions with perfluoro derivatives of benzene, naphthalene, pyridine, thiophene, and furan have been calculated using DFT (B3LYP/6-31++G**) and MP2 (MP2/6-31++G** and MP2/6-311++G**) ab initio methods. The minimum structures show the anion interacting with the pi-cloud of the aromatic compounds. The interaction energies obtained range between -8 and -19 kcal mol(-1). The results obtained at the MP2/6-31++G** and MP2/6-311++G** levels are similar. However, the B3LYP/6-31++G** results provide longer interaction distances and smaller interaction energies than do the MP2 results. The interaction energies have been partitioned using an electrostatic, polarization, and van der Waals scheme. The AIM analysis of the electron density shows a variety of topologies depending on the aromatic system considered.     

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.     

Pi-pi interaction in pyridine

pi-pi Interaction in pyridine dimer and trimer has been investigated in different geometries and orientations at the ab initio (HF, MP2) and DFT (B3LYP) levels of theory using various basis sets (6-31G, 6-31G, 6-311++G) and corrected for basis set superposition error (BSSE). While the HF and DFT calculations show the pyridine dimer and the trimer to be unstable with respect to the monomer, the MP2 calculations show them to be clearly stable, thus emphasizing the need to include electron correlation while determining stacking interaction in such systems. The calculated MP2/6-311++G binding energy (100% BSSE corrected) of the parallel-sandwich, antiparallel-sandwich, parallel-displaced, antiparallel-displaced, T-up and T-down geometries for pyridine dimer are 1.53, 3.05, 2.39, 3.97, 1.91, 1.47 kcal/mol, respectively. The results show the antiparallel-displaced geometry to be the most stable. The binding energies for the trimer in parallel-sandwich, antiparallel-sandwich, and antiparallel-displaced geometry are found to be 3.18, 6.14, and 8.04 kcal/mol, respectively.     

The C7-C10 cycloalkanes revisited

The conformers of cycloheptane through cyclodecane have been examined at the B3LYP/6-311+G* and MP2/6-311+G* theoretical levels, with some additional calculations at the CCD/6-311+G* and CCSD(T)/6-311++G** levels. With cyclooctane, B3LYP predicts that the boat-chair and crown conformers have similar energies, whereas MP2 and CCSD(T) predict that the crown conformer is 2 kcal/mol higher in energy. The latter is in agreement with the electron diffraction data. With cyclononane, B3LYP predicts that two of the higher-energy conformers found in molecular mechanics calculations should convert to one of the lower-energy conformers. However, MP2/6-311+G* optimizations find them to be true minima on the potential energy surface. B3LYP systematically predicts larger C-C-C bond angles for these compounds than either MP2 or CCD. The results of molecular mechanics MM4 calculations are generally in good agreement with those obtained using MP2.     

Conformational energies for 2-substituted butanes

The conformational free energies for some 2-substituted butanes where X = F, Cl, CN, and CCH were calculated using G3-B3, CBS-QB3, and CCSD(T)/6-311++G(2d,p) as well as other theoretical levels. The above methods gave consistent results with free energies relative to the trans conformers as follows: X = CCH, g+ = 0.77 +/- 0.05 kcal/mol. g- = 0.88 +/- 0.05 kcal/mol; X = CN, g+ = 0.85 +/- 0.05 kcal/mol, g- = 0.75 +/- 0.05 kcal/mol; X = Cl, g+ = 0.70 +/- 0.05 kcal/ml, g- = 0.80 +/- 0.05 kcal/mol; and X = F, g+ = 0.53 +/- 0.05 kcal/mol, g- = 0.83 +/- 0.05 kcal/mol. The conformational free energies also were estimated using the observed liquid phase IR spectra and intensities calculated using B3LYP/6-311++G** and MP2/6-311++G**. The rotational free energy profiles for all of the compounds were estimated at the G3-B3 level.     

Interaction of adenine adducts with thymine: a computational study

The existence of DNA adducts bring the danger of carcinogenesis because of mispairing with normal DNA bases. 1,N6-ethenoadenine adducts (epsilonA) and 1,N6-ethanoadenine adducts (EA) have been considered as DNA adducts to study the interaction with thymine, as DNA base. Several different stable conformers for each type of adenine adduct with thymine, [epsilonA(1)-T(I), epsilonA(2)-T(I), epsilonA(3)-T(I) and EA(1)-T(I), EA(2)-T(I), EA(3)-T(I)] and [epsilonA(1)-T(II), epsilonA(2)-T(II), epsilonA(3)-T(II) and EA(1)-T(II), EA(2)-T(II), EA(3)-T(II)], have been considered with regard to their interactions. The differences in their geometrical structures, energetic properties, and hydrogen-bonding strengths have also been compared with Watson-Crick adenine-thymine base pair (A-T). Single-point energy calculations at MP2/6-311++G** levels on B3LYP/6-31+G* optimized geometries have also been carried out to better estimate the hydrogen-bonding strengths. The basis set superposition error corrected hydrogen-bonding strength sequence at MP2/6-311++G**//B3LYP/6-31+G* for the most stable complexes is found to be EA(2)-T(I) (15.30 kcal/mol) > EA(1)-T(II) (14.98 kcal/mol) > EA(3)-T(II) (14.68 kcal/mol) > epsilonA(2)-T(I) (14.54 kcal/mol) > epsilonA(3)-T(II) (14.22 kcal/mol) > epsilonA(3)-T(II) (13.64 kcal/mol) > A-T (13.62 kcal/mol). The calculated reaction enthalpy value for epsilonA(2)-T(I) is 10.05 kcal/mol, which is the highest among the etheno adduct-thymine complexes and about 1.55 kcal/mol more than those obtained for Watson-Crick A-T base pair and the reaction enthalpy value for EA(1)- T(II) is 10.22 kcal/mol, which is highest among the ethano addcut-thymine complexes and about 1.72 kcal/mol more than those obtained for Watson-Crick A-T base pair. The aim of this research is to provide fundamental understanding of adenine adduct and thymine interaction at the molecular level and to aid in future experimental studies toward finding the possible cause of DNA damage.     

Ab initio studies of the conformations of ethynyl formate, cyano formate and related ethynyl and cyano compounds

The optimized geometries, relative free energies and related thermodynamic properties, harmonic frequencies, and dipole moments have been calculated at the HF and MP2 levels for ethynyl formate (1a), ethynyl acetate (1b), cyano formate, HCO₂CN (1c), cyano acetate (1d), S-ethynyl thioformate (2a), S-ethynyl thioacetate (2b), S-cyano thioformate (2c), S-cyano thioacetate (2d), N-ethynylformamide (3a), N-ethynylacetamide (3b), N-cyanoformamide (3c), and N-cyanoacetamide (3d) with the gaussian 98 program. For ethynyl formate, the calculation for 25 °C at the MP2/6-311++G(df,pd) level predicts that the Z isomer is more stable by 1.23 kcal/mol. For S-ethynyl thioformate, calculations at the MP2/6-311++G(2d,2p) level predict that the E isomer is favored by 0.71 kcal/mol at 25 °C. The E isomers of N-ethynylformamide and N-ethynylacetamide were found at all levels to be more stable than the Z isomers at 25 °C. For cyano formate and cyano acetate, calculations at the MP2/6-311++G(df,pd) level predict that the Z isomers are more stable at 25 °C by 1.50 and 2.72 kcal/mol, respectively. At this level and temperature, the Z isomers of 2c, 2d, 3c, and 3d are predicted to have free energies of 0.46, −0.07, 1.22, and 2.28 kcal/mol, respectively, relative to the E conformations. Z to E free-energy barriers at 25 °C of 8.63, 10.64, 17.63, 7.39, and 14.03 kcal/mol were calculated for 1a, 2a, 3a, 1c, and 3c at the HF/6-311G(d,p) level, and at the HF/6-311+G(d,p) level, the free-energy barrier for 2c was 7.08 kcal/mol.     

Cooperative enhancement of water binding to crownophane by multiple hydrogen bonds: analysis by high level ab initio calculations.

The intermolecular interaction energy of the model system of the water-crownophane complex was analyzed. The water molecule has four hydrogen bonds, with the two hydrogen-donating phenolic hydroxy groups and two hydrogen-accepting oxygen atoms of the poly-oxyethylene chain of the crownophane in the complex. The MP2/6-311G(2d,2p) level calculations of the model system of the complex (hydrogen donating unit + hydrogen accepting unit + water) indicate that the binding energy of the water is 21.85 kcal/mol and that the hydrogen bond cooperativity increases the binding energy as much as 3.67 kcal/mol. The calculated interaction energies depend on the basis set, while the basis set dependence of the cooperative increment is negligible. Most of the cooperative increment is covered by the HF level calculation, which suggests that the major source of the hydrogen bond cooperativity in this system has its origin in induction. The BLYP/6-311G** and PW91/6-311G** level interaction energies of the model system are close to the MP2/6-311G** interaction energies, which suggests that the DFT calculations with these functionals are useful methods to evaluated the interactions of hydrogen bonded systems.     

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.     

Structure of the boronic acid dimer and the relative stabilities of its conformers

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order M?ller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.     

Theoretical conformational analysis for neurotransmitters in the gas phase and in aqueous solution. Norepinephrine

The natural neurotransmitter (R)-norepinephrine takes the monocationic form in 93% abundance at the physiological tissue pH of 7.4. Ab initio and DFT/B3LYP calculations were performed for 12 protonated conformers of (R)-norepinephrine in the gas phase with geometry optimizations up to the MP2/6-311++G level, and with single-point calculations up to the QCISD(T) level at the HF/6-31G-optimized geometries. Four monohydrates were studied at the MP2/6-31G//HF/6-31G level. In the gas phase, the G1 conformer is the most stable with phenyl.NH(3)(+) gauche and HO(alc).NH(3)(+) gauche arrangements. A strained intramolecular hydrogen bond was found for conformers (G1 and T) with close NH(3)(+) and OH groups. Upon rotation of the NH(3)(+) group as a whole unit about the C(beta)-C(alpha) axis, a 3-fold potential was calculated with free energies for barriers of 3-12 kcal/mol at the HF/6-31G level. Only small deviations were found in MP2/6-311++G single-point calculations. A 2-fold potential was calculated for the phenyl rotation with free energies of 11-13 kcal/mol for the barriers at T = 310 K and p = 1 atm. A molecular mechanics docking study of (R)-norepinephrine in a model binding pocket of the beta-adrenergic receptor shows that the ligand takes a conformation close to the T(3) arrangement. The effect of aqueous solvation was considered by the free energy perturbation method implemented in Monte Carlo simulations. There are 4-5 strongly bound water molecules in hydrogen bonds to the conformers. Although hydration stabilizes mostly the G2 form with gauche phenyl.NH(3)(+) arrangement and a water-exposed NH(3)(+) group, the conformer population becomes T > G1 > G2, in agreement with the PMR spectroscopy measurements by Solmajer et al. (Z. Naturforsch. 1983, 38c, 758). Solvent effects reduce the free energies for barriers to 3-6 and 9-12 kcal/mol for rotations about the C(beta)-C(alpha) and the C(1)(ring)-C(beta) axes, respectively.     

Multipole electrostatics in hydration free energy calculations

Hydration free energy (HFE) is generally used for evaluating molecular solubility, which is an important property for pharmaceutical and chemical engineering processes. Accurately predicting HFE is also recognized as one fundamental capability of molecular mechanics force field. Here, we present a systematic investigation on HFE calculations with AMOEBA polarizable force field at various parameterization and simulation conditions. The HFEs of seven small organic molecules have been obtained alchemically using the Bennett Acceptance Ratio method. We have compared two approaches to derive the atomic multipoles from quantum mechanical calculations: one directly from the new distributed multipole analysis and the other involving fitting to the electrostatic potential around the molecules. Wave functions solved at the MP2 level with four basis sets (6-311G*, 6-311++G(2d,2p), cc-pVTZ, and aug-cc-pVTZ) are used to derive the atomic multipoles. HFEs from all four basis sets show a reasonable agreement with experimental data (root mean square error 0.63 kcal/mol for aug-cc-pVTZ). We conclude that aug-cc-pVTZ gives the best performance when used with AMOEBA, and 6-311++G(2d,2p) is comparable but more efficient for larger systems. The results suggest that the inclusion of diffuse basis functions is important for capturing intermolecular interactions. The effect of long-range correction to van der Waals interaction on the hydration free energies is about 0.1 kcal/mol when the cutoff is 12?, and increases linearly with the number of atoms in the solute/ligand. In addition, we also discussed the results from a hybrid approach that combines polarizable solute with fixed-charge water in the HFE calculation.     

Estimation on the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides

Computation of accurate intramolecular hydrogen-bonding energies for peptides is of great importance in understanding the conformational stabilities of peptides and developing a more accurate force field for proteins. We have proposed a method to determine the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides. In this article, the method is further applied to evaluate the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in peptides. The optimal structures of the intramolecular 10-membered ring N-H...O=C hydrogen bonds in glycine and alanine tripetide molecules are obtained at the MP2 level with 6-31G(d), 6-31G(d,p), and 6-31+G(d,p) basis sets. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are then evaluated based on our method at the MP2/6-311++G(3df,2p) level with basis set superposition error correction. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are calculated to be in the range of -6.84 to -7.66, -4.44 to -4.98, and -6.95 to -7.88 kcal/mol. The method is also applied to estimate the individual intermolecular hydrogen-bonding energies in the dimers of amino-acetaldehyde, 2-amino-acetamide, formamide, and oxalamide, each dimer having two identical intermolecular hydrogen bonds. According to our method, the individual intermolecular hydrogen-bonding energies in the four dimers are calculated to be -1.77, -1.67, -6.35, and -4.82 kcal/mol at the MP2/6-311++G(d,p) level, which are in good agreement with the values of -1.84, -1.72, -6.23, and -4.93 kcal/mol predicted by the supermolecular method.     

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.     

HArF和HNO二聚体相互作用氢键的理论研究

The hydrogen bonding interaction of 1:1 dimer formed between HNO and HArF molecule has been completely investigated in the present study using Second-order M?ller-Plesset Perturbation (MP2) method in conjunction with 6-311+G**, 6-311++G** and 6-311++G(2d,2p) basis sets. The standard and CP-corrected calculations have been employed to determine the equilibrium structures, the vibrational frequencies and interaction energies. The interaction energies of the dimers were also calculated at G2MP2 level. Two stable structures are found as the minima. Dimer I(H···F)is a five-membered cyclic hydrogen bonded structure and is more stable than the Dimer II(H···O). The blue-shifted N-H···F hydrogen bond is confirmed with standard and CP-corrected calculations by the MP2 and DFT methods in conjunction with different basis sets. The results obtained at MP2 in conjunction with different basis sets show there is a red-shifted hydrogen bond (Ar-H···O) in the Dimer II(H···O). The topological and electronic properties, the origin of red- and blue-shifted hydrogen bonds were investigated at MP2/6-311++G(2d,2p) with CP corrected calculations. From the NBO analysis, the reasonable explanations for the red- and blue-shifted hydrogen bonds were proposed.     

Potential energy scans and vibrational assignments of cyclopropanecarboxylic acid and cyclopropanecarboxamide

The structural stability and internal rotations in cyclopropanecarboxylic acid and cyclopropanecarboxamide were investigated by the DFT-B3LYP and the ab initio MP2 calculations using 6-311G** and 6-311+G** basis sets. The computations were extended to the MP4//MP2/6-311G** and CCSD(T)//MP2/6-311G** single-point calculations. From the calculations the molecules were predicted to exist predominantly in the cis (C=O group eclipses the cyclopropane ring) with a cis-trans barrier of about 4-6kcal/mol. The OCOH torsional barrier in the acid was estimated to be about 12-13kcal/mol while the corresponding OCNH torsional barrier in the amide was calculated to be about 20kcal/mol. The equilibrium constant k for the cis<-->trans interconversion in cyclopropanecarboxylic acid was calculated to be 0.1729 at 298.15K that corresponds to an equilibrium mixture of about 85% cis and 15% trans. The vibrational frequencies were computed at the DFT-B3LYP level. Normal coordinate calculations were carried out and potential energy distributions were calculated for the low energy cis conformer of the molecules. Complete vibrational assignments were made on the basis of normal coordinate calculations and comparison with experimental data of the molecules.     

Conformational studies in the cyclohexane series. 2. Phenylcyclohexane and 1-methyl-1-phenylcyclohexane

The structures and relative energies of the conformers of phenylcyclohexane, and 1-methyl-1-phenylcyclohexane have been calculated at theoretical levels including HF/6-31G, B3LYP/6-311G, MP2/6-311G, MP2/6-311(2df,p), QCISD/6-311G, and QCISD/6-311G(2df,p). The latter gives conformational enthalpy (DeltaH degrees ), entropy (DeltaS degrees ), and free energy (DeltaG degrees ) values for phenylcyclohexane that are in excellent agreement with the experimental data. The calculations for 1-methyl-1-phenylcyclohexane find a free energy difference of 1.0 kcal/mol at -100 degrees C, favoring the conformation having an axial phenyl group, that is in only modest agreement with the experimental value of 0.32 +/- 0.04 kcal/mol. The origin of the phenyl rotational profiles for the conformers of phenylcyclohexane and 1-methyl-1-phenylcyclohexane is discussed.     

Rapid evaluation of the binding energies in hydrogen-bonded amide-thymine and amide-uracil dimers in gas phase

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 48 hydrogen‐bonded amide–thymine and amide–uracil dimers are evaluated from the analytic potential energy function established in our lab recently. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the BSSE correction. For all the 48 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 46 dimers, only two differences are larger than 0.50 kcal/mol and the largest one is only 0.60 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.050 Å for all the 48 dimers. The analytic potential energy function is further applied to four more complicated hydrogen‐bonded amide–base systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded amide–base dimers quickly and accurately. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Computational study of hydrogen-bonded complexes of HOCO with acids: HOCO⋯HCOOH, HOCO⋯H(2)SO(4), and HOCO⋯H(2)CO(3)

Quantum chemistry calculations at the density functional theory (DFT) (B3LYP), MP2, QCISD, QCISD(T), and CCSD(T) levels in conjunction with 6-311++G(2d,2p) and 6-311++G(2df,2p) basis sets have been performed to explore the binding energies of open-shell hydrogen bonded complexes formed between the HOCO radical (both cis-HOCO and trans-HOCO) and trans-HCOOH (formic acid), H(2)SO(4) (sulfuric acid), and cis-cis-H(2)CO(3) (carbonic acid). Calculations at the CCSD(T)∕6-311++G(2df,2p) level predict that these open-shell complexes have relatively large binding energies ranging between 9.4 to 13.5 kcal∕mol and that cis-HOCO (cH) binds more strongly compared to trans-HOCO in these complexes. The zero-point-energy-corrected binding strengths of the cH?Acid complexes are comparable to that of the formic acid homodimer complex (～13-14 kcal∕mol). Infrared fundamental frequencies and intensities of the complexes are computed within the harmonic approximation. Infrared spectroscopy is suggested as a potential useful tool for detection of these HOCO?Acid complexes in the laboratory as well as in various planetary atmospheres since complex formation is found to induce large frequency shifts and intensity enhancement of the H-bonded OH stretching fundamental relative to that of the corresponding parent monomers. Finally, the ability of an acid molecule such as formic acid to catalyze the inter-conversion between the cis- and trans-HOCO isomers in the gas phase is also discussed.     

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.