Theoretical investigations on the mechanistic pathway of the thermal rearrangement of substituted N-acyl-2,2-dimethylaziridines

The mechanism of the thermal rearrangement of substituted N-acyl-2,2-dimethylaziridines 1 has been studied using quantum chemistry methods. Geometries of reactants, transition states and products have been optimized at the B3LYP/6-311++G(2d,2p) level. Relative energies for various stationary points have been determined and reaction identified by IRC calculations. The results show that thermal rearrangements occur in three ways. Firstly, the transition state TS ₁ in which a hydrogen atom of methyl groups migrates from primary carbon to oxygen of amid group to give the N-methallylamide 2. The second is via the transition state TS ₂ in which the attack of oxygen to the tertiary carbon yields the oxazoline 3. The third is via the transition state TS ₃ in which a hydrogen migrate from the secondary carbon to oxygen to give the vinylamide 4. In order to get insights into the factors determining the exact nature of its interactions with electrophiles, the application of reactivity parameters derived from density functional theory in a local sense, in particular the softness and Fukui function, to interpret and predict the mechanisms of the thermal decomposition of the N-acyl-2,2-dimethylaziridines 1, has been discussed.     

Quantum chemistry studies of the catalysis mechanism differences between the two isoforms of glutamic acid decarboxylase

The production of gamma-aminobutyric acid (GABA) is catalyzed by two isoforms of glutamic acid decarboxylase (GAD), using pyridoxal 5′-phosphate (PLP) as the cofactor. Between the two enzymes, GAD67 accounts for normal GABA requirement, while GAD65 stays inactive until emergent demand for GABA. Recent crystal structure findings revealed that the distinct conformation of a common catalytic loop of the enzymes may account for their different functions (Fenalti et al Nat Struct Mol Biol, 14:280-286, 2007). Enlightened by their inferences, we studied the underlying reaction mechanism of the two GAD isoforms using density functional theory (DFT). A rather complete reaction pathway is identified, including nine transition state (TS) structures and 14 intermediate (IM) structures. The rate limiting step occurs early during the reaction and involves a proton transfer. In the late stage, there are two pathways that involve C_4’ and C_α protonation by Tyr or Lys. Our calculations show that the reaction barriers corroborate the conjecture made by Fenalti et al.

Oxidation of pyruvic acid by flavins in the presence of amines

Pyruvic acid or other enolizable α-keto acids in the presence of primary or cyclic secondary amines in aprotic polar solvents efficiently reduce flavins and isoalloxazines to their 1,5-dihydro derivatives. The other product of this reaction is the diamide of citraconic acid (1). In the absence of flavin, the diamide of methylsuccinic acid (2) is obtained as the major product. The specificity for α-keto acid and amine, stoichiometric requirements, rates and quantities of CO₂ evolution, and the kinetics of the reaction were studied. On the basis of these data a mechanism is proposed that involves an enamine condensation followed by decarboxylation and then either reaction with flavin to give ultimately (1) and dihydroflavin, or reaction with an additional amine to give (2). The possible import of this kind of mechanism in biochemical systems is briefly discussed.     

Perspectives on the reaction force constant

A synchronous, concerted chemical process is rigorously divided by the reaction force F(R), the negative gradient of V(R), into “reactant” and “product” regions which are dominated by structural changes and an intervening “transition” region which is electronically intensive. The reaction force constant κ(R), the second derivative of V(R), is negative throughout the transition region, not just at the nominal transition state, at which κ(R) has a minimum. This is consistent with experimental evidence that there is a transition region, not simply a specific point. We show graphically that significant nonsynchronicity in the process is associated with the development of a maximum of κ(R) in the transition region, which increases as the process becomes more nonsynchronous. (We speculate that for a nonconcerted process this maximum is actually positive.) Thus, κ(R) can serve as an indicator of the level of nonsynchronicity.

Theoretical investigations on the synthesis mechanism of cyanuric acid from NH3 and CO2

