Search

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 ω = ωoptl 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. [ABSTRACT FROM AUTHOR]

We mobilize both a generating function approach and the theory of finite Markov processes to compute the probability of irreversible absorption of a randomly diffusing species on a lattice with competing reaction centers. We consider an N-site lattice populated by a single deep trap, and N - 1 partially absorbing traps (absorption probability 0 < s < 1). The influence of competing reaction centers on the probability of reaction at a target site (the deep trap) and the mean walk length of the random walker before localization (a measure of the reaction efficiency) are computed for different geometries. Both analytic expressions and numerical results are given for reactive processes on two-dimensional surfaces of Euler characteristic Ω = 0 and Ω = 2. The results obtained allow a characterization of catalyst deactivation processes on planar surfaces and on catalyst pellets where only a single catalytic site remains fully active (deep trap), the other sites being only partially active as a result of surface poisoning. The central result of our study is that the predicted dependence of the reaction efficiency on system size N and on s is in qualitative accord with previously reported experimental results, notably catalysts exhibiting selective poisoning due to surface sites that have different affinities for chemisorption of the poisoning agent (e.g., acid zeolite catalysts). Deviations from the efficiency of a catalyst with identical sites are quantified, and we find that such deviations display a significant dependence on the topological details of the surface (for fixed values of N and s we find markedly different results for, say, a planar surface and for the polyhedral surface of a catalyst pellet). Our results highlight the importance of surface topology for the efficiency of catalytic conversion processes on inhomogeneous substrates, and in particular for those aimed at industrial applications. From our exact analysis we extract results for the two limiting cases s ≈ 1 and s ≈ 0, corresponding respectively to weak and strong catalyst poisoning (decreasing s leads to a monotonic decrease in the efficiency of catalytic conversion). The results for the s ≈ 0 case are relevant for the dual problem of light-energy conversion via trapping of excitations in the chlorophyll antenna network. Here, decreasing the probability of excitation trapping s at sites other than the target molecule does not result in a decrease of the efficiency as in the catalyst case, but rather in enhanced efficiency of light-energy conversion, which we characterize in terms of N and s. The one-dimensional case and its connection with a modified version of the gambler's ruin problem are discussed. Finally, generalizations of our model are described briefly. [ABSTRACT FROM AUTHOR]

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 ultraslow. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 p(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 p(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. [ABSTRACT FROM AUTHOR]

Recent scanning tunneling hydrogen microscopy (STHM) experiments on PTCDA (perylene-3,4,9,10-tetracarboxylic-3,4,9,10-dianhydride)/Au(lll) have shown unprecedented intramolecular and intermolecular spatial resolution. The origin of this resolution is studied using an accurate STHM theoretical simulation technique that includes a detailed description of the electronic structure of both the tip and sample. Our results show that H2 molecules are dissociated on the Au típ; the adsorbed H atoms change the density of states at the Fermi level (Ep) of the tip, increasing its p-orbital character and reducing the i-orbital contribution. Also, due to the interaction with the H-decorated tip, Ep is shifted to the middle of the PTCDA lowest unoccupied molecular orbital peak, increasing dramatically the density of states of the sample at Ep. These effects give rise to the enhanced STHM resolution. [ABSTRACT FROM AUTHOR]

We consider a mortal random walker on a family of hierarchical graphs in the presence of some trap sites. The configuration comprising the graph, the starting point of the walk, and the locations of the trap sites is taken to be exactly self-similar as one goes from one generation of the family to the next. Under these circumstances, the total probability that the walker hits a trap is determined exactly as a function of the single-step survival probability q of the mortal walker. On the n th generation graph of the family, this probability is shown to be given by the n th iterate of a certain scaling function or map q → f (q). The properties of the map then determine, in each case, the behavior of the trapping probability, the mean time to trapping, the temporal scaling factor governing the random walk dimension on the graph, and other related properties. The formalism is illustrated for the cases of a linear hierarchical lattice and the Sierpinski graphs in two and three Euclidean dimensions. We find an effective reduction of the random walk dimensionality due to the ballistic behavior of the surviving particles induced by the mortality constraint. The relevance of this finding for experiments involving travel times of particles in diffusion-decay systems is discussed. [ABSTRACT FROM AUTHOR]

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". [ABSTRACT FROM AUTHOR]

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 τ =1/a², where a(t) is the expansion scale factor. If the medium expansion is driven by a power law [a(t)∝tγ with γ>0], then 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. [ABSTRACT FROM AUTHOR]

We address the problem of diffusion on a comb whose teeth display varying lengths. Specifically, the length ℓ of each tooth is drawn from a probability distribution displaying power law behavior at large ℓ,P(ℓ)~ ℓ1 (α>0). To start with, we focus on the computation of the anomalous diffusion coefficient for the subdiffusive motion along the backbone. This quantity is subsequently used as an input to compute concentration recovery curves mimicking fluorescence recovery after photobleaching experiments in comblike geometries such as spiny dendrites. Our method is based on the mean-field description provided by the well-tested continuous time random-walk approach for the random-comb model, and the obtained analytical result for the diffusion coefficient is confirmed by numerical simulations of a random walk with finite steps in time and space along the backbone and the teeth. 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, we show that recovery curves obtained with the help of the analytical expression for the anomalous diffusion coefficient 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 comblike systems. [ABSTRACT FROM AUTHOR]

A comparative study of Auger neutralization (AN) of He ions at noble metal surfaces is presented in order to reveal how the electronic structure of the sample influences this charge exchange process. Comparison of calculated ion fractions to experimental data obtained in low energy ion scattering (LEIS) shows that good agreement is achieved only if the relevant aspects of the He-metal interaction are properly taken into consideration. For instance AN depends sensitively on the distance-dependent position of the projectile level, which varies significantly when considering different target materials. [ABSTRACT FROM AUTHOR]

We study normal diffusive and subdiffusive processes in a harmonic potential (Ornstein-Uhlenbeck process) on a uniformly growing or contracting domain. Our starting point is a recently derived fractional Fokker-Planck equation, which covers both the case of Brownian diffusion and the case of a subdiffusive continuous-time random walk (CTRW). We find a high sensitivity of the random walk properties to the details of the domain growth rate, which gives rise to a variety of regimes with extremely different behaviors. At the origin of this rich phenomenology is the fact that the walkers still move while they wait to jump, since they are dragged by the deterministic drift arising from the domain growth. Thus, the increasingly long waiting times associated with the aging of the subdiffusive CTRW imply that, in the time interval between two consecutive jumps, the walkers might travel over much longer distances than in the normal diffusive case. This gives rise to seemingly counterintuitive effects. For example, on a static domain, both Brownian diffusion and subdiffusive CTRWs yield a stationary particle distribution with finite width when a harmonic potential is at play, thus indicating a confinement of the diffusing particle. However, for a sufficiently fast growing or contracting domain, this qualitative behavior breaks down, and differences between the Brownian case and the subdiffusive case are found. In the case of Brownian particles, a sufficiently fast exponential domain growth is needed to break the confinement induced by the harmonic force; in contrast, for subdiffusive particles such a breakdown may already take place for a sufficiently fast power-law domain growth. Our analytic and numerical results for both types of diffusion are fully confirmed by random walk simulations. [ABSTRACT FROM AUTHOR]

The kinetics of encounter-controlled processes in growing domains is markedly different from that in a static domain. Here we consider the specific example of diffusion-limited coalescence and annihilation reactions in one-dimensional space. In the static case, such reactions are among the few systems amenable to exact solution, which can be obtained by means of a well-known method of intervals. In the case of a uniformly growing domain, we show that a double transformation in time and space allows one to extend this method to compute the main quantities characterizing the spatial and temporal behavior. We show that a sufficiently fast domain growth brings about drastic changes in the behavior. In this case, the reactions stop prematurely, as a result of which the survival probability of the reacting particles tends to a finite value at long times and their spatial distribution freezes before reaching the fully self-ordered state. We obtain exact results for the survival probability and for key properties characterizing the degree of self-ordering induced by the chemical reactions, i.e., the interparticle distribution function and the pair correlation function. These results are confirmed by numerical simulations. [ABSTRACT FROM AUTHOR]

Fullerenes on single-layer epitaxial graphene are a model system to study very faint interactions at a molecular level. By a variable temperature scanning tunneling microscope we have been able to study ordered fullerene layers at 40 K, exclusively bound by van der Waals interactions. The experimentally determined adsorption geometry of the molecules is computationally confirmed only if van der Waals interactions are included in the calculation formalism. The relative orientation of fullerenes in their close-packed arrangement is found to be the crucial factor for determining the total energy. Observation of collective movements of fullerene islands points out the weak coupling to the substrate and the important role of the van der Waals cohesion forces within. [ABSTRACT FROM AUTHOR]

The Enskog kinetic theory is applied to compute the mean square displacement of impurities or intruders (modeled as smooth inelastic hard spheres) immersed in a granular gas of smooth inelastic hard spheres (grains). Both species (intruders and grains) are surrounded by an interstitial molecular gas (background) that plays the role of a thermal bath. The influence of the latter on the motion of intruders and grains is modeled via a standard viscous drag force supplemented by a stochastic Langevin-like force proportional to the background temperature. We solve the corresponding Enskog-Lorentz kinetic equation by means of the Chapman-Enskog expansion truncated to first order in the gradient of the intruder number density. The integral equation for the diffusion coefficient is solved by considering the first two Sonine approximations. To test these results, we also compute the diffusion coefficient from the numerical solution of the inelastic Enskog equation by means of the direct simulation Monte Carlo method. We find that the first Sonine approximation generally agrees well with the simulation results, although significant discrepancies arise when the intruders become lighter than the grains. Such discrepancies are largely mitigated by the use of the second Sonine approximation, in excellent agreement with computer simulations even for moderately strong inelasticities and/or dissimilar mass and diameter ratios. We invoke a random walk picture of the intruders' motion to shed light on the physics underlying the intricate dependence of the diffusion coefficient on the main system parameters. This approach, recently employed to study the case of an intruder immersed in a granular gas, also proves useful in the present case of a granular suspension. Finally, we discuss the applicability of our model to real systems in the self-diffusion case. We conclude that collisional effects may strongly impact the diffusion coefficient of the grains.

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&#95;{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&#95;{1}+D&#95;{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&#95;{1}/D&#95;{2}. We also investigate the slow-diffusion limit when D&#95;{2}≪D&#95;{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〉.

We study in this work the dynamics of a collection of identical hollow spheres (ping-pong balls) that rest on a horizontal metallic grid. Fluidization is achieved by means of a turbulent air current coming from below. The upflow is adjusted so that the particles do not levitate over the grid, resulting in quasi-two-dimensional dynamics. We show that the behavior of diffusion and correlations in this system is particularly rich. Noticeably as well (and related to the complex dynamical behavior), a variety of phases appear, with important peculiarities with respect to analogous setups. We observe gas, liquid, glass, and hexagonal crystal phases. Most notably, we show that the melting of the hexagonal crystal occurs in coexistence with a liquid phase. This strikingly differs from the corresponding transition in a purely two-dimensional systems of air-fluidized disks, for which no phase coexistence has been reported in the literature.

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. Inhomogeneously growing domains are frequently encountered in biological phenomena involving stochastic transport, such as tumor growth and morphogen gradient formation. 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 favorably 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. From a more general perspective, our findings help to mitigate the scarcity of analytic results for reaction-diffusion problems in geometries displaying nonuniform growth. They are also expected to pave the way for further results, including the treatment of first-passage problems associated with encounter-controlled reactions in such domains.

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.

We have designed a two-dimensional, fractal-like lattice and explored, both numerically and analytically, the differences between random walks on this lattice and a regular, square-planar Euclidean lattice. We study the efficiency of diffusion-controlled processes for flows from external sites to a centrosymmetric reaction center and, conversely, for flows from a centrosymmetric source to boundary sites. In both cases, we find that analytic expressions derived for the mean walk length on the fractal-like lattice have an algebraic dependence on system size, whereas for regular Euclidean lattices the dependence can be transcendental. These expressions are compared with those derived in the continuum limit using classical diffusion theory. Our analysis and the numerical results quantify the extent to which one paradigmatic class of spatial inhomogeneities can compromise the efficiency of adatom diffusion on solid supports and of surface-assisted self-assembly in metal-organic materials.

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.

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 r at time t given that it started its walk from r&#95;{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 consider various anomalous d -dimensional diffusion problems in the presence of an absorbing boundary with radial symmetry. The motion of particles is described by a fractional diffusion equation. Their mean-square displacement is given by r(2) proportional, variant t(gamma)(00 , the emergence of such series in the long-time domain is a specific feature of subdiffusion problems. We present a method to regularize such series, and, in some cases, validate the procedure by using alternative techniques (Laplace transform method and numerical simulations). In the normal diffusion case, we find that the signature of the initial condition on the approach to the steady state rapidly fades away and the solution approaches a single (the main) decay mode in the long-time regime. In remarkable contrast, long-time memory of the initial condition is present in the subdiffusive case as the spatial part Psi1(r) describing the long-time decay of the solution to the steady state is determined by a weighted superposition of all spatial modes characteristic of the normal diffusion problem, the weight being dependent on the initial condition. Interestingly, Psi1(r) turns out to be independent of the anomalous diffusion exponent gamma .

We investigate the long-time behavior of the survival probability P(t) of a mobile particle in d-dimensional continuous Euclidean media doped with noninteracting mobile traps. The particle is strictly subdiffusive, implying that its mean-square displacement grows as tgamma' with 0or=2 we find that at late times the survival probability is that of the pure target problem (the problem where the particle remains immobile) in agreement with previous studies for the d=1 case. These decay laws remain invariant over the whole gamma range as opposed to the dynamical crossover observed for the case of a purely diffusive particle (gamma'=1) where, for gamma<2/(2+d) , the survival probability becomes that of the so-called trapping problem (the problem where the particle moves in a sea of static traps). This behavior implies that for sufficiently low values of gamma(gamma<2/(2+d)) the survival probability becomes singular in the limit gamma'-->1: trappinglike for gamma'=1 and targetlike for any gamma'<1.

We use a Markov method to study the efficiency of trapping processes involving both a random walker and a deep trap in regular and disordered networks. The efficiency is gauged by the mean absorption time (average of the mean number of steps performed by the random walker before being absorbed by the trap). We compute this quantity in terms of different control parameters, namely, the length of the walker jumps, the mobility of the trap, and the degree of spatial disorder of the network. For a proper choice of the system size, we find in all cases a nonmonotonic behavior of the efficiency in terms of the corresponding control parameter. We thus arrive at the conclusion that, despite the decrease of the effective system size underlying the increase of the control parameter, the efficiency is reduced as a result of an increase of the escape probability of the walker once it finds itself in the interaction zone of the trap. This somewhat anti-intuitive effect is very robust in the sense that it is observed regardless of the specific choice of the control parameter. For the case of a ring lattice, results for the mean absorption time in systems of arbitrary size are given in terms of a two-parameter scaling function. For the case of a mobile trap, we deal with both trapping via a single channel (walker-trap overlap) and via two channels (walker-trap overlap and walker-trap crossing), thereby generalizing previous work. As for the disordered case, our analysis concerns small world networks, for which we see several crossovers of the absorption time as a function of the control parameter and the system size. The methodology used may be well suited to exploring characteristic time scales of encounter-controlled phenomena in networks with a few interacting elements and the effect of geometric constraints in nanoscale systems with a very small number of particles.

We study the efficiency of a two-channel reaction between two walkers on a finite one-dimensional periodic lattice. The walkers perform a combination of synchronous and asynchronous jumps on the lattice and react instantaneously when they meet at the same site (first channel) or upon position exchange (second channel). We develop a method based on a conditional first-passage problem to obtain exact results for the mean number of time steps needed for the reaction to take place as well as for higher order moments. Previous results obtained in the framework of a difference equation approach are fully confirmed, including the existence of a parity effect. For even lattices the maximum efficiency corresponds to a mixture of synchronous events and a small amount of asynchronous events, while for odd lattices the reaction time is minimized by a purely synchronous process. We provide an intuitive explanation for this behavior. In addition, we give explicit expressions for the variance of the reaction time. The latter displays a similar even-odd behavior, suggesting that the parity effect extends to higher order moments.

We study the behavior of the chemical reactions A+A-->A+S and A+A-->S+S (where the reactive species A and the inert species S are both assumed to be immobile) embedded on Bethe lattices of arbitrary coordination number z and on a two-dimensional (2D) square lattice. For the Bethe lattice case, exact solutions for the coverage in the A species in terms of the initial condition are obtained. In particular, our results hold for the important case of an infinite one-dimensional (1D) lattice (z=2). The method is based on an expansion in terms of conditional probabilities which exploits a Markovian property of these systems. Along the same lines, an approximate solution for the case of a 2D square lattice is developed. The effect of dilution in a random initial condition is discussed in detail, both for the lattice coverage and for the spatial distribution of reactants.

We present a stochastic, time-discrete Boolean model that mimics the mesoscopic dynamics of the desorption reactions A+A-->A+S and A+A-->S+S in a one-dimensional lattice. In the continuous-time limit, we derive a hierarchy of dynamical equations for the subset of moments involving contiguous lattice sites. The solution of the hierarchy allows to compute the exact dynamics of the mean coverage for both microscopic and coarse-grained initial conditions, which turn out to be different from the mean field predictions. The evolution equations for the mean coverage and the second-order moments are shown to be equivalent to those provided by a time-continuous master equation. The important role of higher-order fluctuations is brought out by the failure of a truncation scheme retaining only two-particle fluctuation correlations.