Gas-phase Raman spectra and the potential energy function for the internal rotation of 1,3-butadiene and its isotopologues

The gas-phase Raman spectra of 1,3-butadiene and its 2,3-d(2), 1,1,4,4-d(4), and -d(6) isotopologues have been recorded with high sensitivity in the region below 350 cm(-1) in order to investigate the internal rotation (torsional) vibration. Based on more accurate structural information, the internal rotor constants F(n) were calculated as a function of rotation angle (?). The data for all the isotopologues were then fit using a one-dimensional potential energy function of the form V = (1)/(2)∑V(n)(1 - cos ?). Initial V(n) values were based on those generated from theoretical calculations. The agreement between observed and calculated frequencies is very good, although bands not taken into account were present in the spectra. The energy difference between the trans and gauche forms was determined to be about 1030 cm(-1) (2.94 kcal/mol), and the barrier between the two equivalent gauche forms was determined to be about 180 cm(-1) (0.51 kcal/mol), which agrees well with high-level ab initio calculations. An alternative set of assignments also fits the data quite well for all of the isotopologues. For this model, the energy difference between the trans and gauche forms is about 1080 cm(-1) (3.09 kcal/mol), and the barrier between gauche forms is about 405 cm(-1) (1.16 kcal/mol).     

Computational study on kinetics and mechanisms of unimolecular decomposition of succinic acid and its anhydride

The mechanisms and kinetics of unimolecular decomposition of succinic acid and its anhydride have been studied at the G2M(CC2) and microcanonical RRKM levels of theory. It was shown that the ZsgsZ conformer of succinic acid, with the Z-acid form and the gauche conformation around the central C-C bond, is its most stable conformer, whereas the lowest energy conformer with the E-acid form, ECGsZ, is only 3.1 kcal/mol higher in energy than the ZsgsZ. Three primary decomposition channels of succinic acid producing H2O + succinic anhydride with a barrier of 51.0 kcal/mol, H2O + OCC2H3COOH with a barrier of 75.7 kcal/mol and CO2 + C2H5COOH with a barrier of 71.9 kcal/mol were predicted. The dehydration process starting from the ECGCZ-conformer is found to be dominant, whereas the decarboxylation reaction starting from the ZsgsZ-conformer is only slightly less favorable. It was shown that the decomposition of succinic anhydride occurs via a concerted fragmentation mechanism (with a 69.6 kcal/mol barrier), leading to formation of CO + CO2 + C2H4 products. On the basis of the calculated potential energy surfaces of these reactions, the rate constants for unimolecular decomposition of succinic acid and its anhydride were predicted. In addition, the predicted rate constants for the unimolecular decomposition of C2H5COOH by decarboxylation (giving C2H6 + CO2) and dehydration (giving H3CCHCO + H2O) are in good agreement with available experimental data.     

Molecular orbital studies of conformers and the barrier to internal rotation in 1,3-butadiene

The potential curve for rotation around the central bond in 1,3-butadiene has been estimated by ab initio calculations using Gaussian-type basis functions. The calculations, which also include limited geometry variation during rotation, suggest that in the SCF approximation the second stable form of the molecule is a gauche conformation rather than a cis. The predicted energy difference between the planar trans ground state and the stable gauche form is 2.7 kcal/mole and the barrier to internal rotation is found to be 6.0 kcal/mole using a (9,5) basis for carbon and 4s functions on hydrogen.     

The extremely flat torsional potential energy surface of oxalyl chloride

