New mechanism for autoxidation of polyolefins: kinetic Monte Carlo modelling of the role of short-chain branches, molecular oxygen and unsaturated moieties

Using quantum-chemical calculations and kinetic modelling, we reveal the dominant reaction pathway in polymer autoxidation, (i) whether oxygen promotes or hamper degradation and (ii) how defects in the polymer backbone affect degradation.

Computational prediction of the molecular configuration of three-dimensional network polymers

The three-dimensional arrangement of natural and synthetic network materials determines their application range. Control over the real-time incorporation of each building block and functional group is desired to regulate the macroscopic properties of the material from the molecular level onwards. Here we report an approach combining kinetic Monte Carlo and molecular dynamics simulations that chemically and physically predicts the interactions between building blocks in time and in space for the entire formation process of three-dimensional networks. This framework takes into account variations in inter- and intramolecular chemical reactivity, diffusivity, segmental compositions, branch/network point locations and defects. From the kinetic and three-dimensional structural information gathered, we construct structure-property relationships based on molecular descriptors such as pore size or dangling chain distribution and differentiate ideal from non-ideal structural elements. We validate such relationships by synthesizing organosilica, epoxy-amine and Diels-Alder networks with tailored properties and functions, further demonstrating the broad applicability of the platform.

Going Beyond the Carothers, Flory and Stockmayer Equation by Including Cyclization Reactions and Mobility Constraints

A challenge in the field of polymer network synthesis by a step-growth mechanism is the quantification of the relative importance of inter- vs. intramolecular reactions. Here we use a matrix-based kinetic Monte Carlo (kMC) framework to demonstrate that the variation of the chain length distribution and its averages (e.g., number average chain length x_n), are largely affected by intramolecular reactions, as mostly ignored in theoretical studies. We showcase that a conventional approach based on equations derived by Carothers, Flory and Stockmayer, assuming constant reactivities and ignoring intramolecular reactions, is very approximate, and the use of asymptotic limits is biased. Intramolecular reactions stretch the functional group (FG) conversion range and reduce the average chain lengths. In the likely case of restricted mobilities due to diffusional limitations because of a viscosity increase during polymerization, a complex x_n profile with possible plateau formation may arise. The joint consideration of stoichiometric and non-stoichiometric conditions allows the validation of hypotheses for both the intrinsic and apparent reactivities of inter- and intramolecular reactions. The kMC framework is also utilized for reverse engineering purposes, aiming at the identification of advanced (pseudo-)analytical equations, dimensionless numbers and mechanistic insights. We highlight that assuming average molecules by equally distributing A and B FGs is unsuited, and the number of AB intramolecular combinations is affected by the number of monomer units in the molecules, specifically at high FG conversions. In the absence of mobility constraints, dimensionless numbers can be considered to map the time variation of the fraction of intramolecular reactions, but still, a complex solution results, making a kMC approach overall most elegant.

Dehydration of butanol towards butenes over MFI, FAU and MOR: influence of zeolite topology

The effect of zeolite topology on the kinetics, selectivity and catalyst stability for the dehydration of butanol is studied experimentally and through microkinetic modeling.

Alumina-based Coating for Coke Reduction in Steam Crackers

Alumina-based coatings have been claimed as being an advantageous modification in industrial ethylene furnaces. In this work, on-line experimentally measured coking rates of a commercial coating (CoatAlloy™) have pointed out its superiority compared to an uncoated reference material in an electrobalance set-up. Additionally, the effects of presulfiding with 500 ppmw DMDS per H₂O, continuous addition of 41 ppmw S per HC of DMDS, and a combination thereof were evaluated during ethane steam cracking under industrially relevant conditions (T_gasphase = 1173 K, P_tot = 0.1 MPa, X_C2H6 = 70%, dilution δ = 0.33 kg_H2O/kg_HC). The examined samples were further evaluated using online thermogravimetry, scanning electron microscopy and energy diffractive X-ray for surface and cross-section analysis together with X-ray photoelectron spectroscopy and wavelength-dispersive X-ray spectroscopy for surface analysis. The passivating coating illustrated a better performance than the reference Ni-Cr Fe-base alloy after application of an improved pretreatment, followed by piddling changes on the product distribution. Presulfiding of the coating affected negatively the observed coking rates in comparison with the reference alloy, so alternative presulfiding and sulfur addition strategies are recommended when using this barrier coating.

