Towards an Understanding of the Influence of Sedimentation on Colloidal Aggregation by Peclet Number

The Peclet number is a useful index to estimate the importance of sedimentation as compared to the Brownian motion. However, how to choose the characteristic length scale for the Peclet number evaluation is rather critical because the diffusion length increases as the square root of the time whereas the drifting length is linearly related to time. Our Brownian dynamics simulation shows that the degree of sedimentation influence on the coagulation decreases when the dispersion volume fraction increases. Therefore using a fixed length, such as the diameter of particle, as the characteristic length scale for Peclet number evaluation is not a good choice when dealing with the influence of sedimentation on coagulation. The simulations demonstrated that environmental factors in the coagulation process, such as dispersion volume fraction and size distribution, should be taken into account for more reasonable evaluation of the sedimentation influence.

Computer simulation on the collision-sticking dynamics of two colloidal particles in an optical trap

Collisions of a particle pair induced by optical tweezers have been employed to study colloidal stability. In order to deepen insights regarding the collision-sticking dynamics of a particle pair in the optical trap that were observed in experimental approaches at the particle level, the authors carry out a Brownian dynamics simulation. In the simulation, various contributing factors, including the Derjaguin-Landau-Verwey-Overbeek interaction of particles, hydrodynamic interactions, optical trapping forces on the two particles, and the Brownian motion, were all taken into account. The simulation reproduces the tendencies of the accumulated sticking probability during the trapping duration for the trapped particle pair described in our previous study and provides an explanation for why the two entangled particles in the trap experience two different statuses. (c) 2007 American Institute of Physics.

Physical investigations of small particles : (I) aerosol particle charging and flux enhancement and (II) whispering gallery mode sensing

Part I

Particles are a key feature of planetary atmospheres. On Earth they represent the greatest source of uncertainty in the global energy budget. This uncertainty can be addressed by making more measurement, by improving the theoretical analysis of measurements, and by better modeling basic particle nucleation and initial particle growth within an atmosphere. This work will focus on the latter two methods of improvement.

Uncertainty in measurements is largely due to particle charging. Accurate descriptions of particle charging are challenging because one deals with particles in a gas as opposed to a vacuum, so different length scales come into play. Previous studies have considered the effects of transition between the continuum and kinetic regime and the effects of two and three body interactions within the kinetic regime. These studies, however, use questionable assumptions about the charging process which resulted in skewed observations, and bias in the proposed dynamics of aerosol particles. These assumptions affect both the ions and particles in the system. Ions are assumed to be point monopoles that have a single characteristic speed rather than follow a distribution. Particles are assumed to be perfect conductors that have up to five elementary charges on them. The effects of three body interaction, ion-molecule-particle, are also overestimated. By revising this theory so that the basic physical attributes of both ions and particles and their interactions are better represented, we are able to make more accurate predictions of particle charging in both the kinetic and continuum regimes.

The same revised theory that was used above to model ion charging can also be applied to the flux of neutral vapor phase molecules to a particle or initial cluster. Using these results we can model the vapor flux to a neutral or charged particle due to diffusion and electromagnetic interactions. In many classical theories currently applied to these models, the finite size of the molecule and the electromagnetic interaction between the molecule and particle, especially for the neutral particle case, are completely ignored, or, as is often the case for a permanent dipole vapor species, strongly underestimated. Comparing our model to these classical models we determine an “enhancement factor” to characterize how important the addition of these physical parameters and processes is to the understanding of particle nucleation and growth.

Part II

Whispering gallery mode (WGM) optical biosensors are capable of extraordinarily sensitive specific and non-specific detection of species suspended in a gas or fluid. Recent experimental results suggest that these devices may attain single-molecule sensitivity to protein solutions in the form of stepwise shifts in their resonance wavelength, \lambda_{R}, but present sensor models predict much smaller steps than were reported. This study examines the physical interaction between a WGM sensor and a molecule adsorbed to its surface, exploring assumptions made in previous efforts to model WGM sensor behavior, and describing computational schemes that model the experiments for which single protein sensitivity was reported. The resulting model is used to simulate sensor performance, within constraints imposed by the limited material property data. On this basis, we conclude that nonlinear optical effects would be needed to attain the reported sensitivity, and that, in the experiments for which extreme sensitivity was reported, a bound protein experiences optical energy fluxes too high for such effects to be ignored.

