A combined density functional and intrinsic reaction coordinate study on the ground state energy surface of H2CO

Recovering Valuable Chemicals from Polypropylene Waste via a Mild Catalyst-Free Hydrothermal Process

Waste polypropylene (PP) presents a significant environmental challenge, owing to its refractory nature and inert C-C backbone. In this study, we introduce a practical chemical recovery strategy from PP waste using a mild catalyst-free hydrothermal treatment (HT). The treatment converts 64.1% of the processed PP into dissolved organic products within 2 h in an air atmosphere at 160 °C. Higher temperatures increase the PP conversion efficiency. Distinct electron absorption and emission characteristics of the products are identified by spectral analysis. Fourier transform-ion cyclotron resonance-mass spectrometry (FT-ICR-MS) reveals the oxidative cracking of PP into shorter-chain homologues (10-50 carbon atoms) containing carboxylic and carbonyl groups. Density functional theory (DFT) calculations support a reaction pathway involving thermal C-H oxidation at the tertiary carbon sites in the polymer chain. The addition of 1% H2O2 further enhances the oxidation reaction to produce valuable short-chain acetic acids, enabling gram-scale recycling of both pure PP and disposable surgical masks from the real world. Techno-economic analysis (TEA) and environmental life cycle costing (E-LCC) analysis suggest that this hydrothermal oxidation recovery technology is financially viable, which shows significant potential in tackling the ongoing plastic pollution crisis and advancing plastic treatment methodologies toward a circular economy paradigm.

Retro-Cope elimination of cyclic alkynes: reactivity trends and rational design of next-generation bioorthogonal reagents

The retro-Cope elimination reaction between dimethylhydroxylamine (DMHA) and various cyclic alkynes has been quantum chemically explored using DFT at ZORA-BP86/TZ2P. The purpose of this study is to understand the role of the following three unique activation modes on the overall reactivity, that is (i) additional cycloalkyne predistortion via fused cycles, (ii) exocyclic heteroatom substitution on the cycloalkyne, and (iii) endocyclic heteroatom substitution on the cycloalkyne. Trends in reactivity are analyzed and explained by using the activation strain model (ASM) of chemical reactivity. Based on our newly formulated design principles, we constructed a priori a suite of novel bioorthogonal reagents that are highly reactive towards the retro-Cope elimination reaction with DMHA. Our findings offer valuable insights into the design principles for highly reactive bioorthogonal reagents in chemical synthesis.

Reaction Kinetics of the Benzylation of Adenine in DMSO: Regio-Selectivity Guided by Entropy

The factors governing the regio-selectivity of the alkylation of adenine have been of interest for many years due to the biological importance of adenine derivatives, however, no reaction kinetic studies have been conducted. Herein, we report the rate constants and activation parameters of the benzylation of adenine under basic conditions in DMSO in the absence and presence of 15-crown-5 ether using real-time 1H NMR spectroscopy. The reaction is second-order for the formation of the N9- and N3-benzyladenine products, with a regio-selectivity factor 2.3 in favour of the N9-adduct. The Gibbs free energy of activation amounts to 87±2 kJ mol-1 for both reactions. The formation of the N9-adduct is more activated by 7 kJ mol-1, but its effect is offset by a less negative activation entropy, demonstrating that the long-contested reason for the regioselectivity in the benzylation of adenine is dominated by compensation of entropy and enthalpy in the transition state. The kinetic parameters obtained in the presence of the 15-crown-5 ether indicate that the crown ether forms a complex with an adenine-sodium ion-pair, increasing the activation barrier. However, the Gibbs free energy in the absence and presence of the crown ether remains constant.

Substituent effects on the selectivity of ambimodal [6+4]/[4+2] cycloaddition

In this work, we report a density functional theory (DFT) study and a dynamical trajectory study of substituent effects on the ambimodal [6+4]/[4+2] cycloaddition proposed for 1,3,5,10,12-cycloheptadecapentaene, referred to as cycloheptadecapentaene. The proposed cycloaddition proceeds through an ambimodal transition state, which results in both a [6+4] adduct a [4+2] adduct directly. The [6+4] adduct can be readily converted to the [4+2] adduct via a Cope rearrangement. We study the selectivity of the reaction with regard to the position of substituents, steric effects of substituents, and electronic effects of substituents. In the dynamical trajectory study, we find that nitro-substituted reactants lead to a new product from the ambimodal transition state via the hetero Diels-Alder reaction, and this new product can then be converted to a [4+2] adduct by a hetero [3, 3]-sigmatropic rearrangement. These results may provide insights for designing more bridged heterocyclic compounds.

Theoretical Studies on the Isomerization Kinetics of Low-Lying Isomers of the SiC4H2 System