In the synthesis of cyanuric acid from NH₃ and CO₂, urea and isocyanic acid OCNH are two pivotal intermediates. Based on density functional theory (DFT) calculations, the synthesis mechanism of cyanuric acid from NH₃ + CO₂ was investigated systematically. Urea can be synthesized from NH₃ and CO₂, and cyanuric acid can be obtained from urea or NH₃ + CO₂. In the stepwise mechanism of cyanuric acid from urea or NH₃ + CO₂, the energy barriers are relatively high, and the condition of high pressure and temperature does not decrease the energy barriers. Our theoretical model shows that cyanuric acid is actually acquired from OCNH via a one-step cycloaddition reaction.

Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6–311++G(2df,2pd) basis set and CCSD/6–31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

A DFT study on the mechanisms for the cycloaddition reactions between 1−aza-2-azoniaallene cations and acetylenes

The mechanisms of cycloaddition reactions between 1-aza-2-azoniaallene cations 1 and acetylenes 2 have been investigated using the global electrophilicity and nucleophilicity of the corresponding reactants as global reactivity indexes defined within the conceptual density functional theory. The reactivity and regioselectivity of these reactions were predicted by analysis of the energies, geometries, and electronic nature of the transition state structures. The theoretical results revealed that the reaction features a tandem process: an ionic 1,3-dipolar cycloaddition to produce the cycloadducts 3?H-pyrazolium salts 3 followed by a [1,2]-shift affording the thermodynamically more stable adducts 4 or 5. The mechanism of the cycloaddition reactions can be described as an asynchronous concerted pathway with reverse electron demand. The model reaction has also been investigated at the QCISD/6-31++G(d,p) and CCSD(T)/6-31++G(d,p)//B3LYP/6-31++G(d,p) levels as well as by the DFT. The polarizable continuum model, at the B3LYP/6-31++G(d,p) level of theory, was used to study solvent effects on all the studied reactions. In solvent dichloromethane, all the initial cycloadducts 3 were obtained via direct ionic process as the result of the solvent effect. The consecutive [1,2]-shift reaction, in which intermediates 3 are rearranged to the five-membered heterocycles 4/5, is proved to be a kinetically controlled reaction, and the regioselectivity can be modulated by varying the migrant. The LOL function and RDG function based on localized electron analysis were used to analysis the covalent bond and noncovalent interactions in order to unravel the mechanism of the title reactions.     

Mechanisms on electrical breakdown strength increment of polyethylene by acetophenone and its analogues addition: a theoretical study

A theoretical investigation is completed on the mechanism of electrical breakdown strength increment of polyethylene. It is shown that it is one of the most important factors for increasing electrical breakdown strength of polyethylene through keto-enol isomerization of acetophenone and its analogues at the ground state S₀ and the lowest triplet state T₁. The minimum structures and transition states of the keto- and the enol-tautomer of acetophenone and its analogues at the S₀ and T₁ states are obtained at the B3LYP/6-311+G(d,p) level, as well as the harmonic vibration frequencies of the equilibrium geometries and the minimum energy path (MEP) by the intrinsic reaction coordinate (IRC) theory at the same level. The two C–C bond cleavage reaction channels have been identified in acetophenone. The calculated results show that the energy barriers of keto-enol isomerization of acetophenone and its analogues at S₀ and T₁ states are much smaller than the average C-C bond energy of polyethylene, and the acetophenone doping or bond linked into polyethylene can increase the electrical breakdown strength and inhibit polyethylene electrical tree initiation and aging.

Study on the role of SBA-15 in the oxidative dehydrogenation of n-butane over vanadia catalyst using density functional theory

The first step in the mechanism of n-butane oxidative dehydrogenation (ODH) on a V₄O₁₀ cluster and V₄O₁₀ supported SBA-15 is examined using DFT method. The activation and adsorption energies, oxidation state of V atoms are calculated. Over V₄O₁₀ the obtained results indicate that the activation of C-H bond of methylene group can occur at both the terminal and the bridging oxygen atoms with similar barrier (21.5–22.5 kcal mol^?1). The role of SBA-15 (with and without modification by Al) in n-butane adsorption step has been studied in detail. SBA-15 itself has mild effect on the reaction process, but the substitution of silicon atoms by aluminum atoms results in an active supporter for V₂O₅ in ODH reaction. In that, the ratio of Si/Al will decide the direction of initial interaction steps between n-butane and catalyst surface and it will result in the selectivity of the reaction products.

Role of gold in a complex cascade reaction involving two electrocyclization steps

