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.

Vertical growth and resonator properties of hexagonally shaped zinc oxide nanonails

We synthesized vertically aligned nail-shaped ZnO nanocrystal arrays on silicon substrates via a combination of a carbothermal reduction method and textured ZnO seeding layers that were precoated on silicon substrates by thermally decomposing zinc acetate, and studied their optical properties using cathodoluminescence (CL) and photoluminescence techniques. The ZnO nanonails show a sharp band-gap edge UV emission and a defect-related broad green emission. Monochromatic CL images of an individual ZnO nanonail show variations in spatial distributions of respective CL bands that had different origins. We attribute the spatial variation of CL images to an uneven distribution of luminescent defects and/or a structure-related light out-coupling from hexagonal ZnO nanostructures. The most distinct CL feature from the hexagonal head of an individual ZnO nanonail was the occurrence of a series of distinct resonant peaks within the visible wavelength range. It appeared that the head of a nanonail played the role of a hexagonal cavity so that polarizationdependent whispering gallery modes were stimulated by electron beam excitation.

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.

A software tool for simulating spatial separation and kinetics in biological membranes

Self-segregation and compartimentalisation are observed experimentally to occur spontaneously on live membranes as well as reconstructed model membranes. It is believed that many of these processes are caused or supported by anomalous diffusive behaviours of biomolecules on membranes due to the complex and heterogeneous nature of these environments. These phenomena are on the one hand of great interest in biology, since they may be an important way for biological systems to selectively localize receptors, regulate signaling or modulate kinetics; and on the other, they provide an inspiration for engineering designs that mimick natural systems. We present an interactive software package we are developing for the purpose of simulating such processes numerically using a fundamental Monte Carlo approach. This program includes the ability to simulate kinetics and mass transport in the presence of either mobile or immobile obstacles and other relevant structures such as liquid-ordered lipid microdomains. We also present preliminary simulation results regarding the selective spatial localization and chemical kinetics modulating power of immobile obstacles on the membrane, obtained using the program.

Kinetics of electron recombination of dye-senitized solar cells based on TiO2 nanorod arrays sensitized with different dyes

The performance and electron recombination kinetics of dye-sensitized solar cells based on TiO2 films consisting of one-dimensional nanorod arrays (NR-DSSCs) which are sensitized with dye N719, C218 and D205 respectively have been studied. It has been found that the best efficiency is obtained with the dye C218 based NR-DSSCs, benefiting from a 40% higher short-circuit photocurrent density. However, the open circuit photovoltage of the N719 based cell is 40 mV higher than that of the organic dye C218 and D205 based devices. Investigation of the electron recombination kinetics of the NR-DSSCs has revealed that the effective electron lifetime, τn, of the N719 based NR-DSSC is the lowest whereas the τn of the C218 based NR-DSSC is the highest among the three dyes. The higher Voc with the N719 based NR-DSSC is originated from the more negative energy level of the conduction band of the TiO2 film. In addition, in comparison to the DSSCs with conventional nanocrystalline particles based TiO2 films, the NR-DSSCs have shown over two orders of magnitude higher τn when employing N719 as the sensitizer. Nevertheless, the τn of the DSSCs with the C218 based nanorod arrays is only ten-fold higher than the that of the nanoparticles based devices. The remarkable characteristic of the dye C218 in suppressing the electron recombination of DSSCs is discussed.

Investigation of phenol degradation : True reaction kinetics on fixed film titanium dioxide photocatalyst

This study of photocatalytic oxidation of phenol over titanium dioxide films presents a method for the evaluation of true reaction kinetics. A flat plate reactor was designed for the specific purpose of investigating the influence of various reaction parameters, specifically photocatalytic film thickness, solution flow rate (1–8 l min−1), phenol concentration (20, 40 and 80 ppm), and irradiation intensity (70.6, 57.9, 37.1and 20.4 W m−2), in order to further understand their impact on the reaction kinetics. Special attention was given to the mass transfer phenomena and the influence of film thickness. The kinetics of phenol degradation were investigated with different irradiation levels and initial pollutant concentration. Photocatalytic degradation experiments were performed to evaluate the influence of mass transfer on the reaction and, in addition, the benzoic acid method was applied for the evaluation of mass transfer coefficient. For this study the reactor was modelled as a batch-recycle reactor. A system of equations that accounts for irradiation, mass transfer and reaction rate was developed to describe the photocatalytic process, to fit the experimental data and to obtain kinetic parameters. The rate of phenol photocatalytic oxidation was described by a Langmuir–Hinshelwood type law that included competitive adsorption and degradation of phenol and its by-products. The by-products were modelled through their additive effect on the solution total organic carbon.

