Mechanistic study on Mo(CO)6‐catalyzed intramolecular [2 + 2] or [2 + 2 + 1] cycloaddition reaction of 5‐allenyl‐1‐ynes

By means of density functional theory, the Mo(CO)₆‐catalyzed intramolecular [2 + 2] or [2 + 2 + 1] cycloaddition reaction of 5‐allenyl‐1‐ynes was investigated. All the intermediates and transition states were optimized completely at B3LYP/6‐311++G(d,p) level (LANL2DZ(f) for Mo). Calculations indicate that the complexation of 5‐allenyl‐1‐ynes with Mo(CO)₆ occurred preferentially at the triple bond to give the complex M1 and then the complexation with the distal double bond of the allenes generates the complex M5 . In this reaction, Mo(CO)₆‐catalyzed intramolecular [2 + 2] cycloaddition is more favorable than [2 + 2 + 1] cycloaddition. The reaction pathway Mo(CO)₆ + R → M5 → T7 → M12 → M13 → T11 → M18 → P4 is the most favorable one, and the most dominant product predicted theoretically is P4 . The solvation effect is remarkable, and it decreases the reaction energy barriers. Copyright © 2011 John Wiley & Sons, Ltd.     

1,5‐electrocyclization versus 1,5‐proton shift in imidazolium allylides and 2‐phospha‐allylides: a DFT investigation

