Ab initio Mechanism Study on the Reaction of Chlorine Atom with Formic Acid

The potential energy surface(PES) for the reaction of Cl atom with HCOOH is predicted using ab initio molecular orbital calculation methods at UQCIDS(T,full)6-311 G(3df,2p)//UMP2(full)/6-311 G(d,P) level of theory with zero-point vibrational energy (ZPVE) correction.The calculated results show that the reaction mechanism of Cl atom with formic acid is a C-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom with a 3.73kJ/mol reaction barrier height,leading to the formation of cis-HOCO radical which will reacts with Cl atom or other molecules in such a reaction system.Because the reaction barrier height of O-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom which leads to the formation of HCO2 radical is 67.95kJ/mol,it is a secondary reaction channel in experiment,This is in good agreement with the prediction based on the previous experiments.     

An optimum strategy for solution chemistry using semiempirical molecular orbital method: Importance of description of charge distribution

For the purpose to execute direct dynamics calculation in solution chemistry, we propose an optimum strategy for solution chemistry using semiempirical molecular orbital (MO) method with neglect of diatomic differential overlap (NDDO) approximation with specific solution reaction parameters (SSRP), i.e., the NDDO‐SSRP method. In this strategy, the empirical parameters of the semi‐empirical MO method were optimized individually for target molecule or ion by reference to the ab initio MO calculation data for many configurations on the potential energy surface near the reaction path. For demonstration, the NDDO‐SSRP method was applied to two molecules and two ions (OH^?, H₂O, NH₃, NH₄⁺) at their equilibrium states in aqueous solution, respectively. Accordingly, it was verified that both the potential energy surface and the charge distribution of these solutes in aqueous solution are dramatically improved to reproduce themselves accurately at ab initio MO calculation level. In conclusion, it is expected that the NDDO‐SSRP method should become quite useful for dynamic and statistical applications to chemical reaction systems in solution. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Mechanism of cysteine oxidation by a hydroxyl radical: a theoretical study.

Cysteine oxidation by HO(.) was studied at a high level of ab initio theory in both gas phase and aqueous solution. Potential energy surface scans in the gas phase performed for the model system methanethiol+HO(.) indicate that the reactants can form two intermediate states: a sulfur-oxygen adduct and a hydrogen bound reactant complex. However these states appear to play a minor role in the reaction mechanism as long as they are fast dissociating states. Thus the main reaction channel predicted at the QCISD(T)/6-311+G(2df,2pd) level of theory is the direct hydrogen atom abstraction. The reaction mechanism is not perturbed by solvation which was found to induce only small variations in the Gibbs free energy of different reactant configurations. The larger size reactant system cysteine+HO* was treated by the integrated molecular orbital+molecular orbital (IMOMO) hybrid method mixing the QCISD(T)/6-311+G(2df,2pd) and the UMP2/6-311+G(d,p) levels of theory. The calculated potential energy, enthalpy, and Gibbs free energy barriers are slightly different from those of methanethiol. The method gave a rate constant for cysteine oxidation in aqueous solution, k=2.4 x 10(9) mol(-1) dm(3) s(-1), which is in good agreement with the experimental rate constant. Further analysis showed that the reaction is not very sensitive to hydrogen bonding and electrical polarity of the molecular environment.     

The self‐consistent charge density functional tight‐binding theory study of carbon adatoms using tuned Hubbard U parameters

The self‐consistent charge density functional tight‐binding (DFTB) theory is a useful tool for realizing the electronic structures of large molecular complex systems. In this study, the electronic structure of C₆₁ formed by fullerene C₆₀ with a carbon adatom is analyzed, using the fully localized limit and pseudo self‐interaction correction methods of DFTB to adjust the Hubbard U parameter (DFTB + U). The results show that both the methods used to adjust U can significantly reduce the molecular orbital energy of occupied states localized on the defect carbon atom and improve the gap between highest occupied molecular orbital(HOMO) and lowest unoccupied molecular orbital(LUMO) of C₆₁. This work will provide a methodological reference point for future DFTB calculations of the electronic structures of carbon materials.     