A fibrocontractive mechanochemical model of dermal wound closure incorporating realistic growth factor kinetics

Fibroblasts and their activated phenotype, myofibroblasts, are the primary cell types involved in the contraction associated with dermal wound healing. Recent experimental evidence indicates that the transformation from fibroblasts to myofibroblasts involves two distinct processes: the cells are stimulated to change phenotype by the combined actions of transforming growth factor β (TGFβ) and mechanical tension. This observation indicates a need for a detailed exploration of the effect of the strong interactions between the mechanical changes and growth factors in dermal wound healing. We review the experimental findings in detail and develop a model of dermal wound healing that incorporates these phenomena. Our model includes the interactions between TGFβ and collagenase, providing a more biologically realistic form for the growth factor kinetics than those included in previous mechanochemical descriptions. A comparison is made between the model predictions and experimental data on human dermal wound healing and all the essential features are well matched.

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.

Zinc aluminium layered double hydroxides for the removal of iodine and iodide from aqueous solutions

129I is a radioactive isotope of iodine that is readily absorbed by the body. In this paper we investigated the potential of a 3:1 Zn/Al layered double hydroxide (LDH) as a sorbent for the removal of iodine and iodide from water. Synthetic Zn6Al2(OH)16(CO3)∙4H2O was prepared by the co-precipitation before thermal activation. The LDH was treated with solutions containing iodide and iodine. It was found that iodine could be more easily removed from solution than iodide. Powder X-ray diffraction revealed the destruction of the LDH structure during thermal activation and the successful reformation of a similar LDH material after treatment with the iodide or iodine solution. Thermal decomposition of all samples studied by thermogravimetry appeared to be similar. A new decomposition mechanism similar to one previously described in the literature was proposed for the Zn/Al LDH. The total mass loss of samples treated with iodide and iodine was significantly lower than that of the original LDH indicating that iodine species may form non-removable anions when intercalated into the LDH structure. Evolved gas mass spectrometry failed to detect any iodine species lost as gases during the decomposition of iodide treated LDH however, small quantities of iodine species were observed during decomposition of samples treated with iodine solution.

Gold coating of silica and zinc oxide nanoparticles by the surface reduction of gold(I) chloride

