Semiclassical approximation meets Keldysh–Schwinger diagrammatic technique: scalar φ4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi ^4$$\end{document}

We study the evolution of the non-equilibrium quantum fields from a highly excited initial state in two approaches: the standard Keldysh–Schwinger diagram technique and the semiclassical expansion. We demonstrate explicitly that these two approaches coincide if the coupling constant g and the Plank constant ħ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbar $$\end{document} are simultaneously small. Also, we discuss loop diagrams of the perturbative approach, which are summed up by the leading order term of the semiclassical expansion. As an example, we consider shear viscosity for the scalar field theory at the leading semiclassical order. We introduce the new technique that unifies both semiclassical and diagrammatic approaches and open the possibility to perform the resummation of the semiclassical contributions.


Introduction
Highly nonequilibrium dense quantum fields define the initial stage of many physical problems. These include the physics of the early stage of ultrarelativistic heavy ion collisions [1][2][3], cold atomic gases [4,5] and the processes in the early Universe [6][7][8][9][10][11]. At present, there are a variety of approaches used for the description of the quantum field evolution from a highly excited initial state to the quasistationary one, where hydrodynamic equations work effectively.
One of the most advanced approaches is the Keldysh-Schwinger diagram technique which provides a systematic way of studying nonequilibrium phenomena in quantum field theory [3,[12][13][14]. With the help of this technique, one can derive the kinetic equations that describe the evolution of quasiparticle distribution function and observables consequently. Also, this technique can be used for the systematic evaluation of thermodynamical and transport properties of the quantum systems at the thermal equilibrium [3,14]. a e-mail: raan@lpi.ru (corresponding author) b e-mail: semenov@lpi.ru Another way to deal with the nonequilibrium initial state comes from a physical intuition and is based on the assumption that at high energies and/or high occupation numbers the dynamics of the quantum fields is semiclassical, so the classical equations of motion can be used [15][16][17][18][19][20][21][22][23][24][25]. In order to complete this approach, one should make additional assumptions about an ensemble which is used for the averaging of observables. In the literature, this approach is often called the Classical Statistical Approximation (CSA). It is very useful for numerical simulations since it allows numerically extract (nonperturbatively in coupling constant) results for observables [26][27][28] and transport coefficients [29]. In the previous works of one of the authors [20][21][22], it was shown that the CSA arises as the leading order of the semiclassical expansion.
This work aims to demonstrate that the Keldysh-Schwinger diagram technique and the classical statistical approach are two facets of one general way to deal with nonequilibrium quantum fields. It seems that these two approaches are quite different from a practical point of view. The Keldysh-Schwinger diagram technique originates from the perturbative expansion in the coupling constant g. In order to evaluate the observable consistently in this approach, one should (with the help of the Wick theorem and the diagram technique) derive the quantum kinetic equation on the distribution function, solve it, and evaluate observable using the solution. Conversely, the CSA comes from theh expansion. In order to find the observable there, it is necessary to solve the classical equations of motion and to average the observable on the classical trajectories with the ensemble of the initial conditions. In this work, we show that these two approaches can be derived from one general path-integral representation. We demonstrate which assumptions and approximations should be made to obtain one approach or the other. We study the case of both g andh are small, where these two approaches should be consistent, and analyse how the leading order term of the semiclassical expansion (the CSA) sums up certain multi-loop diagrams of the perturbative Keldysh-Schwinger technique.
Previously some efforts were made in this direction. In particular, in works [30,31], the comparison between the Keldysh-Schwinger diagram technique and the classical approximation was performed for the thermal equilibrium, and the agreement of the leading order contributions at high temperatures was demonstrated.
The paper is organised as follows: Starting from the general setup described in Sect. 2, we briefly review the standard Keldysh-Schwinger diagram technique in Sect. 3 and the semiclassical approximation in Sect. 4. In Sect. 5 we perform the detailed comparison of these approaches in the limit where g andh are small simultaneously and analyse the loop contributions. In Sect. 6, on the example of the shear viscosity, we show how to employ the CSA for the calculation of the relevant physical observables and discuss the applicability of the semiclassical approach. Section 7 is devoted to a new diagram technique which naturally combines both approaches described above. The discussion of results and conclusions are in Sect. 8.

Keldysh-Schwinger approach to the non-equilibrium QFT
The standard way to deal with the non-equilibrium quantum field theory includes the Keldysh-Schwinger technique, also known as the closed-time path formalism [3,14,23]. In this approach, averages are calculated as the trace with the density matrix operator. Time evolution of the density matrix is defined by two evolution operators; that is why the doubling of the degrees of freedom occurs. Moreover, the initial density matrix should be additionally defined from the physics of the considered system. In this work we consider the scalar field theory with the action: Here and after we use mostly minus metric convention g μν = (+, −, −, −) and x μ = (t, x). Using Eq. (1) one can calculate the Keldysh action as a difference between actions on the forward and the backward parts of the Keldysh contour. Averages are expressed through the path integrals with the Keldysh action as [3,14,23]: Here and after we keeph explicitly in order to study the semiclassical limit of the theory.
It is convenient to rotate the basis of ϕ F , ϕ B fields to the socalled "classical" and the "quantum" fields (there are equivalent notations for such rotation in the literature φ cl ≡ φ r and φ q ≡ φ a ) [3,14,23]: This new ϕ cl , ϕ q basis has several advantages: the causality of the theory become explicit, the vertices look simpler, and the semiclassical limit is transparent. Then the Keldysh action transforms to (after integration by parts): We keep explicit dependence on the initial time t 0 to take into account highly non-equilibrium initial states. Usually, the initial time moment is set to the past infinity, and the boundary term S init [ϕ cl , ϕ q ] is dropped out. The averages can be calculated by integration over new fields as: Note, ifÔ contains only equal-time operators, then it is sufficient to keep only the ϕ cl component in the integrand of the expression Eq. (5) due to causality. The expressions similar to Eqs. (5) and (2) can be found in many modern textbooks that discuss the Keldysh-Schwinger technique. However, such representation is somewhat misleading since it does not contain information about the initial state of the theory, which makes ϕ F and ϕ B fields correlated. More rigorously the Eq. (5) can be written as [20][21][22]32]: where the integral with i.c. means the initial values for the ϕ cl field, ϕ cl (t 0 , x) = α(x), ∂ t ϕ cl (t 0 , x) = Π(x); whereas the initial values for the ϕ q are not fixed. S init [ϕ cl , ϕ q ] is then taken into account and absorbed by the Wigner functional.
The Wigner functional is related to the initial value of the density matrix operatorρ(t 0 ) as: This functional contains all the information about the initial state of the system. The Eq. (6) represents the general expression from which one can deduce both the perturbative and the semiclassical approaches which were discussed in the introduction. In the next section, we derive the standard Keldysh-Schwinger perturbation technique and discuss its limitations.

