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.     

An exploration of conformational search of leucine molecule and their vibrational spectra in gas phase using ab initio methods

An extensive exploration of the conformational space has been carried out to characterize all possible gas phase structures of leucine. A total of 324 unique trial structures for canonical leucine were generated by considering all possible combinations of single bond rotamers. All trial structures were optimized at the B3LYP/6-311G* level of the DFT method. A total of 77 unique and stationary canonical conformers were found. Further, 15 most stable conformers were reoptimized at B3LYP/6-311++G** level and their respective relative energies, vertical ionization energies, hydrogen bonding patterns, rotational constants and dipole moments were calculated. A single point energy calculations for leucine conformers have also been done at both B3LYP/6-311++G(2df, p) and MP2/6-311++G(2df, p) levels. The good agreement between our estimates of rotational constants for two most stable conformers and available experimental measurements supports the reliability of the B3LYP/6-311++G** level of theory for describing the conformational behavior of leucine molecule. The proton affinity and gas phase basicity were also determined. Using the statistical approach, conformational distributions at various temperatures have also been performed and analyzed. Vibrational spectra were also calculated. It is also observed that zwitterions of leucine are not stable in gas phase.     

Conformation equilibrium of 3-(hydroxymethyl)piperidine in solvents with different polarity

Quantum-chemical calculations of the 3-(hydroxymethyl)piperidine molecule conformers were performed at the B3LYP/6-31+G** level of theory, and four most stable conformations with different relative orientation of CH₂OH and N–H groups were determined. The optimized structures, vibration frequencies, and band intensities in the spectra of the conformers were obtained. The conformational equilibria of the most stable rotational isomers in solvents of different polarity was studied within the polarizable continuum model. According to the results of calculations, the conformational equilibrium in solution is substantially changed on varying the solvent polarity. This conclusion was confirmed by comparison with IR absorption spectra of 3-(hydroxymethyl)piperidine solutions in carbon tetrachloride in the region of ОН-stretchings.

     

Proton and sodium cation affinities of harpagide: a computational study

The aim of this work was to estimate the proton and sodium cation affinities of harpagide (Har), an iridoid glycoside responsible for the antiinflammatory properties of the medicinal plant Harpagophytum. Monte Carlo conformational searches were performed at the semiempirical AM1 level to determine the most stable conformers for harpagide and its protonated and Na+-cationized forms. The 10 oxygen atoms of the molecule were considered as possible protonation and cationization sites. Geometry optimizations were then refined at the DFT B3LYP/6-31G level from the geometries of the most stable conformers found. Final energetics were obtained at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G level. The proton and sodium ion affinities of harpagide have been estimated at 223.5 and 66.0 kcal/mol, respectively. Since harpagide mainly provides HarNa+ ions in electrospray experiments, the DeltarG298 associated with the reaction of proton/sodium exchange between Har and methanol, MeOHNa+ + HarH+ --> MeOH2+ + HarNa+ (1), has been calculated; it has been estimated to be 1.9 kcal/mol. Complexing a methanol molecule to each reagent and product of reaction 1 makes the reaction become exothermic by 1.7 kcal/mol. These values are in the limit of the accuracy of the method and do not allow us to conclude definitely whether the reaction is endo- or exothermic, but, according to these very small values, the cation exchange reaction is expected to proceed easily in the final stages of the ion desolvation process.     

Ab initio molecular orbital study of energies and conformers of 3,4-dihydro-1,2-dithiin, 3,6-dihydro-1,2-dithiin, 4H-1,3-dithiin,and 2,3-dihydro-1,4-dithiin

Optimized geometries and energies for 3,4-dihydro-1,2-dithiin ( 1 ), 3,6-dihydro-1,2-dithiin ( 2 ), 4H-1,3-dithiin ( 3 ), and 2,3-dihydro-1,4-dithiin ( 4 ) were calculated using ab initio 6-31G* and MP2/6-31G*//6-31G* methods. At the MP2/6-31G*//6-31G* level, the half-chair conformer of 4 is more stable than those of 1 , 2 , and 3 by 2.5, 3.5, and 3.6 kcal/mol, respectively. The half-chair conformers of 1 , 2 , 3 , and 4 are 2.9, 7.1, 2.0, and 5.6 kcal/mol, respectively, more stable than their boat conformers. The calculated half-chair structures of 1 – 4 are compared with the calculated chair conformer of cyclohexane and the half-chair structures for cyclohexene, 3,4-dihydro-1,2-dioxin ( 5 ), 3,6-dihydro-1,2-dioxin ( 6 ), 4H-1,3-dioxin ( 7 ), and 2,3-dihydro-1,4-dioxin ( 8 ). © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1064–1071, 1998     