The low-lying isomers of SiC4H2 are investigated to understand the kinetics of isomerization pathways using density functional theory. In our earlier work, we studied the various possible isomers (J. Phys. Chem. A, 2020, 124, 987-1002) and the chemical bonding of low-lying isomers of SiC4H2 (J. Phys. Chem. A, 2022, 126, 9366-9374). Among them, four isomers, 1-ethynyl-3-silacycloprop-1-en-3-ylidene (1), 3-silapent-1,4-diyn-3-ylidene (2), 1-silapent-1,2,3,4-tetraen-1-ylidene (4), and 1-silapent-2,4-diyn-1-ylidene (5) have already been identified in the laboratory. The previously known theoretical isomer 2-methylene-1-silabicyclo[1.1.0]but-1(3)-en-4-ylidene (3) and the newly identified unknown isomer through the present kinetic studies 5-silabicyclo[2.1.0]pent-1(4),2-dien-5-ylidene (N6) remain elusive in the laboratory to date. The isomerization pathways of the low-lying isomers of SiC4H2 are predicted through the transition state structures. Intrinsic reaction coordinate analysis identifies the minimum energy reaction pathways connecting the transition state from one isomer to another of the investigated system. The present kinetic data reveal the isomerization of global minimum energy isomer 1 to thermodynamically stable low-lying isomers, 2 and 5. Interestingly, isomer 3 interconverts to the experimentally known low-energy isomer 4, the second most thermodynamically stable isomer among them. The thermodynamic and kinetic parameters of the low-lying isomers of SiC4H2 are also documented in this work. The rate coefficient and equilibrium constant for isomerization reactions are calculated using the Rice-Ramsperger-Kassel-Marcus theory. The equilibrium constant delineates the difficulties in forming N6 and 3 through the isomerization pathways. Furthermore, ab initio molecular dynamics studies dictate the stability of low-lying isomers of SiC4H2 within the time scale of the simulation.

SN 2 versus E2 Competition of Cyclic Ethers

We have quantum chemically studied the influence of ring strain on the competition between the two mechanistically different SN 2 and E2 pathways using a series of archetypal ethers as substrate in combination with a diverse set of Lewis bases (F- , Cl- , Br- , HO- , H3 CO- , HS- , H3 CS- ), using relativistic density functional theory at ZORA-OLYP/QZ4P. The ring strain in the substrate is systematically increased on going from a model acyclic ether to a 6- to 5- to 4- to 3-membered ether ring. We have found that the activation energy of the SN 2 pathway sharply decreases when the ring strain of the system is increased, thus on going from large to small cyclic ethers, the SN 2 reactivity increases. In contrast, the activation energy of the E2 pathway generally rises along this same series, that is, from large to small cyclic ethers. The opposing reactivity trends induce a mechanistic switch in the preferred reaction pathway for strong Lewis bases from E2, for large cyclic substrates, to SN 2, for small cyclic substrates. Weak Lewis bases are unable to overcome the higher intrinsic distortivity of the E2 pathway and, therefore, always favor the less distortive SN 2 reaction.

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

Rational design of iron catalysts for C-X bond activation

We have quantum chemically studied the iron-mediated CX bond activation (X = H, Cl, CH₃ ) by d⁸ -FeL₄ complexes using relativistic density functional theory at ZORA-OPBE/TZ2P. We find that by either modulating the electronic effects of a generic iron-catalyst by a set of ligands, that is, CO, BF, PH₃ , BN(CH₃ )₂ , or by manipulating structural effects through the introduction of bidentate ligands, that is, PH₂ (CH₂ )_n PH₂ with n = 6-1, one can significantly decrease the reaction barrier for the CX bond activation. The combination of both tuning handles causes a decrease of the CH activation barrier from 10.4 to 4.6 kcal mol^-1 . Our activation strain and Kohn-Sham molecular orbital analyses reveal that the electronic tuning works via optimizing the catalyst-substrate interaction by introducing a strong second backdonation interaction (i.e., "ligand-assisted" interaction), while the mechanism for structural tuning is mainly caused by the reduction of the required activation strain because of the pre-distortion of the catalyst. In all, we present design principles for iron-based catalysts that mimic the favorable behavior of their well-known palladium analogs in the bond-activation step of cross-coupling reactions.

Palladium-catalyzed activation of HnA–AHn bonds (AHn = CH3, NH2, OH, F)

We have quantum chemically studied activation of HnA–AHn bonds (AHn = CH3, NH2, OH, F) by PdLn catalysts with Ln = no ligand, PH3, (PH3)2, using relativistic density functional theory at ZORA-BLYP/TZ2P. The activation energy associated with the oxidative addition step decreases from H3C–CH3 to H2N–NH2 to HO–OH to F–F, where the activation of the F–F bond is barrierless. Activation strain and Kohn–Sham molecular orbital analyses reveal that the enhanced reactivity along this series of substrates originates from a combination of (i) reduced activation strain due to a weaker HnA–AHn bond; (ii) decreased Pauli repulsion as a result of a difference in steric shielding of the HnA–AHn bond; and (iii) enhanced backbonding interaction between the occupied 4d atomic orbitals of the palladium catalyst and σ* acceptor orbital of the substrate.

