skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Relativistic diffusion processes and random walk models

Abstract

The nonrelativistic standard model for a continuous, one-parameter diffusion process in position space is the Wiener process. As is well known, the Gaussian transition probability density function (PDF) of this process is in conflict with special relativity, as it permits particles to propagate faster than the speed of light. A frequently considered alternative is provided by the telegraph equation, whose solutions avoid superluminal propagation speeds but suffer from singular (noncontinuous) diffusion fronts on the light cone, which are unlikely to exist for massive particles. It is therefore advisable to explore other alternatives as well. In this paper, a generalized Wiener process is proposed that is continuous, avoids superluminal propagation, and reduces to the standard Wiener process in the nonrelativistic limit. The corresponding relativistic diffusion propagator is obtained directly from the nonrelativistic Wiener propagator, by rewriting the latter in terms of an integral over actions. The resulting relativistic process is non-Markovian, in accordance with the known fact that nontrivial continuous, relativistic Markov processes in position space cannot exist. Hence, the proposed process defines a consistent relativistic diffusion model for massive particles and provides a viable alternative to the solutions of the telegraph equation.

Authors:
; ;  [1]
  1. Institut fuer Physik, Universitaet Augsburg, Theoretische Physik I, Universitaetstrasse 1, D-86135 Augsburg (Germany)
Publication Date:
OSTI Identifier:
21011033
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physical Review. D, Particles Fields; Journal Volume: 75; Journal Issue: 4; Other Information: DOI: 10.1103/PhysRevD.75.043001; (c) 2007 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; DIFFUSION; GRAPH THEORY; LIGHT CONE; MARKOV PROCESS; MATHEMATICAL SOLUTIONS; PROBABILITY; PROPAGATOR; RANDOMNESS; RELATIVISTIC RANGE; RELATIVITY THEORY; STANDARD MODEL; VISIBLE RADIATION

Citation Formats