Suspension mechanics: I. Inertial and non-Newtonian migration of neutrally buoyant rigid spheres in two-dimensional unidirectional flows; II. The effect of viscoelasticity on the creeping motion of a train of neutrally buoyant Newtonian drops through a circular tube

The lateral migration of neutrally buoyant rigid spheres in two-dimensional unidirectional flows was studied theoretically. The cases of both inertia-induced migration in a Newtonian fluid and normal stress-induced migration in a second-order fluid were considered. Analytical results for the lateral velocities were obtained, and the equilibrium positions and trajectories of the spheres compared favorably with the experimental data available in the literature. The effective viscosity was obtained for a dilute suspension of spheres which were simultaneously undergoing inertia-induced migration and translational Brownian motion in a plane Poiseuille flow. The migration of spheres suspended in a second-order fluid inside a screw extruder was also considered.

The creeping motion of neutrally buoyant concentrically located Newtonian drops through a circular tube was studied experimentally for drops which have an undeformed radius comparable to that of the tube. Both a Newtonian and a viscoelastic suspending fluid were used in order to determine the influence of viscoelasticity. The extra pressure drop due to the presence of the suspended drops, the shape and velocity of the drops, and the streamlines of the flow were obtained for various viscosity ratios, total flow rates, and drop sizes. The results were compared with existing theoretical and experimental data.

Extraction of polydispersity information from photon correlation data

Photon correlation spectroscopy (PCS) is a light-scattering technique for particle size diagnosis. It has been used mainly in the investigation of hydrosol particles since it is based on the measurement of the correlation function of the light scattered from the Brownian motion of suspended particles. Recently this technique also proved useful for studying soot particles in flames and similar aerosol systems. In the case of a polydispersed system the problem of recovering the particle size distribution can be reduced to the problem of inverting the Laplace transform. In this paper we review several methods introduced by the authors for the solution of this problem. We present some numerical results and we discuss the resolution limits characterizing the reconstruction of the size distributions. © 1989.

Localized excitations of charged dust grains in dusty plasma lattices

The nonlinear aspects of charged dust grain motion in a one-dimensional dusty plasma (DP) monolayer are discussed. Both horizontal (longitudinal, acoustic mode) and vertical (transverse, optic mode) displacements are considered, and various types of localized excitations are reviewed, in a continuum approximation. Dust crystals are shown to support nonlinear kink-shaped supersonic longitudinal solitary excitations, as well as modulated envelope (either longitudinal or transverse) localized modes. The possibility for Discrete Breather (DB-) type excitations (Intrinsic Localized Modes, ILMs) to occur is investigated, from first principles. These highly localized excitations owe their existence to lattice discreteness, in combination with the interaction and/or
substrate (sheath) potential nonlinearity. This possibility may open new directions in DP- related research. The relation to previous results on atomic chains as well as to experimental results on strongly-coupled dust layers in gas discharge plasmas is discussed.

LINEAR ACCELERATION EMISSION IN PULSAR MAGNETOSPHERES

