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.

A wavelet-based adaptive technique for adsorption problems involving steep gradients

A general, fast wavelet-based adaptive collocation method is formulated for heat and mass transfer problems involving a steep moving profile of the dependent variable. The technique of grid adaptation is based on sparse point representation (SPR). The method is applied and tested for the case of a gas–solid non-catalytic reaction in a porous solid at high Thiele modulus. Accurate and convergent steep profiles are obtained for Thiele modulus as large as 100 for the case of slab and found to match the analytical solution.

Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

Analysis of multiple gas solid reactions

Multiple gas solid reactions involving one solid and N gaseous reactants are investigated in this study by using a matched asymptotic expansion technique. Two cases are particularly studied. In the first case all N chemical reaction rates are faster than the diffusion rate. While in the second case only M (M < N) chemical reaction rates are faster than the diffusion rate and the rates of the remaining (N-M) chemical reactions are comparable to that of diffusion. For these two cases the solid concentration profile behaves like a travelling wave. In the first case the wave front velocity is contributed linearly by all gaseous reactants (additive law) while in the second case this law does not hold.

On the validity of the shrinking core model in noncatalytic gas solid reaction

By using a matched asymptotic expansion technique, the shrinking core model (SCM) used in non-catalytic gas solid reactions with general kinetic expression is rigorously justified in this paper as a special case of the homogeneous model when the reaction rate is much faster than that of diffusion. The time-pendent velocity of the moving reacted-unreacted interface is found to be proportional to the gas flux at that interface for all geometries of solid particles, and the thickness order of the reaction zone and also the degree of chemical reaction at the interface is discussed in this paper.

Transient response of diffusion cell containing a porous solid

Transient response of an adsorbing or non-adsorbing tracer injected as step or square pulse input in a diffusion cell with two flowing streams across the pellet is theoretically investigated in this paper. Exact solutions and the asymptotic solutions in the time domain and in three different limits are obtained by using an integral transform technique and a singular perturbation technique, respectively. Parametric dependence of the concentrations in the top and bottom chambers can be revealed by investigating the asymptotic solutions, which are far simpler than their exact counterpart. In the time domain investigation, it is found that the bottom-chamber concentration is very sensitive to the value of the macropore effective diffusivity. Therefore this concentration could be used to extract diffusivity by fitting in the time domain. The bottom-chamber concentration is also sensitive to flow rate, pellet length chamber volume and the type of input (step and square input).

Growth of subtropical ECM fungi with different nitrogen sources using a new floating culture technique

Eight species of ectomycorrhizal (ECM) fungi in the genera Amanita. Gymnoboletus, Lactarius, and Russula were isolated from subtropical plant communities in eastem Australia. Two species were isolated from each of rainforest, Nothofagus forest, Eucalyptus forest, and Eucalyptus dominated wallum (heath) forest. These communities differ strongly in their soluble soil nitrogen (N) composition. The ability of the fungi to use inorganic (nitrate, ammonium) and organic (amide, peptide, protein) nitrogen sources was determined. As the fungi did not grow in liquid culture, a 'floating culture' technique was devised that allows hyphal growth on a screen floating on liquid medium. With some exceptions, fungal biomass production in floating culture closely reflected fungal growth on solid media assessed by total colony glucosamine content. Most isolates grown in floating culture had similar glucosamine concentrations on all N sources, with isolate specific concentrations ranging from 6 to 12 mug glucosamine g(-1) DW. However, Russula spp. had up to 1.7-fold higher glucosamine concentrations when growing with glutamine or ammonium compared to nitrate, glutathione or protein. Floating cultures supplied with 0.5, 1.5. 4.5, or 10 mm N mostly produced greatest biomass with 4.5 mM N. In vitro nitrate reductase activity (NRA) ranged from very low (0.03 mumol NO2- g(-1) fw h(-1)) in Russula sp. (wallum) to high (2.16 mumol NO2- g(-1) fw h(-1)) in Gymnoboletus sp. (rainforest) and mirrored the fungi's ability to use nitrate as a N source. All Russula spp. (wallum, Nothofagus and Eucalyptus forests), Lactarills sp, (rainforest) and.4manita sp. (wallum) utilized ammonium and glutamine but had little ability to use other N sources. In contrast,Amanita species (Nothofagus and Eucalyptus forests) grew on all N sources but produced most biomass with ammonium and glutamine. Only Gymnoboletus sp. (rainforest) showed similar growth with nitrate and ammonium as N sources. Fungal N source use was not associated with taxonomic groups, but is discussed in the context of soil N sources in the different habitats.

Axial mixing in rotating drums using magnetic resonance imaging using bran as a model for solid state fermentations

Magnetic resonance imaging (MRI) is an easily automated, reliable technique to investigate axial mixing within rotating drums. Moist bran can be clearly differentiated from dry bran using MRI allowing a non-segregating tracer for axial mixing. For a 20-cm diameter drum, the axial dispersion coefficient in the particle bed was 0.51 cm s(-2). Axial dispersion is scale-dependent.

A method of sampling blood from Australian freshwater turtles

Blood sampling is an essential technique in many herpetological studies. This paper describes a quick and humane technique to collect blood samples from three species of Australian chelid turtles ( Order Pleurodira): Chelodina expansa, Elseya latisternum, and Emydura macquarii signata.

Consequences of incorrect sampling procedures in resonance-based radar target identification

Radar target identification based on complex natural resonances is sometimes achieved by convolving a linear time-domain filter with a received target signature. The filter is constructed from measured or pre-calculated target resonances. The performance of the target identification procedure is degraded if the difference between the sampling rates of the target signature and the filter is ignored. The problem is investigated for the natural extinction pulse technique (E-pulse) for the case of identifying stick models of aircraft.

Clinical pharmacokinetics and pharmacodynamics of tacrolimus in solid organ transplantation

The aim of this review is to analyse critically the recent literature on the clinical pharmacokinetics and pharmacodynamics of tacrolimus in solid organ transplant recipients. Dosage and target concentration recommendations for tacrolimus vary from centre to centre, and large pharmacokinetic variability makes it difficult to predict what concentration will be achieved with a particular dose or dosage change. Therapeutic ranges have not been based on statistical approaches. The majority of pharmacokinetic studies have involved intense blood sampling in small homogeneous groups in the immediate post-transplant period. Most have used nonspecific immunoassays and provide little information on pharmacokinetic variability. Demographic investigations seeking correlations between pharmacokinetic parameters and patient factors have generally looked at one covariate at a time and have involved small patient numbers. Factors reported to influence the pharmacokinetics of tacrolimus include the patient group studied, hepatic dysfunction, hepatitis C status, time after transplantation, patient age, donor liver characteristics, recipient race, haematocrit and albumin concentrations, diurnal rhythm, food administration, corticosteroid dosage, diarrhoea and cytochrome P450 (CYP) isoenzyme and P-glycoprotein expression. Population analyses are adding to our understanding of the pharmacokinetics of tacrolimus, but such investigations are still in their infancy. A significant proportion of model variability remains unexplained. Population modelling and Bayesian forecasting may be improved if CYP isoenzymes and/or P-glycoprotein expression could be considered as covariates. Reports have been conflicting as to whether low tacrolimus trough concentrations are related to rejection. Several studies have demonstrated a correlation between high trough concentrations and toxicity, particularly nephrotoxicity. The best predictor of pharmacological effect may be drug concentrations in the transplanted organ itself. Researchers have started to question current reliance on trough measurement during therapeutic drug monitoring, with instances of toxicity and rejection occurring when trough concentrations are within 'acceptable' ranges. The correlation between blood concentration and drug exposure can be improved by use of non-trough timepoints. However, controversy exists as to whether this will provide any great benefit, given the added complexity in monitoring. Investigators are now attempting to quantify the pharmacological effects of tacrolimus on immune cells through assays that measure in vivo calcineurin inhibition and markers of immuno suppression such as cytokine concentration. To date, no studies have correlated pharmacodynamic marker assay results with immunosuppressive efficacy, as determined by allograft outcome, or investigated the relationship between calcineurin inhibition and drug adverse effects. Little is known about the magnitude of the pharmacodynamic variability of tacrolimus.

Aquatic passive sampling of herbicides on naked particle loaded membranes: Accelerated measurement and empirical estimation of kinetic parameters

Water-sampler equilibrium partitioning coefficients and aqueous boundary layer mass transfer coefficients for atrazine, diuron, hexazionone and fluometuron onto C18 and SDB-RPS Empore disk-based aquatic passive samplers have been determined experimentally under a laminar flow regime (Re = 5400). The method involved accelerating the time to equilibrium of the samplers by exposing them to three water concentrations, decreasing stepwise to 50% and then 25% of the original concentration. Assuming first-order Fickian kinetics across a rate-limiting aqueous boundary layer, both parameters are determined computationally by unconstrained nonlinear optimization. In addition, a method of estimation of mass transfer coefficients-therefore sampling rates-using the dimensionless Sherwood correlation developed for laminar flow over a flat plate is applied. For each of the herbicides, this correlation is validated to within 40% of the experimental data. The study demonstrates that for trace concentrations (sub 0.1 mu g/L) and these flow conditions, a naked Empore disk performs well as an integrative sampler over short deployments (up to 7 days) for the range of polar herbicides investigated. The SDB-RPS disk allows a longer integrative period than the C18 disk due to its higher sorbent mass and/or its more polar sorbent chemistry. This work also suggests that for certain passive sampler designs, empirical estimation of sampling rates may be possible using correlations that have been available in the chemical engineering literature for some time.

Thermal stability analysis of organo-silicates, using solid phase microextraction techniques

An analysis of thermal degradation products evolved during the melt processing of organo-layered silicates (OLS) was carried out via the use of a solid phase microextraction (SPME) technique. Two commerical OLSs and one produced in-house were prepared for comparision. The solid phase microextraction technique proved to be a very effective technique for investigating the degradation of the OLS at a specific processing temperature. The results showed that most available OLSs will degrade under typical conditions required for the melt processing of many polymers, including thermoplastic polyurethanes. It is suggested that these degradation products may lead to changes in the structure and properties of the final polymer, particularly in thermoplastic polyurethanes, which seem significantly succeptable to the presence of these products. It is also suggested that many commercially available OLSs are produced in such a way that results in an excess of unbound organic modifier, giving rise to a greater quantity of degradation products. All OLSs where compared and characterised by TGA and GC-MS. (c) 2004 Elsevier B.V. All rights reserved.