Rhodium catalyzed hydroformylation of olefins

DFT and CCSD(T) methods were used to examine 61 different rhodium catalysts for the hydroformylation of ethylene. The carbon monoxide (CO) stretching frequency was a key electronic parameter to understand the π-accepting nature of the ligand. Normally, π-accepting ligands lead to increased CO stretching frequencies and a reduction in CO dissociation energy. There was no relationship between CO dissociation energy and CO stretching frequency. However, a clear relationship exists between the ethylene insertion barrier (from the rhodium dicarbonyl hydride resting state) and the CO stretching frequency as stronger π-accepting ligands systematically led to a reduction in the barrier. Due to the multistep nature of the rate-limiting step, the overall barrier can be divided into the CO/ethylene equilibrium and an intrinsic ethylene insertion barrier and both are systematically reduced as the π-accepting nature of the ligand is increased. A comparison of the carbonylation transition state (TS) to the ethylene insertion TS allowed us to understand reversibility of olefin insertion. While the ethylene insertion TS systematically decreases with increasing CO stretching frequency, the carbonylation TS is relatively flat. The lines cross at 2156 cm⁻¹ implying a change in the rate-limiting step in this region given a standard set of process conditions. © 2018 Wiley Periodicals, Inc.     

De novo prediction of ground state multiplicity and structural isomerism for transition metal complexes

Prediction of the ground state geometries and multiplicities for 33 transition metal tetrachlorides has been carried out using two different levels of quantum mechanics: semiempirical and density functional theory. All data regarding geometry and spin state provided by both computational methods were compared with experimental data when available. The calculations were performed for all possible spin multiplicities. The most important geometries for coordination number four (tetrahedral, square-planar, dodecahedral, and disphenoidal), as well as less symmetric structural isomers, were evaluated. A match between both computational methods in terms of predicted ground state multiplicity and geometry was found for 26 species, which translated into almost 80% agreement. Even though the PM3(tm) geometry prediction protocol involved more steps for isolating a feasible global minimum, the aggregate of these calculations was still orders of magnitude faster than DFT calculations using extended basis sets. The calculations indicate that caution is needed in the application of the PM3(tm) method to very high-spin transition metal complexes, but point to the suitability of very rapid semiempirical methods for reliable prediction of structural and ‘spin’ isomers, and hence their use in an efficient de novo design protocol for transition metal complexes.     

氧杂环丁烷热解机理的量子化学研究

本文利用半经验分子轨道理论研究了氧杂环丁烷热解为甲醛和乙烯的反应机理计算是采用半经验方法AM1进行的, 各种驻点全部运用Berny梯度方法优化. 同时, 对过渡态的结构进行了振动分析的确证. 计算表明: 1)不存在协同的同面-同面反应途径的过渡态, 其驻点只是一个二级鞍点; 2) 协同的同面-异面反应途径需要经过一个能量很高的过渡态; 3)有利的反应途径是包含了双自由基中间体的分步过程。     

The mechanism of ethylene polymerization by nickel salicylaldiminato catalysts - Agostic interactions and their kinetic isotope effects

The ethylene polymerization reaction of a neutral nickel catalyst was studied by DFT calculations at the Becke3LYP/6-31G(d) level of theory. As in related cases a β-agostic bond stabilizes the nickel alkyl ground states. Transition states for the insertion of the olefin show a distinct α-agostic interaction, which has not been observed for late metal polymerization catalysts before. An ethylene-alkyl complex was identified as the resting state of the reaction. The overall barrier height of the reaction amounts to 17.54 kcal/mol, which slightly increases to 17.60 kcal/mol for the polymerization of deuterated ethylene. Therefore, a small positive kinetic isotope effect (k_H/k_D = 1.09) can be calculated, which is caused by the α-agostic interaction in the transition state. A comparison to other late metal based polymerization systems reveals that the ethylene coordination step of highly active catalysts is significantly lower in energy compared to catalysts which are only moderately active.     

Ab initio study on mechanisms of singlet oxygen reaction with ethylene