How Bases Catalyze Diels-Alder Reactions

We have quantum chemically studied the base-catalyzed Diels-Alder (DA) reaction between 3-hydroxy-2-pyrone and N-methylmaleimide using dispersion-corrected density functional theory. The uncatalyzed reaction is slow and is preceded by the extrusion of CO2 via a retro-DA reaction. Base catalysis, for example, by triethylamine, lowers the reaction barrier up to 10 kcal mol-1 , causing the reaction to proceed smoothly at low temperature, which quenches the expulsion of CO2 , yielding efficient access to polyoxygenated natural compounds. Our activation strain analyses reveal that the base accelerates the DA reaction via two distinct electronic mechanisms: i) by the HOMO-raising effect, which enhances the normal electron demand orbital interaction; and ii) by donating charge into 3-hydroxy-2-pyrone which accumulates in its reactive region and promotes strongly stabilizing secondary electrostatic interactions with N-methylmaleimide.

[3 + 2]-Cycloadditions with Porphyrin β,β'-Bonds: Theoretical Basis of the Counterintuitive meso-Aryl Group Influence on the Rates of Reaction

Removal of a β,β'-bond from meso-tetraarylporphyrin using [3 + 2]-cycloadditions generates meso-tetraarylhydroporphyrins. Literature evidence indicates that meso-tetraphenylporphyrins react more sluggishly with 1,3-dipoles such as ylides and OsO₄ (in the presence of pyridine) than meso-tetrakis(pentafluorophenyl)porphyrin. The trend is counterintuitive for the reaction with OsO₄, as this formal oxidation reaction is expected to proceed more readily with more electron-rich substrates. This work presents a density functional theory-based computational study of the frontier molecular orbital (FMO) interactions and reaction profile thermodynamics involved in the reaction of archetypical cycloaddition reactions (a simple ylide, OsO₄, OsO₄·py, OsO₄·(py)₂, and ozone) with the β,β'-double bonds of variously fluorinated meso-arylporphyrins. The trend observed for the Type I cycloaddition of an ylide is straightforward, as lowering the LUMO of the porphyrin with increasing meso-phenyl-fluorination also lowers the reaction barrier. The corresponding simple FMO analyses of Type III cycloadditions do not correctly model the reaction energetics. This is because increasing fluorination leads to lowering of the porphyrin HOMO-2, thus increasing the reaction barrier. However, coordination of pyridine to OsO₄ preorganizes the transition state complex; lowering of the energy barrier by the preorganization exceeds the increase in repulsive orbital interactions, overall accelerating the cycloaddition and rationalizing the counterintuitive experimental findings.

Rational Tuning of the Reactivity of Three-Membered Heterocycle Ring Openings via SN 2 Reactions

The development of small-molecule covalent inhibitors and probes continuously pushes the rapidly evolving field of chemical biology forward. A key element in these molecular tool compounds is the "electrophilic trap" that allows a covalent linkage with the target enzyme. The reactivity of this entity needs to be well balanced to effectively trap the desired enzyme, while not being attacked by off-target nucleophiles. Here we investigate the intrinsic reactivity of substrates containing a class of widely used electrophilic traps, the three-membered heterocycles with a nitrogen (aziridine), phosphorus (phosphirane), oxygen (epoxide) or sulfur atom (thiirane) as heteroatom. Using quantum chemical approaches, we studied the conformational flexibility and nucleophilic ring opening of a series of model substrates, in which these electrophilic traps are mounted on a cyclohexene scaffold (C6 H10 Y with Y=NH, PH, O, S). It was revealed that the activation energy of the ring opening does not necessarily follow the trend that is expected from C-Y leaving-group bond strength, but steeply decreases from Y=NH, to PH, to O, to S. We illustrate that the HOMONu -LUMOSubstrate interaction is an all-important factor for the observed reactivity. In addition, we show that the activation energy of aziridines and phosphiranes can be tuned far below that of the corresponding epoxides and thiiranes by the addition of proper electron-withdrawing ring substituents. Our results provide mechanistic insights to rationally tune the reactivity of this class of popular electrophilic traps and can guide the experimental design of covalent inhibitors and probes for enzymatic activity.

Rational Design of Bimetallic Metal Chalcogenide Clusters for CO2 Dissociation

Thermochemical dissociation of CO₂ on pure, ligated, and mixed transition metal (W, Cu) chalcogenide clusters are investigated using the first-principles gradient-corrected density functional approach. It is shown that although the pure and ligated metal chalcogenide clusters exhibit significantly high barriers for CO₂ dissociation, the computed barriers for the mixed clusters are relatively lower. The lowest barrier is obtained for the Cu₃W₃Se₈ cluster, which shows a dramatically reduced barrier height of only 0.41 eV. Detailed analysis reveals that the substitution of W by Cu sites leads to a charge transfer from Cu to W sites, resulting in locally active W sites. The lowering of the CO₂ dissociation barriers can be attributed to the facile transfer of charge from the locally active W sites and also due to the alteration of the binding energy of CO₂ to the charged W sites. Our studies provide an alternate strategy to design novel thermochemical catalysts for CO₂ adsorption and subsequent dissociation.

