Dissipative properties and isothermal compressibility of hot and dense hadron gas using non-extensive statistics

We evaluate the transport properties such as shear viscosity (η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}), bulk viscosity (ζ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\zeta $$\end{document}) and their ratios over entropy density (s) for hadronic matter using relativistic non-extensive Boltzmann transport equation (NBTE) in relaxation time approximation (RTA). In NBTE, we argue that the system far from equilibrium may not reach to an equilibrium described by extensive (Boltzmann-Gibbs (BG)) statistics but to a q-equilibrium defined by Tsallis non-extensive statistics after subsequent evolution, where q denotes the degree of non-extensivity. We observe that η/s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta /s$$\end{document} and ζ/s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\zeta /s$$\end{document} decrease rapidly with temperature (T) for various q-values. As q increases, the magnitudes of η/s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta /s$$\end{document} and ζ/s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\zeta /s$$\end{document} decrease with T. We also show the upper mass cutoff dependence of these ratios for a particular q and find that they decrease with the increase in mass cutoff of hadrons. Further, we present the first estimation of isothermal compressibility (κT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa _T$$\end{document}) using non-extensive Tsallis statistics at finite baryon chemical potential (μB\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _B$$\end{document}). It is observed that, κT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa _T$$\end{document} changes significantly with the degree of non-extensivity. We also study the squared speed of sound (cs2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c_{s}^{2}$$\end{document}) as a function of temperature at finite baryon chemical potential for various q and upper mass cutoffs. It is noticed that there is a strong impact of q and mass cutoff on the behaviour of cs2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c_{s}^{2}$$\end{document}.


Introduction
High-energy heavy-ion collisions provide an opportunity to probe the nuclear matter under extreme conditions i.e. at high temperature and/or energy density. The space-time evolution of a system formed in heavy-ion collisions is governed by its dissipative properties such as shear viscosity (η), bulk viscosity (ζ ) etc. The spatial anisotropy created in heavyion collisions gets converted to the momentum anisotropy of the produced particles due to the pressure gradient. The a e-mail: Raghunath.Sahoo@cern.ch equilibration of momentum anisotropy is governed by the shear viscosity. The AdS/CFT calculation shows that there is a lower bound to the value of shear viscosity to entropy density ratio (η/s) for any fluid found in nature. The lower bound (known as the Kovtun-Son-Starinets (KSS) bound) has been set to the value of 1/4π in natural units [1,2]. The elliptic flow measurements at Relativistic Heavy-Ion Collider (RHIC) experiment have found that the η/s is close to the KSS bound, which developed intense interest in this ratio of the strongly interacting matter described by quantum chromodynamics (QCD) [3,4]. Also, a peak in Bulk viscosity to entropy density ratio (ζ /s) is expected near QCD critical temperature (T c ), where conformal symmetry breaking might be significant as expected by various effective models [5][6][7]. Other interesting thermodynamic properties are the isothermal compressibility (κ T ) and speed of sound (c s ), which are used to define the equation of state of the system and quantify the softest point of the phase transition along with the location of the critical point [8].
Due to high multiplicities produced in high-energy collisions, the statistical models are more suitable to describe the particle production mechanism and to study the QCD thermodynamics. Such a statistical description of transverse momentum ( p T ) of final state particles produced in highenergy collisions has been proposed to follow a thermalized Boltzmann-Gibbs (BG) distribution. But, a finite degree of deviation from the equilibrium statistical description of p T spectra has been observed by experiments at RHIC [9,10] and Large Hadron Collider (LHC) [11][12][13][14]. In addition, the matter produced in high energy collisions evolves rapidly in a non-homogeneous way. Hence, the spatial configuration of the matter becomes non-uniform and the global equilibrium is not established [15][16][17][18]. As a consequence of which some of the observables become non-extensive and develop a power-law tail in the spectra instead of exponential distributions. This also happens when there is local temperature fluctuation and long-range correlations in the produced system [19]. In these cases, the use of BG distribution is questionable [20][21][22][23][24]. Systems, where the usual ergodicity is violated, such as, the states in large systems involving longrange forces between particles and metastable states in small systems where the number of particles are relatively smaller, a generalized BG entropy has been introduced by C. Tsallis [25,26]. Recently, a growing attention has been paid towards the possible non-extensive effects in thermodynamics and statistical mechanics [27,28] as well as to explain the particle spectra in high-energy hadronic and heavy-ion collisions [29][30][31][32][33][34]. The nuclear modification factor [35,36] and elliptic flow of identified particles [37] have also been successfully described by using Tsallis non-extensive statistics in Boltzmann Transport equation (BTE) with relaxation time approximation (RTA).
The aim of the present paper is to study various dissipative properties such as shear and bulk viscosities using the relativistic non-extensive Boltzmann transport equation (NBTE), where we employ the RTA for the collision term. In NBTE, we assume that a non-equilibrium system, which dissipates energy and produces entropy, measured by the non-extensive parameter q, relaxes to a local q-equilibrium after a certain relaxation time τ . The hadrons are treated as classical Boltzmann particles and the elastic scattering processes are taken into account. A similar approach has been employed in Ref. [38], where an extensive statistics is used in BTE to calculate the transport coefficients. Here, we extend the approach using Tsallis non-extensive statistics and calculate shear viscosity, bulk viscosity and their ratios over entropy density. First time, we quantify the isothermal compressibility for a hadronic matter using non-extensive statistics. We also calculate the squared speed of sound of hadron resonance gas at finite baryon chemical potential. In view of high-multiplicity events at the LHC and the observation of increase of particle multiplicity driving the system towards thermodynamic equilibrium [39], it becomes important to study these dissipative properties of the systems which are away from equilibrium.
The paper runs as follows. In Sect. 2, the dissipative properties such as shear and bulk viscosities are derived using the relaxation time approximation of the relativistic NBTE. The formulation of isothermal compressibility and squared speed of sound in non-extensive statistics are also given. In Sect. 3, the results and discussions are presented. Finally, we summarize with the findings of this work in Sect. 4.

