Relativistic invariance of the vacuum

Relativistic invariance of the vacuum is (or follows from) one of the Wightman axioms, which is commonly believed to be true. Without these axioms, here we present a direct and general proof of continuous relativistic invariance of all real-time vacuum correlations of fields, not only scattering (forward in time), based on closed time path formalism. The only assumptions are basic principles of relativistic quantum field theories: the relativistic invariance of the Lagrangian, of the form including known interactions (electromagnetic, weak and strong), and standard rules of quantization. The proof is in principle perturbative leaving a possibility of spontaneous violation of invariance. Time symmetry is, however, manifestly violated.


Introduction
Relativistic invariance of the vacuum is a basic property of all quantum field theories based on relativistically invariant dynamics, directly following from one of Wightman axioms [1]. The (third) Wightman postulate assumes the existence of a single ground state with zero energy-momentum eigenvalue [2]. The postulate is generally accepted by high energy community [3] although the problem is more sophisticated than it looks and the validity of the axiom itself has been questioned [4,5]. However, to postulate zero energy is in conflict with renormalization which may add indefinite contribution. Moreover, this postulate should be redundant since the vacuum is already determined by the Lagrangian and the rules of quantization. Indeed, without the Wightman axiom, the invariance has been already proved for special cases of relativistic quantum field theories: single fields (electromagnetic, current, etc.), free theories, Feynman (forward) time-ordered (in-out) correlations of fields, scattering processes, second order and some other classes of correlations [6][7][8][9][10][11], but never in general (for arbitrarily timeordered correlations). A large part of high energy community wrongly claim it is already proved in general [12]. For a e-mail: abednorz@fuw.edu.pl instance, a naive proof by reconstruction of vacuum wavefunction out of Feynman correlations fails because such a reconstruction is incomplete and not unique. A general proof requires to show Lorentz invariance of all arbitrarily ordered correlations at real spacetime points. The invariance could be certainly proved only if the underlying dynamics is invariant and conserved charges are zero. It has been suggested that Nature may violate the invariance directly (by non-invariant dynamics) or spontaneously (by a Higgs-like mechanism) [13]. We do not discuss these two last possibilities here.
In this paper, we confirm the common intuition by a direct proof, not relying on Wightman axioms. Although formally general, it is better to understand it as perturbative. We show that zero-temperature vacuum correlations of fields at real spacetime points are invariant under continuous transformations of a reference frame. We shall use the framework of the closed time path formalism (CTP) [14][15][16] where correlations are defined on the complex path (going downwards with respect to imaginary part). It is defined in a particular reference frame, hence the invariance is not manifest. We only assume basic, accepted principles of relativistic quantum field theories: (i) relativistically invariant Lagrangian, including known interactions (electromagnetic, weak and strong, specified in detail later in the paper) and (ii) rules of quantization (standard construction of field expectation values, consistent with CTP). However, the discrete transformations like time reversal and parity have to be treated separately. If the underlying dynamics is invariant under charge conjugation or parity reflection then the proof remains valid for these transformations [1]. Unfortunately, even if the dynamics is symmetric with respect to time reversal, then CTP is not (except special cases, e.g. at space-like points or in scattering problems). In the same way, CTP is in general not invariant under joint charge-parity-time reversal in contrast to dynamics. It closely related to violation of time symmetry in quantum noninvasive measurements [17].
The paper is organized as follows. We first recall the construction of the complex-time path correlations, defining relevant correlations of fields. Next, we show that the correlations are independent of the particular shape of the path. Then we prove the Lorentz invariance of the zerotemperature no-charge vacuum for continuous transformations. Then we discuss possible pitfalls in other attempts of the proof based e.g. on momentum representation, which is dedicated to all who claim that the invariance is trivial or easy extension of existing proofs [12]. Finally, we will see that time symmetry is absent in CTP. Throughout the paper, we use the convention = c = 1 and β = 1/k B T (inverse temperature).