Dunkel, Joern, Talkner, Peter, and Haenggi, Peter. Relativistic diffusion processes and random walk models. United States: N. p., 2007. Web. doi:10.1103/PHYSREVD.75.043001.
Dunkel, Joern, Talkner, Peter, & Haenggi, Peter. Relativistic diffusion processes and random walk models. United States. doi:10.1103/PHYSREVD.75.043001.
Dunkel, Joern, Talkner, Peter, and Haenggi, Peter. Thu . "Relativistic diffusion processes and random walk models". United States. doi:10.1103/PHYSREVD.75.043001.
@article{osti_21011033,
title = {Relativistic diffusion processes and random walk models},
author = {Dunkel, Joern and Talkner, Peter and Haenggi, Peter},
abstractNote = {The nonrelativistic standard model for a continuous, one-parameter diffusion process in position space is the Wiener process. As is well known, the Gaussian transition probability density function (PDF) of this process is in conflict with special relativity, as it permits particles to propagate faster than the speed of light. A frequently considered alternative is provided by the telegraph equation, whose solutions avoid superluminal propagation speeds but suffer from singular (noncontinuous) diffusion fronts on the light cone, which are unlikely to exist for massive particles. It is therefore advisable to explore other alternatives as well. In this paper, a generalized Wiener process is proposed that is continuous, avoids superluminal propagation, and reduces to the standard Wiener process in the nonrelativistic limit. The corresponding relativistic diffusion propagator is obtained directly from the nonrelativistic Wiener propagator, by rewriting the latter in terms of an integral over actions. The resulting relativistic process is non-Markovian, in accordance with the known fact that nontrivial continuous, relativistic Markov processes in position space cannot exist. Hence, the proposed process defines a consistent relativistic diffusion model for massive particles and provides a viable alternative to the solutions of the telegraph equation.},
doi = {10.1103/PHYSREVD.75.043001},
journal = {Physical Review. D, Particles Fields},
number = 4,
volume = 75,
place = {United States},
year = {Thu Feb 15 00:00:00 EST 2007},
month = {Thu Feb 15 00:00:00 EST 2007}
}
  • A method is presented for calculating exactly the relative width (sigma/sup 2/)/sup 1/2//, the skewness ..gamma../sub 1/, and the kurtosis ..gamma../sub 2/ characterizing the probability distribution function for three random-walk models of diffusion-controlled processes. For processes in which a diffusing coreactant A reacts irreversibly with a target molecule B situated at a reaction center, three models are considered. The first is the traditional one of an unbiased, nearest-neighbor random walk on a d-dimensional periodic/confining lattice with traps; the second involves the consideration of unbiased, non-nearest-neigh bor (i.e., variable-step length) walks on the same d-dimensional lattice; and, the third deals withmore » the case of a biased, nearest-neighbor walk on a d-dimensional lattice (wherein a walker experiences a potential centered at the deep trap site of the lattice). Our method, which has been described in detail elsewhere (P.A. Politowicz and J. J. Kozak, Phys. Rev. B 28, 5549 (1983)) is based on the use of group theoretic arguments within the framework of the theory of finite Markov processes.« less
  • We present computer-simulation results for the trapping rate (rate constant) /ital k/ associated with diffusion-controlled reactions among identical, static spherical traps distributed with an arbitrary degree of impenetrability using a Pearson random-walk algorithm. We specifically consider the penetrable-concentric-shell model in which each trap of diameter sigma is composed of a mutually impenetrable core of diameter lambdasigma, encompassed by a perfectly penetrable shell of thickness (1/minus/lambda)sigma/2: lambda=0 corresponding to randomly centered or ''fully penetrable'' traps and lambda=1 corresponding to totally impenetrable traps. Trapping rates are calculated accurately from the random-walk algorithm at the extreme limits of lambda (lambda=0 and 1) andmore » at an intermediate value (lambda=0.8), for a wide range of trap densities. Our simulation procedure has a relatively fast execution time. It is found that /ital k/ increases with increasing impenetrability at fixed trap concentration. These ''exact'' data are compared with previous theories for the trapping rate. Although a good approximate theory exists for the fully-penetrable-trap case, there are no currently available theories that can provide good estimates of the trapping rate for a moderate to high density of traps with nonzero hard cores (lambda/gt/0).« less
  • The motion of cells and molecules through biological environments is often hindered by the presence of other cells and molecules. A common approach to modeling this kind of hindered transport is to examine the mean squared displacement (MSD) of a motile tracer particle in a lattice-based stochastic random walk in which some lattice sites are occupied by obstacles. Unfortunately, stochastic models can be computationally expensive to analyze because we must average over a large ensemble of identically prepared realizations to obtain meaningful results. To overcome this limitation we describe an exact method for analyzing a lattice-based model of the motionmore » of an agent moving through a crowded environment. Using our approach we calculate the exact MSD of the motile agent. Our analysis confirms the existence of a transition period where, at first, the MSD does not follow a power law with time. However, after a sufficiently long period of time, the MSD increases in proportion to time. This latter phase corresponds to Fickian diffusion with a reduced diffusivity owing to the presence of the obstacles. Our main result is to provide a mathematically motivated, reproducible, and objective estimate of the amount of time required for the transport to become Fickian. Our new method to calculate this crossover time does not rely on stochastic simulations.« less
  • The transport of magnetic field lines is studied numerically where three dimensional (3-D) magnetic fluctuations, with a power law spectrum, and periodic over the simulation box are superimposed on an average uniform magnetic field. The weak and the strong turbulence regime, {delta}{ital B}{similar_to}{ital B}{sub 0}, are investigated. In the weak turbulence case, magnetic flux tubes are separated from each other by percolating layers in which field lines undergo a chaotic motion. In this regime the field lines may exhibit Levy, rather than Gaussian, random walk, changing from Levy flights to trapped motion. The anomalous diffusion laws {l_angle}{Delta}{ital x}{sup 2}{sub {italmore » i}}{r_angle}{proportional_to}{ital s}{sup {alpha}} with {alpha}{gt}1 and {alpha}{lt}1, are obtained for a number of cases, and the non-Gaussian character of the field line random walk is pointed out by computing the kurtosis. Increasing the fluctuation level, and, therefore stochasticity, normal diffusion ({alpha}{congruent}1) is recovered and the kurtoses reach their Gaussian value. However, the numerical results show that neither the quasi-linear theory nor the two dimensional percolation theory can be safely extrapolated to the considered 3-D strong turbulence regime. {copyright} {ital 1995} {ital American} {ital Institute} {ital of} {ital Physics}.« less