A DFT study on the double bond migration of butene catalyzed by ionic pair of 1-ethyl-3-methyl-imidazolium fluoride

The double bond migration of butene catalyzed by 1-ethyl-3-methyl-imidazolium fluoride (EmimF) has been studied using quantum chemical method. The geometries of reactant, transition state and product for the isomerization have been optimized by density functional theory (DFT) at the B3PW91/6-31G(d,p), 6-311++G(d,p) and aug-cc-PVDZ levels. The computed results show that the 4-H atom on imidazole ring of EmimF has a good catalytic activity to the double bond migration of butene and the catalytic reaction of 1-butene to 2-butene is a synergetic and elementary process. The apparent activation energy of isomerization is about 197 kJ/mol.     

DFT Studies on the Isomerization of Butene Double Bond Catalyzed by 1-Butyl-3-methyl-imidazolium in Ionic Liquid

1 INTRODUCTION Butene and its isomers are important petroleum raw materials. Isomerization reaction of butene plays a key role in the course of C4 alkylation and its reaction mechanism has captured the attention of chemists all along[1, 2]. As a green so…     

DFT法研究离子液中EMIM+催化丁烯双键异构反应机理(II)

用密度泛函理论(DFT)在B3LYP/6-31G(d, p)的计算水平上研究了离子液中1-乙基-3-甲基咪唑阳离子(EMIM⁺)的4-H和5-H原子催化丁烯双键异构反应的可能途径，优化了反应体系的平衡态和过渡态的几何构型，分析了反应过程中键参数的变化，通过振动分析对平衡态和过渡态进行了验证. 计算结果表明, 离子液中的EMIM⁺首先通过4-H和5-H原子吸附丁烯, 进而催化丁烯的双键异构反应, EMIM⁺的4-H和5-H催化1-丁烯异构为2-丁烯的正反应活化能分别为204.2和207.3 kJ•mol^-1，逆反应活化能约为220.9和223.8 kJ•mol^-1， 反应为基元反应.     

分子筛催化cis-2-丁烯的双键异构反应机理的DFT研究

基于含有两个Si和一个Al的分子筛3T簇模型, 利用密度泛函方法(DFT)研究了分子筛催化1-丁烯双键异构为cis-2-丁烯的反应机理. 在B3LYP/6-31G(d,p)计算水平上对反应各驻点进行了全优化, 并计算了反应的活化能. 研究发现, 分子筛上的酸性OH基团首先通过物理吸附靠近1-丁烯的双键, 形成了π配位复合物后, 丁烯双键的端基C原子逐渐抽取这个质子, 同时相邻酸性位的一个O原子也抽取丁烯碳链上的一个H原子, 形成吸附态的cis-2-丁烯, 最后通过脱附形成产物, 使分子筛复原, 反应按照协同反应机理发生. 计算得到的表观活化能是55.9 kJ/mol, 与实验结果接近.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

DFT法研究分子筛催化trans-2-丁烯的双键异构

利用一个3T簇模型模拟分子筛催化剂的酸性位, 采用密度泛函理论(DFT)的 B3LYP/6-31G(d, p)方法, 研究了分子筛催化1-丁烯双键异构为trans-2-丁烯的反应机理. 对反应各驻点进行了全优化, 经过零点能校正后, 得到了反应的活化能. 研究表明, 反应分三步进行：物理吸附→化学反应→物理脱附. 分子筛的酸性位OH基团首先吸附1-丁烯的双键形成了π配位复合物, 然后按协同反应机理发生双键异构反应, 生成吸附态的trans-2-丁烯, 最后脱附成产物. 计算得到的表观活化能为57.1 kJ•mol^-1, 与实验结果接近.     

CH3CH2S自由基H迁移异构化及裂解反应的理论研究

采用密度泛函方法(B3LYP)在6-311+G(d,p)基组水平上研究了CH3CH2S自由基H迁移异构化以及裂解反应的微观动力学机理. 在QCISD(T)/6-311++G(d,p)//B3LYP/6-311+G(d,p)+ZPE水平上进行了单点能校正. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了在200~2000 K温度区间内的速率常数kTST和kCVT, 同时获得了经小曲率隧道效应模型(SCT)校正后的速率常数kCVT/SCT. 研究结果表明, CH3CH2S自由基1,2-H迁移、1,3-H迁移、C—C键断裂和β-C—H键断裂反应的势垒ΔE≠分别为149.74, 144.34, 168.79和198.29 kJ/mol. 当温度低于800 K时, 主要发生1,2-H迁移反应, 高于1800 K时, 主要表现为C—C键断裂反应, 在1300—1800 K范围内, 1,3-H迁移反应是优势通道, 在计算的整个温度段内, β-C—H键断裂反应可以忽略.     

