Reactivity of isobutane on zeolites: a first principles study

In this work, ab initio and density functional theory methods are used to study isobutane protolytic cracking, primary hydrogen exchange, tertiary hydrogen exchange, and dehydrogenation reactions catalyzed by zeolites. The reactants, products, and transition-state structures are optimized at the B3LYP/6-31G* level, and the final energies are calculated using the CBS-QB3 composite energy method. The computed activation barriers are 52.3 kcal/mol for cracking, 29.4 kcal/mol for primary hydrogen exchange, 29.9 kcal/mol for tertiary hydrogen exchange, and 59.4 kcal/mol for dehydrogenation. The zeolite acidity effects on the reaction barriers are also investigated by changing the cluster terminal Si-H bond lengths. The analytical expressions between activation barriers and zeolite deprotonation energies for each reaction are proposed so that accurate activation barriers can be obtained when using different zeolites as catalysts.     

Reactivity of alkanes on zeolites: a computational study of propane conversion reactions

In this work, quantum chemical methods were used to study propane conversion reactions on zeolites; these reactions included protolytic cracking, primary hydrogen exchange, secondary hydrogen exchange, and dehydrogenation reactions. The reactants, products, and transition-state structures were optimized at the B3LYP/6-31G level and the energies were calculated with CBS-QB3, a complete basis set composite energy method. The computed activation barriers were 62.1 and 62.6 kcal/mol for protolytic cracking through two different transition states, 30.4 kcal/mol for primary hydrogen exchange, 29.8 kcal/mol for secondary hydrogen exchange, and 76.7 kcal/mol for dehydrogenation reactions. The effects of basis set for the geometry optimization and zeolite acidity on the reaction barriers were also investigated. Adding extra polarization and diffuse functions for the geometry optimization did not affect the activation barriers obtained with the composite energy method. The largest difference in calculated activation barriers is within 1 kcal/mol. Reaction activation barriers do change as zeolite acidity changes, however. Linear relationships were found between activation barriers and zeolite deprotonation energies. Analytical expressions for each reaction were proposed so that accurate activation barriers can be obtained when using different zeolites as catalysts, as long as the deprotonation energies are first acquired.     

H-ZSM-5分子筛上环己烯芳构化反应历程的理论研究

基于76T簇模型,采用量子力学和分子力学联合的ONIOM2(B3LYP/6-31G(d,p):UFF)方法研究了H-ZSM-5分子筛上环己烯芳构化反应历程.结果表明,环己烯首先吸附在分子筛酸性位上,与酸性质子共同脱除一个H2分子后,在分子筛骨架氧上生成烷氧配合物中间体;然后再脱质子得到环己二烯,同时酸性位复原;再经历脱氢和脱质子历程,最后得到产物苯,并吸附在复原的分子筛酸性位上.计算得到脱氢的活化能依次为279.64和260.21kJ/mol,脱质子的活化能依次为74.64和59.14kJ/mol.所有脱氢反应都是吸热过程,生成表面烷氧活性中间体,随后的脱质子反应能垒较低,而且是放热过程.此外,比较了环己烯在分子筛酸性位上的三个竞争反应,即脱氢、质子化和氢交换反应的活化能垒,证明环己烯优先发生脱氢反应.     

NC2S+离子的结构和稳定性的理论研究

在密度泛函和从头算理论水平下计算了单重态的NC2S+离子的结构、能量、光谱以及稳定性. 在B3LYP/6-311G(d)水平下, 得到8个异构体, 它们由15个过渡态相连接. 在CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE水平下, 得到能量最低的异构体是直线型的具有1Σ电子态的NCCS+(1)(0.0 kJ/mol), 其次是直线型的异构体CNCS+(2)(54.8 kJ/mol). 两个低能量的异构体1和2及另外一个高能量的直线型异构体CCNS+(3)(323.8 kJ/mol)都具有相当大的动力学稳定性, 这三个异构体在具备一定条件的实验室和星际条件下是可以进行观测的. 分析了这3个异构体的成键性质.     

亚烷基锡烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet alkylidenestannylene and ethylene has been investigated with MP2/3-21 G^* and B3LYP/3-21 G* methods, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD(T)//MP2/3-2 IG^* and CCSD(T)//B3LYP/3-21G^* methods, respectively. The results show that the dominant reaction pathway of the cycloaddition is that an intermediate (INT) is firstly formed between the two reactants through a barrier-free exothermic reaction of 39.7 kJ/mol, and the intermediate then isomerizes to a four-membered ring product (P2.1) via a transition state TS2.1 with a barrier of 66.8 kJ/mol.     