Activation Energy and Transition State Determination of the Olefin Insertion Process of Metallocene Catalysts Using a Semiempirical Molecular Orbital Calculation

The transition state of the olefin insertion process of metallocene catalysts can be determined by adopting the semiempirical PM3 model. In computational chemistry, the computational methods most employed are the ab initio method and density functional theory, which are very time consuming. The semiempirical molecular orbital method requires much less computational resources than the above methods. However, the accuracy and reliability of the semiempirical molecular orbital method remains to be determined. The PM3 model is the most recently developed the semiempirical molecular orbital method and can also be applied to transition metal calculations. This study is intended to investigate the reliability of computational results determined using semiempirical PM3 model on metallocene catalysts through comparison with published results on the density functional theory (DFT). The saddle point finding procedure is adopted to find the transition state of the ethylene insertion process of metallocene catalysts. Results on the geometry and energy trends of the ethylene insertion process of metallocene catalysts determined using the PM3 model are in good agreement with the DFT results. In addition, the saddle point of the potential energy surface of ethylene insertion is verified in accordance with the eigenvalue of the vibrational frequency spectrum. Correct eigenvalues indicate that the correct saddle point of the potential energy surface of ethylene insertion has been successfully located. Hence, the eigenvalue of the vibrational frequency spectrum is a valuable reference in terms of saddle point justification. Computational results and vibrational frequency spectrum analysis demonstrate that the PM3 model can be used to locate the correct saddle point of the potential energy surface. The results obtained using the PM3 model confirm that the eigenvalue of the transition state lies nearly on the vibrational frequency spectrum. The eigenvalues are also analyzed, providing a valuable reference for further studies of the transition state of olefin insertion of metallocene catalysts. The activation energies for the olefin insertion reaction are also studied for evaluation of the catalyst.     

Density functional study of structural and electronic properties of AlnN (1 ≤ n ≤ 12) clusters

Low‐lying equilibrium geometric structures of Al_nN (n = 1–12) clusters obtained by an all‐electron linear combination of atomic orbital approach, within spin‐polarized density functional theory, are reported. The binding energy, dissociation energy, and stability of these clusters are studied within the local spin density approximation (LSDA) and the three‐parameter hybrid generalized gradient approximation (GGA) due to Becke–Lee–Yang–Parr (B3LYP). Ionization potentials, electron affinities, hardness, and static dipole polarizabilities are calculated for the ground‐state structures within the GGA. It is observed that symmetric structures with the nitrogen atom occupying the internal position are lowest‐energy geometries. Generalized gradient approximation extends bond lengths as compared with the LSDA lengths. The odd–even oscillations in the dissociation energy, the second differences in energy, the highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) gaps, the ionization potential, the electron affinity, and the hardness are more pronounced within the GGA. The stability analysis based on the energies clearly shows the Al₇N cluster to be endowed with special stability. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

The topology of energy hypersurfaces III. The fundamental group of reaction mechanisms on potential energy hypersurfaces

The family of all possible reaction mechanisms on a potential surface has an algebraic structure with potential applications in quantum chemical molecular design and synthesis planning.Transformation properties and equivalence relations of reaction paths on potential energy hypersurfaces lead to a topological definition of reaction mechanisms. The family of all fundamental reaction mechanisms on the hypersurface has a group structure,the fundamental group of an appropriately defined topological space. Isomorphism and homomorphism relations between fundamental groups of reaction mechanisms are used to characterize the chemically important topological properties of various subsets of a hypersurface, or those of different excited state hypersurfaces.     

Realistic Energy Surfaces for Real‐World Systems: An IMOMO CCSD(T):DFT Scheme for Rhodium‐Catalyzed Hydroformylation with the 6‐DPPon Ligand

