Appendix A: Viscosity coefficient ξ ( T ) from LQCD Following the discussion presented in [34], the bulk viscosity of QGP can be calculated from the lattice QCD by Eq. (13) in that paper. We assume that the decay factors for pions and kaons are vanishing above the critical temperature of the phase transition QGP-hadrons. The quark-antiquark condensates can be neglected at temperatures higher than the critical one [35]. Therefore, Eq. (22) of Ref. [34] would be reduced to Our algorithm is the following. Using lattice QCD results on trace anomaly, ( ϵ − 3 p ) /T $^{4}$, and other thermodynamical quantities, we can determine the bulk viscosity. To make use of the lattice QCD results, it is useful to make a suitable fit to the data at high temperatures. Then we obtain the following equations of state where ω = 0 . 319, β = 0 . 718 ± 0 . 054 and r = 0 . 23 ± 0 . 196. Using the equations of state, Eq. (49) in Eq. (47), we obtain For analytical purposes, the function g ( z ), which is defined in z parameter as g ( z ) = F$_{0}$/F$_{1}$ in Eq. (16), can be numerically estimated depending on the parameter z by using the following procedure. First, we plot it parametrically depending on the parameter H , Fig. 17 Appendix B: Estimations of g ( z ) where ρ is the energy density and c 2 s = dp/dρ is the square of the speed of sound. The parameter ω$_{o}$ is a scale depending on the temperature T , and defines the validity of the underlying perturbation theory. In this relation, the viscosity is assumed to have a thermal part which can be determined through lattice calculations, and a vacuum contributing part, which can be fixed using quark and gluon condensates. The vacuum part would take the value 9 ω$_{0}$ ξ = T s ( 1 c 2 s − 3 ) − 4( ρ − 3 p ) + 16 | ϵ$_{v}$ | (47) 16 | ϵ$_{v}$ | (1 + 3 8 · 1 . 6) ≃ (560 MeV) 4 ≃ (3 T$_{c}$ ) 4 (48) ξ ( ϵ ) = 1 9 ω$_{o}$ 9 γ 2 − 24 γ + 16 γ − 1 ρ + 9 ω$_{o}$ T 4 c . (49) p = ωρ, T = βρ $^{r}$, c 2 s = ω