Linear acceleration emission occurs when a charged particle is accelerated parallel to its velocity. We evaluate the spectral and angular distribution of this radiation for several special cases, including constant acceleration (hyperbolic motion) of finite duration. Based on these results, we find the following general properties of the emission from an electron in a linear accelerator that can be characterized by an electric field E acting over a distance L: (1) the spectrum extends to a cutoff frequency (h) over bar omega(c)/mc(2) approximate to L(E/E(Schw))(2)/(lambda) over bar (C), where E(Schw) = 1.3 x 10(18) V m(-1) is the Schwinger critical field and (lambda) over bar (C) = (h) over bar /mc = 3.86 x 10(-13) m is the Compton wavelength of the electron, (2) the total energy emitted by a particle traversing the accelerator is 4/3 alpha(f)(h) over bar omega(c) in accordance with the standard Larmor formula where alpha(f) is the fine-structure constant, and (3) the low frequency spectrum is flat for hyperbolic trajectories, but in general depends on the details of the accelerator. We also show that linear acceleration emission complements curvature radiation in the strongly magnetized pair formation regions in pulsar magnetospheres. It dominates when the length L of the accelerator is less than the formation length rho/gamma of curvature photons, where rho is the radius of curvature of the magnetic field lines and gamma the Lorentz factor of the emitting particle. In standard static models of pair creating regions linear acceleration emission is negligible, but it is important in more realistic dynamical models in which the accelerating field fluctuates on a short length scale.

Irradiation of oxygen and water ices at by 4 keV singly and doubly charged ions; sputtering and molecular synthesis

The icy surfaces of dust grains in the Interstellar Medium and those of comets, satellites and Kuiper Belt Objects are continuously exposed tophoton and charged particle irradiation. These energetic particles maysputter and induce chemical changes in the ices and the underlyingsurfaces.In the present work 258 nm thick O2 and H2O ices were deposited at 10 K with the thickness measured by a laser interferometer method. Asimple model fit to the reflected laser intensity as measured by aphotodiode detector enabled the refractive index of the ices to bedetermined. The ices were then irradiated with various singly and doublycharged ions such as He+, 13C+, N+, O+ , Ar+, 13C2+, N2+ and O2+ at 4keV. The decrease in ice thickness as a function of ion dose wasmonitored by a laser interferometer and the model used to determine thesputtering yield as shown in Figure 1.In the case of O2 ice thesputtering yields increased with increasing ion mass in good agreementwith a model calculation [Fama, J, Shi, R.A Baragiola, Surface Sci.,602, 156 (2007)]. In the case of O2 ice, O2+ has a significant lowersputtering yield when compared to O+. The sputtering yields for O2 icewere found to be at least 9 times larger compared to those for H2O ice.For H2O ice the sputter yields for C, N and O ions were found todecrease with increasing mass. Doubly charged C, N and O ions which werefound to have the same sputtering yield as the singly charged ionswithin the experimental errors. A preliminary TPD study was carried outusing a QMS to detect the desorbed species from water ice afterirradiation by 6 × 10^15 ions of 13C+ and 13C2+. The formation of13CO and 13CO2 was observed with the yield of 13CO almost of a factor of100 larger than of 13CO2. This is in contrast to our earlier work whereonly CO¬2 was observed.

SIMLA: Simulating particle dynamics in intense laser and other electromagnetic fields via classical and quantum electrodynamics

We present the Fortran program SIMLA, which is designed for the study of charged particle dynamics in laser and other background fields. The dynamics can be determined classically via the Lorentz force and Landau–Lifshitz equations or, alternatively, via the simulation of photon emission events determined by strong-field quantum-electrodynamics amplitudes and implemented using Monte-Carlo routines. Multiple background fields can be included in the simulation and, where applicable, the propagation direction, field type (plane wave, focussed paraxial, constant crossed, or constant magnetic), and time envelope of each can be independently specified.

Quantifying diffusion in biofilms : from model hydrogels to living biofilms

