Raman and infrared spectra, conformational stability, barriers to internal rotation, vibrational assignment, and ab initio calculations of 3-methyl-3-butenenitrile

The Raman (3500-30 cm⁻¹) spectra of liquid and solid and the infrared (3500-40 cm⁻¹) spectra of gaseous and solid 3-methyl-3-butenenitrile, CH₂C(CH₃)CH₂CN, have been recorded. Both cis and gauche conformers have been identified in the fluid phases but only the cis form remains in the solid. Variable temperature (−55 to −100 °C) studies of the infrared spectra of the sample dissolved in liquid xenon have been carried out. From these data, the enthalpy difference has been determined to be 163±16 cm⁻¹ (1.20±0.19 kJ mol⁻¹), with the cis conformer the more stable rotamer. It is estimated that there is 48±2% of the gauche conformer present at  25°C. A complete vibrational assignment is proposed for the cis conformer based on infrared band contours, relative intensities, depolarization ratios and group frequencies. Several of the fundamentals for the gauche conformer have also been identified. The vibrational assignments are supported by normal coordinate calculations utilizing ab initio force constants. Complete equilibrium geometries have been obtained for both rotamers by ab initio calculations employing the 6-31G(d), 6-311G(d,p), 6-311+G(d,p) and 6-311+G(2d,2p) basis sets at the levels of restricted Hartree-Fock (HF) and/or Møller-Plesset perturbation theory to the second order (MP2). Only with the 6-311G(2d,2p) and 6-311G(2df,2pd) basis sets with or without diffuse functions is the cis conformer predicted to be more stable than the gauche form. The potential energy terms for the conformational interchange have been obtained at the MP2(full)/6-311+G(2d,2p) level, and compared to those obtained from the experimental data. The results are discussed and compared to the corresponding quantities obtained for some similar molecules.     

Isomerization barriers in bis(4H-thiopyran) and in bithioxanthenes

The inversion and rotation mechanisms for the isomerization of Feringa’s bithioxanthenes existing in two conformations, up/up and up/down, have been calculated at the B3LYP/6-31G(d) and B3LYP/6-311++G(d,p) levels. The inversion mechanism that maintains the double bond nature of the central bond is a classical one but the rotation mechanisms that break the double bond to form a biradical needs to explore the singlet and triplet states. To do this we have removed the four fused phenyl rings of bithioxanthene and calculated at the CASSCF and CASPT2 levels bis(4H-thiopyran) proving that B3LYP calculations yield reasonable results for the rotation barriers.     

Conformational study of the structure of free 12-thiacrown-4 and some of its cation metal complexes

Conformational search of 12-thiacrown-4, 12t4, was performed using the CONFLEX method and the MMFF94S force field whereby 156 conformations were predicted. Optimized geometries of the 156 predicted conformations were calculated at the HF, B3LYP, CAM-B3LYP, M06, M06L, M062x and M06HF levels using the 6-311G** basis set. The correlation energy was recovered at the MP2 level using the same 6-311G** basis set. Optimized geometries at the MP2/6-311G** level and G3MP2 energies were calculated for some of the low energy conformations. The D ₄ conformation was predicted to be the ground state conformation at all levels of theory considered in this work. Comparison between the dihedral angles of the predicted conformations indicated that for the stability of 12t4, a SCCS dihedral angle of 180° requirement is more important than a gauche CSCC dihedral angle requirement. Conformational search was performed also for the 12t4?CAg⁺, Bi³⁺, Cd²⁺, Cu⁺ and Sb³⁺ cation metal complexes using the CONFLEX method and the CAChe-augmented MM3 and MMFF94S force fields. Conformations with relative energies less than 10?kcal/mol at the MP2/6-31+G*//HF/6-31+G* level, with double zeta quality basis set on the metal cations, were considered for computations at the same levels as those used for free 12t4, using also the 6-311G** basis set. The cc-pVTZ-pp basis set was used for the metal cations. The predicted ground state conformations of the 12t4?CAg⁺, Bi³⁺, Cd²⁺, Cu⁺ and Sb³⁺ cation metal complexes are the C ₄, C ₄, C ₄, C _2v and C ₄ conformations, respectively. This is in agreement with the experimental X-ray data for the 12t4?CAg⁺ and Cd²⁺ cation metal complexes, but experimentally by X-ray, the 12t4?CBi³⁺ and Cu⁺ cation metal complexes have C _s and C ₄ structures, respectively.     

