Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

A master equation model for bimolecular reaction via multi-well isomerization intermediates

A reversible linear master equation model is presented for pressure- and temperature-dependent bimolecular reactions proceeding via multiple long-lived intermediates. This kinetic treatment, which applies when the reactions are measured under pseudo-first-order conditions, facilitates accurate and efficient simulation of the time dependence of the populations of reactants, intermediate species and products. Detailed exploratory calculations have been carried out to demonstrate the capabilities of the approach, with applications to the bimolecular association reaction C3H6 + H reversible arrow C3H7 and the bimolecular chemical activation reaction C2H2 +(CH2)-C-1--> C3H3+H. The efficiency of the method can be dramatically enhanced through use of a diffusion approximation to the master equation, and a methodology for exploiting the sparse structure of the resulting rate matrix is established.

Prediction of absolute rate coefficients and product branching ratios for the C(P-3) plus allene reaction system

Complex chemical reactions in the gas phase can be decomposed into a network of elementary (e.g., unimolecular and bimolecular) steps which may involve multiple reactant channels, multiple intermediates, and multiple products. The modeling of such reactions involves describing the molecular species and their transformation by reaction at a detailed level. Here we focus on a detailed modeling of the C(P-3)+allene (C3H4) reaction, for which molecular beam experiments and theoretical calculations have previously been performed. In our previous calculations, product branching ratios for a nonrotating isomerizing unimolecular system were predicted. We extend the previous calculations to predict absolute unimolecular rate coefficients and branching ratios using microcanonical variational transition state theory (mu-VTST) with full energy and angular momentum resolution. Our calculation of the initial capture rate is facilitated by systematic ab initio potential energy surface calculations that describe the interaction potential between carbon and allene as a function of the angle of attack. Furthermore, the chemical kinetic scheme is enhanced to explicitly treat the entrance channels in terms of a predicted overall input flux and also to allow for the possibility of redissociation via the entrance channels. Thus, the computation of total bimolecular reaction rates and partial capture rates is now possible. (C) 2002 American Institute of Physics.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

Dependência das vias de dissociação unimolecular com a energia interna da molécula de éter dimetílico

Internal energy dependence of the competitive unimolecular dissociation channels of dimethyl ether were studied with the statistical RRKM formalism. The C-O and C-H fission reactions and the 1,2-H and 1,3-H shifts, and 1,1-H2 and 1,3-H2 molecular eliminations are discussed as a function of energy dependence of k a(E*), the microcanonical rate constant for production of transition states. C-O fission is the dominant process while reaction channels involving C-H fission, 1,1-H2 and 1,3-H2 elimination and production of MeOH should be competitive at energies around 400 kJ mol-1. The less favorable process is the channel of CH4 formation.

Novel Reactions of TRIS(2-Thienyl)Methane and TRIS(2-Furyal)Methane

A systematic investigation of the reactivity and functionalization of two heterocyclic analogs of triphenylmethane , namely tris(2-thienyl)methane and tris(2-furyl)methane have been carried out and the results are presented in this thesis entitled "NOVEL REACTIONS OF TRIS(2-THIENYL)METHANE AND TRIS(2-FURYL)METHANE.". The history of organic free radicals dates back to Gomberg's monumental discovery of the triphenylmethyl radical in 1900. The heterocyclic analogs of triarylmethane are also interesting from the vantage point of their transformation to the corresponding radicals akin to Gomberg ' s triphenylmethyl radical and also they are prone to further transformation leading to three dimensionally elongated molecules such as dendrimers. Dendritic architectures are one of the most pervasive topologies observed in nature at the macro- and microdimensional length devices. Because of their ability to combine both organic and inorganic compounds and their propensity to either encapsulate or be engineered into unimolecular functional devices , dendrimers are versatile amongst existing nanoscale building blocks and materials.

Structure, energetics and reactions of metal cation complexes of dipeptides in the gas phase

Structure, energetics and reactions of ions in the gas phase can be revealed by mass spectrometry techniques coupled to ions activation methods. Ions can gain enough energy for dissociation by absorbing IR light photons introduced by an IR laser to the mass spectrometer. Also collisions with a neutral molecule can increase the internal energy of ions and provide the dissociation threshold energy. Infrared multiple photon dissociation (IRMPD) or sustained off-resonance irradiation collision-induced dissociation (SORI-CID) methods are combined with Fourier Transform Ion Cyclotron Resonance (FT-ICR) mass spectrometers where ions can be held at low pressures for a long time. The outcome of ion activation techniques especially when it is compared to the computational methods results is of great importance since it provides useful information about the structure, thermochemistry and reactivity of ions of interest. In this work structure, energetics and reactivity of metal cation complexes with dipeptides are investigated. Effect of metal cation size and charge as well as microsolvation on the structure of these complexes has been studied. Structures of bare and hydrated Na and Ca complexes with isomeric dipeptides AlaGly and GlyAla are characterized by means of IRMPD spectroscopy and computational methods. At the second step unimolecular dissociation reactions of singly charged and doubly charged multimetallic complexes of alkaline earth metal cations with GlyGly are examined by CID method. Also structural features of these complexes are revealed by comparing their IRMPD spectra with calculated IR spectra of possible structures. At last the unimolecular dissociation reactions of Mn complexes are studied. IRMPD spectroscopy along with computational methods is also employed for structural elucidation of Mn complexes. In addition the ion-molecule reactions of Mn complexes with CO and water are explored in the low pressures obtained in the ICR cell.