Formulation
We follow the approach mentioned in Ref. [38] to calculate the dissipative properties such as shear and bulk viscosity for a hadronic matter using non-extensive statistics. We start with the BTE given by, where v i p is the velocity of i th particle and F i p is an external force acting on i th particle. I ( f p ) is the collision integral which gives the rate of change of the non-equilibrium distribution function f p when the system approaches qequilibrium.
Assuming no external force and proceeding with the RTA, the collision integral can be approximated as, where τ (E p ) is the relaxation time or collision time. We take non-extensive Tsallis distribution as f 0 p [40] near the local rest frame of the fluid, where the system is described locally by T , μ B and fluid velocity, u, which change slowly in space and time [41]. The thermodynamically consistent Tsallis distribution, ( f 0 p ) [42] in the Boltzmann's approximation is given as, where u is the fluid velocity. T and μ are temperature and chemical potential, respectively. μ = bμ B + sμ s , where b and s are baryon and strangeness quantum numbers, respectively. μ B and μ s are baryon and strange chemical potentials. For the sake of simplicity, we do not consider the strangeness neutrality condition. Now, the stress-energy tensor (T μν ) can be written as, where T μν 0 is the ideal part and T μν dissi is the dissipative part of the stress-energy tensor. In the hydrodynamical description of QCD, shear and bulk viscosities enter in the dissipative part of the stress-energy tensor, which can be written (in the local Lorentz frame) as [41], In terms of distribution function, this can be expressed as, where δ f p is the deviation of the distribution function from the q-equilibrium and is given by (from Eqs. 1 and 2 ), Assuming a steady flow of the form u i = (u x (y), 0, 0) and space-time independent temperature, Eq. 5 simplifies to T xy = −η∂u x /∂ y. Now, from Eqs. 6 and 7 , we get (using μ = 0), Thus the coefficient of shear viscosity for a single component of hadronic matter can be expressed as, The bulk viscosity (ζ ) is related to the dissipation when the system is uniformly compressed. We take the trace of Eq. 5 and get, Using Eq. 6, the above equation modifies as, Now, using the conservation law for energy-momentum i.e. ∂ μ T μν = 0, one can obtain [41], and where, w = ε + P is the enthalpy density, ε is the energy density and P is the pressure. Here, we cannot take steady flow. Thus, Eq. 7 takes the form, where C V is the heat capacity at constant volume. Using Eqs. 11 and 14 and manipulating with the help of conservation laws, the expression for ζ is written as [41], Now, using d P = sdT = wdT /T , where s is the entropy density, the above equation is reduced to, Here, c 2 s = s/C V is the squared speed of sound at constant baryon density. For a multi-component hadron gas at finite chemical potential, the shear and bulk viscosities are written as, and, Here E 2 a = p 2 + m 2 a and the barred quantities represent the antiparticles. f 0 a is the distribution function for a th particle. Now, the energy dependent relaxation time is given as, where W (a, b → c, d) is the transition rate defined as, Here |M| is the transition amplitude. In the center-of-mass frame, Eq. 19 can be simplified as, where v ab is the relative velocity and √ s is the center-of-mass energy. σ ab is the total scattering cross-section in the process For further simplification, τ (E a ) can be approximated to averaged relaxation time ( τ ) [43] and it can be obtained from Eq. 21 by averaging over f 0 a as, is the number density of b th hadronic species. The thermal average for the scattering of same species of particles with constant cross-section at zero baryon density can be calculated as follows [38,44,45].
Here, e (x) q is the q-exponential which is defined as e (x) q = [1+(q −1)x] q/(q−1) . The momentum space volume elements can be written as, The numerator in Eq. 23 is written as, and the denominator is written as, Now, using the generalized Tsallis non-extensive statistics for q-equilibrium, the σ ab v ab takes the following form, Here σ is the cross-section used as a parameter in the calculations. E a and E b are integrated in the limit m a to ∞ and m b to ∞, respectively. The limit of integration for cos θ is -1 to 1. The relaxation time is calculated by using Eqs. 22 and 27 . The other thermodynamical quantities in non-extensive statistics are calculated as [42], n, and P are the number density, energy density and pressure of hadrons, respectively. The non-extensive entropy den-sity, s can be calculated from the above expression as, We use the basic expression for addition of entropies in nonextensive statistics while calculating for the multi-component hadron gas, which is given by, where s(A + B) is the total entropy of A and B. s(A) and s(B) are the entropies of A and B, respectively. We have also calculated the isothermal compressibility for hadron gas using non-extensive statistics. The isothermal compressibility (κ T ) is defined as [46], where V is the volume of the system. Again in terms of fluctuation and average number, isothermal compressibility can be expressed as [46,47], Since we use the thermodynamically consistent Tsallis statistics, the above thermodynamical relation is valid in this case. Using the basic thermodynamical relation where ∂n q ∂μ is given as, In hydrodynamics, speed of sound plays an important role in understanding the Equation of State (EoS) and hence, the associated phase transition. This is because of the fact that speed of sound depends on the properties of the medium. The squared speed of sound (c 2 s ) is defined as, The expression of c 2 s at finite baryon chemical potential is written as [48], The term dμ B dT can be calculated using the constant entropy per particle (s/n) condition [48,49], and given as,