Intramolecular hydrogen bond in the push–pull CF3-aminoenones: DFT and FTIR study,NBO analysis

Postulated conformers of trifluoromethylated β-aminoenones stabilized by intramolecular NH?O and N?HO bonds were studied by IR and NMR spectroscopy and evaluated with quantum chemical calculations (B3LYP/6-311+G(d,p), MP2/6-311+G(d,p)//B3LYP/6-311+G(d,p) and MP2/6-31G(d,p)) and NBO analysis. The influence of the nature of EWG, substituents at the nitrogen atom and double bond, and of orbital interactions of heteroatoms and double bonds in these structures on the proton affinity of basic and acid centers, strength of hydrogen bonds, and the energy of tautomeric transfers is discussed. The theoretical results agree satisfactorily with the experimental observations.     

FT-IR, FT-Raman vibrational spectra and molecular structure investigation of 2-chloro-4-methylaniline: A combined experimental and theoretical study

In this work, the experimental and theoretical vibrational spectra of 2-chloro-4-methylaniline (2Cl4MA, C₇H₈NCl) were studied. FT-IR and FT-Raman spectra of 2Cl4MA in the liquid phase have been recorded in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock (HF) and density functional method (B3LYP) with the 6-31G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p) and 6-311G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p) basis sets. The vibrational frequencies have been calculated and scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The DFT-B3LYP/6-311++G(d,p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 2Cl4MA. The optimized geometric parameters (bond lengths and bond angles) were compared with experimental values of aniline and p-methylaniline molecules.     

Methyl and trifluoromethyl vinyl sulfoxides. Steric and electronic structure

Potential functions of internal rotation about the $C_{sp^2 } $ -S bond in H₂C=CHSOCY₃ (Y = H or F) were determined and stationary points were identified by vibrational analysis at the MP2(full)/6-31+G(d), B3PW91/6-31+G(d), and B3PW91/6-311+G(3df,p) levels. Energetically favorable conformations were established, and rotation barriers and molecular geometry parameters were evaluated. Wave functions [MP2(full)/6-31+G(d)] were analyzed by the NBO method. Energies and dipole moments of bond and lone-pair orbitals, principal types of donor-acceptor interactions, bond orders, and atomic charges were determined.     

Density functional theory studies on the dissociation energies of metallic salts: relationship between lattice and dissociation energies

The formation and physicochemical properties of polymer electrolytes strongly depend on the lattice energy of metal salts. An indirect but efficient way to estimate the lattice energy through the relationship between the heterolytic bond dissociation and lattice energies is proposed in this work. The heterolytic bond dissociation energies for alkali metal compounds were calculated theoretically using the Density Functional Theory (DFT) of B3LYP level with 6‐311+G(d,p) and 6‐311+G(2df,p) basis sets. For transition metal compounds, the same method was employed except for using the effective core potential (ECP) of LANL2DZ and SDD on transition metals for 6‐311+G(d,p) and 6‐311+G(2df,p) calculations, respectively. The dissociation energies calculated by 6‐311+G(2df,p) basis set combined with SDD basis set were better correlated with the experimental values with average error of ca. ±1.0% than those by 6‐311+G* combined with the LANL2DZ basis set. The relationship between dissociation and lattice energies was found to be fairly linear (r>0.98). Thus, this method can be used to estimate the lattice energy of an unknown ionic compound with reasonably high accuracy. We also found that the dissociation energies of transition metal salts were relatively larger than those of alkaline metal salts for comparable ionic radii. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 827–834, 2001     

Reactions of dichlorocarbene,dichlorosilylene, and dichlorogermylene with carboranes(12). A theoretical study

Insertion of CCl₂, SiCl₂, and GeCl₂ into the B-H and C-H bonds of the carborane(12) molecules as well as into the B-Hal, C-Hal bonds of some halogen-substituted ortho-carboranes are studied at the B3LYP/6-311+G(d,p) and M06-2X/6-311+G(d,p) levels of theory. Thermodynamic parameters of the reactions are calculated and their variation for the isomeric carboranes(12) and for a series of carbenoid species is shown. For several insertion reactions the energies of activation are calculated based on the model of three-center synchronous transition state. The energy barriers are regularly varied depending on the position of the B-H bond in the carborane molecule.     

优化几何构型对高级别从头算能量的影响

对《ＣＲＣ物理与化学手册》（第７７版）中第三周期以前的无机双原子分子,当其理论优化构型的相对误差大于２％时,分别在实验构型和最大偏差的理论构型下,计算了ＱＣＩＳＤ（Ｔ）／６－３１１＋Ｇ（３ｄｆ,２ｐ）能量并作了比较。结果表明,大多数能量的差别在４．２ｋＪ．ｍｏｌ＾－１以内。由此说明,目前一般采用的构型优化理论方法,多数情况下不至于明显影响单点高级别从头算的计算精度。同时还发现,Ｇ２（ＱＣＩ）的高级     

Redox potentials of dopamine and its supramolecular complex with aspartic acid

Dopamine (DA) can be oxidized to dopamine quinone (DAquinone) through a one-step, two-electron redox reaction. The electron transfer property of DA and its supramolecular complex with aspartic acid (Asp) has been investigated by the theoretical calculations. We calculated the standard redox potentials (E ^o) of DA/DAquinone at the MP2/6-31G(d,p)//B3LYP/6-31G(d,p), MP2/6-31+G(d,p)//B3LYP/6-31+G(d,p), MP2/6-31G(d,p)//B3LYP/6-311G(d,p), and MP2/6-311+G(d,p)//B3LYP/6-311+G(d,p) levels. Comparing the experimental value, the redox potentials of DA/DAquinone obtained at MP2//B3LYP/6-311G(d,p) and MP2//B3LYP/6-311+G(d,p) levels can be considered as the upper and lower estimates. DA can form supramolecular complex (DA-Asp) with Asp through hydrogen bond (H-bond). Therefore, the values of 0.631 and 0.628 V obtained at MP2//B3LYP/6-311G(d,p) and MP2//B3LYP/6-311+G(d,p) levels for DA-Asp/DAquinone-Asp can be proposed as the upper and lower estimates of a probable (about 0.630 V) value of the corresponding redox potential. The calculated E ^o values of DA-Asp/DAquinone-Asp at the four theoretical levels are upper than those of DA/DAquinone, which indicates that the formation of H-bonds weaken the electron-donating ability of DA.     

Ab initio studies of the allylic hydroxylation: DFT calculation on the reaction of 2-methyl-2-butene with selenium dioxide

The mechanism of the allylic oxidation of 2-methyl-2-butene with selenium dioxide has been investigated by ab initio quantum mechanics. Transition states for two major steps (an ene reaction and a [2,3]-sigmatropic rearrangement) of this reaction have been optimized by the B3LYP/6-311+G(d,p) method. A comparison of the energies of the transition states shows that the anti-endo and syn-endo approaches are the efficient routes in the ene reaction and the methyl (C4) group is sited in a pseudo-equatorial environment in cyclic transition states during the [2,3]-rearrangement. Calculations also show the kind of the terminal alkyl (C4) substituents may control (E)-selectivity in the formation of the allylic alcohols.     

Bond dissociation energies and equilibrium structures of Cu+(MeOH)x, x = 1-6, in the gas phase: competition between solvation of the metal ion and hydrogen-bonding interactions

The solvation of Cu+ by methanol (MeOH) was studied via examination of the kinetic energy dependence of the collision-induced dissociation of Cu+(MeOH)x complexes, where x = 1-6, with Xe in a guided ion beam tandem mass spectrometer. In all cases, the primary and lowest-energy dissociation channel observed is the endothermic loss of a single MeOH molecule. The primary cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) after accounting for the effects of multiple ion-neutral collisions, kinetic and internal energy distributions of the reactants, and lifetimes for dissociation. Density functional theory calculations at the B3LYP/6-31G* level are performed to obtain model structures, vibrational frequencies, and rotational constants for the Cu+(MeOH)x complexes and their dissociation products. The relative stabilities of various conformations and theoretical BDEs are determined from single-point energy calculations at the B3LYP/6-311+G(2d,2p) level of theory using B3LYP/6-31G*-optimized geometries. The relative stabilities of the various conformations of the Cu+(MeOH)x complexes and the trends in the sequential BDEs are explained in terms of stabilization gained from sd hybridization, hydrogen-bonding interactions, electron donor-acceptor natural bond orbital stabilizing interactions, and destabilization arising from ligand-ligand repulsion.     