Standard perturbative approach
The standard Keldysh-Schwinger diagram technique follows naturally from two major assumptions: -The Gaussian form of the initial Wigner functional that allows the Wick theorem to be valid. -The possibility of the perturbative expansion.
Under these assumptions, the Eq. (6) can be rewritten as: where the averaging over the noninteracting fields · · · 0 should be performed with help of the Wick's theorem with four basic contractions [3,14,23]: Here G 0 R(A) is retarded (advanced) Green's functions which can be equivalently defined in the operator formalism as In the absence of the interactions these free correlators Eq. (10) are independent from the initial Wigner functional and solve the equation: with the corresponding boundary conditions (retarded or advanced ones). The solution of the Eq. (11) is, for example: The Keldysh Green's function: solves the equation: This correlator depends crucially on the initial state of the system and cannot be found without specification of the initial density operator. In the simplest case, the initial state of the system is characterised by the one-particle distribution function f p , which is related to the Keldysh Green's function as: [3,14,23] It is necessary to stress here that only for the Gaussian initial state the knowledge of G 0 is enough to perturbatively evaluate the average of any product of the free fields and build up the diagram technique. 1 The basic elements of each diagram are two propagators: and two vertices: Here the "black" and the "white" vertices differ by the power ofh 2 . It is specialised for the exact comparison with the semiclassical approach later. For example, let us draw diagrams for the first two orders of the coupling constant expansion of the full retarded Green's function in the presence of interactions 1 Another way to include the initial conditions is to extend the Keldysh contour onto the imaginary axis to take into account the Matsubara part [14]. However, it works only for the special case of the thermal initial state.
and make some important observations. The first one is related to causality. The zeroth-order retarded Green's func- It means that time increases according to the direction of the arrow on the diagrams and each diagram in this expansion vanishes identically if t ≤ t . So, this diagrammatic expansion respects causality and full Green's function G R (x; x ) = 0 for t ≤ t as expected. The another observation, let us cut all G 0 K lines and inspect what remains. It can be seen that the number of remaining loops exactly equal to theh order of the diagram, i.e. twice of the number of the "white" vertices. We prove this statement Sect. 7. In this place, the Keldysh-Schwinger technique differs drastically from the standard Feynman diagram technique, where theh order of any diagram coincides with the number of the loops. 2 The notation with the arrows, which is useful for the diagrammatic, may seem somewhat misleading in calculating the number of legs with and without arrows in each vertex of a diagram. The important point is that the retarded propagator G R (x; x ) connects two different fields, φ cl and φ q . The arrow in the notation belongs to the φ q field (that is why arrows are moved to one end of the line in the diagrams). In order to check a diagram, it is sufficient to calculate the number of outgoing lines in each vertex. It should be either one (the black vertex) or three (the white vertex). One can find an alternative notation for this theory in the work [31].
As an example, the explicit expression for the "cactus" diagram (the last diagram on the second line of Eq. (17)) can be written: In the next section we demonstrate how this diagram (and all others) originate from the coupling constant expansion of the CSA.