The reaction of singlet oxygen with ethylene has been studied at the ab initio level with both HF/3-21G and HF/6-31G basis sets, fully optimizing the geometries of the critical points. The transition state leading to the intermediate peroxirane from the initial reactants is found,Iocated at 81.71 kJ/mol above the dissociation limit. The vibrational analysis is done with two basis sets. From the results it can be seen that the frequencies have not been made an improvement obviously with augmented 6-31G compared to 3-21G basis; it follows that main reason for too high HF/3-21G frequency could mainly be the vibrational anharmonicity. The eigenvector corresponding to the single imaginary vibrational frequency is dominated by the larger O-O distance. The finding of the transition state confirms that the peroxirane minimum can be reachable passing through a peroxirane-like saddle point. In addition , the mechanisms of the reaction forming oxirane are discussed as well.     

C1桥联的茂金属化合物(Cp)CR1R2(Flu)MCl2的合成及其催化烯烃聚合活性

合成了6种单碳桥联的含芴(Flu)茂(Cp)基B族茂金属催化剂,研究了它们催化烯烃聚合的能力.通过IR,¹HNMR,EI-MS和元素分析对化合物进行了表征.用所合成的茂金属化合物与MAO所组成的催化体系对乙烯、丙烯的聚合进行了研究.其中金属为Ti的催化剂没有聚合活性或活性极低.金属为Zr的催化剂有一定的催化活性,用不同的催化剂得到的聚合物性质有一定的差异.     

A theoretical study of the comonomer effect in the ethylene polymerization with zirconocene catalytic systems

The PM3(tm) semiempirical method has been used to optimize the structures for the reactants and transition states of the first and second ethylene insertion processes into zirconocene catalytic systems. The results obtained for these reactions are compared with calculations published in the literature performed at different ab-initio theoretical levels. The agreement between our calculations and those reported in the literature is satisfactory. Taking advantage of the reduced computational effort required in semiempirical calculations two additional processes related with the so-called comonomer effect were also studied: ethylene/1-hexene copolymerization, and chain termination reaction, both in the homopolymerization and in copolymerization of ethylene with 1-hexene comonomer. The calculated activation energies support some experimental findings such as the higher polymerization activities in the presence of comonomers and also the molecular weight reduction of the copolymers due to the more favorable β-elimination reactions. © 1998 John Wiley & Sons, Inc. J Polym Sci A: Polym Chem 36: 1157–1167, 1998     

Theoretical study of C24N4 molecule

Both ab initio 6-31G, 3-21G and STO-3G basis sets and semiempirical PM3 and AM1 molecular orbital calculations are carried out on the C₂₄N₄ molecule of the T_d symmetry group. Results on the fully optimized structure which constrained T_d symmetry, molecular orbitals and vibrational frequency were obtained by both ab initio and semiempirical methods. The binding energy and various thermodynamic properties were also calculated via the PM3 and AM1 semiempirical methods. All the evidence of this work proves that the C₂₄N₄ molecule is stable and that its four six-membered rings with a remarkable delocalized C…C bond are similar to the related rings in the C₆₀ buckminsterfullerene structure.     

苯酚分子拉曼光谱的密度泛函理论研究

采用B3LYP杂化泛函,6-31++G(d,p)基组,对苯酚分子的结构进行了优化.通过频率计算,获得了苯酚分子的拉曼光谱,并利用势能函数分布(PED)对拉曼光谱进行了指认.绘制了苯酚分子的静电势分布图,讨论了苯酚分子发生化学反应的位置,同时计算了HOMO-LUMO的能级差.通过TD-DFT计算获得了苯酚分子的吸收光谱和激发态.     

Recent progress on olefin polymerization catalysts

Recent progress on metallocene catalysts is reviewed. This consists for the main part of research activities in metallocene catalysts and their polymerization performances (ethylene polymerization, propylene polymerization, styrene polymerization). In addi-tion, the computational design of metallocene catalysts and the commercial status of metallocene technologies in Japan are described.     

Semiempirical configuration interaction calculations for ru‐centered dyes*

Computational investigation of the photochemical properties of transition‐metal‐centered dyes typically involves optimization of the molecular structure followed by calculation of the UV/visible spectrum. At present, these steps are usually carried out using density functional theory (DFT) and time‐dependent DFT calculations. Recently, we demonstrated that semiempirical methods with appropriate parameterization could yield geometries that were in very good agreement with DFT calculations, allowing large sets of molecules to be screened quickly and efficiently. In this article, we modify a configuration interaction (CI) method based on a semiempirical PM6 Hamiltonian to determine the UV/visible absorption spectra of Ru‐centered complexes. Our modification to the CI method is based on a scaling of the two‐center, two‐electron Coulomb integrals. This modified, PM6‐based method shows a significantly better match to the experimental absorption spectra versus the default configuration interaction method (in MOPAC) on a training set of 13 molecules. In particular, the modified PM6 method blue‐shifts the location of the metal‐to‐ligand charge‐transfer (MLCT) peaks, in better agreement with experimental and DFT‐based computational results, correcting a significant deficiency of the unmodified method. Published 2018. This article is a U.S. Government work and is in the public domain in the USA     