Les biofilms sont des communautés de microorganismes incorporés dans une matrice exo-polymérique complexe. Ils sont reconnus pour jouer un rôle important comme barrière de diffusion dans les systèmes environnementaux et la santé humaine, donnant lieu à une résistance accrue aux antibiotiques et aux désinfectants. Comme le transfert de masse dans un biofilm est principalement dû à la diffusion moléculaire, il est primordial de comprendre les principaux paramètres influençant les flux de diffusion. Dans ce travail, nous avons étudié un biofilm de Pseudomonas fluorescens et deux hydrogels modèles (agarose et alginate) pour lesquels l’autodiffusion (mouvement Brownien) et les coefficients de diffusion mutuels ont été quantifiés. La spectroscopie par corrélation de fluorescence a été utilisée pour mesurer les coefficients d'autodiffusion dans une volume confocal de ca. 1 m3 dans les gels ou les biofilms, tandis que les mesures de diffusion mutuelle ont été faites par cellule de diffusion. En outre, la voltamétrie sur microélectrode a été utilisée pour évaluer le potentiel de Donnan des gels afin de déterminer son impact sur la diffusion. Pour l'hydrogel d'agarose, les observations combinées d'une diminution du coefficient d’autodiffusion et de l’augmentation de la diffusion mutuelle pour une force ionique décroissante ont été attribuées au potentiel de Donnan du gel. Des mesures de l'effet Donnan (différence de -30 mV entre des forces ioniques de 10-4 et 10-1 M) et l'accumulation correspondante d’ions dans l'hydrogel (augmentation d’un facteur de 13 par rapport à la solution) ont indiqué que les interactions électrostatiques peuvent fortement influencer le flux de diffusion de cations, même dans un hydrogel faiblement chargé tel que l'agarose. Curieusement, pour un gel plus chargé comme l'alginate de calcium, la variation de la force ionique et du pH n'a donné lieu qu'à de légères variations de la diffusion de sondes chargées dans l'hydrogel. Ces résultats suggèrent qu’en influençant la diffusion du soluté, l'effet direct des cations sur la structure du gel (compression et/ou gonflement induits) était beaucoup plus efficace que l'effet Donnan. De même, pour un biofilm bactérien, les coefficients d'autodiffusion étaient pratiquement constants sur toute une gamme de force ionique (10-4-10-1 M), aussi bien pour des petits solutés chargés négativement ou positivement (le rapport du coefficient d’autodiffusion dans biofilm sur celui dans la solution, Db/Dw ≈ 85 %) que pour des nanoparticules (Db/Dw≈ 50 %), suggérant que l'effet d'obstruction des biofilms l’emporte sur l'effet de charge. Les résultats de cette étude ont montré que parmi les divers facteurs majeurs qui affectent la diffusion dans un biofilm environnemental oligotrophe (exclusion stérique, interactions électrostatiques et hydrophobes), les effets d'obstruction semblent être les plus importants lorsque l'on tente de comprendre la diffusion du soluté. Alors que les effets de charge ne semblaient pas être importants pour l'autodiffusion de substrats chargés dans l'hydrogel d'alginate ou dans le biofilm bactérien, ils ont joué un rôle clé dans la compréhension de la diffusion à travers l’agarose. L’ensemble de ces résultats devraient être très utiles pour l'évaluation de la biodisponibilité des contaminants traces et des nanoparticules dans l'environnement.

Studies on the dynamics and rheology of dilute suspensions of Brownian particles in simple shear flow under the action of an external force

