β-Lactone Synthesis from Epoxide and CO: Reaction Mechanism Revisited

Catalytic Ring-Opening Polymerisation of Cyclic Ethylene Carbonate: Importance of Elementary Steps for Determining Polymer Properties Revealed via DFT-MTD Simulations Validated Using Kinetic Measurements

The production of CO2-containing polymers is still very demanding in terms of controlling the synthesis of products with pre-defined CO2 content and molecular weight. An elegant way of synthesising these polymers is via CO2-containing building blocks, such as cyclic ethylene carbonate (cEC), via catalytic ring-opening polymerisation. However, to date, the mechanism of this reaction and control parameters have not been elucidated. In this work, using DFT-metadynamics simulations for exploiting the potential of the polymerisation process, we aim to shed more light on the mechanisms of the interaction between catalysts (in particular, the catalysts K3VO4, K3PO4, and Na2SnO3) and the cEC monomer in the propagation step of the polymeric chain and the occurring CO2 release. Confirming the simulation results via subsequent kinetics measurements indicates that, depending on the catalyst's characteristics, it can be attached reversibly to the polymeric chain during polymerisation, resulting in a defined lifetime of the activated polymer chain. The second anionic oxygen of the catalyst can promote the catalyst's transfer to another electrophilic cEC monomer, terminating the growth of the first chain and initiating the propagation of the new polymer chain. This transfer reaction is an essential step in controlling the molecular weight of the products.

A DFT-metadynamics study disclosing key properties of ring-opening polymerization catalysts to produce polyethercarbonate polyols from cyclic ethylene carbonate as part of an emerging CCU technology

The ring opening polymerization of cyclic carbonates made from epoxide and CO₂ to CO₂-containing polymers constitutes an emerging technology of particular industrial interest. Considering the reaction of ring-opening polymerization of cyclic ethylene carbonate to produce polyethercarbonate polyols, several types of catalysts were tested experimentally and mechanistic pathways were proposed, but a detailed analysis of structure property relationship including the CO₂-liberation pathways is still lacking. This contribution is using computational methods to investigate reported benchmark catalysts with the lead structure A_xM_yO_z (A: alkali metal or alkyl, M: main group element or transition metal) that are particularly approved as effiecient catalysts for industrial purpose. Employing DFT-metadynamics simulations, free energy surfaces (FESs) for the key-steps in the catalytic polymerization of cyclic ethylene carbonate (cEC) are generated. Important structural criteria and characteristics of the catalysts that influence the catalytic performance and (side)reaction pathways are determined. It turns out that less nucleophilicity of the catalyst anion and more labile cations remain major criteria for prohibiting CO₂ liberation during polymerization. The key learnings of this contribution currently serve as a basis to develop the next generation of catalysts to bring this emerging carbon capture and use (CCU) technology into industrial application.

Cobalt Carbonyls Stabilized by N,P-Ligands: Synthesis, Structure, and Catalytic Property for Ethylene Oxide Hydroalkoxycarbonylation

Reactions of N,P-Ligands as Ph₂ P(o-NMe₂ C₆ H₄ ) (¹ L), 2,6-iPr₂ C₆ H₃ NHC(Ph)=NC₆ H₄ (o-PPh₂ ) (² L), and Ph₂ PN(R)PPh₂ (R=iPr (³ L), cyclo-C₆ H₁₁ (⁴ L), tBu (⁵ L), CH₂ C₄ H₇ O (⁶ L)) each with dicobalt octacarbonyl produced complexes [¹ LCo(CO)₃ ]₂ (1), [² LCo(CO)(μ-CO)₂ Co(CO)₃ ] (2), [³ LCo(CO)₃ ]⁺ [Co(CO)₄ ]^- (3), [³ LCo(CO)₂ ]₂ (4), [⁴ LCo(CO)₂ ]₂ (5), [⁵ LCo(CO)₂ ]⁺ [Co(CO)₄ ]^- (6), and [⁶ LCo(CO)₂ ]⁺ [Co(CO)₄ ]^- (7). Complexes 1-7 have all been structurally characterized by X-ray crystallography, IR and NMR spectroscopies, and elemental analysis. Catalytic tests on transformation of ethylene oxide (EO), CO and MeOH into methyl 3-hydroxypropionate (3-HMP) indicate that complexes 1-7 are active, where ion-pair complexes 3 and 6-7 behave more excellently (by achieving 88.4-93.6% 3-HMP yields) than the neutral species 1-2 and 4-5 (35.0-46.5% 3-HMP yields) when the reactions are all operated at 2 MPa CO pressure and 50 °C in MeOH solvent. Density functional theory (DFT) study by selecting 3 as a model suggests a cooperative catalytic reaction mechanism by [Co(CO)₄ ]^- and its counter cation [³ LCo(CO)₃ ]⁺ . The cobalt-homonuclear ion-pair catalyzed hydroalkoxycarbonylation of EO is present herein.

