Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

Centile reference for total energy expenditure in infants from 1 to 12 months

Background: A knowledge of energy expenditure in infancy is required for the estimation of recommended daily amounts of food energy, for designing artificial infant feeds, and as a reference standard for studies of energy metabolism in disease states. Objectives: The objectives of this study were to construct centile reference charts for total energy expenditure (TEE) in infants across the first year of life. Methods: Repeated measures of TEE using the doubly labeled water technique were made in 162 infants at 1.5, 3, 6, 9 and 12 months. In total, 322 TEE measurements were obtained. The LMS method with maximum penalized likelihood was used to construct the centile reference charts. Centiles were constructed for TEE expressed as MJ/day and also expressed relative to body weight (BW) and fat-free mass (FFM). Results: TEE increased with age and was 1.40,1.86, 2.64, 3.07 and 3.65 MJ/day at 1.5, 3, 6, 9 and 12 months, respectively. The standard deviations were 0.43, 0.47, 0.52, 0.66 and 0.88, respectively. TEE in MJ/kg increased from 0.29 to 0.36 and in MJ/day/kg FFM from 0.36 to 0.48. Conclusions: We have presented centile reference charts for TEE expressed as MJ/day and expressed relative to BW and FFM in infants across the first year of life. There was a wide variation or biological scatter in TEE values seen at all ages. We suggest that these centile charts may be used to assess and possibly quantify abnormal energy metabolism in disease states in infants.

A molecular phylogeny of the Australian skink genera Eulamprus, Gnypetoscincus and Nangura

Skinks from the genera Eulamprus, Gnypetoscincus and Nangura are a prominent component of the reptile fauna of the mesic forests of the east coast of Australia and have been the subject of numerous ecological studies. Highly conserved morphology and the retention of ancestral traits have limited our understanding of the relationships within and among these genera beyond an initial identification of species groups within Eulamprus. To address this deficit and to explore the relationships between Eulamprus and the monotypic genera Nangura and Gnypetoscincus, sections of two mitochondrial genes (ND4 and 16S rRNA) were sequenced and subjected to Bayesian phylogenetic analysis. This phylogenetic analysis supports recognition of the three species groups proposed for Eulamprus (murrayi, quoyii and tenuis) and indicates that this genus is paraphyletic, with Gnypetoscincus and Nangura being proximal to basal lineages of the tenuis group. To resolve these and broader problems of paraphyly, we suggest that each of the species groups from 'Eulamprus' should be recognised as a distinct genus. The phylogenetically and ecologically distinct water skinks of the quoyii group would be retained within Eulamprus and the diverse species of the tenuis group allocated to Concinnia. We suggest placing the monophyletic murrayi group, endemic to the rainforests of central eastern Australia, in a new genus ( yet to be formally described). The sequencing data also revealed the existence of a genetically divergent but morphologically cryptic lineage within E. murrayi and substantial diversity within E. quoyii. There is evidence for two major habitat shifts from rainforest towards drier habitats, one leading to the quoyii group and the second defining a clade of three species within the tenuis complex. These ecological transitions may represent adaptations to general drying across eastern Australia during the late Miocene - Pliocene. Each of the major areas of east coast tropical or subtropical rainforest contains multiple phylogenetically diverse endemic species, reflecting the long-term persistence and high conservation value of wet forest habitats in each area.

A practical method of estimating energy expenditure during tennis play

This study aimed to develop a practical method of estimating energy expenditure (EE) during tennis. Twenty-four elite female tennis players first completed a tennis-specific graded test in which five different intensity levels were applied randomly. Each intensity level was intended to simulate a game of singles tennis and comprised six 14 s periods of activity alternated with 20 s of active rest. Oxygen consumption (VO2) and heart rate (HR) were measured continuously and each player's rate of perceived exertion (RPE) was recorded at the end of each intensity level. Rate of energy expenditure (EEVO2) during the test was calculated using the sum of VO2 during play and the 'O-2 debt' during recovery, divided by the duration of the activity. There were significant individual linear relationships between EEVO2 and RPE, EEVO2 and HR, (rgreater than or equal to0.89 rgreater than or equal to0.93; p

The peroxisome proliferator-activated receptor beta/delta agonist, GW501516, regulates the expression of genes involved in lipid catabolism and energy uncoupling in skeletal muscle cells

