An official website of the United States government
In particle-based stochastic reaction–diffusion models, reaction rates and placement kernels are used to decide the probability per time a reaction can occur between reactant particles and to decide where product particles should be placed. When choosing kernels to use in reversible reactions, a key constraint is to ensure that detailed balance of spatial reaction fluxes holds at all points at equilibrium. In this work, we formulate a general partial-integral differential equation model that encompasses several of the commonly used contact reactivity (e.g., Smoluchowski-Collins-Kimball) and volume reactivity (e.g., Doi) particle models. From these equations, we derive a detailed balance condition for the reversible A + B ⇆ C reaction. In bounded domains with no-flux boundary conditions, when choosing unbinding kernels consistent with several commonly used binding kernels, we show that preserving detailed balance of spatial reaction fluxes at all points requires spatially varying unbinding rate functions near the domain boundary. Brownian dynamics simulation algorithms can realize such varying rates through ignoring domain boundaries during unbinding and rejecting unbinding events that result in product particles being placed outside the domain.
more »
« less
Steady states of active Brownian particles interacting with boundaries
Abstract An active Brownian particle is a minimal model for a self-propelled colloid in a dissipative environment. Experiments and simulations show that, in the presence of boundaries and obstacles, active Brownian particle systems approach nontrivial nonequilibrium steady states with intriguing phenomenology, such as accumulation at boundaries, ratchet effects, and long-range depletion interactions. Nevertheless, theoretical analysis of these phenomena has proven difficult. Here, we address this theoretical challenge in the context of non-interacting particles in two dimensions, basing our analysis on the steady-state Smoluchowski equation for the one-particle distribution function. Our primary result is an approximation strategy that connects asymptotic solutions of the Smoluchowski equation to boundary conditions. We test this approximation against the exact analytic solution in a 2D planar geometry, as well as numerical solutions in circular and elliptic geometries. We find good agreement so long as the boundary conditions do not vary too rapidly with respect to the persistence length of particle trajectories. Our results are relevant for characterizing long-range flows and depletion interactions in such systems. In particular, our framework shows how such behaviors are connected to the breaking of detailed balance at the boundaries.
more »
« less
- PAR ID:
- 10339576
- Date Published:
- Journal Name:
- Journal of Statistical Mechanics: Theory and Experiment
- Volume:
- 2022
- Issue:
- 1
- ISSN:
- 1742-5468
- Page Range / eLocation ID:
- 013208
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We prove the existence and uniqueness of positive analytical solutions with positive initial data to the mean field equation (the Dyson equation) of the Dyson Brownian motion through the complex Burgers equation with a force term on the upper half complex plane. These solutions converge to a steady state given by Wigner’s semicircle law. A unique global weak solution with nonnegative initial data to the Dyson equation is obtained, and some explicit solutions are given by Wigner’s semicircle laws. We also construct a bi-Hamiltonian structure for the system of real and imaginary components of the complex Burgers equation (coupled Burgers system). We establish a kinetic formulation for the coupled Burgers system and prove the existence and uniqueness of entropy solutions. The coupled Burgers system in Lagrangian variable naturally leads to two interacting particle systems, the Fermi–Pasta–Ulam–Tsingou model with nearest-neighbor interactions, and the Calogero–Moser model. These two particle systems yield the same Lagrangian dynamics in the continuum limit.more » « less
-
Many biological processes involve transport and organization of inclusions in thin fluid interfaces. A key aspect of these assemblies is the active dissipative stresses applied from the inclusions to the fluid interface, resulting in long-range active interfacial flows. We study the effect of these active flows on the self-organization of rod-like inclusions in the interface. Specifically, we consider a di- lute suspension of Brownian rods of length L, embedded in a thin fluid interface of 2D viscosity ηm and surrounded on both sides with 3D fluid domains of viscosity ηf . The momentum transfer from the interfacial flows to the surrounding fluids occurs over length l0 = ηm/ηf , known as Saffman- Delbru ̈ck length. We use zeroth, first and second moments of Smoluchowski equation to obtain the conservation equations for concentration, polar order and nematic order fields, and use linear stability analysis and continuum simulations to study the dynamic variations of these fields as a function of L/l0, the ratio of active to thermal stresses, and the dimensionless self-propulsion velocity of the embedded particles. We find that at sufficiently large activities, the suspensions of active extensile stress (pusher) with no directed motion undergo a finite wavelength nematic ordering, with the length of the ordered domains decreasing with increasing L/l0. The ordering transition is hindered with further increases in L/l0. In contrast, the suspensions with active contractile stress (puller) remain uniform with variations of activity. We notice that the self-propulsion velocity results in significant concentration fluctuations and changes in the size of the order domains that depend on L/l0. Our re- search highlights the role of hydrodynamic interactions in the self-organization of active inclusions on biological interfaces.more » « less
-
The phenomenon of the ‘coil-stretch’ (C-S) transition, wherein a long-chain polymer initially in a coiled state undergoes a sudden configuration change to become fully stretched under steady elongational flows, has been widely recognized. This transition can display intricate hysteresis behaviours under specific critical conditions, giving rise to unique rheological characteristics in dilute polymer solutions. Historically, microscopic stochastic models and Brownian dynamics simulations have shed light on the underlying mechanisms of the transition by uncovering bistable configurations of polymer chains. Following the initial work by Cerf (J. Chem. Phys., vol. 20, 1952, pp. 395–402), we introduce a continuum model in this study to investigate the C-S transition in a constant uniaxial elongational flow. Our approach involves approximating the unfolding process of the polymer chain as an axisymmetric deformation of an elastic particle. We make the assumption that the particle possesses uniform material properties, which can be represented by a nonlinear, strain-hardening constitutive equation to replicate the finite extensibility of the polymer chain. Subsequently, we analytically solve for the steady-state deformation using a polarization method. By employing this reduced model, we effectively capture the C-S transition and establish its specific correlations with material and geometric properties. The hysteresis phenomena can be comprehended through a force-balance analysis, which involves comparing the externally applied viscous forces with the intrinsic elastic responsive forces. We demonstrate that our model, while simple, unveils rich elastohydrodynamics of the C-S transition.more » « less
-
We construct a determining form for the 2D Rayleigh-Benard (RB) system in a strip with solid horizontal boundaries, in the cases of no-slip and stress-free boundary conditions. The determining form is an ODE in a Banach space of trajectories whose steady states comprise the long-time dynamics of the RB system. In fact, solutions on the global attractor of the RB system can be further identified through the zeros of a scalar equation to which the ODE reduces for each initial trajectory. The twist in this work is that the trajectories are for the velocity field only, which in turn determines the corresponding trajectories of the temperature.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1742-5468/ac42cf