Results and discussions
In this section, we present the results and discussions of the dissipative and thermodynamic properties of hadronic matter calculated in the NBTE using relaxation time approximation, where non-extensive statistics is used as the q-equilibrium distribution function. The non-extensive parameter, q represents the degree of non-extensivity or in other words the tendency of the system to stay away from the thermodynamic equilibrium. The relaxation time of the system is thus correlated with the q-parameter. The non-extensive parameter, q is related to local temperature fluctuation and possible longrange correlation in the system [19], which is also responsible for the deviation of the system being away from the thermodynamic equilibrium. It has been reported that in the hadronic and elementary collisions where the spectral analysis shows a deviation of the system from the BG equilibrium, the q-value is 1.22 [50]. We study the response of the system as a function of q-parameter. The present study will be very important in view of the multiplicity dependent analysis done in p + p collisions at the LHC. We have seen explicitly from the identified particle spectra in p + p collisions that with higher multiplicity classes, the q-value goes towards one [39]. We consider all hadrons and their resonances with a mass cutoff of 2.0 GeV. First, we estimate various dissipative properties of hadron gas such as shear viscosity (η), bulk viscosity (ζ ) and their ratios to entropy density at zero baryon chemical potential (μ B = 0). We use Eqs. 17 and 31 to estimate η and s, respectively for zero baryon chemical potential. Here, we calculate the relaxation time (τ ) of a th hadron by using Eq. 22, where a constant hadronic collision cross-section (σ ) = 11.3 mb is taken [38]. In order to check the validity of formulations discussed above, we have presented the calculation of η/τ for pion and various mass cutoffs at q = 1.1 in Fig. 1. We observe that η/τ initially varies very slowly with temperature and as temperature increases, it increases rapidly. The results for pion exactly matches with those of presented in Ref. [2], where the degeneracy (g) of pion is taken to be 1. As we add more hadrons and their res-  Fig. 2 (η/s)/τ versus temperature for various q-values at zero baryon chemical potential onances into the system, η/τ increases rapidly at a lower T . These findings suggest that the shear viscosity becomes large with the increasing number of particles into the system. We have also shown the variations of the ratio (η/s)/τ of pion with T at various q-values in Fig. 2. Again, we find the same behaviour as observed in Ref. [2] i.e. there is almost negligible effect of the non-extensive parameter on this ratio particularly at higher temperature. This is due to the fact that entropy density (s) also increases with the same amount as η/τ with q.
The results of shear viscosity (η) of hadrons as a function of temperature (T ) are presented in the top panel of Fig. 3 for various q-values. We find that η increases with temperature in a similar way as observed in Ref. [38], where Excluded-Volume Model with extensive statistics is used. η has smaller values for lower q, which goes in line with the observation that the values of the transport coefficients become smaller as the system goes towards equilibrium. Mid- The general behaviour of η/s is similar to that observed in Ref. [38]. This ratio decreases with T at all q-values due to rapid increase in entropy density.
As the degree of non-extensivity, q increases, the magnitude of η/s decreases and goes below the KSS bound at a higher T . The temperature at which η/s goes below the KSS bound decreases with increasing q. This can be understood as the drastic change in entropy density of multicomponent hadron gas with T , which makes the value of the ratio, η/s very small. This explanation has been given as a counterexample to the KSS bound in Ref. [51]. Also, in certain field theories [52][53][54] KSS bound has been shown to be violated. This is a very important finding while understanding the dissipative properties, non-extensivity and the effect of temperature of the system. For comparison, we have also shown the calculation of η/s for hadron gas at zero baryon chemical potential using Boltzmann-Gibbs distribution function. This is shown by the grey solid line in the middle panel of Fig. 3. η/s as a function of temperature is presented in the lower panel of the figure at q = 1.07 for various upper mass cutoff of hadrons. Study of the response of the system as a function of temperature for various mass cutoffs is a requirement at higher collision energies. With increase of collision energies heavier resonances are produced. We find that the magnitude of η/s decreases with the increase of mass cutoff due to the contributions of heavier resonances in entropy density of the system.
In Fig. 4, we have shown the comparison of η/s as a function of T obtained in different model calculations. Here, the red solid line is η/s for hadron gas calculated using parton-hadron-string dynamics (PHSD) [55], where RTA  with Boltzmann-Gibbs distribution as equilibrium distribution function is used. The blue line is the result taken from Ref. [38], where BTE with RTA is used and BG distribution is taken as the solution of BTE. The red dotted and black dash-dotted lines are the estimations of the present work for q = 1.01 and for extensive case (q = 1), respectively. The difference between the PHSD and our results for q = 1.01(∼ 1) is due to the difference in estimation of relaxation time. Although we have followed a similar procedure as done in Ref. [38] using non-extensive distribution function, for simplicity, we do not include the excluded-volume correction in the estimation of thermodynamical observables such as number density, entropy density etc., as is done in Ref. [38]. We also see from the figure that when q→1, the results are closer to those obtained for extensive case. In Fig. 5, we discuss the results of bulk viscosity (ζ ) and its ratio to entropy density (s) calculated by using Eqs. 18 and 31 at zero baryon chemical potential. Again, we take a constant value of hadronic cross-section (σ ) equal to 11.3 mb while calculating the relaxation time (τ ). The upper panel represents ζ versus T plot for different values of non-extensive parameter, q. We observe that ζ increases slightly with T for a particular q, but as q increases it starts saturating with the increasing T . In the middle panel, we show the variations of bulk viscosity to entropy density ratio (ζ /s) at μ B = 0 for various q-values. We find that ζ /s decreases rapidly with T for all q. The value of ζ /s is lower for higher q-values. We also observe that ζ /s changes rapidly at a particular T , which is found to be lower for higher q-values. These findings suggest that the degree of approach towards an asymptotic value is sensitive to the non-extensive parameter, q. This goes inline with the expectations of dissipative properties of a system approaching thermal equilibrium. In the lower panel, we show the upper mass cutoff dependence of ζ /s for q = 1.07. We observe a similar effect of mass cutoff on this ratio as observed in the case of η/s i.e. as we increase the upper mass cutoff of hadron resonance gas its magnitude decreases, which is due to increase of entropy density after adding the heavier resonances in the system.
Next, we discuss one of the important thermodynamic observables i.e. isothermal compressibility (κ T ), which has drawn considerable interest in these days. It is very interesting to quantify this observable for a hot and dense hadron gas formed in hadronic and heavy-ion collisions. In Ref. [8], κ T is first time estimated for hadrons produced in heavyion collisions in extensive statistics using Hadron Resonance Gas (HRG) model. In this work, we have used non-extensive statistics for the calculation of κ T . In Fig. 6, results of κ T as a function of temperature at μ B = 0 are shown. It is observed that κ T initially decreases rapidly at low temperatures and becomes constant at higher temperatures for all the q-values. The temperature at which κ T becomes constant depends on the q-value. As q increases the value of κ T decrease and  the temperature at which it becomes constant shifts towards the lower values, which suggests that there is a strong dependence of κ T on non-extensivity. Figure 7 describes the dependence of isothermal compressibility on temperature at μ B = 0.5 GeV for various q. We find a similar behaviour of κ T in this case as observed for zero baryon chemical potential but comparatively has smaller values.
Further, in order to see the effect of upper mass cutoff of hadrons on κ T , we proceed as follows. We calculate κ T at μ B = 0.5 GeV for q = 1.07, which is shown in the Fig. 8 as a function of temperature. Again, a similar trend is obtained over all the temperatures as observed in the previous plots and as the upper mass cutoff increases, values of κ T decrease. This effect is more pronounced upto mass cutoff of 1 GeV and after that it gets weaker. κ T becomes independent of system temperature when higher resonances are added to the system. Hence a system is driven towards a critical behaviour, once we have resonances of higher masses as the ingredients of the system. Another important thermodynamical observable is speed of sound. We use Tsallis statistics to calculate the squared speed of sound (c 2 s ) for finite baryon chemical potential. In Ref. [56], the temperature and mass cutoff dependence of c 2 s for various q-values have been studied in Tsallis statistics. It is noticed that the criticality in c 2 s , which is seen as the q-dependent peak in c 2 s , when studied as a function of temperature, disappears after q = 1.13. We have shown c 2 s as a function of temperature at μ B = 0.5 GeV for different qvalues in Fig. 9. We find that c 2 s initially increases with T and after reaching at a maximum value it starts decreasing. The maximum in this observable, which is related to the criticality shifts towards lower T as q increases and disappears after q = 1.13. These findings are same as observed for the case of μ B = 0 in Ref. [56]. In Fig. 10, the mass cutoff dependence of c 2 s is shown for q = 1.13 and μ B = 0.5 GeV. The criticality in c 2 s disappears when lower mass cutoff is reached, which advocates that criticality of the system depends on the number of hadrons present in the system. Also, we find that c 2 s decreases with the mass cutoff, which goes in line with the previous observations [56,57].

