Beyond the relaxation time approximation

The relaxation time approximation (RTA) is a well known method of describing the time evolution of a statistical ensemble by linking distributions of the variables of interest at different stages of their temporal evolution. We show that if all the distributions occurring in the RTA have the same functional form of a quasi-power Tsallis distribution the time evolution of which depends on the time evolution of its control parameter, nonextensivity $q(t)$, then it is more convenient to consider only the time evolution of this control parameter.


Introduction
Many problems in science involve understanding the time evolution of a statistical ensemble. Here we will focus on a system of particles described by the probability distribution f (r, p,t) which depends on position r, momentum p and time t. In general, for an evolving physical system operating irreversibly out of thermodynamical equilibrium f (r, p,t) will differ from that of a Boltzmannian ensemble and its evolution is usually studied using the Boltzmann transport equation (BTE), the general form of which is [1], where F is the external force, u the velocity and C[ f ] the collision term. Narrowing our interest to situations where the considered system is homogeneous (i.e. ∇ r f = 0) and when no external forces are acting (i.e. when F = 0) Eq. (1) simplifies to the form a e-mail: grzegorz.wilk@ncbj.gov.pl b e-mail: zbigniew.wlodarczyk@ujk.edu.pl Because of the freedom in choosing the functional form of C[ f (t)] it is still a very general equation that allows one to deal with a variety of situations. However, usually in many applications it is further simplified to a form called the relaxation time approximation (RTA) which consists in using such a simple form of the collision term [1][2][3][4]: where f eq is the local equilibrium distribution and τ is the relaxation time, understood as the time taken by the nonequilibrium system to reach equilibrium. In this approximation the BTE simplifies further to Solving this equation for the initial conditions such that at t = 0 one has as initial (assumed) distribution, f = f in , and at freeze-out time, t = t f one has a final distribution, f = f f o (which we identify with the distribution we are looking for that actually describes the distribution obtained experimentally) 1 one finds that The continued popularity of such an approach to the analysis of various particle production processes can be proved by the fact that recently the Boltzmann transport equation in the RTA approximation was used to analyze various observables in nucleus-nucleus collisions (in particular, to study the time evolution of temperature fluctuations in a non-equilibrium 1 The statistical system produced in multiparticle production processes quickly reaches an initial distribution (pre-equilibrium state) which slowly evolves to equilibrium but becomes frozen at the freeze-out time (usually the system at freeze-out is not in thermodynamic equilibrium). The experimentally measured spectra of the produced particles reflect the state of the system at freeze-out.

Beyond the RTA -time evolution of internal variable
Note however that a given non-equilibrium process, regardless of the specific dynamics, evolves the probability distribution over the system configurations. The system traverses a probability manifold, which is not a manifold of equilibrium states because the distribution at each point on the manifold need not correspond to a BG distribution [11]. This means that for every t one has some distribution f (z,t) and all objects in the system are described by that (single) distribution. In our case this distribution should smoothly evolve with time from power to exponential (equilibrium) form, i.e., it should be a quasi-power like distribution. The most common distribution of this type is the Tsallis distribution [12][13][14], characterized by the energy dependent parameter q = q(t) (note that for q → 1 the Tsallis distribution becomes the BG from Eq. (6)). To justify such a choice let us recall a unique feature of distribution (9) which distinguishes it from all other distributions used so far in this context. Namely, whereas for z → ∞ (or for z >> 1/(q − 1)) it becomes a power distribution, for z → 0 (or for z << 1/(q − 1)) it goes into an exponential distribution. Thus, for an appropriate selection of the parameter q(t), it can describe all the distributions occurring in the RTA formula: f in (z,t) for q t f = 0 = q in > 1 and f eq for q t f → ∞ = 1. The effectiveness of the Tsallis distribution is best evidenced by the results of works [15][16][17]. In particular, as shown in [16], it nicely describes a wide range of the measured transverse momenta (0.1 < p T < 100 GeV which corresponds in this case to 1 < z < 700) in which the cross section spans a range of ∼ 14 orders of magnitude 3 .
The time-dependent parameter q = q(t) (more specifically, its deviation from q = 1) represents the degree of the non-extensivity or, in other words, the degree of deviation of the system from the thermalized or equilibrated system, which is usually described by the well known BG statistical mechanics. It is also a control parameter that fully defines the shape of the Tsallis distribution, in particular its evolution over time through moments such as the expected value, z(t) and variance, Var[z(t)]: .
In general, the moments < z n >, limited to n+1 < 1/(q−1), are related by the recurrence relation As mentioned above, we assume that the dynamic evolution of a system over time smoothly and monotonically transforms the probability distribution f (z,t) in (z,t)-space and that f (z,t) is a Tsallis distribution fully described by the time-dependent control parameter q(t).
Now note that if we had completely formally used the Tsallis distributions for all distributions in Eq. (5) defining the RTA and comparing f (z = 0) or < z >, we would get such a relationship between the parameters q appearing there (remembering that f eq is assumed as a BG distribution which is equavalent to a Tsallis distribution with q eq = 1): However, if from the very beginning we decide to describe the entire process using only quasi-power Tsallis distributions, the time evolution of which is given only by the time evolution of their control parameters q = q(t), which means that f (t) = f [q(t)], we should go back to Eq. (2), which is now This equation replaces Eq. (2). The form of the function F from Eq. (14) can be deduced by taking f (t) given by the Tsallis distribution with q = q(t) and calculating d f /dt. As a result, we get that . The distributions f in and f eq are the same as in Fig. 1 (given by the Tsallis formula (9) with q in = 1.25 and q eq = 1, respectively). and expressing the dependence of Q(z) on the variable z by To go further, we need to set the time dependence of the parameter q in some way. Note that the nonextensivity parameter q(t) describes deviations of the state of a statistical system from equilibrium and in this sense it plays the role of an internal variable discussed in Refs. [19,20]. Therefore, following such an approach [20] we assume that the equation of the dynamics describing the control parameter q(t) is of the form of the equation of a relaxation: Remembering that we always assume that q eq = 1, the solution of (17) is which coincides with Eq. (13). Fig. 2 shows the resultant schematic distributions f f o for different t f /τ; they all have the form of a Tsallis distribution with q = q t = t f as given by Eq. (18). As one can see the result is now different from that using the RTA approximation shown in Fig. 1. Notice that the relaxation times τ = τ f in Eq.(3) and τ = τ q in Eq.(17) describe the relaxation of different quantities, respectively the entire distribution and its control parameter. Comparing (for the same time) mean values z = 1 + exp − t τ f ( z in − 1) evaluated from the RTA, Eq. (7), and its q version, , Eq. (9) for q(t) given by Eq. (18), we have that and the ratio of relaxation times in both approaches is changing from τ f /τ q = 1/ < z > in for t → 0 to τ f /τ q = 1 for t → ∞ (cf. Fig 3). Note that in a situation where in some isolated system we have a fixed number of particles N and a fixed total energy U, we have a constant average energy E . Therefore, in such a case, the variability in time of z = E /T must be caused by the appropriate variability in time of the scale parameter T (here the temperature). This means that the scale parameter in our scaled variable z = x/x 0 also changes with time: where x 0 = x * 0 at equilibrium (t −→ ∞). Because, as known from [21], fluctuations of the scale parameter x 0 are directly connected with the parameter q: we can write that This therefore means that the relaxation time τ now describes the temporal evolution of the fluctuations of the scale parameter x 0 (in the scaled variable z = x/x 0 ).

Conclusions
A few remarks that may inspire further research in this direction could be of interest here [11]. Notice that in the language of information theory based on Shannon entropy, our time evolution Eq. (14) can be expressed as where the S z = − ln[ f (z)] is suprisal, which measures the information gained by observing the outcome z in the system and the Shannon entropy is its expectation value. Now note that for the Tsallis distribution (9) we have that S z ∼ z and ∂ S z /∂t ∼ z. The entropy rate is given by From [11] we know that the linear relationship between ∂ S z /∂t and z guaranties that distribution (9) saturates the time-information uncertainty bound: More precisely, for the Tsallis distribution because ∂ S z /∂t ·z = 0 and (for any distribution) ∂ S z /∂t = 0. The distance of a given distribution (in our case defined by q(t)) from the equilibrium distribution (defined by q = 1) is given by the difference of the corresponding entropies: From equations (16), (17) and (29) we get that In conclusion, we propose a new, modified form for the relaxation time approximation for the collision term in the Boltzmann equation (2) allowing a smooth transition to the thermalized distribution. It consists of replacing the simple form of this term, given by the Eq. (3), where the relaxation time τ determines how fast the equilibrium state of the studied distribution f is reached, by Eq. (17) describing the time evolution of the most important (control) parameter of the analyzed distribution. The relaxation time τ would now control the rate of change of this parameter from some value to one that corresponds to the equilibrium state. We argue that this is possible if we use for the phenomenological description of the distributions of interest the quasi-power law Tsallis distribution given by Eq. (9) which is able to describe a given process at all stages of time evolution. Its control parameter is the time-dependent non-extensivity parameter q(t) and the relaxation time parameter τ describes its time evolution, as shown in Eq. (17). This single quasi-power law probability distribution (9) smoothly evolving towards thermalization would then replace the two-component distribution given by Eq. (7) which arises from the RTA. The proposed scheme offers multiple applications in situations where one wants to study the time evolution of an ensemble but one does not want to invoke the kinetic theory with complicated collision integrals.