Lipid homeostasis is controlled by the peroxisome proliferator-activated receptors (PPARalpha, -beta/delta, and -gamma) that function as fatty acid-dependent DNA-binding proteins that regulate lipid metabolism. In vitro and in vivo genetic and pharmacological studies have demonstrated PPARalpha regulates lipid catabolism. In contrast, PPARgamma regulates the conflicting process of lipid storage. However, relatively little is known about PPARbeta/delta in the context of target tissues, target genes, lipid homeostasis, and functional overlap with PPARalpha and -gamma. PPARbeta/delta, a very low-density lipoprotein sensor, is abundantly expressed in skeletal muscle, a major mass peripheral tissue that accounts for approximately 40% of total body weight. Skeletal muscle is a metabolically active tissue, and a primary site of glucose metabolism, fatty acid oxidation, and cholesterol efflux. Consequently, it has a significant role in insulin sensitivity, the blood-lipid profile, and lipid homeostasis. Surprisingly, the role of PPARbeta/delta in skeletal muscle has not been investigated. We utilize selective PPARalpha, -beta/delta, -gamma, and liver X receptor agonists in skeletal muscle cells to understand the functional role of PPARbeta/delta, and the complementary and/or contrasting roles of PPARs in this major mass peripheral tissue. Activation of PPARbeta/delta by GW501516 in skeletal muscle cells induces the expression of genes involved in preferential lipid utilization, beta-oxidation, cholesterol efflux, and energy uncoupling. Furthermore, we show that treatment of muscle cells with GW501516 increases apolipoprotein-A1 specific efflux of intracellular cholesterol, thus identifying this tissue as an important target of PPARbeta/delta agonists. Interestingly, fenofibrate induces genes involved in fructose uptake, and glycogen formation. In contrast, rosiglitazone-mediated activation of PPARgamma induces gene expression associated with glucose uptake, fatty acid synthesis, and lipid storage. Furthermore, we show that the PPAR-dependent reporter in the muscle carnitine palmitoyltransferase-1 promoter is directly regulated by PPARbeta/delta, and not PPARalpha in skeletal muscle cells in a PPARgamma coactivator-1-dependent manner. This study demonstrates that PPARs have distinct roles in skeletal muscle cells with respect to the regulation of lipid, carbohydrate, and energy homeostasis. Moreover, we surmise that PPARgamma/delta agonists would increase fatty acid catabolism, cholesterol efflux, and energy expenditure in muscle, and speculate selective activators of PPARbeta/delta may have therapeutic utility in the treatment of hyperlipidemia, atherosclerosis, and obesity.

Rates of electronic energy transfer in conformationally flexible bichromophoric macrocyclic complexes: a combined experimental and molecular modeling study

Electronic energy transfer (EET) rate constants between a naphthalene donor and anthracene acceptor in [ZnL4a](ClO4)(2) and [ZnL4b](ClO4)(2) were determined by time-resolved fluorescence where L-4a and L-4b are the trans and cis isomers of 6-((anthracen-9-yl-methyl)amino)-6,13-dimethyl-13-((naphthalen-1-yl-methyl)amino)-1,4,8,11-tetraazacyclotetradecane, respectively. These isomers differ in the relative disposition of the appended chromophores with respect to the macrocyclic plane. The trans isomer has an energy transfer rate constant (k(EET)) of 8.7 x 10(8) s(-1), whereas that of the cis isomer is significantly faster (2.3 x 10(9) s(-1)). Molecular modeling was used to determine the likely distribution of conformations in CH3CN solution for these complexes in an attempt to identify any distance or orientation dependency that may account for the differing rate constants observed. The calculated conformational distributions together with analysis by H-1 NMR for the [ZnL4a](2+) trans complex in the common trans-III N-based isomer gave a calculated Forster rate constant close to that observed experimentally. For the [ZnL4b](2+) cis complex, the experimentally determined rate constant may be attributed to a combination of trans-Ill and trans-I N-based isomeric forms of the complex in solution.

High energy radiation grafting of fluoropolymers

Fluoropolymers are known as chemically inert materials with good high temperature resistance, so they are often the materials of choice for harsh chemical environments. These properties arise because the carbon-fluorine bond is the strongest of all bonds between other elements and carbon, and, because of their large size, fluorine atoms can protect the carbon backbone of polymers such as poly(tetrafluoroethylene), PTFE, from chemical attack. However, while the carbon-fluorine bond is much stronger than the carbon hydrogen bond, the G values for radical formation on high energy radiolysis of fluoropolymers are roughly comparable to those of their protonated counterparts. Thus, efficient high energy radiation grafting of fluoropolymers is practical, and this process can be used to modify either the surface or bulk properties of a fluoropolymer. Indeed, radiation grafted fluoropolymers are currently being used as separation membranes for fuel cells, hydrophilic filtration membranes and matrix substrate materials for use in combinatorial chemistry. Herein we present a review of recent studies of the high energy radiation grafting of fluoropolymers and of the analytical methods available to characterize the grafts. (C) 2003 Elsevier Ltd. All rights reserved.

Activation energy distribution of thermal annealing of a bituminous coal

A bituminous coal was pyrolyzed in a nitrogen stream in an entrained flow reactor at various temperatures from 700 to 1475 degreesC. Char samples were collected at different positions along the reactor. Each collected sample was oxidized nonisothermally in a TGA for reactivity determination. The reactivity of the coal char was found to decrease rapidly with residence time until 0.5 s, after which it decreased only slightly. On the bases of the reactivity data at various temperatures, a new approach was utilized to obtaining the true activation energy distribution function for thermal annealing without the assumption of any distribution function form or a constant preexponential factor. It appears that the true activation energy distribution function consists of two separate parts corresponding to different temperature ranges, suggesting different mechanisms in different temperature ranges. Partially burnt coal chars were also collected along the reactor when the coal was oxidized in air at various temperatures from 700 to 1475 degreesC. The collected samples were analyzed for the residual carbon content and the specific reaction rate was estimated. The characteristic time of thermal deactivation was compared with that of oxidation under realistic conditions. The characteristic times were found to be close to each other, indicating the importance of thermal deactivation during combustion of the coal studied.