Visualization and design of the functional group distribution during statistical copolymerization

Even though functional copolymers with a low percentage of functional comonomer units (up to 20 mol%) are widely used, for instance for the development of polymer therapeutics and hydrogels, insights in the functional group distribution over the actual chains are lacking and the average composition is conventionally used to describe the functionalization degree. Here we report the visualization of the monomer distribution over the different polymer chains by a synergetic combination of experimental and theoretical analysis aiming at the construction of functionality-chain length distributions (FUNC-CLDs). A successful design of the chemical structure of the comonomer pair, the initial functional comonomer amount (13 mol%), and the temperature (100 °C) is performed to tune the FUNC-CLD of copoly(2-oxazoline)s toward high functionalization degree for both low (100) and high (400) target degrees of polymerization. The proposed research strategy is generic and extendable to a broad range of copolymerization chemistries, including reversible deactivation radical polymerization.

Evaluation of a Ti-Base Alloy as Steam Cracking Reactor Material

Low-coking reactor material technologies are key for improving the performance and sustainability of steam crackers. In an attempt to appraise the coking performance of an alternative Ti-base alloy during ethane steam cracking, an experimental study was performed in a jet stirred reactor under industrially relevant conditions using thermogravimetry (T_gasphase = 1173 K, P_tot = 0.1 MPa, X_C2H6 = 70%, and dilution δ = 0.33 kg_H2O/kg_HC). Initially, a typical pretreatment used for Fe-Ni-Cr alloys was utilized and compared with a pretreatment at increased temperature, aiming at better surface oxidation and thus suppressing coke formation. The results revealed a decrease in coking rates upon high temperature pretreatment of the Ti-base alloy, however, its coking performance was significantly worse compared to the typically used Fe-Ni-Cr alloys, and carbon oxides formation increased by a factor of 30 or more. Moreover, the analyzed coupons showed crack propagation after coking/decoking and cooling down to ambient temperature. Scanning electron microscopy combined with energy-dispersive X-ray spectroscopy indicated that the prompt and unsystematic oxidation of the surface and bulk caused observable crack initiation and propagation due to alloy brittleness. Hence, the tested Ti-base alloy cannot be considered an industrially noteworthy steam cracking reactor alloy.

Modeling of Miniemulsion Polymerization of Styrene with Macro-RAFT Agents to Theoretically Compare Slow Fragmentation, Ideal Exchange and Cross-Termination Cases

A 5-dimensional Smith-Ewart based model is developed to understand differences for reversible addition-fragmentation chain transfer (RAFT) miniemulsion polymerization with theoretical agents mimicking cases of slow fragmentation, cross-termination, and ideal exchange while accounting for chain length and monomer conversion dependencies due to diffusional limitations. The focus is on styrene as a monomer, a water soluble initiator, and a macro-RAFT agent to avoid exit/entry of the RAFT leaving group radical. It is shown that with a too low RAFT fragmentation rate coefficient it is generally not afforded to consider zero-one kinetics (for the related intermediate radical type) and that with significant RAFT cross-termination the dead polymer product is dominantly originating from the RAFT intermediate radical. To allow the identification of the nature of the RAFT retardation it is recommended to experimentally investigate in the future the impact of the average particle size (dp) on both the monomer conversion profile and the average polymer properties for a sufficiently broad dp range, ideally including the bulk limit. With decreasing particle size both a slow RAFT fragmentation and a fast RAFT cross-termination result in a stronger segregation and thus rate acceleration. The particle size dependency is different, allowing further differentiation based on the variation of the dispersity and end-group functionality. Significant RAFT cross-termination is specifically associated with a strong dispersity increase at higher average particle sizes. Only with an ideal exchange it is afforded in the modeling to avoid the explicit calculation of the RAFT intermediate concentration evolution.

