Source author record

E. Abad

E. Abad appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

First-encounter time of two diffusing particles in two- and three-dimensional confinement

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 $\langle \cal{T}\rangle $, 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 \ll D_1$ and 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 $\langle \cal{T}\rangle$.

preprint2020arXiv

Continuous Time Random Walk in a velocity field: Role of domain growth, Galilei-invariant advection-diffusion, and kinetics of particle mixing

We consider the dynamics of a separable Continuous Time Random Walk (CTRW) when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding universe. The considered CTRW can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a Lévy flight. We first consider the case of zero velocity field. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of Lévy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: In one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter-controlled reactions in real systems are discussed.

preprint2020arXiv

First-encounter time of two diffusing particles in confinement

We investigate how confinement may drastically change both the probability density of the first-encounter time and the related 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 over the full time domain. We also evaluate the moments of the first-encounter time when they exist. We then turn to the case when the diffusivities are not equal, 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 translational invariance of such systems is broken.

preprint2020arXiv

Reaction-diffusion and reaction-subdiffusion equations on arbitrarily evolving domains

Reaction-diffusion equations are widely used as the governing evolution equations for modeling many physical, chemical, and biological processes. Here we derive reaction-diffusion equations to model transport with reactions on a one-dimensional domain that is evolving. The model equations, which have been derived from generalized continuous time random walks, can incorporate complexities such as subdiffusive transport and inhomogeneous domain stretching and shrinking. A method for constructing analytic expressions for short time moments of the position of the particles is developed and moments calculated from this approach are shown to compare favourably with results from random walk simulations and numerical integration of the reaction transport equation. The results show the important role played by the initial condition. In particular, it strongly affects the time dependence of the moments in the short time regime by introducing additional drift and diffusion terms. We also discuss how our reaction transport equation could be applied to study the spreading of a population on an evolving interface.

preprint2016arXiv

Anomalous diffusion and FRAP dynamics in the random comb model

We address the problem of diffusion on a comb whose teeth display a varying length. Specifically, the length $\ell$ of each tooth is drawn from a probability distribution displaying the large-$\ell$ behavior $P(\ell) \sim \ell^{-(1+α)}$ ($α>0$). Our method is based on the mean-field description provided by the well-tested CTRW approach for the random comb model, and the obtained analytical result for the diffusion coefficient is confirmed by numerical simulations. We subsequently incorporate retardation effects arising from binding/unbinding kinetics into our model and obtain a scaling law characterizing the corresponding change in the diffusion coefficient. Finally, our results for the diffusion coefficient are used as an input to compute concentration recovery curves mimicking FRAP experiments in comb-like geometries such as spiny dendrites. We show that such curves cannot be fitted perfectly by a model based on scaled Brownian motion, i.e., a standard diffusion equation with a time-dependent diffusion coefficient. However, differences between the exact curves and such fits are small, thereby providing justification for the practical use of models relying on scaled Brownian motion as a fitting procedure for recovery curves arising from particle diffusion in comb-like systems.

preprint2016arXiv

Diffusion in an expanding medium: Fokker-Planck equation, Green's function and first-passage properties

We present a classical, mesoscopic derivation of the Fokker-Planck equation for diffusion in an expanding medium. To this end, we take a conveniently generalized Chapman-Kolmogorov equation as the starting point. We obtain an analytical expression for the Green's function (propagator) and investigate both analytically and numerically how this function and the associated moments behave. We also study first-passage properties in expanding hyperspherical geometries. We show that in all cases the behavior is determined to a great extent by the so-called Brownian conformal time $τ(t)$, which we define via the relation $\dot τ=1/a^2$, where $a(t)$ is the expansion scale factor. If the medium expansion is driven by a power law [$a(t) \propto t^γ$ with $γ>0$], we find interesting crossover effects in the mixing effectiveness of the diffusion process when the characteristic exponent $γ$ is varied. Crossover effects are also found at the level of the survival probability and of the moments of the first passage-time distribution with two different regimes separated by the critical value $γ=1/2$. The case of an exponential scale factor is analyzed separately both for expanding and contracting media. In the latter situation, a stationary probability distribution arises in the long time limit.

preprint2015arXiv

Optimal search strategies of space-time coupled random walkers with finite lifetimes

