Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.

Functionalised zinc oxide nanowire gas sensors : enhanced NO2 gas sensor response by chemical modification of nanowire surfaces

Surface coating with an organic self-assembled monolayer (SAM) can enhance surface reactions or the absorption of specific gases and hence improve the response of a metal oxide (MOx) sensor toward particular target gases in the environment. In this study the effect of an adsorbed organic layer on the dynamic response of zinc oxide nanowire gas sensors was investigated. The effect of ZnO surface functionalisation by two different organic molecules, tris(hydroxymethyl)aminomethane (THMA) and dodecanethiol (DT), was studied. The response towards ammonia, nitrous oxide and nitrogen dioxide was investigated for three sensor configurations, namely pure ZnO nanowires, organic-coated ZnO nanowires and ZnO nanowires covered with a sparse layer of organic-coated ZnO nanoparticles. Exposure of the nanowire sensors to the oxidising gas NO2 produced a significant and reproducible response. ZnO and THMA-coated ZnO nanowire sensors both readily detected NO2 down to a concentration in the very low ppm range. Notably, the THMA-coated nanowires consistently displayed a small, enhanced response to NO2 compared to uncoated ZnO nanowire sensors. At the lower concentration levels tested, ZnO nanowire sensors that were coated with THMA-capped ZnO nanoparticles were found to exhibit the greatest enhanced response. ΔR/R was two times greater than that for the as-prepared ZnO nanowire sensors. It is proposed that the ΔR/R enhancement in this case originates from the changes induced in the depletion-layer width of the ZnO nanoparticles that bridge ZnO nanowires resulting from THMA ligand binding to the surface of the particle coating. The heightened response and selectivity to the NO2 target are positive results arising from the coating of these ZnO nanowire sensors with organic-SAM-functionalised ZnO nanoparticles.

Chemical variability of groundwater samples collected from a Coal Seam Gas exploration well, Maramarua, New Zealand

A pilot study has produced 31 groundwater samples from a coal seam gas (CSG) exploration well located in Maramarua, New Zealand. This paper describes sources of CSG water chemistry variations, and makes sampling and analytical recommendations to minimize these variations. The hydrochemical character of these samples is studied using factor analysis, geochemical modelling, and a sparging experiment. Factor analysis unveils carbon dioxide (CO2) degassing as the principal cause of sample variation (about 33%). Geochemical modelling corroborates these results and identifies minor precipitation of carbonate minerals with degassing. The sparging experiment confirms the effect of CO2 degassing by showing a steady rise in pH while maintaining constant alkalinity. Factor analysis correlates variations in the major ion composition (about 17%) to changes in the pumping regime and to aquifer chemistry variations due to cation exchange reactions with argillaceous minerals. An effective CSG water sampling program can be put into practice by measuring pH at the well head and alkalinity at the laboratory; these data can later be used to calculate the carbonate speciation at the time the sample was collected. In addition, TDS variations can be reduced considerably if a correct drying temperature of 180°C is consistently implemented.

Deformation with coupled chemical diffusion

The deformation of rocks is commonly intimately associated with metamorphic reactions. This paper is a step towards understanding the behaviour of fully coupled, deforming, chemically reacting systems by considering a simple example of the problem comprising a single layer system with elastic-power law viscous constitutive behaviour where the deformation is controlled by the diffusion of a single chemical component that is produced during a metamorphic reaction. Analysis of the problem using the principles of non-equilibrium thermodynamics allows the energy dissipated by the chemical reaction-diffusion processes to be coupled with the energy dissipated during deformation of the layers. This leads to strain-rate softening behaviour and the resultant development of localised deformation which in turn nucleates buckles in the layer. All such diffusion processes, in leading to Herring-Nabarro, Coble or “pressure solution” behaviour, are capable of producing mechanical weakening through the development of a “chemical viscosity”, with the potential for instability in the deformation. For geologically realistic strain rates these chemical feed-back instabilities occur at the centimetre to micron scales, and so produce structures at these scales, as opposed to thermal feed-back instabilities that become important at the 100–1000 m scales.

ortho-Bromo(propa-1,2-dien-1-yl)arenes : substrates for domino reactions

o-Bromo(propa-1,2-dien-1-yl)arenes exhibit novel and orthogonal reactivity under Pd catalysis in the presence of secondary amines to form enamines (concerted Pd insertion, intramolecular carbopalladation, and terminative Buchwald–Hartwig coupling) and of amides to form indoles (addition, Buchwald–Hartwig cyclization, and loss of the acetyl group). The substrates for these reactions can be accessed in a reliable and highly selective two-step process from 2-bromoaryl bromides.