From n-butyl acrylate Arrhenius parameters for backbiting and tertiary propagation to β-scission via stepwise pulsed laser polymerization

A stepwise method to estimate the Arrhenius parameters for backbiting, tertiary propagation, and β-scission in acrylate radical polymerization.

Effect of Long-Term High Temperature Oxidation on the Coking Behavior of Ni-Cr Superalloys

The service time of an industrial cracker is strongly dependent on the long-term coking behavior and microstructure stability of the reactor coil alloy. Super alloys are known to withstand temperatures up to even 1400 K. In this work, several commercially available alloys have been first exposed to a long term oxidation at 1423 K for 500 h, so-called metallurgic aging. Subsequently, their coking behavior was evaluated in situ in a thermogravimetric setup under ethane steam cracking conditions (T_gasphase = 1173 K, P_tot = 0.1 MPa, X_C2H6 = 70%, continuous addition of 41 ppmw S/HC of DMDS, dilution δ = 0.33 kg_H2O/kg_HC) and compared with their unaged coking behavior. The tested samples were also examined using scanning electron microscopy and energy diffractive X-ray for surface and cross-section analysis. The alloys characterized by increased Cr-Ni content or the addition of Al showed improved stability against bulk oxidation and anti-coking behavior after application of metallurgic aging due to the formation of more stable oxides on the top surface.

Ab initio derived group additivity model for intramolecular hydrogen abstraction reactions

A set of group additivity values for intramolecular hydrogen abstraction reactions of alkanes, alkenes and alkynes is reported. Calculating 448 reaction rate coefficients at the CBS-QB3 level of theory for 1-2 up to 1-7 hydrogen shift reactions allowed the estimation of ΔGAV° values for 270 groups. The influence of substituents on (1) the attacking radical, (2) the attacked carbon atom, and (3) the carbon chain between the attacking and attacked reactive atom has been systematically studied. Substituents have been varied between hydrogen atoms and sp3, sp2 and sp hybridized carbon atoms. It has been assumed that substituents further away from the reactive atoms or their connecting carbon chain have negligible influences on the kinetics. This group additivity model is applicable to a wide variety of reactions in the 300-1800 K temperature range. Correlations for tunneling coefficients have been generated which are complementary to the ΔGAV°'s to obtain accurate rate coefficients without the need for imaginary frequencies or electronic energies of activation. These correlations depend on the temperature and activation energy of the exothermic step. The group additivity model has been successfully applied to a test set of reactions also calculated at the CBS-QB3 level of theory. A mean absolute deviation of 1.18 to 1.71 has been achieved showing a good overall accuracy of the model.

Access to the β-scission rate coefficient in acrylate radical polymerization by careful scanning of pulse laser frequencies at elevated temperature

A novel method to estimate the β-scission rate coefficient in acrylate radical polymerization is presented.

Ab-Initio-Based Kinetic Modeling to Understand RAFT Exchange: The Case of 2-Cyano-2-Propyl Dodecyl Trithiocarbonate and Styrene

Ab-initio-calculated rate coefficients for addition and fragmentation in reversible-addition fragmentation chain transfer (RAFT) polymerization of styrene with 2-cyano-2-propyl dodecyl trithiocarbonate initiated by azobisisobutyronitrile allow the reliable simulation of the experimentally observed conversion, number average chain length, and dispersity. The rate coefficient for addition of a macroradical R_i to the macroRAFT agent R_i X at 333 K (6.8 10⁴ L mol^-1 s^-1 ) is significantly lower than to the initial RAFT agent R₀ X (3.2 10⁶ L mol^-1 s^-1 ), mainly due to a difference in activation energy (15.4 vs 3.0 kJ mol^-1 ), which causes the dispersity to spike in the beginning of the polymerization.