How Ionization Catalyzes Diels-Alder Reactions

The catalytic effect of ionization on the Diels-Alder reaction between 1,3-butadiene and acrylaldehyde has been studied using relativistic density functional theory (DFT). Removal of an electron from the dienophile, acrylaldehyde, significantly accelerates the Diels-Alder reaction and shifts the reaction mechanism from concerted asynchronous for the neutral Diels-Alder reaction to stepwise for the radical-cation Diels-Alder reaction. Our detailed activation strain and Kohn-Sham molecular orbital analyses reveal how ionization of the dienophile enhances the Diels-Alder reactivity via two mechanisms: (i) by amplifying the asymmetry in the dienophile's occupied π-orbitals to such an extent that the reaction goes from concerted asynchronous to stepwise and thus with substantially less steric (Pauli) repulsion per reaction step; (ii) by enhancing the stabilizing orbital interactions that result from the ability of the singly occupied molecular orbital of the radical-cation dienophile to engage in an additional three-electron bonding interaction with the highest occupied molecular orbital of the diene.

SN2 versus SN2' Competition

We have quantum chemically explored the competition between the S_N2 and S_N2' pathways for X^- + H₂C═CHCH₂Y (X, Y = F, Cl, Br, I) using a combined relativistic density functional theory and coupled-cluster theory approach. Bimolecular nucleophilic substitution reactions at allylic systems, i.e., C^γ═C^β-C^α-Y, bearing a leaving-group at the α-position, proceed either via a direct attack at the α-carbon (S_N2) or via an attack at the γ-carbon, involving a concerted allylic rearrangement (S_N2'), in both cases leading to the expulsion of the leaving-group. Herein, we provide a physically sound model to rationalize under which circumstances a nucleophile will follow either the aliphatic S_N2 or allylic S_N2' pathway. Our activation strain analyses expose the underlying physical factors that steer the S_N2/S_N2' competition and, again, demonstrate that the concepts of a reaction's "characteristic distortivity" and "transition state acidity" provide explanations and design tools for understanding and predicting reactivity trends in organic synthesis.

Bridging Plastic Recycling and Organic Catalysis: Photocatalytic Deconstruction of Polystyrene via a C–H Oxidation Pathway

Chemical recycling of synthetic polymers represents a promising strategy to deconstruct plastic waste and make valuable products. Inspired by small-molecule C–H bond activation, a visible-light-driven reaction is developed to deconstruct polystyrene (PS) into ∼40% benzoic acid as well as ∼20% other monomeric aromatic products at 50 °C and ambient pressure. The practicality of this strategy is demonstrated by deconstruction of real-world PS foam on a gram scale. The reaction is proposed to proceed via a C–H bond oxidation pathway, which is supported by theoretical calculations and experimental results. Fluorescence quenching experiments also support efficient electron transfer between the photocatalyst and the polymer substrate, providing further evidence for the proposed mechanism. This study introduces concepts from small-molecule catalysis to polymer deconstruction and provides a promising method to tackle the global crisis of plastic pollution.

C(spn )-X (n=1-3) Bond Activation by Palladium

We have studied the palladium-mediated activation of C(spn )-X bonds (n = 1-3 and X = H, CH3 , Cl) in archetypal model substrates H3 C-CH2 -X, H2 C=CH-X and HC≡C-X by catalysts PdLn with Ln = no ligand, Cl- , and (PH3 )2 , using relativistic density functional theory at ZORA-BLYP/TZ2P. The oxidative addition barrier decreases along this series, even though the strength of the bonds increases going from C(sp3 )-X, to C(sp2 )-X, to C(sp)-X. Activation strain and matching energy decomposition analyses reveal that the decreased oxidative addition barrier going from sp3 , to sp2 , to sp, originates from a reduction in the destabilizing steric (Pauli) repulsion between catalyst and substrate. This is the direct consequence of the decreasing coordination number of the carbon atom in C(spn )-X, which goes from four, to three, to two along this series. The associated net stabilization of the catalyst-substrate interaction dominates the trend in strain energy which indeed becomes more destabilizing along this same series as the bond becomes stronger from C(sp3 )-X to C(sp)-X.

Parameter free evaluation of SN2 reaction rates for halide substitution in halomethane

