Stochastic thermodynamics of relativistic Brownian motion

Physical scenarios that require a relativistic treatment are ubiquitous in nature, ranging from cosmological objects to charge carriers in Dirac materials. Interestingly all of these situations have in common that the systems typically operate very far from thermal equilibrium. Therefore, if and how the framework of stochastic thermodynamics applies at relativistic energies is a salient question. In the present work we generalize the notions of stochastic heat and work for the relativistic Langevin equation and derive the fluctuation theorems without and with feedback. For processes with feedback we consider the ramifications of the lack of simultaneity of events in the inertial frames of observer and Brownian particles, and we argue that the framework of absolute irreversibility is instrumental to avoid acausal considerations. The analysis is concluded with a few remarks on potential experimental applications in graphene.


Introduction
The Universe is full of objects that defy the common anthropocentric understanding of nature. A good example are so-called blazars, which are caused by narrow jets of cosmic particles moving at relativistic speeds directly towards Earth [1]. To date, only very little is known about these jets since "most of the action" occurs over length scales that are too small to be resolved even with the best telescopes. Curiously, what is known [1], is that these jets are powered by supermassive black hole "engines" that propel charged particles into the Universe.
From a thermodynamic point of view, these cosmic jets are, thus, nothing else but ensembles of Brownian particles that are driven very far from thermal equilibrium. Therefore, it seems plausible that recent advances in stochastic thermodynamics [2][3][4][5][6][7] should be useful in describing the thermal properties also of cosmic jets. However, with the exception of the fluctuation theorem derived by Fingerle [8], none of the existing results appear directly applicable at relativistic energies.
The purpose of the present paper is to generalize the framework of stochastic thermodynamics to relativistic Brownian motion. For the sake of simplicity, we will focus on the relativistic Ornstein-Uhlenbeck process [40][41][42][43], for which the Langevin equation has been extensively studied [44][45][46][47]. As main results, we will identify the relativistic expression of the stochastic heat, work, and entropy production. For the entropy production we will further prove the detailed as well as the integral fluctuation theorems. In a second part of the analysis, we will study relativistic Brownian motion under feedback. For such a relativistic version of Maxwell's demon, we will derive a statement of the second law that is reminiscent of formulations for absolutely irreversible processes [48]. To this end, we will elucidate how monitored Brownian motion becomes absolutely irreversible due to the relativistic time delay. The analysis will be concluded with a few remarks on potential experimental implementations in novel materials with effectively relativistic dispersion relations.

Stochastic energetics of the relativistic Ornstein-Uhlenbeck process
Consider a Brownian particle in one dimension with rest mass M and in contact with a heat reservoir at temperature T . As usual, the heat bath is taken to be at rest with respect to the inertial frame of the laboratory [40,41]. Then, the dynamics of the Brownian particle is described by the relativistic Langevin equatioṅ where x is the position, p denotes the canonical momentum, and p 0 = p 2 + M 2 is the relativistic kinetic energy. Note that to avoid clutter in the formulas we work in units such the speed of light c = 1. Further, γ is the viscous friction coefficient, and D is the diffusion coefficient for the Gaussian white noise ξ(t), with ξ(t) = 0 and ξ(t)ξ(t ′ ) = δ(t − t ′ ). Finally, U (x, α(t)) is the time-dependent potential and α(t) is an external control parameter. It is then a simple exercise to generalize stochastic energetics [2] to the relativistic Langevin equation (1). To this end, consider the relativistic energy-momentum relation, from which we obtain Inspecting Eq. (3) we immediately recognize the thermodynamic work for a process of length τ [2, 49] which is the identical expression as found in non-relativistic Brownian motion. Accordingly, the heat becomes Note that both heat, Q, and work, W , explicitly depend on the fluctuating position and momentum, and thus also Q and W are stochastic quantities. The relativistic expressions for work and heat in Eqs. (4) and (5), respectively, constitute our first main result. It is interesting to note that the expression for the stochastic work is identical to the non-relativistic case. Work is commonly understood as an "ordered" exchange of internal energy [50]. Thus, from the point of view of the Brownian particle, how to compute work should not depend on how fast the particle is traveling. Heat, on the other hand, is the energy "lost" into the heat bath, which in the present scenario is at rest with respect to the Brownian particle. Microscopically, heat is exchanged through random interaction of the Brownian particle with the degrees of freedom of the bath. Thus, the relativistic description of momentum transfer is pertinent, and one would already intuitively expect that the expression for heat does depend on the relative motion.