The hydroformylation of terminal alkenes is one of the most important homogeneously catalyzed processes in industry, and the atomistic understanding of this reaction has attracted enormous interest in the past. Herein, the whole catalytic cycle for rhodium‐catalyzed hydroformylation with the 6‐diphenylphosphinopyridine‐(2H)‐1‐one (6‐DPPon) ligand 1 was studied. This catalytic transformation is challenging to describe computationally, since two requirements must be met: 1) changes in the hydrogen‐bond network must be modeled accurately and 2) bond‐formation/bond‐breaking processes in the coordination sphere of the rhodium center must be calculated accurately. Depending on the functionals used (BP86, B3LYP), the results were found to differ strongly. Therefore, the complete cycle was calculated by using highly accurate CCSD(T) computations for a PH₃ model ligand. By applying an integrated molecular orbital plus molecular orbital (IMOMO) method consisting of CCSD(T) as high level and DFT as low‐level method, excellent agreement between the two functionals was achieved. To further test the reliability of the calculations, the energetic‐span model was used to compare experimentally derived and computed activation barriers. The accuracy of the new IMOMO method apparently makes it possible to predict the catalytic potential of real‐world systems.     

On the Electrochemical Deposition and Dissolution of Divalent Metal Ions

The deposition of Cu²⁺ and Zn²⁺ from aqueous solution has been investigated by a combination of classical molecular dynamics, density functional theory, and a theory developed by the authors. For both cases, the reaction proceeds through two one‐electron steps. The monovalent ions can get close to the electrode surface without losing hydration energy, while the divalent ions, which have a stronger solvation sheath, cannot. The 4s orbital of Cu interacts strongly with the sp band and more weakly with the d band of the copper surface, while the Zn 4s orbital couples only to the sp band of Zn. At the equilibrium potential for the overall reaction, the energy of the intermediate Cu⁺ ion is only a little higher than that of the divalent ion, so that the first electron transfer can occur in an outer‐sphere mode. In contrast, the energy of the Zn⁺ ion lies too high for a simple outer‐sphere reaction to be favorable; in accord with experimental data this suggests that this step is affected by anions.     

Efficiency of simulated annealing for peptides with increasing geometrical restrictions

Simulated annealing (SA) is a popular global minimizer that can conveniently be applied to complex macromolecular systems. Thus, a molecular dynamics or a Monte Carlo simulation starts at high temperature, which is decreased gradually, and the system is expected to reach the low-energy region on the potential energy surface of the molecule. However, in many cases this process is not efficient. Alternatively, the low-energy region can be reached more effectively by minimizing the energy of selected molecular structures generated along the simulation pathway. The efficiency of SA to locate energy-minimized structures within 5 kcal/mol above the global energy minimum is studied as applied to three peptide models with increasing geometrical restrictions: (1) The linear pentapeptide Leu-enkephalin described by the ECEPP potential, (2) a cyclic hexapeptide described by the GROMOS force field energy E_GRO alone, and (3) the same cyclic peptide with E_GRO combined with a restraining potential based on 31 proton–proton restraints obtained from nuclear magnetic resonance (NMR) experiments. The efficiency of SA is compared to that of the Monte Carlo minimization (MCM) method of Li and Scheraga, and to our local torsional deformations (LTD) method for the conformational search of cyclic molecules. The results for the linear peptide show that SA provides a relatively weak guidance towards the most stable energy region; as expected, this guidance increases for the cyclic peptide and the cyclic peptide with NMR restraints. However, in general, MCM and LTD are significantly more efficient than SA as generators of low-energy minimized structures. This suggests that LTD might provide a better search tool than SA in structure determination of protein regions for which a relatively small number of restraints are provided by NMR. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1659–1670, 1999     

DFT investigation of the mechanism of CH2CO + O(3P) reaction

A detailed theoretical survey of the potential energy surface (PES) for the CH₂CO + O(³P) reaction is carried out at the QCISD(T)/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) level. The geometries, vibrational frequencies, and energies of all stationary points involved in the reaction are calculated at the B3LYP/6‐311+G(d,p) level. More accurate energy information is provided by single‐point calculations at the QCISD(T)/6‐311+G(3df,2p) level. Relationships of the reactants, transition states, intermediates, and products are confirmed by the intrinsic reaction coordinate (IRC) calculations. The results suggest that P1(CH₂+CO₂) is the most important product. This study presents highlights of the mechanism of the title reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Multicomponent Molecular Orbital–Climbing Image–Nudged Elastic Band Method to Analyze Chemical Reactions Including Nuclear Quantum Effect