A theoretical study of neighboring-group participation in thione-to-thiol rearrangement of xanthates. molecular orbital calculation using a conductor-like screening model (COSMO) approach.

The transition structures of the thione-to-thiol rearrangement of O-(2-dimethylaminoethyl and 2-methylthio-) S-methyl xanthates in several solvents were located by semiempirical molecular orbital method (PM3) using the conductor-like screening model (COSMO) approach. Each transition state transforms into the ion-pair intermediate with a three-membered ring structure (aziridinium or thiiranium), indicating that the rearrangement proceeds through an ionic intermediate with the anchimeric assistance of the neighboring group. The intermediary structures in gas phase are also analyzed by ab initio and density functional theory calculations.     

Theoretical study of photoinduced epoxidation of olefins catalyzed by ruthenium porphyrin

Epoxidation of olefin by [Ru(TMP)(CO)(O)](-) (TMP = tetramesitylporphine), which is a key step of the photocatalyzed epoxidation of olefin by [Ru(TMP)(CO)], is studied mainly with the density functional theory (DFT) method, where [Ru(Por)(CO)] is employed as a model complex (Por = unsubstituted porphyrin). The CASSCF method was also used to investigate the electronic structure of important species in the catalytic cycle. In all of the ruthenium porphyrin species involved in the catalytic cycle, the weight of the main configuration of the CASSCF wave function is larger than 85%, suggesting that the static correlation is not very large. Also, unrestricted-DFT-calculated natural orbitals are essentially the same as CASSCF-calculated ones, here. On the basis of these results, we employed the DFT method in this work. Present computational results show characteristic features of this reaction, as follows: (i) The epoxidation reaction occurs via carboradical-type transition state. Neither carbocation-type nor concerted oxene-insertion-type character is observed in the transition state. (ii) Electron and spin populations transfer from the olefin moiety to the porphyrin ring in the step of the C-O bond formation. (iii) Electron and spin populations of the olefin and porphyrin moieties considerably change around the transition state. (iv) The atomic and spin populations of Ru change little in the reaction, indicating that the Ru center keeps the +II oxidation state in the whole catalytic cycle. (v) The stability of the olefin adduct [Ru(Por)(CO)(O)(olefin)](-) considerably depends on the kind of olefin, such as ethylene, n-hexene, and styrene. In particular, styrene forms a stable olefin adduct. And, (vi) interestingly, the difference in the activation barrier among these olefins is small in the quantitative level (within 5 kcal/mol), indicating that this catalyst can be applied to various substrates. This is because the stabilities and electronic structures of both the olefin adduct and the transition state are similarly influenced by the substituent of olefin.     

Novel polystyrene‐supported zirconocene catalyst for olefin polymerization and its catalytic kinetics

Through the Diels–Alder reaction between cyclopentadiene groups attached to polystyrene in the presence of zirconocene, novel polystyrene‐supported metallocene catalysts were prepared. A novel method for immobilizing metallocene catalysts was investigated, and the resultant polystyrene‐supported metallocene for olefin polymerization was studied. The results of olefin polymerization showed that different crosslinking degrees of support in the catalyst system had significant effects on the catalytic behavior. The influence of the [Al]/[Zr] molar ratio and the temperature on the (co)polymerization activity was studied. When 1‐hexene and 1‐dodecene were used for copolymerization with ethylene, an obvious positive comonomer effect was observed. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 43: 2650–2656, 2005     

A sobering assessment of small‐molecule force field methods for low energy conformer predictions

We have carried out a large scale computational investigation to assess the utility of common small‐molecule force fields for computational screening of low energy conformers of typical organic molecules. Using statistical analyses on the energies and relative rankings of up to 250 diverse conformers of 700 different molecular structures, we find that energies from widely used classical force fields (MMFF94, UFF, and GAFF) show unconditionally poor energy and rank correlation with semiempirical (PM7) and Kohn–Sham density functional theory (DFT) energies calculated at PM7 and DFT optimized geometries. In contrast, semiempirical PM7 calculations show significantly better correlation with DFT calculations and generally better geometries. With these results, we make recommendations to more reliably carry out conformer screening.     

