An official website of the United States government
We develop a convergent reaction-drift-diffusion master equation (CRDDME) to facilitate the study of reaction processes in which spatial transport is influenced by drift due to one-body potential fields within general domain geometries. The generalized CRDDME is obtained through two steps. We first derive an unstructured grid jump process approximation for reversible diffusions, enabling the simulation of drift-diffusion processes where the drift arises due to a conservative field that biases particle motion. Leveraging the Edge-Averaged Finite Element method, our approach preserves detailed balance of drift-diffusion fluxes at equilibrium, and preserves an equilibrium Gibbs-Boltzmann distribution for particles undergoing drift-diffusion on the unstructured mesh. We next formulate a spatially-continuous volume reactivity particle-based reaction-drift-diffusion model for reversible reactions of the form A + B <–> C. A finite volume discretization is used to generate jump process approximations to reaction terms in this model. The discretization is developed to ensure the combined reaction-drift-diffusion jump process approximation is consistent with detailed balance of reaction fluxes holding at equilibrium, along with supporting a discrete version of the continuous equilibrium state. The new CRDDME model represents a continuous-time discrete-space jump process approximation to the underlying volume reactivity model. We demonstrate the convergence and accuracy of the new CRDDME through a number of numerical examples, and illustrate its use on an idealized model for membrane protein receptor dynamics in T cell signaling.
more »
« less
Detailed balance for particle models of reversible reactions in bounded domains
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
- Award ID(s):
- 1902854
- PAR ID:
- 10336262
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 20
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 204105
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
An exact pairwise hydrodynamic theory is developed for the flow-induced spatial distribution of particles in dilute polydisperse suspensions undergoing two-dimensional unidirectional flows, including shear and planar Poiseuille flows. Coupled diffusive fluxes and a drift velocity are extracted from a Boltzmann-like master equation. A boundary layer is predicted in regions where the shear rate vanishes with thickness set by the radii of the upstream collision cross-sections for pair interactions. An analysis of this region yields linearly vanishing drift velocities and non-vanishing diffusivities where the shear rate vanishes, thus circumventing the source of the singular particle distribution predicted by the usual models. Outside of the boundary layer, a power-law particle distribution is predicted with exponent equal to minus half the exponent of the local shear rate. Trajectories for particles with symmetry-breaking contact interactions (e.g. rough particles, permeable particles, emulsion drops) are analytically integrated to yield particle displacements given by quadratures of hard-sphere (or spherical drop) mobility functions. Using this analysis, stationary particle distributions are obtained for suspensions in Poiseuille flow. The scale for the particle distribution in monodisperse suspensions is set by the collision cross-section of the particles but its shape is almost universal. Results for polydisperse suspensions show size segregation in the central boundary layer with enrichment of smaller particles. Particle densities at the centreline scale approximately with the inverse square root of particle size. A superposition approximation reliably predicts the exact results over a broad range of parameters. The predictions agree with experiments in suspensions up to approximately 20 % volume fraction without fitting parameters.more » « less
-
Abstract We consider particle-based stochastic reaction-drift-diffusion models where particles move via diffusion and drift induced by one- and two-body potential interactions. The dynamics of the particles are formulated as measure-valued stochastic processes (MVSPs), which describe the evolution of the singular, stochastic concentration fields of each chemical species. The mean field large population limit of such models is derived and proven, giving coarse-grained deterministic partial integro-differential equations (PIDEs) for the limiting deterministic concentration fields’ dynamics. We generalize previous studies on the mean field limit of models involving only diffusive motion, with care to formulating the MVSP representation to ensure detailed balance of reversible reactions in the presence of potentials. Our work illustrates the more general set of PIDEs that arise in the mean field limit, demonstrating that the limiting macroscopic reactive interaction terms for reversible reactions obtain additional nonlinear concentration-dependent coefficients compared to the purely diffusive case. Numerical studies are presented which illustrate that two-body repulsive potential interactions can have a significant impact on the reaction dynamics, and also demonstrate the empirical numerical convergence of solutions to the PBSRDD model to the derived mean field PIDEs as the population size increases.more » « less
-
Ultrafast dynamics of CN radical reactions with chloroform solvent under vibrational strong couplingPolariton chemistry may provide a new means to control molecular reactivity, permitting remote, reversible modification of reaction energetics, kinetics, and product yields. A considerable body of experimental and theoretical work has already demonstrated that strong coupling between a molecular vibrational mode and the confined electromagnetic field of an optical cavity can alter chemical reactivity without external illumination. However, the mechanisms underlying cavity-altered chemistry remain unclear in large part because the experimental systems examined previously are too complex for detailed analysis of their reaction dynamics. Here, we experimentally investigate photolysis-induced reactions of cyanide radicals with strongly-coupled chloroform (CHCl3) solvent molecules and examine the intracavity rates of photofragment recombination, solvent complexation, and hydrogen abstraction. We use a microfluidic optical cavity fitted with dichroic mirrors to facilitate vibrational strong coupling (VSC) of the C–H stretching mode of CHCl3 while simultaneously permitting optical access at visible wavelengths. Ultrafast transient absorption experiments performed with cavities tuned on- and off-resonance reveal that VSC of the CHCl3 C–H stretching transition does not significantly modify any measured rate constants, including those associated with the hydrogen abstraction reaction. This work represents, to the best of our knowledge, the first experimental study of an elementary bimolecular reaction under VSC. We discuss how the conspicuous absence of cavity-altered effects in this system may provide insights into the mechanisms of modified ground state reactivity under VSC and help bridge the divide between experimental results and theoretical predictions in vibrational polariton chemistry.more » « less
-
Abstract Atmospheric aerosol particles impact Earth's radiation balance by acting as seeds for cloud droplet formation. Over half of global cloud seed particles are formed by nucleation, a process where gas‐phase compounds react to form stable particles. Reactions of sulfuric acid (SA) with a wide variety of atmospheric compounds have been previously shown to drive nucleation in the lower troposphere. However, global climate models poorly predict particle nucleation rates since current nucleation models do not describe nucleation for systems containing tens to hundreds of precursor compounds. The nucleation potential model (NPM) was recently developed to model SA nucleation of complex mixtures by measuring an effective base concentration using a 1‐nm condensation particle counter. This technique for estimating particle nucleation rates can be deployed at a much higher spatial and temporal resolution than current methods which require detailed knowledge of all nucleation reactions and measurements, typically using a mass spectrometer, of all nucleation precursor gases. This work expands NPM by showing that this model can capture enhancement and suppression of SA nucleation rates within a complex mixture of organic and inorganic acids, ambient air, and across a range of atmospherically relevant relative humidities. In addition, an expression for calculating atmospheric nucleation rates was also derived from the NPM. Ultimately, NPM provides a simple way to measure and model the extent compounds in a complex mixture enhance SA nucleation rates using a condensation particle counter.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0085296
