1
|
Grebenkov DS. Spectral properties of the Dirichlet-to-Neumann operator for spheroids. Phys Rev E 2024; 109:055306. [PMID: 38907385 DOI: 10.1103/physreve.109.055306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Accepted: 05/01/2024] [Indexed: 06/24/2024]
Abstract
We study the spectral properties of the Dirichlet-to-Neumann operator and the related Steklov problem in spheroidal domains ranging from a needle to a disk. An explicit matrix representation of this operator for both interior and exterior problems is derived. We show how the anisotropy of spheroids affects the eigenvalues and eigenfunctions of the operator. As examples of physical applications, we discuss diffusion-controlled reactions on spheroidal partially reactive targets and the statistics of encounters between the diffusing particle and the spheroidal boundary.
Collapse
|
2
|
Chae SJ, Kim DW, Igoshin OA, Lee S, Kim JK. Beyond microtubules: The cellular environment at the endoplasmic reticulum attracts proteins to the nucleus, enabling nuclear transport. iScience 2024; 27:109235. [PMID: 38439967 PMCID: PMC10909898 DOI: 10.1016/j.isci.2024.109235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 01/03/2024] [Accepted: 02/09/2024] [Indexed: 03/06/2024] Open
Abstract
All proteins are translated in the cytoplasm, yet many, including transcription factors, play vital roles in the nucleus. While previous research has concentrated on molecular motors for the transport of these proteins to the nucleus, recent observations reveal perinuclear accumulation even in the absence of an energy source, hinting at alternative mechanisms. Here, we propose that structural properties of the cellular environment, specifically the endoplasmic reticulum (ER), can promote molecular transport to the perinucleus without requiring additional energy expenditure. Specifically, physical interaction between proteins and the ER impedes their diffusion and leads to their accumulation near the nucleus. This result explains why larger proteins, more frequently interacting with the ER membrane, tend to accumulate at the perinucleus. Interestingly, such diffusion in a heterogeneous environment follows Chapman's law rather than the popular Fick's law. Our findings suggest a novel protein transport mechanism arising solely from characteristics of the intracellular environment.
Collapse
Affiliation(s)
- Seok Joo Chae
- Department of Mathematical Sciences, KAIST, Daejeon 34141, Republic of Korea
- Biomedical Mathematics Group, Pioneer Research Center for Mathematical and Computational Sciences, Institute for Basic Science, Daejeon 34126, Republic of Korea
| | - Dae Wook Kim
- Department of Mathematics, University of Michigan, Ann Arbor, MI 48109, USA
| | - Oleg A. Igoshin
- Department of Bioengineering, Rice University, Houston, TX 77005, USA
- Center for Theoretical Biological Physics, Rice University, Houston, TX 77005, USA
- Department of Chemistry, Rice University, Houston, TX 77005, USA
- Department of Biosciences, Rice University, Houston, TX 77005, USA
| | - Seunggyu Lee
- Biomedical Mathematics Group, Pioneer Research Center for Mathematical and Computational Sciences, Institute for Basic Science, Daejeon 34126, Republic of Korea
- Division of Applied Mathematical Sciences, Korea University, Sejong 30019, Republic of Korea
| | - Jae Kyoung Kim
- Department of Mathematical Sciences, KAIST, Daejeon 34141, Republic of Korea
- Biomedical Mathematics Group, Pioneer Research Center for Mathematical and Computational Sciences, Institute for Basic Science, Daejeon 34126, Republic of Korea
| |
Collapse
|
3
|
Holehouse J, Redner S. First passage on disordered intervals. Phys Rev E 2024; 109:L032102. [PMID: 38632740 DOI: 10.1103/physreve.109.l032102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2023] [Accepted: 02/09/2024] [Indexed: 04/19/2024]
Abstract
We derive unexpected first-passage properties for nearest-neighbor hopping on finite intervals with disordered hopping rates, including (a) a highly variable spatial dependence of the first-passage time, (b) huge disparities in first-passage times for different realizations of hopping rates, (c) significant discrepancies between the first moment and the square root of the second moment of the first-passage time, and (d) bimodal first-passage time distributions. Our approach relies on the backward equation, in conjunction with probability generating functions, to obtain all moments, as well as the distribution of first-passage times. Our approach is simpler than previous approaches based on the forward equation, in which computing the mth moment of the first-passage time requires all preceding moments.
Collapse
Affiliation(s)
- James Holehouse
- The Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501, USA
| | - S Redner
- The Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501, USA
| |
Collapse
|
4
|
Go BG, Jeon E, Kim YW. Active search for a reactive target in thermal environments. J Chem Phys 2024; 160:044103. [PMID: 38258927 DOI: 10.1063/5.0179968] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Accepted: 12/15/2023] [Indexed: 01/24/2024] Open
Abstract
We study a stochastic process where an active particle, modeled by a one-dimensional run-and-tumble particle, searches for a target with a finite absorption strength in thermal environments. Solving the Fokker-Planck equation for a uniform initial distribution, we analytically calculate the mean searching time (MST), the time for the active particle to be finally absorbed, and show that there exists an optimal self-propulsion velocity of the active particle at which MST is minimized. As the diffusion constant increases, the optimal velocity changes from a finite value to zero, which implies that a purely diffusive Brownian motion outperforms an active motion in terms of searching time. Depending on the absorption strength of the target, the transition of the optimal velocity becomes either continuous or discontinuous, which can be understood based on the Landau approach. In addition, we obtain the phase diagram indicating the passive-efficient and the active-efficient regions. Finally, the initial condition dependence of MST is presented in limiting cases.
Collapse
Affiliation(s)
- Byeong Guk Go
- Department of Physics, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea
| | - Euijin Jeon
- Department of Physics, Technion-Israel Institute of Technology, Haifa 3200003, Israel
| | - Yong Woon Kim
- Department of Physics, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea
| |
Collapse
|
5
|
Grebenkov DS. Diffusion-Controlled Reactions: An Overview. Molecules 2023; 28:7570. [PMID: 38005291 PMCID: PMC10674959 DOI: 10.3390/molecules28227570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Revised: 10/20/2023] [Accepted: 10/21/2023] [Indexed: 11/26/2023] Open
Abstract
We review the milestones in the century-long development of the theory of diffusion-controlled reactions. Starting from the seminal work by von Smoluchowski, who recognized the importance of diffusion in chemical reactions, we discuss perfect and imperfect surface reactions, their microscopic origins, and the underlying mathematical framework. Single-molecule reaction schemes, anomalous bulk diffusions, reversible binding/unbinding kinetics, and many other extensions are presented. An alternative encounter-based approach to diffusion-controlled reactions is introduced, with emphasis on its advantages and potential applications. Some open problems and future perspectives are outlined.
Collapse
Affiliation(s)
- Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée, CNRS-Ecole Polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France
| |
Collapse
|
6
|
Niknam M, Bouchard LS. Nuclear induction lineshape modeling via hybrid SDE and MD approach. J Chem Phys 2023; 159:124201. [PMID: 38127390 DOI: 10.1063/5.0163782] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2023] [Accepted: 09/10/2023] [Indexed: 12/23/2023] Open
Abstract
The temperature dependence of the nuclear free induction decay in the presence of a magnetic-field gradient was found to exhibit motional narrowing in gases upon heating, a behavior that is opposite to that observed in liquids. This has led to the revision of the theoretical framework to include a more detailed description of particle trajectories since decoherence mechanisms depend on histories. In the case of free diffusion and single components, the new model yields the correct temperature trends. The inclusion of boundaries in the current formalism is not straightforward. We present a hybrid SDE-MD (stochastic differential equation - molecular dynamics) approach whereby MD is used to compute an effective viscosity and the latter is fed to the SDE to predict the line shape. The theory is in agreement with the experiments. This two-scale approach, which bridges the gap between short (molecular collisions) and long (nuclear induction) timescales, paves the way for the modeling of complex environments with boundaries, mixtures of chemical species, and intermolecular potentials.
Collapse
Affiliation(s)
- Mohamad Niknam
- Department of Chemistry and Biochemistry, University of California Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1059, USA and Center for Quantum Science and Engineering, UCLA, Los Angeles, California 90095-1059, USA
| | - Louis-S Bouchard
- Department of Chemistry and Biochemistry, University of California Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1059, USA and Center for Quantum Science and Engineering, UCLA, Los Angeles, California 90095-1059, USA
| |
Collapse
|
7
|
Kosztołowicz T. Subdiffusion equation with fractional Caputo time derivative with respect to another function in modeling transition from ordinary subdiffusion to superdiffusion. Phys Rev E 2023; 107:064103. [PMID: 37464604 DOI: 10.1103/physreve.107.064103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Accepted: 05/11/2023] [Indexed: 07/20/2023]
Abstract
We use a subdiffusion equation with fractional Caputo time derivative with respect to another function g (g-subdiffusion equation) to describe a smooth transition from ordinary subdiffusion to superdiffusion. Ordinary subdiffusion is described by the equation with the "ordinary" fractional Caputo time derivative, superdiffusion is described by the equation with a fractional Riesz-type spatial derivative. We find the function g for which the solution (Green's function, GF) to the g-subdiffusion equation takes the form of GF for ordinary subdiffusion in the limit of small time and GF for superdiffusion in the limit of long time. To solve the g-subdiffusion equation we use the g-Laplace transform method. It is shown that the scaling properties of the GF for g-subdiffusion and the GF for superdiffusion are the same in the long time limit. We conclude that for a sufficiently long time the g-subdiffusion equation describes superdiffusion well, despite a different stochastic interpretation of the processes. Then, paradoxically, a subdiffusion equation with a fractional time derivative describes superdiffusion. The superdiffusive effect is achieved here not by making anomalously long jumps by a diffusing particle, but by greatly increasing the particle jump frequency which is derived by means of the g-continuous-time random walk model. The g-subdiffusion equation is shown to be quite general, it can be used in modeling of processes in which a kind of diffusion change continuously over time. In addition, some methods used in modeling of ordinary subdiffusion processes, such as the derivation of local boundary conditions at a thin partially permeable membrane, can be used to model g-subdiffusion processes, even if this process is interpreted as superdiffusion.
Collapse
Affiliation(s)
- Tadeusz Kosztołowicz
- Institute of Physics, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| |
Collapse
|
8
|
Bressloff PC. Close encounters of the sticky kind: Brownian motion at absorbing boundaries. Phys Rev E 2023; 107:064121. [PMID: 37464709 DOI: 10.1103/physreve.107.064121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Accepted: 05/26/2023] [Indexed: 07/20/2023]
Abstract
Encounter-based models of diffusion provide a probabilistic framework for analyzing the effects of a partially absorbing reactive surface, in which the probability of absorption depends upon the amount of surface-particle contact time. In this paper we develop a class of encounter-based models that deal with absorption at sticky boundaries. Sticky boundaries occur in a diverse range of applications, including cell biology, colloidal physics, finance, and human crowd dynamics. They also naturally arise in active matter, where confined active particles tend to spontaneously accumulate at boundaries even in the absence of any particle-particle interactions. We begin by constructing a one-dimensional encounter-based model of sticky Brownian motion (BM), which is based on the zero-range limit of nonsticky BM with a short-range attractive potential well near the origin. In this limit, the boundary-contact time is given by the amount of time (occupation time) that the particle spends at the origin. We calculate the joint probability density or propagator for the particle position and the occupation time, and then identify an absorption event as the first time that the occupation time crosses a randomly generated threshold. We illustrate the theory by considering diffusion in a finite interval with a partially absorbing sticky boundary at one end. We show how various quantities, such as the mean first passage time (MFPT) for single-particle absorption and the relaxation to steady state at the multiparticle level, depend on moments of the random threshold distribution. Finally, we determine how sticky BM can be obtained by taking a particular diffusion limit of a sticky run-and-tumble particle (RTP).
Collapse
Affiliation(s)
- Paul C Bressloff
- Department of Mathematics, University of Utah, 155 South 1400 East, Salt Lake City, Utah 84112, USA
| |
Collapse
|
9
|
Grebenkov DS. Encounter-based approach to the escape problem. Phys Rev E 2023; 107:044105. [PMID: 37198799 DOI: 10.1103/physreve.107.044105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2023] [Accepted: 03/23/2023] [Indexed: 05/19/2023]
Abstract
We revise the encounter-based approach to imperfect diffusion-controlled reactions, which employs the statistics of encounters between a diffusing particle and the reactive region to implement surface reactions. We extend this approach to deal with a more general setting, in which the reactive region is surrounded by a reflecting boundary with an escape region. We derive a spectral expansion for the full propagator and investigate the behavior and probabilistic interpretations of the associated probability flux density. In particular, we obtain the joint probability density of the escape time and the number of encounters with the reactive region before escape, and the probability density of the first-crossing time of a prescribed number of encounters. We briefly discuss generalizations of the conventional Poissonian-type surface reaction mechanism described by Robin boundary condition and potential applications of this formalism in chemistry and biophysics.
Collapse
Affiliation(s)
- Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée, CNRS-Ecole Polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France
| |
Collapse
|
10
|
Plunkett CE, Lawley SD. Boundary homogenization for patchy surfaces trapping patchy particles. J Chem Phys 2023; 158:094104. [PMID: 36889970 DOI: 10.1063/5.0135048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Claire E Plunkett
- Department of Mathematics, University of Utah, Salt Lake City, Utah 84112, USA
| | - Sean D Lawley
- Department of Mathematics, University of Utah, Salt Lake City, Utah 84112, USA
| |
Collapse
|
11
|
Scher Y, Lauber Bonomo O, Pal A, Reuveni S. Microscopic theory of adsorption kinetics. J Chem Phys 2023; 158:094107. [PMID: 36889971 DOI: 10.1063/5.0121359] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
Adsorption is the accumulation of a solute at an interface that is formed between a solution and an additional gas, liquid, or solid phase. The macroscopic theory of adsorption dates back more than a century and is now well-established. Yet, despite recent advancements, a detailed and self-contained theory of single-particle adsorption is still lacking. Here, we bridge this gap by developing a microscopic theory of adsorption kinetics, from which the macroscopic properties follow directly. One of our central achievements is the derivation of the microscopic version of the seminal Ward-Tordai relation, which connects the surface and subsurface adsorbate concentrations via a universal equation that holds for arbitrary adsorption dynamics. Furthermore, we present a microscopic interpretation of the Ward-Tordai relation that, in turn, allows us to generalize it to arbitrary dimension, geometry, and initial conditions. The power of our approach is showcased on a set of hitherto unsolved adsorption problems to which we present exact analytical solutions. The framework developed herein sheds fresh light on the fundamentals of adsorption kinetics, which opens new research avenues in surface science with applications to artificial and biological sensing and to the design of nano-scale devices.
Collapse
Affiliation(s)
- Yuval Scher
- School of Chemistry, Center for the Physics and Chemistry of Living Systems, Ratner Institute for Single Molecule Chemistry, and the Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, 6997801 Tel Aviv, Israel
| | - Ofek Lauber Bonomo
- School of Chemistry, Center for the Physics and Chemistry of Living Systems, Ratner Institute for Single Molecule Chemistry, and the Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, 6997801 Tel Aviv, Israel
| | - Arnab Pal
- The Institute of Mathematical Sciences, CIT Campus, Taramani, Chennai 600113, India
| | - Shlomi Reuveni
- School of Chemistry, Center for the Physics and Chemistry of Living Systems, Ratner Institute for Single Molecule Chemistry, and the Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, 6997801 Tel Aviv, Israel
| |
Collapse
|
12
|
Bressloff PC. Renewal equations for single-particle diffusion through a semipermeable interface. Phys Rev E 2023; 107:014110. [PMID: 36797906 DOI: 10.1103/physreve.107.014110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 12/23/2022] [Indexed: 06/18/2023]
Abstract
Diffusion through semipermeable interfaces has a wide range of applications, ranging from molecular transport through biological membranes to reverse osmosis for water purification using artificial membranes. At the single-particle level, one-dimensional diffusion through a barrier with constant permeability κ_{0} can be modeled in terms of so-called snapping out Brownian motion (BM). The latter sews together successive rounds of partially reflected BMs that are restricted to either the left or right of the barrier. Each round is killed (absorbed) at the barrier when its Brownian local time exceeds an exponential random variable parameterized by κ_{0}. A new round is then immediately started in either direction with equal probability. It has recently been shown that the probability density for snapping out BM satisfies a renewal equation that relates the full density to the probability densities of partially reflected BM on either side of the barrier. Moreover, generalized versions of the renewal equation can be constructed that incorporate non-Markovian, encounter-based models of absorption. In this paper we extend the renewal theory of snapping out BM to single-particle diffusion in bounded domains and higher spatial dimensions. In each case we show how the solution of the renewal equation satisfies the classical diffusion equation with a permeable boundary condition at the interface. That is, the probability flux across the interface is continuous and proportional to the difference in densities on either side of the interface. We also consider an example of an asymmetric interface in which the directional switching after each absorption event is biased. Finally, we show how to incorporate an encounter-based model of absorption for single-particle diffusion through a spherically symmetric interface. We find that, even when the same non-Markovian model of absorption applies on either side of the interface, the resulting permeability is an asymmetric time-dependent function with memory. Moreover, the permeability functions tend to be heavy tailed.
Collapse
Affiliation(s)
- Paul C Bressloff
- Department of Mathematics, University of Utah, Salt Lake City, Utah 84112, USA
| |
Collapse
|
13
|
Kay T, Giuggioli L. Subdiffusion in the Presence of Reactive Boundaries: A Generalized Feynman-Kac Approach. JOURNAL OF STATISTICAL PHYSICS 2023; 190:92. [PMID: 37128546 PMCID: PMC10140114 DOI: 10.1007/s10955-023-03105-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/24/2022] [Accepted: 04/03/2023] [Indexed: 05/03/2023]
Abstract
We derive, through subordination techniques, a generalized Feynman-Kac equation in the form of a time fractional Schrödinger equation. We relate such equation to a functional which we name the subordinated local time. We demonstrate through a stochastic treatment how this generalized Feynman-Kac equation describes subdiffusive processes with reactions. In this interpretation, the subordinated local time represents the number of times a specific spatial point is reached, with the amount of time spent there being immaterial. This distinction provides a practical advance due to the potential long waiting time nature of subdiffusive processes. The subordinated local time is used to formulate a probabilistic understanding of subdiffusion with reactions, leading to the well known radiation boundary condition. We demonstrate the equivalence between the generalized Feynman-Kac equation with a reflecting boundary and the fractional diffusion equation with a radiation boundary. We solve the former and find the first-reaction probability density in analytic form in the time domain, in terms of the Wright function. We are also able to find the survival probability and subordinated local time density analytically. These results are validated by stochastic simulations that use the subordinated local time description of subdiffusion in the presence of reactions.
Collapse
Affiliation(s)
- Toby Kay
- Department of Engineering Mathematics, University of Bristol, Bristol, BS8 1UB UK
| | - Luca Giuggioli
- Department of Engineering Mathematics, University of Bristol, Bristol, BS8 1UB UK
- Bristol Centre for Complexity Sciences, University of Bristol, Bristol, BS8 1UB UK
| |
Collapse
|
14
|
Grebenkov DS, Skvortsov AT. Diffusion toward a nanoforest of absorbing pillars. J Chem Phys 2022; 157:244102. [PMID: 36586989 DOI: 10.1063/5.0132197] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Spiky coatings (also known as nanoforests or Fakir-like surfaces) have found many applications in chemical physics, material sciences, and biotechnology, such as superhydrophobic materials, filtration and sensing systems, and selective protein separation, to name but a few. In this paper, we provide a systematic study of steady-state diffusion toward a periodic array of absorbing cylindrical pillars protruding from a flat base. We approximate a periodic cell of this system by a circular tube containing a single pillar, derive an exact solution of the underlying Laplace equation, and deduce a simple yet exact representation for the total flux of particles onto the pillar. The dependence of this flux on the geometric parameters of the model is thoroughly analyzed. In particular, we investigate several asymptotic regimes, such as a thin pillar limit, a disk-like pillar, and an infinitely long pillar. Our study sheds light onto the trapping efficiency of spiky coatings and reveals the roles of pillar anisotropy and diffusional screening.
Collapse
Affiliation(s)
- Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée, CNRS-Ecole Polytechnique, Institut Polytechnique de Paris Paris, 91120 Palaiseau, France
| | - Alexei T Skvortsov
- Maritime Division, Defence Science and Technology Group, 506 Lorimer Street, Fishermans Bend, Port Melbourne, Victoria 3207, Australia
| |
Collapse
|
15
|
Bressloff PC. A probabilistic model of diffusion through a semi-permeable barrier. Proc Math Phys Eng Sci 2022. [DOI: 10.1098/rspa.2022.0615] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Diffusion through semi-permeable structures arises in a wide range of processes in the physical and life sciences. Examples at the microscopic level range from artificial membranes for reverse osmosis to lipid bilayers regulating molecular transport in biological cells to chemical and electrical gap junctions. There are also macroscopic analogues such as animal migration in heterogeneous landscapes. It has recently been shown that one-dimensional diffusion through a barrier with constant permeability
κ
0
is equivalent to snapping out Brownian motion (BM). The latter sews together successive rounds of partially reflecting BMs that are restricted to either the left or the right of the barrier. Each round is killed when its Brownian local time exceeds an exponential random variable parameterized by
κ
0
. A new round is then immediately started in either direction with equal probability. In this article, we use a combination of renewal theory, Laplace transforms and Green’s function methods to show how an extended version of snapping out BM provides a general probabilistic framework for modelling diffusion through a semi-permeable barrier. This includes modifications of the diffusion process away from the barrier (e.g. stochastic resetting) and non-Markovian models of membrane absorption that kill each round of partially reflected BM. The latter leads to time-dependent permeabilities.
Collapse
Affiliation(s)
- Paul C. Bressloff
- Department of Mathematics, University of Utah 155 South1400 East, Salt Lake City, UT 84112, USA
| |
Collapse
|
16
|
Benkhadaj Z, Grebenkov DS. Encounter-based approach to diffusion with resetting. Phys Rev E 2022; 106:044121. [PMID: 36397494 DOI: 10.1103/physreve.106.044121] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2022] [Accepted: 09/27/2022] [Indexed: 06/16/2023]
Abstract
An encounter-based approach consists in using the boundary local time as a proxy for the number of encounters between a diffusing particle and a target to implement various surface reaction mechanisms on that target. In this paper, we investigate the effects of stochastic resetting onto diffusion-controlled reactions in bounded confining domains. We first discuss the effect of position resetting onto the propagator and related quantities; in this way, we retrieve a number of earlier results but also provide complementary insights into them. Second, we introduce boundary local time resetting and investigate its impact. Curiously, we find that this type of resetting does not alter the conventional propagator governing the diffusive dynamics in the presence of a partially reactive target with a constant reactivity. In turn, the generalized propagator for other surface reaction mechanisms can be significantly affected. Our general results are illustrated for diffusion on an interval with reactive end points. Further perspectives and some open problems are discussed.
Collapse
Affiliation(s)
| | - Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91128 Palaiseau, France
| |
Collapse
|
17
|
Bressloff PC. Stochastically switching diffusion with partially reactive surfaces. Phys Rev E 2022; 106:034108. [PMID: 36266901 DOI: 10.1103/physreve.106.034108] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Accepted: 08/19/2022] [Indexed: 06/16/2023]
Abstract
In this paper we develop a hybrid version of the encounter-based approach to diffusion-mediated absorption at a reactive surface, which takes into account stochastic switching of a diffusing particle's conformational state. For simplicity, we consider a two-state model in which the probability of surface absorption depends on the current particle state and the amount of time the particle has spent in a neighborhood of the surface in each state. The latter is determined by a pair of local times ℓ_{n,t}, n=0,1, which are Brownian functionals that keep track of particle-surface encounters over the time interval [0,t]. We proceed by constructing a differential Chapman-Kolmogorov equation for a pair of generalized propagators P_{n}(x,ℓ_{0},ℓ_{1},t), where P_{n} is the joint probability density for the set (X_{t},ℓ_{0,t},ℓ_{1,t}) when N_{t}=n, where X_{t} denotes the particle position and N_{t} is the corresponding conformational state. Performing a double Laplace transform with respect to ℓ_{0},ℓ_{1} yields an effective system of equations describing diffusion in a bounded domain Ω, in which there is switching between two Robin boundary conditions on ∂Ω. The corresponding constant reactivities are κ_{j}=Dz_{j} and j=0,1, where z_{j} is the Laplace variable corresponding to ℓ_{j} and D is the diffusivity. Given the solution for the propagators in Laplace space, we construct a corresponding probabilistic model for partial absorption, which requires finding the inverse Laplace transform with respect to z_{0},z_{1}. We illustrate the theory by considering diffusion of a particle on the half-line with the boundary at x=0 effectively switching between a totally reflecting and a partially absorbing state. We calculate the flux due to absorption and use this to compute the resulting MFPT in the presence of a renewal-based stochastic resetting protocol. The latter resets the position and conformational state of the particle as well as the corresponding local times. Finally, we indicate how to extend the analysis to higher spatial dimensions using the spectral theory of Dirichlet-to-Neumann operators.
Collapse
Affiliation(s)
- Paul C Bressloff
- Department of Mathematics, University of Utah 155 South 1400 East, Salt Lake City, Utah 84112, USA
| |
Collapse
|
18
|
Abstract
Despite having been studied for decades, first passage processes remain an active area of research. In this contribution we examine a particle diffusing in an annulus with an inner absorbing boundary and an outer reflective boundary. We obtain analytic expressions for the joint distribution of the hitting time and the hitting angle in two and three dimensions. For certain configurations we observe a ``diffusive echo", i.e. two well-defined maxima in the first passage time distribution to a targeted position on the absorbing boundary. This effect, which results from the interplay between the starting location and the environmental constraints, may help to significantly increase the efficiency of the random search by generating a high, sustained flux to the targeted position over a short period. Finally, we examine the corresponding one-dimensional system for which there is no well-defined echo. In a confined system, the flux integrated over all target positions always displays a shoulder. This does not, however, guarantee the presence of an echo in the joint distribution.
Collapse
Affiliation(s)
| | - Julian Talbot
- Laboratoire de Physique Théorique de la Matière Condensée (LPTMC), Sorbonne Universités, France
| |
Collapse
|
19
|
Abstract
A probabilistic framework for studying single-particle diffusion in partially absorbing media has recently been developed in terms of an encounter-based approach. The latter computes the joint probability density (generalized propagator) for particle position
X
t
and a Brownian functional
U
t
that specifies the amount of time the particle is in contact with a reactive component
M
. Absorption occurs as soon as
U
t
crosses a randomly distributed threshold (stopping time). Laplace transforming the propagator with respect to
U
t
leads to a classical boundary value problem (BVP) in which the reactive component has a constant rate of absorption
z
, where
z
is the corresponding Laplace variable. Hence, a crucial step in the encounter-based approach is finding the inverse Laplace transform. In the case of a reactive boundary
∂
M
, this can be achieved by solving a classical Robin BVP in terms of the spectral decomposition of a Dirichlet-to-Neumann (D-to-N) operator on
∂
M
. In this paper, we develop the analogous construction in the case of a reactive substrate
M
. In particular, we show that the Laplace transformed propagator can be computed in terms of the spectral decomposition of a pair of D-to-N operators on
∂
M
. However, inverting the Laplace transform with respect to
z
is considerably more involved. We illustrate the theory by considering the D-to-N operators for some simple geometries.
Collapse
Affiliation(s)
- Paul C. Bressloff
- Department of Mathematics, University of Utah, 155 South 1400 East, Salt Lake City, UT 84112, USA
| |
Collapse
|
20
|
Grebenkov DS, Skvortsov AT. Mean first-passage time to a small absorbing target in three-dimensional elongated domains. Phys Rev E 2022; 105:054107. [PMID: 35706289 DOI: 10.1103/physreve.105.054107] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2022] [Accepted: 04/15/2022] [Indexed: 06/15/2023]
Abstract
We derive an approximate formula for the mean first-passage time (MFPT) to a small absorbing target of arbitrary shape inside an elongated domain of a slowly varying axisymmetric profile. For this purpose, the original Poisson equation in three dimensions is reduced to an effective one-dimensional problem on an interval with a semipermeable semiabsorbing membrane. The approximate formula captures correctly the dependence of the MFPT on the distance to the target, the radial profile of the domain, and the size and the shape of the target. This approximation is validated by Monte Carlo simulations.
Collapse
Affiliation(s)
- Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée, UMR No. 7643, CNRS, Ecole Polytechnique, IP Paris, 91120 Palaiseau, France
| | - Alexei T Skvortsov
- Maritime Division, Defence Science and Technology Group, 506 Lorimer Street, Fishermans Bend, Victoria 3207, Australia
| |
Collapse
|
21
|
Chaigneau A, Grebenkov DS. First-passage times to anisotropic partially reactive targets. Phys Rev E 2022; 105:054146. [PMID: 35706315 DOI: 10.1103/physreve.105.054146] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Accepted: 05/05/2022] [Indexed: 06/15/2023]
Abstract
We investigate restricted diffusion in a bounded domain towards a small partially reactive target in three- and higher-dimensional spaces. We propose a simple explicit approximation for the principal eigenvalue of the Laplace operator with mixed Robin-Neumann boundary conditions. This approximation involves the harmonic capacity and the surface area of the target, the volume of the confining domain, the diffusion coefficient, and the reactivity. The accuracy of the approximation is checked by using a finite-elements method. The proposed approximation determines also the mean first-reaction time, the long-time decay of the survival probability, and the overall reaction rate on that target. We identify the relevant lengthscale of the target, which determines its trapping capacity, and we investigate its relation to the target shape. In particular, we study the effect of target anisotropy on the principal eigenvalue by computing the harmonic capacity of prolate and oblate spheroids in various space dimensions. Some implications of these results in chemical physics and biophysics are briefly discussed.
Collapse
Affiliation(s)
- Adrien Chaigneau
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91120 Palaiseau, France
| | - Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91120 Palaiseau, France
| |
Collapse
|
22
|
Grebenkov DS. Depletion of resources by a population of diffusing species. Phys Rev E 2022; 105:054402. [PMID: 35706291 DOI: 10.1103/physreve.105.054402] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2022] [Accepted: 04/19/2022] [Indexed: 06/15/2023]
Abstract
Depletion of natural and artificial resources is a fundamental problem and a potential cause of economic crises, ecological catastrophes, and death of living organisms. Understanding the depletion process is crucial for its further control and optimized replenishment of resources. In this paper, we investigate a stock depletion by a population of species that undergo an ordinary diffusion and consume resources upon each encounter with the stock. We derive the exact form of the probability density of the random depletion time, at which the stock is exhausted. The dependence of this distribution on the number of species, the initial amount of resources, and the geometric setting is analyzed. Future perspectives and related open problems are discussed.
Collapse
Affiliation(s)
- Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91128 Palaiseau, France
| |
Collapse
|
23
|
Le Vot F, Yuste SB, Abad E, Grebenkov DS. First-encounter time of two diffusing particles in two- and three-dimensional confinement. Phys Rev E 2022; 105:044119. [PMID: 35590615 DOI: 10.1103/physreve.105.044119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2022] [Accepted: 03/15/2022] [Indexed: 06/15/2023]
Abstract
The statistics of the first-encounter time of diffusing particles changes drastically when they are placed under confinement. In the present work, we make use of Monte Carlo simulations to study the behavior of a two-particle system in two- and three-dimensional domains with reflecting boundaries. Based on the outcome of the simulations, we give a comprehensive overview of the behavior of the survival probability S(t) and the associated first-encounter time probability density H(t) over a broad time range spanning several decades. In addition, we provide numerical estimates and empirical formulas for the mean first-encounter time 〈T〉, as well as for the decay time T characterizing the monoexponential long-time decay of the survival probability. Based on the distance between the boundary and the center of mass of two particles, we obtain an empirical lower bound t_{B} for the time at which S(t) starts to significantly deviate from its counterpart for the no boundary case. Surprisingly, for small-sized particles, the dominant contribution to T depends only on the total diffusivity D=D_{1}+D_{2}, in sharp contrast to the one-dimensional case. This contribution can be related to the Wiener sausage generated by a fictitious Brownian particle with diffusivity D. In two dimensions, the first subleading contribution to T is found to depend weakly on the ratio D_{1}/D_{2}. We also investigate the slow-diffusion limit when D_{2}≪D_{1}, and we discuss the transition to the limit when one particle is a fixed target. Finally, we give some indications to anticipate when T can be expected to be a good approximation for 〈T〉.
Collapse
Affiliation(s)
- F Le Vot
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, E-06071 Badajoz, Spain
| | - S B Yuste
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, E-06071 Badajoz, Spain
| | - E Abad
- Departamento de Física Aplicada and Instituto de Computación Científica Avanzada (ICCAEx), Centro Universitario de Mérida, Universidad de Extremadura, E-06800 Mérida, Spain
| | - D S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91128 Palaiseau, France
| |
Collapse
|
24
|
Bressloff PC. Narrow capture problem: An encounter-based approach to partially reactive targets. Phys Rev E 2022; 105:034141. [PMID: 35428145 DOI: 10.1103/physreve.105.034141] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2022] [Accepted: 03/09/2022] [Indexed: 06/14/2023]
Abstract
A general topic of current interest is the analysis of diffusion problems in singularly perturbed domains with small interior targets or traps (the narrow capture problem). One major application is to intracellular diffusion, where the targets typically represent some form of reactive biochemical substrate. Most studies of the narrow capture problem treat the target boundaries as totally absorbing (Dirichlet), that is, the chemical reaction occurs immediately on first encounter between particle and target surface. In this paper, we analyze the three-dimensional narrow capture problem in the more realistic case of partially reactive target boundaries. We begin by considering classical Robin boundary conditions. Matching inner and outer solutions of the single-particle probability density, we derive an asymptotic expansion of the Laplace transformed flux into each reactive surface in powers of ε, where ερ is a given target size. In turn, the fluxes determine the splitting probabilities for target absorption. We then extend our analysis to more general types of reactive targets by combining matched asymptotic analysis with an encounter-based formulation of diffusion-mediated surface reactions. That is, we derive an asymptotic expansion of the joint probability density for particle position and the so-called boundary local time, which characterizes the amount of time that a Brownian particle spends in the neighborhood of a point on a totally reflecting boundary. The effects of surface reactions are then incorporated via an appropriate stopping condition for the boundary local time. Robin boundary conditions are recovered in the special case of an exponential law for the stopping local times. Finally, we illustrate the theory by exploring how the leading-order contributions to the splitting probabilities depend on the choice of surface reactions. In particular, we show that there is an effective renormalization of the target radius of the form ρ→ρ-Ψ[over ̃](1/ρ), where Ψ[over ̃] is the Laplace transform of the stopping local time distribution.
Collapse
Affiliation(s)
- Paul C Bressloff
- Department of Mathematics, University of Utah, 155 South 1400 East, Salt Lake City, Utah 84112, USA
| |
Collapse
|
25
|
Grebenkov DS, Kumar A. Reversible target-binding kinetics of multiple impatient particles. J Chem Phys 2022; 156:084107. [DOI: 10.1063/5.0083849] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Certain biochemical reactions can only be triggered after binding a sufficient number of particles to a specific target region such as an enzyme or a protein sensor. We investigate the distribution of the reaction time, i.e., the first instance when all independently diffusing particles are bound to the target. When each particle binds irreversibly, this is equivalent to the first-passage time of the slowest (last) particle. In turn, reversible binding to the target renders the problem much more challenging and drastically changes the distribution of the reaction time. We derive the exact solution of this problem and investigate the short-time and long-time asymptotic behaviors of the reaction time probability density. We also analyze how the mean reaction time depends on the unbinding rate and the number of particles. Our exact and asymptotic solutions are compared to Monte Carlo simulations.
Collapse
Affiliation(s)
- Denis S. Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS–Ecole Polytechnique, IP Paris, 91120 Palaiseau, France
| | - Aanjaneya Kumar
- Department of Physics, Indian Institute of Science Education and Research, Dr. Homi Bhabha Road, Pune 411008, India
| |
Collapse
|
26
|
Rast MP. Contact statistics in populations of noninteracting random walkers in two dimensions. Phys Rev E 2022; 105:014103. [PMID: 35193305 DOI: 10.1103/physreve.105.014103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Accepted: 12/09/2021] [Indexed: 11/07/2022]
Abstract
The interaction between individuals in biological populations, dilute components of chemical systems, or particles transported by turbulent flows depends critically on their contact statistics. This work clarifies those statistics under the simplifying assumptions that the underlying motions approximate a Brownian random walk and that the particles can be treated as noninteracting. We measure the contact-interval (also called the waiting-time or interarrival-time), contact-count, and contact-duration distributions in populations of individuals undergoing noninteracting continuous-space-time random walks on a periodic two-dimensional plane (a torus) as functions of the population number density, walker radius, and random-walk step size. The contact-interval is exponentially distributed for times longer than the ballistic mean-free-collision time but not for times shorter than that, and the contact duration distribution is strongly peaked at the ballistic-crossing time for head-on collisions when the ballistic-crossing time is short compared to the mean step duration. While successive contacts between individuals are independent, the probability of repeat contact decreases with time after a previous contact. This leads to a negative duration dependence of the waiting-time interval and overdispersion of the contact-count probability density function for all time intervals. The paper demonstrates that for populations of small particles (with a walker radius that is small compared to the mean-separation or random-walk step size), the ballistic mean-free-collision interval, the ballistic-crossing time, and the random-walk-step duration can be used to construct temporal scalings which allow for common waiting-time, contact-count, and contact-duration distributions across different populations. Semi-analytic approximations for both the waiting-time and contact-duration distributions are also presented.
Collapse
Affiliation(s)
- Mark Peter Rast
- Department of Astrophysical and Planetary Sciences, Laboratory for Atmospheric and Space Physics, University of Colorado, Boulder, Colorado 80309, USA
| |
Collapse
|
27
|
Multiple Novels and Accurate Traveling Wave and Numerical Solutions of the (2+1) Dimensional Fisher-Kolmogorov- Petrovskii-Piskunov Equation. MATHEMATICS 2021. [DOI: 10.3390/math9121440] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The analytical and numerical solutions of the (2+1) dimensional, Fisher-Kolmogorov-Petrovskii-Piskunov ((2+1) D-Fisher-KPP) model are investigated by employing the modified direct algebraic (MDA), modified Kudryashov (MKud.), and trigonometric-quantic B-spline (TQBS) schemes. This model, which arises in population genetics and nematic liquid crystals, describes the reaction–diffusion system by traveling waves in population genetics and the propagation of domain walls, pattern formation in bi-stable systems, and nematic liquid crystals. Many novel analytical solutions are constructed. These solutions are used to evaluate the requested numerical technique’s conditions. The numerical solutions of the considered model are studied, and the absolute value of error between analytical and numerical is calculated to demonstrate the matching between both solutions. Some figures are represented to explain the obtained analytical solutions and the match between analytical and numerical results. The used schemes’ performance shows their effectiveness and power and their ability to handle many nonlinear evolution equations.
Collapse
|
28
|
Kosztołowicz T. Subdiffusion in a system with a partially permeable partially absorbing wall. Phys Rev E 2021; 103:022117. [PMID: 33736008 DOI: 10.1103/physreve.103.022117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Accepted: 01/14/2021] [Indexed: 11/07/2022]
Abstract
We consider subdiffusion of a particle in a one-dimensional system with a thin partially permeable and partially absorbing wall. The system with the wall can be used to filter diffusing particles. Passing through the wall, the particle can be absorbed with a certain probability. Knowing the Green's functions we derive boundary conditions at the wall. The boundary conditions take a specific form in which fractional time derivatives are involved. The temporal evolution of the probability that a diffusing particle has not been absorbed is also considered.
Collapse
Affiliation(s)
- Tadeusz Kosztołowicz
- Institute of Physics, Jan Kochanowski University, Uniwersytecka 7, 25-406 Kielce, Poland
| |
Collapse
|
29
|
Le Vot F, Yuste SB, Abad E, Grebenkov DS. First-encounter time of two diffusing particles in confinement. Phys Rev E 2020; 102:032118. [PMID: 33076026 DOI: 10.1103/physreve.102.032118] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Accepted: 08/20/2020] [Indexed: 11/07/2022]
Abstract
We investigate how confinement may drastically change both the probability density of the first-encounter time and the associated survival probability in the case of two diffusing particles. To obtain analytical insights into this problem, we focus on two one-dimensional settings: a half-line and an interval. We first consider the case with equal particle diffusivities, for which exact results can be obtained for the survival probability and the associated first-encounter time density valid over the full time domain. We also evaluate the moments of the first-encounter time when they exist. We then turn to the case with unequal diffusivities and focus on the long-time behavior of the survival probability. Our results highlight the great impact of boundary effects in diffusion-controlled kinetics even for simple one-dimensional settings, as well as the difficulty of obtaining analytic results as soon as the translational invariance of such systems is broken.
Collapse
Affiliation(s)
- F Le Vot
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx) Universidad de Extremadura, E-06071 Badajoz, Spain
| | - S B Yuste
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx) Universidad de Extremadura, E-06071 Badajoz, Spain
| | - E Abad
- Departamento de Física Aplicada and Instituto de Computación Científica Avanzada (ICCAEx) Centro Universitario de Mérida Universidad de Extremadura, E-06800 Mérida, Spain
| | - D S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91128 Palaiseau, France and Institute of Physics & Astronomy, University of Potsdam, 14476 Potsdam-Golm, Germany
| |
Collapse
|
30
|
Grebenkov DS. Surface hopping propagator: An alternative approach to diffusion-influenced reactions. Phys Rev E 2020; 102:032125. [PMID: 33075930 DOI: 10.1103/physreve.102.032125] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 08/28/2020] [Indexed: 06/11/2023]
Abstract
Dynamics of a particle diffusing in a confinement can be seen a sequence of bulk-diffusion-mediated hops on the confinement surface. Here, we investigate the surface hopping propagator that describes the position of the diffusing particle after a prescribed number of encounters with that surface. This quantity plays the central role in diffusion-influenced reactions and determines their most common characteristics such as the propagator, the first-passage time distribution, and the reaction rate. We derive explicit formulas for the surface hopping propagator and related quantities for several Euclidean domains: half-space, circular annuli, circular cylinders, and spherical shells. These results provide the theoretical ground for studying diffusion-mediated surface phenomena. The behavior of the surface hopping propagator is investigated for both "immortal" and "mortal" particles.
Collapse
Affiliation(s)
- Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91128 Palaiseau, France
| |
Collapse
|