We present a simple paradigm for detection of an immobile target by a space-time coupled random walker with a finite lifetime. The motion of the walker is characterized by linear displacements at a fixed speed and exponentially distributed duration, interrupted by random changes in the direction of motion and resumption of motion in the new direction with the same speed. We call these walkers "mortal creepers". A mortal creeper may die at any time during its motion according to an exponential decay law characterized by a finite mean death rate $ω_m$. While still alive, the creeper has a finite mean frequency $ω$ of change of the direction of motion. In particular, we consider the efficiency of the target search process, characterized by the probability that the creeper will eventually detect the target. Analytic results confirmed by numerical results show that there is an $ω_m$-dependent optimal frequency $ω=ω_{opt}$ that maximizes the probability of eventual target detection. We work primarily in one-dimensional ($d=1$) domains and examine the role of initial conditions and of finite domain sizes. Numerical results in $d=2$ domains confirm the existence of an optimal frequency of change of direction, thereby suggesting that the observed effects are robust to changes in dimensionality. In the $d=1$ case, explicit expressions for the probability of target detection in the long time limit are given. In the case of an infinite domain, we compute the detection probability for arbitrary times and study its early- and late-time behavior. We further consider the survival probability of the target in the presence of many independent creepers beginning their motion at the same location and at the same time. We also consider a version of the standard "target problem" in which many creepers start at random locations at the same time.

preprint2014arXiv

A reaction-subdiffusion model of fluorescence recovery after photobleaching (FRAP)

Anomalous diffusion, in particular subdiffusion, is frequently invoked as a mechanism of motion in dense biological media, and may have a significant impact on the kinetics of binding/unbinding events at the cellular level. In this work we incorporate anomalous diffusion in a previously developed model for FRAP experiments. Our particular implementation of subdiffusive transport is based on a continuous time random walk (CTRW) description of the motion of fluorescent particles, as CTRWs lend themselves particularly well to the inclusion of binding/unbinding events. In order to model switching between bound and unbound states of fluorescent subdiffusive particles, we derive a fractional reaction-subdiffusion equation of rather general applicability. Using suitable initial and boundary conditions, this equation is then incorporated in the model describing two-dimensional kinetics of FRAP experiments. We find that this model can be used to obtain excellent fits to experimental data. Moreover, recovery curves corresponding to different radii of the circular bleach spot can be fitted by a single set of parameters. While not enough evidence has been collected to claim with certainty that CTRW is the underlying transport mechanism in FRAP experiments, the compatibility of our results with experimental data fuels the discussion as to whether normal diffusion or anomalous diffusion is the appropriate model, and as to whether anomalous diffusion effects are important to fully understand the outcomes of FRAP experiments. On a more technical side, we derive explicit analytic solutions of our model in certain limits.

preprint2013arXiv

Arrival statistics and exploration properties of mortal walkers

We study some of the salient features of the arrival statistics and exploration properties of mortal random walkers, that is, walkers that may die as they move, or as they wait to move. Such evanescence or death events have profound consequences for quantities such as the number of distinct sites visited which are relevant for the computation of encounter-controlled rates in chemical kinetics. We exploit the observation that well-known methods developed decades ago for immortal walkers are widely applicable to mortal walkers. The particular cases of exponential and power-law evanescence are considered in detail. Finally, we discuss the relevance of our results to the target problem with mortal traps and a particular application thereof, namely, the defect diffusion model. Evanescence of defects is postulated as a possible complementary contribution or perhaps even an alternative to anomalous diffusion to explain observed stretched exponential relaxation behavior.

preprint2013arXiv

Evanescent continuous time random walks

We study how an evanescence process affects the number of distinct sites visited by a continuous time random walker in one dimension. We distinguish two very different cases, namely, when evanescence can only occur concurrently with a jump, and when evanescence can occur at any time. The first is characteristic of trapping processes on a lattice, whereas the second is associated with spontaneous death processes such as radioactive decay. In both of these situations we consider three different forms of the waiting time distribution between jumps, namely, exponential, long-tailed, and ultra-slow.

preprint2013arXiv

Exploration and Trapping of Mortal Random Walkers

