Stochastic simulation of chemical reactions in spatially complex media

Discrete stochastic simulations, via techniques such as the Stochastic Simulation Algorithm (SSA) are a powerful tool for understanding the dynamics of chemical kinetics when there are low numbers of certain molecular species. However, an important constraint is the assumption of well-mixedness and homogeneity. In this paper, we show how to use Monte Carlo simulations to estimate an anomalous diffusion parameter that encapsulates the crowdedness of the spatial environment. We then use this parameter to replace the rate constants of bimolecular reactions by a time-dependent power law to produce an SSA valid in cases where anomalous diffusion occurs or the system is not well-mixed (ASSA). Simulations then show that ASSA can successfully predict the temporal dynamics of chemical kinetics in a spatially constrained environment.

Stochastic chemical kinetics and the quasi-steady-state assumption : application to the stochastic simulation algorithm and chemical master equation

Recently the application of the quasi-steady-state approximation (QSSA) to the stochastic simulation algorithm (SSA) was suggested for the purpose of speeding up stochastic simulations of chemical systems that involve both relatively fast and slow chemical reactions [Rao and Arkin, J. Chem. Phys. 118, 4999 (2003)] and further work has led to the nested and slow-scale SSA. Improved numerical efficiency is obtained by respecting the vastly different time scales characterizing the system and then by advancing only the slow reactions exactly, based on a suitable approximation to the fast reactions. We considerably extend these works by applying the QSSA to numerical methods for the direct solution of the chemical master equation (CME) and, in particular, to the finite state projection algorithm [Munsky and Khammash, J. Chem. Phys. 124, 044104 (2006)], in conjunction with Krylov methods. In addition, we point out some important connections to the literature on the (deterministic) total QSSA (tQSSA) and place the stochastic analogue of the QSSA within the more general framework of aggregation of Markov processes. We demonstrate the new methods on four examples: Michaelis–Menten enzyme kinetics, double phosphorylation, the Goldbeter–Koshland switch, and the mitogen activated protein kinase cascade. Overall, we report dramatic improvements by applying the tQSSA to the CME solver.

A multi-scaled approach for simulating chemical reaction systems

In this paper we give an overview of some very recent work, as well as presenting a new approach, on the stochastic simulation of multi-scaled systems involving chemical reactions. In many biological systems (such as genetic regulation and cellular dynamics) there is a mix between small numbers of key regulatory proteins, and medium and large numbers of molecules. In addition, it is important to be able to follow the trajectories of individual molecules by taking proper account of the randomness inherent in such a system. We describe different types of simulation techniques (including the stochastic simulation algorithm, Poisson Runge-Kutta methods and the balanced Euler method) for treating simulations in the three different reaction regimes: slow, medium and fast. We then review some recent techniques on the treatment of coupled slow and fast reactions for stochastic chemical kinetics and present a new approach which couples the three regimes mentioned above. We then apply this approach to a biologically inspired problem involving the expression and activity of LacZ and LacY proteins in E coli, and conclude with a discussion on the significance of this work. (C) 2004 Elsevier Ltd. All rights reserved.

Inducing nonlocal reactions with a local probe

The scanning tunneling microscope (STM) has evolved continually since its invention, as scientists have expanded its use to encompass atomic-scale manipulation, momentum-resolved electronic characterization, localized chemical reactions (bond breaking and bond making) in adsorbed molecules, and even chain reactions at surfaces. This burgeoning field has recently expanded to include the use of the STM to inject hot electrons into substrate surface states; the injected electrons can travel laterally and induce changes in chemical structure in molecules located up to 100 nm from the STM tip. We describe several key demonstrations of this phenomenon, including one appearing in this issue of ACS Nano by Chen et al. Possible applications for this technique are also discussed, including characterizing the dispersion of molecule−substrate interface states and the controlled patterning of molecular overlayers.

Experimental warming shows that decomposition temperature sensitivity increases with soil organic matter recalcitrance

Soil C decomposition is sensitive to changes in temperature, and even small increases in temperature may prompt large releases of C from soils. But much of what we know about soil C responses to global change is based on short-term incubation data and model output that implicitly assumes soil C pools are composed of organic matter fractions with uniform temperature sensitivities. In contrast, kinetic theory based on chemical reactions suggests that older, more-resistant C fractions may be more temperature sensitive. Recent research on the subject is inconclusive, indicating that the temperature sensitivity of labile soil organic matter (OM) decomposition could either be greater than, less than, or equivalent to that of resistant soil OM. We incubated soils at constant temperature to deplete them of labile soil OM and then successively assessed the CO2-C efflux in response to warming. We found that the decomposition response to experimental warming early during soil incubation (when more labile C remained) was less than that later when labile C was depleted. These results suggest that the temperature sensitivity of resistant soil OM pools is greater than that for labile soil OM and that global change-driven soil C losses may be greater than previously estimated.