Mechanistic insights into the formation of butene isomers from 1-butanol in H-ZSM-5: DFT based microkinetic modelling

First principles microkinetic modelling provides in-depth mechanistic insights into the competing reaction pathways for zeolite-catalyzed conversion of 1-butanol to butene isomers.

Ab initio based kinetic Monte Carlo analysis to unravel the propagation kinetics in vinyl acetate pulsed laser polymerization

The radical propagation kinetics of vinyl acetate in pulsed laser polymerization (PLP) is studied by combining ab initio calculated rate coefficients with kinetic Monte Carlo modeling of PLP spectra.

A detailed mechanistic study of bulk MADIX of styrene and its chain extension

By combining experimental and modeling tools, a detailed characterization study of MADIX properties becomes possible.

Thiol-Michael addition in polar aprotic solvents: nucleophilic initiation or base catalysis?

The thiol-Michael addition of ethanethiol to ethyl acrylate, methyl vinylsulfone and maleimide initiated by ethyl-, diethyl-, triethylamine and triethylphosphine in tetrahydrofuran (THF) is investigated at room temperature.

Ab initio coverage-dependent microkinetic modeling of benzene hydrogenation on Pd(111)

Coverage-dependent calculations are required for an accurate DFT-based prediction of the activity and a correct mechanistic understanding of catalytic hydrogenation.

Effect of zeolite confinement on the conversion of 1-butanol to butene isomers: mechanistic insights from DFT based microkinetic modelling

First principles microkinetic modelling shows that, unlike in H-ZSM-5 and H-ZSM-22, trans-2-butene formation in H-FER occurs via direct dehydration of 1-butanol.

Quantitative First-Principles Kinetic Modeling of the Aza-Michael Addition to Acrylates in Polar Aprotic Solvents

This work presents a detailed computational study and kinetic analysis of the aza-Michael addition of primary and secondary amines to acrylates in an aprotic solvent. Accurate rate coefficients for all elementary steps in the various competing mechanisms are calculated using an ONIOM-based approach in which the full system is calculated with M06-2X/6-311+G(d,p) and the core system with CBS-QB3 corrected for solvation using COSMO-RS. Diffusional contributions are taken into account using the coupled encounter pair model with diffusion coefficients calculated based on molecular dynamics simulations. The calculated thermodynamic and kinetic parameters for all forward and reverse elementary reactions are fed to a microkinetic model giving excellent agreement with experimental data obtained using GC analysis. Rate analysis reveals that for primary and secondary amines, the aza-Michael addition to ethyl acrylate occurs preferentially according to a 1,2-addition mechanism, consisting of the pseudoequilibrated formation of a zwitterion followed by a rate controlling amine assisted proton transfer toward the singly substituted product. The alternative 1,4-addition becomes competitive if substituents are present on the amine or double bond of the acrylate. Primary amines react faster than secondary amines due to increased solvation of the zwitterionic intermediate and less sterically hindered proton transfer.

Computational Investigation of the Aminolysis of RAFT Macromolecules

This work presents a detailed computational study and kinetic analysis of the aminolysis of dithioates, dithiobenzoates, trithiocarbonates, xanthates, and thiocarbamates, which are frequently used as chain-transfer agents for reversible addition-fragmentation chain-transfer (RAFT) polymerization. Rate coefficients were obtained from ab initio calculations, taking into account a diffusional contribution according to the encounter pair model. A kinetic model was constructed and reveals a reaction mechanism of four elementary steps: (i) formation of a zwitterionic intermediate, (ii) formation of a complex intermediate in which an assisting amine molecule takes over the proton from the zwitterionic intermediate, (iii) breakdown of the complex into a neutral tetrahedral intermediate with release of the assisting amine molecule, and (iv) amine-assisted breakdown of the neutral intermediate to the products. Furthermore, a comparative analysis indicates that the alkanedithioates and dithiobenzoates react the fastest, followed, respectively, by xanthates and trithiocarbonates, which react almost equally fast, and dithiocarbamates, which are not reactive at typical experimentally used conditions.