A computational study of conformational interconversions in 1,4-dithiacyclohexane (1,4-dithiane)

Ab initio molecular orbital theory with the 6-31G(d), 6-31G(d,p), 6-31+G(d), 6-31+G(d,p), 6-31+G(2d,p), 6-311G(d), 6-311G(d,p), and 6-311+G(2d,p) basis sets and density functional theory (BLYP, B3LYP, B3P86, B3PW91) have been used to locate transition states involved in the conformational interconversions of 1,4-dithiacyclohexane (1,4-dithiane) and to calculate the geometry optimized structures, relative energies, enthalpies, entropies, and free energies of the chair and twist conformers. In the chair and 1,4-twist conformers the C-Hax and C-Heq bond lengths are equal at each carbon, which suggest an absence of stereoelectronic hyperconjugative interactions involving carbon-hydrogen bonds. The 1,4-boat transition state structure was 9.53 to 10.5 kcal/mol higher in energy than the chair conformer and 4.75 to 5.82 kcal/mol higher in energy than the 1,4-twist conformer. Intrinsic reaction coordinate (IRC) calculations showed that the 1,4-boat transition state structure was the energy maximum in the interconversion of the enantiomers of the 1,4-twist conformer. The energy difference between the chair conformer and the 1,4-twist conformer was 4.85 kcal/mol and the chair-1,4-twist free energy difference (deltaG degrees (c-t)) was 4.93 kcal/mol at 298.15 K. Intrinsic reaction coordinate (IRC) calculations connected the transition state between the chair conformer and the 1,4-twist conformer. This transition state is 11.7 kcal/mol higher in energy than the chair conformer. The effects of basis sets on the 1,4-dithiane calculations and the relative energies of saturated and unsaturated six-membered dithianes and dioxanes are also discussed.     