Fine Needle Aspiration Cytology In The Diagnosis Of Perioral Adverse Reactions To Cosmetic Dermal Fillers.

Facial cosmetic procedures are increasingly requested, and dermal filler materials have been widely used as a nonsurgical option since the 1980s. However, injectable fillers have been implicated in local adverse reactions. Therefore, the aim of this article was to describe the use of fine needle aspiration cytology (FNAC) in the diagnosis of foreign-body reactions to the perioral injection of dermal fillers. A 69-year-old woman presented with a painful nodule on her right nasolabial fold. Intraoral FNAC was performed, and cytologic smears were examined under optical and polarized light microscopy, showing birefringent microspheres, confirming the diagnosis of an adverse reaction caused by polymethyl methacrylate filler. FNAC is a less invasive method to confirm the diagnosis of adverse reactions caused by perioral cosmetic dermal fillers.

Impact Of Drug Formulation And Free Platinum/cisplatin Ratio On Hypersensitivity Reactions To Cisplatin: Formulation Matters.

Use of cisplatin can induce type I hypersensitivity reactions that may also be linked to the quality of the drug utilized. We observed cases of hypersensitivity that appeared to be associated with the brand of cisplatin used. The aim of this study was to compare two different brands of cisplatin in relation to type I hypersensitivity reactions. Brand A was used in a tertiary care teaching hospital until 2012, and use of brand B started from January 2013, when the first hypersensitivity cases were observed. Patients were categorized based on symptom. Cisplatin of both brands was analysed by high-performance liquid chromatography (HPLC) and high-resolution electrospray ionization mass spectrometry (ESI-(+)-MS) and characterized according to US Pharmacopeia. There were no cases of hypersensitivity associated with the use of cisplatin brand A, whereas four of 127 outpatients that used cisplatin brand B were affected. The two brands were in accordance with the US Pharmacopeia parameters, and there was no significant difference in the total platinum levels between the two brands when analysed by HPLC. However, high-resolution ESI-(+)-MS analyses show that brand B contains approximately 2.7 times more hydrolysed cisplatin than brand A. The increase in the hydrolysed form of cisplatin found in brand B may be the cause of the hypersensitivity reaction observed in a subset of patients. We present the first study of the quality of drugs by high-resolution ESI-(+)-MS. Drug regulatory agencies and manufacturers should consider including measurement of hydrolysed cisplatin as a quality criterion for cisplatin formulations.

Gas-phase solvolysis type reactions of SiCl3+ cations

Gas-phase SiCl3+ ions undergo sequential solvolysis type reactions with water, methanol, ammonia, methylamine and propylene. Studies carried out in a Fourier Transform mass spectrometer reveal that these reactions are facile at 10-8 Torr and give rise to substituted chlorosilyl cations. Ab initio and DFT calculations reveal that these reactions proceed by addition of the silyl cation to the oxygen or nitrogen lone pair followed by a 1,3-H migration in the transition state. These transition states are calculated to lie below the energy of the reactants. By comparison, hydrolysis of gaseous CCl3+ is calculated to involve a substantial positive energy barrier.

Electrochemical and photocatalytic reactions of polycyclic aromatic hydrocarbons investigated by raman spectroscopy

This paper presents the study of photochemical behavior of polycyclic aromatic hydrocarbons (PAHs), potential pollutants in secondary reactions in aerosols, through Raman spectroscopy compared with its electrochemical behavior. The PAHs studied include pyrene, anthracene, phenanthrene and fluorene. These were adsorbed onto TiO2 and irradiated with ultraviolet light (254 nm). Their electrochemical oxidation was studied by in situ Surface-enhanced Raman Scattering (SERS) and led to the formation of carbonyl-containing products. Oxidized intermediates bearing the C=O group were also formed during photodegradation. The joint analysis of the photodegradation data with those produced by electrochemical means - using spectroscopic techniques for the identification and characterization of the products - revealed the formation of identical products for anthracene, but not for pyrene. A reasonable explanation for this difference in results is that photochemical and electrochemical oxidation reactions proceed via different mechanisms. While photocatalytic degradation over TiO2 is initiated by hydroxyl radicals, electrochemical oxidation is initiated by the direct electron transfer from adsorbed PAH to the electrode, generating PAH cation radicals that undergo subsequent reactions.

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties. (c) 2008 American Institute of Physics.