Structural elucidation of nitro-substituted five-membered aromatic heterocycles utilizing GIAO DFT calculations

The GIAO (Gauge Including Atomic Orbitals) DFT (Density Functional Theory) method is applied at the B3LYP/6-31+G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311+G (2d,p)//B3LYP/6-31+G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-311++G(d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants for 25 nitro-substituted five-membered heterocycles. Difference (1D NOE) spectra in combination with long-range gHMBC experiments were used as tools for the structural elucidation of nitro-substituted five-membered heterocycles. The assigned NMR data (chemical shifts and coupling constants) for all compounds were found to be in good agreement with theoretical calculations using the GIAO DFT method. The magnitudes of one-bond (1JCH) and long-range (nJCH, n>1) coupling constants were utilized for unambiguous differentiation between regioisomers of nitro-substituted five-membered heterocycles.     

3-羰基吡唑质子转移过程的理论研究

采用密度泛函B3LYP／6—311G^＊＊方法,对3-羰基吡唑几何构型进行了全自由度优化,获得了它们的几何结构和电子结构．计算并考察了3-羰基吡唑的两种构象即syn和anti构象的稳定性以及3-羰基吡唑进行结构互变的质子转移过程的四种可能途径：（a）分子内质子转移;（b）水助质子转移;（C）同种二聚体双质子转移;（d）异种二聚体双质子转移．计算结果表明3-羰基吡唑的syn构象中N2-H型的稳定性大于N1-H型,进行质子转移时途径（C）所需要的活化能最小（52．78kJ／mol）,途径（a）所需要的活化能最大（200．59kJ／mol）;3,羰基吡唑的。anti构象中N1-H型的稳定性大于N2-H型,进行质子转移时途径（d）所需要的活化能最小（61．09kJ／mol）,途径（a）所需要的活化能最大（204．15kJ／mol）．     

Synthesis of new cyano-substituted analogues of Tröger's bases from bromo-derivatives. A stereochemical dependence of long-range (nJHH,n = 4, 5, and 6) proton–proton and proton–carbon (nJCH,n = 1, 2, 3, 4, and 5) coupling constants of these compounds

A free-catalyst microwave-assisted cyanation of brominated Tröger's base derivatives ( 2a - f ) is reported. The procedure is simple, efficient, and clean affording the nitrile compounds ( 3a - e, I ) in very good yields. Complete assignment of ¹H and ¹³C chemical shifts of 2a - f, I and 3a - d, I was achieved using gradient selected 1D nuclear magnetic resonance (NMR) techniques (1D zTOCSY, PSYCHE, DPFGSE NOE, and DEPT), homonuclear 2D NMR techniques (gCOSY and zTOCSY), and heteronuclear 2D NMR techniques (gHSQCAD/or pure-shift gHSQCAD, gHMBCAD, bsHSQCNOESY, and gHSQCAD-TOCSY) with adiabatic pulses. Determination of the long-range proton–proton coupling constants ⁿJ_HH (n = 4, 5, 6) was accomplished by simultaneous irradiation of two protons at appropriate power levels. In turn, determined coupling constants were tested by an iterative simulation program by calculating the ¹H NMR spectrum and comparing it to the experimental spectrum. The excitation-sculptured indirect-detection experiment (EXSIDE) and ¹H-¹⁵N CIGARAD-HMBC (constant time inverse-detection gradient accordion rescaled heteronuclear multiple bond correlation) were applied for determination of long-range carbon–proton coupling constants ⁿJ_CH (n = 2, 3, and 4) and for assignment of ¹⁵N chemical shift at natural abundance, respectively. DFT/B3LYP optimization studies were performed in order to determine the geometry of 2c using 6-31G(d,p), 6-311G(d,p), and 6–311 + G(d,p) basis sets. For calculation of ¹H and ¹³C chemical shifts, ⁿJ_HH (n = 2, 3, 4, 5, and 6), and ⁿJ_CH (n = 1, 2, 3, and 4) coupling constants, the GIAO method was employed at the B3LYP/6-31G(d,p), B3LYP/6-31+G(d,p), B3LYP/6-311+G(d,p), B3LYP/6-311++G(2d,2p), B3LYP/cc-pVTZ), and B3LYP/aug-cc-pVTZ) levels of theory. For the first time, a stereochemical dependence magnitude of the long-range ⁿJ_HH (n = 4, 5, and 6) and ⁿJ_CH (n = 1, 2, 3, 4, and 5) have been found in bromo-substituted analogues of Tröger's bases.     

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