Quantum chemical computations (B3LYP/LACVP**) were applied to assess the impact of Au(I) complexation on activation barriers for sequential electrocyclization reactions (one a 1,2-dihydroazete ring-opening and another a pentadienyl cation ring-closure) proposed to occur during a complex reaction cascade that converts alkynes and imines to cyclopentenimines.

Thermal, electronic and ductile properties of lead-chalcogenides under pressure

Fully relativistic pseudo-potential ab-initio calculations have been performed to investigate the high pressure phase transition, elastic and electronic properties of lead-chalcogenides including the less known lead polonium. The calculated ground state parameters, for the rock-salt structure show good agreement with the experimental data. PbS, PbSe, PbTe and PbPo undergo a first-order phase transition from rock-salt to CsCl structure at 19.4, 15.5, 11.5 and 7.3 GPa, respectively. The elastic properties have also been calculated. The calculations successfully predicted the location of the band gap at L-point of Brillouin zone and the band gap for each material at ambient pressure. It is observed that unlike other lead-chalcogenides, PbPo is semi-metal at ambient pressure. The pressure variation of the energy gap indicates that these materials metalize under pressure. The electronic structures of these materials have been computed in parent as well as in high pressure B2 phase.

DFT studies of the adsorption and dissociation of H2O on the Al13 cluster: origins of this reactivity and the mechanism for H2 release

A theoretical study of the chemisorption and dissociation pathways of water on the Al₁₃ cluster was performed using the hybrid density functional B3LYP method with the 6-311+G(d, p) basis set. The activation energies, reaction enthalpies, and Gibbs free energy of activation for the reaction were determined. Calculations revealed that the H₂O molecule is easily adsorbed onto the Al₁₃ surface, forming adlayers. The dissociation of the first H₂O molecule from the bimolecular H₂O structure via the Grotthuss mechanism is the most kinetically favorable among the five potential pathways for O–H bond breaking. The elimination of H₂ in the reaction of an H₂O molecule with a hydrogen atom on the Al cluster via the Eley–Rideal mechanism has a lower activation barrier than the elimination of H₂ in the reaction of two adsorbed H atoms or the reaction of OH and H. Following the adsorption and dissociation of H₂O, the structure of Al₁₃ is distorted to varying degrees.

Reactions of ketones with aromatics in acid media. The effect of trifluoromethyl groups and the acidity media. A theoretical study

The reactions of acetone, 2,2,2-trifluoroacetone and hexafluoroacetone in methanesulfonic (MSA) and triflic acids (TFSA) with benzene have been studied at M06-2X/6-311+G(d,p) level using cluster-continuum model, where the carbonyl group is explicitly solvated by acid molecules. The introduction of a trifluoromethyl group into the ketone structure reduces the activation energy of the tetrahedral intermediates formation due to an increase of the electrophilicity of the carbonyl group and raises the activation and the reaction energies of the C-O bond cleavage in formed carbinol due to the destabilization of the corresponding carbocation. The introduction of the second trifluoromethyl group inhibits the hydroxyalkylation reaction due to a very strong increase of the reaction and activation energies of the C-O bond cleavage which becomes the rate determining step. The most important catalytic effect of TFSA compared to MSA is not the protonation of the ketone carbonyl, but the reduction of the activation and reaction energies of the carbinol C-O bond cleavage due to better protosolvation properties. Even for TFSA no complete proton transfer to carbonyl oxygen has been observed for free ketones. Therefore, the protonation energies of free ketones cannot be considered as a measure of ketone reactivity in the hydroxyalkylation reaction.

Zinc-, cadmium-, and mercury-containing one-dimensional tetraphenylporphyrin arrays: a DFT study

This work describes theoretical and experimental studies on glycerol esterification to obtain acetins focusing on the obtained isomers. The reaction of glycerol with acetic acid was carried out on Amberlyst 36 wet. Density functional theory calculations on the level of M06-2X functional and 6-311+G(d,p) basis set are carried out and the most stable structures of the reactants and products are located by considering a large number of conformers. The thermodynamics is discussed in terms of the calculated reaction Gibbs free energy. The AIM theory was used to characterize reactants and products. The glycerol esterification with acetic acid is found to be thermodynamically favored, with exothermal property. These agree well with experiments and allow us to explain the relative selectivity of products.