Semiempirical and DFT computations of the influence of Tb(III) dopant on unit cell dimensions of cerium(III) fluoride

Several computational methods, both semiempirical and ab initio, were used to study the influence of the amount of dopant on crystal cell dimensions of CeF₃ doped with Tb³⁺ ions (CeF₃:Tb³⁺). AM1, RM1, PM3, PM6, and PM7 semiempirical parameterization models were used, while the Sparkle model was used to represent the lanthanide cations in all cases. Ab initio calculations were performed by means of GGA+U/PBE projector augmented wave density functional theory. The computational results agree well with the experimental data. According to both computation and experiment, the crystal cell parameters undergo a linear decrease with increasing amount of the dopant. The computations performed using Sparkle/PM3 and DFT methods resulted in the best agreement with the experiment with the average deviation of about 1% in both cases. Typical Sparkle/PM3 computation on a 2×2×2 supercell of CeF3:Tb3+ lasted about two orders of magnitude shorter than the DFT computation concerning a unit cell of this material. © 2014 Wiley Periodicals, Inc.     

Mass-analyzed threshold ionization study of vinyl bromide cation in the first excited electronic state using vacuum-ultraviolet radiation generated by four-wave mixing in Hg

The vibrational spectrum of the vinyl bromide cation in the first excited electronic state A 2A' was obtained by one-photon mass-analyzed threshold ionization (MATI) spectroscopy. The use of an improved vacuum-ultraviolet radiation source based on four-wave sum frequency mixing in Hg resulted in excellent sensitivity for MATI signals. From the MATI spectrum, the ionization energy to the A 2A' state of the cation was determined to be 10.9150+/-0.0006 eV. Nearly complete vibrational assignments for the MATI peaks were possible by utilizing the vibrational frequencies and Franck-Condon factors calculated at the density-functional theory (DFT) and time-dependent DFT/B3LYP levels with the 6-311+G(df,p) basis set.     

双组分茂金属催化剂催化乙烯聚合的研究

选择能形成支链的不对称桥联茂金属化合物Me2 C[(Cp) (Ind) ]ZrCl2 和非桥联的不同结构的茂金属化合物二氯二 (烯基取代环戊二烯 )锆如 ( Cp) 2 ZrCl2 ,(Cp) 2 ZrCl2 和 (Cp) 2 ZrCl2 ,以MAO为助催化剂 ,分别组成三组双组分茂金属催化剂的催化体系 ,催化乙烯聚合 .结果表明 ,两类催化剂组成的双组分茂金属催化体系催化乙烯聚合能得到支化的宽分子量分布的聚乙烯 ;聚合温度和改变两种茂金属催化剂的摩尔比对催化活性和分子量有很大影响 .因此可以利用改变双组分茂金属催化剂的摩尔比例和聚合温度来调控聚合物的分子量和分子量分布 .改变两种茂金属催化剂的摩尔比和聚合温度也能使聚合物的结晶度发生改变     

De novo prediction of the ground state structure of transition metal complexes using semiempirical and ab initio quantum mechanics. Coordination isomerism

The ground state coordination isomers for 30 different trigonal bipyramidal transition metal complexes have been predicted using different levels of quantum mechanics: semiempirical (PM3(tm)), ab initio (MP2//HF), pure (BPW91) and hybrid (B3PW91) density functional theory (DFT) methods. For species where these methods failed to reproduce crystallographic data, hybrid quantum mechanics/molecular mechanics (QM/MM) methods were used to study more exact experimental models. Literature deficiencies regarding ground state multiplicity of these species were supplemented by spin predictions using previously tested PM3(tm) methods. Geometry optimization calculations were performed for each possible coordination isomer. The predicted ground state minima provided by the different methods are compared to each other and with crystallographic data. Pure DFT functionals outperformed hybrid functionals and MP2//HF. The very rapid PM3(tm) parameterization method provided accurate predictions in comparison to other levels of theory. An integrated MM/PM3(tm)/DFT de novo scheme accurately reproduced crystallographic data for species where the individual methods failed.     

双核茂金属催化剂催化聚合反应进展

综述了双核茂金属化合物的研究进展及其在烯烃和极性单体聚合方面的应用及双核茂金属化合物的合成及性能研究；并对双核催化剂的作用机理、催化剂结构与性能的关系也做了介绍。