二氯亚甲基锗烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet dichloromethylene germylene and ethylene has been investigated with B3LYP/6-31G＊ method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD（T）//B3LYP/6-31G＊ method. On the basis of the surface energy profile obtained with CCSD（T）// B3LYP/6-31G＊ method for the cycloaddition reaction between singlet dichloromethylene germylene and ethylene, it can be predicted that the dominant reaction pathway is that an intermediate INT1 is firstly formed between the two reactants through a barrier-free exothermic reaction of 61.7 kJ/mol, and the intermediate INT1 then isomerizes to an active four-membered ring product P2.1 via a transition state TS2, an intermediate INT2 and a transition state TS2.1, in which energy barriers are 57.7 and 42.2 kJ/mol, respectively.     

F原子与瞬态自由基CH_2SH反应的理论研究

用量子化学从头算和密度泛函理论（DFT）对F原子与自由基CH_2SH在势能面上 的反应进行了研究。在B3LYP/6-311G水平上计算出了各物种的优化构型、振动频率 和零点振动能(ZPVE);各物种的总能量由B3LYP/6-311 + G(2df, pd)//B3LYP/6- 311G计算,另外对反应物和产物还计算了其G3能量。结果表明：首先F通过与C或S 结合的两种途径与CH_2SH相配位,再通过H（4）原子转移形成甲基,然后甲基再旋 转,甲基中H（4）原子最终与F结合,反应产物为HF和CH_2S。反应为放热反应,分 别为ΔH_r = -370.7 kJ/mol (DFT)和-396.94 kJ/mol (G3)。此外依据计算出的反 应热,可得自由基·CH_2SH的生成热Δ_fH°_(298.15) = 146.44 kJ/mol (DFT), 而Δ_fH°_0 = 167.36 kJ/mol (G3)。它们与以前的实验和理论值是一致的。     

HZSM-5 分子筛上乙烯芳构化过程中 C4 至 C6 中间体的反应机理

 应用量子力学和分子力学联合的 ONIOM2 (B3LYP/6-31G(d,p):UFF) 方法, 采用包含分子筛孔道结构的 78T 簇模型, 对 HZSM-5 分子筛上乙烯芳构化过程中 C4 至 C6 中间体的反应历程进行了研究, 探讨了分子筛的酸催化机理和择形催化作用. 结果表明, 作为乙烯二聚产物的表面正丁基烷氧络合物 (C4) 直接与乙烯作用得到正己基烷氧络合物 (C6), 在分子筛孔穴尺寸的限制下, 很难实现碳链的折叠环化. 按照间歇反应历程, 丁基烷氧络合物先发生 C–O 键断裂, 脱质子生成 1-丁烯, 然后在酸性位上再与乙烯加成, 在分子筛表面生成 3-甲基戊基烷氧络合物. 该烷氧络合物脱除质子给分子筛, 同时环化生成甲基环戊烷, 后者再与分子筛酸性质子共同脱除氢分子, 生成不稳定的碳正离子中间体, 然后重构成环己烷正离子. 丁基烷氧络合物脱质子的活化能为 158.42 kJ/mol; 1-丁烯与乙烯加成反应的活化能为 130.71 kJ/mol; 3-甲基戊基烷氧络合物脱氢环化生成甲基环戊烷的活化能为 122.06 kJ/mol. 由于孔穴的限域作用, 五员环的甲基环戊烷是重要的中间体.     

酸性分子筛催化乙烯二聚反应

应用密度泛函理论(DFT), 采用5T簇模型来模拟分子筛催化剂的酸性位, 在B3LYP/6-311+G(3df, 2p)的条件下通过理论计算研究了乙烯在酸性分子筛上的二聚反应. 对反应各驻点进行了全局优化, 经过零点能校正后, 计算得出乙烯二聚反应的活化能. 研究表明, 乙烯在分子筛上的二聚反应分三步进行: 单个乙烯分子化学吸附→第二个乙烯分子的物理吸附→两乙烯分子二聚反应. 乙烯化学吸附生成的烷氧化合物与物理吸附的乙烯分子发生二聚反应生成新的C—C键同时生成新的烷氧化合物. 计算得到的乙烯化学吸附和二聚反应的反应能垒分别为108和149 kJ·mol-1. 反应的逆过程也就是1-丁烯在酸性分子筛表面的1-丁基烷氧化合物发生β分裂反应, 计算所得相应的1-丁烯β分裂反应的能垒为217 kJ·mol-1, 远高于相应的乙烯二聚反应能垒. 此外还进一步研究了所用基组对计算结果的影响.     

Transition states for deactivation reactions in the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy‐mediated free‐radical polymerization of acrylonitrile

