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 \rho(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 \rho(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

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 $\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_{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.Comment: 18 pages, 7 figures. The title has been changed with respect to the one in the previous versio

Harvesting Thermal Fluctuations: Activation Process Induced by a Nonlinear Chain in Thermal Equilibrium

We present a model in which the immediate environment of a bistable system is a molecular chain which in turn is connected to a thermal environment of the Langevin form. The molecular chain consists of masses connected by harmonic or by anharmonic springs. The distribution, intensity, and mobility of thermal fluctuations in these chains is strongly dependent on the nature of the springs and leads to different transition dynamics for the activated process. Thus, all else (temperature, damping, coupling parameters between the chain and the bistable system) being the same, the hard chain may provide an environment described as diffusion-limited and more effective in the activation process, while the soft chain may provide an environment described as energy-limited and less effective. The importance of a detailed understanding of the thermal environment toward the understanding of the activation process itself is thus highlighted

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.Comment: The previous version was corrupted - some figures misplaced, some strange words that did not belong. Otherwise identica

Bone health in patients with multiple sclerosis relapses

OBJECTIVES: To evaluate the bone health and vitamin D levels of a cohort of patients with relapses of multiple sclerosis (MS) and to propose an algorithm for the management of bone health in this patient group. METHODS: We prospectively studied 56 consecutive patients from our acute relapse clinic. 3 patients were excluded from analysis as they were not deemed to have experienced an acute MS relapse. Bone health was assessed with vitamin D levels and Dual Energy X-ray Absorptiometry (DEXA) scanning (10 patients failed to attend for DEXA). Statistical analyses were used to compare groups and identify predictive variables. A review of the literature led to a proposed management protocol. RESULTS: Pre-relapse the baseline EDSS was ≤6.5 in all subjects, and <4.0 in the majority (66%). Most received corticosteroids. 51% had low bone mineral density (BMD) as defined by a T-score less than −1.0 on DEXA scanning. Three were osteoporotic (T-score less than −2.5). Thirty one of fifty (62%) subjects were Vitamin D deficient (25(OH)D less than 50 nmol/L). A range of variables, including previous corticosteroid usage, were not significantly predictive of reduced BMD. CONCLUSIONS: There was a high frequency of both low BMD and Vitamin D deficiency in this cohort of relatively young and largely ambulatory patients experiencing MS relapses. Current tools, such as the WHO FRAX algorithm, are inadequate in assessing bone status and fracture risk in this patient group, predominantly as they are focused on older age groups. We propose a simple clinical management algorithm

Automated search for galactic star clusters in large multiband surveys: I. Discovery of 15 new open clusters in the Galactic anticenter region

Aims: According to some estimations, there are as many as 100000 open clusters in the Galaxy, but less than 2000 of them have been discovered, measured, and cataloged. We plan to undertake data mining of multiwavelength surveys to find new star clusters. Methods: We have developed a new method to search automatically for star clusters in very large stellar catalogs, which is based on convolution with density functions. We have applied this method to a subset of the Two Micron All Sky Survey catalog toward the Galactic anticenter. We also developed a method to verify whether detected stellar groups are real star clusters, which tests whether the stars that form the spatial density peak also fall onto a single isochrone in the color-magnitude diagram. By fitting an isochrone to the data, we estimate at the same time the main physical parameters of a cluster: age, distance, color excess. Results: For the present paper, we carried out a detailed analysis of 88 overdensity peaks detected in a field of $16\times16$ degrees near the Galactic anticenter. From this analysis, 15 overdensities were confirmed to be new open clusters and the physical and structural parameters were determined for 12 of them; 10 of them were previously suspected to be open clusters by Kronberger (2006) and Froebrich (2007). The properties were also determined for 13 yet-unstudied known open clusters, thus almost tripling the sample of open clusters with studied parameters in the anticenter. The parameters determined with this method showed a good agreement with published data for a set of well-known clusters.Comment: accepted to A&