Investigating the Transformations of Polyoxoanions Using Mass Spectrometry and Molecular Dynamics

The reactions of [γ-SiW10O36](8-) represent one of the most important synthetic gateways into a vast array of polyoxotungstate chemistry. Herein, we set about exploring the transformation of the lacunary polyoxoanion [β2-SiW11O39](8-) into [γ-SiW10O36](8-) using high-resolution electrospray mass spectrometry, density functional theory, and molecular dynamics. We show that the reaction proceeds through an unexpected {SiW9} precursor capable of undertaking a direct β → γ isomerization via a rotational transformation. The remarkably low-energy transition state of this transformation could be identified through theoretical calculations. Moreover, we explore the significant role of the countercations for the first time in such studies. This combination of experimental and the theoretical studies can now be used to understand the complex chemical transformations of oxoanions, leading to the design of reactivity by structural control.

Well-tempered metadynamics as a tool for characterizing multi-component, crystalline molecular machines

The well-tempered, smoothly converging form of the metadynamics algorithm has been implemented in classical molecular dynamics simulations and used to obtain an estimate of the free energy surface explored by the molecular rotations in the plastic crystal, octafluoronaphthalene. The biased simulations explore the full energy surface extremely efficiently, more than 4 orders of magnitude faster than unbiased molecular dynamics runs. The metadynamics collective variables used have also been expanded to include the simultaneous orientations of three neighboring octafluoronaphthalene molecules. Analysis of the resultant three-dimensional free energy surface, which is sampled to a very high degree despite its significant complexity, demonstrates that there are strong correlations between the molecular orientations. Although this correlated motion is of limited applicability in terms of exploiting dynamical motion in octafluoronaphthalene, the approach used is extremely well suited to the investigation of the function of crystalline molecular machines.

The nature of [PdCl(2)(C(2)H(4))(H(2)O)] as an active species in the Wacker process: new insights from ab initio molecular dynamics simulations

First-principles molecular dynamics coupled with metadynamics have been used to gain a deeper insight into the reaction mechanism of the Wacker process by determining the nature of the active species. An explicit and dynamic representation of the aqueous solvent, which was essential for modeling this reaction, was efficiently included into the simulations. Prompted by our earlier results, which showed that the configuration of the catalytically active species [PdCl(2)(H(2)O)(C(2)H(4))] was crucial in the subsequent steps of the Wacker process, herein we focused on the preceding equilibria that led to the formation of both the cis and trans isomers. Starting from the initial catalyst, [PdCl(4)](2-), the free-energy barriers for the forward and backward reactions were calculated. These results confirmed the relevance of the trans intermediate in the reaction mechanism, whilst conversely, they showed that the cis configuration played no role in it. This sole participation of the trans intermediate has some very important implications; besides the mechanistic interpretation of the initial steps in the Wacker reaction mechanism, the analysis of these equilibria provides additional information about the chemical nature of these ligand-substitution processes.

H2CO3 forms via HCO3- in water

According to the generally accepted picture of CO(2) dissolution in water, the formation of H(2)CO(3) proceeds in a single step that involves the attack of a water oxygen on the CO(2) carbon in concert with a proton transfer to a CO(2) oxygen. In the present work, a series of ab initio molecular dynamics simulations have been carried out along with the metadynamics technique which reveals a stepwise mechanism: the reaction of a water molecule with CO(2) yields HCO(3)(-) as an intermediate and a hydronium ion, whereas the protonation of the CO(2) moiety occurs in a separate step representing a well-defined activation barrier toward the H(2)CO(3) molecule. This alternative scenario was already taken into consideration decades ago, but subsequent experiments and calculations have given preference to the concerted mechanism. Employing extended periodic models of the CO(2)-water system that mimic the bulk aqueous environment, the present simulations yield the complete free energy profile of the stepwise mechanism and provide a detailed microscopic mechanism of the elementary steps. HCO(3)(-) formation is found to be the rate-determining step of the entire CO(2) hydration process.