Reactions of simple and peptidic alpha-carboxylate radical anions with dioxygen in the gas phase

alpha-Carboxylate radical anions are potential reactive intermediates in the free radical oxidation of biological molecules (e. g., fatty acids, peptides and proteins). We have synthesised well-defined alpha-carboxylate radical anions in the gas phase by UV laser photolysis of halogenated precursors in an ion-trap mass spectrometer. Reactions of isolated acetate ((center dot)CH(2)CO(2)) and 1-carboxylatobutyl (CH(3)CH(2)CH(2)(center dot)CHCO(2)(-)) radical anions with dioxygen yield carbonate (CO(3)(center dot-)) radical anions and this chemistry is shown to be a hallmark of oxidation in simple and alkyl-substituted cross-conjugated species. Previous solution phase studies have shown that C(alpha)-radicals in peptides, formed from free radical damage, combine with dioxygen to form peroxyl radicals that subsequently decompose into imine and keto acid products. Here, we demonstrate that a novel alternative pathway exists for two alpha-carboxylate C(alpha)-radical anions: the acetylglycinate radical anion (CH(3)C(O)NH(center dot)CHCO(2)(-)) and the model peptide radical anion, YGGFG(center dot-). Reaction of these radical anions with dioxygen results in concerted loss of carbon dioxide and hydroxyl radical. The reaction of the acetylglycinate radical anion with dioxygen reveals a two-stage process involving a slow, followed by a fast kinetic regime. Computational modelling suggests the reversible formation of the C(alpha) peroxyl radical facilitates proton transfer from the amide to the carboxylate group, a process reminiscent of, but distinctive from, classical proton-transfer catalysis. Interestingly, inclusion of this isomerization step in the RRKM/ME modelling of a G3SX level potential energy surface enables recapitulation of the experimentally observed two-stage kinetics.

Gas-phase reactions of aryl radicals with 2-butyne : an experimental and theoretical investigation employing the N-methyl-pyridinium-4-yl radical cation

Aromatic radicals form in a variety of reacting gas-phase systems, where their molecular weight growth reactions with unsaturated hydrocarbons are of considerable importance. We have investigated the ion-molecule reaction of the aromatic distonic N-methyl-pyridinium-4-yl (NMP) radical cation with 2-butyne (CH3C CCH3) using ion trap mass spectrometry. Comparison is made to high-level ab initio energy surfaces for the reaction of NMP and for the neutral phenyl radical system. The NMP radical cation reacts rapidly with 2-butyne at ambient temperature, due to the apparent absence of any barrier. The activated vinyl radical adduct predominantly dissociates via loss of a H atom, with lesser amounts of CH3 loss. High-resolution Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry allows us to identify small quantities of the collisionally deactivated reaction adduct. Statistical reaction rate theory calculations (master equation/RRKM theory) on the NMP + 2-butyne system support our experimental findings, and indicate a mechanism that predominantly involves an allylic resonance-stabilized radical formed via H atom shuttling between the aromatic ring and the C-4 side-chain, followed by cyclization and/or low-energy H atom beta-scission reactions. A similar mechanism is demonstrated for the neutral phenyl radical (Ph center dot)+2-butyne reaction, forming products that include 3-methylindene. The collisionally deactivated reaction adduct is predicted to be quenched in the form of a resonance-stabilized methylphenylallyl radical. Experiments using a 2,5-dichloro substituted methyl-pyridiniumyl radical cation revealed that in this case CH3 loss from the 2-butyne adduct is favoured over H atom loss, verifying the key role of ortho H atoms, and the shuttling mechanism, in the reactions of aromatic radicals with alkynes. As well as being useful phenyl radical analogues, pyridiniumyl radical cations may form in the ionosphere of Titan, where they could undergo rapid molecular weight growth reactions to yield polycyclic aromatic nitrogen hydrocarbons (PANHs).

Ion-molecule reactions reveal facile radical migration in peptides

Ion-molecule reactions between molecular oxygen and peptide radicals in the gas phase demonstrate that radical migration occurs easily within large biomolecules without addition of collisional activation energy.

Base-induced decomposition of alkyl hydroperoxides in the gas phase. Part 3. Kinetics and dynamics in HO-+CH3OOH, C2H5OOH, and tert-C4H9OOH reactions

