An official website of the United States government
Cellular scale decision-making is modulated by the dynamics of signalling molecules and their diffusive trajectories from a source to small absorbing sites on the cellular surface. Diffusive capture problems which model this process are computationally challenging due to their complex geometry and mixed boundary conditions together with intrinsically long transients that occur before a particle is captured. This paper reports on a particle-based kinetic Monte Carlo (KMC) method that provides rapid accurate simulation of arrival statistics for (i) a half-space bounded by a surface with a finite collection of absorbing traps and (ii) the domain exterior to a convex cell, again with absorbing traps. We validate our method by replicating classical results and verifying some newly developed boundary homogenization theories and matched asymptotic expansions on capture rates. In the case of non-spherical domains, we describe a new shielding effect in which geometry can play a role in sharpening cellular estimates on the directionality of diffusive sources.
more »
« less
Boundary homogenization for patchy surfaces trapping patchy particles
Trapping diffusive particles at surfaces is a key step in many systems in chemical and biological physics. Trapping often occurs via reactive patches on the surface and/or the particle. The theory of boundary homogenization has been used in many prior works to estimate the effective trapping rate for such a system in the case that either (i) the surface is patchy and the particle is uniformly reactive or (ii) the particle is patchy and the surface is uniformly reactive. In this paper, we estimate the trapping rate for the case that the surface and the particle are both patchy. In particular, the particle diffuses translationally and rotationally and reacts with the surface when a patch on the particle contacts a patch on the surface. We first formulate a stochastic model and derive a five-dimensional partial differential equation describing the reaction time. We then use matched asymptotic analysis to derive the effective trapping rate, assuming that the patches are roughly evenly distributed and occupy a small fraction of the surface and the particle. This trapping rate involves the electrostatic capacitance of a four-dimensional duocylinder, which we compute using a kinetic Monte Carlo algorithm. We further use Brownian local time theory to derive a simple heuristic estimate of the trapping rate and show that it is remarkably close to the asymptotic estimate. Finally, we develop a kinetic Monte Carlo algorithm to simulate the full stochastic system and then use these simulations to confirm the accuracy of our trapping rate estimates and homogenization theory.
more »
« less
- PAR ID:
- 10421793
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 9
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 094104
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Symmetry-breaking in patch formation on triangular gold nanoparticles by asymmetric polymer graftingAbstract Synthesizing patchy particles with predictive control over patch size, shape, placement and number has been highly sought-after for nanoparticle assembly research, but is fraught with challenges. Here we show that polymers can be designed to selectively adsorb onto nanoparticle surfaces already partially coated by other chains to drive the formation of patchy nanoparticles with broken symmetry. In our model system of triangular gold nanoparticles and polystyrene-b-polyacrylic acid patch, single- and double-patch nanoparticles are produced at high yield. These asymmetric single-patch nanoparticles are shown to assemble into self-limited patch‒patch connected bowties exhibiting intriguing plasmonic properties. To unveil the mechanism of symmetry-breaking patch formation, we develop a theory that accurately predicts our experimental observations at all scales—from patch patterning on nanoparticles, to the size/shape of the patches, to the particle assemblies driven by patch‒patch interactions. Both the experimental strategy and theoretical prediction extend to nanoparticles of other shapes such as octahedra and bipyramids. Our work provides an approach to leverage polymer interactions with nanoscale curved surfaces for asymmetric grafting in nanomaterials engineering.more » « less
-
We use a multiple-scale expansion to average the wave action balance equation over an ensemble of sea-surface velocity fields characteristic of the ocean mesoscale and submesoscale. Assuming that the statistical properties of the flow are stationary and homogeneous, we derive an expression for a diffusivity tensor of surface-wave action density. The small parameter in this expansion is the ratio of surface current speed to gravity wave group speed. For isotropic currents, the action diffusivity is expressed in terms of the kinetic energy spectrum of the flow. A Helmholtz decomposition of the sea-surface currents into solenoidal (vortical) and potential (divergent) components shows that, to leading order, the potential component of the surface velocity field has no effect on the diffusivity of wave action: only the vortical component of the sea-surface velocity results in diffusion of surface-wave action. We validate our analytic results for the action diffusivity by Monte Carlo ray-tracing simulations through an ensemble of stochastic velocity fields.more » « less
-
Abstract Liquid–liquid or liquid–air interfaces provide interesting environments to study colloids and are ubiquitous in nature and industry, as well as relevant in applications involving emulsions and foams. They present a particularly intriguing environment for studying active particles which exhibit a host of phenomena not seen in passive systems. Active particles can also provide on‐demand controllability that greatly expands their use in future applications. However, research on active particles at interfaces is relatively rare compared to those at solid surfaces or in the bulk. Here, magnetically steerable active colloids at liquid–air interfaces that self‐propel by bubble production via the catalytic decomposition of chemical fuel in the liquid medium is presented. The bubble formation and dynamics of “patchy” colloids with a patch of catalytic coating on their surface is investigated and compared to more traditional Janus colloids with a hemispherical coating. The patchy colloids tend to produce smaller bubbles and undergo smoother motion which makes them beneficial for applications such as precise micro‐manipulation. This is demonstrated by manipulating and assembling patterns of passive spheres on a substrate as well as at an air–liquid interface. The propulsion and bubble formation of both the Janus and patchy colloids is characterized and it is found that previously proposed theories are insufficient to fully describe their motion and bubble bursting mechanism. Additionally, the colloids, which reside at the air–liquid interface, demonstrate novel interfacial positive gravitaxis towards the droplet edges which is attributed to a torque resulting from opposing downward and buoyant forces on the colloids.more » « less
-
We consider the approximation of expectations with respect to the distribution of a latent Markov process given noisy measurements. This is known as the smoothing problem and is often approached with particle and Markov chain Monte Carlo (MCMC) methods. These methods provide consistent but biased estimators when run for a finite time. We propose a simple way of coupling two MCMC chains built using Particle Independent Metropolis-Hastings (PIMH) to produce unbiased smoothing estimators. Unbiased estimators are appealing in the context of parallel computing, and facilitate the construction of confidence intervals. The proposed scheme only requires access to off-the-shelf Particle Filters (PF) and is thus easier to implement than recently proposed unbiased smoothers. The approach is demonstrated on a Lévy-driven stochastic volatility model and a stochastic kinetic model.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0135048