Facile synthesis of Lithium Niobate squares by a combustion route

A combustion synthesis of lithium niobate (LN) squares from activated niobium oxide (Nb2 O5.nH2O) and Li2CO3 was studied to understand all the chemical reactions involved, and the nucleation and square-growth mechanisms. It was found that first the lithium ions react with the fuel (urea), then niobium ions of Nb2 O5.nH2O begin a continuous reaction with the fuel to form metal-organic complexes. LN nuclei are formed by the solid-state reaction of Li- and Nb-organic complexes at 430 degrees celcius. Lithium niobate squares are obtained in the crystallization stasge at 700 degrees celcius, which go on the grow into larger squares at 850 degrees celcius because of the agglomeration effect.

Synthesis and characterisation of metal oxyhydroxide and oxide nanomaterials

In this work, a range of nanomaterials have been synthesised based on metal oxyhydroxides MO(OH), where M=Al, Co, Cr, etc. Through a self-assembly hydrothermal route, metal oxyhydroxide nanomaterials with various morphologies were successfully synthesised: one dimensional boehmite (AlO(OH)) nanofibres, zero dimensional indium hydroxide (In(OH)3) nanocubes and chromium oxyhydroxide (CrO(OH)) nanoparticles, as well as two dimensional cobalt hydroxide and oxyhydroxide (Co(OH)2 & CoO(OH)) nanodiscs. In order to control the synthetic nanomaterial morphology and growth, several factors were investigated including cation concentration, temperature, hydrothermal treatment time, and pH. Metal ion doping is a promising technique to modify and control the properties of materials by intentionally introducing impurities or defects into the material. Chromium was successfully applied as a dopant for fabricating doped boehmite nanofibres. The thermal stability of the boehmite nanofibres was enhanced by chromium doping, and the photoluminescence property was introduced to the chromium doped alumina nanofibres. Doping proved to be an efficient method to modify and functionalize nanomaterials. The synthesised nanomaterials were fully characterised by X-ray diffraction (XRD), transmission electron microscopy (TEM) combined with selected area electron diffraction (SAED), scanning electron microscopy (SEM), BET specific surface area analysis, X-ray photoelectron spectroscopy (XPS) and thermo gravimetric analysis (TGA). Hot-stage Raman and infrared emission spectroscopy were applied to study the chemical reactions during dehydration and dehydroxylation. The advantage of these techniques is that the changes in molecular structure can be followed in situ and at the elevated temperatures.

Look before you leap : a confidence-based method for selecting species criticality while avoiding negative populations in τ-leaping

The stochastic simulation algorithm was introduced by Gillespie and in a different form by Kurtz. There have been many attempts at accelerating the algorithm without deviating from the behavior of the simulated system. The crux of the explicit τ-leaping procedure is the use of Poisson random variables to approximate the number of occurrences of each type of reaction event during a carefully selected time period, τ. This method is acceptable providing the leap condition, that no propensity function changes “significantly” during any time-step, is met. Using this method there is a possibility that species numbers can, artificially, become negative. Several recent papers have demonstrated methods that avoid this situation. One such method classifies, as critical, those reactions in danger of sending species populations negative. At most, one of these critical reactions is allowed to occur in the next time-step. We argue that the criticality of a reactant species and its dependent reaction channels should be related to the probability of the species number becoming negative. This way only reactions that, if fired, produce a high probability of driving a reactant population negative are labeled critical. The number of firings of more reaction channels can be approximated using Poisson random variables thus speeding up the simulation while maintaining the accuracy. In implementing this revised method of criticality selection we make use of the probability distribution from which the random variable describing the change in species number is drawn. We give several numerical examples to demonstrate the effectiveness of our new method.

Vibrational spectroscopy of synthetic stercorite H(NH4)Na(PO4)•4H2O – a comparison with the natural cave mineral

