Relativistic Brownian Motion and Diffusion Processes
Relativistic Brownian Motion and Diffusion Processes
Relativistic Brownian Motion and Diffusion Processes
Dissertation
Mathematisch-Naturwissenschaftliche Fakultät
Universität Augsburg
eingereicht von
Jörn Dunkel
Symbols 3
1
2 CONTENTS
Appendices
B Normalization constants 99
B.1 Jüttner function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
B.2 Diffusion propagator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
Bibliography 115
Lebenslauf 151
Danksagung 152
Symbols
3
4 CONTENTS
D noise amplitude
D spatial diffusion constant
B(s) d-dimensional standard Wiener process with time parameter s
P probability measure of the Wiener process
∗ Ito (pre-point) interpretation of the stochastic integral
◦ Stratonovich-Fisk (mid-point) interpretation of the stochastic integral
• backward Ito (post point) interpretation of the stochastic integral
N set of natural numbers 1, 2, . . .
Z set of integer numbers
R set of real numbers
λ Lebesgue measure
µ, ρ measures
hXi expected value of a random variable X
Chapter 1
In his annus mirabilis 1905 Albert Einstein published four manuscripts [1–4] that would
forever change the world of physics. Two of those papers [2, 3] laid the foundations for
the special theory of relativity, while another one [4] solved the longstanding problem of
classical (nonrelativistic) Brownian motion.1 Barring gravitational effects [5, 6], special
relativity has proven to be the correct framework for describing physical processes on
all terrestrial scales [7, 8]. Accordingly, during the past century extensive efforts have
been made to adapt established nonrelativistic theories such as, e.g., thermodynamics,
quantum mechanics or field theories [9] to the requirements of special relativity. Following
this tradition, the present thesis investigates how stochastic concepts such as Brownian
motion may be generalized within the framework of special relativity. The subsequent
chapters intend to provide a cohesive summary of results obtained during the past three
years [10–17], also taking into account important recent contributions by other authors (see,
e.g., [18–24]).
Historically, the term ‘Brownian motion’ refers to the irregular dynamics exhibited by
a test particle (e.g., dust or pollen) in a liquid environment. This phenomenon, already
mentioned by Ingen-Housz [25, 26] in 1784, was first analyzed in detail by the Scottish
botanist Robert Brown [27] in 1827. About 80 years later, Einstein [4], Sutherland [28]
and von Smoluchowski [29] were able to theoretically explain these observations. They
proposed that Brownian motion is caused by quasi-random, microscopic interactions with
molecules forming the liquid. In 1909 their theory was confirmed experimentally by Per-
rin [30], providing additional evidence for the atomistic structure of matter. During the
first half of the 20th century the probabilistic description of Brownian motion processes
was further elaborated in seminal papers by Langevin [31, 32], Fokker [33], Planck [34],
Klein [35], Uhlenbeck and Ornstein [36] and Kramers [37]. Excellent reviews of these early
contributions are given by Chandrasekhar [38] and Wang and Uhlenbeck [39].
1
Einstein’s first paper [1] provided the theoretical explanation for the photoelectric effect.
5
6 CHAPTER 1. INTRODUCTION AND HISTORICAL OVERVIEW
In parallel with the studies in the field of physics, outstanding mathematicians like Bache-
lier [40], Wiener [41–43], Kolmogoroff [44–46], Feller [47], and Lévy [48, 49] provided a
rigorous basis for the theory of Brownian motions and stochastic processes in general.
Between 1944 and 1968 their groundbreaking work was complemented by Ito [50,51], Gih-
man [52–54], Fisk [55, 56] and Stratonovich [57–59], who introduced and characterized dif-
ferent types of stochastic integrals or, equivalently, stochastic differential equations (SDEs).
The theoretical analysis of random processes was further developed over the past decades,
and the most essential results are discussed in several excellent textbook references [60–66]2 .
The modern theory of stochastic processes goes far beyond the original problem considered
by Einstein and his contemporaries, and the applications cover a wide range of different
areas including physics [67–74], biology [75, 76], economy and finance [77–79].
The present thesis is dedicated to the question how SDE-based Brownian motion models
can be generalized within the framework of special relativity. In the physics literature [65],
SDEs are often referred to as Langevin equations [31, 32], and we shall use both terms
synonymously here. From a mathematical point of view, SDEs [64] determine well-defined
models of stochastic processes; from a physicist’s point of view, their usefulness for the
description of a real system is a priori an open issue. Therefore, the derivation of nonrel-
ativistic Langevin equations from microscopic models has attracted considerable interest
over the past 60 years [13,80–86]. Efforts in this direction helped to clarify the applicability
of SDEs to physical problems and led, among others, to the concept of quantum Brownian
motion [82, 87–99].3
If one aims at generalizing the classical Brownian motion concepts to special relativity,
then several elements from relativistic equilibrium thermodynamics and relativistic statis-
tical mechanics play an important role. The first papers on relativistic thermodynamics
were written by Einstein [109] and Planck [110, 111] in 1907. A main objective of their
studies was to clarify the Lorentz transformation laws of thermodynamic variables (tem-
perature, pressure, etc.).4 In 1963 the results of Einstein and Planck were questioned by
Ott [115], whose work initiated an intense debate about the correct relativistic transfor-
mation behavior of thermodynamic quantities [116–160].5 However, as clarified by van
Kampen [137] and Yuen [161], the controversy surrounding relativistic thermodynamics
can be resolved by realizing that thermodynamic quantities can be defined in different,
2
The history of the mathematical literature on Brownian motions and stochastic processes is discussed
extensively in Section 2.11 of Ref. [64]; see also Chapters 2-4 in Nelson [61].
3
The vast literature on classical Brownian motion processes and their various applications in nonrel-
ativistic physics is discussed in several survey articles [68–73, 100–104]. Nonrelativistic generalizations of
the standard theory as, e.g., anomalous diffusion processes have been summarized in [73, 105, 106], while
review articles on nonrelativistic quantum Brownian motion can be found in [97, 98, 107, 108].
4
See also Pauli [112], Tolman [113] and van Dantzig [114] for early discussions of this problem.
5
The pre-1970 literature on this disputed issue has been reviewed by Yuen [161] and Ter Haar and
Wegland [162]; more recent surveys can be found in [163–165].
7
Arens and Babbitt [233], and various semi-relativistic approximations have been discussed,
e.g., in [234–236].10 Another, intensely studied method for describing relativistic interac-
tions is based on the so-called constraint formalism [200, 211, 248–264]. The foundations
of this approach were worked out by Dirac [248] in 1949, who aimed at constructing a
consistent relativistic quantum theory. However, compared with the nonrelativistic case,
it seems fair to say that neither of the various formulations has led to a relativistic statis-
tical many-particle theory that is on the same rigorous, commonly accepted footing as its
nonrelativistic counterpart.11
In spite of the difficulties impeding a rigorous treatment of relativistic many-particle sys-
tems, considerable progress has been made over the past century in constructing an approx-
imate relativistic kinetic theory based on one-particle phase space probability density func-
tions (PDFs). Early pioneering work that paved the way for the relativistic generalization
of the nonrelativistic Boltzmann equation [267,268] was done by Eckart [269], Lichnérowicz
and Marrot [270], Kluitenberg et al. [271], Beliaev and Budker [272], Synge [169], and Is-
rael [273].12 Comprehensive introductions to relativistic kinetic theory can be found in the
textbooks by Stewart [288], de Groot et al. [289], and Cercignani and Kremer [290], or also
in the reviews by Ehlers [291] and Andréasson [292].13
From relativistic kinetic theory [289, 290] it is only a relatively small step to formulating
a theory of relativistic Brownian motion processes in terms of Fokker-Planck equations
(FPEs) and Langevin equations. While the relativistic Boltzmann equation is a nonlinear
partial integro-differential equation, FPEs are linear partial differential equations and,
therefore, can be solved more easily [63]. In this work, we will mostly focus on relativistic
stochastic processes that are characterized by linear evolution equations for their respective
one-particle (transition) PDFs. The research in this direction may be roughly divided into
four different areas, although, of course, there are substantial overlaps and intersections
between them:
a) Relativistic Fokker-Planck equations in phase space. Similar to the relativistic Boltz-
mann equation, relativistic FPEs can be used to model non-equilibrium and relax-
ation phenomena in relativistic many-particle systems. Generally, an FPE can be
derived from a Langevin equation or as an approximation to a more general linear
master equation governing the stochastic process [65, 293]. Yet another way of de-
10
Kerner [237] has edited a reprint collection covering large parts of the pre-1972 literature on relativistic
action-at-a-distance models, and more recent contributions can be found in [199, 200, 238–247].
11
For a more detailed discussion of relativistic many-particle theory, we refer to the insightful consider-
ations in the orginal papers of Van Dam and Wigner [224, 225] and Hakim [192–194, 265] as well as to the
recent review by Hakim and Sivak [266].
12
See also [138, 153, 196, 234, 274–287].
13
Although standard relativistic kinetic theory can be considered as well-established nowadays [289,290],
some authors questioned its validity in recent years and proposed modifications of the relativistic Boltz-
mann equation [206,215–217]. Recent numerical simulations [17,214] support the standard theory [289,290].
9
Rapp et al. [314, 337], who analyzed thermalization effects in quark-gluon plasmas,
and also by Dieckmann et al. [338], who studied the thermalization in ultrarelativistic
plasma beam collisions as common in astrophysical settings.
∂
̺(t, x) = D ∇2̺(t, x). (1.1)
∂t
Concluding this brief historical overview, we may summarize that the theory of relativistic
Brownian motion and diffusion processes has experienced considerable progress during the
past decade, with applications in various areas of high energy physics [308, 310–312, 314,
393, 394, 406] and astrophysics [318, 320, 328, 338]. From a general perspective, relativis-
tic stochastic processes provide a useful approach whenever one has to model the quasi-
random behavior of relativistic particles in a complex environment. Therefore, it may be
expected that relativistic Brownian motion and diffusion concepts will play an increasingly
important role in future investigations of, e.g., thermalization and relaxation processes in
astrophysics or high energy collision experiments. The present work aims to provide a
comprehensive overview of the theory of relativistic Brownian motions with a particular
emphasis on relativistic Langevin equations. For this purpose, the subsequent parts are
organized as follows. Chapter 2 summarizes the Langevin theory of nonrelativistic Brown-
ian motions in phase space. Chapter 3 discusses relevant aspects of relativistic equilibrium
thermostatistics. Relativistic Langevin equations and their associated FPEs, are consid-
ered in Chapter 4. Chapter 5 is dedicated to relativistic diffusion processes in Minkowski
space-time; as outlined above, such processes must necessarily be non-Markovian. The
thesis concludes with a summary of open questions in Chapter 6, which may serve as a
starting point for future investigations and extensions of the theory. In order to present
the most important ideas and concepts in a transparent way, the discussion in the main
text will focus mostly on the simplest case of one space dimension (1D). The generalization
to higher space dimensions is usually straightforward and the corresponding equations are
summarized in the Appendix D.
Chapter 2
In order to briefly introduce the underlying mathematical concepts, we first recall some
basic definitions and results from the Langevin theory [31, 32] of nonrelativistic Brownian
motions. The Langevin and Fokker-Planck equations discussed in this part will be useful
later on, because they represent the nonrelativistic limit case of the relativistic theory,
which will be developed in Sections 4 and 5. The condensed discussion of nonrelativistic
Brownian motion processes in Section 2.1 is primarily based on the papers of Uhlenbeck
and Ornstein [36], Wang and Uhlenbeck [39], and Klimontovich [101]. For further reading
about nonrelativistic stochastic processes and their numerous applications in physics and
mathematics, we refer to the review articles of Chandrasekhar [38], Fox [100], Hänggi
and Thomas [67], Bouchaud and Georges [105], Metzler and Klafter [106], Hänggi and
Marchesoni [71], Frey and Kroy [72], or the textbooks [64–66, 407].
The present chapter is structured as follows. We begin by discussing the linear Langevin
equation of the classical nonrelativistic Ornstein-Uhlenbeck process. Subsequently, nonlin-
ear generalizations of this process will be considered. In this context, we will address the
choice of discretization rules and generalized fluctuation-dissipation theorems. These issues
will become important again later on, when we discuss the Langevin theory of relativistic
Brownian motions in Chapter 4. The last part of this chapter focusses on the question how
stochastic differential equations (SDEs) can be derived from microscopic models. As typ-
ical examples, the well-known harmonic oscillator model [80–86] and a recently proposed
binary collision model [13] will be considered. In contrast to the oscillator model, the
collision model can be generalized to the framework of special relativity, and its relativistic
version will be discussed in Section 4.3.
13
14 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
Free Brownian motion The standard paradigm for a free nonrelativistic Brownian
motion process in the absence of external forces is the Ornstein-Uhlenbeck process. The
Ornstein-Uhlenbeck process is determined by the Langevin equations [31, 32, 36, 38, 39]
dX P
= , (2.1a)
dt M
dP
= −αP + (2D)1/2 ∗ ζ(t). (2.1b)
dt
The first term on the rhs. of Eq. (2.1b) is the linear friction force, where the constant
friction coefficient α > 0 represents an inverse relaxation time. The stochastic Langevin
force L(t) = (2D)1/2 ∗ ζ(t) models the fluctuations in the heat bath.2 In the case of the
Ornstein-Uhlenbeck process, the amplitude of these fluctuations is tuned by the constant
noise parameter D > 0, and the Gaussian white noise process ζ(t) is characterized by:
hζ(t)i = 0, (2.2a)
hζ(t) ζ(s)i = δ(t − s), (2.2b)
with all higher cumulants being zero. In Eqs. (2.2), h · i is understood as an average over
all possible realizations of the noise process ζ(s). We summarize the physical assumptions,
implicitly underlying Eqs. (2.1) and (2.2):
• The heat bath is spatially homogeneous and stationary; i.e., relaxation processes
within the heat bath occur on time scales much shorter than the relevant dynamical
time scales associated with the motion of the heavy Brownian particle.
1
By definition, the mean velocity of the heat bath particles vanishes in Σ.
2
Throughout, the symbol ‘∗’ is used to denote Ito’s stochastic integral definition. A precise specification
of the employed stochastic integral convention (i.e., discretization rule) becomes relevant, if one wishes to
consider a momentum dependent noise amplitude D(P ) and/or nonlinear transfomations of the momentum
process P (t). The most commonly used stochastic integral definitions and their implications are discussed
in App. C. However, for our present discussion in Section 2.1.1, it suffices to read the symbol ‘∗’ in
Eq. (2.1b) as an ordinary multiplication sign.
2.1. LANGEVIN AND FOKKER-PLANCK EQUATIONS 15
• Stochastic impacts between the Brownian particle and the constituents of the heat
bath occur virtually uncorrelated.
• On a macroscopic level, the interaction between Brownian particle and heat bath
is sufficiently well described by the constant viscous friction coefficient α and the
stochastic Langevin force L(t) = (2D)1/2 ∗ ζ(t).
• Eqs. (2.1) hold in the lab frame Σ, corresponding to the specific inertial system,
where the average velocity of the heat bath particles vanishes for all times t.
In Section 2.2 we shall review how stochastic dynamical equations similar to Eqs. (2.1) can
be derived and motivated by means of specific microscopic models.
In the mathematical literature [64, 66], SDEs like the Langevin Eq. (2.1b) are usually
written in the differential notation
Here, dX(t) := X(t + dt) − X(t) denotes the position increment, dP (t) := P (t + dt) − P (t)
the momentum increment; B(t) is a standardized 1D Brownian motion or, equivalently, a
standard Wiener process [41, 64, 66, 67], whose increments
i.e., the increments dB(t) are independent random numbers drawn from a normal distri-
bution with variance dt. The two different representations (2.1) and (2.3) of the Ornstein-
Uhlenbeck process may be connected by formally identifying
In the remainder, SDEs will primarily be written in the differential notation of Eq. (2.3),
which may also be viewed as a simple numerical integration scheme, cf., e.g., Ref. [79, 336]
3
This means that the joint probability density of an arbitrary collection of subsequent increments dB(ti )
is a product of the Gaussians P[dB(ti )]; see, e.g., [64, 66] for a precise mathematical definition.
16 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
and App. C. From Eq. (2.3d) and the independence of the increments at different times
s 6= t, it follows that
(
0, t 6= s
hdB(t)i = 0, hdB(t) dB(s)i = (2.5)
dt, t = s,
where now the expectation h · i is taken with respect to the probability measure of the
Wiener process B(t).
In order for Eqs. (2.3) to define a well-posed problem, they must be complemented by
initial conditions.4 Generally, one could consider either deterministic initial conditions by
fixing X(0) = x0 and P (0) = p0 or probabilistic initial conditions by specifying initial
distributions for X(0) and P (0). Then, the solutions of Eqs. (2.3) read explicitly
Z t
X(t) = X(0) + ds P (s)/M, (2.6a)
0
Z t
−αt 1/2 −αt
P (t) = P (0) e + (2D) e eαs ∗ dB(s). (2.6b)
0
In the remainder, we primarily refer to deterministic initial conditions, assuming that the
initial position X(0) = x0 and the initial momentum P (0) = p0 of the Brownian particle
are known exactly. Combining the solution (2.6) with Eq. (2.3d), one finds for the first
two moments of the momentum coordinate [36, 38]
P (0)
hX(t) − X(0)i = (1 − e−αt ), (2.8a)
αM
2
2D t P (0) −αt 2
[X(t) − X(0)]2
= + 1 − e +
(αM)2 αM
D −αt −2αt
−3 + 4e − e . (2.8b)
α3 M 2
The asymptotic spatial diffusion constant D∞ , not to be confused with the noise amplitude
D, is usually defined by
1
2D∞ := lim [X(t) − X(0)]2 . (2.9)
t→∞ t
4
Without loss of generality we fix the initial time t0 = 0.
2.1. LANGEVIN AND FOKKER-PLANCK EQUATIONS 17
From Eq. (2.8b) we find for the Ornstein-Uhlenbeck process the classical result
D∞ := D/(αM)2 . (2.10)
When studying SDEs of the type (2.3b), one is typically interested in the probability
f (t, x, p) dx dp
of finding the Brownian particle at time t in the infinitesimal phase space interval [x, x +
dx] × [p, p + dp]. The non-negative phase space PDF f (t, x, p) ≥ 0 of the Brownian particle
is normalized at all times, i.e.
Z
1 = dxdp f (t, x, p) , ∀ t > 0; (2.11)
where, here and below, unspecified integrals range over the full phase space, position space,
or momentum space, respectively. Given the phase space PDF f (t, x, p), the marginal
momentum PDF φ(t, p) and the marginal position PDF ̺(t, x) are defined by
Z
φ(t, p) = dx f (t, x, p), (2.12a)
Z
̺(t, x) = dp f (t, x, p). (2.12b)
Deterministic initial data X(0) = x0 and P (0) = p0 translate into the initial conditions
For the Ornstein-Uhlenbeck process from Eq. (2.3b), the Fokker-Planck equation (FPE)
governing the momentum PDF φ(t, p) reads [220]
∂φ ∂ ∂φ
= αpφ + D . (2.14)
∂t ∂p ∂p
Adopting the deterministic initial condition (2.13b), the time-dependent solution of
Eq. (2.14) is given by [36, 220]
1/2
α[p − p0 exp(−αt)]2
α
φ(t, p) = exp − . (2.15)
2πD[1 − exp(−2αt)] 2D[1 − exp(−2αt)]
In the limit t → ∞ this solution reduces to the stationary Gaussian distribution
αp2
α 1/2
φ∞ (p) = exp − . (2.16)
2πD 2D
18 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
For a given momentum distribution φ(t, p) of the Brownian particle, the corresponding
velocity PDF ψ(t, v) is defined by
dp
ψ(t, v) := φ(t, p(v)), (2.17)
dv
where p = Mv in the nonrelativistic case. Hence, by imposing the Einstein relation
D = αMkB T, (2.18)
the stationary momentum PDF (2.16) is seen to be equivalent to Maxwell’s velocity dis-
tribution
1/2
Mv 2
M
ψM (v) = exp − , (2.19)
2πkB T 2kB T
where T is the temperature of the heat bath and kB the Boltzmann constant. Moreover,
the asymptotic spatial diffusion constant from Eq. (2.10) takes the form
D∞ := kB T/(αM). (2.20)
The Einstein relation (2.18) represents the simplest example of a fluctuation-dissipation
relation (FDR) by linking the noise amplitude D and the friction coefficient α to the
temperature T of the heat bath. On the level of the Langevin description, this relation is
motivated by the plausible assumption that, after a certain relaxation time, the Brownian
particle will be in thermodynamic equilibrium with the surrounding bath. In Section 2.2
it will be discussed how generalized FDRs may arise from specific microscopic models for
the interaction between Brownian particle and heat bath. Before doing this, however, we
briefly address a few generalizations of the free Ornstein-Uhlenbeck process (2.3b).
For arbitrary time and position dependent force fields F(t, x) it is generally very difficult,
and in many cases even impossible, to find exact time-dependent solutions of the Fokker-
Planck equation (2.22). In the simpler case of a time-independent, conservative force field
F(t, x) ≡ F (x) with confining6 potential Φ(x), i.e.
∂
F (x) = − Φ(x), (2.23)
∂x
one can determine the stationary solution attained in the limit t → ∞. Imposing as above
the Einstein relation D = αMkB T, the stationary solution of Eq. (2.22) is given by the
Maxwell-Boltzmann distribution [220, 221]
2
−1 p
f (x, p) = Z exp −β + Φ(x) , β := (kB T)−1 , (2.24)
2M
Here, the symbol ‘•’ signals the post-point discretization interpretation [60] of the
SDE (2.25), which means that the coefficient function D(p) is evaluated at the post-
point P (t + dt). A stochastic force with nonlinearly momentum dependent noise amplitude
function D(p) as in Eq. (2.25) is usually referred to as ‘multiplicative’ noise, in contrast
6
Conventionally, a potential Φ(x) is called ‘confining’ if it increases sufficiently fast for |x| → ∞ so that
the phase space PDF f is normalizable.
20 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
to the ‘additive’ noise encountered in Eqs. (2.3) and (2.21). When considering SDEs that
contain multiplicative noise terms, the specification of the discretization rule is necessary
because of the fact that, for fixed functions α(p) and D(p), different discretization schemes
in general lead to nonequivalent stochastic processes; put differently, the values of the
stochastic integral P (t) defined by Eq. (2.25) depend on the choice of discretization rule.
This is the most essential difference compared with ordinary differential equations, whose
integral curves (i.e., solutions) are independent of the underlying discretization scheme
when taking the continuum limit dt → 0.
In Eq. (2.25) we opted for the post-point rule; in principle, other discretization rules can
be used as well [65, 418, 419]. The most prominent alternatives are Ito’s [50, 51] pre-point
discretization (∗), corresponding to computing function D(p) at P (t), and the mid-point
rule (◦) of Stratonovich [57–59] and Fisk [55, 56], where D is evaluated at the mean value
[P (t) + P (t + dt)]/2. From the mathematical point of view, the choice of the discretization
rule reduces to a matter of convenience due to the following fact: For each pair of sufficiently
smooth functions (α• (p), D(p)), one can determine a pair of functions (α◦|∗ (p), D(p)) which
describes exactly the same stochastic dynamics when combined with another discretization
rule ◦ and ∗, respectively. The corresponding conversion formulae are summarized in
App. C.
From the practical point of view, each of the three above mentioned discretization meth-
ods possesses its own merits and drawbacks: Ito’s pre-point rule (∗) is particularly conve-
nient for numerical simulations, but care is required when considering nonlinear transfor-
mations G(P ) of the momentum coordinate due to modifications of the differential calculus,
cf. App. C. By contrast, if one adopts the Stratonovich-Fisk mid-point rule (◦), then the
transformation rules from ordinary differential calculus carry over, but it becomes more
difficult to implement this mid-point rule in numerical simulations. The latter disadvan-
tage also applies to the post-point rule employed in Eq. (2.25). However, as we shall see
next, the post-point rule (•) leads to a particularly simple form of the FDR.
Adopting the post-point rule, the Fokker-Planck equation for the momentum PDF φ(t, p)
of the stochastic process defined by Eq. (2.25) reads
∂φ ∂ ∂φ
= α• (p)p φ + D(p) . (2.26)
∂t ∂p ∂p
where N is a normalization constant, and p∗ some arbitrary constant such that the integral
in the exponential exists.
As follows from the general form (2.27) of the stationary solution, one may generate ar-
bitrary momentum distributions (e.g., Maxwell, Bose, Fermi or power law distributions)
by choosing the friction and noise amplitude functions α• (p) and D(p) in a suitable man-
ner [12, 420]. To briefly illustrate this, consider some normalized target PDF φ∗ (p) ≥ 0.
We would like to fix the relation between α• and D such that the stationary solution φ∞ (p)
coincides with φ∗ (p). Equating φ∗ (p) with φ∞ (p) from Eq. (2.27), taking the logarithm
and differentiating with respect to p we find the condition
α• (p) d
p = − log φ∗ (p). (2.28)
D(p) dp
α• (p)
= (MkB T)−1 , (2.29)
D(p)
the stationary distribution φ∞ (p) reduces to the Maxwell distribution from Eq. (2.16). It
should be stressed, however, that the FDRs (2.28) and (2.29) do fix only one of the two
coefficients α• (p) and D(p). Put differently, one is still free to adapt, e.g., the function
α• (p) such that the stochastic process (2.25) exhibits the correct relaxation behavior. This
freedom is a main reason why the Langevin approach is successfully applicable to a wide
range of thermalization processes [101]. Physically reasonable expressions for α• (p) may
be obtained from kinetic theory [308, 310, 421–424] or microscopic Hamiltonian models
that take into account the interactions as well as the statistical properties of the heat
bath [13, 80–83, 86, 425]. Examples will be discussed in Section 2.2.
Langevin equations of the type (2.25) and nonlinear friction effects [426] have been
studied extensively in various contexts during the past decades (see, e.g., the review by
Klimontovich [101]). The applications cover a wide range of different areas including laser
physics [101, 102, 427], optical lattices [428, 429], plasma physics [430–433], high energy
physics [308, 310], biologically and chemically motivated population and reaction dynam-
ics [434], active Brownian motion models [424, 435–440], or theoretical and experimental
studies of excitation and transition phenomena in nonlinear systems [441–444].
However, with regard to our subsequent discussion of relativistic Brownian motions, it will
be most important to keep in mind that the nonlinear Langevin equation (2.25) provides
a tool for constructing Brownian motion processes with arbitrary stationary velocity and
momentum distributions [12, 420].
22 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
1
hζ(t) ζ(s)i = exp(−|t − s|/τn ), (2.30)
τn
where the parameter τn is the relaxation time of the driving noise ζ. The mathematical
analysis of processes driven by colored noise is considerably more complicated than in the
case of δ-correlated white noise; for a detailed discussion we refer to the review by Hänggi
and Jung [103].
b) Alternatively, one can try to motivate and derive Langevin equations from micro-
scopic models. If successful, this approach yields explicit expressions for the friction
and noise functions in terms of the microscopic model parameters.
The remainder of this section addresses the latter problem, which has attracted consider-
able interested over the past decades [13,80–86,294,326,423,425,452]. From the physicist’s
point of view, Langevin equations provide an approximate description of the ‘exact’ micro-
scopic dynamics. Hence, in order to derive SDEs like (2.1) or (2.25) from, e.g., microscopic
Hamiltonian mechanics one has to impose certain approximations. These approximations
determine the range of applicability of the Langevin approach. Generally, one can pursue
at least two different routes for deriving SDEs of the type (2.1) and (2.25) from more
precise models:
(i) Starting from a Boltzmann equation [267, 268, 290] or master equation [65] for the
one-particle probability density of the Brownian particle, one can try to reduce these
integro-differential equations to a Fokker-Planck equation by performing suitable approxi-
mations [294,308,310,326,422–424,452]. Once the Fokker-Planck equation has been found,
it is straightforward to write down a corresponding Langevin dynamics [65, 66]. The mi-
croscopic collision dynamics is then encoded in the scattering cross-sections appearing in
the collision integral of the Boltzmann equation [268, 290].
(ii) Alternatively, one may start from a microscopic (e.g., Hamiltonian) model describing
the interaction between Brownian particle and heat bath. After eliminating the heat bath
degrees of freedom from the equations of motion for the Brownian particle [13,80–86,425],
one obtains a generalized Langevin equation which may be reduced to the form (2.1)
and (2.25) in certain limit cases. As a byproduct, the FDRs arise naturally within this
approach upon assuming a probability distribution for the (initial) bath configuration. To
illustrate this by example, we next consider the oscillator model [80–86] and the elastic
binary collision model [13]. In contrast to the more frequently studied harmonic oscillator
model, the collision model from Section 2.2.2 can be extended to the relativistic case; cf.
discussion in Section 4.3.
Here M, X and P are the mass, position and momentum of the Brownian particle and Φ(x)
is an external potential field; xr and pr denote the position and momentum of a heat bath
24 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
particle with mass mr , oscillator frequency ωr and coupling constant cr . Equation (2.31)
yields the following Hamilton equations of motions:
X cr
M Ẋ = P, Ṗ = F (X) + cr xr − 2
X ; (2.32a)
r
mr ω r
where, for given initial values X(0), P (0), {xr (0), pr (0)}, the memory friction kernel ν(t−s)
and the Langevin noise force L(t) are given by [86]
1 X c2r
ν(t − s) := cos[ωr (t − s)], (2.33c)
M r mr ωr2
X cr
pr (0)
L(t) := cr xr (0) − 2
X(0) cos(ωr t) + sin(ωr t) . (2.33d)
r
m r ω r m r ω r
In order to be able to characterize the properties of the noise force L(t) by means of an FDR,
one still needs to impose a distribution for the initial conditions {xr (0), pr (0)} of the bath
variables. In principle, this initial distribution can be chosen arbitrarily. Of particular in-
terest in canonical thermostatistics are equilibrium distributions of the Maxwell-Boltzmann
type. In the case of the generalized Langevin equation (2.33), a plausible choice for the
initial bath distribution corresponds to the PDF
fb ({xr (0), pr (0)} | X(0) ) = Z−1 ×
X pr (0)2 mr ω 2 2
r cr
exp −β + xr (0) − X(0) ; (2.34)
r
2mr 2 mr ωr2
β = (kB T)−1 is the inverse thermal energy, T the temperature, and Z the normalization
constant. The initial position X(0) of the Brownian particle enters in Eq. (2.34) as an
independent parameter, i.e., averages with respect to fb ({xr (0), pr (0)} | X(0)) are condi-
tional on the initial Brownian particle position X(0). Averaging the stochastic force L(t)
with respect to fb from Eq. (2.34), one finds
hL(t)ib = 0, (2.35a)
hL(t)L(s)ib = MkB T ν(t − s). (2.35b)
2.2. MICROSCOPIC MODELS 25
Equation (2.35b) represents the FDR for the generalized Langevin equation (2.33) given
the initial bath distribution (2.34). The generalized Langevin equation (2.33) differs from
Eqs. (2.1), (2.21) and (2.25) through the memory friction ν(t − s). The SDE (2.21),
which describes the Ornstein-Uhlenbeck process in an external force field, is recovered
from Eqs. (2.33) in the limit case9
where α is a constant friction coefficient. The limit case (2.36) can be illustrated by
rewriting the friction kernel (2.33c) in the more general form
Z ∞
ν(t − s) = dω C(ω) cos[ω(t − s)]. (2.37)
0
1 X c2r
C(ω) = δ(ω − ωr ), (2.38)
M r mr ωr2
one recovers the memory friction (2.33c) as a special case of Eq. (2.37). In order to obtain
the limit case (2.36) from Eq. (2.37), one can use the cosine-decomposition of the Dirac
δ-function
Z ∞
1
δ(t − s) = dω eiω(t−s)
2π −∞
1 ∞
Z
= dω cos[ω(t − s)]. (2.39)
π 0
Hence, upon comparing Eqs. (2.39) and (2.37), the white noise limit (2.36) corresponds to
the choice
The harmonic oscillator model10 provides a useful microscopic justification for the
Langevin equations (2.1) and (2.21) of the Ornstein-Uhlenbeck process. Unfortunately,
this model cannot be transferred to special relativity, as it is based on instantaneous har-
monic interactions-at-distance which violate fundamental relativistic principles. Therefore,
in the last part of this chapter we shall consider another microscopic model which is based
on strictly localized elastic binary collisions and, thus, can be extended to special relativity.
9
Rt
The prefactor 2 is required in Eq. (2.36) because of the convention 0 ds δ(t − s)P (s) = P (t)/2.
10
The quantum mechanical generalization of the harmonic oscillator model represents a paradigm for
quantum Brownian motions and has been studied, e.g., in [82, 87–93, 99]; see also the reviews by Grabert
et al. [107] and Hänggi and Ingold [98].
26 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
Collision kinematics We consider the elastic collision of the Brownian particle (mo-
mentum P , kinetic energy E) with a heat bath particle (momentum pr , kinetic energy ǫr ).
The collision process is governed by the energy-momentum conservation laws
E + ǫr = Ê + ǫ̂r , (2.41a)
P + pr = P̂ + p̂r , (2.41b)
where hat-symbols refer to the state after the collision. In the nonrelativistic case, we have,
e.g., before the collision
with V and vr denoting the velocities. Taking into account the kinematic conservation laws
(2.41), we find that the momentum gain ∆Pr of the Brownian particle per single collision
is given by
2m 2M
∆Pr := P̂ − P = − P+ pr . (2.43)
M +m M +m
In order to construct a Langevin-like equation from Eqs. (2.41)–(2.43), we consider the
total momentum change δP (t) of the Brownian particle within the time interval [t, t + δt],
assuming that:
• collisions occurring within [t, t + δt] can be viewed as independent events;
holds true. In particular, δt is supposed to be so small that there occurs at most only
one collision between the Brownian particle and a specific heat bath particle r; on
the other hand, δt should still be large enough, so that the total number of collisions
within δt is larger than 1. These requirements can be fulfilled simultaneously only if
m ≪ M holds.
11
Similar approaches are known from unimolecular rate theory, see, e.g., Section V in [68]. In the context
of quantum Brownian motions, a quantum-mechanical version of the collision model was proposed and
studied by Pechukas [95], and Tsonchev and Pechukas [96].
2.2. MICROSCOPIC MODELS 27
where Ir (t, δt) ∈ {0, 1} is the indicator function for a collision with the heat bath particle r
during the interval [t, t + δt]; i.e.
(
1 if a collision has occurred,
Ir (t, δt) = (2.45)
0 otherwise.
Evidently, the collision indicators depend on the position and velocity coordinates of the
collision partners. In the 1D case, Ir (t, δt) can be expressed in the form12
where X = X(t), xr = xr (t) are the ‘initial’ positions of the colliding particles at time t,
and
X ′ = X + V δt, x′r = xr + vr δt
their projected positions at time t + δt. The collision indicator from Eq. (2.46) is charac-
terized by
Ir (t, 0) = 0, (2.47a)
j
[Ir (t, δt)] = Ir (t, δt), j = 1, 2, . . . ; (2.47b)
δt
Ir (t, δt) ≈ |vr − V | δ(xr − X). (2.48a)
2
Combining Eqs. (2.43), (2.44) and (2.48a) yields
" N # N
X m X M
δP (t) ≈ −2 Ir (t, δt) P (t) + 2 pr Ir (t, δt), (2.48b)
r=1
M +m r=1
M + m
12
The Heaviside-function Θ(x) is defined as the integral over the Dirac δ-function, i.e., Θ(x) := 0, x < 0;
Θ(0) := 1/2; Θ(x) := 1, x > 0. When considering higher space dimensions, the expression (2.46) for the
indicator function has to be modified accordingly, e.g., by taking into account the geometric shape of the
Brownian particle.
28 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
where, additionally, it was assumed that for each collision occurring within [t, t + δt], the
momentum of the Brownian particle before the collision is approximately equal to the
‘initial’ value P (t). In view of m ≪ M, Eq. (2.48b) can be simplified further to give
" N # N
Xm X
δP (t) ≈ −2 Ir (t, δt) P (t) + 2 pr Ir (t, δt). (2.48c)
r=1
M r=1
A comparison with the Langevin equations (2.3b) and (2.25) suggests that, heuristically,
the first term on the rhs. of Eq. (2.48c) can be interpreted as a ‘friction’ term, while
the second term represents ‘noise’. On the other hand, although looking quite similar
to a Langevin equation, Eq. (2.48c) is still considerably more complicated than, e.g., the
nonlinear Langevin equation (2.25). This is due to the fact that the collision indicators
Ir (t, δt) from Eq. (2.48a) depend not only on the Brownian particle’s position and velocity
but also on the stochastic bath variables {xr , vr }. Nevertheless, it is possible to calculate
the statistical properties of the momentum increments δP (t) from Eqs. (2.48), provided
one specifies a distribution for the heat bath particles.
Bath distribution and drift In principle, one can use Eqs. (2.48) to calculate the
statistical moments h(δP )j ib for an arbitrarily given heat bath PDF fbN ({xr , pr }). Here,
we shall focus on the situation where the (infinitely large) heat bath is given by a quasi-ideal
gas which is in thermal equilibrium with its environment. In this case, the one-particle
PDF fb1 (xr , pr ) is given by the spatially homogeneous Maxwell distribution
p2r
1 −1/2 −1
fb (xr , pr ) = (2πmkB T) L exp − , (2.49)
2mkB T
where xr ∈ [0, L] with L being the 1D container volume. Moreover, we will assume that:
• the collisions with the Brownian particle do not significantly alter the bath distribu-
tion, so that the total energy of the bath particles remains constant.
The above assumptions can be justified for a sufficiently large bath, if collisions between
the bath particles rapidly reestablish a spatially homogeneous bath distribution.
In order to calculate the moments h(δP )j ib , we note that, for a spatially uniform bath
distribution as in Eq. (2.49), the one-particle expectation value hG(xr , vr ) [Ir (t, δt)]j ib is
given by
(2.47a)
G(xr , vr ) [Ir (t, δt)]j b
= hG(xr , vr ) Ir (t, δt)ib
δt ∞
Z
(2.48a)
≈ dvr G(X, vr ) |vr − V | ψb (vr ), (2.50)
2L −∞
2.2. MICROSCOPIC MODELS 29
with ψb (vr ) denoting the one-particle velocity PDF of the heat bath particles. For the
canonical heat bath distribution from Eq. (2.49), the velocity PDF ψb (v) corresponds to
the Maxwellian
−1/2
vB2 π exp −vr2 /vB2 , vB := (2kB T/m)1/2 .
ψb (vr ) = (2.51)
In particular, we obtain
2
nb vB −1/2 P P P
N hIr (t, δt)ib = π exp − + erf δt,
2 pB pB pB
nb vB m P
N hpr Ir (t, δt)ib = − pB erf δt, (2.52)
4 M pB
2
nb vB 2 m 2 P 1 P P
N p2r Ir (t, δt) b
= pB π −1/2
exp − + erf δt,
2 M pB 2 pB pB
where nb = N/L is the number density of the heat bath particles, pB := MvB =
M(2kB T/m)1/2 , and the error function erf(z) is defined by
Z z
2 2
erf(z) := √ dx e−x .
π 0
By making use of Eqs. (2.52), we find the mean drift of the collision model:
N
mX N
X
hδP (t)ib = −2 hIr (t, δt)ib P + 2
hpr Ir (t, δt)ib
M r=1 r=1
m
= −2N hIr (t, δt)ib P + 2N hpr Ir (t, δt)ib
M
2
(2.52)
−1/2 P P P P
≈ −2nb kB T π exp − + erf δt −
pB pB pB pB
P
nb kB T erf δt. (2.53)
pB
Figure 2.1 depicts the mean drift force hδP (t)/δtib , obtained from Eq. (2.53). The absolute
value of this drift force grows linearly for small momentum values (Ornstein-Uhlenbeck
regime) and quadratically for large momentum values.
We still consider the second moment:
N X
m 2 X N N X
mX N
2 2
[δP (t)] b
= 4 hIr Is ib P − 8 hpr Ir Is ib P +
M s=1 r=1
M s=1 r=1
N X
X N
4 hpr ps Ir Is ib .
r=1 s=1
30 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
20
hδP/δtib / (nb kB T)
10
−10
−20
−4 −2 0 2 4
P / pB
Figure 2.1: Nonrelativistic binary collision model. Mean drift force hδP (t)/δtib from Eq. (2.53)
with nb = N/L denoting the number density of the heat bath particles and pB := M (2kB T/m)1/2
the characteristic momentum of a Brownian particle (mass M ), surrounded by heat bath particles
of mass m and temperature T.
Neglecting contributions of order δt2 , only the terms with s = r remain, and we obtain
m 2 m
[δP (t)]2 b ≈ 4N hIr ib P 2 − 8N hpr Ir ib P + 4N p2r Ir b ,
M M
where the averages are determined by Eqs. (2.52). Evidently, higher moments
j
h[δP (t)] ib , j ≥ 3 can be estimated in the same manner.
In the remainder of this section, however, we are going to discuss a systematic proce-
dure for approximating the model equations (2.48) by a nonlinear SDE of the Langevin-
type (2.25).
Langevin approximation As stated before, the Eq. (2.48c) for the momentum incre-
ments δP (t) in the binary collision model is not yet a Langevin equation. Therefore, we
conclude this section by discussing how one can approximate Eqs. (2.48) by a nonlinear
SDE of the form [cf. Eq (2.25)]
Langevin equations of the type (2.54) are phenomenological model equations that provide
a simplified description of the ‘exact’ microscopic dynamics. In order to obtain a useful
Langevin model for a given microscopic dynamics, the coefficient functions α and D in
Eq. (2.54) have to be chosen such that they yield the best possible approximation within
this class of SDEs. Here, we define the ‘best approximation’ by means of the following two
criteria: The stochastic process described by Eq. (2.54) should
2.2. MICROSCOPIC MODELS 31
• approach the correct stationary momentum distribution for the Brownian particle;
According to the discussion in Section 2.1.2, this implies that α and D must be coupled
by the Einstein condition
The rhs. of Eq. (2.57) may be determined from Eq. (2.53), yielding the mean drift force
δP (t)
g(p) := P (t) = p
δt b
2
−1/2 p p
= −2 nb kB T π exp − +
pB pB
" #
2
p 1 p
+ erf . (2.58a)
pB 2 pB
In order to evaluate the lhs. of Eq. (2.57), we note that for the post-point discretization
rule (•) it is known that [cf. (C.25)]
where the prime denotes the derivative with respect to the momentum variable. Substi-
tuting the Einstein relation (2.56), i.e., D(p) = α(P )MkB T, we obtain
Combining this with the conditional expectation for the friction term in Eq. (2.54),
Hence, by virtue of Eqs. (2.58), we see that the drift criterion (2.57) is equivalent to the
following ordinary differential equation (ODE) for α(p):
With respect to the two criteria formulated above, the solution of this ODE yields the
friction function α that provides the ‘best’ Langevin approximation to the binary collision
model. Information about the collision model and the bath distribution is encoded in the
mean drift force g(p). Evidently, the procedure leading to Eq. (2.59) can be generalized to
other interaction models/bath distributions as well – provided the stationary distribution
of the Brownian particle is known. Other types of interactions (e.g., nonelastic) would
result in another function g(p). A non-Maxwellian bath distribution would affect not only
the rhs. of Eq. (2.59) but also its lhs. due to a modified fluctuation-dissipation relation.
In order to analyze Eq. (2.59) for our specific collision model, it is useful to define rescaled
dimensionless quantities
The general solution of Eq. (2.61) with initial value α̃∗ (0) reads
π 1/2 1/2
p2∗ /µ∗ −1/2 −1/2
α̃∗ (p∗ ) = |p∗ | + e α̃∗ (0) − µ∗ + µ∗ erf µ∗ |p∗ | . (2.62)
2
The asymptotic solution (2.62) implies that the friction coefficient α would grow or decrease
exponentially unless one chooses a particular initial condition. Hence, the physically correct
2.2. MICROSCOPIC MODELS 33
0
−4 −2 0 2 4
P / pB
Figure 2.2: Nonrelativistic binary collision model. Nonlinear friction coefficent α(P ) from
Eq. (2.64), with pB = M (2kB T/m)1/2 denoting the characteristic thermal momentum of the
Brownian particle (mass M ) and nb = N/L the number density of the heat bath particles
(mass m).
asymptotic behavior, which is given by α∗ (p∗ ) ≃ |p∗ | for p∗ → ±∞, suffices to single out
the specific initial condition that had to be used in the general solution of Eq. (2.59).
Unfortunately, it seems to be very difficult or perhaps even impossible to find the exact
analytical solution of the ODE (2.59).
Therefore, for practical purposes, one could use the simpler ad hoc approximation
N
2 Xm
α(P ) ≈ hIr (t, δt)ib , (2.63)
δt r=1 M
which reflects the earlier heuristic interpretation of Eq. (2.48) in terms of ‘friction’ and
‘noise’ contributions. One then finds
2
kB T −1/2 P P P
α(P ) ≈ 2 nb π exp − + erf . (2.64)
pB pB pB pB
Figure 2.2 depicts the nonlinear friction coefficient function α(P ) from Eq. (2.64). A
Langevin equation based on α(P ) from Eq. (2.64) and D(P ) = α(P )MkB T yields the
correct stationary momentum distribution (2.55), but does not exhibit exactly the same
mean relaxation behavior as Eq. (2.48). Nonetheless, even such an approximate Langevin
equation will provide a considerably more accurate description of the Brownian motion in
a gaseous heat bath than, e.g., the classical Ornstein-Uhlenbeck process, which assumed
constant friction and noise coefficients, cf. Eq. (2.3b). For instance, an Ornstein-Uhlenbeck
(or Stokes-like) approximation could be obtained by replacing α(P ) through its minimum
34 CHAPTER 2. NONRELATIVISTIC BROWNIAN MOTION
value
1/2
m 2kB T
α(0) = nb , (2.65a)
M πm
D(0) = α(0) M kB T. (2.65b)
Relativistic equilibrium
thermostatistics
The preceding discussion of the nonrelativistic Brownian motion theory has shown that
equilibrium thermostatistics plays an important role in constraining the relation between
friction and noise coefficients in Langevin equations. Therefore, the present chapter intends
to clarify several aspects of relativistic equilibrium thermostatistics. These considerations
will become relevant in Chapter 4, when we will have to specify the fluctuation-dissipation
relations for the relativistic Langevin theory.
Evidently, the nonrelativistic Brownian motion models from Chapter 2 are in conflict
with special relativity because they do not prevent particles from moving faster than the
speed of light c. For example, the stationary velocity distribution of the classical, e.g.,
Ornstein-Uhlenbeck process (2.1) is given by a Maxwell distribution that is non-zero for
velocities |v| > c. Hence, if one wants to construct a relativistic version of the Ornstein-
Uhlenbeck process then one has to know the relativistic generalization of the Maxwell
distribution first. The recent literature has seen considerable debate about the correct gen-
eralization of Maxwell’s velocity distribution in special relativity [14, 16, 203, 211, 215, 216].
In Sec. 3.2 we shall present results of fully relativistic 1D molecular dynamics simula-
tions [17] which clearly favor a distribution that was proposed by Jüttner [167] in 1911,
i.e., six years after Einstein had formulated his theory of special relativity [2,3]. Moreover,
as discussed in the last part of this chapter, our simulations can also be used to illustrate
the meaning of temperature and thermal equilibrium in special relativity [17].
35
36 CHAPTER 3. RELATIVISTIC EQUILIBRIUM THERMOSTATISTICS
3.1 Preliminaries
Section 3.1.1 contains a brief summary of the definitions and the notation that will be used
in the remainder. Subsequently, the ‘peculiar’ transformation behavior of one-particle
phase space PDFs under Lorentz transformations will be addressed. For simplicity, we
restrict ourselves to discussing the 1D case (i.e., one time dimension and one space dimen-
sion); the corresponding generalization to higher space dimensions is straightforward.
with M > 0 denoting the particle rest mass and Lorentz factor
−1/2
γ(v) := 1 − v 2 /c2 . (3.2)
To keep the subsequent formulas as simple as possible, we will from now on adopt a natural
unit system with c = 1 yielding, e.g., the simplified relations
Time and position, and energy and momentum can be combined into the contravariant
four-vectors1
Here, ηαβ denote the components of the Minkowski metric tensor, and Einstein’s summation
convention has been used, i.e.,
X
ηαβ X β := ηαβ X β .
β
The above definitions refer to the inertial lab frame Σ. Now consider a second inertial
frame Σ′ , moving with velocity w relative to Σ. According to Einstein’s theory of special
1
We shall use the term ‘four-vector’ regardless of the actual number of spatial dimensions. Contravariant
(covariant) four-vectors will be treated as column (row) vectors.
3.1. PRELIMINARIES 37
The inverse Λ−1 (w) is equal to Λ(−w). The transformation law of an arbitrary four-vector
(Aα ) reads A′ α = Λ(w)α β Aβ , leading in the case of (X α ) and (P α) to the explicit results
′ ′
t t − wX E E − wP
= γ(w) , = γ(w) . (3.7)
X′ −wt + X P′ −wE + P
The rest mass M of a particle is Lorentz invariant, i.e.,
α 2
M = (E 2 − P 2 )1/2 = (−Pα P α )1/2 = (−Pα′ P ′ )1/2 = (E ′2 − P ′ )1/2 = M ′ .
Throughout, we will assume that the rest mass of a particle is not changed by interactions.
Finally, considering a particle with velocity curve V (t), its Lorentz invariant proper time
span ∆τ , elapsing between coordinate times t1 and t2 , is given by
Z t2
∆τ := dt [1 − V (t)2 ]1/2 ⇔ dτ := dt [1 − V (t)2 ]1/2 . (3.8)
t1
N f (t, x, p) dxdp
particles in the phase space interval [x, x + dx] × [p, p + dp]. Alternatively, when con-
sidering the random motion of a single Brownian particle in a fluctuating medium, the
quantity f (t, x, p)dxdp gives the probability of finding the Brownian particle at lab time t
in [x, x + dx] × [p, p + dp]. In either case, the function f is subject to the t-simultaneous
normalization condition
Z
1 = dxdp f (t, x, p). (3.9)
38 CHAPTER 3. RELATIVISTIC EQUILIBRIUM THERMOSTATISTICS
Now consider a second observer O′ , moving with velocity w 6= 0 relative to Σ. The ob-
server O′ will measure another distribution f ′ (t′ , x′ , p′ ) and one is led to wonder how the two
functions f ′ (t′ , x′ , p′ ) and f (t, x, p) are related to each other. In the nonrelativistic theory,
the change from one inertial system to another does not affect the time coordinate; hence,
one can use the standard transformation laws for PDFs in that case [see, e.g., Eq. (2.17)].
By contrast, in the relativistic case the situation becomes more complicated, because now
the definition of f and f ′ is based on an observer-dependent notion of simultaneity. Put dif-
ferently, the measurements of O and O′ refer to the two different hyperplanes “t=constant”
and “t′ =constant” in Minkowski space, respectively. Van Kampen [138] was able to prove
that the one-particle phase space density f transforms as a Lorentz scalar, i.e.,2
where (t, x, p) and (t′ , x′ , p′ ) are connected by a Lorentz transformation with velocity param-
eter w. Moreover, he showed that the function f ′ satisfies the t′ -simultaneous normalization
condition
Z
1 = dx′ dp′ f (t′ , x′ , p′ ). (3.11)
Van Kampen’s proof [138] of Eq. (3.10) uses an assumption about the uniqueness of particle
trajectories and a reparameterization of the particles trajectories in terms of their invariant
proper times. As a consequence, Eq. (3.10) represents a generic kinematical result and
applies to a broad class of interaction models.
We next summarize several useful implications of Eq. (3.10). For this purpose, we define
the marginal densities
Z Z
φ(t, p) = dx f (t, x, p), φ (t , p ) = dx′ f ′ (t′ , x′ , p′ ),
′ ′ ′
(3.12a)
Z Z
̺(t, x) = dp f (t, x, p), ̺′ (t′ , x′ ) = dp′ f ′ (t′ , x′ , p′ ). (3.12b)
(ii) Ideal gas in a container. Consider a spatially homogeneous gas enclosed in a box
that rests in the lab frame Σ. Assume that the gas is in equilibrium and can be described
by a time-independent one-particle phase space density
φ(p)
I(x)
f (x, p) =
V
in Σ, where V is the rest volume of the box in Σ and I(x) the indicator for the box, i.e.,
I(x) = 1 if x is within the box and I(x) = 0, otherwise. In this case, Eq. (3.10) implies
that3
φ(p) φ′ (p′ )
= , (3.14)
V V′
where V ′ = V /γ(w) is the Lorentz-contracted box volume in the moving frame Σ′ . At first
sight it is surprising that the presence of the box alters the transformation properties of the
momentum distribution. However, this can be explained by the fact that the observations
by O and O′ are not synchronous, and that in the time between their observations some
particles collide with the container walls.4
(iii) Current-density vector. The quantities
Z Z
̺(t, x) = dp f (t, x, p) , j(t, x) = dp v f (t, x, p) (3.15)
can be combined into a current-density four-vector (j α ) = (̺, j), since they transform as
̺′ (t′ , x′ ) = γ(w) ̺(t, x) − γ(w) w j(t, x), (3.16a)
j ′ (t′ , x′ ) = −γ(w) w ̺(t, x) + γ(w) j(t, x). (3.16b)
Furthermore, it can be shown that j α satisfies the continuity equation
∂̺
∂α j α = + divj = 0. (3.17)
∂t
Equation (3.10) and statements (i) − (iii) remain valid in higher space dimensions d > 1
upon replacing x → x, p → p, j → j, etc. [138]. Moreover, Eq. (3.10) can be generalized
to the case of N-particle phase space PDFs fN , yielding for arbitrary space dimensions
fN (t1 , x1 , p1 ; . . . ; tN , xN , pN ) = fN′ (t′1 , x′1 , p′1 ; . . . ; t′N , x′N , p′N ), (3.18)
where for n = 1, . . . , N the coordinates (t′n , x′n , p′n ) and (tn , xn , pn ) are connected by a
Lorentz-transformation, and fN is the multiple-time probability density for lab observer O
to observe particle 1 at time t1 near (x1 , p1 ), and particle 2 at time t2 near (x2 , p2 ), etc..
The above results clarify the transformation behavior of PDFs in special relativity, but
they do not yet answer the question as to which PDF provides the correct description for
a given physical system as, e.g., a relativistic gas in equilibrium. The latter problem will
be addressed in the next part.
3
This result was already derived by Dirac [453] in 1924.
4
Cf. discussion in Section 6 of van Kampen’s paper [138].
40 CHAPTER 3. RELATIVISTIC EQUILIBRIUM THERMOSTATISTICS
space PDF5
where V is the volume of the gas container, E(m, p) = (m2 +p2 )1/2 = mγ(v) the relativistic
particle energy, and p = mvγ(v) the relativistic momentum; the corresponding marginal
momentum PDF reads6
φJ (p) = Z−1
J exp[−βJ E(p)]. (3.20c)
Equations (3.20) refer to a lab frame Σ where the box, enclosing the gas, is at rest.
Jüttner’s proposal (3.20) became widely accepted among theorists during the first three
quarters of the last century [112,162,169,273,289] – although a rigorous microscopic deriva-
tion remained lacking due to the difficulty of formulating a relativistically consistent Hamil-
ton mechanics of interacting particles [223, 224, 226, 228, 252]. Doubts about the Jüttner
function fJ began to arise in the 1980s, when Horwitz et al. [215, 216] suggested a ‘man-
ifestly covariant’ relativistic Boltzmann equation, whose stationary solution differs from
Eq. (3.20) and, in particular, predicts a different mean energy-temperature relation in the
ultra-relativistic limit [206]. Since then, partially conflicting results and proposals from
other authors [14,16,203,207,211] have led to an increasing confusion as to which distribu-
tion actually represents the correct generalization of the Maxwellian (3.19). For example,
another recently discussed alternative to Eq. (3.20), which also reduces to the Maxwell
distribution in the nonrelativistic limit case, is the ‘modified’ Jüttner function [14,16,211]
md γ(v)2+d
ψMJ (v; m, βMJ , d) = exp[−βMJ mγ(v)]. (3.21a)
ZMJ mγ(v)
For a spatially homogeneous gas, Eq. (3.21a) corresponds to the phase space PDF
(ZMJ )−1
φMJ (p) = exp[−βMJ E(p)], (3.21c)
E(p)
Compared with the Jüttner distribution (3.20) at the same parameter values
βJ = βMJ . 1/m, the modified PDFs (3.21) exhibits a significantly lower particle popu-
lation in the high energy tail because of the additional 1/E-prefactor.
5
Here and below, it will be assumed that f vanishes outside the gas container.
6
Jüttner [167] derived the distribution (3.20) from a maximum entropy principle; cf. Section 3.2.1. An
alternative derivation based on the microscopic ensemble is given by Matolcsi et al. [198].
42 CHAPTER 3. RELATIVISTIC EQUILIBRIUM THERMOSTATISTICS
It is worthwhile to note that the relative entropy (3.23) is manifestly invariant under coor-
dinate transformations [16]. In the mathematical literature, the definition of the relative
entropy is often given in the form
dµφ
Z
Ŝ[µφ |µρ ] := − dµφ ln , (3.24)
dµρ
which upon identifying
dµφ φ(p)
dµφ = dp φ(p) , dµρ = dp ρ(p) , =
dµρ ρ(p)
becomes equivalent to Eq. (3.23). The relative entropy is also known as the Kullback-
dµ
Leibler entropy [461], and the non-negative function dµφρ (p) is the Radon-Nikodym den-
sity [457] of µφ with respect to µρ . In the remainder, we are going to work with the density
representation (3.23).
7
See Ref. [16] for a discussion of the d-dimensional case.
8
Since we restrict ourselves to the spatially homogeneous case, it suffices to focus on the momentum
distribution φ.
9
In particular, if µφ (RM1 ) = 1 holds true then φ ≥ 0 is a PDF and µφ a probability measure on RM1 ;
at this point, however, we do not require normalization of φ or ρ.
3.2. THERMOSTATISTICS OF A RELATIVISTIC GAS 43
The relative entropy S[φ|ρ] from Eq. (3.23) constitutes the basis of the maximum rela-
tive entropy principle. This principle formalizes the idea that, for a broad class of physical
systems, the canonical equilibrium distribution can be obtained by maximizing an appro-
priate entropy functional in the presence of constraints. The constraints encode a priori
knowledge about the system. In our case, we are interested in maximizing S[φ|ρ] with
respect to φ under the conditions
Z
1 = dp φ(p), (3.25a)
Z
ǫ = dp E(p) φ(p). (3.25b)
The first constraint (3.25a) ensures that φ is a PDF. The second constraint (3.25b) reflects
the assumption that the mean energy per particle, ǫ, is known. These constraints may be
incorporated into the maximum entropy principle via the method of Lagrangian multipli-
ers [462]. Denoting the Lagrangian multipliers by (α, β), the maximum entropy principle
results the necessary condition
δ
Z Z
0 ≡ S[φ|ρ] + α 1 − dp φ(p) + β ǫ − dp E(p) φ(p)
δφ φ=φ∗
φ∗ (p)
= − 1 + ln − α − βE(p), (3.26)
ρ(p)
The Lagrangian multipliers (α, β) can be determined from the two conditions (3.25) yield-
ing, e.g., β as a function of the given mean value ǫ. Moreover, it is evident now that
normalization of the reference density ρ is irrelevant, since constant prefactors will be
absorbed by the multiplier α. We next discuss three specific cases:
b) Jüttner distribution. Fixing again a constant reference density, ρ(p) ≡ ρ0 , and consid-
ering the relativistic energy E(p) = (m2 + p2 )1/2 , one recovers Jüttner’s momentum
PDF (3.20c).
10
Due to the appearance of the logarithm, the maximization of the entropy functional (3.23) gives rise
to an ‘exponential’ distribution (3.27). By considering other non-logarithmic ‘entropies’ like, e.g., Tsallis’
or Kaniadakis’ entropy, one may construct other forms of distributions (e.g., power law distributions).
44 CHAPTER 3. RELATIVISTIC EQUILIBRIUM THERMOSTATISTICS
This shows that the two candidate distribution (3.20) and (3.21) may be derived from a
common maximum relative entropy principle, but they refer to different reference densities
(i.e., reference measures), respectively. It is worthwhile to clarify the difference between
the two reference measures by analyzing their respective symmetries:
The constant reference density ρ(p) ≡ ρ0 , underlying the Jüttner distribution (3.20),
corresponds to the Lebesgue measure λ on RM1 , i.e., in this case we have
φ(p)
Z
S[φ|ρ0 ] := − dp φ(p) ln (3.29)
ρ0
coincides with the usual Shannon-Boltzmann-Gibbs entropy on RM1 . The Lebesgue mea-
sure λ on RM1 is distinguished by the fact that it is the only translation invariant measure
on the relativistic momentum space. Here, translation invariance of λ means that
holds for all a, b, x ∈ RM1 . In more mathematical terms, λ represents the Haar measure11
of the (additive) momentum translation group (RM1 , +).
11
In a seminal paper [463] published in 1933, the Hungarian mathematician Alfred Haar studied the
possibility to introduce a measure µ◦ on a continuous group (G, ◦) such that µ◦ is invariant under the group
multiplication ‘◦’. To briefly sketch this idea, consider a subset A of the group G and some arbitrary, fixed
group element g ∈ G. By multiplying each element a ∈ A with g, the subset A is mapped onto another
subset of G, denoted by
g ◦ A := {g ◦ a |a ∈ A } .
Now consider a measure µ◦ on G that assigns to A ⊆ G some non-negative real number µ◦ (A). The
measure µ◦ is said to be group invariant, if
µ◦ (g ◦ A) = µ◦ (A)
holds for any g ∈ G and A ⊆ G. In the case of non-commutative (i.e., non-Abelian) groups, one may
distinguish invariance under multiplications from the right or left. Haar [463] was able to prove the
existence of an invariant measure µ◦ , and its uniqueness apart from an irrelevant multiplicative constant
for locally compact, topological groups. Such group invariant measures µ◦ are referred to as Haar measures
nowadays [457]. They give a mathematically precise meaning to the notion ‘uniform distribution’ by
combining measure and group theoretical concepts.
3.2. THERMOSTATISTICS OF A RELATIVISTIC GAS 45
Now consider the reference density ρ(p) = E(p)−1 , which yields the modified momentum
PDF (3.21c). We may define the Lorentz transformation L[A] of a set A ⊂ RM1 by
dp dp′
= (3.31)
E(p) E(p′ )
That is, the measure µρ induced by ρ(p) = E(p)−1 is invariant under Lorentz transforma-
tions.
To briefly summarize, adopting the Lebesgue measure in momentum space as reference
measure yields the Jüttner function (3.20c), whereas the Lorentz invariant reference mea-
sure yields the modified distribution (3.21c). Hence, at this stage, both distribution func-
tions appear to be plausible candidates,12 and it remains the question which reference
measure is the physically correct one. This question cannot be answered on the level of
the maximum entropy principle and, therefore, other approaches have to be explored. To
identify the physically correct distribution, we have performed numerical simulations of a
fully relativistic dynamical 1D gas model that will be discussed in the next section.
The basic time step of the algorithm involves three partial tasks:
(i) determine the next collision event (xc , tc );
(ii) evolve the system up to time tc ;
(iii) calculate the momenta after the collision.
The third task is solved as follows: If two particles A and B meet at the space-time point
(xc , tc ), then they exchange momentum according to the relativistic energy momentum
conservation laws
where E(m, p) = (m2 + p2 )1/2 is the relativistic particle energy. Given the momenta
(pA , pB ) before the collision, these conservation laws determine the momenta (p̂A , p̂B ) after
the collision by [14]
where
microcanonical ensemble, since the total initial energy Etot in Σ is conserved in the mi-
croscopic collision processes. The above conventions define the simplest interacting model
system that
a) complies with all principles of special relativity,
Assuming that (i) an equilibrium state exists where both species can be described by the
same value β, and that (ii) for a gas in equilibrium the mean energy per particle is the
same for particles of the same species, the total energy can be expressed as
Etot = N1 ǫ(m1 , β) + N2 ǫ(m2 , β). (3.36)
48 CHAPTER 3. RELATIVISTIC EQUILIBRIUM THERMOSTATISTICS
In our case, the energy mean values ǫJ/MJ of the two 1D candidate PDFs ψJ and ψMJ can
be calculated analytically, yielding [cf. App. B.1]
K0 (βJ m) + K2 (βJ m)
ǫJ (m, βJ ) = m , (3.37a)
2K1 (βJ m)
K1 (βMJ m)
ǫMJ (m, βMJ ) = m , (3.37b)
K0 (βMJ m)
with Kn denoting modified Bessel functions of the second kind [466]. For each simulation
run the parameter tuple (Etot , N1 , N2 , m1 , m2 ) is known. Hence, upon inserting them into
[c−1 ]
1 heavy particles ψJ
ψMJ
velocity PDF
0.5
ǭ = 2.5 m1 c2 w = 0.0
0
[c−1 ]
light particles ψJ
2 ψMJ
velocity PDF
0
−1 −0.5 0 0.5 1
v/c
Figure 3.1: Equilibrium velocity PDFs in the lab frame Σ: Numerically obtained one-particle
velocity PDFs (◦) based on a simulation with N1 = 1000 light particles of mass m1 and N2 = 1000
heavy particles with mass m2 = 2m1 . The mean energy per particle is ǭ = Etot /(N1 + N2 ) =
2.5m1 c2 . The solid curves in the upper and lower panel correspond to Jüttner functions (3.20a)
with same inverse temperature parameter βJ = 0.701 (m1 c2 )−1 , but different particle masses,
respectively. Dashed lines show the corresponding modified distribution (3.21a) with βMJ =
0.402 (m1 c2 )−1 . The simulation data is consistent with the standard Jüttner distribution ψJ
(solid lines), and thus provides evidence against the modified distribution (3.21).
3.2. THERMOSTATISTICS OF A RELATIVISTIC GAS 49
Eqs. (3.36) and (3.37), these parameters uniquely determine the parameter value βJ/MJ
that is consistent with the chosen velocity PDF fJ/MJ .
However, for this concept to be meaningful, a restriction of the accessible spatial volume
is required – be it by means of periodic boundary conditions, or by imposing reflecting
walls.14 Otherwise, it cannot be expected that a many-particle system approaches a uni-
versal stationary state which is independent of the specific initial conditions. This obser-
vation has an important implication: Any (relativistic or non-relativistic) Boltzmann-type
equation [206, 215, 216, 273, 289, 290, 292] that gives rise to a universal stationary velocity
PDF implicitly assumes the presence of a spatial confinement, thus singling out a preferred
frame of reference.
14
The critical role of the boundary conditions in relativistic systems has been emphasized by
Sinyukov [467] and van Kampen [137,138]. Loosely speaking, if a many-particle system has reached a uni-
versal stationary equilibrium state, then each particle ‘knows’ about the presence of the confinement/walls
because equilibration typically requires momentum reversal at the walls, e.g., in order to maintain a uni-
form density and a well-balanced average collision frequency. In particular, if the walls are considered to
be stationary objects then they single out a preferred frame of reference. The relevance of the boundary
conditions is even more obvious in quantum mechanics/statistics due to their direct effect on the energy
spectra and, thus, on the density of states. In fact, quantum mechanical arguments [170] suggest ρ(p) ≡ ρ0
in Eq. (3.23).
50 CHAPTER 3. RELATIVISTIC EQUILIBRIUM THERMOSTATISTICS
Moving observers From our simulations we may further determine the equilibrium
velocity distributions as seen from another frame Σ′ moving with velocity w relative to the
lab frame Σ. Figure 3.2 depicts results for w = 0.2 and same simulation parameters as
in Fig. 3.1. In contrast to Fig. 3.1, the numerical data points in Fig. 3.2 were obtained
by measuring velocities Σ′ -simultaneously. The solid curves in Fig. 3.2 correspond to the
PDF
m γ(v ′ )3
ψJ′ (v ′ ; m, βJ , w) = exp[−βJ γ(w)mγ(v ′) (1 + wv ′)] (3.39)
ZJ γ(w)
[v ′ is the particle velocity in the moving frame Σ′ ]. The velocity PDF (3.39) is obtained by
making use of Eq. (3.14). For w = 0, the PDF ψJ′ reduces to the Jüttner function (3.20).
Due to the excellent agreement between the numerical simulations and Eq. (3.39), we may
state more precisely: Two relativistic gas components are in ‘thermodynamic equilibrium’
for any observer if their one-particle velocity PDFs are given by generalized Jüttner func-
tions (3.39) with same parameters βJ and w. Only in this case the net energy transfer
between the different gas components in the container vanishes. Last but not least, the
above results shed light on a longstanding, highly debated question [127, 129–131, 158]
originally posed by Landsberg [126, 128]:
Does a moving body appear cool? Evidently, the answer depends on the ther-
mometers employed by different observers. Adopting, for the reasons discussed above,
T := (kB βJ )−1 as a reasonable temperature definition, a moving observer with rest frame
Σ′ can measure T by exploiting the Lorentz invariant equipartition theorem15 [128]
′
kB T = mγ(w)3 γ(v ′ ) (v ′ + w)2 , (3.40a)
where
w = −hv ′ i′ (3.40b)
is the mean velocity of the gas measured by the moving observer, and the averages h · i′
are taken Σ′ -simultaneously. We verified the validity of Eq. (3.40a) explicitly by using
simulation data obtained for different values of w, see Fig. 3.3. Hence, Eq. (3.40a) defines
a Lorentz invariant statistical gas thermometer. Put differently, this intrinsic statistical
thermometer determines the proper temperature of the gas by making use of simulta-
neously measured particle velocities only; thus, by adopting the statistical thermometer
definition (3.40a), moving bodies appear neither hotter nor colder.16
15
Equation (3.40a) is obtained by combining the microcanonical equipartition theorem for a Hamiltonian
PN1 PN2
H = i=1 E(m1 , pi ) + j=1 E(m2 , pj ) with the Lorentz invariance of the relativistic phase space PDF f .
16
The mean value from Eq. (3.40a) can be used to measure the rest temperature, which plays a cen-
tral role in van Kampen’s [137] approach to relativistic thermodynamics. Evidently, upon multiplying
3.2. THERMOSTATISTICS OF A RELATIVISTIC GAS 51
1.5
[c−1 ]
heavy particles ψJ′
1
velocity PDF
0.5
ǭ = 2.5 m1 c2 w = 0.2
0
4
[c−1 ]
0
−1 −0.5 0 0.5 1
v/c
Figure 3.2: Equilibrium velocity PDFs in a moving frame Σ′ : Velocity PDFs as measured by
an observer who moves with velocity w = 0.2c relative to the lab frame Σ. Parameter values
and initial conditions are the same as those in Fig. 3.1. The solid lines correspond to Jüttner
functions ψJ′ from Eq. (3.39) with the same inverse temperature parameter βJ = 0.701 (m1 c2 )−1
as in Fig. 3.1 and different masses m1 and m2 , respectively.
To briefly summarize the results of this part: Our fully relativistic MD simulations confirm
the Jüttner distribution (3.20) as the correct relativistic one-particle equilibrium velocity
distribution. Furthermore, our simulations corroborate van Kampen’s [137] and Lands-
berg’s view [126,128] that the temperature of classical gaseous systems can be defined and
measured in a Lorentz invariant way.
The extension of the MD approach to higher space dimensions is nontrivial, due to the
aforementioned difficulties of treating 2D and 3D two-body collisions in a relativistically
consistent manner [223, 224, 226, 228, 252]. In order to be fully consistent with the require-
Eq. (3.40a) by factors γ(w)α , α 6= 0 one can construct thermometers that measure ‘other’ temperatures;
e.g., a = −1 would correspond to Planck’s [111] formulation of relativistic thermodynamics and a = 1 to
proposals made by Eddington [468] and Ott [115].
52 CHAPTER 3. RELATIVISTIC EQUILIBRIUM THERMOSTATISTICS
0 1.6
kB T / (m1 c2 )
hv ′ i′ / c
−0.2
1.4
−0.4
1.2
−0.6
0 0.2 0.4 0.6 0 0.2 0.4 0.6
w/c w/c
Figure 3.3: Measured mean particle velocity (◦, left digram) and estimated temperature (◦, right
diagram) based on Eqs. (3.40) as a function of the observer velocity w. Solid lines correspond to
the theoretically expected values, respectively, using the same simulation parameters and initial
conditions as in Fig. 3.2.
17
Typical examples are unconfined systems with a limited number of collisions per particle as, e.g., in
an expanding dilute gas.
Chapter 4
The present chapter discusses the generalization of the Langevin theory of Brownian mo-
tions to the framework of special relativity [10, 11, 13, 14, 16, 18–24, 331, 333–335, 374, 470].
More precisely, we will consider stochastic differential equations (SDEs) that describe
Markov processes in relativistic one-particle phase space.1 Relativistic Langevin equa-
tions present a useful tool for modeling the dynamics of relativistic particles in a ran-
dom environment. Recently, for example, SDEs have been applied to analyze thermal-
ization effects in quark-gluon plasmas produced at the Relativistic Heavy Ion Collider
(RHIC) [312–314, 337].2 The subsequent sections intend to give a comprehensive discus-
sion and illustration of the underlying mathematical theory.
The phenomenological or axiomatic Langevin approach to relativistic Brownian motion
was initiated by Debbasch et al. [18], who in 1997 proposed a simple relativistic generaliza-
tion of the classical Ornstein-Uhlenbeck process [36]. The relativistic Ornstein-Uhlenbeck
process (ROUP) of Debbasch et al. [18] is obtained by postulating additive white noise for
the particle’s momentum change in the rest frame of the bath. During the past decade
various properties of the ROUP were studied by Debbasch and Rivet [19, 20], Barbachoux
et al. [21,22], and Zygadlo [334]. An alternative approach to relativistic Langevin equations
was pursued by us in Refs. [10, 11]. By starting from a ‘Newtonian’ Ornstein-Uhlenbeck-
type Langevin equation in the comoving rest frame of a Brownian particle, we obtained
a modified relativistic Brownian motion (RBM) process, whose relaxation behavior differs
from that of the ROUP. As we are going to illustrate below, the two different processes
represent special limit cases of a larger class of relativistic Langevin models [12], which
may be used to describe the random motions of relativistic particles.
The content of the present chapter is structured as follows. Section 4.1 discusses the
axiomatic Langevin approach, i.e., suitable SDEs and fluctuation-dissipation relations are
1
As discussed by Lopuszaǹski [339], Dudley (Theorem 11.3 in [341]) and Hakim (Proposition 2 in [346]),
it is impossible to define nontrivial relativistic Markov processes in position space.
2
Potential applications in high energy astrophysics are discussed by Dieckmann et al. [338] and [471].
53
54 CHAPTER 4. RELATIVISTIC BROWNIAN MOTION
postulated in order to provide a simplified model of the complex interaction between the
Brownian particle and its environment (heat bath). After having outlined the general
conceptual foundations (Section 4.1.1), specific examples are considered (Section 4.1.2).
In Section 4.1.3, we will analyze the temperature dependence of the asymptotic mean
square displacement for different example processes, using general analytic formulas re-
cently derived by Lindner [24], and Angst and Franchi [23]. In particular, we are going
to demonstrate that the diffusion constant is sensitive with respect to variations of the
friction coefficients. This result implies that measurements of the diffusion constant may
reveal information about the underlying microscopic interactions. Section 4.2 discusses
relativistic Brownian motion processes from the viewpoint of a moving observer. In the
last part of this chapter (Section 4.3), we will generalize the binary collision model from
Section 2.2.2 to the relativistic case in order to obtain a simple microscopic model for
relativistic Brownian motions.
When considering Langevin equations as models of Brownian motion, one implicitly as-
sumes that it is possible and reasonable to separate the degrees of freedom of the Brownian
particle from those of the environment (heat bath). Adopting this point of view, one can
specify two distinguished frames of reference: the stationary inertial rest frame Σ of the
heat bath, and the inertial frame Σ∗ that is comoving with the Brownian particle at a
given instant of time.3 As before, Σ is referred to as lab frame. In the present section, we
focus on describing relativistic Brownian motion processes with respect to the space-time
coordinates of Σ.
3
Apart from an irrelevant shift of the origin, the inertial lab frame Σ is uniquely determined by the
requirement that the mean velocity of the heat bath particles, which is assumed to be constant in any
inertial frame, must vanish in Σ. Similarly, the instantaneously comoving frame Σ∗ is determined by
the condition that the Σ∗ -velocity of the Brownian particle is equal to zero at the given instant of time.
Generally, we assume that the time coordinates t and t∗ can be measured, e.g., by using atomic clocks
that are at rest in Σ or Σ∗ , respectively.
4.1. LANGEVIN AND FOKKER-PLANCK EQUATIONS 55
|P | |P |
|V | = =
P 0
(M + P 2 )1/2
2
is always less than c = 1 (throughout, M > 0 denotes the rest mass of the Brownian
particle). However, before we can actually write down specific Langevin equations for
the relativistic momentum components P i, an important question needs to be addressed,
namely, the choice of the time parameter in relativistic Langevin equations.
relation
where B(t) is a standard Wiener process with increment PDF (2.3d). If the background
medium (heat bath) is stationary and spatially homogeneous, then both the friction co-
efficient α and the noise amplitude D should depend on the relativistic particle energy
E(P ) = (M 2 + P 2 )1/2 only, i.e.,
In the remainder, we will always assume that Eqs. (4.2b) hold true. An additional con-
straint results from the requirement that Eq. (4.2a) should yield the correct stationary
momentum distribution.
The FPE for the phase space density f (t, x, p) of a relativistic Brownian particle, governed
by Eq. (4.1) with d = 1 and Eq. (4.2a), reads
∂f p ∂f ∂ ∂f
+ = α• (p) p f + D(p) . (4.2c)
∂t E ∂x ∂p ∂p
Analogous to Eq. (2.27), the corresponding stationary distribution is given by6
Z p ′
′ α• (p ) ′
f∞ (x, p) = N exp − dp p . (4.2d)
p∗ D(p′ )
5
In principle, one could also consider other driving process (as, e.g., Levy or Poisson processes) and/or
coefficient functions α(t, X, P ) and D(t, X, P ).
6
Here, we assume the presence of a spatial confinement or periodic boundary conditions.
4.1. LANGEVIN AND FOKKER-PLANCK EQUATIONS 57
α• (p) ! d d
p ≡ − log φJ (p) = βE(p). (4.2e)
D(p) dp dp
α• (p) β
≡ . (4.2f)
D(p) E(p)
This is the relativistic fluctuation-dissipation relation, also referred to as the relativistic
Einstein relation. Compared with the nonrelativistic Einstein relation (2.29), the mass has
been replaced with energy on the rhs. of Eq. (4.2f).
Furthermore, it is straightforward to derive from Eq. (4.2a) the corresponding SDE for the
relativistic energy P 0 = E(P ) = (M 2 +P 2)1/2 . Applying the (backward) Ito formula (C.27)
with Y = P , Z = P 0 and G(p) = (M 2 + p2 )1/2 , we obtain the following SDE for the
relativistic energy P 0 = E(P ):
( " 2 # 2 )
0
M D̂(P ) M
dP 0 (t) = −α̂• (P 0) P 0 1 − 0
− 0
dt +
P P P0
( " 2 #)1/2
M
2 D̂(P 0 ) 1 − • dB(t), (4.2g)
P0
where the friction coefficient functions α◦ and α• are related by [cf. Eq. (C.43)]7
and D ′ (p) := dD(p)/dp. The relativistic Einstein relation (4.2f), reexpressed in terms of
α◦ (p), reads
βD(p) 1 D ′ (p)
α◦ (p) ≡ − , (4.3c)
E(p) 2 p
7
See also the corresponding discussion by Hänggi [473], and Hänggi and Thomas (page 293 of Ref. [67]).
58 CHAPTER 4. RELATIVISTIC BROWNIAN MOTION
i.e., only if α◦ (p) and D(p) satisfy the criterion (4.3c), then the stationary momentum
distribution is given by Jüttner’s PDF. Moreover, if D and α• depend on the relativistic
particle energy P 0 = E only, i.e., if D(p) = D̂(E) and α• (p) = α̂• (E) hold true, then we
may write
The rules of ordinary differential calculus are preserved for Stratonovich-Fisk SDEs. Con-
sequently, we find the following Stratonovich-Fisk SDE for the energy P 0 = E(P ):
" 2 #
0 0 0 M
dP (t) = −α̂◦ (P ) P 1 − dt +
P0
( " 2 #)1/2
M
2 D̂(P 0) 1 − ◦ dB(t). (4.3f)
P0
Pre-point discretization The relativistic Langevin equations (4.2a) and (4.3a) can also
be rewritten in terms of the equivalent Ito SDE
where the friction coefficients α∗ and α• are related by [cf. Eq. (C.43)]
Compared with Eqs. (4.2a) and (4.3a), the Ito SDE (4.4) is most convenient for numerical
simulations. For a homogeneous isotropic bath with D(p) = D̂(E) and α• (p) = α̂• (E),
Eq. (4.4b) is equivalent to
Applying the Ito formula (C.10), we obtain the following Ito SDE for the energy P 0
( " 2 # 2 )
D̂(P 0 ) M
M
dP 0 (t) = −α̂∗ (P 0) P 0 1− + dt +
P0 P0 P0
( " 2 #)1/2
M
2 D̂(P 0 ) 1 − ∗ dB(t). (4.4e)
P0
Having outlined the general ideas underlying the axiomatic Langevin approach to relativis-
tic Brownian motions in the lab frame, we next consider several example processes.
4.1.2 Examples
We discuss three specific 1D relativistic Langevin models whose stationary momentum
distributions are given by the Jüttner function φJ (p) = exp[−βE(p)]. In this case, the
relativistic Einstein relation (4.2f) implies that only one of the two functions α• (p) and
D(p) can be chosen arbitrarily.
where αc > 0 is a constant friction parameter. From the relativistic Einstein relation (4.2f),
one then finds
i.e., the ROUP corresponds to the limit case of constant noise amplitude. The associated
Langevin equation reads
1/2
M 2αc M
dP (t) = −αc 0 P dt + • dB(t) (4.5c)
P β
1/2
M 2αc M
= −αc 0 P dt + ∗ dB(t). (4.5d)
P β
The discretization rule is irrelevant here, because the noise amplitude Dc = αc M /β does
not depend on the momentum P for this particular case. However, the rules of stochastic
calculus have to be specified, if one wishes to write down the SDE for the associated velocity
60 CHAPTER 4. RELATIVISTIC BROWNIAN MOTION
process V (t) := P/P 0 . For example, adopting the post-point discretization, the Langevin
equation for V (t) reads
2 3/2 3kB T 2 2
dV (t) = −αc (1 − V ) − (1 − V ) V dt +
M
1/2
2αc kB T 2 3
(1 − V ) • dB(t). (4.5e)
M
For comparison, the corresponding Ito SDE is given by
2 3/2 3kB T 2 2
dV (t) = −αc (1 − V ) + (1 − V ) V dt +
M
1/2
2αc kB T 2 3
(1 − V ) ∗ dB(t). (4.5f)
M
Adopting the post-point discretization scheme, the relativistic Langevin equation of this
model reads
1/2
2α† P 0
dP (t) = −α† P dt + • dB(t). (4.6b)
β
The corresponding SDE for the velocity process V (t) = P/P 0 is given by
2 3kB T 2 3/2
dV (t) = −α† (1 − V ) − (1 − V ) V dt +
M
1/2
2αc kB T 2 5/2
(1 − V ) • dB(t). (4.6c)
M
Recently, various properties of the RBM process (4.6b) have been analyzed by Fa [331],
Lindner [24], Fingerle8 [470], and Angst and Franchi [23].
With regard to numerical simulations,9 the Ito form of Eqs. (4.6b) is often more conve-
nient. By making use of Eq. (4.4c), the equivalent Ito SDEs are obtained as
0 1/2
2α† P 0
βP − 1
dP (t) = −α† P dt + ∗ dB(t), (4.6d)
βP 0 β
8
Fingerle [470] discusses a fluctuation theorem for this process; see also Cleuren et al. [218].
9
Cf. App. C.
4.1. LANGEVIN AND FOKKER-PLANCK EQUATIONS 61
and
2 2kB T 2 3/2
dV (t) = −α† (1 − V ) + (1 − V ) V dt +
M
1/2
2αc kB T 2 5/2
(1 − V ) ∗ dB(t). (4.6e)
M
Constant friction coefficient in the Ito-SDE The RBM process defined by Eq. (4.6b)
is characterized by a constant friction coefficient α† , when adopting the post-point dis-
cretization rule (•). Another model, referred to as RBM(I) hereafter, is obtained by con-
sidering a constant friction coefficient α∗ in the Ito-Langevin equation
1/2
2α∗
dP (t) = −α∗ P dt + (1 + βP 0 ) ∗ dB(t), (4.7a)
β2
where the noise amplitude is chosen such that the Einstein relation (4.4d) is satisfied. For
completeness, the Ito-SDE of the associated velocity process V (t) := P/P 0 is given by
2
2 kB T 2 3/2 kB T 2 2
dV (t) = −α∗ (1 − V ) + 3 (1 − V ) + 3 (1 − V ) V dt +
M M
1/2
kB T 2 5/2 kB T 2 3
2α∗ (1 − V ) + (1 − V ) ∗ dB(t). (4.7b)
M M
The three model processes (4.5d), (4.6b) and (4.7a) give rise to the same stationary
momentum PDF φJ (p), but their respective relaxation behavior differs strongly. This
is illustrated in Fig. 4.1, which depicts the time evolution of the spatial mean square
displacement divided by time,
for all three models at same temperature T. The curves in Fig. 4.1 were calculated numer-
ically from Eqs. (4.5d), (4.6d) and (4.7a), respectively, using an algorithm similar to those
described in [23, 24], see also App. C.5.
1.5
M c2 /(kB T) = 1.0
1
Dt / (c2 /α)
0.5
0
0 100 200 300 400 500
−1
t/α
Figure 4.1: Time evolution of the spatial mean square displacement Dt := [X(t) − X(0)]2 /(2t)
for the ROUP [18] model (solid line) from Eq. (4.5d), the RBM [10] model (dotted) from Eq. (4.6b)
and the RBM(I) model (dashed) from Eq. (4.7a) at same temperature kB T/(M c2 ) = 1. The plots
are based on a simulation with N = 1000 trajectories, initial conditions X(0) = 0, P (0) = 0 for
each trajectory, and discretization time step ∆t = 10−4 α−1c/† .
The asymptotic diffusion constant D∞ may be expressed in terms of the velocity correlation
function hV (t)V (s)i by virtue of
1 d
D∞ = lim [X(t) − X(0)]2
t→∞ 2 dt
1 d t
Z Z t
= lim ds ds′ hV (s)V (s′ )i
t→∞ 2 dt 0 0
Z t
= lim ds hV (t)V (s)i . (4.10)
t→∞ 0
Assuming that the velocity process V (t) is (approximately) stationary, which means that
hV (t)V (s)i = hV (t − s)V (0)i holds (at least in good approximation), and substituting
u = t − s, we recover Kubo’s formula
Z t
D∞ = lim du hV (u)V (0)i . (4.11)
t→∞ 0
with symmetric coefficient functions, a• (v) = a• (−v) and b(v) = b(−v), the Kubo for-
mula (4.11) gives rise to the following integral representation for the asymptotic diffusion
constant:10
R v+ hR i2
U (y) v+ −U (x)
0
dy e y
dx e x/b(x)
D∞ = R v+ . (4.13)
0
dz e−U (z) /b(z)
Here, v+ ∈ [0, ∞] represents the upper bound for the velocity range, and
Z v
U(v) := dw µ∗ (w)/b(w) (4.14a)
0
In general, the formula (4.13) has to be integrated numerically, but for the first two models
from Section 4.1.2 the integrals may also be evaluated analytically.
For example, upon comparing with Eq. (4.6e), we see that the RBM model from Eq. (4.6b)
is described by
2 2kB T 2 3/2
µ∗ (v) = α† (1 − v ) + (1 − v ) v, (4.15a)
M
kB T
b(v) = α† (1 − v 2 )5/2 , (4.15b)
M
with an upper velocity bound v+ = c = 1. In this case, the rhs. of Eq. (4.13) can be
calculated analytically by making use of the identity
ν Z 1
ν d χ
(−1) K0 (χ) = dv exp − √ (1 − v 2 )−(ν+2)/2 , (4.16)
dχν 0 1 − v 2
where, for ν = 0, 1, 2 . . ., Kν (z) denotes the modified Bessel function of the second
kind [466]. One then finds that11
K0 (βM)
DRBM
∞ = (α† βM)−1 . (4.17)
K1 (βM)
At low temperatures β := (kB T)−1 → ∞, Eq. (4.17) reduces to the well-known classical
result Dclass
∞ = kB T/(Mα† ), cf. Eq. (2.20). In the opposite limit of very high temperatures,
i.e., for βM ≪ 1, one finds a logarithmic dependence [23]
RBM −1
2 2
D∞ = (α† M) −γǫ + ln + O[(βM) ] , (4.18)
βM
10
A d-dimensional generalization of Eq. (4.13) was recently derived by Angst and Franchi [23].
11
Equation (4.17) is an equivalent, more compact representation of Lindner’s result Eq. (10) in [24].
64 CHAPTER 4. RELATIVISTIC BROWNIAN MOTION
where γǫ ≃ 0.577216 is the Euler constant. However, it should be kept in mind that, due to
the increasing importance of particle annihilation/creation at high energies, classical non-
quantum theories become invalid in the high temperature limit βM ≪ 1, and, therefore,
the asymptotic expansion (4.18) appears to be of limited practical use.
For comparison, the ROUP (4.5d) corresponds to [cf. Eq. (4.5f)]
2 3/2 3kB T 2 2
µ∗ (v) = αc (1 − v ) + (1 − v ) v, (4.19a)
M
kB T
b(v) = αc (1 − v 2 )3 . (4.19b)
M
Remarkably, in this case the general integral formula (4.13) and the identity (4.16) yield
the ‘classical’ result
DROUP
∞ = kB T/(Mαc ) = (αc βM)−1 (4.20)
for all parameter values (αc , T, M). Moreover, we observe that DRBM
∞ ≤ DROUP
∞ holds true
for same values of the friction coefficients αc = α† . Intuitively, this can be explained by
the fact that, for the ROUP, the absolute value of the friction force is bounded by αc M,
cf. Eq. (4.5d), whereas the friction force is unbounded for the RBM model (4.6b), thereby
suppressing spatial diffusion more strongly in the latter case.
As the last example, we consider the RBM(I) model, defined in Eq. (4.7a) and described
by [cf. Eq. (4.7b)]
2
2 kB T 2 3/2 kB T 2 2
µ∗ (v) = α∗ (1 − v ) + 3 (1 − v ) + 3 (1 − v ) v, (4.21a)
M M
kB T 2 5/2 kB T 2 3
b(v) = α∗ (1 − v ) + (1 − v ) . (4.21b)
M M
The remaining integral can be evaluated numerically. As illustrated in Fig. 4.2, the the-
oretical predictions from Eqs. (4.17), (4.20) and (4.22) are in good agreement with the
numerically obtained estimates of the asymptotic diffusion constant.
4.1. LANGEVIN AND FOKKER-PLANCK EQUATIONS 65
The ROUP and the two RBM models considered in this part represent special limit
cases of the general Langevin equation (4.2a) with arbitrarily chosen friction coefficient
functions. Nonetheless, they yield useful insights: As evident from Fig. 4.2, at moderate-
to-high temperatures the diffusion constant can vary significantly for different friction
models. For realistic systems, the exact functional shape (i.e., energy dependence) of
the friction function α is determined by the microscopic interactions. This result implies
that simultaneous measurements of the temperature and the diffusion constants can reveal
information about the underlying microscopic forces. Below, in Section 4.3, we will outline
a general procedure for deducing more realistic friction coefficients α from microscopic
models.
10 ROUP
RBM
RBM(I)
D∞ / (c2 /α)
0.1
0.01
0.01 0.1 1 10
kB T / (M c2 )
Figure 4.2: Temperature dependence of the asymptotic diffusion constant D∞ for the ROUP [18]
from Eq. (4.5d), the RBM model [10] from Eq. (4.6b), and the RBM(I) model from Eq. (4.7a).
Symbols ’◦’, ’+’ and ’×’ represent the results of computer simulations (N = 100 trajec-
tories, initial conditions X(0) = 0, P (0) = 0 for each trajectory, discretization time step
∆t = 10−4 in units of α−1 c/†/∗ ), obtained by averaging over the numerically determined values
{D100 , D110 , D120 , . . . , D500 }. Solid, dashed and dotted lines correspond to the theoretical predic-
tions from Eqs. (4.20), (4.17) and (4.22), respectively.
66 CHAPTER 4. RELATIVISTIC BROWNIAN MOTION
and demand that α and D satisfy the relativistic Einstein relation βD(p) = α(p)E(p) with
E(P ) = (M 2 + P 2 )1/2 = P 0 . Then, upon multiplying by p0 = (M 2 + p2 )1/2 , the FPE (4.2c)
for the associated PDF f (t, x, p) can be written in the form
0 ∂ ∂ 0 ∂ ∂
p +p f =p α(p) p f + D(p) f . (4.24)
∂t ∂x ∂p ∂p
As before, we will assume deterministic initial conditions X(0) = x0 , P (0) =
p0 for the stochastic process (4.23), corresponding to a localized initial distribution
f (0, x, p) = δ(x − x0 ) δ(p − p0 ) in Σ.
In order to characterize the process (4.23) from the viewpoint of a moving inertial
frame Σ′ , one can proceed as follows: According to van Kampen [138], the one-particle
phase space PDF f transforms as a Lorentz scalar, i.e., we have
and, conversely,
with w denoting the velocity of Σ′ relative to Σ, γ(w) = (1−w 2 )−1/2 and p′ 0 = (M 2 +p′ 2 )1/2 .
Thus, in order to find f ′ , one merely needs to solve the Fokker-Planck equation (4.24) in
the lab frame Σ and, subsequently, insert the solution into Eq. (4.25a).
4.2. MOVING OBSERVER 67
where hat-symbols refer to the state after the collision. Inserting Eqs. (4.36) into the
conservation laws (4.37), and solving for the momentum of Brownian particle after the
collision, P̂ , we obtain [13]
∆Pr := P̂ − P
of the Brownian particle in a single collision with the heat bath particle ‘r’ is given by
ǫr E
∆Pr = −2γ(ur )2 P + 2γ(ur )2 pr , (4.39)
E + ǫr E + ǫr
where ur := u(pr , P ) and ǫr := ǫ(pr ). In the non-relativistic limit case, where u2r ≪ 1,
E ≃ M and ǫr ≃ m, Eq. (4.39) reduces to Eq. (2.43).
Furthermore, by making the same assumptions as in Section 2.2.2, we find that the
momentum change δP (t) := P (t + δt) − P (t) of the Brownian particle during a small-but-
sufficiently-long time interval [t, t + δt] can be approximated by
N
X
δP (t) ≈ ∆Pr Ir (t, δt)
r=1
N N
X ǫr2
X E
≈ −2 γ(ur ) P (t) Ir (t, δt) + 2 γ(ur )2 pr Ir (t, δt). (4.40a)
r=1
E + ǫr r=1
E + ǫr
Formally, the collision indicator Ir (t, τ ) is again given by [cf. Eq. (2.48a)]
δt
Ir (t, δt) ≈ |vr − V | δ(xr − X), (4.40b)
2
but now we have to use the relativistic velocities
Equation (4.40a) is the relativistic counterpart of Eq. (2.48b). Heuristically, the first
term on the rhs. of Eq. (4.40a) can again be interpreted as ‘friction’, while the second
contribution may be viewed as ‘noise’.
Bath distribution and drift Similar to the nonrelativistic case, Eqs. (4.40) can be used
to calculate the statistical moments h[δP (t)]j ib of the momentum increments – provided
one specifies the phase space distribution of the heat bath particles. We will assume here
that the heat bath is in a thermal equilibrium state, so that the one-particle phase space
PDF is given by a spatially homogeneous Jüttner function
(m2 + p2r )1/2
1 −1
fb (xr , pr ) = (ZJ L) exp − , (4.41a)
kB T
where L is the 1D volume, T = (βkB )−1 the temperature, and
ZJ = 2m K1 (βm), (4.41b)
4.3. RELATIVISTIC BINARY COLLISION MODEL 71
with K1 (z) denoting the modified Bessel function. With regard to our subsequent discus-
sion, we are interested in calculating the mean drift force g, defined by15
δP (t)
g(P ) := . (4.42)
δt b
where nb = N/L is the number density of the bath particles. In order to determine g(P ),
we note that for some arbitrary function G(p, P ), we have
Z ∞
(m2 + p2r )1/2
LIr (t, δt) −1
G(pr , P ) = (2ZJ ) dp G(pr , P ) exp − ×
δt b −∞ kB T
pr P
2 2 1/2
− . (4.44)
(m + pr ) (M + P 2)1/2
2
The first term on the rhs. of Eq. (4.43) involves the function
ǫ(p)
G1 (p, P ) := 2γ(ur )2 , (4.45a)
E(P ) + ǫ(p)
The considerations from Chapter 3 imply that the stationary momentum distribution of
the Brownian particle in the binary collision model is given by the Jüttner function
1
m/M = 0.01
hδP/δtib / (nb M c2 )
0.5
0
kB T/(M c2 ) = 0.1
−0.5 kB T/(M c2 ) = 0.01
kB T/(M c2 ) = 0.001
−1
−4 −2 0 2 4
P / (M c)
Figure 4.3: Relativistic binary collision model. Mean drift force g(P ) := hδP (t)/δtib numerically
evaluated from Eq. (4.43) for different values kB T, with nb = N/L denoting the number density
of the heat bath particles.
Hence, in order for Eq. (4.46) to yield the correct stationary distribution, the functions α
and D must be coupled by the relativistic Einstein relation (4.2f), reading
where E = (p2 + M 2 )1/2 . In order to determine the function α, we demand that the
Langevin equation yields the same mean drift force g as the collision model, i.e.,
dP (t) ! δP (t)
P (t) = p = P (t) = p . (4.49)
dt δt b
For the post-point (backward-Ito) Langevin equation (4.46) we know that [cf. Eq. (C.25)]
dP (t) d
P (t) = p = −α(p) p + D(p). (4.50)
dt dp
Thus, by means of the Einstein relation (4.48), the lhs. of Eq. (4.49) is given by
dP (t) d
P (t) = p = −α(p) p + β −1 [α(p) E(p)], (4.51)
dt dp
and the condition (4.49) becomes equivalent to the differential equation [cf. Eq. (2.59)]
d
−α(p) p + β −1 [α(p) E(p)] = g(p). (4.52)
dp
4.3. RELATIVISTIC BINARY COLLISION MODEL 73
In the case of the collision model, where the function g(p) is not exactly known, one could,
e.g., try to fit g(p) by a simple analytic expression and, subsequently, use this approximation
in Eq. (4.52).16
16
Alternatively, one can also try to derive a Fokker-Planck equation from, e.g., a relativistic Boltzmann
equation [278, 294].
74 CHAPTER 4. RELATIVISTIC BROWNIAN MOTION
Chapter 5
The preceding chapter was dedicated to relativistic Brownian motions in phase space. In
the remainder we will discuss relativistic diffusion models in Minkowski space-time, i.e.,
continuous relativistic processes that do not explicitly depend on the momentum coor-
dinate. On the one hand, such space-time processes may be constructed, e.g., from a
Brownian motion processes in phase space by integrating out the momentum coordinates.
As a result of this averaging procedure, the reduced process for the position coordinate will
be non-Markovian. Alternatively, one can try to derive or postulate a relativistic diffusion
equation and/or diffusion propagators by means of microscopic models [140, 375, 377, 381]
or plausibility considerations [15]. Regardless of the approach adopted, in order to comply
with the principles of special relativity, the resulting space-time process must be non-
Markovian, as rigorously proven by Dudley (Theorem 11.3 in [341]) and Hakim (Propo-
sition 2 in [346]). Put differently, any relativistic generalization of the nonrelativistic
diffusion equation (1.1) with constant coefficients should be of at least second order in the
time coordinate.
The construction and analysis of relativistic diffusion models in Minkowski space-time
poses an interesting problem in its own right. In addition, the investigation of these
processes becomes relevant in view of potential analogies with relativistic quantum the-
ory [329, 398], similar to the analogy between Schrödinger’s equation and the diffusion
equation (1.1) in the nonrelativistic case [474, 475]. The present section intends to provide
an overview over classical relativistic diffusion models that have been discussed in the lit-
erature [15, 329, 345, 375–377, 381, 390–392, 395, 396]. For this purpose, we first recall basic
properties of the Wiener (Gaussian) process, which constitutes the standard paradigm for
nonrelativistic diffusions in position space (Section 5.1). Subsequently, relativistic general-
izations of the nonrelativistic diffusion equation (1.1) and/or the nonrelativistic Gaussian
diffusion propagator will be discussed [15].
75
76 CHAPTER 5. NON-MARKOVIAN RELATIVISTIC DIFFUSION
∂
̺(t, x) = D ∇2 ̺(t, x), t ≥ 0, (5.1)
∂t
where D > 0 is the spatial diffusion constant, and ∇2 = ∂ 2 /∂x2 in the 1D case. It is well-
known for a long time that the diffusion equation (5.1) is in conflict with the postulates of
special relativity. A simple way of seeing this is to consider the propagator of Eq. (5.1),
which for d = 1 space dimensions is given by
1/2
(x − x0 )2
1
p(t, x|t0 , x0 ) = exp − , t > t0 . (5.2)
4π D(t − t0 ) 4D(t − t0 )
The propagator (5.2) represents the solution of Eq. (5.1) for the initial condition
̺(t0 , x) = δ(x − x0 ).
That is, if X(t) denotes the random path of a particle with fixed initial position X(t0 ) = x0 ,
then p(t, x|t0 , x0 )dx gives the probability that the particle is found in the infinitesimal
volume element [x, x + dx] at time t > t0 . As evident from Eq. (5.2), for each t > t0 there
is a small, but non-vanishing probability that the particle may be observed at distances
|x − x0 | > c(t − t0 ), where c = 1 is the speed of light in natural units. The evolution of the
nonrelativistic Gaussian PDF from Eq. (5.2) is depicted in Fig. 5.1.
0.3
t=1
̺(t, x) / (D−1 c)
t=2
t=4
0.2 t = 10
0.1
0
−10 −5 0 5 10
−1
x / (Dc )
Figure 5.1: Spreading of the Gaussian PDF ̺(t, x) = p(t, x|0, 0) from Eq. (5.2) at different
times t, where t is measured in units of D/c2 . At initial time t = t0 = 0, the PDF corresponds to
a δ-function centered at the origin.
5.1. REMINDER: NONRELATIVISTIC DIFFUSION EQUATION 77
We next summarize basic properties of Eqs. (5.1) and (5.2): Equation (5.1) is a linear
parabolic partial differential equation. Due to the linearity, more general solutions may
be constructed by superpositioning, i.e., by integrating the solution (5.2) over some given
initial PDF ̺0 (x0 ). Equation (5.1) describes a Markov process which means that the
transition PDF (5.2) satisfies the Chapman-Kolmogoroff criterion
Z
p(t, x|t0 , x0 ) = dx1 p(t, x|t1 , x1 ) p(t1 , x1 |t0 , x0 ) (5.3)
R
for all t1 ∈ (t0 , t). The corresponding diffusion process X(t) can be characterized in terms
of the following SDE:
where B(t) is a standard Wiener process as defined in Section 2.1.1. Formally, Eq. (5.4)
may be obtained from the Langevin equation (2.3b) of the classical Ornstein-Uhlenbeck
process as follows: First we rewrite Eq. (2.3b) as
1/2
dV (t) 2D
= −V dt + ∗ dB(t). (5.5)
Mα M 2 α2
Upon letting (Mα) → ∞ and D → ∞ such that D = D/(αM)2 remains constant, the
lhs. of Eq. (5.5) should become negligible. Then, by making use of dX = V dt, Eq. (5.4) is
recovered. This limiting procedure defines the so-called overdamped regime of the Ornstein-
Uhlenbeck process. The mean square displacement of the overdamped process (5.4) is given
by [220]
Z
2
[X(t) − X(t0 )] := dx (x − x0 )2 p(t, x|t0 , x0 )
= 2D (t − t0 ), (5.6)
qualitatively similar to that of classical Ornstein-Uhlenbeck process; cf. Eq. (2.8). Finally,
we note that the solution of Eq. (5.1) with initial condition
where h · i indicates an average with respect to the standard Wiener process B(t) with
initial condition B(t0 ) = 0. Equation (5.8) yields an efficient Monte-Carlo simulation
scheme for computing the solutions of the diffusion equation (5.1).
78 CHAPTER 5. NON-MARKOVIAN RELATIVISTIC DIFFUSION
∂2
∂
τv 2 + ̺(t, x) = D ∇2 ̺(t, x). (5.9)
∂t ∂t
Here, D > 0 plays again the role of a diffusion constant, while τv > 0 is an additional
relaxation time parameter. Similar to Eq. (5.1), the telegraph equation (5.9) refers to a
special frame where the background medium, causing the random motion of the diffusing
test particle, is at rest (on average). The ‘nonrelativistic limit’ corresponds to letting τv → 0
in Eq. (5.9), which leads back to Eq. (5.1). For τv > 0, Eq. (5.9) is a linear hyperbolic
partial differential equation. Because of the second order time derivative in Eq. (5.9), one
now also has to specify the first order time derivative of the initial distribution at time t0 .
Considering particular initial conditions
∂
̺(t0 , x) = δ(x − x0 ), ̺(t0 , x) ≡ 0, (5.10)
∂t
one finds that the corresponding solution of Eq. (5.9) is given by [375, 377]
e−(t−t0 )/(2τv )
p(t, x|t0 , x0 ) = δ[|x − x0 | − v(t − t0 )] +
2
Θ(ξ 2 )
t I1 (ξ)
I0 (ξ) + . (5.11a)
2τv v 2τv ξ
Here, we have abbreviated
" 2 2 #1/2
1 t − t0 x − x0
ξ := − , v := (D/τv )1/2 , (5.11b)
2 τv τv v
and the modified Bessel functions of the first kind, Iν (z), are defined by
∞
X 1 z 2k+ν
Iν (z) :=
k=0
Γ(k + ν + 1) k! 2
5.2. TELEGRAPH EQUATION 79
with Γ(z) denoting the Euler gamma function. According to our knowledge, the solu-
tion (5.11) was first obtained by Goldstein in 1938/1939. Actually, Goldstein derived the
result (5.11) by considering the continuum limit of a persistent random walk model [379];
subsequently, he proved that this function satisfies the telegraph equation (5.9), cf. Sec-
tion 8 of his paper [375].
The propagator (5.11) is characterized by two salient features:
a) As evident from the δ-function term, the solution exhibits two singular diffusion
fronts traveling at absolute velocity v := (D/τv )1/2 to the left and right, respectively;
b) due to the appearance of the Heaviside Θ-function, the solution is non-zero only
within the region |x − x0 | ≤ v(t − t0 ), i.e., upon fixing τv such that v = c = 1 the
solution vanishes outside the light cone.
Thus, in contrast to the nonrelativistic propagator (5.2), Eqs. (5.11) define a relativisti-
cally acceptable diffusion model. Because of the second order time derivative, the telegraph
equation (5.9) describes a non-Markovian process, in accordance with the aforementioned
theorems of Dudley [341] and Hakim [346]. The non-Markovian character of the propa-
gator (5.11) can also be proven directly by verifying that this solution does not fulfill the
condition (5.3).
The linearity of Eq. (5.9) implies that more general solutions can be obtained by in-
tegrating the propagator (5.11) over some given initial distribution ̺0 (x0 ). In principle,
one may also construct other classes of solutions with ∂̺(t0 , x)/∂t 6≡ 0, e.g., by applying
a Laplace-Fourier transformation [377, 478, 479] to Eq. (5.9). We note, however, that in
order for the solution ̺(t, x) to remain normalized at all times t > t0 , it is required that
Z +∞
∂
dx ̺(t0 , x) ≡ 0, (5.12a)
−∞ ∂t
and
∂
̺(t, ±∞) ≡ 0. (5.12b)
∂x
In the case of Eq. (5.10), these conditions are automatically satisfied. Various solutions and
extensions of the telegraph equation (5.9), including different types of boundary conditions,
additional external sources, etc., have been discussed, e.g., by Goldstein [375], Masoliver
et al. [478, 479], Foong and Kanno [480], Renardy [481], and Dorogovtsev [482].
Similar to the nonrelativistic diffusion equation (5.1), the telegraph equation (5.9) may
be derived and/or motivated in many different ways. A detailed overview is given by
Masoliver and Weiss [377], who discuss four different possibilities of deducing Eq. (5.9)
from underlying models; see also Koide [395, 396]. During the past decades, the telegraph
equation (5.9) has been used to describe a number of different phenomena. The applications
include:
80 CHAPTER 5. NON-MARKOVIAN RELATIVISTIC DIFFUSION
• Transmission of electrical signals. According to Masoliver and Weiss (5.9), the earli-
est derivation of the telegraph equation is based on a paper by William Thomson [483]
(who later became Lord Kelvin), published in 1855. He considered the problem of
how to transmit electrical signals without distortion, a question closely related to the
design of the first transatlantic cable.
• Heat transport and propagation of heat waves. In this case, the function ̺(t, x) in
Eq. (5.9) is interpreted as a temperature field and the normalization condition is
usually dropped. For a detailed account of the vast literature on heat waves we refer
to the review article of Joseph and Preziosi [388, 389].
• High energy ion collision experiments. In recent years, the telegraph equation has
been used to estimate the dissipation of net charge fluctuations, which may obliterate
signals of QCD phase transitions in nuclear collisions [393, 394]. In this context,
however, the coordinate x in Eq. (5.9) is interpreted as a rapidity variable.
Another interesting aspect of the telegraph equation is elucidated in a paper by Kac [376].
He observed that the solutions of Eq. (5.9) with initial conditions
∂
̺(t0 , x) ≡ ̺0 (x), ̺(t0 , x) ≡ 0, (5.13)
∂t
1
Boguñá et al. [381] discuss persistent random walks in higher space dimensions.
5.2. TELEGRAPH EQUATION 81
where v = (D/τv )1/2 , and h · i indicates an average with respect to the τv -parameterized
Poisson process N(t); i.e., for any given time t > t0 we have
k
e−(t−t0 )/(2τv )
t − t0
Prob{N(t) = k} = , k = 0, 1, 2, . . . ; (5.15a)
k! 2τv
and for any finite sequence t0 < t1 < . . . < tn the increments
are independent. Equation (5.14) is the direct counterpart of Feynman-Kac formula (5.8).
Equations (5.14) and (5.15) together provide a very efficient Monte-Carlo simulation
scheme for computing solutions of the telegraph equation (5.9). Moreover, the Poisson
path integral representation (5.14) discloses an interesting correspondence between the
free-particle Dirac equation [397] and the telegraph equation (5.9), which was first pointed
out by Gaveau et al. [398] in 1984: The solutions of both equations may be linked by means
of an analytic continuation. This connection is similar to the relation between the diffusion
equation (5.1) and the free particle Schrödinger equation in the nonrelativistic case.3 The
crucial difference is given by the fact that the measures of the functional integration refer
to different underlying processes, respectively.
However, the telegraph equation (5.9) is not the only possible relativistic generalization of
the nonrelativistic diffusion equation (1.1) and, recently, there has been some controversy
about its applicability and validity [390–392]. An early critical discussion of Eq. (5.9) in
the context of relativistic heat transport was given by van Kampen [140] in 1970. Starting
from a simple microscopic model, consisting of a cloud of material particles that exchange
electromagnetic radiation, van Kampen derived an integral equation for the temperature
of the particles as function of time and space. He then showed how the telegraph equa-
tion (5.9) can be recovered as an approximation to the more precise integral equation, but
that the validity of this approximation breaks down in the vicinity of the diffusion fronts.
2
The result (5.14) may be generalized to an arbitrary number of space dimensions; cf. pp. 500 in Kac’s
paper [376].
3
For further reading about path integral representations of the Dirac propagator we refer to the papers of
Ichinose [399,400], Jacobson and Schulman [401], Barut and Duru [402], and Gaveau and Schulman [403];
see also footnote 7 in Gaveau et al. [398] and problem 2-6, pp. 34-36 in Feynman and Hibbs [404].
82 CHAPTER 5. NON-MARKOVIAN RELATIVISTIC DIFFUSION
Similarly, the singular diffusion fronts predicted by Eq. (5.11) represent a source of con-
cern if one wishes to adopt the telegraph equation (5.9) as a model for particle transport
in a random medium. While these singularities may be acceptable in the case of photon
diffusion [382–386], they seem unrealistic for massive particles, because such fronts would
imply that a finite fraction of particles carries a huge amount of kinetic energy (much
larger than mc2 ). In view of these shortcomings, it appears reasonable to explore other
constructions of relativistic diffusion processes [15, 405]. In the next part we will propose
a different approach [15] that may provide a viable alternative to the solutions of the
telegraph equation.
1 t ′
Z
a(x̄|x̄0 ) = dt v(t′ )2 , (5.16)
2 t0
where the velocity v(t′ ) is a piecewise constant function, satisfying
Z t
x = x0 + dt′ v(t′ ). (5.17)
t0
Clearly, the nonrelativistic action (5.16) becomes minimal for the deterministic (direct)
path, i.e., if the particle does not collide at all. In this case, it moves with constant
velocity v(t′ ) ≡ (x − x0 )/(t − t0 ) for all t′ ∈ [t0 , t], yielding the smallest possible action
value
(x − x0 )2
a− (x̄|x̄0 ) = . (5.18)
2(t − t0 )
On the other hand, to match the boundary conditions it is merely required that the mean
velocity equals (x − x0 )/(t − t0 ). Consequently, in the nonrelativistic case, the absolute
5.3. RELATIVISTIC DIFFUSION PROPAGATOR 83
velocity of a particle may become arbitrarily large during some intermediate time inter-
val [t′ , t′′ ] ⊂ [t0 , t]. Hence, the largest possible action value is a+ (x̄, x̄0 ) = +∞. These
considerations put us in the position to rewrite the Wiener propagator (5.2) as an integral-
over-actions:
Z a+ (x̄|x̄0 )
a
p(x̄|x̄0 ) ∝ da exp − , (5.19a)
a− (x̄|x̄0 ) 2D
Analogous to the nonrelativistic case, the relativistic action (5.20) assumes its minimum
a− for the deterministic (direct) path from x0 to x, characterized by a constant velocity
v(t′ ) ≡ (x − x0 )/(t − t0 ). One explicitly obtains
1/2
a− (x̄, x̄0 ) = − (t − t0 )2 − (x − x0 )2
, (5.21a)
i.e., a− is the negative Minkowski distance of the two space-time events x̄0 and x̄. The
maximum action value is realized for particles moving at light speed, yielding a+ = 0.
Hence, the transition PDF for the relativistic generalization of the Wiener process reads
−1 a− (x̄, x̄0 )
p(x̄|x̄0 ) = N exp − −1 , (5.21b)
2D
0.6 t=1
̺(t, x) / (D−1 c)
t=2
t=4
0.4 t = 10
0.2
0
−10 −5 0 5 10
x / (Dc−1 )
Figure 5.2: Transition PDF ̺(t, x) = p(t, x|0, 0) for the one-dimensional (d = 1) relativistic
diffusion process (5.21) at different times t (measured in units of D/c2 ). At time t = t0 = 0,
the function ̺(t, x) reduces to a δ-function centered at x0 = 0. In contrast to the nonrelativistic
diffusion propagator, cf. Fig. 5.1, the PDF (5.21) vanishes outside of the light cone.
It is interesting to note that the PDF (5.19) is a special case of a larger class of diffusion
processes, defined by
Z a+ (x̄|x̄0 )
pw (x̄|x̄0 ) = Nw−1 da w(a), (5.22)
a− (x̄|x̄0 )
This can be achieved, e.g., by choosing the integral boundaries as ã− = (x − x0 )2 /(t − t0 )α , α 6= 1 and
4
a+ = ∞, but then the variable a may not be interpreted as a conventional action anymore.
5.3. RELATIVISTIC DIFFUSION PROPAGATOR 85
1
hX 2 (t)i / (2Dt)
0.75
0.5
0.25 nonrelativistic
relativistic
0
0 25 50 75 100
t / (Dc−2 )
Figure 5.3: Comparison of the mean square displacements X 2 (t) , divided by 2Dt, for the one-
dimensional (d = 1) nonrelativistic Wiener process (5.2) and its relativistic generalization from
Eq. (5.21) with initial condition (t0 , x0 ) = (0, 0).
86 CHAPTER 5. NON-MARKOVIAN RELATIVISTIC DIFFUSION
Chapter 6
In this thesis, we have investigated how the concepts of Brownian motion and diffusion [4]
may be incorporated into the framework of special relativity [2, 3]. As the starting point,
we summarized in Chapter 2 the foundations of the nonrelativistic Langevin theory [31,32]
of Brownian motions. Langevin equations are stochastic differential equations (SDEs),
and provide a simplified description of the complex microscopic interactions between a
Brownian particle and its environment [64,65]. Nonrelativistic stochastic processes as, e.g.,
the classical Ornstein-Uhlenbeck process [36, 39] can be used as guidance for constructing
relativistic Langevin equations. In particular, these nonrelativistic processes should also be
recovered from the relativistic theory in the limit of an infinite light speed c → ∞. Since, in
the past, nonrelativistic SDEs have proven to be very useful for describing a wide range of
different phenomena [66, 79, 336], relativistic Langevin equations can be expected to play
a similarly successful role in the modeling of thermalization and relaxation processes in
relativistic and astrophysical system [312–314, 337, 338, 471].
When adopting a phenomenological or axiomatic Langevin equation approach, the fluctu-
ation-dissipation relations of the friction and noise coefficients have to be chosen such that
the correct equilibrium distributions are recovered in the stationary state. For example,
in the nonrelativistic case it is commonly accepted that the stationary velocity distribu-
tion of a free1 Brownian particle in a thermally equilibrated environment (heat bath) is
given by a Maxwell distribution [168]. As well-known, Maxwell’s distribution is in conflict2
with special relativity and, therefore, has to be replaced in a relativistic theory. In recent
years, there has been considerable debate in the literature about the correct generalization
of Maxwell’s distribution in special relativity [14, 206, 208, 209, 211, 215–217]. In order to
elucidate the differences between two of the most commonly considered candidate distri-
butions [14, 167, 211], we considered a maximum principle for the relative entropy [16],
1
By ‘free’ we mean the absence of external fields.
2
The Maxwell distribution, representing a Gaussian velocity distribution, assigns a non-vanishing prob-
ability to superluminal velocities.
87
88 CHAPTER 6. SUMMARY AND OUTLOOK
effects. On the other hand, it would be desirable to identify classes of physical sys-
tems that can be described by these processes. For example, it seems likely that that
the quasi-random particle acceleration [495] in various astrophysical systems may be
efficiently modeled by non-Gaussian driving processes.
This appendix summarizes basic definitions and a few properties of Lorentz transforma-
tions. For a more detailed introduction to special relativity we refer to Weinberg [6] or
Sexl and Urbantke [8].
By definition, the components of the covariant vector (xα ) are obtained by contracting the
contravariant vector (xα ) with ηαβ , i.e.,1
d
X
xα := ηαβ xβ =: ηαβ xβ , (A.2)
β=0
1
Sometimes, it is convenient to interpret the contravariant vectors (aα ) as column vectors and the
corresponding covariant vectors (aα ) as row vectors.
93
94 APPENDIX A. SPECIAL RELATIVITY (BASICS)
yielding explicitly
The rhs. of Eqs. (A.2) defines the Einstein’s summation convention. The vectors (xα )
and (xα ) will be called four-vectors, regardless of the number of space dimensions. The
corresponding four-vector gradients are defined by
∂ ∂ ∂ ∂ ∂
(∂α ) := = , ,..., d = ,∇ , (A.4a)
∂xα ∂t ∂x1 ∂x ∂t
α ∂ ∂ ∂ ∂ ∂
(∂ ) := = − , 1,..., d = − ,∇ . (A.4b)
∂xα ∂t ∂x ∂x ∂t
The components η αβ of the inverse metric tensor are determined by the condition
!
xα = η αβ xβ = η αβ ηβγ xγ ∀ (xα ), (A.5a)
or, equivalently, by
!
η αβ ηβγ = δ α γ , (A.5b)
η αβ = ηαβ . (A.6)
The Minkowski space-time distance between two events x̄A = (xαA ) = (tA , xA ) and
x̄B = (xαB ) = (tB , xB ) is defined by
In special relativity, events with time-like separation can be causally connected by (a series
of) signals travelling slower or equal to the speed of light. Events with light-like separation
can be causally related only by undisturbed signals travelling at the speed of light. Events
with space-like separation are causally disconnected.
A.2. LORENTZ-POINCARÉ TRANSFORMATIONS 95
uα uα = −1. (A.9b)
where m > 0 is the rest mass of the particle. Upon comparing with (A.8), one finds for a
particle with velocity v in Σ
x′α = Λα β xβ + aα , (A.12a)
mapping the ‘old’ Σ-space-time-coordinates xα onto the ‘new’ Σ′ -coordinates x′α . The
constant four-vector aα shifts the origins of time and space, while the constant Lorentz
matrix (Λα β ) may account for a spatial rotation, a change of orientation and/or a relative
velocity between the two frames Σ and Σ′ , cf. Chapter 1 of Weinberg’s book [6]. The
matrix components Λα β are determined by the condition
! !
d(x̄′A , x̄′B )2 = d(x̄A , x̄B )2 ⇔ Λα γ Λβ δ ηαβ = ηγδ . (A.12b)
96 APPENDIX A. SPECIAL RELATIVITY (BASICS)
The condition (A.12b) means that causal relations remain preserved during transitions
between inertial systems.
The LPTs (A.12) form a group.2 Of particular interest for our purpose, is the subgroup
of the proper LPTs, defined by imposing the additional constraints
The requirements (A.13) exclude time reversal and space inversion. Examples are pure
rotations
where (Rij ) is a rotation matrix [i.e., det(Rij ) = 1 and Rij Rkj = δij ], and Lorentz boosts [6]
(γ − 1)
Λ0 0 = γ, Λi 0 = Λ0 i = γw i , Λi j = δ i j + w i w j . (A.15)
w2
with velocity w = (w 1 , . . . , w d ) and Lorentz factor
γ := (1 − w 2 )−1/2 . (A.16)
To illustrate the effect of a boost, we first consider a particle that rests at the spatial
origin of Σ and, therefore, being described by (xα ) = (t, 0) in Σ. By applying the Lorentz
boost (A.15) to (xα ) = (t, 0), we find
which means that the particle moves at constant velocity w through Σ′ . Similarly, we may
consider a particle that moves at velocity −w through Σ. If the particle at Σ-time t = 0
was located at x = 0 then it is described by the world-line (xα ) = (t, −tw) in Σ. By
applying the Lorentz boost (A.15) to (xα ) = (t, −tw), we obtain
p′α = Λα β pβ . (A.19)
2
A detailed discussion can be found in [6, 8, 9].
A.2. LORENTZ-POINCARÉ TRANSFORMATIONS 97
Combining Eqs. (A.19) and (A.12b), one can verify the mass-shell condition
which means that the rest mass m is a Lorentz invariant. In particular, the mass shell
condition (A.20) implies that Eq. (A.19) is equivalent to the nonlinear restricted transfor-
mation p 7→ p′ , given by
and we find
∂p′i pk
k
= Λi 0 2 + Λi k . (A.21b)
∂p (m + p2 )1/2
98 APPENDIX A. SPECIAL RELATIVITY (BASICS)
Appendix B
Normalization constants
Z1 = 2m K1 (βm), (B.3a)
K2 (βm)
Z3 = 4πm3 , (B.3b)
βm
with Kn (z) denoting the modified Bessel functions of the second kind [466]. The energy
mean values are obtained by logarithmic differentiation,
∂
hEid = − ln Zd , (B.4)
∂β
K0 (βm) + K2 (βm)
hEi1 = m , (B.5a)
2K1 (βm)
3 K1 (βm)
hEi3 = +m . (B.5b)
β K2 (βm)
99
100 APPENDIX B. NORMALIZATION CONSTANTS
in terms of modified Bessel functions of the first kind and modified Struve functions [466].
Introducing z := x − x0 and u := t − t0 , we have to calculate
2
(u − z 2 )1/2
Z
d
Nd = d z Θ(u − |z|) exp −1 .
Rd 2D
where Od = 2π d/2 /Γ(d/2) is the surface area of the d-dimensional unit-sphere. It is conve-
nient to split the integral (B.7) in the form
ud
Nd = Nd′ − Od , (B.8)
d
where
u
(u2 − |z|2 )1/2
Z
Nd′ = Od d|z| |z| d−1
exp , (B.9)
0 2D
Next we substitute |z| = u sin φ, where φ ∈ [0, π/2]. Then Eq. (B.9) takes the form
π/2
u cos φ
Z
Nd′ d
= u Od dφ cos φ sin d−1
φ exp .
0 2D
Based on this integral representation, Nd′ can be expressed in terms of modified Bessel
functions of the first kind In , and modified Struve functions Lk [466], and one finds
where χ = u/(2D).
Appendix C
In this appendix we briefly summarize the most commonly considered stochastic integral
definitions and the corresponding rules of stochastic calculus. For a more rigorous and
more comprehensive introduction, we refer to Refs. [63, 64, 66, 407].
We consider a Wiener process (standard Brownian motion) B(t) as defined in Sec-
tion 2.1.1; i.e., the increments dB(t) := B(t + dt) − B(t) are stochastically indepen-
dent [64, 66] and characterized by the Gaussian PDF
where f (y) is some real-valued function and Y (s) a real-valued time-dependent process. If
B(s) were some ordinary differentiable function of s ∈ [0, t], then the integral in Eq. (C.2)
would simply be given by1
Z t
I= f (Y (s)) Ḃ(s) ds, (C.3)
0
where Ḃ = dB/ds. Unfortunately, Ḃ(s) is not well-defined for the Wiener process [64,407],
but it is possible to generalize the concept of integration to also include the Wiener process
and other stochastic processes [64,66,407]. However, in contrast to the standard Riemann-
Stieltjes integral (C.3), the integral with respect to a stochastic process may depend on
the choice of the discretization scheme and, in particular, also require modifications of
differential calculus.
1
By writing Eq. (C.3), it is implictly assumed that f (y), Y and Ḃ are sufficiently smooth functions so
that this integral exists in the sense of Riemann-Stieltjes; in this case, the value of the integral (C.3) is
independent of the underlying discretization scheme [66].
101
102 APPENDIX C. STOCHASTIC INTEGRALS AND CALCULUS
To illustrate these aspects for the most commonly considered stochastic integral defini-
tions, we will always consider the following equidistant partition {t0 , t1 . . . , tN } of the time
interval [0, t]:
where the partition {t0 , . . . , tN } is given by (C.4). The peculiar, defining feature of this
integral is that, on the rhs. of Eq. (C.5), the argument of the function f must be evaluated
at the lower boundary points tk of the discrete intervals [tk , tk+1 ]; i.e., the definition of the
Ito integral is non-anticipating. Accordingly, the Ito discretization scheme is also known
as the pre-point rule.
Now consider a stochastic process Y (t) which, for two given functions A(y) and C(y), is
defined by
Z t Z t
Y (t) = Y (0) + A(Y (s)) ds + C(Y (s)) ∗ dB(s), (C.6)
0 0
and where the last term is interpreted as an Ito integral (C.5). Stochastic integral equations
like Eq. (C.6) are usually abbreviated by rewriting them as an Ito stochastic differential
equation (I-SDE)
complemented by the initial condition Y (0). From the non-anticipating definition (C.5) of
the Ito integral and the properties of the Wiener process it follows that [407]2
The Fokker-Planck equation for the PDF f (t, y) of the stochastic process defined by
Eq. (C.7) reads
∂f ∂ 1 ∂ 2
= −Af + (C f ) , (C.9)
∂t ∂y 2 ∂y
where A = A(y) and C = C(y). The deterministic initial condition Y (0) = y0 translates
into f (0, y) = δ(y − y0 ).
Finally, an important pecularity arises when one considers nonlinear transformations G of
the stochastic process Y (t). More precisley, assuming that Y is defined by the I-SDE (C.7),
then the differential change of the process Z(t) := G(Y (t)) is given by (see, e.g., Section
4.3.2 in [407])
1
dZ(t) = G′ (Y ) ∗ dY + C(Y )2 G′′ (Y ) dt
2
1
′
= A(Y ) G (Y ) + C(Y ) G (Y ) dt + C(Y ) G′ (Y ) ∗ dB(t),
2 ′′
2
(C.10)
where G′ (y) = dG(y)/dy and G′′ (y) = d2 G(y)/dy 2. Within ordinary differential calculus,
the term containing G′′ is absent. Equation (C.10) is usually referred to as Ito formula.
where ∂i := ∂/∂y i .
104 APPENDIX C. STOCHASTIC INTEGRALS AND CALCULUS
t N −1
1
Z X
f (Y (s)) ◦ dB(s) := lim [f (Y (tk+1 )) + f (Y (tk ))] ×
0 N →∞
k=0
2
[B(tk+1 ) − B(tk )] (C.14)
where the partition {t0 , . . . , tN } is given by (C.4). In contrast to Ito’s integral (C.5), the
SF definition (C.14) uses the mean of the boundary values of f on the intervals [tk , tk+1 ];
i.e., the definition of the SF integral is semi-anticipating. This discretization scheme is also
known as the mid-point rule.
Similar to Eq. (C.6), we may consider a stochastic process Y (t) defined by
Z t Z t
Y (t) = Y (0) + A(Y (s)) ds + C(Y (s)) ◦ dB(s), (C.15)
0 0
where now the last term is interpreted as an SF integral (C.14). The integral equa-
tion (C.15) can be abbreviated in terms of the equivalent SF stochastic differential equation
(SF-SDE)
with initial condition Y (0). From the semi-anticipating definition (C.14) of the SF integral
and the properties of the Wiener process it follows that [407]
1
hC(Y ) ◦ dB(t) | Y (t) = yi = C(y) C ′(y) dt, (C.17)
2
where C ′ (y) = dC(y)/dy. The Fokker-Planck equation for the PDF f (t, y) of the stochastic
process (C.16) reads
∂f ∂ 1 ∂
= −Af + C (Cf ) (C.18)
∂t ∂y 2 ∂y
C.2. STRATONOVICH-FISK INTEGRAL 105
where A = A(y), C = C(y) and C ′ = C ′ (y). The deterministic initial condition Y (0) = y0
translates into f (0, y) = δ(y − y0 ).
It can be shown [66, 407] that the SF integral definition preserves the rules of ordinary
stochastic calculus; i.e., if Y (t) is defined by the SF-SDE (C.16), then the differential
change of the process Z(t) := G(Y (t)) is given by (see, e.g., Section 4.3.2 in [407])
dZ(t) = G′ (Y ) ◦ dY
= A(Y ) G′ (Y ) dt + C(Y ) G′ (Y ) ◦ dB(t),
(C.19)
Consider the n-dimensional stochastic process Y (t) = (Y 1 (t), . . . , Y n (t)), defined by the
following n-dimensional generalization of Eq. (C.16):
∂f ∂ i 1 i ∂ j
= i −A f + C r j (C r f ) , (C.21)
∂t ∂y 2 ∂y
where the partition {t0 , . . . , tN } is given by (C.4). On the rhs. of Eq. (C.22), in contrast
to the Ito and SF integrals, the argument of the function f must be evaluated at the upper
boundary points tk+1 of the discrete intervals [tk , tk+1 ]; i.e., the definition of this integral
is anticipating. This discretization scheme is also known as the post-point rule.
Similar to above, we may consider a stochastic process Y (t) which, for two given functions
A(y) and C(y), is defined by
Z t Z t
Y (t) = Y (0) + A(Y (s)) ds + C(Y (s)) • dB(s), (C.23)
0 0
and where the last term is now interpreted as a BI integral (C.22). Equation (C.6) can be
abbreviated by rewriting it as a backward Ito stochastic differential equation (BI-SDE)
complemented by the deterministic initial condition Y (0). From the anticipating defini-
tion (C.22) of the BI integral and the properties of the Wiener process it follows that [407]
The Fokker-Planck equation for the PDF f (t, y) of the stochastic process defined by
Eq. (C.24) reads
∂f ∂ 1 2∂
= −Af + C f , (C.26)
∂t ∂y 2 ∂y
where A = A(y) and C = C(y). The deterministic initial condition Y (0) = y0 translates
into f (0, y) = δ(y − y0 ).
C.4. COMPARISON OF STOCHASTIC INTEGRALS 107
It can be shown that, similar to the Ito integral, also the BI integral requires a modification
of differential calculus. More precisley, assuming that Y is defined by the BI-SDE (C.24),
the differential change of the process Z(t) := G(Y (t)) is given by
1
dZ(t) = G′ (Y ) • dY − C(Y )2 G′′ (Y ) dt
2
1
= A(Y ) G (Y ) − C(Y ) G (Y ) dt + C(Y ) G′ (Y ) • dB(t),
′ 2 ′′
2
(C.27)
Thus, upon comparing Eqs. (C.37) and (C.38) with Eq. (C.36), we see that Eqs. (C.31a)
and (C.31c) describe the same process like Eq. (C.31b), if we fix
1 j 1 i
Ai∗ = Ai◦ + C r ∂j C i r , Ai• = Ai◦ − C r ∂j C j r . (C.39)
2 2
1 i j
∂t f = ∂i −Ai• f + C r C r ∂j f . (C.40)
2
Next we rewrite the Ito Fokker-Planck equation (C.12), corresponding to Eq. (C.31a), as
i 1 i j 1 i j
∂t f = ∂i − A∗ − ∂j (C r C r ) f + C r C r ∂j f , (C.41)
2 2
Similarly, the SF Fokker-Planck equation (C.21), corresponding to Eq. (C.31c), can also
be written as
i 1 i j 1 i j
∂t f = ∂i − A◦ − C r ∂j C r f + C r C r ∂j f . (C.42)
2 2
Thus, upon comparing Eqs. (C.41) and (C.42) with Eq. (C.40), we see that Eqs. (C.31a)
and (C.31b) describe the same process like Eq. (C.31c), if we fix
1 1 i
Ai∗ = Ai• + ∂j (C i r C j r ) , Ai◦ = Ai• + C r ∂j C j r . (C.43)
2 2
To summarize, by means of Eqs. (C.35), (C.39) and (C.43) one can change between the
different forms of stochastic integration and stochastic differential calculus, respectively.
Each SDE type has advantages and disadvantages: The Ito formalism is well suited for
numerical simulations [63, 79, 336] and yields a vanishing noise contribution to conditional
expectations of the form (C.8). The Stratonovich-Fisk approach is more difficult to imple-
ment numerically, but preserves the rules of ordinary differential calculus (in contrast to
Ito/backward Ito integration). Finally, within the backward Ito scheme, fluctuation dissi-
pation relations take a particularly elegant form (cf. Sec. 6.2 in Ref. [67], and Ref. [101]).
110 APPENDIX C. STOCHASTIC INTEGRALS AND CALCULUS
Here, the ∆B r (t) are random numbers, sampled from a Gaussian normal distribution with
density
1/2
[∆B r (t)]2
r 1
P[∆B (t)] = exp − . (C.45)
2π∆t 2∆t
As evident from Eqs. (C.44), for given functions Ai∗ and C i r , the discretized Ito SDE (C.44a)
allows for calculating the values Y i (t + ∆t) directly from the preceding values Y i (t). By
contrast, the discretized SF SDEs (C.44b) and BI SDEs (C.44c) are implicit equations,
which must be solved for Y i (t + ∆t). The latter difficulty can be avoided by transforming a
given SF/BI SDE to the corresponding Ito SDE by means of Eqs. (C.35), (C.39) and (C.43).
Appendix D
for i = 1, . . . , d. From Eq. (C.29), the Fokker-Planck equation corresponding to the stochas-
tic process (D.1) is obtained as
pi ∂
∂ ∂ i j 1 i k ∂
+ f = a j p f + c rc r k f (D.3)
∂t p0 ∂xi ∂pi 2 ∂p
with f (t, x, p) denoting the phase space PDF of the relativistic Brownian particle in Σ,
and p0 = (M 2 + p2 )1/2 the relativistic energy.
111
112 APPENDIX D. HIGHER SPACE DIMENSIONS
In general, the friction and noise coefficients, ai j and ci j may be functions of the momen-
tum and position coordinates, P = (P i) and X = (X i ), and of the lab time t. Restricting
the discussion to a heat bath that is stationary, isotropic and homogeneous in the lab frame
Σ, the coefficient matrices take the simplified diagonal form
ai j = α δ i j , ci j = (2D)1/2 δ i j , (D.4a)
where the functions α and D depend only on the Brownian particles’ absolute momentum –
or, equivalently, on its relativistic energy P 0 = (M 2 + P 2 )1/2 . Denoting by U β = (1, 0) the
mean four velocity vector of the heat bath in Σ, we have P 0 = −U β Pβ and may therefore
write
For an isotropic background modelled by Eqs. (D.4), the Langevin equations (D.1) can be
rewritten as
pi ∂
∂ ∂ i ∂
+ f = αp f + D i f , (D.6)
∂t p0 ∂xi ∂pi ∂p
and if the motion is restricted to a finite volume V, then the stationary solution f∞ (x, p)
of Eq. (D.6) is given by a spatially homogeneous Jüttner distribution [167]
and, conversely,
Specifically, for an observer moving at constant velocity w through the lab frame, the
corresponding Lorentz boost matrix elements read [cf. Eq. (A.15)]
(γ − 1)
Λ0 0 = γ, Λi0 = Λ0 i = −γw i , Λi j = δ i j + w i w j . (D.11)
w2
Hence, in order to find f ′ , it suffices to solve the the Fokker-Planck equations (D.3) or (D.6)
in lab frame Σ, e.g., for a given t-simultaneous initial condition f (0, x, p), and to insert
the solution into (D.9a).
114 APPENDIX D. HIGHER SPACE DIMENSIONS
Bibliography
[1] A. Einstein. Über einen die Erzeugung und Verwandlung des Lichtes betreffenden
heuristischen Gesichtspunkt. Ann. Phys. (Leipzig), 17(4):132–148, 1905.
[2] A. Einstein. Zur Elektrodynamik bewegter Körper. Ann. Phys. (Leipzig), 17(4):891–
921, 1905.
[3] A. Einstein. Ist die Trägheit eines Körpers von seinem Energiegehalt abhängig? Ann.
Phys. (Leipzig), 18(4):639–641, 1905.
[4] A. Einstein. Über die von der molekularkinetischen Theorie der Wärme geforderte
Bewegung von in ruhenden Flüssigkeiten suspensierten Teilchen. Ann. Phys.
(Leipzig), 17(4):549–560, 1905.
[6] S. Weinberg. Gravitation and Cosmology. John Wiley & Sons, 1972.
[7] C. Lämmerzahl. Special Relativity and Lorentz Invariance. Ann. Phys. (Leipzig),
14(1-3):71–102, 2005.
[9] S. Weinberg. The Quantum Theory of Fields, volume 1. Cambridge University Press,
Cambridge, 2002.
[10] J. Dunkel and P. Hänggi. Theory of the relativistic Brownian motion: The (1+1)-
dimensional case. Phys. Rev. E, 71:016124, 2005.
[11] J. Dunkel and P. Hänggi. Theory of the relativistic Brownian motion: The (1+3)-
dimensional case. Phys. Rev. E, 72:036106, 2005.
115
116 BIBLIOGRAPHY
[13] J. Dunkel and P. Hänggi. Relativistic Brownian motion: From a microscopic binary
collision model to the Langevin equation. Phys. Rev. E, 74:051106, 2006. Erratum,
Phys. Rev. E 74:069902(E), 2006.
[15] J. Dunkel, P. Talkner, and P. Hänggi. Relativistic diffusion processes and random
walk models. Phys. Rev. D, 75:043001, 2007.
[16] J. Dunkel, P. Talkner, and P. Hänggi. Relative entropy, Haar measures and relativistic
canonical velocity distributions. New J. Phys., 9:144, 2007.
[19] F. Debbasch and J. P. Rivet. A diffusion equation from the relativistic Ornstein-
Uhlenbeck process. J. Stat. Phys, 90:1179–1199, 1998.
[20] F. Debbasch. Relativistic stochastic processes: A review. AIP Conf. Proc., 861:488–
493, 2006.
[23] J. Angst and J. Franchi. Central limit theorem for a class of relativistic diffusions.
J. Math. Phys., 48(8):083101, 2007.
[24] B. Lindner. Exact solution for the diffusion coefficient of nonlinear brownian motion.
New J. Phys., 9:136, 2007.
BIBLIOGRAPHY 117
[27] R. Brown. A brief account of microscopical observations made in the months of June,
July and August, 1827, on the particles contained in the pollen of plants; and on the
general existence of active molecules in organic and inorganic bodies. Philos. Mag,
4:161–173, 1828. reprinted in Edinburgh New Philos. J. 5:358-371, 1928.
[28] W. Sutherland. A dynamical theory of diffusion for non-electrolytes and the molec-
ular mass of albumin. Philos. Mag, 9:781–785, 1905.
[29] A. Einstein and M. von Smoluchowski. Untersuchungen über die Theorie der Brown-
schen Bewegung/Abhandlungen über die Brownsche Bewegung und verwandte Er-
scheinungen, volume 199. Harri Deutsch, Frankfurt, 3 edition, 1999.
[30] J. Perrin. Mouvement brownien et réalité moléculaire. Ann. Chim. Phys., 18:5–114,
1909.
[31] P. Langevin. Sur la théorie du mouvement brownien. Comptes rendus Acad. Sci.
(Paris), 146:530–533, 1908.
[32] D. S. Lemons and A. Gythiel. Paul Langevins’s 1908 paper “On the Theory of
Brownian Motion”. Am J. Phys., 65(11):1079–1081, 1997.
[34] M. Planck. Über einen Satz der statistischen Dynamik und seine Erweiterung in der
Quantentheorie. Sitzungsber. Preuss. Akad. Wiss. Berlin, 1917.
[35] O. Klein. Zur statistischen Theorie der Suspensionen und Lösungen. PhD thesis,
Stockholm Högskola, 1921. published in Ark. Mat. Astron. Fys. 16(5):1-51, 1922.
Almqvist & Wiksell, Uppsala.
[36] G. E. Uhlenbeck and L. S. Ornstein. On the Theory of the Brownian Motion. Phys.
Rev., 36:823–841, 1930.
[37] H. A. Kramers. Brownian motion in a field of force and the diffusion model of
chemical reactions. Physica, 7(4):284–304, 1940.
118 BIBLIOGRAPHY
[38] S. Chandrasekhar. Stochastic Problems in Physics and Astronomy. Rev. Mod. Phys.,
15(1):1–89, 1943.
[39] M. C. Wang and G. E. Uhlenbeck. On the Theory of the Brownian Motion II. Rev.
Mod. Phys., 17:323–342, 1945.
[45] A. Kolmogoroff. Zur Theorie der stetigen zufälligen Prozesse. Math. Annalen,
108(1):149–160, 1933.
[47] W. Feller. Zur Theorie der stochastischen Prozesse (Existenz- und Ein-
deutigkeitssätze). Ann. Math., 113(1):113–160, 1936.
[48] P. Lévy. Sur certains processus stochastique homogènes. Compositio Math., 7:283–
339, 1939.
[50] K. Ito. Stochastic integral. Proc. Imp. Acad. Tokyo, 20:519–524, 1944.
[51] K. Ito. On stochastic differential equations. Mem. Amer. Mathem. Soc., 4:51–89,
1951.
[53] I. I. Gihman. Certain differential equations with random functions. Ukrain. Math.
Zh., 2(3):45–69, 1950. In Russian.
BIBLIOGRAPHY 119
[55] D. Fisk. Quasi-martingales and stochastic integrals. PhD thesis, Michigan State
University, Dept. of Statistics, 1963.
[59] R. L. Stratonovich. Conditional Markov processes and their application to the theory
of optimal control. American Elsevier Pub. Co., 1968.
[60] H. P. MacKean Jr. Stochastic Integrals. Academic Press, New York, 1969.
[62] D. R. Cox and H. D. Miller. The Theory of Stochastic Processes. Chapman and Hall,
Oxford, 1972.
[63] H. Risken and T. Frank. The Fokker-Planck Equation: Methods of Solutions and
Applications. Springer Series in Synergetics. Springer, Berlin, 1996.
[64] I. Karatzas and S. E. Shreve. Brownian Motion and Stochastic Calculus. Number
113 in Graduate Texts in Mathematics. Springer, New York, Berlin, 2 edition, 1991.
[67] P. Hänggi and H. Thomas. Stochastic processes: Time evolution, symmetries and
linear response. Phys. Rep., 88(4):207–319, 1982.
[68] P. Hänggi, P. Talkner, and M. Borkovec. Reaction rate theory: fifty years after
Kramers. Rev. Mod. Phys., 62(2):251–341, 1990.
[69] P. Reimann. Brownian motors: noisy transport far from equilibrium. Phys. Rep.,
361(2):57–265, 2002.
120 BIBLIOGRAPHY
[71] P. Hänggi and F. Marchesoni. 100 Years of Brownian Motion. Chaos, 15:026101,
2005.
[72] E. Frey and K. Kroy. Brownian motion: a paradigm of soft matter and biological
physics. Ann. Phys. (Leipzig), 14(1-3):20–50, 2005.
[73] I. M. Sokolov and J. Klafter. From diffusion to anomalous diffusion: A century after
Einstein’s Brownian motion. Chaos, 15:026103, 2005.
[78] J.-P. Bouchaud and M. Potters. Theory of Financial Risks: From Statistical Physics
to Risk Management. Cambridge University Press, Cambridge, 2001.
[81] V. B. Magalinskii. Dynamical model in the theory of the Brownian motion. Sov.
Phys. JETP, 9:1381–1382, 1959.
[84] E. Pollak. Theory of activated rate processes: A new derivation of Kramers’ expres-
sion. J. Chem. Phys., 85:865–867, 1986.
[85] A. M. Levine, M. Shapiro, and E. Pollak. Hamiltonian theory for vibrational dephas-
ing rates of small molecules in liquids. J. Chem. Phys., 88(3):1959–1966, 1988.
[86] P. Hänggi. Generalized Langevin equations: A useful tool for the perplexed modeller
of non-equilibrium fluctuations. In Stochastic dynamics, volume 484 of Lecture Notes
in Physics, pages 15–22. Springer, Heidelberg, 1997.
[90] P. Talkner. The failure of the quantum regression hypothesis. Ann. Phys. (N.Y.),
167(2):390–436, 1986.
[91] H. Grabert, U. Weiss, and P. Talkner. Quantum theory of the damped harmonic
oscillator. Z. Phys. B Cond. Mat., 55:87–94, 1984.
[93] G. W. Ford and M. Kac. On the Quantum Langevin Equation. J. Stat. Phys.,
46(5/6):803–810, 1987.
[96] S. Tsonchev and P. Pechukas. Binary collision model for quantum Brownian motion.
Phys. Rev. E, 61:6171–6182, 2000.
[97] S. Kohler, J. Lehmann, and P. Hänggi. Driven quantum transport on the nanoscale.
Phys. Rep., 406:379–443, 2005.
122 BIBLIOGRAPHY
[99] P. Hänggi and G. L. Ingold. Quantum Brownian motion and the third law of ther-
modynamics. Acta Phys. Polonica B, 37(5):1537–1549, 2006.
[100] R. F. Fox. Gaussian stochastic processes in physics. Phys. Rep., 48(3):179–283, 1978.
[102] Y. L. Klimontovich. Statistical Theory of Open Systems. Kluwer Acad. Publ., Dor-
drecht, 1995.
[103] P. Hänggi and P. Jung. Colored Noise in Dynamical Systems. Adv. Chem. Phys.,
89:239–326, 1995.
[105] J.-P. Bouchaud and A. Georges. Anomalous diffusion in disordered media: Statistical
mechanisms, models and physical applications. Phys. Rep., 195(4-5):127–293, 1990.
[106] R. Metzler and J. Klafter. The random walk’s guide to anomalous diffusion: a
fractional dynamics approach. Phys. Rep., 339(1):1–77, 2000.
[107] H. Grabert, P. Schramm, and G.-L. Ingold. Quantum brownian motion: The func-
tional integral approach. Phys. Rep., 168(3):115–207, 1988.
[108] M. Grifoni and P. Hänggi. Driven quantum tunneling. Phys. Rep., 304:229–354,
1998.
[109] A. Einstein. Über das Relativitätsprinzip und die aus demselben gezogenen Folgerun-
gen. Jahrbuch der Radioaktivität und Elektronik, 4:411–462, 1907.
[110] M. Planck. Zur Dynamik bewegter Systeme. Sitzungsber. Preuss. Akad. Wiss. Berlin,
pages 542–570, 1907.
[111] M. Planck. Zur Dynamik bewegter Systeme. Ann. Phys. (Leipzig), 26:1–34, 1908.
[115] H. Ott. Lorentz-Transformation der Wärme und der Temperatur. Z. Phys., 175:70–
104, 1963.
[116] A. Børs. Note on the temperature of a gas moving at relativistic speed. Proc. Roy.
Soc. (London), 86:1141–1143, 1965.
[120] H. Arzelies. Comment on dr. kibble’s article. Nuovo Cimento B, 41(1):81–82, 1966.
[121] F. Rohrlich. True and Apparent Transformations, Classical Electrons and Relativistic
Thermodynamics. Nuovo Cimento B, 45:76–83, 1966.
[126] P. T. Landsberg. Does a moving body appear cool? Nature, 212:571–572, 1966.
[127] J. H. Fremlin. Does a moving body appear cool? Nature, 213:277, 1967.
[128] P. T. Landsberg. Does a moving body appear cool? Nature, 214:903–904, 1967.
[129] P. D. Noerdlinger. A moving body must “appear” cool. Nature, 213:1117, 1967.
124 BIBLIOGRAPHY
[131] I. P. Williams. Does a moving body appear cooler. Nature, 214:1105, 1967.
[135] A. G. Georgiu. Special Relativity and Thermodynamics. Proc. Cambridge Phil. Soc.,
66:423–429, 1966.
[139] N. G. van Kampen. Relativistic thermodynamics. J. Phys. Soc. Jap. Suppl., 26:316–
321, 1969.
[140] N. G. van Kampen. A model for relativistic heat transport. Physica, 46:315–332,
1970.
[141] I. Müller. Toward relativistic thermodynamics. Arch. Rat. Mech. Anal., 34:259–282,
1969.
[146] P. T. Landsberg and K. A. Johns. The Lorentz Transformation of Heat and Work.
Ann. Phys., 56:299–318, 1970.
BIBLIOGRAPHY 125
[156] I-Shih Liu, I. Müller, and T. Ruggeri. Relativistic thermodynamics of gases. Ann.
Physics (N.Y.), 169:191–219, 1986.
[158] P. T. Landsberg and G. E. A. Matsas. Laying the ghost of the relativistic temperature
transformation. Phys. Lett. A, 223:401–403, 1996.
[162] D. Ter Haar and H. Wegeland. Thermodynamics and statistical mechanics on the
special theory of relativity. Phys. Rep., 1(2):31–54, 1971.
126 BIBLIOGRAPHY
[166] K. von Mosengeil. Theorie der stationären Strahlung in einem gleichförmig bewegten
Hohlraum. Ann. Phys. (Leipzig), 22:867–904, 1907.
[168] J. C. Maxwell. On the dynamical theory of gases. Phil. Trans. Roy. Soc. London,
157:49–88, 1867.
[170] F. Jüttner. Die relativistische Quantentheorie des idealen Gases. Z. Phys., 47:542–
566, 1928.
[173] John E. Krizan. Matter and radiation in equilibrium. Phys. Rev. D, 31(8):1968–1971,
1985.
[174] W. Lucha and F. F. Schöberl. Relativistic virial theorem. Phys. Rev. Lett.,
64(23):2733–2735, 1990.
[179] N. A. Chernikov. Equilibrium distribution of the relativistic gas. Acta Phys. Pol.,
26:1069–1092, 1964.
[180] C. Møller. Gibbs’ Statistical Mechanics in the Theory of Relativity. Kgl. Dan.
Vidensk. Selsk., Mat.-Fys. Medd., 36(16):1–44, 1968.
[185] R. Balescu, T. Kotera, and E. Piña. Lorentz transformations in phase space and in
physical space. Physica, 33(3):581–594, 1967.
[186] R. Hakim. Relativistic perfect gas in an external force field. Phys. Rev., 162:128–133,
1967.
[193] R. Hakim. Remarks on Relativistic Statistical Mechanics. II. Hierarchies for the
Reduced Densities. J. Math. Phys., 8(7):1397–1400, 1965.
[194] R. Hakim. Sur la mécanique statistique relativiste. Ann. Inst. Henri Poincaré (A),
6(3):225–244, 1967.
128 BIBLIOGRAPHY
[203] G. Kaniadakis. Statistical mechanics in the context of special relativity. Phys. Rev.
E, 66:056125, 2002.
[207] R. Silva and J. A. S. Lima. Relativity, nonextensivity and extended power law
distributions. Phys. Rev. E, 72:057101, 2005.
[208] G. Kaniadakis. Statistical mechanics in the context of special relativity. ii. Phys.
Rev. E, 72:036108, 2005.
[210] R. Silva. The relativistic statistical theory and Kaniadakis entropy: an approach
through a molecular chaos hypothesis. Eur. Phys. J. B, 54(4):499–502, 2006.
[212] F. Becattini and L. Ferroni. The microcanonical ensemble of the ideal relativistic
quantum gas. arXiv:nucl-th/0704.1967v1.
[213] G. Lina, D. Jiulin, and L. Zhipeng. The property of kappa-deformed statistics for a
relativistic gas in an electromagnetic field: kappa parameter and kappa-distribution.
arXiv:0704.1437v1.
[218] B. Cleuren, K. Willaert, A. Engel, and C. Van den Broeck. Fluctuation theorem
for entropy production during effusion of a relativistic ideal gas. Phys. Rev. E,
77(2):022103, 2008.
[222] J. A. Wheeler and R. P. Feynman. Interaction with the Absorber as the Mechanism
of Radiation. Rev. Mod. Phys., 17(2-3):157–181, 1945.
[224] H. Van Dam and E. P. Wigner. Classical Relativistic Mechanics of Interacting Point
Particles. Phys. Rev., 138(6B):B1576–B1582, 1965.
[225] H. Van Dam and E. P. Wigner. Instantaneous and Asymptotic Conservation Laws
for Classical Relativistic Mechanics of Interacting Point Particles. Phys. Rev.,
142(4):838–843, 1966.
[229] A. D. Fokker. Ein invarianter Variationssatz für die Bewegung mehrerer elektrischer
Massenteilchen. Z. Physik A: Hadrons and Nuclei, 58(5-6):386–393, 1929.
[230] M. H. L. Pryce. The mass-centre in the restricted theory of relativity and its con-
nexion with the quantum theory of elementary particles. Proc. Roy. Soc. London,
195(1040):62–81, 1948.
[239] T. Gill and J. Lindesay. Canonical proper time formulation of relativistic particle
dynamics. Int. J. Theor. Phys., 32(11):2087–2098, 1993.
[240] P. Ramond. Action-at-a - Distance Theories and Dual Models. Phys. Rev. D,
7(2):449–458, 1973.
[255] F. Rohrlich. Many-body forces and the cluster decomposition. Phys. Rev. D,
23(6):1305–1315, 1981.
[258] L. P. Horwitz and F. Rohrlich. Constraint relativistic quantum dynamics. Phys. Rev.
D, 24(6):1528–1542, 1981.
[261] G. Longhi and L. Lusanna. Bound-state solutions, invariant scalar products, and
conserved currents for a class of two-body relativistic systems. Phys. Rev. D,
34(12):3707–3731, 1986.
[262] M. Montesinos, C. Rovelli, and T. Thiemann. SL(2, R) model with two Hamiltonian
constraints. Phys. Rev. D, 60(4):044009, 1999.
[263] Cheuk-Yin Wong and Horace W. Crater. Relativistic N-body problem in a separable
two-body basis. Phys. Rev. C, 63(4):044907, 2001.
BIBLIOGRAPHY 133
[264] V. Gerdt, A. Khvedelidze, and Yu. Palii. Towards an algorithmisation of the dirac
constraint formalism. In J. Calmet, W. M. Seiler, and R. W. Tucker, editors, Global
Integrability of Field Theories. Proceedings of GIFT 2006., pages 135–154, Daresbury,
UK, 2006. Cocroft Institute.
[265] R. Hakim. Statistical mechanics of relativistic dense matter systems. Riv. Nuovo
Cim., 1(6):1–52, 1978.
[266] R. Hakim and H. D. Sivak. Relativistic statistical mechanics, a brief overview. AIP
Conf. Proc., 841:63–99, 2006. A Century of Relativity Physics: ERE 2005; XXVIII
Spanish Relativity Meeting.
[268] R. L. Liboff. Kinetic Theory. Prentice Hall, Inc., Englewood Cliffs, New Jersey
07632, 1990.
[272] S. T. Beliaev and G. I. Budker. Relativistic kinetic equation. Dokl. Akad. Nauk
SSSR, 6:807–810, 1956. English transl.: Soviet Phys. Doklady 1, 218 (1956).
[273] W. Israel. Relativistic kinetic theory of a simple gas. J. Math. Phys., 4(9):1163–1181,
1963.
[278] H. Akama. Relativistic Boltzmann equation for plasmas. J. Phys. Soc. Japan,
28(2):478–488, 1970.
[280] R. J. Gould. Kinetic theory of relativistic plasmas. Phys. Fluids, 24(1), 1981.
[287] L. Hsiao and H. Yu. Asymptotic stability of the relativistic Maxwellian. Math. Meth.
Appl. Sci., 29(13):1481 – 1499, 2006.
[291] J. Ehlers. General relativity and kinetic theory. In R. K. Sachs, editor, General
Relativity and Cosmology, pages 1–70. Academic Press, New York, 1971. Proc. of
the Int. School of Physics ”Enrico Fermi”, Course XLVII.
[293] G. Schay. The equations of diffusion in the special theory of relativity. PhD thesis,
Princeton University, 1961. available through University Microfilms, Ann Arbor,
Michigan, https://wwwlib.umi.com.
[296] H. van Erkelens. Relativistic Boltzmann theory for a plasma. VI. The relativistic
Landau equation. Physica, 107A:48–70, 1981.
[298] C. F. F. Karney and N. J. Fisch. Efficiency of current drive by fast waves. Phys.
Fluids, 28(1):116–126, 1985.
[300] Ph. de Gottal and J. Gariel. Test-particle motion in a relativistic plasma. II. The
Lorentz gas. Physica A, 163:615–624, 1990.
[301] Ph. de Gottal and J. Gariel. Test-particle motion in a relativistic plasma. III. Col-
lective effects. Physica A, 186:534–548, 1992.
[308] B. Svetitsky. Diffusion of charmed quarks in the quark-gluon plasma. Phys. Rev. D,
37(9):2484–2491, 1988.
[309] P. Roy, J. E. Alam, S. Sarkar, B. Sinha, and S.Raha. Quark-gluon plasma diagnostics
in a successive equilibrium scenario. Nucl. Phys. A, 624(4):687–705, 1997.
[310] H. van Hees and R. Rapp. Thermalization of heavy quarks in the quark-gluon plasma.
Phys. Rev. C, 71(3):034907, 2005.
[312] H. van Hees, V. Greco, and R. Rapp. Heavy-quark probes of the quark-gluon plasma
and interpretation of recent data taken at the BNL Relativistic Heavy Ion Collider.
Phys. Rev. C, 73:034913, 2006.
[313] H. van Hees, V. Greco, and R. Rapp. Thermalization and flow of heavy quarks in
the quark-gluon plasma. AIP Conf. Proc., 842:77–79, 2006.
[314] R. Rapp, V. Greco, and H. van Hees. Heavy-Quark Spectra at RHIC and Resonances
in the QGP. Nucl. Phys. A, 774:685–688, 2006.
[315] P. J. D. Mauas and D. O. Gomez. Kinetic description of electron beams in the solar
chromosphere. Astophys. J., 398:682–691, 1997.
[318] B. Wolfe and F. Melia. Covariant kinetic theory with an application to the coma
cluster. Astrophys. J., 638, 2006.
[320] E. Bertschinger. Brownian Motion of Stars, Dust, and Invisible Matter. AIP Conf.
Proc., 861(1):97–105, 2006.
[321] C. F. F. Karney. Fokker-Planck and quasilinear codes. Comp. Phys. Rep., 4(3-4):183–
244, 1986.
[322] M. Shoucri, V. Fuchs, and A. Bers. The application of the finite-element library
TWODEPEP for the numerical solution of the relativistic Fokker-Planck equation.
Comp. Phys. Comm., 46(3):337–344, 1987.
[324] T. H. Boyer. Inconsistency between the Boltzmann distribution for relativistic free
particles and the Rayleigh-Jeans law for thermal radiation within classical electro-
dynamics. Phys. Rev. D, 19:1112–1115, 1979.
[325] T. H. Boyer. Inconsistency between the Boltzmann distribution for relativistic free
particles and the Planck spectrum for thermal radiation in quantum theory. Phys.
Rev. D, 19:3635–3640, 1979.
[326] U. Ben-Ya’acov. Relativistic Brownian motion and the spectrum of thermal radiation.
Phys. Rev. D, 23(6):1441–1450, 1981.
[327] R. Blanco, L. Pesquera, and E. Santos. Equilibrium between radiation and mat-
ter for classical relativistic multiperiodic systems. Derivation of Maxwell-Boltzmann
distribution from Rayleigh-Jeans spectrum. Phys. Rev. D, 27(6):1254–1287, 1983.
[329] O. Oron and L. P. Horwitz. Relativistic Brownian Motion and Gravity as an Eikonal
Approximation to a Quantum Evolution Equation. Found. Phys., 35(7):1181 – 1203,
2005.
[330] F. Castejón and S. Eguilior. Langevin equations for quasi-linear wave–particle inter-
action. Plasma Phys. Control. Fusion, 45(2):159–167, 2003.
[331] K. S. Fa. Analysis of the relativistic brownian motion in momentum space. Braz. J.
Phys., 36(3A):777–782, 2006.
[333] R. Zygadlo. Two Markovian models leading to the Pascalian (1D ultrarelativistic)
stationary distribution. Phys. Lett. A, 329:459–463, 2004.
[334] R. Zygadlo. Free Brownian motion approach to the ideal gas of relativistic particles.
Phys. Lett. A, 345:323–329, 2005.
[337] R. Rapp and H. van Hees. Heavy quark diffusion as a probe of the quark-gluon
plasma. arXiv:0803.0901v1, 2008.
[339] J. Lopuszaǹski. Relativisierung der theorie der stochastischen prozesse. Acta Phys.
Polon., 12:87–99, 1953.
[340] H. Rudberg. On the Theory of Relativistic Diffusion. Almquist and Wiksells, Upp-
sala, 1957.
[343] R. M. Dudley. Asymptotics of Some Relativistic Markov Processes. Proc. Nat. Acad.
Sci. USA, 70(12):3551–3555, 1973.
[347] E. Nelson. Derivation of the Schrödinger Equation from Newtonian Mechanics. Phys.
Rev., 150(4):1079–1085, 1966.
BIBLIOGRAPHY 139
[351] J. P. Caubet. Relativistic Brownian motion. Adv. Appl. Prob., 8(2):230–231, 1976.
[354] J. P. Caubet. Brownian Motion, Minkowski Space and Principle of Special Relativity.
Comp. Rend. Sci. Acad. Serie A, 284(5):339–341, 1977.
[356] F. Guerra and P. Ruggiero. A note on relativistic Markov processes. Lett. Nouvo
Cimento, 23:529–534, 1978.
[357] F. Guerra and P. Ruggiero. New Interpretation of the Euclidean-Markov Field in the
Framework of Physical Minkowski Space-Time. Phys. Rev. Lett., 31(16):1022–1025,
1978.
[359] S. Roy. Relativistic Brownian motion and the space–time approach to quantum
mechanics. J. Math. Phys., 21:71–76, 1980.
[360] S. Roy. Relativistic Brownian motion and the theory of measurement in quantum
mechanics. Nuovo Cim. B, 64:81–93, 1981.
[361] J. P. Caubet. The Brownian Theory of Light. Stoch. Anal. Appl., 2(1):55–85, 1984.
[362] J. P. Caubet. The great fugue of the Brownian theory of light. Stoch. Anal. Appl.,
3(2):119–151, 1985.
140 BIBLIOGRAPHY
[365] L. M. Morato and L. Viola. Markov diffusions in comoving coordinates and stochastic
quantization of the free relativistic particle. J. Math. Phys., 36:4691–4710, 1995.
Erratum, ibid. 37, 4769 (1996).
[367] M. Nagasawa. Quantum theory, theory of Brownian motions, and relativity theory.
Chaos, Solitons & Fract., 7(5):631–643, 1996.
[368] A. Posilicano. Poincareé invariant Markov processes and Gaussian random fields on
relativistic phase-spaces. Lett. Math. Phys., 42:85–93, 1997.
[369] M. Pavon. On the stochastic mechanics of the free relativistic particle. J. Math.
Phys., 42(10):4946–4856, 2001.
[374] J. Franchi and Y. Le Jan. Relativistic diffusions and Schwarzschild geometry. Comm.
Pure Appl. Math., 60:187–251, 2006.
[376] M. Kac. A stochastic model related to the telegrapher’s equation. Rocky Mountain
J. Math., 4:497–509, 1974.
[379] R. Fürth. Die Brownsche Bewegung bei Berücksichtigung einer Persistenz der Bewe-
gungsrichtung. Mit Anwendungen auf die Bewegung lebender Infusorien. Z. Physik,
2:244–256, 1920.
[380] G. I. Taylor. Diffusion by continuous movements. Proc. London Math. Soc., 20:196–
212, 1922.
[382] M. Boguñá, J. M. Porrà, and J. Masoliver. Persistent random walk model for trans-
port through thin slabs. Phys. Rev. E, 59(6):6517–6526, 1999.
[383] A. Ishimaru. Diffusion of light in turbid material. Appl. Opt., 28:2210–2215, 1989.
[384] M. Miri and H. Stark. Persistent random walk in a honeycomb structure: Light
transport in foams. Phys. Rev. E, 68(3):031102, 2003.
[385] M. Miri, E. Madadi, and H. Stark. Diffusive transport of light in a Kelvin foam.
Phys. Rev. E, 72(3):031111, 2005.
[388] D. D. Joseph and L. Preziosi. Heat waves. Rev. Mod. Phys., 61(1):41–73, 1989.
[389] D. D. Joseph and L. Preziosi. Addendum to the paper ”Heat waves” [Rev. Mod.
Phys. 61, 41 (1989)]. Rev. Mod. Phys., 62(2):375–391, 1990.
[390] P. Kostädt and M. Liu. Causality and stability of the relativistic diffusion equation.
Phys. Rev. D, 62:023003, 2000.
142 BIBLIOGRAPHY
[391] P. Kostädt and M. Liu. Alleged acausality of the diffusion equations: A reply. Phys.
Rev. D, 64:088504, 2001.
[392] L. Herrera and D. Pavón. Why hyperbolic theories of dissipation cannot be ignored:
Comment on a paper by Kostädt and Liu. Phys. Rev. D, 64:088503, 2001.
[393] M. Abdel-Aziz and S. Gavin. Causal diffusion and the survival of charge fluctuations
in nuclear collisions. Phys. Rev. C, 70(3):034905, 2004.
[394] M. Abdel-Aziz and S. Gavin. Causal diffusion at RHIC. J. Phys. G: Nucl. Phys.,
31:S77–S84, 2005.
[395] T. Koide. Microscopic derivation of causal diffusion equation using the projection
operator method. Phys. Rev. E, 72:026135, 2005.
[399] T. Ichinose. Path integral for the Dirac equation in two space-time dimensions. Proc.
Japan Acad., 58:290–293, 1982.
[400] T. Ichinose. Path integral formulation of the propagator for a two-dimensional Dirac
particle. Physica A, 124(1-3):419–425, 1984.
[401] T. Jacobson and L. S. Schulman. Quantum stochastics: the passage from a relativistic
to a non-relativistic path integral. J. Phys. A: Math. Gen., 17:375–383, 1984.
[402] A. O. Barut and I. H. Duru. Path-Integral Derivation of the Dirac Propagator. Phys.
Rev. Lett., 53(25):2355–2358, 1984.
[403] B. Gaveau and L. S. Schulman. Grassmann-valued processes for the Weyl and the
Dirac equations. Phys. Rev. D, 36(4):1135–1140, 1987.
[404] R. P. Feynman and A. R. Hibbs. Quantum Mechanics and Path Integrals. McGraw-
Hill, New York, 1965. problem 2-6, pp. 34-36.
[406] G. Wolschin. Diffusion and local deconfinement in relativistic systems. Phys. Rev.
C, 69(2):024906, February 2004.
[410] C. Nicolis and G. Nicolis. Stochastic aspects of climatic transitions - additive fluc-
tuations. Tellus, 33:225–234, 1981.
[411] C. Nicolis. Long-term climatic transitions and stochastic resonance. J. Stat. Phys.,
70(1/2):3–13, 1993.
[417] W. Ebeling. Nonlinear Brownian motions – mean square displacement. Cond. Mat.
Phys., 7(3(39)):539–550, 2004.
[418] D. Ryter and U. Deker. Properties of the noise-induced ”spurious” drift. I. J. Math.
Phys., 21:2662–2665, 1980.
[419] D. Ryter and U. Deker. Properties of the noise-induced (”spurious”) drift. II. Sim-
plifications of Langevin Equations. J. Math. Phys., 21:2666–2669, 1980.
[422] M. R. Hoare. The linear gas. Adv. Chem. Phys., 20:135–214, 1971.
[423] M. R. Hoare, S. Raval, and M. Rahman. The Rayleigh Model: Singular Transport
Theory in One Dimension. Philos. Trans. R. Soc. Lond. A, 305:383–440, 1982.
[426] J. W. Rayleigh. The Theory of Sound, volume 1. Dover, New York, 2 edition, 1945.
[428] S. Marksteiner, K. Ellinger, and P. Zoller. Anomalous diffusion and Lévy walks in
optical lattices. Phys. Rev. A, 53:3409–3430, 1996.
[429] E. Lutz. Anomalous diffusion and Tsallis statistics in an optical lattice. Phys. Rev.
A, 67:051402, 2003.
[435] M. Schienbein and H. Gruler. Langevin Equation, Fokker-Planck Equation and Cell
Migration. Bull. Math. Biology, 55:585–608, 1993.
BIBLIOGRAPHY 145
[436] W. Ebeling, F. Schweitzer, and B. Tilch. Active Brownian particles with energy
depots modelling animal mobility. BioSystems, 49:17–29, 1999.
[442] W. Ebeling, P. S. Landa, and V. Ushakov. Self-oscillations in Toda ring chains with
negative friction. Phys. Rev. E, 63:046601, 2001.
[445] S. Jespersen, R. Metzler, and H. C. Fogedby. Lévy flights in external force fields:
Langevin and fractional Fokker-Planck equations and their solutions. Phys. Rev. E,
59(3):2736–2745, 1999.
[447] A. Weron, K. Burnecki, S. Mercik, and K. Weron. Complete description of all self-
similar models driven by Lévy stable noise. Phys. Rev. E, 71(1):016113, 2005.
[453] P. A. M. Dirac. The conditions for statistical equilibrium between atoms, electrons
and radiation. Proc. Roy. Soc. London A, 106:581–596, 1924.
[456] M. Jones et al. An image of the Sunyaev-Zel’dovich effect. Nature, 365:320 – 323,
1993.
[460] A. Wehrl. General properties of entropy. Rev. Mod. Phys., 50(2):221–260, 1978.
[461] S. Kullback and R. A. Leibler. On information and sufficiency. Ann. Math. Statistics,
22:79–86, 1951.
[463] A. Haar. Der Massbegriff in der Theorie der kontinuierlichen Gruppen. Ann. Math.,
34(1):147–169, 1933.
[465] J. Masoliver and J. Marro. Kinetics of a finite one-dimensional mixture of hard rods
with different masses. J. Stat. Physics, 31(3):565–575, 1983.
[467] Yu. M. Sinyukov. On the relativistic partition function of ideal gases. Phys. Lett. B,
127(6):443–447, 1983.
[470] A. Fingerle. Relativistic fluctuation theorems: Theory and explicit examples. Compt.
Rend. Acad. Sci. Physique, 8(5-6):696–713, 2007. arXiv:0705.3115v1.
[471] L. Zaninetti. X-ray profiles in symmetric and asymmetric supernova remnants. Baltic
Astron., 16:251–285, 2007.
[472] M. Rigotti. A relativistic stochastic process. Master’s thesis, ETH Zürich, Institute
for Theoretical Physics, 2005.
[483] W. Thomson. On the Theory of the Electric Telegraph. Proc. Roy. Soc. London,
Series I, 7:382–399, 1854-1855.
[485] P. A. M. Dirac. Classical theory of radiating electrons. Proc. Roy. Soc. London Ser.
A, 167(929):148–169, 1938.
[486] C. J. Eliezer. On the classical theory of particles. Proc. Royal Soc. London, 194:543–
555, 1948.
[488] F. Rohrlich. The self-force and radiation reaction. Amer. J. Phys., 68(12):1109–1112,
2000.
[489] H. Spohn. The Critical Manifold of the Lorentz-Dirac Equation. Europhys. Lett.,
50:287–292, 2000.
[490] G. Ares de Parga, R. MAres, and S. Dominguez. An unphysical result for the Landau-
Lifshitz equation of motion for a charged particle. Rev. Mexicana Fis., 52(2):139–142,
2006.
[495] E. Fermi. On the origin of the cosmic radiation. Phys. Rev., 75(8):1169–1174, 1949.
[496] C. Jarzynski. Nonequilibrium Equality for Free Energy Differences. Phys. Rev. Lett.,
78(14):2690–2693, 1997.
[497] U. Seifert. Entropy Production along a Stochastic Trajectory and an Integral Fluc-
tuation Theorem. Phys. Rev. Lett., 95:040602, 2005.
[500] P. Talkner, E. Lutz, and P. Hänggi. Fluctuation theorems: Work is not an observable.
Phys. Rev. E, 75(5):050102, 2007.
[502] H.-P. Breuer and F. Petruccine. A Lorentz covariant stochastic wave function dy-
namics for open systems. Phys. Lett. A, 242(4-5):205–210, 1998.
[503] H.-P. Breuer and F. Petruccine. Relativistic formulation of quantum state diffusion.
J. Phys. A: Math. Gen., 31(1):33–52, 1998.
[504] H.-P. Breuer and F. Petruccine. Reply to ‘Relativistic formulation of quantum state
diffusion?’. J. Phys. A: Math. Gen., 31(47):9605–9612, 1998.
150 BIBLIOGRAPHY
[509] M. Rigotti and F. Debbasch. An H-theorem for the general relativistic Ornstein-
Uhlenbeck process. J. Math. Phys., 46:103303, 2005.
Lebenslauf
151
152 BIBLIOGRAPHY
Danksagung
Ich möchte mich bei allen Kollegen und Freunden bedanken, die im Laufe der vergangen
Jahre durch Vorlesungen, Bemerkungen und Kommentare direkt oder indirekt zur der
vorliegenden Dissertation beigetragen haben. Besonders großer Dank gilt dabei meinem
Betreuer Peter Hänggi, ohne dessen Anregungen, Unterstützung und Enthusiasmus diese
Arbeit nicht möglich gewesen wäre. Ein weiteres herzliches Dankeschön richtet sich an
Peter Talkner und Gert Ingold für ihre bewundernswerte Geduld beim Beantworten meiner
Fragen. Also, many thanks to David and Jesus (und Stefan) for the fruitful collaboration.
Zurückblickend bin ich sehr glücklich darüber, daß ich während der vergangenen drei
Jahre in einer so angenehmen, freundlichen und hilfsbereiten Umgebung arbeiten und ler-
nen durfte. Dafür möchte ich allen ehemaligen und gegenwärtigen Mitarbeitern, Postdocs,
Doktoranden und Zimmerkollegen aus theo1 nocheinmal ganz herzlich danken. Ihr alle
habt dafür gesorgt, daß ich das tägliche Pendeln zwischen Augsburg und München gerne
in Kauf genommen habe!
153