用密度函数理论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计算值是至今与实验值符合得最好的理论计算结果。     

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     

Predicting a New Δ-Proton Sponge-Base of 4,12-Dihydrogen-4,8,12-triazatriangulene through Proton Affinity,Aromatic Stabilization Energy,and Aromatic Magnetism

Herein, we report designing a new Δ (delta-shaped) proton sponge base of 4,12-dihydrogen-4,8,12-triazatriangulene (compound 1 ) and calculating its proton affinity (PA), aromatic stabilization, natural bond orbital (NBO), electron density ρ(r), Laplacian of electron density ∇²ρ(r), (2D-3D) multidimensional off-nucleus magnetic shielding (σ_zz(r) and σ_iso(r)), and scanning nucleus-independent chemical shift (NICS_zz and NICS). Density functional theory (DFT) at B3LYP/6-311+G(d,p), ωB97XD/6-311+G(d,p), and PW91/def2TZVP were used to compute the magnetic shielding variables. In addition, relevant bases like pyridine, quinoline, and acridine were also studied and compared. The protonation of compound 1 yields a highly symmetric carbocation of three Hückel benzenic rings. Comparing our findings of the studied molecules showed that compound 1 precedes others in PA, aromatic isomerization stabilization energy, and basicity. Therefore, the basicity may be enhanced when a conjugate acid gains higher aromatic features than its unprotonated base. Both multidimensional σ_zz(r) and σ_iso(r) off-nucleus magnetic shieldings outperformed electron-based techniques and can visually monitor changes in aromaticity that occur by protonation. The B3LYP/6-311+G(d,p), ωB97XD/6-311+G(d,p), and PW91/def2TZVP levels showed no significant differences in detailing isochemical shielding surfaces.     

3-卤代吡唑质子转移过程的理论研究

采用密度泛函B3LYP/6-311G**方法,对3-卤(-F、-Cl、-Br)代吡唑几何构型进行了全自由度优化,获得了它们的几何结构和电子结构。计算结果显示,N1-H型的稳定性大于N2-H型。计算并考察了3-卤代吡唑进行结构互变的质子转移过程的四种可能途径:(a)分子内质子转移;(b)水助质子转移;(c)同种二聚体双质子转移;(d)异种二聚体双质子转移。计算结果表明(以3-氟代吡唑为例),途径d所需要的活化能最小(54.89 kJ/mol),而途径a所需要的活化能最大(198.83kJ/mol),途径b和c的活化能居中间分别为(104.05 kJ/mol和69.05 kJ/mol)。研究还表明氢键在降低活化能方面起着重要的作用,卤素(-F、-Cl、-Br)对活化能的影响不大。     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

Ab initio Study on Structures and Isomerization of Magnesium Fluorosilylenoid H2SiFMgF

The structures and isomerization of magnesium fluorosilylenoid H2SiFMgF were investigated by ab initio molecular orbital theory. Four equilibrium structures and three isomeric transition states were located and fully optimized at the B3LYP/6-31G(d,p) and G3MP2B3 levels, respectively. Based on the B3LYP/6-31G(d,p) optimized geometries, harmonic frequencies of various structures were obtained and 29Si chemical shifts were calculated. The solvent effects were investigated by means of the polarizable continuum model using THF as a solvent at B3LYP/6-31G(d,p) level. Isomerization paths for isomers were confirmed by intrinsic reaction coordinate calculations. The calculated results show that tetrahedral structure has the lowest energy and is the most stable; tetrahedral, three-membered ring, and p-complex structures are suggested to be the experimentally detectable ones; and σ-complex structure has the highest energy and will not exist.     

Proton affinities of ketones, vicinal diketones and α-keto esters: a computational study