The geometries and energetics of transition states (TS) for radical deactivation reactions, including competitive combination and disproportionation reactions, have been studied for the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy (TEMPO)‐mediated free‐radical polymerization of acrylonitrile with quantum mechanical calculations at the DFT/UB3‐LYP/6‐311+G(3df,2p)//(U)AM1 level of theory (where DFT is density functional theory, AM1 is Austin model 1, and UAM1 is unrestricted Austin model 1). A method providing reasonable starting geometries for an effective search for TS between the TEMPO radical and 1‐cyanopropyl radical mimicking the growing polyacrylonitrile macroradical is shown. For the hydrogen atom abstraction reaction by the TEMPO radical from the 1‐cyanopropyl radical, practically one TS has been found, whereas for the combination reaction of the radicals, several TS have been found, mainly differing in out‐of‐plane angle α of the N? O bond in the TEMPO structure. α in the TS is correlated with the activation energy, ΔE, determined from the single‐point calculation at the DFT UB3‐LYP/6‐311+G(3df, 2p)//UAM1 level for the combination reaction of CH₃AN· with the TEMPO radical. The theoretical activation energy for the coupling reaction from DFT UB3‐LYP/6‐311+G(3df, 2p)//UAM1 calculations has been estimated to be 11.6 kcal mol^?1, that is, only about 4.5 times smaller than ΔE for the disproportionation reaction obtained with the DFT UB3‐LYP/6‐311+G(3df, 2p)//(U)AM1 approach. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 914–927, 2006     

Absorption spectra of azobenzenes simulated with time‐dependent density functional theory

Using time‐dependent density functional theory and the polarizable continuum model, we have simulated the absorption spectra of an extended series of azobenzene dyes. First, we have determined a theoretical level optimal for this important class of dyes, and it turned out that a C‐PCM‐CAM‐B3LYP/6‐311+G(d,p)//C‐PCM‐B3LYP/6‐311G(d,p) approach represents an effective compromise between chemical accuracy and computational cost. In a second stage, we have compared the theoretical and experimental transition energies for 46 n → π^☆ and 141 π → π^☆ excitations. For the full set, that spans over a 302–565 nm domain, we obtained a mean absolute deviation of 13 nm (0.10 eV) and a linear correlation coefficient of 0.95, illustrating the accuracy of our approach, though some significant outliers pertained. In a last step, the impact of several modifications, that is, trans/cis isomerization, variation of the acidity of the medium and azo/hydrazo tautomerism have been modeled with two functionals. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Ketene-ketenimine rearrangements in the gas phase and in polar media. 1,3-Migration intermediates and sequential transition states

[reaction: see text] Calculations of the activation barrier for the 1,3-shifts of substituents X in alpha-imidoylketenes 1 (HN=C(X)-CH=C=O), which interconverts them with alpha-oxoketenimines 3 (HN=C=CH-C(X)=O) via a four-membered cyclic transition state TS2 have been performed at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G* level. Substituents with accessible lone pairs have the lowest activation barriers for the 1,3-shift (halogens, OR, NR2). The corresponding activation barriers for the alpha-oxoketene-alpha-oxoketene rearrangement of 4 via TS5 are generally lower by 1-30 kJ/mol. A polar medium (acetonitrile, epsilon = 36.64) was simulated using the polarizable continuum (PCM) solvation model. The effect of the solvent field is a reduction of the activation barrier by an average of 12 kJ/mol. In the cases of 1,3-shifts of amino and dimethylamino groups, the stabilization of the transition state TS2 in a solvent field is so large that it becomes an intermediate, Int2, flanked by transition states (TS2' and TS2') that are due primarily to internal rotation of the amine functions, and secondarily to the 1,3-bonding interaction. In the case of the alpha-oxoketene-alpha-oxoketene rearrangement of 4, there is a corresponding intermediate Int5 for the 1,3-amine shift already in the gas phase.     

Transition structures and reaction barriers in the styrene–acrylonitrile copolymerization system according to quantum mechanical calculations