To analyze the H/D isotope effects on hydrogen transfer reactions in XHCHCHCHY?XCHCHCHYH (X, Y=O, NH, or CH₂) including the nuclear quantum effect of proton and deuteron, we propose a multicomponent molecular orbital‐climbing image‐nudged elastic band (MC_MO–CI–NEB) method. We obtain not only transition state structures but also minimum‐energy paths (MEPs) on the MC_MO effective potential energy surface by using MC_MO–CI–NEB method. We find that nuclear quantum effect affects not only stationary‐point geometries but also MEPs and electronic structures in the reactions. We clearly demonstrate the importance of including nuclear quantum effects for H/D isotope effect on rate constants (k_H/k_D).     

Reaction Mechanism of the Multi‐channel Decomposition Reactions of 1‐Pentenyl Free Radicals

The reactions of 1‐pentenyl decomposition system have been studied extensively at the B3LYP/6‐311++G?? ?? level with Gaussion 98 package. The potential energy surface with zero‐point energy correction was drawn. All reaction channels were fully investigated with the vibrational mode analysis, frontier orbital analysis and electron population analysis to confirm the transition states and reveal the reaction mechanism.     

Effect of a complex formation on the calculated low-pressure rate constant of a bimolecular gas-phase reaction governed by tunneling

Many important bimolecular hydrogen-transfer processes that take place in the atmosphere proceed via a potential energy minimum (hydrogen-bonded complex) that precedes along the minimum energy path the unique saddle point of the reaction, the one corresponding to the hydrogen transfer. It is clear that the one-step low-pressure rate constant of such a reaction does not depend on the existence of any complex along the minimum energy path below the reactant if the reaction takes place by thermal activation over a transition state that lies quite above the reactants (for instance 10 kcal/mol). However, we have quantitatively shown in this article that the scenario notoriously changes if the reaction involves significant tunneling. In this work, we have theoretically calculated the rate constants and their temperature dependence for the reaction HO+HOH→HOH+OH by means of a canonical variational transition state theory and a canonical unified statistical theory (when necessary). Multidimensional tunneling effects have been included with a semiclassical transmission coefficient. Two kinds of modified potential energy surfaces (PESs), obtained from an original ab initio potential energy surface, previously calculated by us, have been used. The Eckart-modified PESs serve to model the hydrogen-abstraction profiles with no complexes along the path, while the Gaussian-modified PESs model the energy profiles with two complexes along the path symmetrically distributed at each side of the abstraction saddle point. Our results show that the existence of those complexes reduces the thickness of the classically forbidden region for energies below the adiabatic barrier, and then tunneling is promoted and the reaction is accelerated. The effect of the complex formation in several kinetic magnitudes, as the Arrhenius parameters and the kinetic isotope effect has also been analyzed. © 1999 John Wiley & Sons, Inc. J Comput Chem 20: 1685–1692, 1999     

New hybrid method for reactive systems from integrating molecular orbital or molecular mechanics methods with analytical potential energy surfaces

A computational approach to calculating potential energy surfaces for reactive systems is presented and tested. This hybrid approach is based on integrated methods where calculations for a small model system are performed by using analytical potential energy surfaces, and for the real system by using molecular orbital or molecular mechanics methods. The method is tested on a hydrogen abstraction reaction by using the variational transition-state theory with multidimensional tunneling corrections. The agreement between the calculated and experimental information depends on the quality of the method chosen for the real system. When the real system is treated by accurate quantum mechanics methods, the rate constants are in excellent agreement with the experimental measurements over a wide temperature range. When the real system is treated by molecular mechanics methods, the results are still good, which is very encouraging since molecular mechanics itself is not at all capable of describing this reactive system. Since no experimental information or additional fits are required to apply this method, it can be used to improve the accuracy of molecular orbital methods or to extend the molecular mechanics method to treat any reactive system with the single constraint of the availability of an analytical potential energy surface that describes the model system.     

Ab initio study of C + H3+ reactions

