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.

Starting from a continuous-time random-walk (CTRW) model of particles that may evanesce as they walk, our goal is to arrive at macroscopic integrodifferential equations for the probability density for a particle to be found at point $\mathbf{r}$ at time $t$ given that it started its walk from ${\mathbf{r}}_0$ at time $t=0$. The passage from the CTRW to an integrodifferential 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.

We study gradient formation of subdiffusive morphogens. The morphogens are produced at a source point at a constant rate. From there they move subdiffusively and are also subject to degradation at a rate that may depend on location and on time. Our analysis is based on a reaction-subdiffusion equation obtained from a continuous time random-walk model with a long-tailed waiting time distribution that also incorporates an evanescence process. Spatially uniform degradation at a constant rate leads to an exponentially decreasing stationary concentration profile hardly distinguishable from that obtained with normal diffusion. On the other hand, with location-dependent degradation we find a rich gamut of profiles, some qualitatively quite different from those occurring with normal diffusion. We conclude that long-time morphogen concentration profiles are very sensitive to the spatial dependence of the reactivity and may also serve as a sensitive measure of the occurrence of anomalous diffusion.

Expanding media are typical in many different fields, e.g., in biology and cosmology. In general, a medium expansion (contraction) brings about dramatic changes in the behavior of diffusive transport properties such as the set of positional moments and the Green's function. Here, we focus on the characterization of such effects when the diffusion process is described by the continuous-time random-walk (CTRW) model. As is well known, when the medium is static this model yields anomalous diffusion for a proper choice of the probability density function (pdf) for the jump length and the waiting time, but the behavior may change drastically if a medium expansion is superimposed on the intrinsic random motion of the diffusing particle. For the case where the jump length and the waiting time pdfs are long-tailed, we derive a general bifractional diffusion equation which reduces to a normal diffusion equation in the appropriate limit. We then study some particular cases of interest, including Lévy flights and subdiffusive CTRWs. In the former case, we find an analytical exact solution for the Green's function (propagator). When the expansion is sufficiently fast, the contribution of the diffusive transport becomes irrelevant at long times and the propagator tends to a stationary profile in the comoving reference frame. In contrast, for a contracting medium a competition between the spreading effect of diffusion and the concentrating effect of contraction arises. In the specific case of a subdiffusive CTRW in an exponentially contracting medium, the latter effect prevails for sufficiently long times, and all the particles are eventually localized at a single point in physical space. This “big crunch” effect, totally absent in the case of normal diffusion, stems from inefficient particle spreading due to subdiffusion. We also derive a hierarchy of differential equations for the moments of the transport process described by the subdiffusive CTRW model in an expanding medium. From this hierarchy, the full time evolution of the second-order moment is obtained for some specific types of expansion. In the case of an exponential expansion, exact recurrence relations for the Laplace-transformed moments are obtained, whence the long-time behavior of moments of arbitrary order is subsequently inferred. Our analytical and numerical results for both Lévy flights and subdiffusive CTRWs confirm the intuitive expectation that the medium expansion hinders the mixing of diffusive particles occupying separate regions. In the case of Lévy flights, we quantify this effect by means of the so-called “Lévy horizon.”

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 ${\omega }_m$. While still alive, the creeper has a finite mean frequency $\omega$ 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 ${\omega }_m$-dependent optimal frequency $\omega ={\omega }_{\mathrm{opt}}$ that maximizes the probability of eventual target detection. We work primarily in one-dimensional $\left(d=1\right)$ 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.