The possibility of a surface inner sphere electron transfer mechanism leading to the coating of gold via the surface reduction of gold(I) chloride on metal and semi-metal oxide nanoparticles was investigated. Silica and zinc oxide nanoparticles are known to have very different surface chemistry, potentially leading to a new class of gold coated nanoparticles. Monodisperse silica nanoparticles were synthesised by the well known Stöber protocol in conjunction with sonication. The nanoparticle size was regulated solely by varying the amount of ammonia solution added. The presence of surface hydroxyl groups was investigated by liquid proton NMR. The resultant nanoparticle size was directly measured by the use of TEM. The synthesised silica nanoparticles were dispersed in acetonitrile (MeCN) and added to a bis acetonitrile gold(I) co-ordination complex [Au(MeCN)2]+ in MeCN. The silica hydroxyl groups were deprotonated in the presence of MeCN generating a formal negative charge on the siloxy groups. This allowed the [Au(MeCN)2]+ complex to undergo ligand exchange with the silica nanoparticles, which formed a surface co-ordination complex with reduction to gold(0), that proceeded by a surface inner sphere electron transfer mechanism. The residual [Au(MeCN)2]+ complex was allowed to react with water, disproportionating into gold(0) and gold(III) respectively, with gold(0) being added to the reduced gold already bound on the silica surface. The so-formed metallic gold seed surface was found to be suitable for the conventional reduction of gold(III) to gold(0) by ascorbic acid. This process generated a thin and uniform gold coating on the silica nanoparticles. This process was modified to include uniformly gold coated composite zinc oxide nanoparticles (Au@ZnO NPs) using surface co-ordination chemistry. AuCl dissolved in acetonitrile (MeCN) supplied chloride ions which were adsorbed onto ZnO NPs. The co-ordinated gold(I) was reduced on the ZnO surface to gold(0) by the inner sphere electron transfer mechanism. Addition of water disproportionated the remaining gold(I) to gold(0) and gold(III). Gold(0) bonded to gold(0) on the NP surface with gold(III) was reduced to gold(0) by ascorbic acid (ASC), which completed the gold coating process. This gold coating process of Au@ZnO NPs was modified to incorporate iodide instead of chloride. ZnO NPs were synthesised by the use of sodium oxide, zinc iodide and potassium iodide in refluxing basic ethanol with iodide controlling the presence of chemisorbed oxygen. These ZnO NPs were treated by the addition of gold(I) chloride dissolved in acetonitrile leaving chloride anions co-ordinated on the ZnO NP surface. This allowed acetonitrile ligands in the added [Au(MeCN)2]+ complex to surface exchange with adsorbed chloride from the dissolved AuCl on the ZnO NP surface. Gold(I) was then reduced by the surface inner sphere electron transfer mechanism. The presence of the reduced gold on the ZnO NPs allowed adsorption of iodide to generate a uniform deposition of gold onto the ZnO NP surface without the use of additional reducing agents or heat.

Intrinsic kinetics of titania photocatalysis : simplified models for their investigation

Even though titanium dioxide photocatalysis has been promoted as a leading green technology for water purification, many issues have hindered its application on a large commercial scale. For the materials scientist the main issues have centred the synthesis of more efficient materials and the investigation of degradation mechanisms; whereas for the engineers the main issues have been the development of appropriate models and the evaluation of intrinsic kinetics parameters that allow the scale up or re-design of efficient large-scale photocatalytic reactors. In order to obtain intrinsic kinetics parameters the reaction must be analysed and modelled considering the influence of the radiation field, pollutant concentrations and fluid dynamics. In this way, the obtained kinetic parameters are independent of the reactor size and configuration and can be subsequently used for scale-up purposes or for the development of entirely new reactor designs. This work investigates the intrinsic kinetics of phenol degradation over titania film due to the practicality of a fixed film configuration over a slurry. A flat plate reactor was designed in order to be able to control reaction parameters that include the UV irradiance, flow rates, pollutant concentration and temperature. Particular attention was paid to the investigation of the radiation field over the reactive surface and to the issue of mass transfer limited reactions. The ability of different emission models to describe the radiation field was investigated and compared to actinometric measurements. The RAD-LSI model was found to give the best predictions over the conditions tested. Mass transfer issues often limit fixed film reactors. The influence of this phenomenon was investigated with specifically planned sets of benzoic acid experiments and with the adoption of the stagnant film model. The phenol mass transfer coefficient in the system was calculated to be km,phenol=8.5815x10-7Re0.65(ms-1). The data obtained from a wide range of experimental conditions, together with an appropriate model of the system, has enabled determination of intrinsic kinetic parameters. The experiments were performed in four different irradiation levels (70.7, 57.9, 37.1 and 20.4 W m-2) and combined with three different initial phenol concentrations (20, 40 and 80 ppm) to give a wide range of final pollutant conversions (from 22% to 85%). The simple model adopted was able to fit the wide range of conditions with only four kinetic parameters; two reaction rate constants (one for phenol and one for the family of intermediates) and their corresponding adsorption constants. The intrinsic kinetic parameters values were defined as kph = 0.5226 mmol m-1 s-1 W-1, kI = 0.120 mmol m-1 s-1 W-1, Kph = 8.5 x 10-4 m3 mmol-1 and KI = 2.2 x 10-3 m3 mmol-1. The flat plate reactor allowed the investigation of the reaction under two different light configurations; liquid and substrate side illumination. The latter of particular interest for real world applications where light absorption due to turbidity and pollutants contained in the water stream to be treated could represent a significant issue. The two light configurations allowed the investigation of the effects of film thickness and the determination of the catalyst optimal thickness. The experimental investigation confirmed the predictions of a porous medium model developed to investigate the influence of diffusion, advection and photocatalytic phenomena inside the porous titania film, with the optimal thickness value individuated at 5 ìm. The model used the intrinsic kinetic parameters obtained from the flat plate reactor to predict the influence of thickness and transport phenomena on the final observed phenol conversion without using any correction factor; the excellent match between predictions and experimental results provided further proof of the quality of the parameters obtained with the proposed method.

Vibrational spectroscopic characterization of the phosphate mineral phosphophyllite – Zn2Fe(PO4)2·4H2O, from Hagendorf Süd, Germany and in comparison with other zinc phosphates

This research was undertaken on phosphophyllite sample from the Hagendorf Süd pegmatite, Bavaria, Germany. Chemical analysis was carried out by Scanning Electron Microscope in the EDS mode and indicates a zinc and iron phosphate with partial substitution of manganese, which partially replaced iron. The calculated chemical formula of the studied sample was determined to be: Zn2(Fe0.65, Mn0.35)P1.00(PO4)2- �4(H2O). The intense Raman peak at 995 cm�1 is assigned to the m1 PO3� 4 symmetric stretching mode and the two Raman bands at 1073 and 1135 cm�1 to the m3 PO3� 4 antisymmetric stretching modes. The m4 PO3� 4 bending modes are observed at 505, 571, 592 and 653 cm�1 and the m2 PO3� 4 bending mode at 415 cm�1. The sharp Raman band at 3567 cm�1 attributed to the stretching vibration of OH units brings into question the actual formula of phosphophyllite. Vibrational spectroscopy enables an assessment of the molecular structure of phosphophyllite to be assessed.

Aspects of the kinetics and solubility of silica and calcium oxalate composites in sugar solutions

Fouling of industrial surfaces by silica and calcium oxalate can be detrimental to a number of process streams. Solution chemistry plays a large roll in the rate and type of scale formed on industrial surfaces. This study is on the kinetics and thermodynamics of SiO2 and calcium oxalate composite formation in solutions containing Mg2+ ions, trans-aconitic acid and sucrose, to mimic factory sugar cane juices. The induction time (ti) of silicic acid polymerization is found to be dependent on the sucrose concentration and SiO2 supersaturation ratio (SS). Generalized kinetic and solubility models are developed for SiO2 and calcium oxalate in binary systems using response surface methodology. The role of sucrose, Mg, trans-aconitic acid, a mixture of Mg and trans-aconitic acid, SiO2 SS ratio and Ca in the formation of com- posites is explained using the solution properties of these species including their ability to form complexes.

Plasma zinc and immune markers in runners in response to a moderate increase in training volume

Changes in plasma zinc concentration and markers of immune function were examined in a group of 10 male runners (n = 10) following a moderate increase in training over four weeks. Seven sedentary males acted as controls. Fasting blood samples were taken at rest, before (T0) and after (T4) four weeks of increased (+ 16 %) training and after two weeks of reduced (-31 %) training (T6). Blood was analysed for plasma zinc concentration, differential leucocyte counts, lymphocyte subpopulations and lymphocyte proliferation using incorporation of 3H-thymidine. The runners increased their training volume by 16 % over the four weeks. When compared with the nonathletes, the runners had lower concentrations of plasma zinc (p = 0.012), CD3 + (p = 0.042) and CD19 + lymphocytes (p = 0.010) over the four weeks. Lymphocyte proliferation in response to Concanavalin A stimulation was greater in the runners (p = 0.0090). Plasma zinc concentration and immune markers remained constant during the study. Plasma zinc concentration correlated with total leucocyte counts in the athletes at T6 (r = -0.72, p < 0.05) and with Pokeweed mitogen stimulation in the nonathletes at T6 (r = -0.92, p < 0.05). Therefore, athletes are unlikely to benefit from zinc supplementation during periods of moderately increased training volume.