The proton affinities of seven different ketones, vicinal diketones, and α-keto esters (acetophenone, 2,2,2-trifluoroacetophenone, 2,3-butanedione, 1-phenyl-1,2-propanedione, methyl pyruvate, ethyl benzoylformate, and ketopantolactone) have been evaluated theoretically using the conventional ab initio HF and several post-HF methods (MP2, MP4, CCSD), density functional methods with the B3LYP hybrid functional, as well as some ab initio model chemistries [CBS-4M, G2(MP2), and G3(MP2)//B3LYP]. The chemical compounds studied are frequently used substrates in the asymmetric hydrogenation over chirally modified platinum catalysts where the protonation properties of the chiral modifier and the substrates are of great interest. In most cases, the proton affinities (PAs) evaluated with the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) methods are in good agreement with the existing experimental ones. However, the previously reported PA of 2,3-butanedione seems to be too high by 10-15 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP model chemistries predict proton affinities that are systematically higher and lower than the experimental PAs, respectively. If proton affinities are evaluated as the average of the PAs calculated with these two theoretical methods a very good agreement with the experimental results is obtained. The mean absolute deviation (MAD) from experiment of this combination method for the PAs of 13 test molecules is 4.0 kJ mol⁻¹. For 9 molecules composed only of first-row atoms the MAD is 2.5 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP methods provide significant savings in computational time and disk space compared to the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) models. Therefore, it is suggested that if no experimental or highly accurate theoretical data is available (due to computational cost), the proton affinities of similar compounds as investigated in this paper, can be evaluated with the combination method. For the studied molecules, this method gives the following PAs (in kJ mol⁻¹): 788 (2,3-butanedione, exptl 802); 798 (2,2,2-trifluoroacetophenone, exptl 799); 811 (ketopantolactone); 813 (methyl pyruvate); 825 (1-phenyl-1,2-propanedione); 862 (acetophenone, exptl 861); 865 (ethyl benzoylformate).     

2-(3-巯基-2-吡啶基)苯并咪唑分子内质子转移及溶剂化效应

在B3LYP/6-31G(d,p)和TD B3LYP/6-31++G(d,p)//CIS/6-31G(d,p)水平上,研究了2-(3-巯基-2-吡啶基)苯并咪唑(MPyBI)在气相和七种溶剂(环己烷、苯、三氯甲烷、乙醇、乙腈、二甲亚砜和水)中基态和激发态的分子内质子转移(GSIPT和ESIPT)过程.在基态势能面的研究中发现,该化合物存在分子内双质子转移,其中分步的双质子转移在动力学上具有优势.同时对激发态质子转移势能面及激发态转移过程中的光物理现象进行了研究,结果表明该化合物存在快速的无能垒的激发态分子内质子转移,随着溶剂极性的增强,可以降低基态过渡态的能垒,改变硫醇式与硫酮式互变异构体的比例,从而灵敏地控制荧光的强度.     

The rate enhancement for prolyl cis-to-trans isomerization of cyclic CPFC peptide is caused by an increase in the vibrational entropy of the transition state

The conformational preferences and prolyl cis-trans isomerization of oxidized and reduced Ac-Cys-Pro-Phe-Cys-NH2 (CPFC peptides) have been carried out using the ab initio HF/6-31+G(d) and hybrid density functional B3LYP/6-311++G(d,p) levels of theory. The most preferred conformations of oxidized and reduced CPFC peptides with the trans prolyl peptide bond have a type-I beta-turn for the Pro-Phe sequence in common. In particular, the transition states for both forms are stabilized by the intramolecular hydrogen bonds between the prolyl nitrogen and the N-H group of the Phe3 residue. The rotational barrier DeltaGct to the cis-to-trans isomerization for the oxidized CPFC peptide is calculated to be 19.37 kcal/mol at the B3LYP/6-311++G(d,p)//HF/6-31+G(d) level of theory, which is lower by 0.88 kcal/mol than that of the reduced CPFC peptide. This may indicate that the rate constant kc-->t of the prolyl cis-to-trans isomerization for the oxidized form is about 4 times larger than that of the reduced form, which is reasonably consistent with the value deduced from NMR experiments. In particular, the increase in vibrational entropy for the transition state of the oxidized form over that of the reduced form contributes to enhance the rate constant for the prolyl cis-to-trans isomerization of the oxidized form.     

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.