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.

Effect of vehicle pretreatment on the flux, retention, and diffusion of topically applied penetrants in vitro

Purpose. The flux of a topically applied drug depends on the activity in the skin and the interaction between the vehicle and skin. Permeation of vehicle into the skin can alter the activity of drug and the properties of the skin barrier. The aim of this in vitro study was to separate and quantify these effects. Methods. The flux of four radiolabeled permeants (water, phenol, diflunisal, and diazepam) with log K-oct/water values from 1.4 to 4.3 was measured over 4 h through heat-separated human epidermis pretreated for 30 min with vehicles having Hildebrand solubility parameters from 7.9 to 23.4 (cal/cm(3))(1/2). Results. Enhancement was greatest after pretreatment with the more lipophilic vehicles. A synergistic enhancement was observed using binary mixtures. The flux of diazepam was not enhanced to the same extent as the other permeants, possibly because its partitioning into the epidermis is close to optimal (log K-oct 2.96). Conclusion. An analysis of the permeant remaining in the epidermis revealed that the enhancement can be the result of either increased partitioning of permeant into the epidermis or an increasing diffusivity of permeants through the epidermis.

Surface diffusion of strongly adsorbing vapors in activated carbon by a differential permeation method

Conventional methods to determine surface diffusion of adsorbed molecules are proven to be inadequate for strongly adsorbing vapors on activated carbon. Knudsen diffusion permeability (B-k) for strongly adsorbing vapors cannot be directly estimated from that of inert gases such as helium. In this paper three models are considered to elucidate the mechanism of surface diffusion in activated carbon. The transport mechanism in all three models is a combination of Knudsen diffusion, viscous flow and surface diffusion. The collision reflection factor f (which is the fraction of molecules undergoing collision to the solid surface over reflection from the surface) of the Knudsen diffusivity is assumed to be a function of loading. It was found to be 1.79 in the limit of zero loading, and decreases as loading increases. The surface diffusion permeability increases sharply at very low pressures and then starts to decrease after it has reached a maximum (B(mum)s) at a threshold pressure. The initial rapid increase in the total permeability is mainly attributed to surface diffusion. Interestingly the B(mum)s for all adsorbates appear at the same volumetric adsorbed phase concentration, suggesting that the volume of adsorbed molecules may play an important role in the surface diffusion mechanism in activated carbon. (C) 2003 Elsevier Ltd. All rights reserved.

Hydrodynamic origin of diffusion in nanopores

We study the transport of a subcritical Lennard-Jones fluid in a cylindrical nanopore, using a combination of equilibrium and nonequilibrium as well as dual control volume grand canonical molecular dynamics methods. We show that all three techniques yield the same value of the transport coefficient for diffusely reflecting pore walls, even in the presence of viscous transport. We also demonstrate that the classical Knudsen mechanism is not manifested, and that a combination of viscous flow and momentum exchange at the pore wall governs the transport over a wide range of densities.

Radial Distortion Self-Calibration

In cameras with radial distortion, straight lines in space are in general mapped to curves in the image. Although epipolar geometry also gets distorted, there is a set of special epipolar lines that remain straight, namely those that go through the distortion center. By finding these straight epipolar lines in camera pairs we can obtain constraints on the distortion center(s) without any calibration object or plumbline assumptions in the scene. Although this holds for all radial distortion models we conceptually prove this idea using the division distortion model and the radial fundamental matrix which allow for a very simple closed form solution of the distortion center from two views (same distortion) or three views (different distortions). The non-iterative nature of our approach makes it immune to local minima and allows finding the distortion center also for cropped images or those where no good prior exists. Besides this, we give comprehensive relations between different undistortion models and discuss advantages and drawbacks.

Agronomic performance of 27 cocoa progenies and plant selection based on productivity, self-compatibility and disease resistance

The objective of the present work was to evaluate 27 progenies of cocoa crosses considering the agronomic traits and select F1 plants within superior crosses. The experiment was installed in March 2005, in the Experimental Station Joaquim Bahiana (ESJOB), in Itajuipe, Bahia. The area of the experiment is of approximately 3 ha, with a total of 3240 plants. Thirteen evaluations of vegetative brooms, five of cushion brooms and 15 of number of pods per plant were accomplished. Thirty pollinations were made for each selected plant to test for self-compatibility. The production, based on the number of pods per plant, and resistance to witches´ broom indicated CEPEC 94 x CCN 10, RB 39 x CCN 51 and CCN 10 x VB 1151 as superior progenies. All selections tested were self-compatible. The analyses of progenies and individual tree data, associated to visual field observations, allowed the selection of 17 plants which were included in a network of regional tests to determine the phenotypic stability.

Analysis of mediation and reciprocal effects of changes in weight, self-esteem, and exercise in a behavioral weight loss treatment