The geometries and electronic properties of substrates, transition structures (TS), and product radicals in modeled elementary propagation reactions were studied for the styrene–acrylonitrile monomer system by use of quantum‐mechanical calculations: (DFT/B3‐LYP/6–31G(d), ROMP2/6–311+G(3df,2p)//DFT/B3‐LYP/6–31G(d), and DFT/B3‐LYP/6–311+G(3df,2p)//DFT/B3‐LYP/6–31G(d)) and for some parameters, the high‐level composite method G3 (Gaussian‐3, G3/MP2). Activation enthalpies (ΔH^act) and reaction enthalpies (ΔH_r) for modeled propagation reactions at 298.15 K were evaluated. The enthalpy of activation energy (ΔH^act, kJ/mol) for the investigated elementary reactions rises for the B3‐LYP calculation in the following order: (CH₃A?+S) < (CH₃A?+A) < (CH₃S?+A) < (CH₃S?+S). For three propagation reactions, (CH₃A?+A), (CH₃A?+S), and (CH₃S?+A), correlation between reaction enthalpy and enthalpy of activation suggests weak or negligible polar effects reflecting the Evans–Polanyi relation. However, from the electron affinities and ionization energies values data, it is not excluded that at least for [CH₃A?+S[b]] and [CH₃S?+A[b]] reactions, nucleophilic and electrophilic polar effects, respectively, can also be expected. The dependencies between TS geometries, electronic parameters, and enthalpic effects suggest the presence of a steric factor in all TS, including its exceptionally high contribution to the activation enthalpy for the CH₃S?+S addition. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 43: 1827–1844, 2005     

Theoretical investigation on multinuclear NMR chemical shifts of some tris(trifluoromethyl)boron complexes

Tris(trifluoromethyl)boron complexes have unusual properties and may find applications in many fields of chemistry, biology, and physics. To gain insight into their NMR properties, the isotropic ¹¹B, ¹³C, and ¹⁹F NMR chemical shifts of a series of tris(trifluoromethyl)boron complexes were systematically studied using the gauge‐included atomic orbitals (GIAO) method at the levels of B3LYP/6‐31 + G(d,p)//B3LYP/6‐31G* and B3LYP/6‐311 + G(d,p)//B3LYP/6‐311 + G(d,p). Solvent effects were taken into account by polarizable continuum models (PCM). The calculated results were compared with the experimental values. The reason that the structurally inequivalent fluorine atoms in a specific species give a same chemical shift in experimental measurements is attributed to the fast rotation of CF₃ group around the B? C(F₃) bond because of the low energy barrier. The calculated ¹¹B, ¹³C(F₃), and ¹⁹F chemical shifts are in good agreement with the experimental measurements, while the deviations of calculated ¹³C(X, X = O, N) chemical shifts are slightly large. For the latter, the average absolute deviations of the results from B3LYP/6‐311 + G(d,p)//B3LYP/6‐311 + G(d,p) are smaller than those from B3LYP/6‐31 + G(d,p)//B3LYP/6‐31G*, and the inclusion of PCM reduces the deviation values. The calculated ¹⁹F and ¹¹B chemical shieldings of (CF₃)₃BCO are greatly dependent on the optimized structures, while the influence of structural parameters on the calculated ¹³C chemical shieldings is minor. Copyright © 2009 John Wiley & Sons, Ltd.     

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-     

分子筛催化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, 与实验结果接近.     

亚甲基硅烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet methylidenesilene and ethylene has been investigated with MP2/6-31G^＊ and B3LYP/6-31G^＊ methods, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies of the involved conformers were calculated by CCSD（T）//MP2/6-31G＊ and CCSD（T）//B3LYP/6-31 G＊ methods, respectively. The results show that the dominant reaction pathway of the cycloaddition reaction is that a complex intermediate is firstly formed between the two reactants through a barrier-free exothermic reaction of 13.3 kJ/mol, and the complex is then isomefized to a four-membered ring product P2,1 via a transition state TS2.1 with a barrier of 32.0 kJ/mol.     

An investigation on the reaction mechanism of the F2+Cl2→2ClF using the B3LYP method

The reaction mechanism of F₂+Cl₂→2ClF has been investigated with the density functional theory at the B3LYP/6‐311G* level. Six transition states have been found for the three possible reaction paths and verified by the normal mode vibrational and IRC analyses. Ab initio MP2/6‐311G* geometry optimizations and CCSD(T)/6‐311G(2df)//MP2/6‐311G* single‐point energy calculations have been performed for comparison. It is found that when the F₂ (or Cl₂) molecule decomposes into atoms first and then the F (or Cl) atom reacts with the molecule Cl₂ (or F₂) nearly along the molecular axis, the energy barrier is very low. The calculated energy barrier of F attacking Cl₂ is zero and that of Cl attacking F₂ is only 15.57 kJ?mol^?1 at the B3LYP level. However, the calculated dissociation energies of F₂ and Cl₂ are as high as 145.40 and 192.48 kJ?mol^?1, respectively. When the reaction proceeds through a bimolecular reaction mechanism, two four‐center transition states are obtained and the lower energy barrier is 218.69 kJ?mol^?1. Therefore, the title reaction F₂+Cl₂→2ClF is most probably initiated from the atomization of the F₂ molecule and terminated by the reaction of F attacking Cl₂ nearly along the Cl? Cl bond. MP2 calculations lead to the same conclusion, but the geometry of TS and the energy barrier are somewhat different. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem, 2002