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.

Biomolecular engineering at interfaces

A broad review of technologically focused work concerning biomolecules at interfaces is presented. The emphasis is on developments in interfacial biomolecular engineering that may have a practical impact in bioanalysis, tissue engineering, emulsion processing or bioseparations. We also review methods for fabrication in an attempt to draw out those approaches that may be useful for product manufacture, and briefly review methods for analysing the resulting interfacial nanostructures. From this review we conclude that the generation of knowledge and-innovation at the nanoscale far exceeds our ability to translate this innovation into practical outcomes addressing a market need, and that significant technological challenges exist. A particular challenge in this translation is to understand how the structural properties of biomolecules control the assembled architecture, which in turn defines product performance, and how this relationship is affected by the chosen manufacturing route. This structure-architecture-process-performance (SAPP) interaction problem is the familiar laboratory scale-up challenge in disguise. A further challenge will be to interpret biomolecular self- and directed-assembly reactions using tools of chemical reaction engineering, enabling rigorous manufacturing optimization of self-assembly laboratory techniques. We conclude that many of the technological problems facing this field are addressable using tools of modem chemical and biomolecular engineering, in conjunction with knowledge and skills from the underpinning sciences. (c) 2005 Elsevier Ltd. All rights reserved.

Degradation of azo-dye Orange II by a photoassisted Fenton reaction using a novel composite of iron oxide and silicate nanoparticles as a catalyst

A novel nanocomposite of iron oxide and silicate, prepared through a reaction between a solution of iron salt and a dispersion of Laponite clay, was used as a catalyst for the photoassisted Fenton degradation of azo-dye Orange II. This catalyst is much cheaper than the Nafion-based catalysts, and our results illustrate that it can significantly accelerate the degradation of Orange II under the irradiation of UV light (lambda = 254 nm). An advantage of the catalyst is its long-term stability that was confirmed through using the catalyst for multiple runs in the degradation of Orange II. The effects of the H2O2 molar concentration, solution pH, wavelength and power of the LTV light, catalyst loading, and initial Orange II concentration on the degradation of Orange 11 were studied in detail. In addition, it was also found that discoloration of Orange 11 undergoes a faster kinetics than mineralization of Orange II and 75% total organic carbons of 0.1 mM Orange II can be eliminated after 90 min in the presence of 1.0 g of Fe-nanocomposite/L, 4.8 mM H2O2, and 1 x 8W UVC.

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.

Kinetics of reduction of lead smelting slags with solid carbon

Research techniques and a methodology have been developed that enable the reduction kinetics of molten lead smelting slags with solid carbon to be studied. The rates of reduction of PbO-FeO-Fe2O3-CaO-SiO2 slags with carbon have been measured for a range of slag compositions for PbO concentrations between 3 and 100 weight percent, and temperatures between 1423 and 1573 K. The reduction rates were determined for both graphite and coke. Within the range of process conditions examined, it has been shown that the reaction rates are almost independent of carbon reactivity, SiO2/CaO and SiO2/Fe ratio in the range of compositions investigated and are not influenced by the presence of sulphur in the slag.The apparent first order rate constants for oxygen removal increase with increasing PbO concentration and oxygen activity in the slag. The data indicate that the rate limiting reaction step for the reduction of lead slags with solid carbon is the chemical reaction at the gas/slag interface.

Kinetics of reduction of FeO from slag by graphite and coal chars

The rates of reduction of FeO from iron-saturated FeO-CaO-Al2O3-SiO2 slags by graphite, coke, bituminous coal and anthracitic coal chars at temperatures in the range 1 673-1873 K have been measured using a sessile drop technique. The extents of reaction were determined using EPMA analysis of quenched samples, and on line gas analysis using a quadrupole mass spectrometer. The reaction rates have been shown to be dependent critically on carbon type. For the reaction geometry used in this investigation the reduction rates of graphite and coke are observed to be faster than with coal chars. This unexpected finding is shown to be associated with differences in the dominant chemical and mass transfer mechanisms occurring at the reaction interface. High reaction rates are observed to occur with the formation of liquid Fe-C alloy product and the associated gasification of carbon from the alloy. The rates of reduction by coal chars are determined principally by the chemical reaction at the carbon/gas interface and slag phase mass transfer.

Structural and mechanical properties of nanostructured granular alumina catalysts

