Source author record

S. B. Yuste

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

ResearcherUnclaimed source record

cond-mat.stat-mech physics.chem-ph cond-mat.soft Biological Physics cond-mat.dis-nn math-ph math.MP physics.class-ph Subcellular Processes

Catalog footprint

What is connected

23works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

preprint2020arXiv

Structural properties of additive binary hard-sphere mixtures

An approach to obtain the structural properties of additive binary hard-sphere mixtures is presented. Such an approach, which is a nontrivial generalization of the one recently used for monocomponent hard-sphere fluids [S. Pieprzyk, A. C. Brańka, and D. M. Heyes, Phys. Rev. E 95, 062104 (2017)], combines accurate molecular-dynamics simulation data, the pole structure representation of the total correlation functions, and the Ornstein-Zernike equation. A comparison of the direct correlation functions obtained with the present scheme with those derived from theoretical results stemming from the Percus-Yevick (PY) closure and the so-called rational-function approximation (RFA) is performed. The density dependence of the leading poles of the Fourier transforms of the total correlation functions and the decay of the pair correlation functions of the mixtures are also addressed and compared to the predictions of the two theoretical approximations. A very good overall agreement between the results of the present scheme and those of the RFA is found, thus suggesting that the latter (which is an improvement over the PY approximation) can safely be used to predict reasonably well the long-range behavior, including the structural crossover, of the correlation functions of additive binary hard-sphere mixtures.

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.

preprint2016arXiv

Theoretical approaches to the structural properties of the square-shoulder fluid

A comparison of simulation results with the prediction of the structural properties of square-shoulder fluids is carried out to assess the performance of three theories: Tang--Lu's first-order mean spherical approximation, the simplified exponential approximation of the latter and the rational-function approximation. These three theoretical developments share the characteristic of being analytical in Laplace space and of reducing in the proper limit to the Percus--Yevick result for the hard-sphere fluid. Overall, the best agreement with the simulation data is obtained with the simplified exponential approximation.

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

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.

preprint2012arXiv

Rational-function approximation for fluids interacting via piece-wise constant potentials

The structural properties of fluids whose molecules interact via potentials with a hard-core plus n piece-wise constant sections of different widths and heights are derived using a (semi-analytical) rational-function approximation method. The results are illustrated for the cases of a square-shoulder plus square-well potential and a shifted square-well potential and compared both with simulation data and with those that follow from the (numerical) solutions of the Percus-Yevick integral equation.

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.

preprint2011arXiv

Asymptotic solutions of decoupled continuous-time random walks with superheavy-tailed waiting time and heavy-tailed jump length distributions

We study the long-time behavior of decoupled continuous-time random walks characterized by superheavy-tailed distributions of waiting times and symmetric heavy-tailed distributions of jump lengths. Our main quantity of interest is the limiting probability density of the position of the walker multiplied by a scaling function of time. We show that the probability density of the scaled walker position converges in the long-time limit to a non-degenerate one only if the scaling function behaves in a certain way. This function as well as the limiting probability density are determined in explicit form. Also, we express the limiting probability density which has heavy tails in terms of the Fox $H$-function and find its behavior for small and large distances.

preprint2011arXiv

Structure of the square-shoulder fluid

The structural properties of square-shoulder fluids are derived from the use of the rational function approximation method. The computation of both the radial distribution function and the static structure factor involves mostly analytical steps, requiring only the numerical solution of a single transcendental equation. The comparison with available simulation data and with numerical solutions of the Percus-Yevick and hypernetted-chain integral equations shows that the present approximation represents an improvement over the Percus-Yevick theory for this system and a reasonable compromise between accuracy and simplicity.

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

Coagulation reaction in low dimensions: Revisiting subdiffusive A+A reactions in one dimension

We present a theory for the coagulation reaction A+A -> A for particles moving subdiffusively in one dimension. Our theory is tested against numerical simulations of the concentration of $A$ particles as a function of time (``anomalous kinetics'') and of the interparticle distribution function as a function of interparticle distance and time. We find that the theory captures the correct behavior asymptotically and also at early times, and that it does so whether the particles are nearly diffusive or very subdiffusive. We find that, as in the normal diffusion problem, an interparticle gap responsible for the anomalous kinetics develops and grows with time. This corrects an earlier claim to the contrary on our part.

