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.

A Two-Step Sol–Gel Method for Synthesis of Nanoporous TiO2

This article reports a study of the effects of synthesis parameters on the preparation and formation of mesoporous titania nanopowders by employing a two-step sol-gel method. These materials displayed crystalline domains characteristic of anatase. The first step of the process involved the hydrolysis of titanium isopropoxide in a basic aqueous solution mediated by neutral surfactant. The solid product obtained from step 1 was then treated in an acidified ethanol solution containing the same titanium precursor to thicken the pore walls. Low pH and higher loading of the Ti precursor in step 2 produced better mesoporosity and crystallinity of titanium dioxide polymorphs. The resultant powder exhibited a high surface area (73.8 m(2)/g) and large pore volume (0.17 cm(3)/g) with uniform mesopores. These materials are envisaged to be used as precursors for mesoporous titania films as a wide band gap semiconductor in dye-sensitized nanocrystalline TiO2 solar cells.

A quasi-chemical viscosity model for fully liquid slags in the Al2O3-CaO-'FeO'-SiO2 system

A structurally based viscosity model for fully liquid silicate slags has been proposed and applied to the Al2O3-CaO-'FeO'-SiO2 system at metallic iron saturation. The model links the slag viscosity to the internal structure of melts through the concentrations of various anion/cation structural units (SUs). The concentrations of structural units are equivalent to the second nearest neighbor bond concentrations calculated by the quasi-chemical thermodynamic model. This viscosity model describes experimental data over the entire temperature and composition range within the Al2O3-CaO-'FeO'-SiO2 system at metallic iron saturation and can be extended to other industrial slag systems.

Development of a quasi-chemical viscosity model for fully liquid slags in the Al2O3-CaO-'FeO'-MgO-SiO2 system. Part 1. Description of the model and its application to the MgO, MgO-SiO2, Al2O3-MgO and CaO-MgO sub-systems

A structurally-based quasi-chemical viscosity model for fully liquid slags in the Al2O3 CaO-'FeO'-MgO-SiO2 system has been developed. The model links the slag viscosities to the internal structures of the melts through the concentrations of various Si0.5O, Me2/nn+O and Me1/nn+Si0.25O viscous flow structural units. The concentrations of these structural units are derived from a quasi-chemical thermodynamic model of the system. The model described in this series of papers enables the viscosities of liquid slags to be predicted within experimental uncertainties over the whole range of temperatures and compositions in the Al2O3 CaOMgO-SiO2 system.

Development of a quasi-chemical viscosity model for fully liquid slags in the Al2O3-CaO-'FeO'-MgO-SiO2 system. Part 2. A review of the experimental data and the model predictions for the Al2O3-CaO-MgO, CaO-MgO-SiO2 and Al2O3-MgO-SiO2 systems

A structurally-based quasi-chemical viscosity model has been developed for the Al2O3 CaO-'FeO'-MgO-SiO2 system. The model links the slag viscosity to the internal structure of melts through the concentrations of various anion/cation Si0.5O, Me2/nn+O and Me1/nn+Si0.25O viscous flow structural units. The concentrations of structural units are derived from the quasi-chemical thermodynamic model. The focus of the work described in the present paper is the analysis of experimental data and the viscosity models for fully liquid slags in the Al2O3-CaO-MgO, Al2O3 MgO-SiO2 and CaO-MgO-SiO2 systems.

Development of a quasi-chemical viscosity model for fully liquid slags in the Al2O3-CaO-'FeO'-MgO-SiO2 system. Part 3. Summary of the model predictions for the Al2O3-CaO-MgO-SiO2 system and its sub-systems

A structurally-based quasi-chemical viscosity model for fully liquid slags in the Al2O3 CaO-'FeO'-MgOSiO2 system has been developed. The focus of the work described in the present paper is the analysis of the experimental data and viscosity models in the quaternary system Al2O3 CaO-MgO-SiO2 and its subsystems. A review of the experimental data, viscometry methods used and viscosity models available in the Al2O3 CaO-MgO-SiO2 and its sub-systems is reported. The quasi-chemical viscosity model is shown to provide good agreement between experimental data and predictions over the whole compositional range.