Granular gamma-Al2O3 support and 8 wt % CuO/gamma-Al2O3 catalyst were synthesized by a sol-gel granulation method. The pore structure, crush strength, hardness, and elasticity of these sol-gel-derived catalysts were studied and compared with similar commercial catalysts prepared by non-sol-gel methods. Alumina and CuO-coated alumina granular particles prepared by different methods have different macro- and microstructure. The sol-gel-derived granular gamma-alumina and CuO-coated gamma-alumina granular particles have a structure defined by compact packing of uniform, nanosized gamma-alumina crystallites. They are characterized by a more uniform pore size distribution and larger surface area as compared to similar commercial samples with a structure defined by packing of aggregates consisting of nonuniform gamma-alumina crystallites. Because of the differences in the macro- and microstructure, the sol-gel-derived granular samples offer higher crush strength and greater hardness than the commercial samples.

Symmetry contaminations in reactive scattering through long-lived collision complexes

We extend our Lanczos subspace time-independent wave packet method [J. Chem. Phys. 116 (2002) 2354] to investigate the issue of symmetry contaminations for the challenging deep-well H + O-2 reaction. Our central objective is to address the issue of whether significant symmetry contamination can occur if a wavepacket initially possessing the correct O-O exchange symmetry is propagated over tens of thousands of recursive steps using a basis which does not explicitly enforce the correct symmetry, and if so how seriously this affects the results. We find that symmetry contamination does exist where the symmetry constraint is not explicitly enforced in the basis. While it affects individual resonances and the associated peak amplitudes, the overall shape of the more averaged quantities such as total reaction probabilities and vibrational branching ratios are not seriously affected. (C) 2004 Elsevier B.V. All rights reserved.

Modern environmental improvement pathways for the coal power generation industry in Australia

In this work we assess the pathways for environmental improvement by the coal utilization industry for power generation in Australia. In terms of resources, our findings show that coal is a long term resource of concern as coal reserves are likely to last for the next 500 years or more. However, our analysis indicates that evaporation losses of water in power generation will approach 1000 Gl (gigalitres) per year, equivalent to a consumption of half of the Australian residential population. As Australia is the second driest continent on earth, water consumption by power generators is a resource of immediate concern with regards to sustainability. We also show that coal will continue to play a major role in energy generation in Australia and, hence, there is a need to employ new technologies that can minimize environmental impacts. The major technologies to reduce impacts to air, water and soils are addressed. Of major interest, there is a major potential for developing sequestration processes in Australia, in particular by enhanced coal bed methane (ECBM) recovery at the Bowen Basin, South Sydney Basin and Gunnedah Basin. Having said that, CO2 capture technologies require further development to support any sequestration processes in order to comply with the Kyoto Protocol. Current power generation cycles are thermodynamic limited, with 35-40% efficiencies. To move to a high efficiency cycle, it is required to change technologies of which integrated gasification combined cycle plus fuel cell is the most promising, with efficiencies expected to reach 60-65%. However, risks of moving towards an unproven technology means that power generators are likely to continue to use pulverized fuel technologies, aiming at incremental efficiency improvements (business as usual). As a big picture pathway, power generators are likely to play an increasing role in regional development; in particular EcoParks and reclaiming saline water for treatment as pressures to access fresh water supplies will significantly increase.

Reconsidering the use of photosynthetic bacteria for removal of sulfide from wastewater

The feasibility of using photosynthetic sulfide-oxidizing bacteria to remove sulfide from wastewater in circumstances where axenic cultures are unrealistic has been completely reconsidered on the basis of known ecophysiological data, and the principles of photobioreactor and chemical reactor engineering. This has given rise to the development of two similar treatment concepts relying on biofilms dominated by green sulfur bacteria (GSB) that develop on the exterior of transparent surfaces suspended in the wastewater. The GSB are sustained and selected for by radiant energy in the band 720 - 780 nm, supplied from within the transparent surface. A model of one of these concepts was constructed and with it the reactor concept was proven. The dependence of sulfide-removal rate on bulk sulfide concentration has been ascertained. The maximum net areal sulfide removal rate was 2.23 g m(-2) day(-1) at a bulk sulfide concentration of 16.5 mg L-1 and an incident irradiance of 1.51 W m(-2). The system has a demonstrated capacity to mitigate surges in sulfide load, and appears to use much less radiant power than comparable systems. The efficacy with which this energy was used for sulfide removal was 1.47 g day(-1) W-1. The biofilm was dominated by GSB, and evidence gathered indicated that other types of phototrophs were not present. (C) 2004 Wiley Periodicals, Inc.

Diffusion of linear paraffins in nanoporous silica