preprint2009arXiv

Mean Field Model of Coagulation and Annihilation Reactions in a Medium of Quenched Traps: Subdiffusion

We present a mean field model for coagulation ($A+A\to A$) and annihilation ($A+A\to 0$) reactions on lattices of traps with a distribution of depths reflected in a distribution of mean escape times. The escape time from each trap is exponentially distributed about the mean for that trap, and the distribution of mean escape times is a power law. Even in the absence of reactions, the distribution of particles over sites changes with time as particles are caught in ever deeper traps, that is, the distribution exhibits aging. Our main goal is to explore whether the reactions lead to further (time dependent) changes in this distribution.

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.

preprint2001arXiv

Comment on "Theory and computer simulation for the equation of state of additive hard-disk fluid mixtures"

A flaw in the comparison between two different theoretical equations of state for a binary mixture of additive hard disks and Monte Carlo results, as recently reported in C. Barrio and J. R. Solana, Phys. Rev. E 63, 011201 (2001), is pointed out. It is found that both proposals, which require the equation of state of the single component system as input, lead to comparable accuracy but the one advocated by us [A. Santos, S. B. Yuste, and M. López de Haro, Mol. Phys. 96, 1 (1999)] is simpler and complies with the exact limit in which the small disks are point particles.

preprint2001arXiv

Multiparticle trapping problem in the half-line

A variation of Rosenstock's trapping model in which $N$ independent random walkers are all initially placed upon a site of a one-dimensional lattice in the presence of a {\em one-sided} random distribution (with probability $c$) of absorbing traps is investigated. The probability (survival probability) $Φ_N(t)$ that no random walker is trapped by time $t$ for $N \gg 1$ is calculated by using the extended Rosenstock approximation. This requires the evaluation of the moments of the number $S_N(t)$ of distinct sites visited in a {\em given} direction up to time $t$ by $N$ independent random walkers. The Rosenstock approximation improves when $N$ increases, working well in the range $Dt\ln^2(1-c) \ll \ln N$, $D$ being the diffusion constant. The moments of the time (lifetime) before any trapping event occurs are calculated asymptotically, too. The agreement with numerical results is excellent.

preprint2000arXiv

Demixing in binary mixtures of hard hyperspheres

The phase behavior of binary fluid mixtures of hard hyperspheres in four and five dimensions is investigated. Spinodal instability is found by using a recent and accurate prescription for the equation of state of the mixture that requires the equation of state of the single component fluid as input. The role played by the dimensionality on the possible metastability of the demixing transition with respect to a fluid-solid transition is discussed. The binodal curves in the pressure--chemical potential representation are seen to lie on a common line, independent of the size ratio

S. B. Yuste

What is connected

Connect this record

See the researcher in context

Building this map preview

23 published item(s)

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

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

First-encounter time of two diffusing particles in confinement

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

Structural properties of additive binary hard-sphere mixtures

Anomalous diffusion and FRAP dynamics in the random comb model

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

Theoretical approaches to the structural properties of the square-shoulder fluid

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

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

Evanescent continuous time random walks

Exploration and Trapping of Mortal Random Walkers

Rational-function approximation for fluids interacting via piece-wise constant potentials

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

Asymptotic solutions of decoupled continuous-time random walks with superheavy-tailed waiting time and heavy-tailed jump length distributions

Structure of the square-shoulder fluid

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

Coagulation reaction in low dimensions: Revisiting subdiffusive A+A reactions in one dimension

Mean Field Model of Coagulation and Annihilation Reactions in a Medium of Quenched Traps: Subdiffusion

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

Comment on "Theory and computer simulation for the equation of state of additive hard-disk fluid mixtures"

Multiparticle trapping problem in the half-line

Demixing in binary mixtures of hard hyperspheres