The reaction C + H3+ --> CH(+) + H2 is frequently used in models of dense interstellar cloud chemistry with the assumption that it is fast, i.e. there are no potential energy barriers inhibiting it. Ab initio molecular orbital study of the triplet CH3+ potential energy surface (triplet because the reactant carbon atom is a ground state triplet) supports this hypothesis. The reaction product is 3 pi CH+; the reaction is to exothermic even though the product is not in its electronic ground state. No path has been found on the potential energy surface for C + H3+ --> CH2(+) + H reaction.     

Chemistry of novolac resins. X. Polymerization studies of HMTA and strategically synthesized model compounds

The reactions between hexamethylenetetramine (HMTA) and compounds which model novolac resins have been studied by ¹³C‐ and ¹⁵N‐NMR techniques. The dimer and tetramer compounds vary in molecular size and structure and react with HMTA to yield benzylamines and benzoxazine as the major initial‐formed intermediates and convert to methylene linked compounds at increased temperatures. The reaction of the compounds with only ortho reactive sites paralleled the 2,4‐xylenol–HMTA case reported by us previously; however, increasing molecular weight favored the formation of benzylamines and not benzoxazines. Those compounds with only para reactive sites paralleled the 2,6‐xylenol–HMTA case. The reactivity of the systems containing both ortho and para reactive sites depends on the ratio of ortho/para sites and various aspects such as the chemical structure and molecular weight of the compound, the HMTA level, and the melting point and pH of the system. These results parallel those obtained from novolac/HMTA systems. The xylenol/HMTA reactions formed similar products but showed quite different relative reaction rates by varying the HMTA ratio and structures of the materials. The importance of careful selection of model systems is also discussed. © 1999 John Wiley & Sons, Inc. J Polym Sci A: Polym Chem 37: 1347–1355, 1999     

The Rearrangement of 2,2‐Diphenyl‐1‐[(E)‐2‐phenylethenyl]cyclopropane to 3,4,4‐Triphenylcyclopent‐1‐ene: a DFT Analysis

A computational study on the rearrangement of 2,2‐diphenyl‐1‐[(E)‐2‐phenylethenyl]cyclopropane ( 1 ) is presented, using density functional theory (DFT), (U)B3LYP with the 6‐31G* basis set (DFT1) and (U)M05‐2X with the 6‐311+G** basis set (DFT2). In agreement with a biradical character of the transition structure (TS) or intermediate, the potential‐energy hypersurface is lowered by the influence of three conjugated Ph groups. Surprisingly, two conformations of the geminal diphenyl group (different twist angles) induce two different minimum‐energy pathways for the rearrangement. Independent of the functional used, the first hypersurface harbors true biradical intermediates, whereas the second energy surface is a flat, slightly ascending slope from the starting material to the TS. The functional (U)M05‐2X with the basis set 6‐311+G** provides realistic energies which seem to be close to experiment. The activation energy for racemization of enantiomers of 1 is lower than that of rearrangement by 2.5 kcal mol^?1, in agreement with experiment.     

O2 dissociative adsorption on the Cu‐, Ag‐, and W‐doped Al(111) surfaces from DFT computation

The O₂ adsorption and dissociation on M‐doped (M = Cu, Ag, W) Al(111) surface were studied by density functional theory. The adsorption energy of adsorbate, the average binding energy and surface energy of Al surface, and the doping energy of doping atom were calculated. All the doped atoms can be stably combined with Al atoms, while being slightly embedded in the surface to a certain depth. The TOP‐type surfaces are the most stable doped surfaces for O₂ adsorption, which is related to the orbital hybridization between the adsorbate and the surface atoms, the electronegativity, and the orbital energy level of the doping atoms. Moreover, the O atoms and doping atoms contribute significantly to the density of states (DOS), especially the O‐p orbital electrons and the d orbital electrons of doping atoms. The degree of O₂ dissociation is related to the doping atoms on Al surfaces, and the doping atoms actually resist the dissociation of O₂. W atoms have the best resistance effect on the O₂ dissociation as compared with Cu and Ag atoms, especially W‐1NN surface, which has both large barrier energy and reaction energy.