Conformational analysis and vibrational assignments of 2,2,3,3-tetrafluoro-1-propanol CHF2CF2CH2OH as three-rotors system

The conformational stability of 2,2,3,3-tetrafluoro-1-propanol was investigated by DFT-B3LYP/6-311+G** and ab initio MP2/6-311+G** calculations. The calculated potential energy curves of the molecule at DFT-B3LYP level were consistent with five distinct minima that correspond to gauche(-)-gauche-gauche (G1gg), trans-trans-gauche (Ttg), trans-gauche-gauche (Tgg), trans-gauche-gauche(-) (Tgg1) and gauche(-)-gauche-trans (G1gt) conformers in the order of decreasing relative stability. The equilibrium constants for the conformational interconversion of 2,2,3,3-tetrafluoro-1-propanol were calculated and found to correspond to an equilibrium mixture of about 38% G1gg, 28% Ttg, 13% Tgg, 11% Tggt and 10% G1gt conformations at 298.15K. The vibrational frequencies of 2,2,3,3,-tetrafluoro-1-propanol in its five stable forms were computed at B3LYP level and complete vibrational assignments were made based on normal coordinate calculations and comparison with experimental data of the molecule.     

The Role of Intramolecular Hydrogen Bonds vs. Other Weak Interactions on the Conformation of Hyponitrous Acid and Its Mono- and Dithio-Derivatives

High level ab initio and density functional theory calculations have been carried out to investigate the relative stability of the different conformers of hyponitrous acid and its mono- and dithio-derivatives. Geometries and vibrational frequencies were obtained at the B3LYP/6-311+G(d,p) level and final energies through B3LYP/6-311++G(3df,2pd) single point calculations. The reliability of this theoretical scheme has been assessed by comparing these DFT results with those obtained at the G3 level of theory, for some suitable cases. The cis conformers of hyponitrous acid and its mono- and dithio-derivatives are systematically more stable than the trans ones because in the cis conformation a dative interaction between the nitrogen-lone pairs and the σ_NX^* antibonding orbital is significantly favored. Quite interestingly, in general, the conformers presenting an intramolecular hydrogen bond (IHB) are not the global minima of the corresponding potential energy surfaces and only for hyponitrous acid the conformer with a OH ⋅s O IHB is slightly more stable than the cis conformer without IHB. The low stability of the tautomers with IHB is closely related with another weak intramolecular interaction which involves the lone-pairs of the chalcogen atoms and the π_NN^* antibondig orbital, and which is significantly perturbed when the IHB is formed.     

CF3CH2CF2CH3（HFC-365mfc）与Cl原子反应的微观机理及动力学性质

采用密度泛函理论方法 M06-2X结合6-31+G(d,p)基组研究了CF3CH2CF2CH3与Cl原子反应的反应机理.计算获得了CF3CH2CF2CH3的两种可区分的稳定几何构象RC1和RC2以及与它们相对应的8条氢提取反应通道和2条取代反应通道.运用改进的正则变分过渡态理论(ICVT)并结合小曲率隧道效应校正(SCT),在M06-2X/6-31+G(d,p)水平上计算了各氢提取通道的速率常数,并由Boltzmann配分函数得到总包反应的速率常数kT(cm3.molecule-1.s-1).计算结果表明,体系的总反应速率常数与已有实验值相吻合,进而给出了该反应在200～1000 K温度区间内反应速率常数kT的三参数表达式kT=1.88×10-22T3.76.exp(-1780.69/T),并讨论了两种构象RC1和RC2对总反应的贡献及各构象中氢提取发生在—CH3或—CH2—基团上的位置选择性.此外,由于缺少相关反应物及产物自由基标准生成焓ΔHf,298 K的数据,利用等化学键法估算了在上述物种的标准生成焓.     

Molecular structure,conformational stability,energetic and intramolecular hydrogen bonding in ground,and electronic excited state of 3-mercapto propeneselenal

In the present work, a conformational analysis of 3-mercapto propeneselenal is performed using several computational methods, including DFT (B3LYP), MP2, and G2MP2. At the DFT and G2MP2 levels the most stable conformers of title compound are characterized by an extended backbone structure, minimizing the steric repulsions between the sulfur and selenium lone pairs. Two conformers exhibit hydrogen bonding. This feature, although not being the dominant factor in energetic terms, appears to be of foremost importance to define the geometry of the molecule. The influence of the solvent on the stability order of conformers and the strength of intramolecular hydrogen bonding was considered using the PCM, SCI–PCM, and IEF–PCM methods. The results of analysis by quantum theory of “Atoms in Molecules” and natural bond orbital method fairly support the DFT results. The calculated HOMO and LUMO energies showed that charge transfer occurs within the molecule. Further verification of the obtained transition state structures was implemented via intrinsic reaction coordinate analysis. Calculations of the ¹H NMR chemical shift at GIAO/B3LYP/6–311++G** levels of theory are also presented. The excited-state properties of intramolecular hydrogen bonding in hydrogen-bonded systems have been investigated theoretically using the time-dependent density functional theory method.     

