Particle shape strongly influences the diffusion charging of aerosol particles exposed to bipolar/unipolar ions and accurate modeling is needed to predict the charge distribution of non-spherical particles. A prior particle-ion collision kernel β_i model including Coulombic and image potential interactions for spherical particles is generalized for arbitrary shapes following a scaling approach that uses a continuum and free molecular particle length scale and Langevin dynamics simulations of non-spherical particle-ion collisions for attractive Coulomb-image potential interactions. This extended β_i model for collisions between unlike charged particle-ion (bipolar charging) and like charged particle-ion (unipolar charging) is validated by comparing against published experimental data of bipolar charge distributions for diverse shapes. Comparison to the bipolar charging data for spherical particles shows good agreement in air, argon, and nitrogen, while also demonstrating high accuracy in predicting charge states up to ±6. Comparisons to the data for fractal aggregates reveal that the LD-based β_i model predicts within overall ±30% without any systematic bias. The mean charge on linear chain aggregates and charge fractions on cylindrical particles is found to be in good agreement with the measurements (~±20% overall). The comparison with experimental results supports the use of LD-based diffusion charging models to predict the bipolarmore »
Modeling nanoparticle charge distribution in the afterglow of non-thermal plasmas and comparison with measurements
Particle charging in the afterglows of non-thermal plasmas typically take place in a non-neutral space charge environment. We model the same by incorporating particle-ion collision rate constant models, developed in prior work by analyzing particle-ion trajectories calculated using Langevin Dynamics simulations, into species transport equations for ions, electrons and charged particles in the afterglow. A scaling analysis of particle charging and additional Langevin Dynamics calculations of the particle-ion collision rate constant are presented to extend the range of applicability to ion electrostatic to thermal energy ratios of 300 and diffusive Knudsen number (that scales inversely with gas pressure) up to 2000. The developed collision rate constant models are first validated by comparing predictions of particle charge against measured values in a stationary, non-thermal DC plasma from past PK-4 campaigns published in Phys. Rev. Lett. 93(8): 085001 and Phys. Rev. E 72(1): 016406). The comparisons reveal excellent agreement within ±35% for particles of radius 0.6,1.0,1.3 μm in the gas pressure range of ~20-150 Pa. The experiments to probe particle charge distributions by Sharma et al. (J. Physics D: Appl. Phys. 53(24): 245204) are modeled using the validated particle-ion collision rate constant models and the calculated charge fractions are compared with measurements. more »
- Award ID(s):
- 1903432
- Publication Date:
- NSF-PAR ID:
- 10228833
- Journal Name:
- Journal of physics
- Volume:
- 54
- Issue:
- 27
- Page Range or eLocation-ID:
- 275205
- ISSN:
- 0022-3727
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Based on the prior work of Chahl and Gopalakrishnan (2019) to infer particle-ion collision time distributions using a Langevin Dynamics (LD) approach, we develop a model for the non-dimensional diffusion charging collision kernel β_i or H that is applicable for 0≤Ψ_E≤60,0≤Ψ_I/Ψ_E ≤1,Kn_D≤2000 (defined in the main text). The developed model for β_i for attractive Coulomb and image potential interactions, along with the model for β_i for repulsive Coulomb and image potential interactions from Gopalakrishnan et al. (2013b), is tested against published diffusion charging experimental data. Current state of the art charging models, Fuchs (1963) and Wiedensohler (1988) regression for bipolar charging, are also evaluated and discussed. Comparisons reveal that the LD-based model accurately describes unipolar fractions for 10 – 100 nm particles measured in air (Adachi et al., 1985), nitrogen and argon but not in helium (Adachi et al., 1987). Fuchs model and the LD-based model yield similar predictions in the experimental conditions considered, except in helium. In the case of bipolar charging, the LD-based model captures the experimental trends quantitatively (within ±20%) across the entire size range of 4 – 40 nm producing superior agreement than Wiedensohler’s regression. The latter systematically underpredicts charge fraction below ~20 nm in airmore »
-
Radiative double-electron capture (RDEC), in which two-electron capture is accompanied by simultaneousemission of a single photon, was investigated for fully stripped and one-electron projectiles colliding withgaseous and thin-foil targets. RDEC can be considered the inverse of double photoionization by a single photon.For the gaseous targets, measurements were done for 2.11 MeV/uF9+and F8+ions interacting with N2and Ne,while for the thin-foil target the measurements were done for 2.11 MeV/uF9+and F8+and 2.19 MeV/uO8+andO7+ions striking thin C targets. Reports on this work were already published separately in shorter accounts by LaMantiaet al.[Phys. Rev. Lett.124, 133401 (2020)for the gas targets andPhys.Rev.A102, 060801(R) (2020)forthe thin-foil targets]. The gas targets were studied under single-collision conditions, while the foil targets sufferedunavoidable multiple collisions. The measurements were carried out by detecting x-ray emission from the targetat 90◦to the beam direction in coincidence with outgoing ions undergoing double, single, and, in the caseof the foil targets, no charge change inside the target. Striking differences between the gaseous and foil targetswere found from these measurements, with RDEC for the gaseous targets occurring only in coincidence with q-2outgoing projectiles as expected, while RDEC for the foil targets was seen in each of the outgoing q-2, q-1, and nocharge-change states. The no charge-changemore »
-
Abstract The use of silicon nanoparticles for lithium-ion batteries requires a precise control over both their average size and their size distribution. Particles larger than the generally accepted critical size of 150 nm fail during lithiation because of excessive swelling, while very small particles (<10 nm) inevitably lead to a poor first cycle coulombic efficiency because of their excessive specific surface area. Both mechanisms induce irreversible capacity losses and are detrimental to the anode functionality. In this manuscript we describe a novel approach for enhanced growth of nanoparticles to ∼20 nm using low-temperature flow-through plasma reactors via pulsing. Pulsing of the RF power leads to a significant increase in the average particle size, all while maintaining the particles well below the critical size for stable operation in a lithium-ion battery anode. A zero-dimensional aerosol plasma model is developed to provide insights into the dynamics of particle agglomeration and growth in the pulsed plasma reactor. The accelerated growth correlates with the shape of the particle size distribution in the afterglow, which is in turn controlled by parameters such as metastable density, gas and electron temperature. The accelerated agglomeration in each afterglow phase is followed by rapid sintering of the agglomerates intomore »
-
The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes:more »