Exploration and trapping properties of random walkers that may evanesce at any time as they walk have seen very little treatment in the literature, and yet a finite lifetime is a frequent occurrence, and its effects on a number of random walk properties may be profound. For instance, whereas the average number of distinct sites visited by an immortal walker grows with time without bound, that of a mortal walker may, depending on dimensionality and rate of evanescence, remain finite or keep growing with the passage of time. This number can in turn be used to calculate other classic quantities such as the survival probability of a target surrounded by diffusing traps. If the traps are immortal, the survival probability will vanish with increasing time. However, if the traps are evanescent, the target may be spared a certain death. We analytically calculate a number of basic and broadly used quantities for evanescent random walkers.

preprint2013arXiv

Quasi-resonant neutralization of He+ ions at a germanium surface

When low-energy He ions are scattered from a Ge surface, the fraction of positive ions exhibits characteristic oscillations as a function of ion energy. These oscillations are caused by quasi-resonant neutralization (qRN), a process which can be observed for materials with a narrow band nearly resonant with the unperturbed He-1s level. In the present manuscript, we employ recently developed theoretical methods to extract quantitative information on the efficiency of qRN. Our evaluation reveals that qRN is a highly efficient process leading to ion fractions two orders of magnitude lower than in systems for which neutralization is only due to Auger processes.

preprint2012arXiv

Survival probability of an immobile target in a sea of evanescent diffusive or subdiffusive traps: a fractional equation approach

We calculate the survival probability of an immobile target surrounded by a sea of uncorrelated diffusive or subdiffusive evanescent traps, i.e., traps that disappear in the course of their motion. Our calculation is based on a fractional reaction-subdiffusion equation derived from a continuous time random walk model of the system. Contrary to an earlier method valid only in one dimension (d=1), the equation is applicable in any Euclidean dimension d and elucidates the interplay between anomalous subdiffusive transport, the irreversible evanescence reaction and the dimension in which both the traps and the target are embedded. Explicit results for the survival probability of the target are obtained for a density ρ(t) of traps which decays (i) exponentially and (ii) as a power law. In the former case, the target has a finite asymptotic survival probability in all integer dimensions, whereas in the latter case there are several regimes where the values of the decay exponent for ρ(t) and the anomalous diffusion exponent of the traps determine whether or not the target has a chance of eternal survival in one, two and three dimensions.

preprint2010arXiv

Application of Fractional Calculus to Reaction-Subdiffusion Processes and Morphogen Gradient Formation

It is a well known fact that subdiffusion equations in terms of fractional derivatives can be obtained from Continuous Time Random Walk (CTRW) models with long-tailed waiting time distributions. Over the last years various authors have shown that extensions of such CTRW models incorporating reactive processes to the mesoscopic transport equations may lead to non-intuitive reaction-subdiffusion equations. In particular, one such equation has been recently derived for a subdiffusive random walker subject to a linear (first-order) death process. We take this equation as a starting point to study the developmental biology key problem of morphogen gradient formation, both for the uniform case where the morphogen degradation rate coefficient (reactivity) is constant and for the non-uniform case (position-dependent reactivity). In the uniform case we obtain exponentially decreasing stationary concentration profiles and we study their robustness with respect to perturbations in the incoming morphogen flux. In the non-uniform case we find a rich phenomenology at the level of the stationary profiles. We conclude that the analytic form of the long-time morphogen concentration profiles is very sensitive to the spatial dependence of the reactivity and the specific value of the anomalous diffusion coefficient.

preprint2009arXiv

Reaction-subdiffusion and reaction-superdiffusion equations for evanescent particles performing continuous time random walks

Starting from a continuous time random walk (CTRW) model of particles that may evanesce as they walk, our goal is to arrive at macroscopic integro-differential equations for the probability density for a particle to be found at point r at time t given that it started its walk from r_0 at time t=0. The passage from the CTRW to an integro-differential equation is well understood when the particles are not evanescent. Depending on the distribution of stepping times and distances, one arrives at standard macroscopic equations that may be "normal" (diffusion) or "anomalous" (subdiffusion and/or superdiffusion). The macroscopic description becomes considerably more complicated and not particularly intuitive if the particles can die during their walk. While such equations have been derived for specific cases, e.g., for location-independent exponential evanescence, we present a more general derivation valid under less stringent constraints than those found in the current literature.