The conformational behavior of oxalyl chloride has been investigated using ab initio Hartree-Fock (HF) and second-order Moller-Plesset (MP2) perturbation theories, and the coupled-cluster singles and doubles method appended with a perturbative inclusion of connected triple excitations [CCSD(T)]. Correlation consistent polarized valence quadruple-zeta (cc-pVQZ) and quintuple-zeta (cc-pV5Z) basis sets were used in this research. At the cc-pVQZ and cc-pV5Z HF levels, there is no stationary point corresponding to a stable gauche conformer. On the other hand, at the cc-pVQZ and cc-pV5Z MP2 levels and with the cc-pVQZ CCSD(T) method, the gauche conformer of oxalyl chloride was found at O[Double Bond]C-C[Double Bond]O dihedral angles of 81.9 degrees , 79.4 degrees , and 83.4 degrees , respectively. At the cc-pV5Z MP2 level, the energy barrier from trans to gauche was predicted to be 0.74 kcal mol(-1) and that from gauche to trans to be 0.09 kcal mol(-1). Thus, the potential-energy surface along the O[Double Bond]C-C[Double Bond]O torsional mode is exceedingly flat. The existence of the gauche conformation is mainly due to the minimization of steric repulsion.     

Conformational equilibrium of 1,2-dichloroethane in water: comparison of PCM and RISM-SCF methods

The RISM-SCF and polarizable continuum model (PCM) approaches have been applied to study the conformational equilibrium of 1,2-dichloroethane (DCE) in water. Both the electron correlation effect and basis sets play an important role in the relative energies of the gauche and trans conformers in gas and solution phases. Both PCM and RISM-MP2 methods resulted in a consistent trend with the previous experimental and theoretical studies that the population of the gauche conformer increases in going from the gas phase to the aqueous solution. However, the PCM treatment could not describe the solvent effect completely in that the sign of the relative free energy of the gauche and trans forms is opposite to the most recent experimental and theoretical data, while the RISM-MP2 gives the right sign in the free energy difference. We found that the larger excess chemical potential gain (by ca. -4.1 kcal/mol) for the gauche conformer is large enough to result in the gauche preference of DCE in water, though it has to compensate for more solute reorganization energy (approximately 1.6 kcal/mol) and overcome the energy difference (approximately 1.6 kcal/mol) in the gas phase. The radial distribution functions between DCE and the nearest water shows that the electrostatic repulsion between chlorine and oxygen atoms is higher in the trans conformer than in the gauche one, while the attractive interaction between chlorine and hydrogen of water is higher in the gauche conformer.     

13C NMR, infrared, solvation and theoretical investigation of the conformational isomerism in 1-haloacetones (X = Cl, Br and I)

The solvent dependence of the 13C NMR spectra of chloroacetone (CA), bromoacetone (BA) and iodoacetone (IA) are reported and the 3J(CH) couplings analysed using ab initio calculations and solvation theory. In CA the energy difference (E(cis) - E(gauche)) between the cis (Cl-C-C=O 0 degrees) and gauche (Cl-C-C=O 155 degrees) conformers is 1.7 kcal mol(-1) in the vapour, decreasing to 0.8 kcal mol(-1) in CCl4 solution and to -1.0 kcal mol(-1) in the pure liquid. The conformational equilibrium, in BA, is between the more polar cis (Br-C-C=O 0 degrees) and gauche (Br-C-C=O 132 degrees) conformations. The energy difference (E(cis) - E(gauche)) is 1.8 kcal mol(-1) in the vapour, decreasing to 0.9 kcal mol(-1) in CCl4 solution and to -0.4 kcal mol(-1) in the pure liquid. The energy difference (E(cis) - E(gauche)), in IA, between the cis (I-C-C=O 0 degrees) and gauche (I-C-C=O 104 degrees) conformers is 1.1 kcal mol(-1) in the vapour phase, decreasing to 0.5 kcal mol(-1) in CCl4 solution and to -0.5 kcal mol(-1) in the pure liquid. The vapour state energy difference for BA [1.4 kcal mol(-1) at B3LYP/6-311++G(d,p)] and for IA [1.6 kcal mol(-1) at B3LYP/6-311++G(d,p)/LANL2DZ)] are in very good agreement with the above values. For CA the agreement is also satisfactory [1.4 kcal mol(-1) at B3LYP/6-311++G(d,p)].     

Computational mutation analysis of hydrogen abstraction and radical rearrangement steps in the catalysis of coenzyme B12-dependent diol dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.     

Conformational Study of 8—C—glucosyl—prunetin by Dynamic NMR Spectroscopy