Semiclassical approach
In order to construct the semiclassical expansion, we add an auxiliary source J (x) to the theory described by Eq. (1): The source J (x) is used for the intermediate steps only and should be set to zero at the end of the calculations. Let us rewrite the Keldysh action (Eq. (4)) in a more convenient form: There are two key features of this action: -A[ϕ cl ] = 0 corresponds to projecting onto the classical equation of motion of the Lagrangian (Eq. (19)). -As far as we explicitly keeph-dependence, it is clear that the semiclassical approach is, in fact, the expansion of the last term:

Classical statistical approximation
The Leading Order term of the semiclassical expansion (Eq. (21)) is also known as the Classical Statistical Approximation, or the classical approach. In this case, the integration over ϕ q and ϕ cl fields can be done, and the Eq. (6) reproduce the well-known result [15][16][17][18][19][20][21][22]: where φ c is the solution of the classical equation of motion: with the initial conditions given by: and at zero axillary source J (t, x). Hence, the recipe for the CSA is the following: -find the classical trajectory as a function of the initial conditions. -calculate observables on this trajectory. -average over the initial conditions with the Wigner functional corresponding the considered problem.
Let us introduce new notation for averaging over initial conditions with the Wigner functional as: Then the definition of the CSA approximation (Eq. (22)) can be rewritten in this notation as It may seem that in the semiclassical expansion there are no linear inh contributions. However, it is not the case since the Wigner functional may depend on theh explicitly and averaging over the initial conditions may produce these terms [24,25]. For example, for the initial thermal state with the Bose distribution function and in the absence of the interactions the Kedlysh Green function, which is G 0 K ∼ φ c φ c i.c. , contains the combination (see Eq. (16)): In the zero temperature limit G 0 K ∼h/ω p which is linear inh, whereas for high temperature G 0 K ∼ T /ω 2 p and this contribution is purely classical and independent fromh.

Quantum corrections
The quantum corrections to the CSA (or the next-to-leading order of the semiclassical expansion) can be found with the help of the second term of the expansion (Eq. (21)). The integration over ϕ q can not be performed straightforwardly because of the new ϕ 3 q term. 3 However, each ϕ q can be replaced by the functional derivative over the source J due to ϕ q J term in the Keldysh action (Eq. (19)) as: 3 The same problem arise during the calculation of the correlation function like G R (x, x ).
Then the quantum corrections to the CSA averages are 4 The recipe of Eq. (28) similar to the recipe for the CSA: -find the classical trajectory as a function of the initial conditions. -perform three variations over the auxiliary source (not really needed). It is easy to recast all terms of the semiclassical approximation to the following general form: HereT denotes the anti-time ordering which is required to recover exponential form. The Eq. (29) shows that the building block of the semiclassical expansion is the full nonperturbative solution of the classical EoM φ c (x) and its variations over the additional source J (x).
It turns out that it is not necessary to calculate the variations of the classical solution explicitly. Let us define n-th variation as: Φ n (x; x 1 , x 2 , . . . x n ) can be calculated by variation of the classical equation of motion: Hence, to calculate the quantum correction to the CSA one need to find the solution of the n coupled differential equa-tions without knowledge of the exact dependence of the classical solution φ c (x) from the auxiliary source J (x). The initial conditions for these equations are zero because of the causality (φ c (x) depends on the source J only at the preceding times).

