Modeling the thermochemistry of nitrogen-containing compounds group additivity.

First-principles based kinetic modeling is essential to gain insight into the governing chemistry of nitrogen-containing compounds over a wide range of technologically important processes, pyrolysis, oxidation and combustion. It also enables the development of predictive, fundamental models key to improving understanding of the influence of nitrogen-containing compounds present as impurities or process additives, considering safety, operability and quality of the product streams. A prerequisite for the generation of detailed fundamental kinetic models is the availability of accurate thermodynamic properties.

Modeling the kinetics of hydrogen abstraction reactions in nitrogen-containing compounds group additivity.

New group additivity values are presented to enable the modeling of a broad range of intermolecular hydrogen abstraction reactions involving nitrogen-containing compounds. From a dataset of 316 reaction rate coefficients calculated at the CBS-QB3 level of theory in the high-pressure limit, 76 group additivity values and 14 resonance corrections have been estimated. The influence of substituents on both the attacked hydrogen and attacking radical, being a carbon or nitrogen atom, has been investigated systematically.

Effect of Phosphine on Coke Formation during Steam Cracking of Propane.

In conventional steam cracking feedstocks, contaminants such as sulfur, phosphine, and heavy metal components, present in trace levels, are believed to affect coke formation on high temperature alloys. To gain an understanding of the role of phosphine coking rates on 25/35, CrNi and Al-containing reactor materials were determined in a plug flow reactor during cracking of a propane feedstock doped with ppb levels of PH in the presence of DMDS. The presence of phosphine decreased the asymptotic coking rates by more than 20%, while it had a smaller influence on the catalytic coking rate.

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 (MC) framework to demonstrate that the variation of the chain length distribution and its averages (e.

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.

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 HO, 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 = 1173 K, P = 0.

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.

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 = 1173 K, P = 0.1 MPa, X = 70%, and dilution δ = 0.

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.

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.

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.

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 to the macroRAFT agent R X at 333 K (6.8 10 L mol s ) is significantly lower than to the initial RAFT agent R X (3.

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.

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.

Insights into the Reaction Mechanism of Ethanol Conversion into Hydrocarbons on H-ZSM-5.

Ethanol dehydration to ethene is mechanistically decoupled from the production of higher hydrocarbons due to complete surface coverage by adsorbed ethanol and diethyl ether (DEE). The production of C3+ hydrocarbons was found to be unaffected by water present in the reaction mixture. Three routes for the production of C3+ hydrocarbons are identified: the dimerization of ethene to butene and two routes involving two different types of surface species categorized as aliphatic and aromatic.

Improved Livingness and Control over Branching in RAFT Polymerization of Acrylates: Could Microflow Synthesis Make the Difference?

The superior capabilities of structured microreactors over batch reactors are demonstrated for reversible addition-fragmentation chain transfer (RAFT) solution polymerization of n-butyl acrylate with the aid of simulations, explicitly accounting for the chain length distribution of all macrospecies types. Since perfect isothermicity can be established in a microreactor, less side products due to backbiting and β-scission are formed compared to the batch operation in which ineffective heat removal leads to an undesirable temperature spike. For a given RAFT chain transfer agent (CTA), additional microstructural control results under microflow conditions by optimizing the reaction temperature, lowering the dilution degree, or decreasing the initial molar ratio of monomer to RAFT CTA.

Computational Study and Kinetic Analysis of the Aminolysis of Thiolactones.

The aminolysis of three differently α-substituted γ-thiolactones (C4H5OSX, X = H, NH2, and NH(CO)CH3) is modeled based on CBS-QB3 calculated free energies corrected for solvation using COSMO-RS. For the first time, quantitative kinetic and thermodynamic data are provided for the concerted path and the stepwise path over a neutral tetrahedral intermediate. These paths can take place via an unassisted, an amine-assisted, or a thiol-assisted mechanism.

Group Additive Kinetics for Hydrogen Transfer Between Oxygenates.

Hydrogen abstraction reactions involving oxygenates in gaseous phase play an important role in many biomass-related conversion processes. In this work, group additivity is used to provide Arrhenius parameters in a temperature range of 300-2500 K for hydrogen abstractions between oxygenate compounds such as alcohols, ethers, esters, acids, ketones, diketones, aldehydes, hydroxyperoxides, alkyl peroxides, and unsaturated ethers and ketones. The group additive values for Arrhenius parameters of hydrogen transfer reactions of the type O--H--C and O--H--O are derived from CBS-QB3 calculations in the high-pressure limit.

Symmetry calculation for molecules and transition states.

The symmetry of molecules and transition states of elementary reactions is an essential property with important implications for computational chemistry. The automated identification of symmetry by computers is a very useful tool for many applications, but often relies on the availability of three-dimensional coordinates of the atoms in the molecule and hence becomes less useful when these coordinates are a priori unavailable. This article presents a new algorithm that identifies symmetry of molecules and transition states based on an augmented graph representation of the corresponding structures, in which both topology and the presence of stereocenters are accounted for.

Thermodynamic study of benzene and hydrogen coadsorption on Pd(111).

Periodic density functional theory (DFT) has been used to study the coadsorption of hydrogen and benzene on Pd(111). The most stable coverages are predicted by constructing the thermodynamic phase diagram as a function of gas-phase temperature and pressure. The common approximation that neglects vibrational contributions to the surface Gibbs free energy, using the PW91 functional, is compared to the one that includes vibrational contributions.

Kinetic modeling of α-hydrogen abstractions from unsaturated and saturated oxygenate compounds by hydrogen atoms.

Hydrogen-abstraction reactions play a significant role in thermal biomass conversion processes, as well as regular gasification, pyrolysis, or combustion. In this work, a group additivity model is constructed that allows prediction of reaction rates and Arrhenius parameters of hydrogen abstractions by hydrogen atoms from alcohols, ethers, esters, peroxides, ketones, aldehydes, acids, and diketones in a broad temperature range (300-2000 K). A training set of 60 reactions was developed with rate coefficients and Arrhenius parameters calculated by the CBS-QB3 method in the high-pressure limit with tunneling corrections using Eckart tunneling coefficients.

Experimental and theoretical methods in kinetic studies of heterogeneously catalyzed reactions.

This review aims to illustrate the potential of kinetic analysis in general and microkinetic modeling in particular for rational catalyst design. Both ab initio calculations and experiments providing intrinsic kinetic data allow us to assess the effects of catalytic properties and reaction conditions on the activity and selectivity of the targeted reactions. Three complementary approaches for kinetic analysis of complex reaction networks are illustrated, using select examples of acid zeolite-catalyzed reactions from the authors' recent work.

Kinetic modeling of α-hydrogen abstractions from unsaturated and saturated oxygenate compounds by carbon-centered radicals.

Hydrogen abstractions are important elementary reactions in a variety of reacting media at high temperatures in which oxygenates and hydrocarbon radicals are present. Accurate kinetic data are obtained from CBS-QB3 ab initio (AI) calculations by using conventional transition-state theory within the high-pressure limit, including corrections for hindered rotation and tunneling. From the obtained results, a group-additive (GA) model is developed that allows the Arrhenius parameters and rate coefficients for abstraction of the α-hydrogen from a wide range of oxygenate compounds to be predicted at temperatures ranging from 300 to 1500 K.

Kinetic modeling of hydrogen abstractions involving sulfur radicals.

One of the requisites for the development of detailed reaction networks is the availability of accurate kinetic data. Group additivity based models linking the Arrhenius parameters to structural characteristics of the transition state have proven to be a valuable tool to obtain those data. In this work, group additivity values are presented to allow a broad range of CH and SH hydrogen abstraction reactions by S radicals to be modeled.

Group additive values for the gas-phase standard enthalpy of formation, entropy and heat capacity of oxygenates.

A complete and consistent set of 60 Benson group additive values (GAVs) for oxygenate molecules and 97 GAVs for oxygenate radicals is provided, which allow to describe their standard enthalpies of formation, entropies and heat capacities. Approximately half of the GAVs for oxygenate molecules and the majority of the GAVs for oxygenate radicals have not been reported before. The values are derived from an extensive and accurate database of thermochemical data obtained by ab initio calculations at the CBS-QB3 level of theory for 202 molecules and 248 radicals.