By means of variable temperature NMR spectra,conformation of 8-C-glucosyl prunetin, isolated from the leaves of Dalbergia hainanensis (Leguminosae), was studied. The restricted rotation around the C(sp^3)-C(sp^2) bond in the C-glucosides isoflavonoid results in two main conformers (syn and anti). With the help of MM calculation, the preferred conformation A has H-1″ gauche to the 7-OCH3. The barrier to rotation was 18.1 kcal/mol. This result agrees with the calculated value 16.2 kcal/mol of free energy of activation for the interconversion between the conformers.     

Adjacent gauche stabilization in linear alkanes: implications for polymer models and conformational analysis

High-level ab initio quantum mechanical calculations are used to study various gauche conformational energies of n-pentane to n-decane. The destabilizing "pentane effect" (adjacent gauche states of opposite sign) for alkanes is confirmed, but the energies were found to depend slightly on chain length. In contrast, introducing an adjacent gauche of the same sign requires only 0.22-0.37 kcal/mol, approximately half of the single gauche state energy. This adjacent gauche stabilization should be taken into account when formulating or analyzing rotational isomeric models, carrying out conformational analysis, and developing force fields for alkanes, lipids, and related polymers.     

2-Silyl group effect on the reactivity of cyclopentane-1,3-diyls. Intramolecular ring-closure versus silyl migration

Generation of singlet and triplet 2-silylcyclopentane-1,3-diyls and their reactivity have been investigated in the thermal and photochemical denitrogenation of 2,3-diaza-7-silylbicyclo[2.2.1]hept-2-ene. 5-Silylcyclopentene (silyl migration product) is quantitatively obtained, while 5-silylbicyclo[2.1.0]pentane (intramolecular ring-closure product) is not detected in the denitrogenation reactions. Deuterium labeling studies clarify that 5-silylcyclopentene is formed by a suprafacial [1,2] silyl migration in singlet 2-silylcyclopentane-1,3-diyl. UDFT calculations closely reproduce the observed reactivity of the singlet diradical: The enthalpic barriers of the intramolecular ring-closure are calculated to be DeltaH++exo468 = 5.8 kcal/mol and DeltaH++endo468 = 6.7 kcal/mol, which are much higher than the energy barrier for the [1,2] silyl migration, DeltaH++468 = 2.7 kcal/mol. The notable effect of the silyl group on raising the energy barrier of the intramolecular cyclization is rationalized by an electronic configuration of the lowest singlet state of 2-silylcyclopentane-1,3-diyls.     

Ab initio study of the mechanism of singlet-dioxygen addition to hydroxyaromatic compounds: negative evidence for the involvement of peroxa and endoperoxide intermediates

In this article we report our study of two possible mechanisms of photooxidation of hydroxyaromatic compounds, involving the intermediacy of zwitterionic peroxa intermediates or 1,4-endoperoxides. To study the pathway of the first of them, as yet unexplored by theoretical methods, a simpler system composed of 1,3-butadiene-1-ol and singlet ((1)Delta(g)) dioxygen was considered first, for which calculations were carried out at the CASSCF/MCQDPT2 ab initio level, mostly with the 6-31G* basis set. The cumulative activation barrier to this reaction was found to be 20 kcal/mol and corresponded to a proton transfer (from the hydroxy oxygen atom to the attached oxygen molecule) in the cyclic zwitterionic peroxacyclopenta-3-ene-2-ol intermediate. This intermediate and the proton-transfer transition state were found to have a closed-shell character, which enabled us to estimate the corresponding activation barrier for the phenol-dioxygen system by carrying out optimization at the RHF level and single-point calculations at the MP2, CASSCF, and MCQDPT2 levels of theory. The energy barrier to the reaction was estimated to at least about 40 kcal/mol, rendering this mechanism for the phenol-oxygen system unlikely for nonpolar solvents. Similarly, calculations of the barrier to proton transfer from the 1,4-endoperoxide of phenol to its hydroperoxide were found to exceed 60 kcal/mol, eliminating such a mechanism too, which leaves only the earlier postulated mechanisms involving an initial charge or hydrogen-atom transfer to dioxygen as probable.     