The competitive 1,5‐electrocyclization versus intramolecular 1,5‐proton shift in imidazolium allylides and imidazolium 2‐phosphaallylides has been investigated theoretically at the DFT (B3LYP/6‐311 + +G**//B3LYP/6‐31G**) level. 1,5‐Electrocyclization follows pericyclic mechanism and its activation barrier is lower than that for the pseudopericyclic mechanism by ～5–6 kcal mol^?1. The activation barriers for 1,5‐electrocyclization of imidazolium 2‐phosphaallylides are found to be smaller than those for their nonphosphorus analogues by ～3–5 kcal mol^?1. There appears to be a good correlation between the activation barrier for intramolecular 1,5‐proton shift and the density of the negative charge at C8, except for the ylides having fluorine substituent at this position ( 7b and 8b ). The presence of fluorine atom reduces the density of the negative charge at C8 (in 7b it becomes positively charged) and thus raises the activation barrier. The ylides 7f and 8f having CF₃ group at C8, in preference to the 1,5‐proton shift, follow an alternative route leading to different carbenes which is accompanied by the loss of HF. The carbenes Pr _{7 , 8b – e} resulting from intramolecular 1,5‐proton shift have a strong tendency to undergo intramolecular S_N2 type reaction, the activation barrier being 7–28 kcal mol^?1. Copyright © 2011 John Wiley & Sons, Ltd.     

A complete active-space self-consistent-field study on cubic N8

Summary Complete Active-Space Self-Consistent-Field (CAS-SCF) calculations for cubic N₈ are presented. We studied the N₈↔4N₂ reaction inD _4h symmetry and found its energy release and activation barrier with three different atomic basis sets. The energy release for this reaction is predicted to be around 526 kcal/mol, while the energy barrier to dissociation is estimated about 159 kcal/mol. These results are in substantial agreement with previousab initio estimates. The authors of this paper have agreed to not receive the proofs for correction.     

Stereoelectronic effects in Menshutkin‐type SN2 reactions: theoretical study based on through‐space/bond orbital interaction analysis

Through‐space/bond orbital interaction analysis has been applied to investigate the stereoelectronic effects on stabilizing the transition state of Menshutkin‐type S_N2 reactions. The mechanism of how the substituent effects work on accelerating the reactions has been demonstrated from orbital interaction perspective. The geometrical structures and Mulliken charge distributions have been compared to elucidate the substituent effects for the S_N2 reaction center. It is found that the substituents lower the activation energies by strengthening the orbital interactions in the S_N2 reaction process. When electron‐donating and electron‐accepting substituents (–C₆H₅ and –CHO) are introduced to the same central carbon at the reaction center, the symmetry allows the π–π* interactions among the donor and acceptor in the transition state. It stabilizes the transition state much more than the reactant complex. And the π–π* interactions are estimated to decrease about 2.28 kcal/mol of the energy for transition state. The σ‐like orbitals of the partial bond around the central carbon are reactive, and the σ–π* orbital interactions stabilize the reactant complex much more than the π–σ* interaction. When the σ–π* and π–σ* interactions are deleted from the system, the activation energy increases and turns close to the values of the systems which are without such substituents. It can be concluded that the π–π*, σ–π*, and π–σ* interactions cooperatively accelerates the S_N2 reaction by stabilizing its transition state. Copyright © 2013 John Wiley & Sons, Ltd.     

Reactions C<Subscript>2</Subscript>H<Subscript>2</Subscript> + OH and C<Subscript>2</Subscript> + H<Subscript>2</Subscript>O: Ab initio study of the potential energy surfaces

The stationary points of the potential energy surfaces for the reactions C₂H₂ + OH and C₂ + H₂O are calculated using density functional theory and the coupled cluster method. The relative energies and geometric parameters of the stable intermediates and transition states are in good agreement with the results of independent studies. In most cases, the relative energies differ from the earlier published values by no more than 3 kcal/mol, whereas the rotational constants, by 1–2%. The mechanism of the reaction CCOH₂ → C₂ + H₂O is studied in detail. The possible sources of errors in the calculation methods are examined.     

Theoretical study on the mechanism and kinetics of atmospheric reactions NH2OH + OOH and NH2CH3 + OOH

Mechanism and kinetics of NH₂OH + OOH and NH₂CH₃ + OOH reactions were studied at the B3LYP and M062X levels of theory using the 6-311++G(3df, 3pd) basis set. The NH₂OH + OOH and NH₂CH₃ + OOH reactions proceed through different paths which lead to different products. Transition state structure and activation energy of each path were calculated. The calculated activation energies of hydrogen abstraction reactions were smaller than 25 kcal/mol and of substitution reactions are in the range of 50–70 kcal/mol. The rate constants were calculated using transition state theory (TST) modified for tunneling effect at 273–2000 K.     

Understanding the molecular mechanism of the [3 + 2] cycloaddition reaction of benzonitrile oxide toward electron‐rich N‐vinylpyrrole: a DFT study

The [3 + 2] cycloaddition (32CA) reaction of benzonitrile oxide, BNO 2 , with an electron‐rich N‐vinylpyrrole derivative, NVP 3a , in the presence of dichloromethane, has been theoretically studied using density functional theory (DFT) methods at the B3LYP/6‐31G(d) level. This 32CA reaction presents a relatively high activation Gibbs free energy as a result of the low polar character of this zwitterionic‐type (zw‐type) reaction. Analyses of the calculated relative Gibbs free energies and transition state geometries indicate that the studied 32CA reaction, in excellent agreement with experimental outcomes, takes place in a complete regioselective manner as a consequence of the steric repulsions that appear at the most unfavorable transition state. An electron localization function (ELF) topological analysis of the bonding changes along this 32CA reaction supports a non‐concerted two‐stage one‐step molecular mechanism in which the formation of the O3‐ C5 single bond takes place at the end of the reaction after the complete formation of the C1‐C4 one. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamics of [n.3]paracyclophanes (n = 2 – 4) as studied by NMR. Obtaining separate Arrhenius parameters for two dynamic processes in [4.3]paracyclophane

Extending our earlier findings for [3.3]paracyclophane, NMR line shape studies of the conformational dynamics in [3.2] and [4.3]paracyclophanes are reported, of which the former is conformationally homogeneous and the latter occurs in two enantiomeric forms. For [3.2]paracyclophane, the Arrhenius activation energy E_a = 11.6 ± 0.1 kcal/mol and preexponential factor log (A/s^?1) = 12.92 ± 0.07 were found. In [4.3]paracyclophane, the conformational dynamics are quite complicated because, apart from interconversions of each enantiomer into itself proceeding via inversion of the propano bridge with rate constant k₁, the enantiomers mutually rearrange with rate constant k₂ due to inversion of the butano bridge. The determination of Arrhenius parameters from dynamic ¹H spectra of the aromatic protons for these two conformational processes (E_a = 11.2 ± 0.5 kcal/mol and log (A/s^?1) = 13.6 ± 0.5 for the former, and E_a = 9.7 ± 0.4 kcal/mol and log (A/s^?1) = 13.2 ± 0.4 for the latter) is the highlight of this work. In the investigated temperature range, in [4.3]paracyclophane, the occurrence of other conformational processes beyond those mentioned above can be excluded, because they would produce different line shape patterns than those actually observed. Copyright © 2013 John Wiley & Sons, Ltd.     

Conformational analysis of 4,4‐dimethyl‐4‐silathiane and its S‐oxides

4,4‐Dimethyl‐4‐silathiane and its S‐oxides [n = 0 ( 1 ), 1 ( 2 ), 2 ( 3 )] were studied experimentally by variable temperature dynamic NMR spectroscopy down to 103 K and the frozen ring inversion was revealed for all three compounds. The barriers for the degenerate ring inversion in 1 and 3 were measured to be 4.8 and 5.0 kcal/mol at the coalescence temperatures of 111 and 116 K, respectively, and practically coincide with the calculated barriers of 4.60 kcal/mol in 1 and 4.46 kcal/mol in 3 . The frozen equilibrium mixture 2‐ax/2‐eq contains 37% of the 2‐ax and 63% of the 2‐eq conformer. The ring inversion barrier proved to be ca. 4.8 kcal/mol. Calculations at the B3LYP/6‐311+G(d,p) level of theory showed the 2‐ax conformer to be 0.90 kcal/mol more stable than the 2‐eq conformer in the gas phase whereas in solution the relative stability of the conformers calculated using the PCM model at the same level of theory is inverted to become 0.19 (in CHCl₃) or 0.36 kcal/mol (in DMSO) in favor of the 2‐eq conformer. The chair–chair interconversion mechanism of sulfoxide 2 includes two intermediate energetically equivalent 1,4‐twist forms and the 2,5‐boat transition state: 2‐ax (chair) ? 2 (1,4‐twist) ? [ 2 (2,5‐boat)]^≠ ? 2 (1,4‐twist) ? 2‐eq (chair). The calculated ring inversion barriers are 5.1 ( 2‐ax → 2‐eq ) and 4.2 kcal/mol ( 2‐eq → 2‐ax ) in the gas phase, and 4.03 and 4.22 kcal/mol, respectively, in chloroform. Copyright © 2011 John Wiley & Sons, Ltd.     

Breathing viability into cyclonona‐3,5,7‐trienylidenes via α‐dimethyl and ά‐moieties at DFT

Cyclonona‐3,5,7‐trienylidene ( 1 ) changes from being a transition state (TS) to minimum states when substituted by α‐methyl groups and ?‐X, where X = CMe₂, NMe, PMe, O, S, cyclopropyl, and SiMe₂ ( 2 , 3 , 4 , 5 , 6 , 7 , 8 , respectively) at density functional theory. Specifically, the parent carbene 1 exhibits a negative vibrational force constant and proves to be an unreachable electrophilic TS while shows C_s symmetry with an NBO atomic charge of +0.70 on its carbenic center. It has a triplet ground state with a rather small singlet‐triplet energy gap (ΔE_s–t = ?4.1 kcal/mol). In contrast, all of its seven scrutinized derivatives enjoy reachable global minima, with C₁ symmetry, desired nucleophilicity, and singlet closed shell (S_cs) ground states (for all but 8 which remains triplet). Stability is indicated by relative ΔE_s–t values: 2 > 3 > 4 > 5 > 6 > 7 > 1 > 8 . The highest ΔE_s–t as well as NBO carbenic atomic negative charge (?0.74) are displayed by 2 . Our carbenes ( 2 , 3 , 4 , 5 , 6 , 7 ) appear more nucleophilic than the synthesized N‐heterocyclic carbenes (imidazol‐2‐ylidenes). Copyright © 2013 John Wiley & Sons, Ltd.     

The production of prop-2-ynylideneamine by thermolysis of N-chloropropargylamine,N-fluoropropargylamine and N-hydroxypropargylamine: a computational study

Ab initio molecular orbital calculations using MP2 and DFT/B3LYP methods at the 6-311++G(d,p) and aug-cc-pvdz basis sets were applied to characterise the kinetics of the thermal dissociation of HC≡CCH₂NHX [where X = OH(I), F(II) and Cl(III)] to produce Z- and E-prop-2-ynylideneamines (HC≡CCH=NH) (IV and V, respectively), which tautomerise to vinyl cyanide (CH₂=CHC≡N) (VI). The optimised geometries and electronic energies of reactants, transition states and products were estimated and discussed. A concerted proton migration and HX abstraction mechanism was proposed for the imine formation. The reliance of these properties on the elected levels of theory was discussed. The activation energies and barrier heights for the Z- and E-forms and their vinyl cyanide tautomers were estimated and analysed. The Z-form was computed to be more stable than the E-form. Using natural bond orbital calculations, the origin of the preference of the Z-form was attributed mainly to the N lone pair delocalisations. Vinyl cyanide was located to have a lower energy (33–35 kcal/mol) than prop-2-ynylideneamine. The provenance of the preference of the former and its tautomerisation mechanism will be addressed in a separate publication.     

Effects of α‐cyclopropyl on heterocyclic carbenes stability at DFT

α‐Cyclopropyl stability impacts on singlet and triplet heterocyclic carbenes with acyclic, cyclic, and cyclic‐unsaturated structures are compared and contrasted to di‐t‐butyl as well as t‐butylcyclopropylcarbenes through appropriate isodesmic reactions at B3LYP/AUG‐cc‐pVTZ level. Substitution of one of the t‐butyl groups of di‐t‐butylcarbene with a cyclopropyl alters the ground state multiplicity from triplet to singlet with a singlet–triplet energy separation (ΔE_s–t) of 7.2 kcal/mol. Additional heteroatom substitution increases ΔE_s–t values for the resulting α‐heteroatom cyclopropylcarbenes in the following order: amino > oxy > thio > phophino. α‐Cyclopropyl group stabilizes singlet states of all our carbenes two to three times more than their corresponding triplet states. The ΔE_s–t values of all the carbenes are increased through cyclization, while the introduction of unsaturation in the rings causes small and rather random changes. To probe the kinetic stability of the species, we calculated the transition states for the opening of cyclopropyl through 1,2‐C shift. Interestingly, the 4.1 kcal/mol energy barrier in cyclopropylcarbene is significantly increased in the presence of heteroatoms to 31.2 kcal/mol for aminocyclopropylcarbene. The reactivity of the species is discussed in terms of nucleophilicity and electrophilicity issues showing our carbenes, especially acyclic ones, more nucleophilic than the common N‐heterocyclic carbenes. Copyright © 2010 John Wiley & Sons, Ltd.     

Stereoselectivity and kinetics of [4 + 2] cycloaddition reaction of cyclopentadiene to para‐substituted E‐2‐arylnitroethenes

In spite of diversified electrophilicity of E‐2‐arylnitroethenes, their [4 + 2] cycloaddition reactions with cyclopentadiene leads to the corresponding 6‐endo‐aryl‐5‐exo‐nitronorbornenes and 6‐exo‐aryl‐5‐endo‐nitronorbornenes as the only reaction products. Stereoselectivity, substituent and solvent effects, and activation parameters, suggest that these reactions occur via a synchronous concerted mechanism on both competing pathways. The experimental results obtained are consistent with the data from B3LYP/6‐31G(d) calculations. Due to high electrophilicity of E‐2‐arylnitroethenes, the reactions studied should be considered as polar [4 + 2] cycloadditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Theoretical investigation on the excited‐state intramolecular proton transfer mechanism of 2‐(2′‐benzofuryl)‐3‐hydroxychromone

Spectroscopic studies on excited‐state proton transfer of a new chromophore 2‐(2′‐benzofuryl)‐3‐hydroxychromone (BFHC) have been reported recently. In the present work, based on the time‐dependent density functional theory (TD‐DFT), the excited‐state intramolecular proton transfer (ESIPT) of BFHC is investigated theoretically. The calculated primary bond lengths and angles involved in hydrogen bond demonstrate that the intramolecular hydrogen bond is strengthened. In addition, the phenomenon of hydrogen bond reinforce has also been testified based on infrared (IR) vibrational spectra as well as the calculated hydrogen bonding energies. Further, hydrogen bonding strengthening manifests the tendency of excited state proton transfer. Our calculated results reproduced absorbance and fluorescence emission spectra of experiment, which verifies that the TD‐DFT theory we used is reasonable and effective. The calculated Frontier Molecular Orbitals (MOs) further demonstrate that the excited state proton transfer is likely to occur. According to the calculated results of potential energy curves along O―H coordinate, the potential energy barrier of about 14.5 kcal/mol is discovered in the S₀ state. However, a lower potential energy barrier of 5.4 kcal/mol is found in the S₁ state, which demonstrates that the proton transfer process is more likely to happen in the S₁ state than the S₀ state. In other words, the proton transfer reaction can be facilitated based on the photo‐excitation effectively. Moreover, the phenomenon of fluorescence quenching could be explained based on the ESIPT mechanism. Copyright © 2015 John Wiley & Sons, Ltd.     

过氧硝酸分解势能面的理论研究

使用密度泛函理论B3LYP/6-311+ G(2d,2p)研究了过氧硝酸的最低能量结构．采用耦合簇方法CCSD(T)/aug-cc-pVDZ首次分别扫描了过氧硝酸沿氧-氮和氧-氧键的分解势能面．计算结果表明在氧-氮势能面上,当O3—N4键长是2.82 ?时,对应的疏松过渡态的能垒是25.6 kcal/mol;在氧$-$氧键的势能面上,当O2—O3键长是2.35 ?时,对应的疏松过渡态的能垒是37.4 kcal/mol．这表明过氧硝酸更容易分解为HO₂和NO₂．     

Conformational analysis of 3,3‐dimethyl‐3‐silathiane, 2,3,3‐trimethyl‐3‐silathiane and 2‐trimethylsilyl‐3,3‐dimethyl‐3‐silathiane—preferred conformers,barriers to ring inversion and substituent effects

The first conformational analysis of 3‐silathiane and its C‐substituted derivatives, namely, 3,3‐dimethyl‐3‐silathiane 1 , 2,3,3‐trimethyl‐3‐silathiane 2 , and 2‐trimethylsilyl‐3,3‐dimethyl‐3‐silathiane 3 was performed by using dynamic NMR spectroscopy and B3LYP/6‐311G(d,p) quantum chemical calculations. From coalescence temperatures, ring inversion barriers ΔG^≠ for 1 and 2 were estimated to be 6.3 and 6.8 kcal/mol, respectively. These values are considerably lower than that of thiacyclohexane (9.4 kcal/mol) but slightly higher than the one of 1,1‐dimethylsilacyclohexane (5.5 kcal/mol). The conformational free energy for the methyl group in 2 (?ΔG° = 0.35 kcal/mol) derived from low‐temperature ¹³C NMR data is fairly consistent with the calculated value. For compound 2 , theoretical calculations give ΔE value close to zero for the equilibrium between the 2 ‐Me_ax and 2 ‐Me_eq conformers. The calculated equatorial preference of the trimethylsilyl group in 3 is much more pronounced (?ΔG° = 1.8 kcal/mol) and the predominance of the 3 ‐SiMe_{3 eq} conformer at room temperature was confirmed by the simulated ¹H NMR and 2D NOESY spectra. The effect of the 2‐substituent on the structural parameters of 2 and 3 is discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

A DFT study on the (2 + 3) cycloaddition reactions of 2‐nitropropene‐1 with Z‐C,N‐diarylnitrones

B3LYP/6‐31G* calculations for competing (2 + 3)‐cycloaddition pathways for 2‐nitropropene‐1 (1) to Z‐C, N‐diarylnitrones ( 2a – e ) suggest a concerted reaction mechanism. However, the results point to the strongly asymmetric nature of transition complexes. Increasing polarity of the reaction environment and presence of electron‐donating substituents in the nitrone phenyl rings contribute to the higher asymmetry of these structures. Copyright © 2009 John Wiley & Sons, Ltd.     

P‐Heterocyclic silylenes: a survey of stability with density functional theory

The effects of phosphorous atom on the stability, multiplicity, and reactivity of six‐member cyclic silylenes are investigated at B3LYP/AUG‐cc‐pVTZ//B3LYP/6‐31+G* and MP2/6‐311++G**//B3LYP/6‐31+G* coupled with appropriate isodesmic reactions. From a thermodynamic point of view, 1H‐2‐silaphosphinine‐2‐ylidene ( 1_a ) and 1H‐4‐silaphosphinine‐4‐ylidene ( 2_a ) are relatively the most stable with singlet–triplet energy gaps (ΔE_S–T) of 37.0 and 28.1 kcal/mol, respectively. The calculated energy barrier for the 1,2‐H shift of 1_a to the corresponding 2‐silapyridine ( 1 ) is 26.5 kcal/mol, which is lower than the 28.8 kcal/mol required for the 1,4‐H shift of 2_a to the corresponding 4‐silapyridine ( 2 ). In contrast to the previous reports, isodesmic reactions indicate that π‐donor/σ‐donor phosphorous destabilizes the singlet while stabilizes the triplet state. Both 1_a and 2_a silylenes appear invulnerable to the head‐to‐head as well as the head‐to‐tail dimerization, inviting experimental explorations. Copyright © 2011 John Wiley & Sons, Ltd.     

Palladium‐catalyzed cyclization of 2‐alkynyl‐N‐ethanoyl anilines to indoles: synthesis,structural, spectroscopic,and mechanistic study

This study reports a facial regio‐selective synthesis of 2‐alkyl‐N‐ethanoyl indoles from substituted‐N‐ethanoyl anilines employing palladium (II) chloride, which acts as a cyclization catalyst. The mechanistic trait of palladium‐based cyclization is also explored by employing density functional theory. In a two‐step mechanism, the palladium, which attaches to the ethylene carbons, promotes the proton transfer and cyclization. The gas‐phase barrier height of the first transition state is 37 kcal/mol, indicating the rate‐determining step of this reaction. Incorporating acetonitrile through the solvation model on density solvation model reduces the barrier height to 31 kcal/mol. In the presence of solvent, the electron‐releasing (–CH₃) group has a greater influence on the reduction of the barrier height compared with the electron‐withdrawing group (–Cl). These results further confirm that solvent plays an important role on palladium‐catalyzed proton transfer and cyclization. For unveiling structural, spectroscopic, and photophysical properties, experimental and computational studies are also performed. Thermodynamic analysis discloses that these reactions are exothermic. The highest occupied molecular orbital?lowest unoccupied molecular orbital gap (4.9–5.0 eV) confirms that these compounds are more chemically reactive than indole. The calculated UV–Vis spectra by time‐dependent density functional theory exhibit strong peaks at 290, 246, and 232 nm, in good agreement with the experimental results. Moreover, experimental and computed ¹H and ¹³C NMR chemical shifts of the indole derivatives are well correlated. Copyright © 2015 John Wiley & Sons, Ltd.     

Luminescence of CaGa<Subscript>2</Subscript>Se<Subscript>4</Subscript>:Eu crystals

We have studied photoluminescence and thermoluminescence (PL and TL) in CaGa₂Se₄:Eu crystals in the temperature range 77–400 K. We have established that broadband photoluminescence with maximum at 571 nm is due to intracenter transitions 4f⁶ 5d–4f⁷ (⁸S_7/2) of the Eu²⁺ ions. From the temperature dependence of the intensity (log I–10³/T), we determined the activation energy (E _a = 0.04 eV) for thermal quenching of photoluminescence. From the thermoluminescence spectra, we determined the trap depths: 0.31, 0.44, 0.53, 0.59 eV. The lifetime of the excited state 4f6 5d of the Eu²⁺ ions in the CaGa₂Se₄ crystal found from the luminescence decay kinetics is 3.8 μsec. Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 76, No. 1, pp. 112–116, January–February, 2009.