Leading coordinate analysis of reaction pathways in proton chain transfer: Application to a two-proton transfer model for the green fluorescent protein

The ‘leading coordinate’ approach to computing an approximate reaction pathway, with subsequent determination of the true minimum energy profile, is applied to a two-proton chain transfer model based on the chromophore and its surrounding moieties within the green fluorescent protein (GFP). Using an ab initio quantum chemical method, a number of different relaxed energy profiles are found for several plausible guesses at leading coordinates. The results obtained for different trial leading coordinates are rationalized through the calculation of a two-dimensional relaxed potential energy surface (PES) for the system. Analysis of the 2-D relaxed PES reveals that two of the trial pathways are entirely spurious, while two others contain useful information and can be used to furnish starting points for successful saddle-point searches. Implications for selection of trial leading coordinates in this class of proton chain transfer reactions are discussed, and a simple diagnostic function is proposed for revealing whether or not a relaxed pathway based on a trial leading coordinate is likely to furnish useful information.

Binomial leap methods for simulating stochastic chemical kinetics

This paper discusses efficient simulation methods for stochastic chemical kinetics. Based on the tau-leap and midpoint tau-leap methods of Gillespie [D. T. Gillespie, J. Chem. Phys. 115, 1716 (2001)], binomial random variables are used in these leap methods rather than Poisson random variables. The motivation for this approach is to improve the efficiency of the Poisson leap methods by using larger stepsizes. Unlike Poisson random variables whose range of sample values is from zero to infinity, binomial random variables have a finite range of sample values. This probabilistic property has been used to restrict possible reaction numbers and to avoid negative molecular numbers in stochastic simulations when larger stepsize is used. In this approach a binomial random variable is defined for a single reaction channel in order to keep the reaction number of this channel below the numbers of molecules that undergo this reaction channel. A sampling technique is also designed for the total reaction number of a reactant species that undergoes two or more reaction channels. Samples for the total reaction number are not greater than the molecular number of this species. In addition, probability properties of the binomial random variables provide stepsize conditions for restricting reaction numbers in a chosen time interval. These stepsize conditions are important properties of robust leap control strategies. Numerical results indicate that the proposed binomial leap methods can be applied to a wide range of chemical reaction systems with very good accuracy and significant improvement on efficiency over existing approaches. (C) 2004 American Institute of Physics.

Modelling of entrainment in industrial flotation cells: the effect of solids suspension

Solids concentration and particle size distribution gradually change in the vertical dimension of industrial flotation cells, subject primarily to the flotation cell size and design and the cell operating conditions. As entrainment is a two-step process and involves only the suspended solids in the top pulp region near the pulp-froth interface, the solids suspension characteristics have a significant impact on the overall entrainment. In this paper, a classification function is proposed to describe the state of solids suspension in flotation cells, similar to the definition of degree of entrainment for classification in the froth phase found in the literature. A mathematical model for solids suspension is also developed, in which the classification function is expressed as an exponential function of the particle size. Experimental data collected from three different Outokumpu tank flotation cells in three different concentrators are well fitted by the proposed exponential model. Under the prevailing experimental conditions, it was found that the solids content in the top region was relatively independent of cell operating conditions such as froth height and air rate but dependent on the cell size. Moreover, the results obtained from the total solids tend to be similar to those from a particular gangue mineral and hence may be applied to all minerals in entrainment calculation. (C) 2004 Elsevier Ltd. All rights reserved.

Chemical and structural disorder in eumelanins: A possible explanation for broadband absorbance

We report the results of an experimental and theoretical study of the electronic and structural properties of a key eumelanin precursor-5,6,-dihydroxyindole-2-carboxylic acid ( DHICA) - and its dimeric forms. We have used optical spectroscopy to follow the oxidative polymerization of DHICA to eumelanin and observe red shifting and broadening of the absorption spectrum as the reaction proceeds. First principles density functional theory calculations indicate that DHICA oligomers ( possible reaction products of oxidative polymerization) have the gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital red-shifted gaps with respect to the monomer. Furthermore, different bonding configurations ( leading to oligomers with different structures) produce a range of gaps. These experimental and theoretical results lend support to the chemical disorder model where the broadband monotonic absorption characteristic of all melanins is a consequence of the superposition of a large number of nonhomogeneously broadened Gaussian transitions associated with each of the components of a melanin ensemble. These results suggest that the traditional model of eumelanin as an amorphous organic semiconductor is not required to explain its optical properties and should be thoroughly reexamined. These results have significant implications for our understanding of the physics, chemistry, and biological function of these important biological macromolecules. Indeed, one may speculate that the robust functionality of melanins in vitro is a direct consequence of its heterogeneity, i.e., chemical disorder is a "low cost" natural resource in these systems

Finite mixture regression model with random effects: application to neonatal hospital length of stay

A two-component mixture regression model that allows simultaneously for heterogeneity and dependency among observations is proposed. By specifying random effects explicitly in the linear predictor of the mixture probability and the mixture components, parameter estimation is achieved by maximising the corresponding best linear unbiased prediction type log-likelihood. Approximate residual maximum likelihood estimates are obtained via an EM algorithm in the manner of generalised linear mixed model (GLMM). The method can be extended to a g-component mixture regression model with the component density from the exponential family, leading to the development of the class of finite mixture GLMM. For illustration, the method is applied to analyse neonatal length of stay (LOS). It is shown that identification of pertinent factors that influence hospital LOS can provide important information for health care planning and resource allocation. (C) 2002 Elsevier Science B.V. All rights reserved.

A novel laponite clay-based Fe nanocomposite and its photo-catalytic activity in photo-assisted degradation of Orange II

A novel laponite RD clay-based Fe nanocomposite (Fe-Lap-RD) has been successfully synthesized through a reaction between a solution of iron salt and an aqueous dispersion of laponite RD clay. The X-ray diffraction (XRD) results reveal that the Fe-Lap-RD mainly consists of Fe2O3 (maghemite) and Fe2Si4O10(OH)2 (iron silicate hydroxide), which have tetragonal and monoclinic structures, respectively, and has a high specific surface area as well as a high pore volume. The photo-catalytic activity of the Fe-Lap-RD was examined in the photo-assisted degradation of an organic azo dye Orange II. It was found that the mineralization of Orange 11 undergoes a slower kinetics than discoloration, and 70% total organic carbon (TOC) of 0.2 mM Orange 11 can be removed in 90 min, implying that the Fe-Lap-RD exhibited a high photo-catalytic activity in the presence of H2O2 and UV light (254 nm) in the photo-assisted degradation of Orange II. In addition, our experiments also illustrate that the Fe-Lap-RD has a long-term stability but is of low cost. This study illustrates the possibility of photo-assisted degradation of organic compounds without the requirements to remove the Fe ions after reaction. Two possible catalytic reaction mechanisms are also proposed. (C) 2003 Elsevier Science Ltd. All rights reserved.

General practitioners' role in preventive medicine: scenario analysis using alcohol as a case study

The purpose of this analysis is threefold: first, to extract from the literature, current levels of GP detection of at-risk drinking by their patients, rates at which general practitioners (GPs) offer an intervention; and the effectiveness of these interventions; secondly, to develop a model based on this literature to be used in conjunction with scenario analysis; and thirdly, to consider the cost implications of current efforts and various scenarios. This study deals specifically with Australian general practice. A two-step procedure is used in the scenario analysis, which involves identifying opportunities for detection, intervention, effectiveness and assigning probabilities to outcomes. The results suggest that increasing rates of GP intervention achieves greatest benefit and return on resource use. For every 5% point increase in the rate of GP intervention, an additional 26 754 at-risk drinkers modify their drinking behaviour at a cost of $231.45 per patient. This compares with a cost per patient modifying drinking behaviour of $232.60 and $208.31 for every 5% point increase in the rates of detection and effectiveness, respectively. The knowledge, skill and attitude of practitioners toward drinking are significant, and they can be the prime motivators in persuading their patients to modify drinking behaviour.