Nickel/zinc-catalyzed decarbonylative addition of anhydrides to alkynes: A DFT study

Density functional theory (DFT) was used to investigate the nickel- or nickel(0)/zinc- catalyzed decarbonylative addition of phthalic anhydrides to alkynes. All intermediates and transition states were optimized completely at the B3LYP/6-31+G(d,p) level. Calculated results indicated that the decarbonylative addition of phthalic anhydrides to alkynes was exergonic, and the total free energy released was ?87.6 kJ mol^?1. In the five-coordinated complexes M4a and M4b, the insertion reaction of alkynes into the Ni-C bond occurred prior to that into the Ni-O bond. The nickel(0)/zinc-catalyzed decarbonylative addition was much more dominant than the nickel-catalyzed one in whole catalytic decarbonylative addition. The reaction channel CA→M1'→T1'→M2'→T2'→M3a'→M4a'→T3a1'→M5a1' →T4a1'→M6a'→P was the most favorable among all reaction pathways of the nickel- or nickel(0)/zinc- catalyzed decarbonylative addition of phthalic anhydrides to alkynes. And the alkyne insertion reaction was the rate-determining step for this channel. The additive ZnCl₂ had a significant effect, and it might change greatly the electron and geometry structures of those intermediates and transition states. On the whole, the solvent effect decreased the free energy barriers.

Effect of surface hydroxyls on dimethyl ether synthesis over the γ-Al2O3 in liquid paraffin: a computational study

In a recent paper (Zuo et al., Appl Catal A 408:130–136, 2011), the mechanism of dimethyl ether (DME) synthesis from methanol dehydration over γ-Al₂O₃ (110) was studied using density functional theory (DFT). Using the same method, the effect of surface hydroxyls on γ-Al₂O₃ in liquid paraffin during DME synthesis from methanol dehydration is investigated. It is found that DME is mainly formed from two adsorbed CH₃O groups via methanol dehydrogenation on both dehydrated and hydrated γ-Al₂O₃ in liquid paraffin. No close correlation between catalytic activity and acid intensity was found. Before and after water adsorption at typical catalytic conditions (e.g., 553 K), the reaction rate is not obviously changed on γ-Al₂O₃(100) surface in liquid paraffin, but the reaction rate decreases by about 11 times on the (110) in liquid paraffin. Considering the area of the (110) and (100) surfaces under actual conditions, the catalytic activity decreased mainly because the Al3 sites on the (110) surface gradually become inactive. Catalytic activity decreased mainly due to surface hydrophilicity. The calculated results were consistent with the experiment.

Variation of the electronic dipole polarizability on the reaction path

The reaction force and the electronic flux, first proposed by Toro-Labbé et al. (J Phys Chem A 103:4398, 1999) have been expressed by the existing conceptual DFT apparatus. The critical points (extremes) of the chemical potential, global hardness and softness have been identified by means of the existing and computable energy derivatives: the Hellman-Feynman force, nuclear reactivity and nuclear stiffness. Specific role of atoms at the reaction center has been unveiled by indicating an alternative method of calculation of the reaction force and the reaction electronic flux. The electron dipole polarizability on the IRC has been analyzed for the model reaction HF + CO→HCOF. The electron polarizability determined on the IRC α _e (ξ) was found to be reasonably parallel to the global softness curve S(ξ). The softest state on the IRC (not TS) coincides with zero electronic flux.

Excited-state relaxation paths of oxo/hydroxy and N9H/N7H tautomers of guanine: a CC2 theoretical study

We performed a theoretical investigation, at the CC2/aug-cc-pVDZ level, of the ring-deformation mechanisms of four guanine tautomers (oxo, hydroxy, N9H, and N7H). The study showed that the optimized conical intersections S₀/S₁ are accessible through the ¹ππ* excited states of tautomers. The optimized conical intersections S₀/S₁, which show deformation at the pyrimidine ring, have high energies. This means that the relaxations of the ¹ππ* excited states via internal conversion are disfavored. For two tautomers we found crossing points ¹ππ*/¹πσ* of the excited-state reaction paths, revealing the possibility of a population of the ¹πσ* excited state by the ¹ππ* excited state.