Comparison g 2 andh 2 expansions
Now we are ready to compare the perturbative and the semiclassical approaches up to two loops. For this purpose, we perform the semiclassical expansion of the G R (x 1 , x 2 ) up tō h 2 terms and show how this result, being decomposed further up to g 2 , reproduces the perturbative answer of the Sect. 3. Let us consider the full retarded Green's function and expand it according to Eqs. (21), (27), and (30) Let us denote the Leading Order retarded Green function as G C S A R (x 1 , x 2 ), then from Eq. (33) it is obvious that The result of Eq. (33) presents the Leading and Nextto-Leading orders of the semiclassical expansion; however, it is still the full nonperturbative answer in the sense of the coupling constant. In order to perform expansion in g we need to express φ c (x) and Φ 1 (x 1 ; x 2 ) through the non-interacting counterparts: where φ 0 (x) and G 0 R (x 1 , x 2 ) are solutions of the free differential equation (Eq. (11)) The iterative expansion of the Eq. (36) up to g 2 is the following: The higher variations Φ 3 and Φ 4 can be rewritten through φ c and Φ 1 with the help of the integration representations of the differential equations of Eq. (31). However, it is enough to expand the higher variation only up to g, because of the addition power of g in the second term of the Eq. (33). Moreover, the contribution of the Φ 4 vanishes, because the lowest term in this variation proportional to g 2 . The remaining h 2 term of the Eq. (33) is:

Let us draw the contributions to the full retarded Green's function
All lines and vertices have the same meaning as in Sect. 3. The only new element -the grey blob -denotes the free field φ 0 (x). Since only φ 0 (x) depends on the initial conditions in the above expansion, it is straightforward to perform averaging according to the rule: Since we consider the Gaussian form of the Wigner functional (to satisfy the demands of the perturbative approach), the Eq. (41) represents the basic element of the Wick's theorem -the contraction of two φ 0 (x). For example, the "cactus" diagram, that we mention earlier in Eq. (18), is recovered from the fourth term of the expansion (Eq. (5)), or the last line of Eq. (39).
i.c. = 3 All other contributions of the expansion (17) are recovered correspondingly.
One can see that the Leading Order semiclassical term (the CSA) reproduces all the contributions of the g 2 terms of the perturbative approach except the last one. However, this term is subleading for a highly occupied initial state as we discuss in the next section.