We estimate the kinetic constants of a series of archetypal S_N2 reactions, i.e., the nucleophilic substitutions of halides in halomethane. A parameter free, multiscale approach recently developed [Campeggio et al., Phys. Chem. Chem. Phys., 2020, 22, 3455] is employed. The protocol relies on quantum mechanical calculations for the description of the energy profile along the intrinsic reaction coordinate, which is then mapped onto a reaction coordinate conveniently built for the reactive process. A Kramers-Klein equation is used to describe the stochastic time evolution of the reaction coordinate and its velocity; friction is parameterized using a hydrodynamic model and Kramers theory is used to derive the rate constant of the reaction. The method is here applied to six S_N2 reactions in water at 295.15 K, which differ in the nucleophile and the leaving group. The computed reaction rates are in good agreement with the experimental data and correlate well with the trends observed for the activation energies.

How Solvation Influences the SN2 versus E2 Competition

We have quantum chemically investigated how solvation influences the competition between the S_N2 and E2 pathways of the model F^– + C₂H₅Cl reaction. The system is solvated in a stepwise manner by going from the gas phase, then via microsolvation of one to three explicit solvent molecules, then last to bulk solvation using relativistic density functional theory at (COSMO)-ZORA-OLYP/QZ4P. We explain how and why the mechanistic pathway of the system shifts from E2 in the gas phase to S_N2 upon strong solvation of the Lewis base (i.e., nucleophile/protophile). The E2 pathway is preferred under weak solvation of the system by dichloromethane, whereas a switch in reactivity from E2 to S_N2 is observed under strong solvation by water. Our activation strain and Kohn–Sham molecular orbital analyses reveal that solvation of the Lewis base has a significant impact on the strength of the Lewis base. We show how strong solvation furnishes a weaker Lewis base that is unable to overcome the high characteristic distortivity associated with the E2 pathway, and thus the S_N2 pathway becomes viable.

Theoretical studies on triplet-state driven dissociation of formaldehyde by quasi-classical molecular dynamics simulation on machine-learning potential energy surface

The H-atom dissociation of formaldehyde on the lowest triplet state (T1) is studied by quasi-classical molecular dynamic simulations on the high-dimensional machine-learning potential energy surface (PES) model. An atomic-energy based deep-learning neural network (NN) is used to represent the PES function, and the weighted atom-centered symmetry functions are employed as inputs of the NN model to satisfy the translational, rotational, and permutational symmetries, and to capture the geometry features of each atom and its individual chemical environment. Several standard technical tricks are used in the construction of NN-PES, which includes the application of clustering algorithm in the formation of the training dataset, the examination of the reliability of the NN-PES model by different fitted NN models, and the detection of the out-of-confidence region by the confidence interval of the training dataset. The accuracy of the full-dimensional NN-PES model is examined by two benchmark calculations with respect to ab initio data. Both the NN and electronic-structure calculations give a similar H-atom dissociation reaction pathway on the T1 state in the intrinsic reaction coordinate analysis. The small-scaled trial dynamics simulations based on NN-PES and ab initio PES give highly consistent results. After confirming the accuracy of the NN-PES, a large number of trajectories are calculated in the quasi-classical dynamics, which allows us to get a better understanding of the T1-driven H-atom dissociation dynamics efficiently. Particularly, the dynamics simulations from different initial conditions can be easily simulated with a rather low computational cost. The influence of the mode-specific vibrational excitations on the H-atom dissociation dynamics driven by the T1 state is explored. The results show that the vibrational excitations on symmetric C-H stretching, asymmetric C-H stretching, and C=O stretching motions always enhance the H-atom dissociation probability obviously.

Origin of asynchronicity in Diels-Alder reactions

Asynchronicity in Diels-Alder reactions plays a crucial role in determining the height of the reaction barrier. Currently, the origin of asynchronicity is ascribed to the stronger orbital interaction between the diene and the terminal carbon of an asymmetric dienophile, which shortens the corresponding newly formed C-C bond and hence induces asynchronicity in the reaction. Here, we show, using the activation strain model and Kohn-Sham molecular orbital theory at ZORA-BP86/TZ2P, that this rationale behind asynchronicity is incorrect. We, in fact, found that following a more asynchronous reaction mode costs favorable HOMO-LUMO orbital overlap and, therefore, weakens (not strengthens) these orbital interactions. Instead, it is the Pauli repulsion that induces asynchronicity in Diels-Alder reactions. An asynchronous reaction pathway also lowers repulsive occupied-occupied orbital overlap which, therefore, reduces the unfavorable Pauli repulsion. As soon as this mechanism of reducing Pauli repulsion dominates, the reaction begins to deviate from synchronicity and adopts an asynchronous mode. The eventual degree of asynchronicity, as observed in the transition state of a Diels-Alder reaction, is ultimately achieved when the gain in stability, as a response to the reduced Pauli repulsion, balances with the loss of favorable orbital interactions.

Lewis Acid-Catalyzed Diels-Alder Reactions: Reactivity Trends across the Periodic Table