Analysis of the Electronic Structure of 1-(4-Aminobutyl)-3-methylimidazolium Chloride via DFT Method

Density functional theory and natural bond orbital analysis are applied to study the microelectronic structural properties of 1-(4-aminobutyl)-3-methylimidazolium chloride(C_4NH_2C_1imCl). We obtain 48 stable cation conformers at the B3LYP/6-311++G(d, p) level of theory. When the Cl anion is located at different positions around the two stable cation conformers,optimization and frequency are calculated at the same level of theory, and seven stable C_4NH_2C_1imCl structures are obtained. Ion-pair association energies, natural population analysis(NPA) charges, and the second-order stabilization energies are also calculated. Results show that three energy-degenerated states are present in seven conformers. The calculated ion-pair association energy, NPA charges, and second-order stabilization energy values are the same as those with the same energy level conformers. Finally, H bond formation in the molecule is analyzed by the interactions of frontier molecular orbitals. The lone pair of N(11) atom in the amino group exhibits higher reactivity than that of the N(4) atom in the ring.     

Conformational analysis and comparison between theoretical and experimental vibrational spectra for chloroacetyl isocyanate

The conformational stability and vibrational infrared and Raman spectra of chloroacetyl isocyanate (CH2ClCONCO) were investigated by ab initio MP2 and density functional B3LYP calculations using the 6-311 + + G** basis set. From the potential energy scans of the internal rotations of both the halomethyl and the isocyanate rotors, chloroacetyl isocyanate was predicted to exist predominantly in a mixture of the cis-cis (chlorine atom and NCO group eclipse C=O bond) and the gauche-cis (one hydrogen atom and NCO group eclipse C=O bond) conformations with a comparable relative stability. The vibrational wavenumbers of each of the two conformers of the molecule were computed at DFT-B3LYP/6-311 + + G** level. Normal coordinate calculations were carried out to obtain the potential energy distributions (PED) among the symmetry coordinates of the normal modes for each of the stable conformers of chloroacetyl isocyanate. The theoretical vibrational assignments are compared with experimental ones and a ratio of observed/calculated wavenumbers of about 0.97-1.04 was obtained.     

Ab initio studies of aspartic acid conformers in gas phase and in solution

Systematic and extensive conformational searches of aspartic acid in gas phase and in solution have been performed. For the gaseous aspartic acid, a total of 1296 trial canonical structures and 216 trial zwitterionic structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were optimized at the B3LYP/6-311G* level and then subjected to further optimization at the B3LYP/6-311++G** level. A total of 139 canonical conformers were found, but no stable zwitterionic structure was found. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, and vertical ionization energies of the canonical conformers were determined. Single-point energies were also calculated at the MP2/6-311++G** and CCSD/6-311++G** levels. The equilibrium distributions of the gaseous conformers at various temperatures were calculated. The proton affinity and gas phase basicity were calculated and the results are in excellent agreement with the experiments. The conformations in the solution were studied with different solvation models. The 216 trial zwitterionic structures were first optimized at the B3LYP/6-311G* level using the Onsager self-consistent reaction field model (SCRF) and then optimized at the B3LYP/6-311++G** level using the conductorlike polarized continuum model (CPCM) SCRF theory. A total of 22 zwitterions conformers were found. The gaseous canonical conformers were combined with the CPCM model and optimized at the B3LYP/6-311++G** level. The solvated zwitterionic and canonical structures were further examined by the discrete/SCRF model with one and two water molecules. The incremental solvation of the canonical and zwitterionic structures with up to six water molecules in gas phase was systematically examined. The studies show that combining aspartic acid with at least six water molecules in the gas phase or two water molecules and a SCRF solution model is required to provide qualitatively correct results in the solution.     

Vibrational spectroscopic and conformational studies of 1-(4-pyridyl)piperazine

Possible stable conformers of the 1-(4-pyridyl)piperazine (1-4pypp) molecule were experimentally and theoretically studied by FT-IR and Raman spectroscopy. FT-IR and Raman spectra were recorded in the region of 4000–200 cm⁻¹. Optimized geometric structures related to the minimum on the potential energy surface were investigated by the B3LYP hybrid density functional theory method using the 6-31G(d) basis set. Comparison of the experimental and theoretical results indicates that the density functional B3LYP method provides satisfactory results for the prediction of vibrational wavenumbers and structural parameters and equatorial-equatorial (e-e) isomer is supposed to be the most stable form of the 1–4pypp molecule.     