Closed time path formalism and definitions
The complex closed time path C is a specially chosen curve of time t in the complex plane as a function of a real parameter s. The curve is continuous with nonincreasing Im t (important for convergence of integrals, involving energies bounded from below and open from above) except a jump of iβ, which becomes imaginary infinity at zero temperature, see Fig. 1. The physical part of the contour is located on the real axis (both upper and lower). In particular, the curve in Fig. 1a The generalized contour in Fig. 1b contains multiple real parts. Fig. 1 The closed time path in the complex-time plane. All the path cannot go upwards in the imaginary direction. The finite temperature implies the jump iβ, which extends to infinity at zero temperature. Apart from the above conditions, the shape of the path is arbitrary. We stress that the path is ordered by the real parameter s. All the important times have been marked. The shift between real parts is infinitesimal only for better visibility but in fact can be zero. (a) The simple in-out-in (shortly in-in) contour with a single pair of real parts (between 0 and t 0 ). (b) The generalized contour for arbitrarily ordered correlations To facilitate an easier notation throughout the paper we define Heaviside and Dirac functions for complex-time arguments t, where, as usual, θ(s − s ) = 1 for s > s and 0 otherwise while dsδ(s − s ) = 1 and δ(s − s ) = 0 for s = s . Therefore θ defines a special order of times along the path, and it may be converse to the usual order. This is important since only s is the real, ordering parameter, while t is complex and thus has no natural ordering. Note that in all contour time-integrals we essentially integrate over s. Let us take the generating functional for quantum correlations (often referred to as Green functions) with respect to the action, making use of the Lagrangian approach, namely [11, 15, 18- where with auxiliary commuting (c-) field χ (χχ = χ χ ) and anticommuting Grassmann field η (ηη = −η η). Here t k can be located everywhere on the path (upper real, lower real, left complex part). In the above, A and φ are generic bosonic and fermionic (Grassmann) fields, respectively, and the Lagrangian is a time-local function of fields and their time derivatives, namely L = L[A,Ȧ, φ,φ]. The integral C dt is over the just defined time path, and should be read as an integral over real s, ds(dt/ds). The boundary conditions for fields are A(t i ) = A(t f ) and φ(t i ) = −φ(t f ), for bosonic and fermionic field, respectively (Kubo-Martin-Schwinger condition [32][33][34]). Standard (in-out) correlations are obtained for the auxiliary fields nonzero only on the upper real part of the contour. General correlations involve both parts of the contour (in-out-in or in-in) or even more (in-in-in, etc.). The path-integral formulation of CTP is equivalent to Hamiltonian formulation, as we show in Appendix A. The standard issues such as Wick theorem, conservation laws and renormalization are resolved in CTP analogously to inout correlations, see Appendix B.

Contour shape independence
The fact that we can freely wiggle the contour is reasonably believed for CTP [24,[35][36][37], but a direct proof is very instructive [38,39]. Let us prove that (4) is independent of a particular shape of CTP, because the proof of Lorentz invariance will be analogous. It follows from the invariance of DA and Dφ under infinitesimal transformation δt (s) parametrized by real s (the variation may change endpoints but the jump must be kept constant) which gives L(t + δt) = exp(δt∂ t )L(t). This means that the effect of the contour changes can be absorbed into a shift under time integral and a such it disappears (analogously to an integral over a total derivative). We stress that this proof is valid only if do not change the shape of the real part of the contour, except shifting t 0 if it lies beyond all relevant times. It is essential that the fields are smooth. An alternative, more explicit proof relies on perturbative approach. One starts with Gaussian, diagonalized Lagrangian In the above, we have denoted complex ψ instead of real φ and distinguished the linear time derivative bosonic field A from the quadratic B. The two-point Green function for B reads where For A we get (A.7) and for ψ (A.8).
If we add interaction term, L = L 0 + L I , and L I depends on fields appearing in L 0 , the generating function (4)-by virtue of the Wick theorem (B.1)-will be a sum of integrals represented by Feynman diagrams, with lines (propagators) corresponding to free Green functions (A.7), (A.8) and (7) and vertices corresponding to L I . Suppose that t (s) is varied by infinitesimal δt. Let us consider the variation of Green functions. Note that δ X(t)X * t 0 = δt∂ t + δt ∂ t X(t)X * t 0 (8) where X = A, B, ψ and subscript 0 denotes averaging (3) with respect to L 0 . The free Green function depends only on endpoint times t, t , not the shape of CTP between them. Now, each physical diagram corresponds to average of the type where X = ψ, A, B and t k lie in the real (right hand) part of the contour for k > N. So we get where we used Wick-decoupling into products of free Green functions, Leibniz rule δ(AB) = BδA + AδB, δ dt = dδt = ds∂ s δt, (8), inverse Leibniz rule applied to ∂ s and back Wick-coupling. The final result of the shape variation integral is zero, which completes the proof.

Lorentz invariance
It is natural to expect that Green functions are invariant under relativistic Lorentz transformations at zero temperature if the Lagrangian is Lorentz-invariant. However, starting from CTP makes it not obvious because the contour prefers some time direction. Intuitively, we expect that this should not bring about any problem but the warning light comes already from the non-invariance due to finite temperatures, which enters only as a jump in the contour. Below, we show that the Green functions at real times in zerotemperature vacuum are indeed invariant under continuous Lorentz transformations in Lorentz-invariant field theories but not necessarily under time reversal.

Relativistic notation in CTP
Let us briefly recall relativistic notation. The components of position four-vector are denoted x μ , μ = 0, 1, 2, 3, with time x 0 = t and three spatial coordinates x k , k = 1, 2, 3. Fields, functions of position will be denoted as . The derivatives are denoted by ∂ μ = ∂/∂x μ . From now on, Greek indices are reserved for fourvectors and higher Lorentz-like tensors and structures. It is extremely important for the purpose of this paper that the time x 0 (s) is here a complex function of the real parameter s. In particular the time derivative ∂ 0 is translated into the real derivative ∂ s by Similarly, we introduce Dirac four-delta, including complex time, and four-integral C d 4 x with d 4 x = dx 0 dx 1 dx 2 dx 3 and dx 0 = (dt/ds) ds. For a compact notation of relative four position including the complex-time component of the fourvector on CTP we will denote |x − where θ(x 0 − x 0 ) = θ(s − s ) applies to complex x 0 and x 0 along CTP. The rest of conventions, including Lorentz generators J , spinors and Lagrangian density L (becomes Lagrangian L when space-integrated) are consistent with textbooks [7,8,10,40], which we recall in Appendix C for completeness.
The essential postulate of QFT is that the Lagrangian density is a Lorentz-invariant. It means that the generating functional (4) can be written now as where all fields are functions of spacetime, e.g. η k (x) with complex time x 0 = t, L transforms as a scalar field and normalization imposes constraint S[0, 0, 0] = 0. Here A k denotes Lorentz scalars and B μ k -four-vectors with the dot · denoting Minkowski scalar product. In our convention, assuming that L is a scalar we also assumed that average charge densities are zero. The first Minkowski coordinate (time) has the jump x 0 i − x 0 f = iβ. We stress that the charge conservation, Ward-Takahashi identity, holds in CTP if the Lagrangian does not mix different classes of fermions [10] The great advantage of Lagrangian and path-integral formalism is that it seems to be perfectly Lorentz-invariant at first sight. However, the possible caveats are hidden in the shape of CTP, which is bent into imaginary time but not spatial coordinates. The common sense intuition tells us that the Lorentz invariance must be broken by finite temperatures, which essentially dictate the contour jump. Moreover, even at zero temperature, the contour retains preferred time direction so Lorentz invariance is not manifest and self-evident.

Formal proof of invariance
We shall now present a brief formal proof of Lorentz invariance of vacuum generating function (4), with the auxiliary field nonzero only for real times (right hand part of the contour) at zero temperature k B T = 0 or β = +∞. Note that a zero temperature makes the jump on CTP between t i and t f infinite. Due to the freedom of choice of the shape, we assume the contour in Fig. 1 where we have used the fact that both the measure d 4 x and the integration boundary is invariant under Lorentz trans- All the exponent in the generating functional (14) is a Lorentz scalar so the above reasoning remains valid. We write Zero-temperature CTP. The boundary times t i and t f are sent to imaginary infinity, i∞ and −i∞, respectively. Similarly to Fig. 1 we can take the simple contour (a) and the generalized one (b) which proves the invariance. However, we had to assume that Λ acts on functions of complex arguments, not just real domain, as time is complex. To be correct we should work only in the real domain (s, x 1 , x 2 , x 3 ) and define Λ by Taylor expansion of generator exponential This is equivalent to real Lorentz transformation for real times, which is the case we need (Green functions at real times). The above proof, although formally correct, is hardly acceptable because one uses ill-defined mathematical operations. The vague issues are: -Lorentz transformations are analytically continued to complex spacetime. -We ignored possible problems for x 0 → ±i∞.
-Convergence and renormalization is not discussed.
To give more insight to the proof, dispel possible objections and make it robust, we will now reformulate it in terms of perturbative diagrams showing that each diagram separately is Lorentz-invariant.
It is enough to prove that J G = 0 for all Lorentz generators J and all Green functions G-defined by (3) at real times, because every finite continuous transformation can be written as Λ = exp wJ and (d/dw)ΛG = J ΛG = J G w where ΛG = G w is always some Green function at real times.

Free theories
We shall first prove Lorentz invariance of vacuum zerotemperature Green functions in free theories, where they can be found exactly. The Lagrangian density for free scalar theory reads where m is the mass of the field. The Green function (propagator) satisfies Klein-Gordon equation where ∂ acts on x (we shall write ∂ if acting on x ). The solution reads where p μ is the four momentum with p 0 = p 2 + m 2 . It is important that (21) is valid for time components, x 0 and x 0 along the CTP contour as defined in Fig. 1 and (13). For the Dirac field, we have with the summation over c and m ab = mδ ab . Now, we get two cases. If X = A (bosonic field, unphysical, usually ap- where For vector fields B μ , the internal Lorentz structure is simply added by metric tensor Possible problems occur for expressions of the type ∂ μ ∂ ν G 0 because the derivative produces unphysical contact terms ∼ δ(x 0 − x 0 ) [40]. Fortunately, they appear always in pairs with a renormalization field, which cancels the singularity. Nonzero temperature β < ∞ will break Lorentz invariance, which means that there is a four-vector related to temperature. There is no unique choice. One can define β μ [41,42] with β 0 = β and β k = 0, k = 1, 2, 3 in a particular reference frame. However, good alternatives are k B T μ = β μ /β · β or u μ = β μ / √ β · β, see also a general discussion [43]. However, one has to remember that the zero temperature limit means infinitely time-like β μ . The zero four-vector β limit, when Lorentz invariance might be obvious, corresponds unfortunately to infinite temperature so this is not what we are looking for.
In the zero-temperature limit β → +∞, (21) and (25) take the form Applying Lorentz generator acting on x and x by J 0x and J 0x , respectively, one can check explicitly (a lengthy textbook exercise) that J G = 0 (so ΛG = G) when applied to (26), so the free theory is indeed Lorentz-invariant (there are no additional problems with transformation of internal vector or spinor structure). Remember that the time derivative in J is always taken along the contour, namely with complex-time derivative ∂ 0 (11).

Interaction
The rest of the proof is similar to the shape-independence proof. Again, it is important that all time components below, x 0 k , run along the CTP contour as defined in Fig. 1 (in particular Fig. 2 at zero temperature) and (13). The Green functions k X k (x k ) (28) are, by virtue of the Wick theorem, represented by sums of diagrams with external vertices k (not integrated), internal vertices (integrated) corresponding to the interaction part of the Lagrangian density L I and connecting lines corresponding to free Green functions (21) and (25). We want to show that J vanishes when applied to the whole diagram. Suppose that the diagram has fixed number N of internal vertices 1..N and K external ones N + 1..N + K. then we have to show that where X k represent fields whose expectation value define the Green function, J m acts on all fields meeting at the vertex m. The average is taken with respect to free Lagrangian density L 0 . We decouple (29) Since L I is a Lorentz scalar we can change n J n → n J n 0 . Finally, we perform integration over each vertex n ≤ N (x = x n ), which yields only boundary values due to dt x k ∂ 0 = ds x k ∂ s , see Eq. (C.7). However, free propagators vanish exponentially at large distances and, at zero temperature, with x 0 → ±i∞. Hence (31) vanishes along with the left hand side of (29), which completes the proof. This is our main result. The proof has been obtained in the time domain, without stretching the real time to infinities. Note the analogy to contour shape independence. The heart of the proof is the fact that analytic continuation of time does not hurt the generating function. It is quite intuitive but never before shown explicitly.

Incomplete proofs
It is tempting to ask: Why not to prove Lorentz invariance using simply energy-momentum representation [44][45][46][47][48]? This approach seems attractive and practical but it has to be supplemented by arguments following from spacetime CTP representation. The main problems arise for offshell Green functions (not describing scattering processes) and field mass/strength renormalization. They cannot be resolved without returning to spacetime representation.
The energy-momentum representation can be used only on the real part of the contour stretched to the whole real axis (t 0 → +∞ and t (s ± ) → −∞), taking separately its upper and lower branch, as Fourier transform is only defined in real axis. The Green function gets additional indices to keep track of the branch, where, in the space (+, −), where δ ± (p 2 − m 2 ) = θ(±p 0 )δ(p 2 − m 2 ), the upper index +/− denotes bosonic/fermionic field with n ± (q) = (e βq ∓ 1) −1 . We denoted ε → 0 + defined for momenta to distinguish it from defined in spacetime. In the zero-temperature limit n → 0 and the free Green functions are manifestly Lorentz-invariant. The problems arise for interactions. Now, the interacting Green function takes on renormalization corrections, where Z(p) is regular at p 2 = M 2 . The renormalization is obtained by summation of Dyson series of self-energy contributions. However, such step is ill-defined for p 2 = M 2 because ε dominates in this regime. This can be cured by a slight tilt of the real-time contour as in Fig. 3 [10]. In fact, the new contour contains the essential feature of Fig. 2the endpoints in ±i∞. However, for G ±∓ the trick with the contour does not help. We are forced to simply accept because this part is only defined for p 2 = M 2 . One can circumvent the problem by ad hoc defining extrapolation of G ±± to G ±∓ and the p 2 = M 2 case. Another argument is the fluctuation dissipation theorem, which imposes relations between elements of G [49,50], or mass derivative [36]. However, this is only the case of elementary particles like electron and photon. Composed bound states like nucleons, atoms, and molecules get their renormalization corrections, too, but they cannot be simply reduced to a single dressed propagator line. In the time domain, bound states do not lead to any problems because the poles/deltas in propagators appear only in the infinite time limit and we do not have to sum infinite series of diagrams. Further problems may arise from pinch singularities [51], which may occur for products of propagators [(p 2 − m 2 + iε)(p 2 − m 2 − iε)] −1 . Fortunately, a careful analysis makes them always cancel [46,47]. Having the invariance of in-out correlations, it would be easy to prove Lorentz invariance if all the Green functions factorized at zero temperature [41,48], namely where k± denotes fields defined on the upper(+)/ lower(−) branch by the Fourier transform X ± (p) = Fig. 3 The tilted contour which circumvents the problem of pole renormalization p 2 ∼ m 2 → M 2 d 4 x ± e ip μ x μ± X ± (x ± ). An attempt to prove the factorization requires generalization of CTP to the contours like that in Fig. 4 with varying σ [41,44,45,48]. If σ = β/2 then the factorization occurs indeed for β → ∞. However, not every Green function is σ -independent [52,53]. From translational invariance the nonvanishing Green functions satisfy p = 0. The new vertical part of the contour gives additional factors e ±σp 0 to each propagator from lower to upper branch (+) or vice versa (−). These factors cancel only if which means that the frequencies/energies of upper and lower branch separately add up to zero as seen in Fig. 5. This situation is common to scattering processes but not in general. It holds e.g. for in-out correlations when all times are on the upper real part due to time shift invariance but not necessarily for in-in correlations when times are located on both real parts (upper and lower).
Having presented above the complications in energymomentum space, we conclude that the proof is fully correct in the spacetime domain, while in the momentum/Fourier space it would require assuming at best several additional technical rules and may not be general (for very complicated diagrams). Still, most of practical calculations are easier performed in the latter case, with spacetime arguments used only when the momentum rules are ambiguous.
Another attempt of proof requires reconstruction of vacuum wavefunction out of in-out Green functions. This is possible in simple quantum models, including free theories, finite and harmonic systems but not in fully interacting relativistic theory because renormalizing fields, especially for fermions, make the information in in-out functions incomplete and insufficient to uniquely reconstruct vacuum wavefunction. On the other hand, dimensional regularization is useless when proving Lorentz invariance because Lorentz transformations are defined only for integer dimensions. Lastly, an often heard argument is analyticity of Green functions [12], but how to define it on CTP? The condition of nonincreasing imaginary part makes it impossible to define analytic conditions (Cauchy-Riemann) because for free real times the middle one is pinched-the imaginary part of the derivative is not allowed. One cannot allow the contour to go upwards in imaginary direction because the energy spectrum is bounded from below and open form above, making integrals nonrenormalizably divergent. On the other hand analyticity with respect to real s is questionable because of turning points t 0 (s 0 ) and crossings (t n t m ). Such a proof will be restricted to a particular branch of CTP, e.g. upper/lower real or vertical imaginary part, but not general. One cannot also use Gell-Mann and Low theorem [54], which could simply extend the proof from free theories to interacting ones by adiabatic switching, because it works itself in a preferred frame. Any attempt to prove Lorentz invariance in general will ultimately fail or arrive at the proof here presented.

Absence of time-reversal symmetry
The momentum representation hides the fact that the Lorentz invariance is not necessarily valid for time reversal, which is only true for second order correlations or if the factorization (36) holds for higher order correlations. Note that the reversal p 0 → −p 0 break the time symmetry because of θ(±p 0 ) factors in G ±∓ . Trying to reverse branches + ↔ − does not help either, because then G ±± gets conjugated, so there is no way to restore original Green function. If we both reverse p 0 and make conjugation then we will return to the original Green function but the time will be also twice reversed which is again not what we wanted to have. The time (generally charge-parity-time) symmetry is valid only for scattering processes or space-like separated points. This lack of symmetry is fundamental and especially surprising in context of noninvasive measurements which are naturally defined on CTP [17].

Conclusions
We have shown that the closed time path formalism is consistent with relativistic invariance at zero temperature. Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

Appendix A: Hamiltonian formulation of CTP
The traditional starting point of every quantum-mechanical problem is the Hermitian Hamiltonian operatorĤ in a given Hilbert space. The physical quantities are described by Hermitian operatorsX k corresponding to their classical counterparts X k (numbers). Depending on the measurement scheme, there exists a correspondence between operators and measurement outcomes,X k → X k . Evolution and state is described by HamiltonianĤ (t) →Ĥ +Ĥ 1 (t), wherê H 1 (t) is the time dependent part, due to changing external forces, which are assumed to be absent for Re t < 0. We shall assume that t lies on CTP as in Fig. 1. Every conserved quantity can be incorporated into the definition of the Hamiltonian,Ĥ →Ĥ − μ QQ . For the main purpose of this paper, H 1 andQ can be omitted. However, the whole CTP formalism works perfectly also for nonzeroĤ 1 , which is necessary in problems involving external (changing) forces and charges [55]. The generating functional can be written in the form where the subscript C denotes integration over the CTP as defined in Sect. 2 and the time order T C denotes ordering of operators in the Taylor expansion according to CTP, namely T CŶ1 (t 1 ) · · ·Ŷ n (t n ) = sgn(σ f )Ŷ σ (1) (t σ (1) ) · · ·Ŷ σ (n) (t σ (n) ), where σ is a permutation such that s σ (1) ≥ · · · ≥ s σ (n) (in the case ofX k we take the time of χ ). Here sgn(σ f ) denotes the sign of the sub-permutation fermionic operators. IfŶ k is a fermionic operator for k ∈ {1..n f }, then σ f is defined as the map between the ordered sequences (σ (k)) n f k=1 and (σ (k)), ordered according to k and σ (k), respectively.
The correlations (Green functions) are given by derivatives of the generating functional, Performing a Taylor expansion of exponentials in (A.1) one can explicitly check that the result does not depend on a particular shape of CTP [38,39], as long as (1) is satisfied and t 0 is greater than all interesting times. In such a check it is important to note thatĤ is time-independent along the wiggly part of the contour, so the ordering does not matter there. On the other hand (A.1) is useful also for deriving the thermodynamic functions and transport coefficients, e.g. the Kubo formula [32][33][34].
Usually one separates free harmonic HamiltonianĤ 0 writingĤ (t) =Ĥ 0 +Ĥ I (t). The remaining calculations are usually performed perturbatively, with the help of the free two-point Green's functions and Wick's theorem [56,57], which allows a decoupling of the free many-point Green functions into products of two-point functions, see below.

A.1 Free Green functions
The most convenient free Hamiltonian is the quadratic form of bosonic and fermionic operatorŝ wherex andφ are Hermitian operators with bosonic and fermionic commutation relations, respectively, [x k ,x l ] = ig kl1 , {φ k , φ l } = h kl1 , with real g and h.
One can always diagonalizeĤ 0 so that where H vac is the vacuum energy (can be ignored),Â and ψ are linear combinations ofx andφ, respectively, with the property It is especially simple and instructive to find Green functions for the above Hamiltonian, according to (A.3) and (A.1), extended to the whole CTP. In the case of bosonic operators, For fermionic operators, it is important that φ (ψ , ψ * ) is never a c-number but a Grassmann number with the property φ a φ b = −φ b φ a , which follows from anticommutation relations. The anticommutation makes it necessary to use different boundary conditions, due to Kubo, Martin and Schwinger [32][33][34], imposing the sign reversal on the jump t f → t i , so the Green functions read ψ(t) 0 = ψ * (t) 0 = 0 and In both cases, the Green functions satisfy the equation

A.2 Wick theorem
The many-point Green functions are obtained from Wick theorem [56,57]. For products of odd number of operators, Green function vanishes while for the even number 2n, 2n k=1 x k (t k ) (A.10) The Wick theorem states essentially that every manypoint Green function for a quadratic Hamiltonian splits into products of two-point functions. This fact is analogous to the property of Fourier transforms of Gaussian functions, which are again Gaussian. Wick theorem is useful for interactions-one proceeds perturbatively, expanding (A.1) in interaction strength.

A.3 Hamiltonian-Lagrangian equivalence for bosonic fields
Let us begin with the Lagrangian for generalized coordinates q andq = dq/dt. The classical Hamiltonian is obtained by Legendre transformation, defining momenta, so that the Hamiltonian reads .
We assume quantum evolution of the wavefunction Ψ (q) given by the path integral where Dq includes all necessary normalization factors. To find the operator form of the Hamiltonian, we expand the above equation for small times, The ignored terms are of the order ∼ t 2 . We now use the fact that Let us briefly comment the last term in bracket. It has been obtained by assuming θ(t)δ(t) → δ(t)/2 or symmetrizing the last integral in (A.15) [58]. However, in the fundamental theories, like quantum electrodynamics, weak, strong interactions and generally Standard Model this term is absent. It follows from the fact that the even the most dangerous terms, due to non-Abelian gauge fields, have ∂f k /∂q k = 0. Nevertheless, this indicates potential problems for completely general Lagrangians, especially for gravity. Finally, we get which proves the equivalence between bosonic path integrals and Hamiltonian picture, with the commutation rule [q k ,p l ] = iδ kl .

A.4 Hamiltonian-Lagrangian equivalence for fermionic fields
Fermionic path integrals are more complicated due to anticommutation rules. We have to define abstract anticommuting Grassmann numbers, φ k , with the property φ k φ l = −φ l φ k and the integrals [10], One can check that Here φ = φ r + iφ i and φ * = φ r − iφ i are independent fields and dφ * dφ = i dφ i dφ r /2.
Let us now take an abstract Lagrangian According to the classical rule we can construct an abstract Hamiltonian In the continuous case for CTP, which is the equivalent to (A.8). The equivalence between bosonic-fermionic path integrals and Hamiltonian operator can be now easily extended to general family of Lagrangians of the type where f k (q, ψ) and g(q, ψ) are real-valued and even in ψ, and Hamiltonian with commutation relations [p k ,q l ] = i and {ψ k ,ψ * l } = δ kl . We emphasize, however, that the equivalence does not generalize to completely arbitrary Lagrangians, containing higher powers of time derivatives, e.g.q 4 orψ 2 .

B.1 Wick theorem for path integrals
The Wick theorem is now a straightforward consequence of Gaussian integrals DqDφe − kl (q k g kl q l −φ k b kl φ l )/2+ k (χ k q k +η k φ k ) ∝ e kl (χ k g −1 kl χ l −η k b −1 kl η l )/2 , (B.1) where g and b are real symmetric and antisymmetric matrices, respectively, while g −1 and b −1 are their inverses.

B.2 Conservation laws
As noticed already by Kadanoff and Baym [15], the CTP technique maintains classical conservation laws. Suppose that the fermionic terms in Lagrangian (A.25) appear in separate sets, where L I is built from bosonic field and fermionic quadratic forms Sets N form a family of disjoint sets and L I may contains arbitrary products and sums of (B.3) for different N s together with bosonic fields. It is, however, forbidden to include terms like ψ * k ψ k where k and k belong to different sets N and N , respectively. Each set N defines one conserved quantity, Q N (t) = n∈N ψ * n (t)ψ n (t). The conservation law states that 3) so that they are invariant under ψ → ψ . We could certainly repeat the proof perturbatively, as we did to prove shape independence in Sect. 3 (some issues are resolved by renormalization [10]).

B.3 Hubbard-Stratonovich transformation and renormalization
Here we discuss the two techniques often used in context of path integrals. The Hubbard-Stratonovich transformation [59,60] allows to reduce order of terms in path integral by an auxiliary bosonic field X(t), Note that √ a is imaginary if a < 0 and in this case X does not correspond to any physical field although path integrals can be still performed. The renormalization is used to kill divergences appearing in quantum field theory [7,8]. It requires introducing unphysical fields, too. The most common form of renormalizing terms in Lagrangian is whereX m is the unphysical field as it appears with i. Both techniques are fully consistent with CTP formalism. For the purpose of this paper-the proof of Lorentz invariance, it is more convenient to use Pauli-Villars renormalization scheme, with up to several renormalizing fields with large masses (e.g. two bosonic and two fermions in the case of vacuum polarization) [61], because dimensional regularization makes Lorentz transformations undefined. For details of these techniques we refer readers to textbooks [7,8,10,40].

Appendix C: Lorentz transformations
For general multi-indexed relativistic structures we adopt Einstein convention A · B = A μ B μ = μ=0,1,2,3 A μ B μ and A 2 = A · A. We denote Minkowski metric tensor g μν = g μν equal +1 for μ = ν = 0, −1 for μ = ν = 1, 2, 3 and 0 otherwise. Next, g is used to lower/raise an index, A μ = g μν A ν , A μ = g μν A ν . In this notation g μ ν = δ μ ν is equal +1 if μ = ν and 0 otherwise and A 0 = A 0 while A 1,2,3 = −A 1,2,3 . The four-vector measure is denoted d 4 x = dx 0 dx 1 dx 2 dx 3 . We extract the usual three-vectors x = (x 1 , x 2 , x 3 ), its square x 2 = (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 3 and measure d 3 x = dx 1 dx 2 dx 3 . The natural Lorentz transformation is defined on four-vectors, The transformation must satisfy conservation of metric tensor In general, the continuous Lorentz transformation can be written in the form where J is the generator of the Lorentz group, satisfying J μν , J αβ = i g να J μβ − g μα J νβ − g νβ J μα + g μβ J να (C.4) which yields J μν = −J νμ . The generator for (C.1) reads J μν α β = i g μα g ν β − g να g μ α . (C.5) For the scalar field with the corresponding generator If the field has vector structure, e.g. A μ (x) then the generator reads J + J 0 , acting both on argument x and four-vector (indices μ). Apart from continuous transformations there exist two special discrete transformations, time reversal x 0 → −x 0 and parity (mirror) inversion x k → −x k , k = 1, 2, 3. One can also make charge conjugation e → −e (the interaction strength). The proof of Lorentz invariance of zero temperature vacuum will be valid only for continuous transformations. If Lagrangian is invariant with respect to parity inversion or charge conjugation then we can include it, too (which is e.g. not the case for weak interactions).
Instead of the Lagrangian, in relativistic spacetime we will rather use its density L with the property  and the list may be incomplete. When considering weak interactions, one should include also pseudoscalarsψ k γ 5 ψ l , A μψk γ 5 γ μ ψ l where γ 5 = iγ 0 γ 1 γ 2 γ 3 , which transforms as scalar under continuous transformations but changes sign under parity inversion.