Barriers to rotation in methyl formate by dynamic NMR spectroscopy and barriers to 1,3 oxygen-to-oxygen migration in methyl formate and trifluoromethyl formate by ab initio calculations

Free-energy barriers of 9.85 and 11.91 +/- 0.15 kcal/mol at -70.8 degrees C were found by dynamic NMR spectroscopy for the E-to-Z and Z-to-E conversions, respectively, of methyl formate (1) enriched in 13C to 99% for the carbonyl carbon [methyl formate 13C (2)]. These barriers are higher than the literature values reported for -53 degrees C. The free-energy barrier to 1,3 oxygen-to-oxygen migration of the methyl group in methyl formate was determined by ab initio calculations at several levels. The value of 58.7 kcal/mol obtained at the MP2/6-311+G (df,pd) level was compared to a literature barrier for this process (MINDO/3) and to barriers for related compounds. A free-energy barrier of 63.0 kcal/mol for the oxygen - to - oxygen migration of the CF3 group in trifluoromethyl formate (3) was calculated at the MP2/6-31+G level.     

Beta-trimethylsilyl cyclopropylcarbenes

Thermal decomposition of the in situ generated lithium salt of the tosylhydrazone derivative of cyclopropyl trimethylsilylmethyl ketone gave 1-cyclopropyl-1-trimethylsilylethylene, a product of exclusive silyl migration. Thermal decomposition of the sodium salts of tosylhydrazone derivatives of 1-trimethylsilylcyclopropyl alkyl ketones also gave methylenecyclopropane products derived from trimethylsilyl migration. These reactions were interpreted in terms of rapid trimethylsilyl migration to carbene-like centers that compete effectively with ring expansion processes of cyclopropylcarbenes. Computational studies (B3LYP/6-31G) suggest that cyclopropyl stabilization of carbenes is more effective than beta-trimethylsilyl stabilization. However, beta-trimethylsilyl stabilized conformations are easily attained, and these conformations can lead to silyl migrations. There are two minimum energy conformations of methyl-1-trimethylsilylcyclopropylcarbene, 27, and the rotational barrier to interconversion of these conformations (5.4 kcal/mol) is substantially lower than in the parent cyclopropylcarbene (15 kcal/mol). The onset of a stabilizing interaction in the transition state between the carbene vacant orbital with the adjacent Si-C sigma-orbital is proposed. Computational studies also show a very small (2.0 kcal/mol) barrier for trimethylsilyl migration in trimethylsilylmethyl cyclopropylcarbene, 11.     

Catalytic mechanism and product specificity of the histone lysine methyltransferase SET7/9: an ab initio QM/MM-FE study with multiple initial structures

Histone lysine methylation is emerging as an important mechanism to regulate chromatin structure and gene activity. To provide theoretical understanding of its reaction mechanism and product specificity, ab initio quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations and molecular dynamics simulations have been carried out to investigate the histone lysine methyltransferase SET7/9. It is found that the methyl-transfer reaction catalyzed by SET7/9 is a typical in-line S(N)2 nucleophilic substitution reaction with a transition state of 70% dissociative character. The calculated average free energy barrier at the MP2(6-31+G) QM/MM level is 20.4 +/- 1.1 kcal/mol, consistent with the activation barrier of 20.9 kcal/mol estimated from the experimental reaction rate. The barrier fluctuation has a strong correlation with the nucleophilic attack distance and angle in the reactant complex. The calculation results show that the product specificity of SET7/9 as a monomethyltransferase is achieved by disrupting the formation of near-attack conformations for the dimethylation reaction.     

Theoretical study of 5-aminolevulinic acid tautomerization: a novel self-catalyzed mechanism

