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
- NSF-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
-
-
more » « less
Abstract Active particle systems can vary greatly from one-component systems of spheres to mixtures of particle shapes at different composition ratios. We investigate computationally the combined effect of anisotropy and stoichiometry on the collective behavior of two-dimensional active colloidal mixtures of polygons. We uncover three emergent phenomena not yet reported in active Brownian particle systems. First, we find that mixtures containing hexagons exhibit micro-phase separation with large grains of hexagonal symmetry. We quantify a measurable, implicit ‘steric attraction’ between the active particles as a result of shape anisotropy and activity. This calculation provides further evidence that implicit interactions in active systems, even without explicit attraction, can lead to an effective preferential attraction between particles. Next, we report stable fluid clusters in mixtures containing one triangle or square component. We attribute the fluidization of the dense cluster to the interplay of cluster destabilizing particles, which introduce grain boundaries and slip planes into the system, causing solid-like clusters to break up into fluid clusters. Third, we show that fluid clusters can coexist with solid clusters within a sparse gas of particles in a steady state of three coexisting phases. Our results highlight the potential for a wide variety of behavior to be accessible to active matter systems and establish a route to control active colloidal systems through simple parameter designs.
-
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
-
We investigate the effect of constant-vorticity background shear on the properties of wavetrains in deep water. Using the methodology of Fokas ( A Unified Approach to Boundary Value Problems , 2008, SIAM), we derive a higher-order nonlinear Schrödinger equation in the presence of shear and surface tension. We show that the presence of shear induces a strong coupling between the carrier wave and the mean-surface displacement. The effects of the background shear on the modulational instability of plane waves is also studied, where it is shown that shear can suppress instability, although not for all carrier wavelengths in the presence of surface tension. These results expand upon the findings of Thomas et al. ( Phys. Fluids , vol. 24 (12), 2012, 127102). Using a modification of the generalized Lagrangian mean theory in Andrews & McIntyre ( J. Fluid Mech. , vol. 89, 1978, pp. 609–646) and approximate formulas for the velocity field in the fluid column, explicit, asymptotic approximations for the Lagrangian and Stokes drift velocities are obtained for plane-wave and Jacobi elliptic function solutions of the nonlinear Schrödinger equation. Numerical approximations to particle trajectories for these solutions are found and the Lagrangian and Stokes drift velocities corresponding to these numerical solutions corroborate the theoretical results. We show that background currents have significant effects on the mean transport properties of waves. In particular, certain combinations of background shear and carrier wave frequency lead to the disappearance of mean-surface mass transport. These results provide a possible explanation for the measurements reported in Smith ( J. Phys. Oceanogr. , vol. 36, 2006, pp. 1381–1402). Our results also provide further evidence of the viability of the modification of the Stokes drift velocity beyond the standard monochromatic approximation, such as recently proposed in Breivik et al. ( J. Phys. Oceanogr. , vol. 44, 2014, pp. 2433–2445) in order to obtain a closer match to a range of complex ocean wave spectra.more » « less
-
Stimulated by the effect of the nearest neighbor interactions in vehicular traffic and motor proteins, we study a 1D driven lattice gas model, in which the nearest neighbor particle interactions are taken in accordance with the thermodynamic concepts. The non-equilibrium steady-state properties of the system are analyzed under both open and periodic boundary conditions using a combination of cluster mean-field analysis and Monte Carlo simulations. Interestingly, the fundamental diagram of current versus density shows a complex behavior with a unimodal dependence for attractions and weak repulsions that turns into the bimodal behavior for stronger repulsive interactions. Specific details of system-reservoir coupling for the open system have a strong effect on the stationary phases. We produce the steady-state phase diagrams for the bulk-adapted coupling to the reservoir using the minimum and maximum current principles. The strength and nature of interaction energy has a striking influence on the number of stationary phases. We observe that interactions lead to correlations having a strong impact on the system dynamical properties. The correlation between any two sites decays exponentially as the distance between the sites increases. Moreover, they are found to be short-range for repulsions and long-range for attractions. Our results also suggest that repulsions and attractions asymmetrically modify the dynamics of interacting particles in exclusion processes.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1742-5468/ac42cf