A Computational Study of Chair–Twist Energy Differences and the Chair–Chair Conformational Interconversion in Tetrahydro-2H-thiopyran (Thiacyclohexane,Thiane)

Ab initio molecular orbital theory with the LANL2DZ, 3-21G, 6-31G(d), 6-31+G(d), 6-31+G(d,p), 6-311+G(d,p),6-31G(2d), 6-31G(3d), and 6-311G(d,p) basis sets and density functional theory (B3P86, B3LYP, B3PW91) have been used to calculate the structures, relative energies, enthalpies, entropies, and free energies of the chair, 1,4-twist, and 2,5-twist conformers of tetrahydro-2H-thiopyran (tetrahydrothiopyran, thiacyclohexane, thiane, pentamethylene sulfide). All levels of theory calculated similar energy values and the effect of basis sets on the calculated energies was small. The chair conformer of tetrahydro-2H-thiopyran was 5.27 kcal/mol more stable than the 1,4-twist conformer, which was slightly more stable (0.81 kcal/mol) than the 2,5-twist conformer. The chair–1,4-twist and chair–2,5-twist free energy differences ( G°_{c – t}) were 5.44 and 5.71 kcal/mol, respectively. Intrinsic reaction coordinate [IRC, minimum-energy path (MEP)] calculations connected the transition state between the chair and the 2,5-twist conformers. This transition state is 9.73 kcal/mol higher in energy than the chair conformer and the energy differences between the chair and the 1,4-boat and 2,5-boat transition states were 8.07 and 6.38 kcal/mol, respectively. Stereoelectronic hyperconjugative interactions were observed in the chair, 1,4-twist, and 2,5-twist conformers of tetrahydro-2H-thiopyran. The stereoelectronic hyperconjugative effects in the chair conformer of tetrahydro-2H-thiopyran have been compared to those in the respective chair conformers of tetrahydro-2H-pyran, tetrahydro-2H-selenane, and tetrahydro-2H-tellurane.     

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

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

On the conformation of the propranolol molecule

The structure of the propranolol molecule has been optimized within the AM1 and PM3 semiempirical framework followed by ab initio HF/6-31G^* refinement. On each calculation level the conformational space was sampled to search for the lowest-energy conformer(s) from among a few hundreds of conformers at the semiempirical step and next from among a few dozens of conformers at the ab initio level. Finally, five stable conformers were found; each stabilized by one or two of the three possible hydrogen bonds. The geometrical and electronic parameters were established and found to differ only slightly in the structures with the hydrogen bond either present or not.     

Ab initio study on the internal rotation of five π-conjugated hydrocarbons at MP2 level

The potential-energy curves of internal rotation were calculated for 1,3-butadiene at the MP2/6-311G** level, for isoprene and 1,3-pentadiene at the MP2/6-311G* level, and for 2,3-dimethyl-1,3-butadiene and styrene at the MP2/6-31G* level. The geometries of the energy minima (stable conformers) and maxima (transition states) on the curves are completely optimized. For butadiene and its methyl derivatives, two stable rotamers, s-trans and gauche conformers, are obtained. s-trans forms have the lowest energies and gauche conformers twisted by 39.9°–48.3° around the central bond of the butadiene skeleton are, on average, 9.8 kJ/mol above the trans forms. s-cis forms are rotational transition states. The computed gauche–cis barriers range from 4.30 to 11.70 kJ/mol. The regular effects of methyl substitutions at the end and central carbons are found. For styrene, the planar form is calculated to be a saddle point which is only about 1 kJ/mol higher in total energy than a twisted minimum, in which the torsional angle between the phenyl and vinyl planes is 27.4°. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 659–667, 1998     

Ab initio study of the conformational equilibrium of ethylene glycol

Summary The conformational equilibrium of ethylene glycol (CH₂OHCH₂OH) has been examined by performing geometry optimizations at the 6-31G*, MP2/6-31G* and 6-31G** levels. Final energies have been calculated at the MP3 level with the optimized geometries. The two most stable conformers are atGg andgGg but it is verified that the inclusion of electronic correlations reduces their energy difference of 0.6 kcal/mol at the HF level to less than 0.2 kcal/mol. The possible coexistence of the two most stable conformers is in agreement with some previous studies of Frei et al. For thetXg conformer a detailed analysis of the intramolecular potential as a function of rotation around the C-C bond is also reported.