The catalytic effect of various weakly interacting Lewis acids (LAs) across the periodic table, based on hydrogen (Group 1), pnictogen (Group 15), chalcogen (Group 16), and halogen (Group 17) bonds, on the Diels-Alder cycloaddition reaction between 1,3-butadiene and methyl acrylate was studied quantum chemically by using relativistic density functional theory. Weakly interacting LAs accelerate the Diels-Alder reaction by lowering the reaction barrier up to 3 kcal mol-1 compared to the uncatalyzed reaction. The reaction barriers systematically increase from halogen Collapse

Key Words

• Diels-Alder reaction
• Lewis acids
• activation strain model
• density functional calculations
• reactivity

Collapse

MESH Headings Collapse

Grants

• Netherlands Organization for Scientific Research (NWO)
• Dutch Astrochemistry Network (DAN)
• PID2019-106184GB-I00 Spanish MICIIN
• RED2018-102387-T Spanish MICIIN
• Nederlandse Organisatie voor Wetenschappelijk Onderzoek
• Dutch Astrochemistry Network
• Nederlandse Organisatie voor Wetenschappelijk Onderzoek

Collapse

Affiliation(s)

Pascal Vermeeren Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands
Marco Dalla Tiezza Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands
Michelle van Dongen Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands
Israel Fernández Departamento de Química Orgánica I and Centro de Innovación en Química Avanzada (ORFEO-CINQA)Facultad de Ciencias QuímicasUniversidad Complutense de Madrid28040MadridSpain
F. Matthias Bickelhaupt Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands Institute for Molecules and MaterialsRadboud UniversityHeyendaalseweg 1356525 AJNijmegen (TheNetherlands
Trevor A. Hamlin Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdam (TheNetherlands

Collapse

Chemical reactivity from an activation strain perspective

Chemical reactions are ubiquitous in the universe, they are at the core of life, and they are essential for industrial processes. The drive for a deep understanding of how something occurs, in this case, the mechanism of a chemical reaction and the factors controlling its reactivity, is intrinsically valuable and an innate quality of humans. The level of insight and degree of understanding afforded by computational chemistry cannot be understated. The activation strain model is one of the most powerful tools in our arsenal to obtain unparalleled insight into reactivity. The relative energy of interacting reactants is evaluated along a reaction energy profile and related to the rigidity of the reactants' molecular structure and the strength of the stabilizing interactions between the deformed reactants: ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ). Owing to the connectedness between the activation strain model and Kohn-Sham molecular orbital theory, one is able to obtain a causal relationship between both the sterics and electronics of the reactants and their mutual reactivity. Only when this is accomplished one can eclipse the phenomenological explanations that are commonplace in the literature and textbooks and begin to rationally tune and optimize chemical transformations. We showcase how the activation strain model is the ideal tool to elucidate fundamental organic reactions, the activation of small molecules by metallylenes, and the cycloaddition reactivity of cyclic diene- and dipolarophiles.

Substituent effect on ESIPT mechanisms and photophysical properties of HBT derivatives

Substituent effects on excited-state intramolecular proton transfer (ESIPT) and photophysical properties of 2-(2-Hydroxyphenyl) benzothiazole (HBT) derivatives have been theoretically unveiled via the density functional theory (DFT) and time-dependent DFT (TDDFT). The optimized geometrical configurations and normal mode analyses confirm that the proton transfer processes are more reactive in excited state. Through calculating the activation energies and rate constants of ESIPT processes, finding that the processes are increasingly inactive when substituent group changes from -CN, -CO₂Me, -Cl, -Me, -NMe₂ to -NO₂. In addition, the photophysical properties analyses indicate the vertical transition energies are in good agreement with those observed in experiment. Note that all the absorption and emission maxima of enol and keto forms have the significant red-shift. In order to clarify the substituent effect on ESIPT and photophysical properties, we draw the frontier molecular orbitals (FMOs) isosurfaces and calculate the distances of electrons and holes and atomic charges. It follows that the intramolecular charge transfer (ICT) degrees are increasingly enlarged as substituting from -CN, -CO₂Me, -Cl, -Me, -NMe₂ to -NO₂ groups, which not only causes the red-shift of absorption and emission of enol and keto forms, but also affects the charge distribution of proton donor and acceptor, inhibiting the occurrence of ESIPT processes.

How Oriented External Electric Fields Modulate Reactivity

A judiciously oriented external electric field (OEEF) can catalyze a wide range of reactions and can even induce endo/exo stereoselectivity of cycloaddition reactions. The Diels-Alder reaction between cyclopentadiene and maleic anhydride is studied by using quantitative activation strain and Kohn-Sham molecular orbital theory to pinpoint the origin of these catalytic and stereoselective effects. Our quantitative model reveals that an OEEF along the reaction axis induces an enhanced electrostatic and orbital interaction between the reactants, which in turn lowers the reaction barrier. The stronger electrostatic interaction originates from an increased electron density difference between the reactants at the reactive center, and the enhanced orbital interaction arises from the promoted normal electron demand donor-acceptor interaction driven by the OEEF. An OEEF perpendicular to the plane of the reaction axis solely stabilizes the exo pathway of this reaction, whereas the endo pathway remains unaltered and efficiently steers the endo/exo stereoselectivity. The influence of the OEEF on the inverse electron demand Diels-Alder reaction is also investigated; unexpectedly, it inhibits the reaction, as the electric field now suppresses the critical inverse electron demand donor-acceptor interaction.

SN2 versus E2 Competition of F- and PH2- Revisited

We have quantum chemically analyzed the competition between the bimolecular nucleophilic substitution (S_N2) and base-induced elimination (E2) pathways for F^- + CH₃CH₂Cl and PH₂^- + CH₃CH₂Cl using the activation strain model and Kohn-Sham molecular orbital theory at ZORA-OLYP/QZ4P. Herein, we correct an earlier study that intuitively attributed the mechanistic preferences of F^- and PH₂^-, i.e., E2 and S_N2, respectively, to a supposedly unfavorable shift in the polarity of the abstracted β-proton along the PH₂^--induced E2 pathway while claiming that ″...no correlation between the thermodynamic basicity and E2 rate should be expected.″ Our analyses, however, unequivocally show that it is simply the 6 kcal mol^-1 higher proton affinity of F^- that enables this base to engage in a more stabilizing orbital interaction with CH₃CH₂Cl and hence to preferentially react via the E2 pathway, despite the higher characteristic distortivity (more destabilizing activation strain) associated with this pathway. On the other hand, the less basic PH₂^- has a weaker stabilizing interaction with CH₃CH₂Cl and is, therefore, unable to overcome the characteristic distortivity of the E2 pathway. Therefore, the mechanistic preference of PH₂^- is steered to the S_N2 reaction channel (less-destabilizing activation strain).

How Alkali Cations Catalyze Aromatic Diels-Alder Reactions

We have quantum chemically studied alkali cation-catalyzed aromatic Diels-Alder reactions between benzene and acetylene forming barrelene using relativistic, dispersion-corrected density functional theory. The alkali cation-catalyzed aromatic Diels-Alder reactions are accelerated by up to 5 orders of magnitude relative to the uncatalyzed reaction and the reaction barrier increases along the series Li+ < Na+ < K+ < Rb+ < Cs+ < none. Our detailed activation strain and molecular-orbital bonding analyses reveal that the alkali cations lower the aromatic Diels-Alder reaction barrier by reducing the Pauli repulsion between the closed-shell filled orbitals of the dienophile and the aromatic diene. We argue that such Pauli mechanism behind Lewis-acid catalysis is a more general phenomenon. Also, our results may be of direct importance for a more complete understanding of the network of competing mechanisms towards the formation of polycyclic aromatic hydrocarbons (PAHs) in an astrochemical context.

Sequential oxidations of phenylchalcogenides by H2O2: insights into the redox behavior of selenium via DFT analysis

The biological activity of sulfur and selenium, despite their similarity, shows some remarkable differences that have been recognized in many different scenarios.

PyFrag 2019-Automating the exploration and analysis of reaction mechanisms

We present a substantial update to the PyFrag 2008 program, which was originally designed to perform a fragment-based activation strain analysis along a provided potential energy surface. The original PyFrag 2008 workflow facilitated the characterization of reaction mechanisms in terms of the intrinsic properties, such as strain and interaction, of the reactants. The new PyFrag 2019 program has automated and reduced the time-consuming and laborious task of setting up, running, analyzing, and visualizing computational data from reaction mechanism studies to a single job. PyFrag 2019 resolves three main challenges associated with the automated computational exploration of reaction mechanisms: it (1) computes the reaction path by carrying out multiple parallel calculations using initial coordinates provided by the user; (2) monitors the entire workflow process; and (3) tabulates and visualizes the final data in a clear way. The activation strain and canonical energy decomposition results that are generated relate the characteristics of the reaction profile in terms of intrinsic properties (strain, interaction, orbital overlaps, orbital energies, populations) of the reactant species. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Dual Activation of Aromatic Diels-Alder Reactions

The unusually fast Diels-Alder reactions of [5]cyclophanes were analyzed by DFT at the BLYP-D3(BJ)/TZ2P level of theory. The computations were guided by an integrated activation-strain and Kohn-Sham molecular orbital analysis. It is revealed why both [5]metacyclophane and [5]paracyclophane exhibit a significant rate enhancement compared to their planar benzene analogue. The activation strain analyses revealed that the enhanced reactivity originates from 1) predistortion of the aromatic core resulting in a reduced activation strain of the aromatic diene, and/or 2) enhanced interaction with the dienophile through a distortion-controlled lowering of the HOMO-LUMO gap within the diene. Both of these physical mechanisms and thus the rate of Diels-Alder cycloaddition can be tuned through different modes of geometrical distortion (meta versus para bridging) and by heteroatom substitution in the aromatic ring. Judicious choice of the bridge and heteroatom in the aromatic core enables effective tuning of the aromatic Diels-Alder reactivity to achieve activation barriers as low as 2 kcal mol-1 , which is an impressive 35 kcal mol-1 lower than that of benzene.

Rocha MV, Hamlin TA, Poater J, Bickelhaupt FM. Understanding the differences between iron and palladium in cross-coupling reactions

We aim at developing design principles, based on quantum chemical analyses, for a novel type of iron-based catalysts that mimic the behavior of their well-known palladium analogs in the bond activation step of cross coupling reactions. To this end, we have systematically explored C-X bond activation via oxidative addition of CH3X substrates (X = H, Cl, CH3) to model catalysts mFe(CO)4q (q = 0, -2; m = singlet, triplet) and, for comparison, Pd(PH3)2 and Pd(CO)2, using relativistic density functional theory at the ZORA-OPBE/TZ2P level. We find that the neutral singlet iron catalyst 1Fe(CO)4 activates all three C-X bonds via barriers that are lower than those for Pd(PH3)2 and Pd(CO)2. This is a direct consequence of the capability of the iron complex to engage not only in π-backdonation, but also in comparably strong σ-donation. Interestingly, whereas the palladium complexes favor C-Cl activation, 1Fe(CO)4 shows a strong preference for activating the C-H bond, with a barrier as low as 10.4 kcal mol-1. Our results suggest a high potential for iron to feature in palladium-type cross-coupling reactions.

Arylic C-X Bond Activation by Palladium Catalysts: Activation Strain Analyses of Reactivity Trends

We have quantum chemically explored arylic carbon-substituent bond activation via oxidative insertion of a palladium catalyst in C₆H₅X + PdL_n model systems (X = H, Cl, CH₃; L_n = no ligand, PH₃, (PH₃)₂, PH₂C₂H₄PH₂) using relativistic density functional theory at ZORA-BLYP/TZ2P. Besides exploring reactivity trends and comparing them to aliphatic C-X activation, we aim at uncovering the physical factors behind the activity and selectivity. Our results show that barriers for arylic C-X activation are lower than those for the corresponding aliphatic C-X bonds. However, trends along bonds or upon variation of ligands are similar. Thus, bond activation barriers increase along C-Cl < C-H < C-C and along Pd < Pd(PH₃) or Pd(PH₂C₂H₄PH₂) < Pd(PH₃)₂. Activation strain analyses in conjunction with quantitative molecular orbital theory trace these trends to the rigidity and bonding capability of the various C-X bonds, model catalysts, and ligands.

Fluorine substituent effect on the stereochemistry of catalyzed and non-catalyzed Diels–Alder reactions. The case of R-butenone with cyclopentadiene: a computational assessment of the mechanism

A semiquantitative prediction of stereoselectivity due to substitutions of dienophile is obtained.

Reaction Coordinates and the Transition-Vector Approximation to the IRC

The appearance of a reaction profile or potential energy surface (PES) associated with the reaction path (defined as the path of steepest descent from the saddle point) depends on the choice of reaction coordinate onto which the intrinsic reaction coordinate is projected. This provides one with the freedom, but also the problem, of choosing the optimal perspective (i.e., the optimal reaction coordinate) for revealing what is essential for understanding the reaction. Here, we address this issue by analyzing a number of different reaction coordinates for the same set of model reactions, namely, prototypical oxidative addition reactions of C-X bonds to palladium. We show how different choices affect the appearance of the PES, and we discuss which qualities make a particular reaction coordinate most suitable for comparing and analyzing the reactions. Furthermore, we show how the transition vector (i.e., the normal mode associated with a negative force constant that leads from the saddle point to the steepest descent paths) can serve as a useful and computationally much more efficient approximation (designated TV-IRC) for full IRC computations, in the decisive region around the transition state.

Selective C-H and C-C Bond Activation: Electronic Regimes as a Tool for Designing d(10) MLn Catalysts

We wish to understand how a transition-metal catalyst can be rationally designed so as to selectively activate one particular bond in a substrate, herein, C-H and C-C bonds in ethane. To this end, we quantum chemically analyzed the activity and selectivity of a large series of model catalysts towards ethane and, for comparison, methane, by using the activation strain model and quantitative molecular orbital theory. The model catalysts comprise d(10) MLn complexes with coordination numbers n=0, 1, and 2; metal centers M=Co(-), Rh(-), Ir(-), Ni, Pd, Pt, Cu(+), Ag(+), and Au(+); and ligands L=NH3, PH3, and CO. Our analyses reveal that rather subtle electronic differences between bonds can be exploited to induce a lower barrier for activating one or the other, depending, among other factors, on the catalysts electronic regime (i.e., s-regime versus d-regime catalysts). Interestingly, the concepts and design principles emerging from this work can also be applied to the more challenging problem of differentiating between activation of the C-H bonds in ethane versus those in methane.