Summary
In summary, we have studied the transport coefficients using the relativistic non-extensive Boltzmann transport equation (NBTE), where we employ the relaxation time approximation for the collision integral. Here, we have extended the approach using Tsallis non-extensive statistics as qequilibrium distribution function and calculated shear and bulk viscosity to the entropy density ratios. We have also calculated isothermal compressibility and squared speed of sound using Tsallis statistics at finite baryon chemical potential for different q. The important findings of this work are summarized as follows: 1. We have found that the ratios of shear and bulk viscosity to entropy density have a strong dependence on nonextensivity. Their values decrease with the increasing qparameter with the temperature. 2. We have observed the effect of mass cutoff on these ratios and as mass cutoff increases the magnitudes of these ratios decrease. These effects are more pronounced at higher values of temperature. 3. We have first time studied the effect of non-extensive parameter, q on isothermal compressibility. We have studied the temperature dependence of isothermal compressibility for various q-values and observed that it decreases rapidly with temperature and becomes constant at higher temperatures. The temperature at which the isothermal compressibility shows a saturation, may indicate an appearance of criticality in the system, which is dependent on q and shifts towards lower values of temperatures for higher q-values. 4. While studying the temperature dependence of isothermal compressibility for various q-values, it has been found that the values of κ T are smaller for higher q. It has also been found that, its values decrease at finite baryon chemical potential in comparison to that obtained at zero baryon chemical potential. 5. We have found that, isothermal compressibility strongly depends on the upper mass cutoff of hadron resonances and decreases with the increasing mass cutoffs. It becomes independent of system temperature when higher resonances are added to it. Hence a system is driven towards a critical behaviour, once we have resonances of higher masses. 6. It is observed that as the non-extensivity of the system increases the criticality observed as a peak in the squared speed of sound versus T plot shifts towards lower temperatures and disappears after q = 1.13. We have also noticed a strong dependence of criticality on mass cutoff. As upper mass cutoff decreases, the criticality of the system disappears.