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.

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 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.

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.

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.

New calibration technique for multiple-component stress wave force balances

Force measurement in hypervelocity expansion tubes is not possible using conventional techniques. The stress wave force balance technique can be applied in expansion tubes to measure forces despite the short test times involved. This paper presents a new calibration technique for multiple-component stress wave force balances where an impulse response created using a load distribution is required and no orthogonal surfaces on the model exist.. This new technique relies on the tensorial superposition of single-component impulse responses analogous to the vectorial superposition of the calibration loads. The example presented here is that of a scale model of the Mars Pathfinder, but the technique is applicable to any geometry and may be useful for cases where orthogonal loads cannot be applied.

High-Frequency Turbulence and Sediment Flux Measurements in an Upper Estuarine Zone

In natural estuaries, the predictions of scalar dispersion are rarely predicted accurately because of a lack of fundamental understanding of the turbulence structure in estuaries. Herein detailed turbulence field measurements were conducted continuously at high frequency for 50 hours in the upper zone of a small subtropical estuary with semi-diurnal tides. Acoustic Doppler velocimetry was deemed the most appropriate measurement technique for such shallow water depths (less than 0.4 m at low tides), and a thorough post-processing technique was applied. In addition, some experiments were conducted in laboratory under controlled conditions using water and soil samples collected in the estuary to test the relationship between acoustic backscatter strength and suspended sediment load. A striking feature of the field data set was the large fluctuations in all turbulence characteristics during the tidal cycle, including the suspended sediment flux. This feature was rarely documented.

Isotopes, wool, and rangeland monitoring: Let the sheep do the sampling

Stable carbon isotope analyses of wool staples provided insight into the vegetation consumed by sheep at a temporal resolution not previously studied. Contemporary Australian and historic South African samples dating back to 1916 were analyzed for their stable carbon isotope ratio, a proxy for the proportion of C-3 and C-4 plant species consumed by animals. Sheep sample vegetation continuously throughout a year, and as their wool grows it integrates and stores information about their diet. In subtropical and tropical rangelands the majority of grass species are C-4. Since sheep prefer to graze, and their wool is an isotopic record of their diet, we now have the potential to develop a high resolution index to the availability of grass from a sheep's perspective. Isotopic analyses of wool suggest a new direction for monitoring grazing and for the reconstruction of past vegetation changes, which will make a significant contribution to traditional rangeland ecology and management. It is recommended that isotopic and other analyses of wool be further developed for use in rangeland monitoring programs to provide valuable feedback for land managers.

CO2 reforming of methane on Ni catalysts: Effects of the support phase and preparation technique

The effects of the support phase and catalyst preparation methods on catalytic activity and carbon deposition were systematically investigated over nickel catalysts supported on Al2O3, SiO2 and MgO for the reforming reaction of methane with carbon dioxide. It is found that the pore structure of the support and metal-support interaction significantly affected the catalytic activity and coking resistance. Catalyst with well-developed porosity exhibited higher catalytic activity. Strong interaction between metal and the support made the catalyst more resistant to sintering and coking, thus resulting in a longer time of catalyst stability. (C) 1998 Elsevier Science B.V.

Effect of load distance on self-selected manual lifting technique

The aim of this research is to determine the effects of constraining the horizontal distance of the feet from the load on the posture adopted at the start of the lift. Kinematic data were collected while each of 24 subjects lifted 3, 6, and 9 kg loads from a starting height 18 cm above the ground. The position of the feet was controlled relative to the load such that the horizontal distance from the hand to the ankle at the start of extension was either 20, 40, or 60 cm. Subjects performed 20 trials in each of six combinations of load and ankle-load distance chosen to provide three sets of equivilent load moment pairs. The initial horizontal distance from the load to the ankle had a large influence on the posture adopted to lift the load. Ankle and knee flexion, in particular, were reduced when the ankle-load distance was smaller, and particularly so when the distance was reduced to 20 cm. Hip flexion was reduced to a smaller extent, while lumbar vertebral flexion remained relatively unchanged. The inclination of the trunk at the start of the lift was unchanged when the ankle-load distance was 60 or 40 cm, but was 10 degrees greater when the load was 20 cm from the ankles, indicating that subjects adopted a posture closer to a stoop when the ankle-load distance was small. Comparison of conditions of equal load moment (but different load mass and ankle-load distance) revealed differences which mirrored the effects of ankle-load distance alone, suggesting that the effects of ankle-load distance on the posture adopted at the start of extension were largely independent of the load moment. While the forces and torques required to lift a load must be to some extent dependent on the load moment, rather than load or ankle-load distance per se, the posture adopted to lift the load is not.

Bacteria in shrimp pond sediments: their role in mineralizing nutrients and some suggested sampling strategies

Strategies for sampling sediment bacteria were examined in intensive shrimp, Penaeus monodon (Fabricius), ponds in tropical Australia. Stratified sampling of bacteria at the end of the production season showed that the pond centre, containing flocculated sludge, had significantly higher bacterial counts (15.5 X 10(9) g(-1) dw) than the pond periphery (8.1 X 10(9) g(-1) dw), where the action of aerators had swept the pond floor. The variation in bacterial counts between these two zones within a pond was higher than that between sites within each zone or between ponds. Therefore, sampling effort should be focused within these zones: for example, sampling two ponds at six locations within each of the two zones resulted in a coefficient of variation of approximate to 5%. Bacterial numbers in the sediment were highly correlated with sediment grain size, probably because eroded soil particles and organic waste both accumulated in the centre of the pond. Despite high inputs of organic matter added to the ponds, principally as pelleted feeds, the mean bacterial numbers and nutrient concentrations (i.e. organic carbon, nitrogen and phosphorus) in the sediment were similar to those found in mangrove sediments. This suggests that bacteria are rapidly remineralizing particulates into soluble compounds. Bacterial numbers were highly correlated with organic carbon and total kjeldahl nitrogen in the sediment, suggesting that these were limiting factors to bacterial growth.