The diffusion of hexane, heptane, octane, and decane in nanoporous MCM-41 silica at various temperatures is investigated by the zero-length-column method. The diffusion coefficients are derived by a complete-time-range analysis of desorption curves at different purge flow rates and temperatures. The results show that the calculated low-coverage diffusivity values decrease monotonically, and the derived Henry's law constants increase, as the carbon number of paraffins increases. The study reveals that transport is strongly influenced by intracrystalline diffusion and dominated by the sorbate-sorbent interaction. The diffusion activation energy and adsorption isosteric heat at zero loading increase monotonically with the carbon number of linear paraffins, but their ratio is essentially constant for each adsorbate compound.

Effect of metal chlorides on the sintering and densification of hydroxyapatite adsorbent

This work is part of a series of studies dealing with the evaluation of the effects of major elements of solid waste, especially metallic oxides, nitrates, sulfates, and chlorides, on the sintering and the densification of calcium hydroxyapatite (Ca-HAP) adsorbent. The effects of chloride salts of potassium (KCl) and zinc (ZnCl2) on sintering and densification of Ca-HAP were studied using surface area reduction and shrinkage measurements. The addition of KCl (2% w/w) activated the sintering process by bringing a swift reduction in surface area and lowering the densification temperature. However, a low final densification was achieved. Increasing the amount of this additive to 10% w/w further lowered the final densification and lowered the densification temperature of hydroxyapatite by 150 degrees C. On the other hand, the addition of 2 wt % of ZnCl2 deactivated the sintering process by slowing down the densification process and raising the densification temperature. However, the reduction of surface area was comparable to that of Ca-HAP. The densification rate contained two or more rate maxima indicating the additives (salts) bring multiple speeds in the densification process.

Tetra- and tribromopbenoxyanisoles in Marine Samples from Oceania

Some methoxylated polybrominated diphenyl ethers (MeO-BDEs) are known halogenated natural products (HNPs) and are frequently detected in higher organisms of the marine environment. In this study we demonstrate that a prominent MeO-BDE, previously detected in marine mammals from Australia, is identical to 3,5-dibromo-2-(2',4'-dibromo)phenoxyanisole(BC-3,6-MeO-BDE47). Up to 1.9mg/ kg of 6-MeO-BDE 47 was present in cetaceans from Australia, 0.2-0.3 mg/kg in two crocodile eggs from Australia, but concentrations of 1 or 2 orders of magnitude lower were found in shark liver oil from New Zealand and in marine mammals from Africa and the Antarctic. Concentrations of 6-MeO-BDE47 in samples from Australia were in the same range as anthropogenic pollutants such as PCB 153 and p,p'-DDE. Along with 6-MeO-BDE 47 and the known HNP 4,6-dibromo-2-(2',4'-dibromo)phenoxyanisole (BC-2,2'-MeO-BDE 68), several tribromophenoxyanisoles (MeO-triBDE) were present in tissue of Australian cetaceans. To determine their structure, abiotic debromination experiments were performed using 6-MeO-BDE 47 and 2'-MeO-BDE 68 and superreduced di cyanocobalamine. These experiments resulted in formation of eight MeO-triBDEs, all of which were detected in the cetacean samples. Five of these eight MeO-triBDEs could be identified based on two standard compounds as well as gas chromatographic and mass spectrometric features. It was also shown that the first eluting isomer (compound 1), 6-MeO-BDE 17 (compound 2), and 2-MeO-BDE 39 (compound 5) were the most prominent MeO-triBDEs in the Australian cetacean samples. The concentrations of the MeO-triBDEs in two cetacean samples were 0.20 and 0.36 mg/kg, respectively. Although the reductive debromination with dicyanocobalamine resulted in a different congener pattern than was found in the marine mammals, it could not be excluded that the tribromo congeners of 6-MeO-BDE 47 and 2'-MeO-BDE 68 in the samples were metabolites of the latter.

Bounding the stability and rupture condition of emulsion and foam films

A scaling law is presented that provides a complete solution to the equations bounding the stability and rupture of thin films. The scaling law depends on the fundamental physicochemical properties of the film and interface to calculate bounds for the critical thickness and other key film thicknesses, the relevant waveforms associated with instability and rupture, and film lifetimes. Critical thicknesses calculated from the scaling law are shown to bound the values reported in the literature for numerous emulsion and foam films. The majority of critical thickness values are between 15 to 40% lower than the upper bound critical thickness provided by the scaling law.