5-Aminolevulinic acid (5ALA) is the key synthetic building block in protoporphyrin IX (PpIX), the heme chromophore in mitochondria. In this study density functional theory calculations were performed on the tautomers of 5ALA and the tautomerization reaction mechanism from its enolic forms (5-amino-4-hydroxypent-3-enoic acid and 5-amino-4-hydroxypent-4-enoic acid) to the more stable 5ALA. The hydrated form 5-amino-4,4-dihydroxypentanoic acid was also studied. The lowest energy pathway of 5ALA tautomerization is by means of autocatalysis, in that an oxygen of the carboxylic group transfers the hydrogen atom as a "crane", with an activation energy of approximately 15 kcal/mol. This should be compared to the barriers of about 35 kcal/mol for water assisted tautomerization, and 60 kcal/mol for direct hydrogen transfer. For hydration of 5ALA, the water catalyzed activation barrier is found to be approximately 35 kcal/mol, approximately 5 kcal/mol lower than direct hydration.     

Tautomeric and conformational equilibrium of 2-nitrosophenol and 9,10-phenanthrenequinonemonooxime: ab initio and NMR study

The tautomeric and conformational equilibrium of 2-nitrosophenol and 9,10-phenanthrenequinonemonooxime was studied by ab initio methods. The geometry optimizations of the structures investigated were done without any geometrical restrictions at HF/6-31G** and MP2/6-31G** levels of theory. The transition structures for tautomeric and rotameric conversions were located. To correct for electron correlation, single-point calculations were carried out up to MP4/6-311G*//MP2/6-31G* level of theory.

Ab initio calculations for 2-nitrosophenol in agreement with the available experimental data define the nitroso form as more stable. It was found that the influence of the correlation energy on the relative stabilities is smaller for the rotamers of the nitroso tautomer but substantially (4–6 kcal/mol) for the oxime forms. It was found that the barrier height of tautomerization reaction is 10.24 kcal/mol.

The structure of the 9,10-phenanthrenequinonemonooxime was studied by solid and liquid state NMR spectroscopy. Ab initio calculations in agreement with our experimental data predict that the compound exists as oxime tautomer and the syn-oxime is most stable. It was found that the solvent influence on the relative stabilities of both isomers: syn- and anti-oxime. While in chloroform solution the syn-oxime is preferred but in DMSO anti-oxime is more stable in energy.

At the MP4/6-311G*//MP2/6-31G**+ZPE level of theory the barrier of tautomerization was predicted to be 10.96 kcal/mol and the rotational barrier around the single C–O bond in the syn-oxime was found to be 7.57 kcal/mol. The rotation is facile and this explains the absence of nitroso tautomers in solution.     

Theoretical study on singlet oxygen adsorption onto surface of graphene-like aromatic hydrocarbon molecules

Recently, graphene sheet is one of interesting systems to realize novel electronic properties. Especially, interaction between graphene and adsorbed oxygen molecule is very important to control electronic condition. In this paper, we employed some aromatic hydrocarbons as simple systems of graphene sheet and ab initio MO calculations were carried out to investigate inter-molecular interaction. It is found that not triplet but singlet O₂ molecule have potential of chemisorption onto graphene surface. From the calculated potential energy surface (PES) for distance between benzene and O₂ molecules, meta-stable structure is found at about 1.5 Å with potential barrier. In the optimized structure of its meta-stable state, structural strain can be relaxed through bending of planer benzene ring. Its energy is estimated at 70.10 kcal/mol for benzene. We also estimated the strain effects for naphthalene and pyrene molecules as larger case of graphene and they were 80.85 and 72.45 kcal/mol, respectively.     

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.     

Theoretical studies on the stability of the H3O radical based on ab initio UHF-CI calculations

An UHF-CI investigation of parts of the energy surface for the H₃O radical is reported. Several types of basis sets have been used and the CI expansion included all singly and doubly replaced configurations using an UHF determinant as the reference state. H₃O, constrained-to C_3v symmetry is in the best approximation found to be 20.5 kcal/mol less stable than H₂O + H. A small local barrier of 4.6 kcal/mol for dissociation is found on the UHF level of approximation. Correlation effects lower this barrier to 3.4 kcal/mol making the existence of a quasibound state with a measurable lifetime improbable. The height of the barrier was found to be very sensitive to the detailed form of the diffuse singly occupied orbital.