Shear viscosity and the CSA applicability
In this section we consider the implementation of the semiclassical formalism to the numerical simulation. The CSA is useful for numerical calculations since the path-integrals over initial conditions can be done with the help of the Monte-Carlo approach [26,29]. Another advantage of the CSA is the possibility to take into account the strongly correlated initial conditions (non-Gaussian ones). Let us consider the shear viscosity as an example. In order to evaluate it one can use the Kubo linear response theory [29,[33][34][35][36][37][38][39], where transport coefficients can be expressed through the retarded correlator R μν αβ of two components of the stress-energy tensor T μν as: The Kubo theory is valid if a system is in a (quasi)stationary state when hydrodynamical description [40] is applicable and R μν αβ (x; x ) depends only on x −x . In the rest frame, the shear viscosity can be expressed through the Fourier transform of the (12-12) correlation function: as: However, for the real physical systems, it is useful to consider the shear viscosity in the general frame where the energy flow can exist. One can define the flow velocity u μ for energy current as the only time-like eigenvector of the average stressenergy tensor T μν with eigenvalue equal to energy density ε. We normalize the flow velocity as u μ u μ = 1. In this case, one can rewrite the expression for shear viscosity in the covariant form through the retarded correlator of the traceless part of the stress tensor. The final expression is: where: The Eq. (44) describes the quasistationary state where the viscosity can be a slow function of the space-time point x, rather than the constant value of Eq. (43). It is likely to obtain such a physical system in the later stages of the nonequilibrium matter evolution. Moreover, the Eq. (44) allows evaluation of the shear viscosity in the system with the nonzero energy flow without transformation to the rest frame of the medium.
Note that the Eq. (44) can be considered as the shear viscosity only in the range of validity of the hydrodynamic description. For the correlator of the Eq. (44) it means that αβ μν R μν αβ (x + y; x) weakly changes as a function of x and fast decays as a function of y. In other words, all the microscopic dynamics should enter into the large scale behaviour only through the number of the transport coefficients.
For the system under consideration the stress energy tensor equals to: and only the first part Θ μν (x) = ∂ μ ϕ(x)∂ ν ϕ(x) contribute to the shear viscosity. It means that we need to evaluate retarded correlator which is proportional to Θ μν cl (x)Θ αβ q (x ) where the definition of the "classical" and the "quantum" components are the same as before: In terms of ϕ cl(q) the result is: where ∂ α = ∂ ∂ x α . Now we can apply the semiclassical approach and derive the leading order contribution for R(x; x ). We obtain the CSA answer: Here φ c (x) is the solution of the classical EoM and Φ 1 (x; x ) -the solution of the linear differential equation (32). The averaging over the initial conditions · · · i.c. is done with the Wigner functional Eq. (25) defined by a physical problem under consideration. The dynamics in this approximation is classical; the onlyh -corrections can be obtained from the Wigner functional. In principle, by taking variational derivatives of the Eq. (47), one can obtain the next-to-leading order semiclassical terms (h 2 ) or the quantum corrections to the CSA. Another way to obtain these quantum corrections is to use the technique described in Sect. 7 below. Let us explicitly consider the two similar diagrams Eq. (48) which are included to the Eq. (46). The first diagram contributes to the leading order shear viscosity, and it is taken into account by the CSA, whereas the second one proportional toh 2 and belongs to the NLO semiclassical term.
Here the squared vertices arise from the value we consider -the shear viscosity Eq. (46) as: One can observe that the difference in (48) comes only from the central loop. In the first case, it contains the product of two Keldysh Green functions ∼ G 0 K G 0 K , whereas the second one has ∼ G 0 R G 0 R insertion. Every Keldysh Green function has 2 f p + 1 multiplier in contrast to the retarded one (see Eqs. (16), (13)). If the initial state is highly occupied (i.e. the system is almost classical), then f p 1, and we can neglect the second contribution. This analysis can be extended to any diagram, contributing to the viscosity or any other observable. For each diagram ∼h 2n there is the diagram ∼h 0 which differs by 2n times substitutions of G 0 R by G 0 K . In other words, in this diagram n "white" vertices are changed to the "black" ones. The resulting diagram is greater due to 2 f p + 1 factors and is already included in the CSA. It explains why the CSA works well for the highly excited initial state and sums up all leading contributions. Hence, the results of works [27,31] are clarified.

7h 2 diagram technique
In previous sections, we consider two different approaches. The first one is the Keldysh-Schwinger diagram technique with two propagators G R/K (x, x ) and two vertices as building blocks. The other one is the semiclassical approach where the main objects are the classical solution φ c (x) and the classical response functions Φ n (x; x 1 , . . . , x n ) for a given initial condition. In order to obtain observables in the latter approach, one should perform the averaging over all possible initial conditions. In Sect. 5, we show that both approaches are equivalent for small g andh. Also, we demonstrate that the leading contributions of the semiclassical expansion (the CSA) represent a sum of an infinite number of diagrams in the perturbative approach. In this section, we develop a new approach that combines the advantages of both the CSA and the diagram technique and permits systematic improvement of the CSA and analysis of the higher-order quantum corrections. Let us start again from the general expression for the some observable (Eq. (6)) and shift the integration variable is the solution of the classical equation of motion with the corresponding boundary conditions. After thatφ cl obeys zero boundary conditions and all dependence of α(x) and Π(x) enters into the path integral only through φ c (x). This trick can be done for any quantity of interest. For example, the full retarded Green's function can be represented as: Now we perform the perturbative expansion of the above expression (without averaging over the initial conditions yet). One can check that due to the zero initial conditions forφ cl there is only one non-zero contraction: Here, by double angle brackets . . . we denote the functional integration without averaging over the initial conditions. This contraction is represented by the dashed line in diagrams: The "price" for the absence of the φ clφcl propagator is the presence of four different vertices in the theory: With the help of the above definitions, the perturbative expansion can be performed in a standard manner.
The advantage of thish 2 expansion (which is absent in the standard Keldysh-Schwinger technique) is that the number of diagrams for a givenh 2n order is finite. Moreover, the diagrams of theh 2n order have exactly n loops. Let us prove these facts. Due to causality Φ 1 (x; x ) is nonzero only for t > t . It means that any loop of the form y) is zero, and one can fix the direction of a time flow for each diagram. The proof is based on the observation that the "white" vertices Eq. (53) increase the number of lines in time, whereas the "black" vertices Eq. (52) decrease. However, for each fixedh 2n order, there can be only n "white" vertices, so we have only the finite number of possibilities to insert the "black" vertices. The statement about the number of the loops is based on the observation that it is impossible to construct a loop using only black vertices because there is only one outgoing leg in each vertex and all lines should be oriented along the time direction due to causality. Hence, every white vertex is associated with two loops in the diagram.
For example, let us consider the two-loop diagram: x x x 1 Up to a numerical prefactor it is equal to: By following the time direction of this diagram, one can observe that from t to t 1 there is only one line. The "white" vertex at t 1 creates three lines and two loops, which are closed by the "black" vertices at times t 2 and t 3 . The overall strategy for the usage of this new technique consists of three main steps: Let us note that the last two steps (solution of the equations of motion and the averaging) can be performed numerically. Usually, the Wigner functional decays fast for the large initial values of the field and the momenta. Since the Wigner functional is normalised: the "sign problem" does not occur for the quantum corrections, like in pure CSA simulations. Moreover, the precise form of the initial Wigner functional is not specified. Hence, with the help of this technique, one can work with correlated and highly nonequilibrium initial states. All complexity of the initial state appears only on the last step -the averaging over the initial conditions (which are included to the observables through the classical trajectory φ c (x)).
It is interesting to note that the separation of the path integral evaluation and the averaging over initial conditions in some sense are similar to the method used in work [41]. In this work, authors attribute one Lefschetz thimble to every initial condition and average with the Wigner function at the end of the calculation.
It is important to mention the physical interpretation of φ c and Φ 1 (x; x ). Since the initial conditions are random, then the value of φ c (x) at a given space-time point is also random in some sense. Therefore, one can think of φ c (x) as of some noisy background fluctuating in space and time. At the same time the equation for Φ 1 (x; x ) describes the retarded propagator for a scalar particle in that noisy background. It looks similar to the Langevene equation for a Brownian particle interacting with the environment (see, for example, [42] and references therein). However, here noise is generated not by external bath as in the case of Brownian particle but by the scalar field itself. From this point of view, vertices in our diagrammatic technique correspond to the creation or annihilation of additional scalar particles. Contribution from each event of a particle creation contains an additional smallh 2 factor. It means that in the semiclassical limit, only contribution from noisy background survives, and the CSA works well.

Conclusions
In this work, we compare two approaches to the descriptions of the nonequilibrium quantum scalar fields: -The standard Keldysh-Schwinger diagram technique, which requires the Gaussian initial conditions and the small coupling constant. -The semiclassical expansion, which works with the arbitrary coupling constant, but valid for highly excited (or highly occupied) initial states only.
We analyse these two expansions at the limit where the coupling constant g and the Plank constanth are simultaneously small. We prove the consistency of these approaches and explicitly demonstrate that already the first term of the semiclassical expansion (the Classical Statistical Approximation) includes almost all two loop-diagrams of the standard perturbative approach. We show that the only remaining g 2h2 diagram is small if the initial conditions are overoccupied i.e. the one-particle distribution function f p 1. In practice, this condition defines the applicability of the CSA.
As an example of the usefulness of the semiclassical approach, we evaluate the shear viscosity in a more general case of nonzero energy flow.
Also, we present a new diagram technique that combines both the advantages of the semiclassical and the Keldysh-Schwinger diagrammatic approaches. We believe that this technique allows to perform the resummation of the next-toleading order semiclassical contributions and to improve the CSA.
Acknowledgements This work was supported by the RFBR project 18-02-40131.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This work does not contain experimental data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .