Preferred conformers and photochemical (lambda > 200 nm) reactivity of serine and 3,3-dideutero-serine in the neutral form

A systematic investigation of the conformational potential energy surface of neutral serine [HOCH2CHNH2COOH] and 3,3-dideutero-serine [HOCD2CHNH2COOH] was undertaken, revealing the existence of 61 different minima. The structures and vibrational spectra of the most stable conformers, which were estimated to have relative energies within 7 kJ mol(-1) and account for ca. 93% of the total conformational population at room temperature, were calculated at both the MP2 and DFT/BLYP levels of theory with the 6-311++G(d,p) basis-set and used to interpret the spectroscopic data obtained for the compounds isolated in low-temperature inert matrixes. The assignment of the main spectral infrared features observed in the range 4000-400 cm(-1) to the most stable conformers of serine was undertaken. In addition, UV irradiation (lambda > 200 nm) of the matrix-isolated compounds was also performed, leading to decarboxylation, which was found to be strongly dependent on the conformation assumed by the reactant molecule.     

Conformational study of monomeric 2,3-butanediols by matrix-isolation infrared spectroscopy and DFT calculations

The FT-IR spectra of two diastereomers of 2,3-butanediol, (R,S) and (S,S), isolated in low-temperature argon and xenon matrixes were studied, allowing the identification of two different conformers for each compound. These conformers were characterized by a +/-gauche arrangement around the O-C-C-O dihedral angle, thus enabling the establishment of a very weak intramolecular hydrogen bond of the O...H-O type. No other forms of these compounds were identified in matrixes, despite the fact that these four conformers had calculated relative energies from 0 to 5.1 kJ mol(-1) and were expected to be thermally populated from 50 to 6% in the gaseous phase of each compound. The nonobservation of additional conformers was explained in terms of low barriers to intramolecular rotation, resulting in the conformational relaxation of the compounds during deposition of the matrixes. The barriers to internal rotation of the OH groups were computed to be less than 4 kJ mol(-1) and are easily overcome in matrixes within the family of conformers with the same heavy atom backbone. The barriers for intramolecular rearrangement of the O-C-C-O dihedral angle in both diastereomers were calculated to range from 20 to 30 kJ mol(-1). Interconversions between the latter conformers were not observed in matrixes, even after annealing up to 65 K. Energy calculations, barriers, and calculated infrared spectra were carried out at the DFT(B3LYP)/6-311++G theory. Additional MP2/6-311++G calculations of energies and vibrational frequencies were performed on the most relevant conformers. Finally, independent estimations of the hydrogen-bond enthalpy in the studied molecules were also obtained based on theoretical structural data and from vibrational frequencies (using well-established empirical correlations). The obtained values for -DeltaH for both diastereomers of 2,3-butanediol amount to ca. 6-8 kJ mol(-1).     

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.     

Conformational stability from variable temperature infrared spectra of krypton solutions, ab initio calculations, vibrational assignment, and r0 structural parameters of 1,3-difluoropropane

The infrared spectra (3200-400 cm(-1)) of krypton solutions of 1,3-difluoropropane, FCH2CH2CH2F, at variable temperatures (-105 to -150 degrees C) have been recorded. Additionally, the infrared spectra (3200-50 cm(-1)) of the gas and solid have been recorded as well as the Raman spectrum of the liquid. From a comparison of the spectra of the fluid phases with that in the solid, all of the fundamental vibrations of the C2 conformer (gauche-gauche) where the first gauche indicates the form for one of the CH2F groups and the second gauche the other CH2F, and many of those for the C1 form (trans-gauche) have been identified. Tentative assignments have been made for a few of the fundamentals of the other two conformers, i.e. C2v (trans-trans) and Cs (gauche-gauche'). By utilizing six pairs of fundamentals for these two conformers in the krypton solutions, an enthalpy difference of 277 +/- 28 cm(-1) (3.31 +/- 0.33 kJ mol(-1)) has been obtained for the C2 versus C1 conformer with the C2 conformer the more stable form. For the C2v conformer, the enthalpy difference has been determined to be 716 +/- 72 cm(-1) (8.57 +/- 0.86 kJ mol(-1)) and for the Cs form 971 +/- 115 cm(-1) (11.6 +/- 1.4 kJ mol(-1)). It is estimated that there is 64 +/- 3% of the C2 form, 34 +/-3% of the C1 form, 1% of the C2v form and 0.6% of the Cs conformer present at ambient temperature. Equilibrium geometries and total energies of the four stable conformers have been determined from ab initio calculations with full electron correlation by the perturbation method to second order as well as by hybrid density functional theory calculations with the B3LYP method using a number of basis sets. The MP2 calculations predict the C1 conformer stability to be slightly higher than the experimentally determined value whereas for the C2v and Cs conformers the predicted energy difference is much larger than the experimental value. The B3LYP calculations predict a better energy difference for both the C1 and C2v as well as for the Cs conformers than the MP2 values. A complete vibrational assignment is proposed for the C2 conformer and many of the fundamentals have been identified for the C1 form based on the force constants, relative intensities and rotational-vibrational band contours obtained from the predicted equilibrium geometry parameters. By combining previously reported rotational constants for the C2 and C1 conformers with ab initio MP2/6-311 + G(d, p) predicted parameters, adjusted r0 parameters have been obtained for both conformers. Comparisons are made with the parameters obtained for some other molecules containing the FCH2 group. The spectroscopic and theoretical results are compared to the corresponding properties for some similar molecules.     

Dissociation of carbonic acid: gas phase energetics and mechanism from ab initio metadynamics simulations

A comprehensive metadynamics study of the energetics, stability, conformational changes, and mechanism of dissociation of gas phase carbonic acid, H2CO3, yields significant new insight into these reactions. The equilibrium geometries, vibrational frequencies, and conformer energies calculated using the density functional theory are in good agreement with the previous theoretical predictions. At 315 K, the cis-cis conformer has a very short life time and transforms easily to the cis-trans conformer through a change in the O=C-O-H dihedral angle. The energy difference between the trans-trans and cis-trans conformers is very small (approximately 1 kcal/mol), but the trans-trans conformer is resistant to dissociation to carbon dioxide and water. The cis-trans conformer has a relatively short path for one of its hydroxyl groups to accept the proton from the other end of the molecule, resulting in a lower activation barrier for dissociation. Comparison of the free and potential energies of dissociation shows that the entropic contribution to the dissociation energy is less than 10%. The potential energy barrier for dissociation of H2CO3 to CO2 and H2O from the metadynamics calculations is 5-6 kcal/mol lower than in previous 0 K studies, possibly due to a combination of a finite temperature and more efficient sampling of the energy landscape in the metadynamics calculations. Gas phase carbonic acid dissociation is triggered by the dehydroxylation of one of the hydroxyl groups, which reorients as it approaches the proton on the other end of the molecule, thus facilitating a favorable H-O-H angle for the formation of a product H2O molecule. The major atomic reorganization of the other part of the molecule is a gradual straightening of the O=C=O bond. The metadynamics results provide a basis for future simulation of the more challenging carbonic acid-water system.     

Molecular structures of benzoic acid and 2-hydroxybenzoic acid, obtained by gas-phase electron diffraction and theoretical calculations

The structures of benzoic acid (C6H5COOH) and 2-hydroxybenzoic acid (C6H4OHCOOH) have been determined in the gas phase by electron diffraction using results from quantum chemical calculations to inform restraints used on the structural parameters. Theoretical methods (HF and MP2/6-311+G(d,p)) predict two conformers for benzoic acid, one which is 25.0 kJ mol(-1) (MP2) lower in energy than the other. In the low-energy form, the carboxyl group is coplanar with the phenyl ring and the O-H group eclipses the C=O bond. Theoretical calculations (HF and MP2/6-311+G(d,p)) carried out for 2-hydroxybenzoic acid gave evidence for seven stable conformers but one low-energy form (11.7 kJ mol(-1) lower in energy (MP2)) which again has the carboxyl group coplanar with the phenyl ring, the O-H of the carboxyl group eclipsing the C=O bond and the C=O of the carboxyl group oriented toward the O-H group of the phenyl ring. The effects of internal hydrogen bonding in 2-hydroxybenzoic acid can be clearly observed by comparison of pertinent structural parameters between the two compounds. These differences for 2-hydroxybenzoic acid include a shorter exocyclic C-C bond, a lengthening of the ring C-C bond between the substituents, and a shortening of the carboxylic single C-O bond.     

Isoelectronic homologues and isomers: tropolone, 5-azatropolone, 1-H-azepine-4,5-dione, saddle points, and ions

Computational studies of 12 64-electron homologues and isomers of tropolone in the S(0) electronic ground state are reported. Three minimum-energy structures, tropolone (Tp), 5-azatropolone (5Azt), and 5-H-5-azatropolonium (5AztH(+)), have an internal H-bond and planar C(s)) geometry, and three, tropolonate (TpO(-)), 5-azatropolonate (5AzO(-)), and 1-H-azepine-4,5-dione (45Di), lack the H-bond and have twisted C(2) geometry. All 6 substances have an equal double-minimum potential energy surface and a saddle point with planar C(2)(v) geometry. The energy for the gas-phase isomerization reaction 45Di --> 5Azt is near +4 kJ mol(-1) at the MP4(SDQ)/6-311++G(df,pd)//MP2/6-311++G(df,pd) (energy//geometry) theoretical level and around -20 kJ mol(-1) at lower theoretical levels. The dipole moments computed for 45Di and 5Azt are 9.6 and 2.1 D, respectively, and this large difference contributes to MO-computed free energies of solvation that strongly favor--as experimentally observed--45Di over 5Azt in chloroform solvent. The MO-computed energy for the gas-phase protonation reaction 45Di + H(+) --> 5AztH(+) is -956.4 kJ mol(-1), leading to 926.8 kJ mol(-1) as the estimated proton affinity for 45Di at 298 K and 1 atm. The intramolecular dynamical properties predicted for 5Azt and 5AztH(+) parallel those observed for tropolone. They are therefore expected to exhibit spectral tunneling doublets. Once they are synthesized, they should contribute importantly to the understanding of multidimensional intramolecular H transfer and dynamical coupling processes.     

Infrared multiple photon dissociation spectroscopy of potassiated proline

The structure of proline in [proline + K]+ has been investigated in the gas phase using high level DFT and MP2 calculations and infrared photo dissociation spectroscopy with a free electron laser (FELIX). The respective FELIX spectrum of [proline + K]+ matches convincingly the calculated spectra of two structurally closely related and nearly iso-energetic zwitterionic salt bridge (SB) structures. An additional unresolved band at approximately 1725 cm(-1) matching with the characteristic CO stretching mode of charge solvation (CS) structures points toward the presence of a minor population of these conformers of proline in [proline + K]+. However, theory predicts a significant energy gap of 18.9 kJ mol(-1) (B3LYP/6-311++G(2d,2p)) or 15.6 kJ mol(-1) (MP2) between the lowest CS conformer of proline and the clearly favored SB structure.     

Matrix isolation FTIR spectroscopic and theoretical study of dimethyl sulfite

The preferred conformations of dimethyl sulfite and their vibrational spectra were studied by matrix-isolation Fourier transform infrared spectroscopy and theoretical methods (density functional theory (DFT) and Moller-Plesset (MP2), with basis sets of different sizes, including the quadruple-zeta, aug-cc-pVQZ basis). Five minima were found at these levels of theory. At the MP2/6-31++G(d,p) and DFT/B3LYP/aug-cc-pVQPZ levels, the GG conformer (where the O-S-O-C dihedral angles are 73.2 and 70.8 degrees ) resulted in the conformational ground state. At the highest level of theory used, the GT conformer (O-S-O-C = +68.5 and -173.2 degrees ) is 0.83 kJ mol(-1) higher in energy than the GG form, while conformer GG' (O-S-O-C = +85.7 and -85.7 degrees ) has a relative energy of 1.18 kJ mol(-1). The remaining two conformers (G'T and TT) are high-energy forms and not experimentally relevant. In consonance with the theoretical predictions, conformer GG was found to be the most stable conformer in the gaseous phase as well as in the low-temperature matrices. Annealing of the argon matrices first promotes the GG'-->GT isomerization, which is followed by conversion of GT into the most stable conformer. There is no evidence of occurrence of GG'-->GG direct conversion in the low-temperature matrices. On the other hand, during deposition of the xenon matrices conformer GG' totally converts to conformer GT. Two observations demonstrated this fact: no evidence of bands corresponding to GG' were observed in xenon matrices and the GG/GT intensity ratio became similar to the GG/(GT + GG') intensity ratio observed in argon matrices. All these results could be explained by taking into account the relative values of the theoretically predicted energy barriers for the different isomerization processes: GG'-->GT, 1.90 kJ mol(-1); GT-->GG, 9.64 kJ mol(-1); and GG'-->GG, 19.46 kJ mol(-1).     

Unusual photochemical C-N bond cleavage in the novel methyl 2-chloro-3-methyl-2H-azirine-2-carboxylate

The structure, preferred conformers, vibrational spectrum, and photochemical behavior of the novel azirine, methyl 2-chloro-3-methyl-2H-azirine-2-carboxylate (MCMAC) were investigated in low-temperature matrixes and in the neat solid amorphous state by infrared spectroscopy and quantum-chemical calculations. Two conformers of the compound were observed in argon, krypton, and xenon matrixes, in agreement with the DFT(B3LYP)/6-311++G(d,p) and MP2/6-311++G(d,p) theoretical calculations. Both conformers were found to exhibit the carboxylic ester group in the cis conformation, differing in the arrangement defined by the O=C-C-Cl dihedral angle (cis and trans, for Ct and Cc forms, respectively). The Ct conformer was found to be the most stable conformer in the gaseous phase as well as in both argon and krypton matrixes, whereas the more polar Cc conformer became the most stable form in the xenon matrix and in the neat solid amorphous phase. In situ broadband UV (lambda > 235 nm) excitation of matrix-isolated MCMAC led to azirine ring C-C and C-N bond cleavages, the latter process corresponding to the most efficient reaction channel. The photochemical cleavage of the C-N bond had never been previously observed in the case of aliphatic 2H-azirines. Two electron withdrawing substituents (methoxycarbonyl group and chlorine atom) are connected to the azirine ring in the novel MCMAC azirine. The simultaneous presence of these two groups accelerates intersystem crossing toward the triplet state where cleavage of the C-N bond takes place. The primary photoproducts resulting from the C-N and C-C ring-opening reactions were also found to undergo further photochemical decarbonylation or decarboxylation reactions.     

Acid-base thermochemistry of gaseous oxygen and sulfur substituted amino acids (Ser, Thr, Cys, Met)

Acid-base thermochemistry of isolated amino acids containing oxygen or sulfur in their side chain (serine, threonine, cysteine and methionine) have been examined by quantum chemical computations. Density functional theory (DFT) was used, with B3LYP, B97-D and M06-2X functionals using the 6-31+G(d,p) basis set for geometry optimizations and the larger 6-311++G(3df,2p) basis set for energy computations. Composite methods CBS-QB3, G3B3, G4MP2 and G4 were applied to large sets of neutral, protonated and deprotonated conformers. Conformational analysis of these species, based on chemical approach and AMOEBA force field calculations, has been used to identify the lowest energy conformers and to estimate the population of conformers expected to be present at thermal equilibrium at 298 K. It is observed that G4, G4MP2, G3B3, CBS-QB3 composite methods and M06-2X DFT lead to similar conformer energies. Thermochemical parameters have been computed using either the most stable conformers or equilibrium populations of conformers. Comparison of experimental and theoretical proton affinities and Δ(acid)H shows that the G4 method provides the better agreement with deviations of less than 1.5 kJ mol(-1). From this point of view, a set of evaluated thermochemical quantities for serine, threonine, cysteine and methionine may be proposed: PA = 912, 919, 903, 938; GB = 878, 886, 870, 899; Δ(acid)H = 1393, 1391, 1396, 1411; Δ(acid)G = 1363, 1362, 1367, 1382 kJ mol(-1). This study also confirms that a non-negligible ΔpS° is associated with protonation of methionine and that the most acidic hydrogen of cysteine in the gas phase is that of the SH group. In several instances new conformers were identified thus suggesting a re-examination of several IRMPD spectra.     

Ab initio- and density-functional studies of conformational behaviour of N-formylmethionine in gaseous phase

The current work is a study of the conformational space of the non-ionic N-formylmethionine molecule around its seven structurally significant internal backbone torsional angles at B3LYP/6-31++G(d,p) levels of theory in the gaseous phase. The potential energy surface exploration reveals that a total of 432 different conformers would result if all the possible combinations of the internal rotations were to be considered. A set of twelve conformers of the N-formylmethionine molecule are then further analysed in terms of their relative stabilities, theoretically predicted harmonic vibrational frequencies, HOMO-LUMO energy gaps, ESP charges, rotational constants and dipole moments calculated using MP2/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels. The calculated relative energy-range of the conformers at the MP2 level is 11.08 kcal mol^?1 (1 kcal = 4.1868 kJ), whereas the same obtained at the B3LYP level is 10.02 kcal mol^?1. The results of this study provide a good account of the role of four types of intramolecular H-bonds, namely O…H—O, O…H—N, O…H—C and N…H—C, in influencing the energies of the conformers as well as their conformational and vibrational spectroscopic aspects. The relative stability order of the conformers appears to depend on the level of theory used while the vibrational frequencies calculated at the B3LYP level are in better agreement with the experimental values.     

Ab initio and density functional theory studies of the structure,gas-phase acidity and aromaticity of tetraselenosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure,acidity (free energy of deprotonation,G0) and aro-maticity of tetraselenosquaric acid (3,4-diselenyl-3-cy-dobutene-1,2-diselenone,H2C4Se4) are reported.The global minimum found on the potential energy surface of tetraselenosquaric acid presents a planar conformation.The ZZ iso-mer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very dose in energy.The optimized geometric parameters exhibit a bond length equalization relative to reference compounds,cyclobu-tanediselenone,and cydobutenediselenol.The computed aromatic stabilization energy (ASE) by homodesmotic reaction is -77.4 (MP2(fu)/6 - 311 G //RHF/6 - 311 G) and - 54.8 kJ/mol (B3LYP/6 - 311 G //B3LYP/6 -311 G).The aromaticity of tetraselenosquaric add is indicated by the calculated diamagnetic susceptibility exaltation (A) - 19.13 (CSGT(IGAEM) - RHF/6 - 311 G// RHF/6-     

Importance of entropy in the conformational equilibrium of phenylalanine: a matrix-isolation infrared spectroscopy and density functional theory study

The conformational behavior and infrared spectrum of l-phenylalanine were studied by matrix-isolation infrared spectroscopy and DFT [B3LYP/6-311++G(d,p)] calculations. The fourteen most stable structures were predicted to differ in energy by less than 10 kJ mol(-1), eight of them with abundances higher than 5% at the temperature of evaporation of the compound (423 K). Experimental results suggest that six conformers contribute to the spectrum of the isolated compound, whereas two conformers (IIb(3) and IIIb(3)) relax in matrix to a more stable form (IIb(2)) due to low energy barriers for conformational isomerization (conformational cooling). The two lowest-energy conformers (Ib(1), Ia) differ only in the arrangement of the amino acid group relative to the phenyl ring; they exhibit a relatively strong stabilizing intramolecular hydrogen bond of the O-H...N type and the carboxylic group in the trans configuration (O=C-O-H dihedral angle ca. 180 degrees ). Type II conformers have a weaker H-bond of the N-H...O=C type, but they bear the more favorable cis arrangement of the carboxylic group. Being considerably more flexible, type II conformers are stabilized by entropy and the relative abundances of two conformers of this type (IIb(2) and IIc(1)) are shown to significantly increase with temperature due to entropic stabilization. At 423 K, these conformers are found to be the first and third most abundant species present in the conformational equilibrium, with relative populations of ca. 15% each, whereas their populations could be expected to be only ca. 5% if entropy effects were not taken into consideration. Indeed, phenylalanine can be considered a notable example of a molecule where entropy plays an essential role in determining the relative abundance of the possible low-energy conformational states and then, the thermodynamics of the compound, even at moderate temperatures. Upon UV irradiation (lambda > 235 nm) of the matrix-isolated compound, unimolecular photodecomposition of phenylalanine is observed with production of CO(2) and phenethylamine.     

Theoretical Gas Phase Study of the Gauche and Trans Conformers of 1-Bromo-2-Iodoethane and Solvent Effects

This is a gas-phase study of the gauche and trans conformers of 1-bromo-2-iodoethane. The methods used are the second-order Møller-Plesset theory (MP2) and density functional theory (DFT). The functional used for the DFT method is B3LYP and the basis sets used are 6-311++G(d,p) for all atoms except that different basis sets, namely 3-21G, LANECP, CRENBL ECP, Stuttgart RLC ECP and 6-311G(d,p), have been explored for the iodine atom. The results indicate that the trans conformer is preferred. The energy difference between the gauche and trans conformers (ΔE _g?t) and related thermodynamic parameters are reported. The ΔE _g?t values are 12.50 kJ?mol^?1 (B3LYP) and 10.00 kJ?mol^?1 (MP2) with the basis sets being 6-311++G(d,p)[C,H,Br]/6-311G(d,p)[I]. The conformers of 1-bromo-2-iodoethane have also been subjected to vibrational analysis. The results from the two theoretical levels are in good agreement but they are not much affected by the basis set of the iodine atom. The study has been extended to explore solvent effects using Self-Consistent Reaction Field methods. The structural parameters of the conformers are little affected by the polarity of the solvent but ΔE _g?t decreases and the solvation Gibbs energy increases with increasing polarity of the solvent.     

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.     

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.     

Influence of H2 on the gas-phase decomposition of formic acid: a theoretical study

Gas-phase decomposition of formic acid results in final products CO + H2O and CO2 + H2. Experimentally, the CO/CO2 ratio tends to be large, in contradiction with mechanism studies, which show almost equal activation energies for dehydration and decarboxylation. In this work, the influence of H2 on the decomposition mechanism of HCOOH was explored using ab initio calculations at the CCSD(T)/6-311++G**//MP2/6-311++G** level. It was found that, in the presence of H2, the reaction channels leading to CO + H2O are more than those leading to CO2 + H2. With competitive energy, H2 addition to HCOOH can reduce the latter into HCHO, which then dissociates into CO + H2 catalyzed by H2O. Compared to trans-HCOOH, cis-HCOOH and cis-C(OH)2, conformers required for decarboxylation, are less populated due to interactions with H2.     

Theoretical investigation of hydrogen bonds between CO and HNF2, H2NF, and HNO

Ab initio quantum mechanics methods were applied to investigate the hydrogen bonds between CO and HNF2, H2NF, and HNO. We use the Hartree-Fock, MP2, and MP4(SDQ) theories with three basis sets 6-311++G(d,p), 6-311++G(2df,2p), and AUG-cc-pVDZ, and both the standard gradient and counterpoise-corrected gradient techniques to optimize the geometries in order to explore the effects of the theories, basis sets, and different optimization methods on this type of H bond. Eight complexes are obtained, including the two types of C...H-N and O...H-N hydrogen bonds: OC...HNF2(C(s)), OC...H2NF(C(s) and C1), and OC...HNO(C(s)), and CO...HNF2(C(s)), CO...H2NF(C(s) and C1), and CO...HNO(C(s)). The vibrational analysis shows that they have no imaginary frequencies and are minima in potential energy surfaces. The N-H bonds exhibit a small decrease with a concomitant blue shift of the N-H stretch frequency on complexation, except for OC...HNF2 and OC...H2NF(C1), which are red-shifting at high levels of theory and with large basis sets. The O...H-N hydrogen bonds are very weak, with 0 K dissociation energies of only 0.2-2.5 kJ/mol, but the C...H-N hydrogen bonds are stronger with dissociation energies of 2.7-7.0 kJ/mol at the MP2/AUG-cc-pVDZ level. It is notable that the IR intensity of the N-H stretch vibration decreases on complexation for the proton donor HNO but increases for HNF2 and H2NF. A calculation investigation of the dipole moment derivative leads to the conclusion that a negative permanent dipole moment derivative of the proton donor is not a necessary condition for the formation of the blue-shifting hydrogen bond. Natural bond orbital analysis shows that for the C...H-N hydrogen bonds a large electron density is transferred from CO to the donors, but for the O...H-N hydrogen bonds a small electron density transfer exists from the proton donor to the acceptor CO, which is unusual except for CO...H2NF(C(s)). From the fact that the bent hydrogen bonds in OC(CO)...H2NF(C(s)) are quite different from those in the others, we conclude that a greatly bent H-bond configuration shall inhibit both hyperconjugation and rehybridization.     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.