A cluster computing implementation of a Monte Carlo simulation of a particle growth mechanism

Clusters of computers can be used together to provide a powerful computing resource. Large Monte Carlo simulations, such as those used to model particle growth, are computationally intensive and take considerable time to execute on conventional workstations. By spreading the work of the simulation across a cluster of computers, the elapsed execution time can be greatly reduced. Thus a user has apparently the performance of a supercomputer by using the spare cycles on other workstations.

Advancements in the estimation of ice particle fall speeds using laboratory and field measurements

Accurate estimates for the fall speed of natural hydrometeors are vital if their evolution in clouds is to be understood quantitatively. In this study, laboratory measurements of the terminal velocity vt for a variety of ice particle models settling in viscous fluids, along with wind-tunnel and field measurements of ice particles settling in air, have been analyzed and compared to common methods of computing vt from the literature. It is observed that while these methods work well for a number of particle types, they fail for particles with open geometries, specifically those particles for which the area ratio Ar is small (Ar is defined as the area of the particle projected normal to the flow divided by the area of a circumscribing disc). In particular, the fall speeds of stellar and dendritic crystals, needles, open bullet rosettes, and low-density aggregates are all overestimated. These particle types are important in many cloud types: aggregates in particular often dominate snow precipitation at the ground and vertically pointing Doppler radar measurements. Based on the laboratory data, a simple modification to previous computational methods is proposed, based on the area ratio. This new method collapses the available drag data onto an approximately universal curve, and the resulting errors in the computed fall speeds relative to the tank data are less than 25% in all cases. Comparison with the (much more scattered) measurements of ice particles falling in air show strong support for this new method, with the area ratio bias apparently eliminated.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KMGAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KMGAP is based on the PRA model framework (P¨oschl-Rudich- Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modeled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system and the computational constraints, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmospheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270K is close to unity (Winkler et al., 2006). Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for efficient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Coupling an aerosol box model with one-dimensional flow: a tool for understanding observations of new particle formation events

Field observations of new particle formation and the subsequent particle growth are typically only possible at a fixed measurement location, and hence do not follow the temporal evolution of an air parcel in a Lagrangian sense. Standard analysis for determining formation and growth rates requires that the time-dependent formation rate and growth rate of the particles are spatially invariant; air parcel advection means that the observed temporal evolution of the particle size distribution at a fixed measurement location may not represent the true evolution if there are spatial variations in the formation and growth rates. Here we present a zero-dimensional aerosol box model coupled with one-dimensional atmospheric flow to describe the impact of advection on the evolution of simulated new particle formation events. Wind speed, particle formation rates and growth rates are input parameters that can vary as a function of time and location, using wind speed to connect location to time. The output simulates measurements at a fixed location; formation and growth rates of the particle mode can then be calculated from the simulated observations at a stationary point for different scenarios and be compared with the ‘true’ input parameters. Hence, we can investigate how spatial variations in the formation and growth rates of new particles would appear in observations of particle number size distributions at a fixed measurement site. We show that the particle size distribution and growth rate at a fixed location is dependent on the formation and growth parameters upwind, even if local conditions do not vary. We also show that different input parameters used may result in very similar simulated measurements. Erroneous interpretation of observations in terms of particle formation and growth rates, and the time span and areal extent of new particle formation, is possible if the spatial effects are not accounted for.

Pseudoclassical description of scalar particle in non-abelian background and path-integral representations

Path-integral representations for a scalar particle propagator in non-Abelian external backgrounds are derived. To this aim, we generalize the procedure proposed by Gitman and Schvartsman of path-integral construction to any representation of SU(N) given in terms of antisymmetric generators. And for arbitrary representations of SU(N), we present an alternative construction by means of fermionic coherent states. From the path-integral representations we derive pseudoclassical actions for a scalar particle placed in non-Abelian backgrounds. These actions are classically analyzed and then quantized to prove their consistency.

Self-adjoint extensions and spectral analysis in the Calogero problem

In this paper, we present a mathematically rigorous quantum-mechanical treatment of a one-dimensional motion of a particle in the Calogero potential alpha x(-2). Although the problem is quite old and well studied, we believe that our consideration based on a uniform approach to constructing a correct quantum-mechanical description for systems with singular potentials and/or boundaries, proposed in our previous works, adds some new points to its solution. To demonstrate that a consideration of the Calogero problem requires mathematical accuracy, we discuss some `paradoxes` inherent in the `naive` quantum-mechanical treatment. Using a self-adjoint extension method, we construct and study all possible self-adjoint operators (self-adjoint Hamiltonians) associated with a formal differential expression for the Calogero Hamiltonian. In particular, we discuss a spontaneous scale-symmetry breaking associated with self-adjoint extensions. A complete spectral analysis of all self-adjoint Hamiltonians is presented.

Self-adjoint extensions and spectral analysis in the generalized Kratzer problem

We present a mathematically rigorous quantum-mechanical treatment of a one-dimensional non-relativistic motion of a particle in the potential field V(x) = g(1)x(-1) + g(2)x(-2), x is an element of R(+) = [0, infinity). For g(2) > 0 and g(1) < 0, the potential is known as the Kratzer potential V(K)(x) and is usually used to describe molecular energy and structure, interactions between different molecules and interactions between non-bonded atoms. We construct all self-adjoint Schrodinger operators with the potential V(x) and represent rigorous solutions of the corresponding spectral problems. Solving the first part of the problem, we use a method of specifying self-adjoint extensions by (asymptotic) self-adjoint boundary conditions. Solving spectral problems, we follow Krein`s method of guiding functionals. This work is a continuation of our previous works devoted to the Coulomb, Calogero and Aharonov-Bohm potentials.

Tightened Lieb-Oxford Bound for Systems of Fixed Particle Number

The Lieb-Oxford bound is a constraint upon approximate exchange-correlation functionals. We explore a nonempirical tightening of that bound in both universal and electron number-dependent form. The test functional is PBE. Regarding both atomization energies (slightly worsened) and bond lengths (slightly improved), we find the PBE functional to be remarkably insensitive to the value of the Lieb-Oxford bound. This both rationalizes the use of the original Lieb-Oxford constant in PBE and suggests that enhancement factors more sensitive to sharpened constraints await discovery.

Effects of corn particle size and physical form of the diet on the gastrointestinal structures of broiler chickens

The objective of this experiment was to investigate the effects of different particle sizes, expressed as Geometric Mean Diameter (GMD) of corn (0.336mm, 0.585mm, 0.856mm and 1.12mm) of mash and pelleted broiler chicken diets on the weight of the gizzard, duodenum and jejunum+ileum; on the pH of the gizzard and small intestine and on the characteristics of the duodenal mucous layer (number and height of villi and crypt depth) in 42-day-old broilers. The physical form and the particle size of the diet had no significant effect on gizzard and intestine pH (p > 0.05). A greater gizzard weight was seen in the birds receiving pelleted diet and particle size of 0.336mm (p < 0.008). However, for the particle sizes of 0.856 and 1.12 mm, a greater weight was found in birds that received mash diet (p < 0.039 and p < 0.006, respectively). Also, gizzard weight was greater with increasing corn GMD independent of the physical form of the diet. In the mash diet, the increase in particle size promoted a quadratic response in the weight of duodenum and jejunum + ileum. The pelleted diet promoted a greater number of villi per transverse duodenum cut (p < 0.007) and greater crypt depth (p < 0.05). As the particle size increased, there was a linear increase of villus height and crypt depth in the duodenum, irrespective of the physical form of the diet.

Effects of the particle size and nucleation temperature on tellurite 20Li(2)O-80TeO(2) glass crystallization

20Li(2)O-80TeMO(2) glasses were heat annealed at different temperatures between T-g and T-x and studied by using XRD, FTIR spectroscopy and DSC techniques to understand the crystallization kinetics in this glass matrix. The infrared band structure of this glass is similar to what was observed in glassy TeO2. XRD results reveal the presence of three distinct crystalline gamma-TeO2, alpha-TeO2 and Li2Te2O5 phases during the crystallization process. This is a first report of gamma-TeO2 phase crystallization in this glass matrix. DSC results confinn the crystallization of three distinct structures in the glass. In summary, our results suggest a crystallization hierarchy on this glass matrix since the gamma-TeO2 and alpha-TeO2 phases crystallization occurs before the Li2Te2O5 phase crystallization. (c) 2006 Elsevier B.V. All rights reserved.

Effect of Airborne-Particle Abrasion and Mechanico-Thermal Cycling on the Flexural Strength of Glass Ceramic Fused to Gold or Cobalt-Chromium Alloy

Purpose: To evaluate the effect of airborne-particle abrasion and mechanico-thermal cycling on the flexural strength of a ceramic fused to cobalt-chromium alloy or gold alloy.Materials and Methods: Metallic bars (n = 120) were made (25 mm x 3 mm x 0.5 mm): 60 with gold alloy and 60 with Co-Cr. At the central area of the bars (8 mm x 3 mm), a layer of opaque ceramic and then two layers of glass ceramic (Vita VM13, Vita Zahnfabrick) were fired onto it (thickness: 1 mm). Ten specimens from each alloy group were randomly allocated to a surface treatment [(tungsten bur or air-particle abrasion (APA) with Al(2)O(3) at 10 mm or 20 mm away)] and mechanico-thermal cycling (no cycling or mechanically loaded 20,000 cycles; 10 N distilled water at 37 degrees C and then thermocycled 3000 cycles; 5 degrees C to 55 degrees C, dwell time 30 seconds) combination. Those specimens that did not undergo mechanico-thermal cyclingwere stored inwater (37 degrees C) for 24 hours. Bond strength was measured using a three-point bend test, according to ISO 9693. After the flexural strength test, failure types were noted. The data were analyzed using three factor-ANOVA and Tukey's test (alpha = 0.05).Results: There were no significant differences between the flexural bond strength of gold and Co-Cr groups (42.64 +/- 8.25 and 43.39 +/- 10.89 MPa, respectively). APA 10 and 20 mm away surface treatment (45.86 +/- 9.31 and 46.38 +/- 8.89 MPa, respectively) had similar mean flexural strength values, and both had significantly higher bond strength than tungsten bur treatment (36.81 +/- 7.60 MPa). Mechanico-thermal cycling decreased the mean flexural strength values significantly for all six alloy-surface treatment combinations tested when compared to the control groups. The failure type was adhesive in the metal/ceramic interface for specimens surface treated only with the tungsten bur, and mixed for specimens surface treated with APA 10 and 20 mm.Conclusions: Considering the levels adopted in this study, the alloy did not affect the bond strength; APA with Al(2)O(3) at 10 and 20 mm improved the flexural bond strength between ceramics and alloys used, and the mechanico-thermal cycling of metal-ceramic specimens resulted in a decrease of bond strength.

CRITICAL EXPONENTS and SCALING PROPERTIES FOR THE CHAOTIC DYNAMICS of A PARTICLE IN A TIME-DEPENDENT POTENTIAL BARRIER

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)