In order to mimic the chemical reactions in cave systems, the analogue of the mineral stercorite H(NH4)Na(PO4)•4H2O has been synthesised. X-ray diffraction of the stercorite analogue matches the stercorite reference pattern. A comparison is made with the vibrational spectra of synthetic stercorite analogue and the natural Cave mineral. The mineral in nature is formed by the reaction of bat guano chemicals on calcite substrates. A single Raman band at 920 cm-1 (Cave) and 922 cm-1 (synthesised) defines the presence of hydrogen phosphate in the mineral. In the synthetic stercorite analogue, additional bands are observed and are attributed to the dihydrogen and phosphate anions. The vibrational spectra of synthetic stercorite only partly match that of the natural stercorite. It is suggested that natural stercorite is more pure than that of synthesised stercorite. Antisymmetric stretching bands are observed in the infrared spectrum at 1052, 1097, 1135 and 1173 cm-1. Raman spectroscopy shows the stercorite mineral is based upon the hydrogen phosphate anion and not the phosphate anion. Raman and infrared bands are found and assigned to PO43-, H2O, OH and NH stretching vibrations. Raman spectroscopy shows the synthetic analogue is similar to the natural mineral. A mechanism for the formation of stercorite is provided.

Noble metal photocatalysts under visible light and UV light irradiation

One of the greatest challenges for the study of photocatalysts is to devise new catalysts that possess high activity under visible light illumination. This would allow the use of an abundant and green energy source, sunlight, to drive chemical reactions. Gold nanoparticles strongly absorb both visible light and UV light. It is therefore possible to drive chemical reactions utilising a significant fraction of full sunlight spectrum. Here we prepared gold nanoparticles supported on various oxide powders, and reported a new finding that gold nanoparticles on oxide supports exhibit significant activity for the oxidation of formaldehyde and methanol in the air at ambient temperature, when illuminated with visible light. We suggested that visible light can greatly enhance local electromagnetic fields and heat gold nanoparticles due to surface plasmon resonance effect which provides activation energy for the oxidation of organic molecules. Moreover, the nature of the oxide support has an important influence on the activity of the gold nanoparticles. The finding reveals the possibility to drive chemical reactions with sunlight on gold nanoparticles at ambient temperature, highlighting a new direction for research on visible light photocatalysts. Gold nanoparticles supported on oxides also exhibit significant dye oxidation activity under visible light irradiation in aqueous solution at ambient temperature. Turnover frequencies of the supported gold nanoparticles for the dye degradation are much higher than titania based photocatalysts under both visible and UV light. These gold photocatalysts can also catalyse phenol degradation as well as selective oxidation of benzyl alcohol under UV light. The reaction mechanism for these photocatalytic oxidations was studied. Gold nanoparticles exhibit photocatalytic activity due to visible light heating gold electrons in 6sp band, while the UV absorption results in electron holes in gold 5d band to oxidise organic molecules. Silver nanoparticles also exhibit considerable visible light and UV light absorption due to surface plasmon resonance effect and the interband transition of 4d electrons to the 5sp band, respectively. Therefore, silver nanoparticles are potentially photocatalysts that utilise the solar spectrum effectively. Here we reported that silver nanoparticles at room temperature can be used to drive chemical reactions when illuminated with light throughout the solar spectrum. The significant activities for dye degradation by silver nanoparticles on oxide supports are even better than those by semiconductor photocatalysts. Moreover, silver photocatalysts also can degrade phenol and drive the oxidation of benzyl alcohol to benzaldehyde under UV light. We suggested that surface plasmon resonance effect and interband transition of silver nanoparticles can activate organic molecule oxidations under light illumination.

Reaction products between Bi-Sr-Ca-Cu-oxide thick films and alumina substrates

The structure and composition of reaction products between Bi-Sr-Ca-Cu-oxide (BSCCO) thick films and alumina substrates have been characterized using a combination of electron diffraction, scanning electron microscopy and energy dispersive X-ray spectrometry (EDX). Sr and Ca are found to be the most reactive cations with alumina. Sr4Al6O12SO4 is formed between the alumina substrates and BSCCO thick films prepared from paste with composition close to Bi-2212 (and Bi-2212 + 10 wt.% Ag). For paste with composition close to Bi(Pb)-2223 + 20 wt.% Ag, a new phase with f.c.c. structure, lattice parameter about a = 24.5 A and approximate composition Al3Sr2CaBi2CuOx has been identified in the interface region. Understanding and control of these reactions is essential for growth of high quality BSCCO thick films on alumina. (C) 1997 Elsevier Science S.A.

Titanium-oxide magneli phases in four chondritic porous interplanetary dust particles