Theoretical and experimental study of the absolute configuration of helical structure of (2R,3S)-Rubiginone A2 analog

A novel analogue of (2R,3S)-Rubiginone A₂ was synthesized as a chiral helical model compound via an eight-step procedure (2.7% overall yield). Quantum methods, such as density functional theory (DFT) at different basis sets of 6-311+(d), 6-311++G(2d,p), were used to compute its optical rotation and electronic circular dichroism at the B3LYP/6-311++G(2d,p) level in the gas phase and in solution using PCM model, respectively. UV corrections were performed in electronic circular dichroism (ECD) simulations to match the experimental ECD well. The suitable computational methods, e.g., B3LYP/6-311++G(2d,p)//B3LYP/6-311++G(2d,p) in the gas phase using zero-point energy in Boltzmann statistics, were found and suggested for optical rotation and circular dichroism computations that can be used for absolute configuration determination of chiral helical compounds.     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

Neutralization-Reionization of alkenylammonium cations: An experimental and ab initio study of intramolecular N-H … C=C interactions in cations and hypervalent ammonium radicals

A series of isomeric hexenylammonium and hexenyldimethylammonium cations were neutralized by collisional electron transfer in the gas phase in an attempt to generate hypervalent ammonium radicals. The radicals dissociated completely on the 4.8–5.4 µs time scale. Radicals in which the hexene double bond was in the 3-, 4-, and 5-positions dissociated by competitive N-H and N=C bond cleavages. Allylic 2-hexen-1-ylammonium and 2-hexen-1-yldimethylammonium radicals underwent predominant cleavages of allylic N-C bonds. Deuterium labeling experiments revealed no intramolecular hydrogen transfer from the hypervalent ammonium group to the hexene double bond. Ab initio and density functional theory calculations showed that alkenylammonium and alkenylmethyloxonium ions preferred hydrogen bonded structures in the gas phase. The stabilization through intramolecular H bonding in 3-buten-1-ylammonium and 3-buten-1-yl methyloxonium ions was calculated by B3LYP/6-311G(2d,p) at 26 and 18 kJ mol^?1, respectively. No intramolecular hydrogen bonding was found for the allylammonium ion. The hypervalent 3-buten-1-yl-methyloxonium radical was calculated to be unbound and predicted to dissociate exothermically by O-H bond cleavage. This dissociation may provide kinetic energy for the hydrogen atom to overcome a small energy barrier for exothermic addition to the double bond. The 3-butten-1-ylammonium and allylammonium radicals were found to be bound and preferred gauche conformations without intramolecular hydrogen bonding. Vertical neutralization of alkenylammonium ions was accompanied by small Franck-Condon effects. The failure to detect stable or metastable hypervalent alkenylammonium radicals was ascribed to the low activation barriers to exothermic dissociations by N-H and N-C bond cleavages.     