We present a novel approach to computing the orientation moments and rheological properties of a dilute suspension of spheroids in a simple shear flow at arbitrary Peclct number based on a generalised Langevin equation method. This method differs from the diffusion equation method which is commonly used to model similar systems in that the actual equations of motion for the orientations of the individual particles are used in the computations, instead of a solution of the diffusion equation of the system. It also differs from the method of 'Brownian dynamics simulations' in that the equations used for the simulations are deterministic differential equations even in the presence of noise, and not stochastic differential equations as in Brownian dynamics simulations. One advantage of the present approach over the Fokker-Planck equation formalism is that it employs a common strategy that can be applied across a wide range of shear and diffusion parameters. Also, since deterministic differential equations are easier to simulate than stochastic differential equations, the Langevin equation method presented in this work is more efficient and less computationally intensive than Brownian dynamics simulations.We derive the Langevin equations governing the orientations of the particles in the suspension and evolve a procedure for obtaining the equation of motion for any orientation moment. A computational technique is described for simulating the orientation moments dynamically from a set of time-averaged Langevin equations, which can be used to obtain the moments when the governing equations are harder to solve analytically. The results obtained using this method are in good agreement with those available in the literature.The above computational method is also used to investigate the effect of rotational Brownian motion on the rheology of the suspension under the action of an external force field. The force field is assumed to be either constant or periodic. In the case of con- I stant external fields earlier results in the literature are reproduced, while for the case of periodic forcing certain parametric regimes corresponding to weak Brownian diffusion are identified where the rheological parameters evolve chaotically and settle onto a low dimensional attractor. The response of the system to variations in the magnitude and orientation of the force field and strength of diffusion is also analyzed through numerical experiments. It is also demonstrated that the aperiodic behaviour exhibited by the system could not have been picked up by the diffusion equation approach as presently used in the literature.The main contributions of this work include the preparation of the basic framework for applying the Langevin method to standard flow problems, quantification of rotary Brownian effects by using the new method, the paired-moment scheme for computing the moments and its use in solving an otherwise intractable problem especially in the limit of small Brownian motion where the problem becomes singular, and a demonstration of how systems governed by a Fokker-Planck equation can be explored for possible chaotic behaviour.

Law of large numbers for super-brownian motions with a single point source

We investigate the super-Brownian motion with a single point source in dimensions 2 and 3 as constructed by Fleischmann and Mueller in 2004. Using analytic facts we derive the long time behavior of the mean in dimension 2 and 3 thereby complementing previous work of Fleischmann, Mueller and Vogt. Using spectral theory and martingale arguments we prove a version of the strong law of large numbers for the two dimensional superprocess with a single point source and finite variance.

Restoring the azimuthal symmetry of lateral distributions of charged particles in the range of the KASCADE-Grande experiment

The reconstruction of Extensive Air Showers (EAS) observed by particle detectors at the ground is based on the characteristics of observables like the lateral particle density and the arrival times. The lateral densities, inferred for different EAS components from detector data, are usually parameterised by applying various lateral distribution functions (LDFs). The LDFs are used in turn for evaluating quantities like the total number of particles or the density at particular radial distances. Typical expressions for LDFs anticipate azimuthal symmetry of the density around the shower axis. The deviations of the lateral particle density from this assumption arising from various reasons are smoothed out in the case of compact arrays like KASCADE, but not in the case of arrays like Grande, which only sample a smaller part of the azimuthal variation. KASCADE-Grande, an extension of the former KASCADE experiment, is a multi-component Extensive Air Shower (EAS) experiment located at the Karlsruhe Institute of Technology (Campus North), Germany. The lateral distributions of charged particles are deduced from the basic information provided by the Grande scintillators - the energy deposits - first in the observation plane, then in the intrinsic shower plane. In all steps azimuthal dependences should be taken into account. As the energy deposit in the scintillators is dependent on the angles of incidence of the particles, azimuthal dependences are already involved in the first step: the conversion from the energy deposits to the charged particle density. This is done by using the Lateral Energy Correction Function (LECF) that evaluates the mean energy deposited by a charged particle taking into account the contribution of other particles (e.g. photons) to the energy deposit. By using a very fast procedure for the evaluation of the energy deposited by various particles we prepared realistic LECFs depending on the angle of incidence of the shower and on the radial and azimuthal coordinates of the location of the detector. Mapping the lateral density from the observation plane onto the intrinsic shower plane does not remove the azimuthal dependences arising from geometric and attenuation effects, in particular for inclined showers. Realistic procedures for applying correction factors are developed. Specific examples of the bias due to neglecting the azimuthal asymmetries in the conversion from the energy deposit in the Grande detectors to the lateral density of charged particles in the intrinsic shower plane are given. (C) 2011 Elsevier B.V. All rights reserved.

Transverse-Momentum and Pseudorapidity Distributions of Charged Hadrons in pp Collisions at root s=7 TeV

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