Improving the treatment of obesity remains a critical challenge. Several health behaviour change models, often based on a social-cognitive framework, have been used to design weight management interventions (Baranowski et al., 2003). However, most interventions have only produced modest weight reductions (Wadden et al., 2002) and socialcognitive variables have shown limited power to predict weight outcomes (Palmeira et al., 2007). Other predictors, and possibl alte nati e e planatory models, are needed to better understand the mechanisms by which weight loss and other obesity treatment-outcomes are brought about (Baranowski, 2006). Self-esteem is one of these possible mechanisms, because is commonly reported to change during the treatment, although these changes are not necessarily associated with weight loss (Blaine et al., 2007; Maciejewski et al., 2005). This possibility should be more evident if the program integrates regular exercise, as it promotes improvements in subjective well-being (Biddle & Mutrie, 2001), with possible influences on long-term behavioral adherence (e.g. diet, exercise). Following the reciprocal effects model tenets (Marsh & Craven, 2006), we expect that the influences between changes in weight, selfesteem and exercise to be reciprocal and might present one of the mechanisms by which obesity treatments can be improved.

Modeling the Diffusion of Energy Performance Contracting

The Portuguese Energy Policy considers the development of a commercially viable and competitive market for energy performance contracting (EPC) as a main mechanism to achieve the objectives of energy efficiency improvement. This paper proposes a study to investigate how to achieve widespread adoption of energy performance contracting by means of system dynamics modelling and simulation. To explore and gather insights on this question, a system dynamics model representing the system of the Portuguese EPC market at industry level will be created. The simulation of that model will provide a helpful basis for analysing and explaining the development of key variables, and accelerating learning on the managerial, organizational and political adaptation processes that foster the diffusion of EPC adoption. The first phase of this research project aims at identifying and analysing the key factors and critical cause-effect relations that drive the adoption of EPC. With this purpose, a qualitative content analysis on relevant documents was performed and a set of interviews was conducted. That data was analysed to capture the critical variables and its interrelation to formulate a preliminary representation of the system structure as stock and flow diagrams.

Self-interacting scalar field cosmologies: unified exact solutions and symmetries

We investigate a mechanism that generates exact solutions of scalar field cosmologies in a unified way. The procedure investigated here permits to recover almost all known solutions, and allows one to derive new solutions as well. In particular, we derive and discuss one novel solution defined in terms of the Lambert function. The solutions are organised in a classification which depends on the choice of a generating function which we have denoted by x(phi) that reflects the underlying thermodynamics of the model. We also analyse and discuss the existence of form-invariance dualities between solutions. A general way of defining the latter in an appropriate fashion for scalar fields is put forward.

Equilibrium self-assembly of colloids with distinct interaction sites: Thermodynamics, percolation, and cluster distribution functions

We calculate the equilibrium thermodynamic properties, percolation threshold, and cluster distribution functions for a model of associating colloids, which consists of hard spherical particles having on their surfaces three short-ranged attractive sites (sticky spots) of two different types, A and B. The thermodynamic properties are calculated using Wertheim's perturbation theory of associating fluids. This also allows us to find the onset of self-assembly, which can be quantified by the maxima of the specific heat at constant volume. The percolation threshold is derived, under the no-loop assumption, for the correlated bond model: In all cases it is two percolated phases that become identical at a critical point, when one exists. Finally, the cluster size distributions are calculated by mapping the model onto an effective model, characterized by a-state-dependent-functionality (f) over bar and unique bonding probability (p) over bar. The mapping is based on the asymptotic limit of the cluster distributions functions of the generic model and the effective parameters are defined through the requirement that the equilibrium cluster distributions of the true and effective models have the same number-averaged and weight-averaged sizes at all densities and temperatures. We also study the model numerically in the case where BB interactions are missing. In this limit, AB bonds either provide branching between A-chains (Y-junctions) if epsilon(AB)/epsilon(AA) is small, or drive the formation of a hyperbranched polymer if epsilon(AB)/epsilon(AA) is large. We find that the theoretical predictions describe quite accurately the numerical data, especially in the region where Y-junctions are present. There is fairly good agreement between theoretical and numerical results both for the thermodynamic (number of bonds and phase coexistence) and the connectivity properties of the model (cluster size distributions and percolation locus).

Validity of self-reported weight: a study of urban brazilian adults

In order to evaluate the validity of self-reported weight for use in obesity prevalence surveys, self-reported weight was compared to measured weight for 659 adults living in the Porto Alegre county, RS Brazil in 1986-87, both weights being obtained by a technician in the individual's home on the same visit. The mean difference between self-reported and measured weight was small (-0.06 +/- 3.16 kg; mean +/- standard deviation), and the correlation between reported and measured weight was high (r=0.97). Sixty-two percent of participants reported their weight with an error of < 2 kg, 87% with an error of < 4 kg, and 95% with an error of < 6 kg. Underweight individuals overestimated their weight, while obese individuals underestimated theirs (p<0.05). Men tended to overestimate their weight and women underestimate theirs, this difference between sexes being statistically significant (p=0.04). The overall prevalence of underweight (body mass index < 20) by reported weight was 11%, by measured weight 13%; the overall prevalence of obesity (body mass index > 30) by reported weight was 10%, by measured weight 11%. Thus, the validity of reported weight is acceptable for surveys of the prevalence of ponderosity in similar settings.