Stochastic thermodynamics and fluctuation theorems
We now continue the analysis by focusing on the relativistic entropy production. To this end, we will first identify the internal and external contributions by comparing with the stochastic heat (5), before we will proof the fluctuation theorems.

Relativistic entropy production
The stochastic entropy production is more conveniently described through the timeevolution of the phase space density ρ(Γ, t) [3,[51][52][53]. Thus, we consider the relativistic Fokker Planck equation [44] where Γ = (x, p) is a point in phase space, ∇ ≡ (∂ x , ∂ p ), and j(Γ, t) is the probability current. Under time-reversal, the position, x, is an even function whereas the momentum, p, is an odd function. Thus, it has proven convenient to separate j into two contributions, j ≡ j r + j d [53]. We have which is symmetric under time-reversal, and a dissipative part that is asymmetric under time-reversal. It is this second contribution that is the origin of irreversible entropy production [52,53].

3.1.1.
Relativistic fluctuation-dissipation theorem and detailed balance Before analyzing the entropy production we need to consider the proper definition of thermodynamic temperature. To this end, consider that in equilibrium for a specific value of α the dissipative part of the current (8) vanishes, j d = 0, and therefore we also have ∇ · j r = 0. Accordingly, the stationary solution of (6) can be written as, where E(x, α) is again the relativistic energy (2). Now demanding that the equilibrium state (9) is given by the Maxwell-Jüttner distribution [44,54] we obtain the corresponding fluctuation-dissipation theorem which is identical to the non-relativistic case.
3.1.2. Second law of relativistic stochastic thermodynamics Now, consider the nonequilibrium Shannon entropy of the Brownian particle whose rate can be written aṡ Using standard arguments [53] it is easy to see that the contribution from the timereversal current, j r vanishes. Thus, we have after integrating by parts twice and simplifying expressionṡ and where we introduced the notation j p d = −γM p/p 0 ρ − D ∂ p ρ. To continue we now need to identify the contribution to the entropy production that is solely due to the exchange of heat between the Brownian particle and the thermal reservoir. Thus, inspecting the second term in Eq. (13) we define, Now re-writing the momentum component of the total probability current, j(Γ, t) = (j x , j p ), as we have where we also used thatṗ = j p /ρ. Finally, employing the relativistic Langevin equation (1) we can write, where we identified the stochastic heat (5) substituted the fluctuation-dissipation theorem (10). Therefore, the total irreversible entropy production becomeṡ which is always non-negative and of the standard form [53]. Equation (18) is a statement of the second law of thermodynamics for relativistic Brownian motion. It states that the irreversible entropy production is always nonnegative, and thus Eq. (18) can also be considered a relativistic generalization of the Clausius inequality.

Fluctuation theorem for the irreversible entropy production
Finally, we need to show that the fluctuation theorem also remains valid at relativistic energies. To this end, we will generalize the path integral formalism pioneered in Refs. [55,56]. As before we consider a relativistic Brownian particle that is driven by the time-dependent external control parameter α(t). Moreover, we denote a trajectory in phase space by X t ≡ {Γ(t)}. Now, consider a sequence of time steps of length δt with N δt = τ , such that the trajectory X t is separated into N infinitesimal segments. Then, ρ + i (Γ i+1 , t i + δt|Γ i , t i ) is the transition probability between two consecutive points along X t . It can be written as where ... ξi is an average over the noise with probability distribution and J + ξi,pi is the Jacobian. In the present case, the Jacobian becomes, which follows from straight forwardly evaluating terms.
Using the explicit expression for the noise distribution (20) the transition probability for the forward process reads The corresponding transition probability for the reverse process, i.e., the transition probability for the time-reversed driving becomes where J − ξi,pi is the reverse Jacobian. Note that the reverse Jacobian is identical to the Jacobian of the forward process, since it only contains even functions of the momentum [57].
We are now equipped with everything we need to compute the ratio of the transition probabilities for forward and reverse process. We have, where we used the definition of the integral in the first line, and identified the ∆s res in the second line. Further assuming that both, forward as well as reverse process start in the respective equilibrium state ρ s in and ρ s fin , respectively, the system entropy change is ∆s = k B ln(ρ s in /ρ s fin ). Therefore, using ρ F = ρ s in ρ + and ρ R = ρ s fin ρ − we finally obtain, where as before ∆s tot = ∆s + ∆s res . Equation (25) is the detailed fluctuation theorem for the total entropy production in relativistic Brownian motion. Accordingly, we also have an integral fluctuation theorem exp (−∆s tot /k B ) = 1. As a main result, we have shown that with the proper identification of the relativistic work (4) and heat (5) the major results of stochastic thermodynamics hold without modification also at relativistic energies.

Maxwell's demon at relativistic energies
Initiated by the seminal work of Sagawa and Ueda [26,27] the study of Maxwell's demon as experienced a renaissance. In the remainder of the analysis, we will now show that the fundamental postulates of special relativity inhibit the common working principles of the demon.
Einstein's central insight was that no object can be travel faster than the speed of light [58], and hence no information can be transferred faster than light. This observation is at variance with the usual set-up of Maxwell's gedankenexperiment, in which a neat fingered being observes the physical state of gas molecules and takes immediate action [59]. If the observed particles are moving at relativistic speed relative to the observer, any feedback will be subject to the relativistic time delay. This is nothing else but a direct consequence of lack of simultaneity of events in the inertial frames of the observer and the Brownian particle. Consequently, the natural question arises how this fact manifests itself in the statements of the second law.
To gain insight, we now consider a relativistic Brownian particle, whose dynamics is described bẏ where F fb [X t−τ0 ] is the force arising from autonomous or non-autonomous feedback [60], which depends on the microscopic trajectory upto time (t − τ 0 ). For instance, if F fb [X t−τ0 ] is proportional to the position x at time (t − τ 0 ), then τ 0 is just given by standard time dilation. Interestingly, the mathematical description of relativistic Brownian motion (26) is analogous to physical systems with response lags [61]. Accordingly, the energy balance, i.e., the first law of thermodynamics (3) now readsĖ =Q +Ẇ +Ẇ fb (27) whereẆ fb = (p/p 0 )F fb [X t−τ0 ] is the power input due to the feedback force. Equation (27) constitutes the generalization of relativistic stochastic energetics to systems with feedback. More interesting are the ramifications for the second law and the fluctuation theorems. To this end, we re-consider the transition probability between two increment (22). We immediately conclude that where as before the Jacobian is given by Eq. (21).
The problem now is how to identify the correct reverse or adjoint process, which should be compared in the fluctuation theorem [62]. Whereas already in the case of instantaneous feedback the question is non-trivial [49], the situation becomes significantly more involved at relativistic energies. Due to the relativistic time delay any direct time-reversal of the trajectories would necessarily lead to acausal events, and as such naive time-reversal cannot be considered thermodynamically sound.
A possible solution can be found by considering the recently developed notion of absolute irreversibility [48]. Consider a system that undergoes a thermodynamic process with microscopic trajectory X F that occurs with probability ρ F [X F ]. A time-reversed trajectory X R may then occur with probability ρ R [X R ]. The systems exhibits absolute irreversibility if either (i) ρ F [X F ] > 0 and ρ R [X R ] = 0, or (ii) ρ F [X F ] = 0 and ρ R [X R ] > 0. In other words, a process becomes absolutely irreversible if the probability to observe the reverse trajectory vanishes.
Realizing that the probability of acausal events occurring has to be zero the notion of absolute irreversibility appears well-suited to study relativistic Brownian motion. In complete analogy to the non-relativistic case [48] we therefore consider the "reverse" dynamics to be described bẏ which is nothing else but the time-reversed dynamics of Eq. (26) in the absence of any feedback. For this scenario the ratio of forward (26) and reverse (29) transition probabilities (24) becomes which differs from the case without feedback (24) only by one additional term. In complete analogy to before the probability distribution of all forward trajectories is then given by ρ F = ρ + ρ F in , where ρ in is again the initial distribution. Correspondingly we also have for the reverse process ρ R = ρ − ρ R in . Hence the detailed fluctuation theorem becomes where we introduced the dissipation function R, which reads We immediately recognize the first two terms as the stochastic entropy production ∆s tot , and thus we can also write which is the total dissipation function comprised of entropy production and irreversible work input due to the feedback force. Note that R in Eq. (33) contains only those pairs of trajectories, X F τ and X R τ , that are the time-reversal of each other. In particular, there exist some trajectories that are generated by Eq. (26), whose time-reversed version cannot be generated by the Langevin equation (29). Similarly, Eq. (29) generates trajectories that do not have a conjugate in the forward process.
Let Ξ F denote the set of all trajectories generated by the forward dynamics Eq. (26) that do have a time-reversed version and Ξ F c be the set of trajectories that do not have a time-reversed version, i.e., the complement of Ξ F . Similarly, Ξ R denotes the set of all trajectories generated by reverse dynamics Eq. (29) that have a conjugate in the forward process and Ξ R c is the complementary set consisting of all the other trajectories that are generated by reverse dynamics. Then, we can write where we used the detailed fluctuation theorem (31) in the second line. Now, denoting the total probability of reverse trajectories that do not have any counterpart in the forward process by Equation (35) constitutes a modified version of the integral fluctuation theorem due to the underlying absolutely irreversible processes, see also Ref. [48]. Further noting that λ and Λ are both positive numbers, we can write with the help of Jensen's inequality which is a generalized Clausius inequality for relativistic Brownian motion under feedback.

Concluding remarks
In the present analysis we have comprehensively shown how the framework of stochastic thermodynamics applies to Brownian motion at relativistic energies. Whereas cosmic gases and jets [1] may one day become the physical systems to be studied with this framework, there are experimental platforms that are somewhat "closer to home".

Potential experimental significance
In a rather recent experiment Pototsky et al. [63] demonstrated that relativistic Brownian motion can be studied in graphene. Graphene is essentially a single layer of graphite, i.e., carbon atoms in a honeycomb lattice. Its charge carriers have the remarkable property that they obey a linear dispersion relation, which means that the charge carriers behave like relativistic, massless Dirac fermions [64]. Pototsky et al. [63] then demonstrated that for strong electric fields and at high temperatures the quantum signatures are suppressed, and that the dynamics of the charge carriers can be effectively described by the relativistic Langevin equation (1). Our present results now open an avenue for the experimental study of relativistic stochastic thermodynamics in an already existing system.

Summary
In summary, we demonstrated that the hallmark results from stochastic thermodynamics, namely the fluctuation theorems with and without feedback hold at relativistic energies. To this end, we generalized stochastic energetics to the relativistic Langevin equation. We concluded that the expression for the stochastic work remains the same, whereas the stochastic heat accounts for the relativistic corrections for energy and momentum transfer between the traveling particle and a heat reservoir at rest. With this identification obtaining the fluctuation theorems without feedback turned out to be a straight forward exercise. On the other hand, relativistic systems experiencing feedback are more interesting, as due to the relativistic time dilation measurement and feedback action are not necessarily simultaneous for observer and relativistic particle. Thus, we argued that the framework of absolute irreversibility needs to be invoked to avoid formulating statements of the second law from acausal, i.e., unphysical events.