The E-CO(2) elimination reactions of alkyl hydroperoxides proceed via abstraction of an (x-hydrogen by a base: X- + (RRHCOOH)-R-1-H-2 -> HX + (RRC)-R-1-C-2=O + HO-. Efficiencies and product distributions for the reactions of the hydroxide anion with methyl, ethyl, and tert-butyl hydroperoxides are studied in the gas phase. On the basis of experiments using three isotopic analogues, HO- + CH3OOH, HO- + CD3OOH, and H18O- + CH3OOH. the overall intrinsic reaction efficiency is determined to be 80% or greater. The E(CO)2 decomposition is facile for these methylperoxide reactions, and predominates over competing proton transfer at the hydroperoxide moiety. The CH3CH2OOH reaction displays a similar E(CO)2 reactivity, whereas proton transfer and the formation of HOO- are the exclusive pathways observed for (CH3)(3)COOH, which has no (x-hydrogen. All results are consistent with the E-CO(2) mechanism, transition state structure, and reaction energy diagrams calculated using the hybrid density functional B3LYP approach. Isotope labeling for HO- + CH3OOH also reveals some interaction between H2O and HO- within the E(CO)2 product complex [H2O center dot center dot center dot CH2=O center dot center dot center dot HO-]. There is little evidence, however. for the formation of the most exothermic products H2O + CH2(OH)O-, which would arise from nuclephilic condensation of CH2=O and HO-. The results suggest that the product dynamics are not totally statistical but are rather direct after the E-CO(2) transition state. The larger HO- + CH3CH2OOH system displays more statistical behavior during complex dissociation.

Reactions of the hydroperoxide anion with dimethyl methylphosphonate in an ion trap mass spectrometer : evidence for a gas phase alpha-effect

The gas phase degradation reactions of the chemical warfare agent (CWA) simulant, dimethyl methylphosphonate (DMMP), with the hydroperoxide anion (HOO(-)) were investigated using a modified quadrupole ion trap mass spectrometer. The HOO(-) anion reacts readily with neutral DMMP forming two significant product ions at m/z 109 and m/z 123. The major reaction pathways correspond to (i) the nucleophilic substitution at carbon to form \[CH(3)P(O)(OCH(3))O](-) (m/z 109) in a highly exothermic process and (ii) exothermic proton transfer. The branching ratios of the two reaction pathways, 89% and 11% respectively, indicate that the former reaction is significantly faster than the latter. This is in contrast to the trend for the methoxide anion with DMMP, where proton transfer dominates. The difference in the observed reactivities of the HOO(-) and CH(3)O(-) anions can be considered as evidence for an a-effect in the gas phase and is supported by electronic structure calculations at the B3LYP/aug-cc-pVTZ//B3LYP/6-31+G(d) level of theory that indicate the S(N)2(carbon) process has an activation energy 7.8 kJ mol(-1) lower for HOO(-) as compared to CH(3)O(-). A similar alpha-effect was calculated for nucleophilic addition-elimination at phosphorus, but this process an important step in the perhydrolysis degradation of CWAs in solution - was not observed to occur with DMMP in the gas phase. A theoretical investigation revealed that all processes are energetically accessible with negative activation energies. However, comparison of the relative Arrhenius pre-exponential factors indicate that substitution at phosphorus is not kinetically competitive with respect to the S(N)2(carbon) and deprotonation processes.

Ion-molecule reactions of O,S-Dimethyl methylphosphonothioate : Evidence for intramolecular sulfur oxidation during VX perhydrolysis

The alkaline perhydrolysis of the nerve agent O-ethyl S-[2-(diisopropylamino)ethyl] methylphosphonothioate (VX) was investigated by studying the ion-molecule reactions of HOO(-) with O,S-dimethyl methylphosphonothioate in a modified linear ion-trap mass spectrometer. In addition to simple proton transfer, two other abundant product ions are observed at m/z 125 and 109 corresponding to the S-methyl methylphosphonothioate and methyl methylphosphonate anions, respectively. The structure of these product ions is demonstrated by a combination of collision-induced dissociation and isotope-labeling experiments that also provide evidence for their formation by nucleophilic reaction pathways, namely, (i) S(N)2 at carbon to yield the S-methyl methylphosphonothioate anion and (ii) nucleophilic addition at phosphorus affording a reactive pentavalent intermediate that readily undergoes internal sulfur oxidation and concomitant elimination of CH(3)SOH to yield the methyl methylphosphonate anion. Consistent with previous Solution phase observations of VX perhydrolysis, the toxic P-O cleavage product is not observed in this VX model system and theoretical calculations identify P-O cleavage to be energetically uncompetitive. Conversely, intramolecular sulfur oxidation is calculated to be extremely exothermic and kinetically accessible explaining its competitiveness with the facile gas phase proton transfer process. Elimination of a sulfur moiety deactivates the nerve agent VX and thus the intramolecular sulfur oxidation process reported here is also able to explain the selective perhydrolysis of the nerve agent to relatively nontoxic products.

Control and diagnostics of inductively coupled plasmas for chemical vapour deposition of nanocomposite carbon nitride-based films

Control and diagnostics of low-frequency (∼ 500 kHz) inductively coupled plasmas for chemical vapor deposition (CVD) of nano-composite carbon nitride-based films is reported. Relation between the discharge control parameters, plasma electron energy distribution/probability functions (EEDF/EEPF), and elemental composition in the deposited C-N based thin films is investigated. Langmuir probe technique is employed to monitor the plasma density and potential, effective electron temperature, and EEDFs/EEPFs in Ar + N2 + CH4 discharges. It is revealed that varying RF power and gas composition/pressure one can engineer the EEDFs/EEPFs to enhance the desired plasma-chemical gas-phase reactions thus controlling the film chemical structure. Auxiliary diagnostic tools for study of the RF power deposition, plasma composition, stability, and optical emission are discussed as well.

Visible light-driven cross-coupling reactions at lower temperatures using a photocatalyst of palladium and gold alloy nanoparticles

Palladium (Pd)-catalyzed cross-coupling reactions are among the most important methods in organic synthesis. We report the discovery of highly efficient and green photocatalytic processes by which cross-coupling reactions, including Sonogashira, Stille, Hiyama, Ullmann, and Buchwald–Hartwig reactions, can be driven with visible light at temperatures slightly above room temperature using alloy nanoparticles of gold and Pd on zirconium oxide, thus achieving high yields. The alloy nanoparticles absorb visible light, and their conduction electrons gain energy, which is available at the surface Pd sites. Results of the density functional theory calculations indicate that transfer of the light excited electrons from the nanoparticle surface to the reactant molecules adsorbed on the nanoparticle surface activates the reactants. When the light intensity was increased, a higher reaction rate was observed, because of the increased population of photoexcited electrons. The irradiation wavelength also has an important impact on the reaction rates. Ultraviolet irradiation can drive some reactions with the chlorobenzene substrate, while visible light irradiation failed to, and substantially improve the yields of the reactions with the bromobenzene substrate. The discovery reveals the possibility of using low-energy and -density sources such as sunlight to drive chemical transformations.

Rapid, selective, and reversible Nitroxide Radical Coupling (NRC) reactions at amibient temperature

High activation of polystyrene with bromine end groups (PSTY-Br) to their incipient radicals occurred in the presence of Cu(I)Br, Me6TREN, and DMSO solvent. These radicals were then trapped by nitroxide species leading to coupling reactions between PSTY-Br and nitroxides that were ultrafast and selective in the presence of a diverse range of functional groups. The nitroxide radical coupling (NRC) reactions have the attributes of a “click” reaction with near quantitative yields of product formed, but through the reversibility of this reaction, it has the added advantage of permitting the exchange of chemical functionality on macromolecules. Conditions were chosen to facilitate the disproportionation of Cu(I)Br to the highly activating nascent Cu(0) and deactivating Cu(II)Br2 in the presence of DMSO solvent and Me6TREN ligand. NRC at room temperature gave near quantitative yields of macromolecular coupling of low molecular weight polystyrene with bromine chain-ends (PSTY-Br) and nitroxides in under 7 min even in the presence of functional groups (e.g., −≡, −OH, −COOH, −NH2, =O). Utilization of the reversibility of the NRC reaction at elevated temperatures allowed the exchange of chain-end groups with a variety of functional nitroxide derivatives. The robustness and orthogonality of this NRC reaction were further demonstrated using the Cu-catalyzed azide/alkyne “click” (CuAAC) reactions, in which yields greater than 95% were observed for coupling between PSTY-N3 and a PSTY chain first trapped with an alkyne functional TEMPO (PSTY-TEMPO-≡).