The effects of water on beta-D-xylose condensation reactions

Car-Parrinello-based ab initio molecular dynamics simulations (CPMD) combined with metadynamics (MTD) simulations were used to determine the reaction energetics for the beta-D-xylose condensation reaction to form beta-1,4-linked xylobiose in a dilute acid solution. Protonation of the hydroxyl group on the xylose molecule and the subsequent breaking of the C-O bond were found to be the rate-limiting step during the xylose condensation reaction. Water and water structure was found to play a critical role in these reactions due to the proton's high affinity for water molecules. The reaction free energy and reaction barrier were determined using CPMD-MTD. We found that solvent reorganization due to proton partial desolvation must be taken into account in order to obtain the correct reaction activation energy. Our calculated reaction free energy and reaction activation energy compare well with available experimental results.

Towards a Rational Design of Ruthenium CO2 Hydrogenation Catalysts by Ab Initio Metadynamics

Complete reaction pathways relevant to CO2 hydrogenation by using a homogeneous ruthenium dihydride catalyst ([Ru(dmpe)2H2], dmpe=Me2PCH2CH2PMe2) have been investigated by ab initio metadynamics. This approach has allowed reaction intermediates to be identified and free-energy profiles to be calculated, which provide new insights into the experimentally observed reaction pathway. Our simulations indicate that CO2 insertion, which leads to the formation of formate complexes, proceeds by a concerted insertion mechanism. It is a rapid and direct process with a relatively low activation barrier, which is in agreement with experimental observations. Subsequent H2 insertion into the formate--Ru complex, which leads to the formation of formic acid, instead occurs via an intermediate [Ru(eta2-H2)] complex in which the molecular hydrogen coordinates to the ruthenium center and interacts weakly with the formate group. This step has been identified as the rate-limiting step. The reaction completes by hydrogen transfer from the [Ru(eta2-H2)] complex to the formate oxygen atom, which forms a dihydrogen-bonded Ru--HHO(CHO) complex. The activation energy for the H2 insertion step is lower for the trans isomer than for the cis isomer. A simple measure of the catalytic activity was proposed based on the structure of the transition state of the identified rate-limiting step. From this measure, the relationship between catalysts with different ligands and their experimental catalytic activities can be explained.

Ab initio molecular dynamics study of the keto-enol tautomerism of acetone in solution

We have studied the keto-enol interconversion of acetone to understand the mechanism of tautomerism relevant to numerous organic and biochemical processes. Applying the ab initio metadynamics method, we simulated the keto-enol isomerism both in the gas phase and in the presence of water. For the gas-phase intramolecular mechanism we show that no other hydrogen-transfer reactions can compete with the simple keto-enol tautomerism. We obtain an intermolecular mechanism and remarkable participation of water when acetone is solvated by neutral water. The simulations reveal that C deprotonation is the kinetic bottleneck of the keto-enol transformation, in agreement with experimental observations. The most interesting finding is the formation of short H-bonded chains of water molecules that provide the route for proton transfer from the carbon to the oxygen atom of acetone. The mechanistic picture that emerged from the present study involves proton migration and emphasizes the importance of active solvent participation in tautomeric interconversion.

Catalytic Double Carbonylation of Epoxides to Succinic Anhydrides: Catalyst Discovery, Reaction Scope, and Mechanism

The first catalytic method for the efficient conversion of epoxides to succinic anhydrides via one-pot double carbonylation is reported. This reaction occurs in two stages: first, the epoxide is carbonylated to a beta-lactone, and then the beta-lactone is subsequently carbonylated to a succinic anhydride. This reaction is made possible by the bimetallic catalyst [(ClTPP)Al(THF)2]+[Co(CO)4]- (1; ClTPP = meso-tetra(4-chlorophenyl)porphyrinato; THF = tetrahydrofuran), which is highly active and selective for both epoxide and lactone carbonylation, and by the identification of a solvent that facilitates both stages. The catalysis is compatible with substituted epoxides having aliphatic, aromatic, alkene, ether, ester, alcohol, nitrile, and amide functional groups. Disubstituted and enantiomerically pure anhydrides are synthesized from epoxides with excellent retention of stereochemical purity. The mechanism of epoxide double carbonylation with 1 was investigated by in situ IR spectroscopy, which reveals that the two carbonylation stages are sequential and non-overlapping, such that epoxide carbonylation goes to completion before any of the intermediate beta-lactone is consumed. The rates of both epoxide and lactone carbonylation are independent of carbon monoxide pressure and are first-order in the concentration of 1. The stages differ in that the rate of epoxide carbonylation is independent of substrate concentration and first-order in donor solvent, whereas the rate of lactone carbonylation is first-order in lactone and inversely dependent on the concentration of donor solvent. The opposite solvent effects and substrate order for these two stages are rationalized in terms of different resting states and rate-determining steps for each carbonylation reaction.

Carbonylation of heterocycles by homogeneous catalysts

This article summarizes the recent developments (particularly the uses of homogeneous organometallic catalysts) in ring-opening carbonylations, ring-opening carbonylative polymerizations and ring-expansion carbonylations of heterocycles such as epoxides, aziridines, lactones and oxazolines.

Ab Initio Exploration of Rearrangement Reactions: Intramolecular Hydrogen Scrambling Processes in Acetone

The recently developed metadynamics method is applied to the intramolecular hydrogen migration reactions of acetone in the gas phase. Comparison of different sets of collective coordinates allows efficient description of the underlying free energy surface. The simulations yielded numerous reactions: the enol-oxo tautomerism, the decomposition of acetone to various products, and rearrangement reactions. On the basis of the calculated activation barriers it is concluded that the enol-oxo tautomerism is the most frequent intramolecular proton-exchange process the acetone undergoes in the gas phase.

Ab initio Simulation of the Grafting of Phenylacetylene on Hydrogenated Surfaces of Crystalline Silicon Catalyzed by a Lewis Acid

Car-Parrinello simulations have been carried out to identify the grafting mechanism of phenylacetylene, a prototypical alkyne, on the hydrogenated surfaces of crystalline silicon, catalyzed by a Lewis acid (AlCl3). To this purpose, we have made use of a new technique, metadynamics, devised recently to deal with complex chemical reactions in first principles simulations. The reaction mechanism, leading to a styrenyl-terminated surface, turns out to be equivalent to the corresponding gas-phase hydrosilylation reaction by silanes that we have identified in a previous work. The activation energies for the surface reactions (0.43, 0.42, 0.35 eV, for H-Si(111), H-Si(100)2 x 1, and H-Si(100)1 x 1, respectively) are very close to that of the corresponding gas-phase reaction (0.37 eV). The estimated activation free energy at room temperature is sufficiently low for the grafting reaction to be viable at normal conditions and at low coverage on the crystalline silicon surfaces, as already well documented to occur on the surface of porous silicon. However, the conformation of the transition state shadows a large area of the surface, which might contribute to making the grafting process self-limiting.

The Mechanism of Epoxide Carbonylation by [Lewis Acid]+[Co(CO)4]- Catalysts

A detailed mechanistic investigation of epoxide carbonylation by the catalyst [(salph)Al(THF)2]+ [Co(CO)4]- (1, salph = N,N'-o-phenylenebis(3,5-di-tert-butylsalicylideneimine), THF = tetrahydrofuran) is reported. When the carbonylation of 1,2-epoxybutane (EB) to beta-valerolactone is performed in 1,2-dimethoxyethane solution, the reaction rate is independent of the epoxide concentration and the carbon monoxide pressure but first order in 1. The rate of lactone formation varies considerably in different solvents and depends primarily on the coordinating ability of the solvent. In mixtures of THF and cis/trans-2,5-dimethyltetrahydrofuran, the reaction is first order in THF. From spectroscopic and kinetic data, the catalyst resting state was assigned to be the neutral (beta-aluminoxy)acylcobalt species (salph)AlOCH(Et)CH2COCo(CO)4 (3a), which was successfully trapped with isocyanates. As the formation of 3a from EB, CO, and 1 is rapid, lactone ring closing is rate-determining. The favorable impact of donating solvents was attributed to the necessity of stabilizing the aluminum cation formed upon generation of the lactone.