skip to main content


Title: Evolution of coagulation-fragmentation stochastic processes using accurate chemical master equation approach
Coagulation and fragmentation (CF) is a fundamental process in which smaller particles attach to each other to form larger clusters while existing clusters break up into smaller particles. It is a ubiquitous process that plays important roles in many physical and biological phenomena. CF is typically a stochastic process that often occurs in confined spaces with a limited number of available particles . Here, we study the CF process formulated with the discrete Chemical Master Equation (dCME). Using the newly developed Accurate Chemical Master Equation (ACME) method, we examine calculate the time-dependent behavior of the CF system. We investigate the effects of a number of important factors that influence the overall behavior of the system, including the dimensionality, the ratio of attachment to detachment rates among clusters, and the initial conditions. By comparing CF in one and three dimensions, we conclude that systems in three dimensions are more likely to form large clusters. We also demonstrate how the ratio of the attachment to detachment rates affects the dynamics and the steady-state of the system. Finally, we demonstrate the relationship between the formation of large clusters and the initial condition.  more » « less
Award ID(s):
1814364
PAR ID:
10092709
Author(s) / Creator(s):
Date Published:
Journal Name:
Communications in information and systems
Volume:
19
Issue:
1
ISSN:
1526-7555
Page Range / eLocation ID:
37-55
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We developed a physical model to fundamentally understand the conductive filament (CF) formation and growth behavior in the switching layer during electroforming process in the metal-oxide-based resistive random-access memories (RRAM). The effects of the electrode and oxide layer properties on the CF morphology evolution, current-voltage characteristic, local temperature, and electrical potential distribution have been systematically explored. It is found that choosing active electrodes with lower oxygen vacancy formation energy and oxides with small Lorenz number (ratio of thermal and electrical conductivity) enables CF formation at a smaller electroforming voltage and creates a CF with more homogeneous morphology. This work advances our understanding of the kinetic behaviors of the CF formation and growth during the electroforming process and could potentially guide the oxide and electrode materials selection to realize a more stable and functional RRAM.

     
    more » « less
  2. Using the example of metal clusters, an experimental setup and procedure is presented, which allows for the generation of size and charge-state selected polyanions from monoanions in a molecular beam. As a characteristic feature of this modular setup, the further charging process via sequential electron attachment within a three-state digital trap takes place after mass-selection. In contrast to other approaches, the rf-based concept permits access to heavy particles. The procedure is highly flexible with respect to the preparation process and potentially suitable for a wide variety of anionic species. By adjusting the storage conditions, i.e., the radio frequency, to the change in the mass-to-charge ratio, we succeeded in producing clusters in highly negative charge states, i.e., [Formula: see text]. The capabilities of the setup are demonstrated by experiments extracting electronic and optical properties of polyanionic metal clusters by analyzing the corresponding photoelectron spectra.

     
    more » « less
  3. Anisotropic potentials in Janus colloids provide additional freedom to control particle aggregation into structures of different sizes and morphologies. In this work, we perform Brownian dynamics simulations of a dilute suspension of magnetic spherical Janus colloids with their magnetic dipole moments shifted radially towards the surface of the particle in order to gain valuable microstructural insight. Properties such as the mean cluster size, orientational ordering, and nucleation and growth are examined dynamically. Differences in the structure of clusters and in the aggregation process are observed depending on the dipolar shift ( s )—the ratio between the displacement of the dipole and the particle radius—and the dipolar coupling constant ( λ )—the ratio between the magnetic dipole–dipole and Brownian forces. Using these two dimensionless quantities, a structural “phase” diagram is constructed. Each phase corresponds to unique nucleation and growth behavior and orientational ordering of dipoles inside clusters. At small λ , the particles aggregate and disaggregate resulting in short-lived clusters at small s , while at high s the particles aggregate in permanent triplets (long-lived clusters). At high λ , the critical nuclei formed during the nucleation process are triplets and quadruplets with unique orientational ordering. These small clusters then serve as building blocks to form larger structures, such as single-chain, loop-like, island-like, worm-like, and antiparallel-double-chain clusters. This study shows that dipolar shifts in colloids can serve as a control parameter in applications where unique size, morphology, and aggregation kinetics of clusters are required. 
    more » « less
  4. The role of negative hydroxyl ions in liquid-phase plasma discharge formation is investigated using an inhouse modeling framework. Two tunneling sources for electrons are considered—tunneling ionization of water molecules and tunneling detachment of negative hydroxyl ions together with additional reaction steps. The simulations are conducted for a needle-like powered electrode with two different nanosecond rise time voltage profiles—a linear and an exponential rise. Both the profiles have a maximum voltage of 15 kV. The predictions show that the electron detachment, which has a much lower threshold energy requirement, provides a stream of electrons at low applied voltage during the initial rise time. The electrical forces from the electron detachment process generate stronger compression but a weaker expansion regime in the liquid resulting in ∼40% increase in the density and only ∼1% decrease. The electron detachment tunneling process is found to be not limited by the electric field, but rather by the availability of negative hydroxyl ions in the system and ceases when these ions are depleted. The tunnel ionization of water molecules forms the electron wave at a higher applied voltage, but the resulting peak electron number density is typically six orders of magnitude larger than the detachment tunneling. The higher electron number density allows the recycling of depleted negative hydroxyl ions in the system and can reestablish tunneling detachment. In addition, the system experiences a larger variation in density; specifically, a decrease in density due to tunnel ionization. The prediction also shows that irrespective of the initial electron sources (i.e. tunnel ionization or tunnel detachment) the reduced electric field is not sufficient enough to allow electron impact ionization to be active and make a significant contribution. Path flux analysis is conducted to determine the kinetics responsible for the recycling of the negative hydroxyl ions. 
    more » « less
  5. We develop a physical and computational model for performing fully coupled, grain-resolved direct numerical simulations of cohesive sediment, based on the immersed boundary method. The model distributes the cohesive forces over a thin shell surrounding each particle, thereby allowing for the spatial and temporal resolution of the cohesive forces during particle–particle interactions. The influence of the cohesive forces is captured by a single dimensionless parameter in the form of a cohesion number, which represents the ratio of cohesive and gravitational forces acting on a particle. We test and validate the cohesive force model for binary particle interactions in the drafting–kissing–tumbling (DKT) configuration. Cohesive sediment grains can remain attached to each other during the tumbling phase following the initial collision, thereby giving rise to the formation of flocs. The DKT simulations demonstrate that cohesive particle pairs settle in a preferred orientation, with particles of very different sizes preferentially aligning themselves in the vertical direction, so that the smaller particle is drafted in the wake of the larger one. This preferred orientation of cohesive particle pairs is found to remain influential for systems of higher complexity. To this end, we perform large simulations of 1261 polydisperse settling particles starting from rest. These simulations reproduce several earlier experimental observations by other authors, such as the accelerated settling of sand and silt particles due to particle bonding, the stratification of cohesive sediment deposits, and the consolidation process of the deposit. They identify three characteristic phases of the polydisperse settling process, viz. (i) initial stir-up phase with limited flocculation, (ii) enhanced settling phase characterized by increased flocculation, and (iii) consolidation phase. The simulations demonstrate that cohesive forces accelerate the overall settling process primarily because smaller grains attach to larger ones and settle in their wakes. For the present cohesive number values, we observe that settling can be accelerated by up to 29 %. We propose physically based parametrization of classical hindered settling functions introduced by earlier authors, in order to account for cohesive forces. An investigation of the energy budget shows that, even though the work of the collision forces is much smaller than that of the hydrodynamic drag forces, it can substantially modify the relevant energy conversion processes. 
    more » « less