CH_2CO+CH_2的多通道反应的理论研究

利用密度泛函(DFT)和自然键轨道理论(NBO)及高级电子耦合簇[CCSD(T)]和电子密度拓扑(AIM)方法,对单重态和三重态CH2与CH2CO反应的微观机理进行了研究.在B3LYP/6-311+G(d,p)水平上优化了反应通道各驻点的几何构型.在CCSD(T)/6-311+G(d,p)水平上计算了各物种的单点能量,并对总能量进行了校正.计算表明,单重态CH2与CH2CO的C-H键可发生插入反应,与C=C、C=O可发生加成反应,存在三条反应通道,产物为CO和C2H4,从能量变化和反应速控步骤能垒两方面考虑,反应II更容易发生.对反应通道中的关键点进行了自然键轨道及电子密度拓扑分析.三重态CH2与CH2CO的反应存在三条反应通道,一条是与C-H键的插入反应,另一条是三重态CH2与C=C发生加成反应,产物为CO和三重态C2H4,通道II势垒较低,更容易发生.最后一条涉及双自由基的反应活化能最大,最难发生.     

Methyl and trifluoromethyl vinyl sulfones: Steric and electronic structure

The potential functions of internal rotation around the $C_{sp^2 } $ -S bond in the compounds H₂C= CHSO₂CY₃ (Y = H or F) were obtained; the stationary points were revealed and identified by solving vibrational problems [MP2(full)/6-31+G(d), B3PW91/6-31+G(d), B3PW91/6-311+G(3df,p)]. The energetically favorable conformations were determined, and the internal rotation barriers and geometric parameters of the molecules were estimated. The wave functions [MP2(full)/6-31+G(d)] were analyzed by the NBO method. The energies, populations, and dipole moments of natural orbitals of bonds and lone electron pairs, the main types of donor-acceptor interactions, the bond orders (Wiberg indices), and the natural charges on atoms were estimated. The steric and electronic structures of vinyl sulfoxides and sulfones are compared.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out