Detailed analytical electron microscope analyses of four fine-grained chondritic porous interplanetary dust particles (IDPs)reveal the presence of titanium oxide Magneli phases, TinO2n-1 (n=4,5,6), and rare Ti-metal. The titanium minerals are indigenous to these chondritic IDPs. The association of Magneli phases, Ti-metal, and carbonaceous material in chondritic IDPs, along with the grain size distributions support in situ solid carbon gasification in these extraterrestrial particles. The active catalyst in this process is titanium metal that we infer may be of interstellar origin. This favorable catalysis uniquely leads to the formation of Magneli phases. As chondritic IDPs may be solid debris of short-period comets, our data indicate that nuclei of short-period comets may show distinctive chemical reactions that lead to Ti-mineral assemblages that typically include Magneli phases. The proposed model provides a plausible mechnism to explain the higher solid carbon content of chondritic IDPs relative to bulk carbon abundances typical for carbonaceous chondrite matrices that represent another type of more evolved, that is, metamorphosed, undifferentiated solar system bodies.

Recent advances in catalytic/biocatalytic conversion of greenhouse methane and carbon dioxide to methanol and other oxygenates

Methane gas has been identified as the most destructive greenhouse gas (Liu et al., 2004). It was reported that the global warming potential of methane per molecule relative to CO2 is approximately 23 on a 100-year timescale or 62 over a 20-year period (IPCC, 2001). Methane has high C-H bond energy of about 439 kJ/mol and other higher alkanes (or saturated hydrocarbons) also have a very strong C-C and C-H bonds, thus making their molecules to have no empty orbitals of low energy or filled orbitals of high energy that could readily participate in chemical reactions as is the case with unsaturated hydrocarbons such as olefins and alkynes (Crabtree, 1994; Labinger & Bercaw, 2002)...

Development of modified inorganic adsorbents for radioactive iodine removal and biomolecule adsorption

Development and application of inorganic adsorbent materials have been continuously investigated due to their variability and versatility. This Master thesis has expanded the knowledge in the field of adsorption targeting radioactive iodine waste and proteins using modified inorganic materials. Industrial treatment of radioactive waste and safety disposal of nuclear waste is a constant concern around the world with the development of radioactive materials applications. To address the current problems, laminar titanate with large surface area (143 m2 g−1) was synthesized from inorganic titanium compounds by hydrothermal reactions at 433 K. Ag2O nanocrystals of particle size ranging from 5–30 nm were anchored on the titanate lamina surface which has crystallographic similarity to that of Ag2O nanocrystals. Therefore, the deposited Ag2O nanocrystals and titanate substrate could join together at these surfaces between which there forms a coherent interface. Such coherence between the two phases reduces the overall energy by minimizing surface energy and maintains the Ag2O nanocrystals firmly on the outer surface of the titanate structure. The combined adsorbent was then applied as efficient adsorbent to remove radioactive iodine from water (one gram adsorbent can capture up to 3.4 mmol of I- anions) and the composite adsorbent can be recovered easily for safe disposal. The structure changes of the titanate lamina and the composite adsorbent were characterized via various techniques. The isotherm and kinetics of iodine adsorption, competitive adsorption and column adsorption using the adsorbent were studied to determine the iodine removal abilities of the adsorbent. It is shown that the adsorbent exhibited excellent trapping ability towards iodine in the fix-bed column despite the presence of competitive ions. Hence, Ag2O deposited titanate lamina could serve as an effective adsorbent for removing iodine from radioactive waste. Surface hydroxyl group of the inorganic materials is widely applied for modification purposes and modification of inorganic materials for biomolecule adsorption can also be achieved. Specifically, γ-Al2O3 nanofibre material is converted via calcinations from boehmite precursor which is synthesised by hydrothermal chemical reactions under directing of surfactant. These γ-Al2O3 nanofibres possess large surface area (243 m2 g-1), good stability under extreme chemical conditions, good mechanical strength and rich surface hydroxyl groups making it an ideal candidate in industrialized separation column. The fibrous morphology of the adsorbent also guarantees facile recovery from aqueous solution under both centrifuge and sedimentation approaches. By chemically bonding the dyes molecules, the charge property of γ-Al2O3 is changed in the aim of selectively capturing of lysozyme from chicken egg white solution. The highest Lysozyme adsorption amount was obtained at around 600 mg/g and its proportion is elevated from around 5% to 69% in chicken egg white solution. It was found from the adsorption test under different solution pH that electrostatic force played the key role in the good selectivity and high adsorption rate of surface modified γ-Al2O3 nanofibre adsorbents. Overall, surface modified fibrous γ-Al2O3 could be applied potentially as an efficient adsorbent for capturing of various biomolecules.