Thermal out-of-time-order correlators, KMS relations, and spectral functions

We describe general features of thermal correlation functions in quantum systems, with specific focus on the fluctuation-dissipation type relations implied by the KMS condition. These end up relating correlation functions with different time ordering and thus should naturally be viewed in the larger context of out-of-time-ordered (OTO) observables. In particular, eschewing the standard formulation of KMS relations where thermal periodicity is combined with time-reversal to stay within the purview of Schwinger-Keldysh functional integrals, we show that there is a natural way to phrase them directly in terms of OTO correlators. We use these observations to construct a natural causal basis for thermal n-point functions in terms of fully nested commutators. We provide several general results which can be inferred from cyclic orbits of permutations, and exemplify the abstract results using a quantum oscillator as an explicit example.


Introduction
Understanding thermal equilibrium in quantum systems is a necessary precursor to developing a comprehensive physical picture of out-of-equilibrium dynamics. Fortunately, the simplicity of the thermal state, wherein the Gibbs density matrix is a simple functional of the quantum Hamiltonian,ρ T = e −β H , allows one significant insight. By analytic continuation to Euclidean time one can study an equivalent classical statistical mechanical problem; a compact thermal circle allows decomposition into discrete Matsubara modes and thus one can analyze various equilibrium properties. However, it remains of interest to ask how deviations from equilibrium can be physically quantified. These are captured by response functions, which by virtue of causal ordering of events, necessarily involve real-time ordering. A remarkable fact about thermal states is that the real-time response functions are related to fluctuations about equilibrium thanks to the fluctuation-dissipation relations, which in turn follow from the Kubo-Martin-Schwinger (KMS) relations [1,2]. These conditions have played a key role in not only in thermal field theory (cf., [3] for a review), but have also been useful in the context of axiomatic formulations of QFT [4].
Traditionally, the fluctuation-dissipation relations (abbreviated FDT) are characterized by relating the retarded Green's function in thermal equilibrium (causal response), to the symmetrized two-point function (the fluctuation). This well-celebrated generalization of the Einstein relation clearly admits generalization to higher point functions. Over the years many authors have attempted to construct a full set of relations, translating the information into a series of statements about higher point spectral functions in the field theory. Some preliminary attempts are for instance summarized in the excellent review [3]; further attempts to understand these relations within the real-time Schwinger-Keldysh formalism can be found in [5][6][7][8][9]. In addition, several groups have tried to ascertain the relation between the real-time and imaginary time formalisms, cf., [10][11][12][13][14][15][16][17]. These works demonstrate that an effective implementation of KMS relations can be very useful in simplifying the finite temperature Schwinger-Keldysh formalism as applied to φ 4 theories as well as gauge theories.
Much of the aforementioned analysis was inspired by the need to better understand the connection between the imaginary time formalism (or statistical field theory) and the real time Schwinger-Keldysh formalism for finite temperature QFT. Another motivation was the need to compute thermal observables in hot QCD plasma. More recently, there is growing interest in quantum systems driven out of equilibrium such as e.g., quantum quenches, owing to our improved ability to experimentally simulate time-dependent many-body quantum Hamiltonians.
The current work is aimed at synthesizing these developments by attempting to form a coherent picture that transcends the limitations of the Schwinger-Keldysh construction. To appreciate this perspective, let us first note that the general FDT arises from the fact that thermal correlators are trace class observables in the Gibbs density matrix. As such the cyclicity inherent in such observables descends directly onto the correlation functions. This is very clear if we think of the cyclic structure made explicit in the imaginary time formalism. Whilst simple, this cyclicity does not play well with causal ordering -cyclic transpositions of ordered operators are no longer ordered in a similar manner. In the Schwinger-Keldysh formalism we are able to compute correlators which have a restricted form of out-of-time ordering, but not the most general OTO correlator. Despite being reasonably intuitive, it was only recently argued in [18] (see also preliminary observations in [19]), that of the set of real-time Wightman n-point functions, characterized by n! orderings, only a small subset are computed by the real-time Schwinger-Keldysh contour. Most of the Wightman functions, are out-of-time-order, and are computed by more general functional integral contours, which involve many instances of forward/backward evolution. Such out-of-time order (OTO) functional integral contours, and correlation functions computed therefrom, have been the focus of recent interest owing to the intricate connection between chaos, ergodicity, thermalization, and black hole physics [20]. 1 Traditionally, the approach in the analysis of fluctuation-dissipation theorems has been to eschew the occurrence of OTO-correlators by suitably combining the KMS transformation with other discrete symmetries such as time-reversal or CPT, cf., [24]. This is not necessary, and indeed we will show that one can formulate the general set of KMS relations quite simply once we enlarge our considerations to include OTO correlation functions. One of the aims of the current discussion is to give a unified and general picture of such thermal relations from a real-time perspective and elucidate the simplicity gained by moving to a framework involving OTO correlators. For a special class of OTO 2k-point functions, fluctuation-dissipation relations for suitably regulated "bipartite" correlators have been investigated in [25]. 2 Summary of results: We can distill the essential features of thermal correlation functions into the following set of statements: • The KMS relations and general fluctuation-dissipation theorems for higher-point functions can be formulated for any time-ordering of the operators. One does not need to restrict the relations to only involve correlation functions obtained from the Schwinger-Keldysh path integral contour. More specifically, the n! Wightman functions can be 1 Various authors have also explored the OTO correlators from a quantum information perspective, cf., [21,22] for connections to quantum channels, and [23] for connections to weak measurements and quasiprobabilities. 2 It was also argued in [26] that the 4-point chaos correlator could be given an interpretation of a generalized FD relation along the lines of the Jarzynski relation [27,28]. We however believe that the relation derived there is better viewed as a suitable writing of the generating function after fusing two operators into a single composite, which is somewhat different from the more standard notion of fluctuation-dissipation relations.
partitioned into (n−1)! equivalence classes after taking into account the KMS relations, which act by cyclic permutations and imaginary time shifts of the operator insertions.
• Once one allows for this general perspective, one does not need to invoke any form of Z 2 involution originating from time-reversal or CPT. The rationale for doing so in more traditional presentations of KMS relations originates from the desire to relate correlators computed within the Schwinger-Keldysh contour (which in our classification would be a 1-OTO contour). There is no need for such a restriction when using the k-OTO path integral contours and the general perspective espoused herein.
• The KMS relations end up relating proper q-OTO Wightman functions with proper (q + 1)-OTO correlators. The proper OTO number refers to the minimal number of forward-backward evolutions necessary to account for the time-ordering in the correlation function. In particular, this set of relations has the effect of reducing the switchbacks in the path integral contours necessary for computing thermal Wightman functions.
• In particular, while it is well known that all time orderings in two-point functions are captured by the Schwinger-Keldysh 1-OTO contour, for three-and four-point functions we need to employ 2-OTO contours for generic initial states. However, in the thermal state the 2-OTO three-point functions can be related to 1-OTO 3-point functions via the KMS relations, thereby reducing the required proper OTO number. The first time we encounter a genuinely 2-OTO thermal correlation is for four-point functions, as for example illustrated by the now familiar chaos correlator [20].
• More generally, an n-point thermal Wightman function can be computed using at most an n 2 -OTO contour, whereas in the absence of the KMS relations one would have to go up to n+1 2 -OTO contours. Consequently, one may be tempted to speculate, as we do in §6, that even-point functions where we first encounter new OTO contour order (e.g., 2-OTO for four-point, 3-OTO for six-point etc.) may be the natural place to look for detailed features of how systems thermalize.
While the general set of statements above holds for any thermal Wightman functions, it is useful to express the relations in terms of nested commutators and anti-commutators. As is well known from the Keldysh construction, fully nested commutators capture causal response, and fully nested anti-commutators encode fluctuations. One expects based on the two-point FDT that there would be general relations between such objects and this is indeed borne out. We will in particular argue that n-point function generalized FDTs (which constrain the physics of thermal OTO correlators) can be described as follows: • The set of nested commutators with the innermost operator held fixed provides a complete basis of thermal correlation functions (written down in (3.22)). That is, there are (n − 1)! such correlators and they are not related by either generalized Jacobi or KMS relations. We refer to this set of correlators as the causal basis.
• The KMS relations take a very simple form when expressed in terms of fully nested correlators. This statement can be argued for in a couple of different ways. The first involves writing down thermal Jacobi type operator identities, tJacobi relations, which can be proven using thermally deformed commutators and anti-commutators (see Appendix B). Alternately, one can recurse the KMS relations taking into account the generalized Jacobi relations of [18] to obtain a sequence of iterated KMS relations.
• We solve the KMS relations explicitly by giving formulae that express classes of correlation functions (such as thermal Wightman functions, nested correlators, and advanced/retarded Green functions) in terms of our causal basis (see §3.4).
• All told, the causal basis captures the essence of the KMS relations most efficiently. We argue that it is a very useful alternative to a basis of Wightman functions, since it implements nice causal properties and conforms with various analyses of thermal spectral functions in the literature.
The outline of the paper is as follows: In §2 we explore general features of thermal Wightman functions exemplifying the cyclic structure and construct the master Wightman function which generates all the KMS relations. 3 In §3 we describe various features of nested commutators, construct the causal basis, and explore iterated KMS relations. §4 uses the simple harmonic oscillator to illustrate the general features of our the discussion in a somewhat explicit manner. Finally, in §5 we turn to the OTO characterization of the KMS relations. We have collected in the appendices several technical details which prove useful to verify the statements in the main text. In addition, for completeness we also check that the thermal identities we derive are consistent with the relations obtained in the thermal quantum field theory literature (Appendix A).

Wightman correlators and the KMS condition
The Wightman correlation functions of interest are correlation functions of Heisenberg operators O(t) = U (t 0 , t) † O U (t 0 , t) with no prescribed time-ordering. We will consider generic n-point functions, and w.l.o.g. fix the temporal insertion points to be ordered t 1 > t 2 > t 3 > · · · > t n and simply permute the operators of interest. The Wightman basis of n-point functions is then given by where S n denotes the group of permutations of n objects and O σ(i) ≡ O σ(i) (t σ(i) ). This accounts for the n! possibilities of time-ordering of n-operators.
In what follows we will find it convenient to simplify notation -we refer to the operators by their temporal insertion points thereby abbreviating, O j (t j ) ≡ j, which helps declutter formulae below.
While Wightman correlation functions can be studied in any given state, in this paper we will focus exclusively on thermal correlation functions in equilibrium. We therefore take our quantum system to be in a thermal density matrix The correlation functions of interest are then

3)
A special feature of the thermal density matrix is the fact that it involves the Hamiltonian, and can be interpreted as evolving the system by an imaginary amount t = −i β. This fact can be encoded in the Kubo-Martin-Schwinger (KMS) conditions [1,2]. Formally, one may state the KMS conditions in terms of the Schwinger functions of the Euclidean theory.
Two-point functions are required to be analytic in a strip in the complex time plane: with t C = t σ(1) − t σ(2) − i t E the KMS condition requires that the Green's functions are analytic in the strip {t C ∈ C : 0 < t E < β}. Hence for two Heisenberg operators A(t) and B(t) which are elements of the algebra of observables, the two-point thermal correlator obeys the periodicity condition for bosonic operators 4 A and B. We used here conjugation of A by the density matrix operator ρ T and cyclicity of the trace. This motivates us to define the KMS conjugate of an operator: We are now in a position to state the general KMS condition for n-point functions. Consider a correlation function of the form We have allowed complex time shifts by k i β being mindful of the analyticity of the Euclidean correlator to only continue the imaginary argument into the lower half plane. This analytic continuation is valid provided as we read from left to right in the above correlator, the imaginary part of the time argument always decreases, viz., k n + 1 ≥ k 1 ≥ k 2 ≥ . . . ≥ k n , independent of the real time ordering, cf., [13]. 5 4 The generalization to fermionic operators is straightforward. We will exclusively work with bosonic operators for simplicity. 5 This can be seen by inserting a complete set of states between operators and demanding the complex phase factors thence generated are exponentially damped at high energies.
where we have assumed that all the time arguments fall into the admissible domain.
More generally, we can consider the correlation function defined by the permutation σ ∈ S n with a sequence of imaginary time excursions characterized by k i for each of the operators., viz., The KMS condition relates correlators within a cyclic orbit of the permutation σ giving relations of the form: for p = 1, · · · , (n − 1). (2.9) Thus given a particular time-ordering (specified by the permutation σ) the KMS conditions relate it to its cyclic permutations, with operators that pass through the density matrix being shifted in imaginary time by an extra unit along the thermal circle. This clearly, breaks up the n! Wightman correlators into (n−1)! equivalence class each comprising of n correlators related by a KMS relation and imaginary time shifted arguments. 6 To make these statements more explicit, we will give some examples below. To be concise, we will limit ourselves to the case where all the correlators we will consider have the complex time shifts either 0 or −iβ, i.e., k i ∈ {0, 1}.

KMS relations in time domain
The KMS relations can be written out explicitly for various low-point Wightman functions. For instance, we have for n = 2 and n = 3 the following decomposition of the correlators: Two-point functions: At the two-point level there is the standard relation which we express as 12 = 2 β 1 = 1 β 2 β (2.10) Note that this is the same as the relation quoted originally in (2.4), and in particular means that once we take into account the KMS condition there is a single two-point function in thermal equilibrium.
Three-point functions: One can similarly carry out the exercise for three-point functions.

KMS relations in frequency domain
In frequency space, KMS relations take a much simpler form since we can write: 7 By virtue of the fact that the imaginary time shift factors out, we end up with simple forms for the KMS relations in frequency space: for p = 1, · · · , (n − 1) . (2.13) This is self consistent since n i=1 ω i = 0 owing to time translational invariance. These relations are again best illustrated with examples for low values of n.
Two-point functions: The n = 2 KMS relations (2.10) take the form: (2.14) Three-point functions: The n = 3 KMS relations in frequency space can again be immediately written down from the earlier expressions in (2.11). They take the much simpler compact form (2.15) 7 We will abuse notation and not distinguish between the time and frequency domain incarnation of the operators. By convention we take O(ω) =´dt e i ω t O(t).
Four-point functions: Finally, we record the n = 4 KMS relations in frequency space. To wit, where σ ∈ S 3 is any permutation acting on 2, 3, 4.

The generator of KMS relations
The cyclic symmetry inherent in the KMS condition is very clearly visible in the frequency domain, where associated with a given cyclic permutation we have associated a prefactor involving the frequencies of operators that have been conjugated through the density matrix measured in thermal units. This simple transformation law suggests that we can recover the KMS relations from a set of master functions that diagonalize the cyclic permutations. Consider a Wightman correlation function specified by a given permutation σ ∈ S n . We consider the action of the cyclic group Z n on the sequence σ(1)σ(2) · · · σ(n). Denoting by π an element of Z n , we would end up with map Taking a weighted average of the cyclic permutations we construct the master function for the permutation σ ∈ S n : Ψ σ (1, · · · , n) ≡ π∈Zn πσ(1) πσ(2) · · · πσ(n) π∈Zn e − β n j j ω πσ(j) , (2.18) where π ∈ Z n is a cyclic permutation. Then the KMS relations can be stated in terms of recovering a given Wightman correlator in frequency space directly from Ψ σ . One finds: σ(1) · · · σ(n) = Ψ σ (1, · · · , n) e − β n j j ω σ(j) . (2.19) Note that the master function Ψ σ is indexed by a permutation σ ∈ S n . Since we are averaging over the Z n orbit of σ, we only need to consider (n − 1)! elements of S n . Combined with the action of π we would end up covering all n! permutations.
There is a physically interesting way to write Ψ σ . We claim that which makes the cyclic invariance manifest. For n = 2 this spectral function is naturally interpreted in the thermofield double construction with factors ofρ 1 2 T appearing between the forward and backward evolutions. Similarly, the chaos correlator in [20] was computed using an Euclidean regulator which inserts powers ofρ 1 4 T between the operators in a 4-point function. More generally, inspired by these developments, [25] considered arbitrary combinations of commutators and anti-commutators of two operators, [A, B] and {A, B}, to derive a set of bipartite fluctuation-dissipation relations. Their results follow from the statements mentioned above as the reader can quickly infer.
The statement of (2.19) is that Ψ σ is the natural object that symmetrizes the KMS relations amongst the (n − 1)! independent Wightman functions, and should be viewed as providing a suitable basis for the spectral functions. 8 Once again it is helpful to view these relations through some examples of low point correlation functions.

Generalized fluctuation-dissipation relations
While Wightman functions encode all the information inherent in the KMS relations quite succinctly, it is often useful in physical applications to consider correlation functions of various combinations involving nested sequence of commutators and anti-commutators. This basis was referred to as the nested basis in [18] and its physical utility is that causal response functions are determined in terms of nested commutators. Unlike the Wightman basis the set of nested correlators has redundancies. For instance even without including imaginary time shifted operators, it was noted in [18] that one naively has 2 n−2 n! nested correlators, since each of the n! permutations allows for (n − 1) binary choices of either a commutator or an anti-commutator (and we only need to consider alternating permutations owing to an overall sign flip resulting from swapping order of operators in a commutator). In this section we resolve this redundancy, leading to a minimal set of (n − 1)! spectral functions. Generalized fluctuation-dissipation relations encode how to express other correlators (such as Wightman functions, or generic nested correlators) in terms of these spectral functions.

Notation
To simplify some of the formulae that appear in the sequel, we introduce some new notation which we collect here for quick reference: • We will be interested in nested correlators with both commutators and anti-commutators.
To avoid clutter of brackets, we will write out a single unbroken string of operators with the innermost operator in the sequence being the left-most entry in the sequence. We specify anti-commutators by marking the outer-most entry with a '+' subscript, so that where we remind the reader of our shortcut notation O j (t j ) ≡ j. For example for 3-point correlators we would have combinations of the form among others. In particular, note that a fully nested commutator will have no + subscript markings, i.e., • Various thermal factors which appear in the formulae will also be abbreviated: Thermal factors : 3) The following relation between Bose-Einstein factors in correlators with n insertions is very often useful to simplify various expressions: This is a consequence of energy conservation, ω 1 + · · · + ω n = 0. For example, in 4-point functions we have f 1,2,4 = −(1 + f 3 ) etc.. An analogous statement for thermal factors can be written as • In the context of fluctuation-dissipation relations, we will use a shortcut for the combination that appears in the n = 2 FDT: • When necessary we capture both commutators and anti-commutators by the definition • Inspired by various relations which we encounter, we also introduce a new bracket notation. Define the thermally deformed commutator and anti-commutator as follows: Note that the definition singles out the inner-most operator in the (anti-)commutator and reduces to the standard commutator and anti-commutator when β → 0. This thermal commutator is used in Appendix B and provides an alternative route to deriving the fluctuation-dissipation relations for nested thermal correlators.

Relations among nested correlators: examples
It is helpful to first see some explicit examples (which are well known) to understand the rationale behind the use of the nested correlators.
Two-point functions: The usual presentation of the KMS relation is in the form of a fluctuation-dissipation (FD) condition. Recall that response functions, which involve commutators, naturally diagnose dissipation (and transport), while the symmetrized Green's functions capture fluctuations [3]. Starting with (2.14) one can by taking suitable linear combinations arrive at the standard form of the FD relation, viz., We have chosen here to use the KMS relation to solve for the fluctuation measure in terms of the commutator. Whilst it is conventional to do so, this relation is singular when β ω 1 → 0 for N 1 diverges. What this really means is that in this limit the two operators simply commute, as would have been manifest had we solved the KMS relation to give instead In what follows we will typically solve the KMS relations and present the results in terms of nested commutators. The reader should be aware then that the relations hold modulo contact-terms, which could arise when the statistics factors N i diverge.
Three-point functions: The analogous story for 3-point functions is a bit more involved. We have 6 permutations of the operators, and two possible nestings of each ordering into a commutator/anti-commutator. Of these 12 nested correlators, 6 can be eliminated by use of generalized Jacobi relations, called sJacobi relations in [18]. 9  to always move the operator 1 to the left-most position (a pattern which we will repeat below). While we have written only four of the remaining 10 nested correlators, the remaining six may be obtained by permutations of the above. It is also useful to record here the expression for the thermal 3-point Wightman correlator in terms of the two independent nested basis elements identified above. To wit, (3.13) The remaining Wightman functions can be similarly expressed after permuting the operators on the l.h.s and using the sJacobi relation in Eq. (3.12) to eliminate [ [231] ]. For example, we can write the following set of relations: , , which involves all possible positions of the operator 1, so the remaining 3 Wightman functions are obtained by permuting 2 and 3, which doesn't require any use of Jacobi identities.
Four-point functions: Let us now turn to 4-point functions. There are 96 possible correlators, obtained by taking the 24 permutations of operators and inserting them into commutators and anti-commutators. We now show that these are determined by 6 independent structures after accounting for the sJacobi and KMS relations. We will choose to express all the correlators in terms of the basis generated by [ [1σ(2)σ(3)σ(4)] ] with σ ∈ S 3 (we show that this is a basis in Appendix C). 10 Explicitly, all commutator/anti-commutator nestings can be written in terms of this basis using the KMS relations (cf., the caveat around (3.10)):  Chaos correlator: The chaos correlator 2 is an out-of-time-order correlator as a measure of chaotic behaviour in a thermal quantum system [20,29,30]. We consider the following parameterization: (3.18) Using the aforementioned results we can express this correlator in terms of the thermal spectral functions in the Wightman basis or the nested commutator basis, respectively as: where C(t) = C(0, t, 0, t). As explained in [20] it is useful to consider the regulated correlator which is simply the spectral Wightman 4-point function. Letting F (0, t, 0, t) = F (t) (3.21)

Towards a causal basis
After having seen various examples at low-point order, let us finally explain the general construction for n-point functions. To understand the situation better we should first enumerate the space of nested correlators and ascertain the relations we expect amongst them. Once we do this we will be in a position to construct a useful basis in terms of which to express the full set of correlators. As the reader might anticipate from the results in §3.2, we will now argue that the appropriate basis is provided by a subset of fully nested commutators (which therefore pick out a basis with nice causality properties). We only give the general arguments below, and relegate most of the details to Appendices B and C, respectively.

Counting
Any n-point nested correlator should be expressible as a linear combination of (n−1)! spectral functions after using all sJacobi relations and KMS relations. First, recall from the beginning of this section that there are 2 n−2 n! nested correlators involving (anti-)commutators [18], but our analysis in terms of Wightman correlators reveals that these should be expressible in terms of (n − 1)! spectral functions for a thermal state.
Of the set of nested correlators (which allow for both commutators and anti-commutators), one has 1 2 n! = 1 2 n (n − 1)! nested commutators. It would already be sufficient to pick out from amongst these (n − 1)! nested commutators and declare them to be our basis. For this to be true, we need to show that the remaining nested correlators are related to our basis choice by either sJacobi relations or through KMS relations.
We recall that amongst the 2 n−2 n! nested correlators, there are (2 n−2 − 1)n! sJacobi relations. 11 They are operator relations in that they hold by virtue of the products of operators generating a free algebra with a graded commutator. The detailed proof of this statement and examples of relations can be found in [18]. Accounting for these sJacobi relations we would end up with a basis of n! correlators which agrees with the Wightman count. However, for thermal correlators, the KMS conditions further reduce the number of independent correlators to (n − 1)!. Thus there should be n! − (n − 1)! further relations. The quest now is to find a useful way to characterize the full set of relations in one go.

The causal thermal commutator basis
In §3.2 we gave various examples, expressing low-point Wightman functions and nested correlators in terms of nested commutators [ [1 σ(2)σ(3) · · · σ(n)] ] . We now make a completeness statement about this representation: Lemma 1: A useful basis of thermal n-point correlators is given by the nested commutators for all σ ∈ S n−1 .
We refer to this as the causal basis. Consequently, for a thermal density matrix, any nested correlator (and therefore any Wightman function) can always be expressed as a linear combination of causal spectral functions involving only nested commutators of the form (3.22), which we will make explicit below.
The equations expressing any thermal Wightman function in terms of the causal basis are equivalent to the Fluctuation-Dissipation (FD) relations for n-point functions (see below). For correlators that can be computed by the Schwinger-Keldysh path integral these were given by [7] using the av(erage)-dif(ference) basis of contour correlators. Since the latter can be expressed in terms of the nested correlators using the Keldysh rules, this is an equivalent presentation. We indeed check in Appendix A that their expressions match with ours for low-point functions. Later in §5 we will see that for n ≥ 3 there are non-trivial KMS relations relating out-of-time order correlators which transcend the Schwinger-Keldysh ordering (cf., also [25]). There are a couple of different ways to illustrate the veracity of our claims above: • The most straightforward way is to explicitly demonstrate how every Wightman function can be expressed in terms of the above causal basis. This is done in §3.4 and can be seen as a generalization of FD relations. A related approach is to write any one of the 2 n−2 n! nested correlators in terms of the causal basis (see §3.4.1).
• Another approach is to work directly with generalized Jacobi relations to show completeness of the basis (3.22). One can first account for the sJacobi relations of [18] (which are valid in any state), and then subsequently isolate the extra relations implied by the thermal nature of the density matrix. We illustrate the latter method in §3.5.
• Finally one can mimic the strategy followed in [18] more closely and encapsulate the set of relations implied by the sJacobi relations and thermal KMS conditions in one swoop. This leads to the notion of tJacobi operator relations, which are thermal deformations of Jacobi relations and their generalizations. We explore these in Appendix B.

Wightman functions in terms of causal basis
It is possible to generalize (3.13), (3.17) to n-point functions and give a general formula that expresses Wightman functions in terms of the causal basis of correlators The result can be stated as follows (the proof can be found in Appendix C.2): 12 Lemma 2: An n-point thermal Wightman correlator can be written in terms of nested commutators as 12 · · · n = ρ∈S n−1 where ρ ≡ (ρ(2), . . . , ρ(n)) are understood to be permutations acting on the set {2, . . . , n} and the numberss (ρ) i ∈ {0, 1} describe the run structure of the permutation ρ. We define the run structure of a permutation ρ as follows. We call the index i an ascent if ρ(i) < ρ(i + 1) (and a descent otherwise). The run structure is then the collection of ascents and descents, which culminate in peaks and valleys, respectively, along the string of elements being permuted (see [31] for a recent attempt to enumerate permutations by their run structure).
The numberss are defined in terms of the Heaviside step-function scanning the run structure of the permutation ρ: For reasons discussed around Eq. (3.10) one has to pay closer attention to the loci where the Bose-Einstein functions conspire to diverge, viz., in the limit β ω i → 0. One can check that the linear map between the Wightman correlators and the nested commutator basis has an eigenspectrum built from the set of products of the statistical factors. More precisely, we find that the eigenvalues of the transformation (3.23) are given by 13 for all permutations σ ∈ S n−1 and all indices 2 ≤ i 1 ≤ . . . ≤ i N < n for any 1 ≤ N ≤ n − 2. In words: the set of (n − 1)! eigenvalues is constructed by distributing all frequencies ω 2 , . . . , ω n , 12 We thank Simon Caron-Huot for sharing valuable insights from his unpublished work which was extremely useful in obtaining this result. 13 We have inferred these eigenvalues by analyzing low-point examples (up to and including n = 5). In the formula below we do not give the degeneracies. Of the (n − 1)! eigenvalues, we empirically find and so on. We will not here attempt to prove these statements, since we are going to eschew these singular loci in our discussion.
arbitrarily over an arbitrary number of f-factors. There exist divergent eigenvalues whenever any subset of the frequencies adds up to zero, i.e., whenever there exists a permutation σ such that β k≤n i=2 ω σ(i) = 0. At these points the relation given in (3.23) breaks down and one would have to account for additional contact terms on the right hand side. In the discussion below we will assume that we are always away from these singular points in frequency space and refrain from writing down the contact terms. We note however that the origin of such terms was already inferred in the analysis of [13] based on the analyticity domains of the Euclidean correlator.
Lemma 2 can be seen as one way to state n-point function fluctuation relations. It provides an explicit and complete implementation of KMS conditions: a trivial permutation of labels {2, . . . , n} on both sides of (3.23) gives a set of (n−1)! relations of the same form, which then provide the matrix that transforms between Wightman functions and the causal basis (in Appendix C.2 we give an explicit formula for these permutations). Note also that this matrix is obviously invertible -the converse problem of expressing nested commutators in terms of Wightman functions is rather trivial to solve (one simply expands out the commutators), and an explicit formula can be found in Lemma 5.

Application: nested correlators in terms of causal basis
As an application of Lemma 2, we can give the relations that give the complete set of 2 n−2 n! nested correlators in terms of the causal basis. This again amounts to a complete set of generalized fluctuation-dissipation relations. To achieve this, we expand a nested correlator [ [1 2 ε 2 · · · n εn ] ] in terms of Wightman functions, bring them to canonical form (i.e., with O 1 being the left-most operator) using the KMS condition, and finally apply Lemma 2 to each term. The result is Lemma 3: An n-point nested correlator can be written in terms of nested commutators as with thermal factors given by For more details on the derivation of (3.26), we refer to Appendix C.2. Rewriting the Bose-Einstein factors in terms of N i , one can use this formula to verify, e.g., the four-point function fluctuation relations stated in (3.15).
Lemma 3 can again be interpreted as solving the generalized n-point sJacobi relations and the KMS relations at the same time to explicitly write all nested correlators in terms of our causal basis. Let us also note that the Keldysh rules (and their generalization to k-OTO contours of [18]) can be used to write standard retarded and advanced thermal n-point functions in terms of the objects appearing on the left hand side of Lemma 3. Therefore, (3.26) allows us to express standard response functions in terms of the causal basis. We refer to Appendix A for a more detailed consistency check with existing literature on this topic.

Iterated KMS relations
Let us now turn to an alternative approach to verify Lemma 1, which entails an implementation of KMS relations at the level of generalized Jacobi relations.
We will start with known KMS relations and employ the sJacobi relations to whittle down the space of nested correlators. We will see that a naive application leads to an over-complete set of relations, which then implies that there ought to be non-trivial identities amongst the various relations we construct. These turn out to the interesting identities amongst thermal correlators, which is the raison d'etre for the construction we describe herein.
The general discussion is best illustrated by an example. Consider the n = 2 FD theorem (3.9) for two operators X and A which follows from the KMS relation, viz., 14 (3.28) As this holds for any X in the operator algebra it in particular holds for the choices X = [C, B] or X = {C, B}. This then implies These can be further simplified using the n = 3 sJacobi relations [18], This exercise shows that all the n = 3 nested correlators can be expressed in terms of n = 2 nested commutators and upon using the sJacobi relations, one isolates the n = 3 FD relation.
One therefore suspects that all the n-point fluctuation-dissipation relations can be derived this way, i.e., using (3.9), iterating, and applying the appropriate set of sJacobi identities. This makes clear the primacy of the n = 2 FDT (3.9), since the sJacobi identities are but operator identities owing to the presence of a graded commutator (as explained in [18] this simply follows from thinking about the operator algebra as a free Lie algebra). In particular, all the higher point relations (for example those obtained in [7]) are generated from the n = 2 FDT algebraically, with no further dynamical input. 14 We emphasize again that the expressions we write down have potential ambiguities owing to the divergence of the statistics factor N as βω → 0. See the discussion around Eq. (3.10) and Eq. (3.25).
Let us make this precise. Consider the n = 2 FDT written in the form (3.28). Now take X to be a nested commutator involving n − 1 operators O i with i ∈ {1, · · · , n − 1}. The number of such identities we can write down is n 2 (n−1)−2 (n − 1)! = 2 n−3 n!, the count following from the number of nested structures and cyclic permutations of the resulting n objects. However, these cannot all be independent. For one we have to account for sJacobi relations: if X = j Y j = 0 is an sJacobi relation between some set of nested correlators Y j , then we trivially satisfy j Y j A + +N A Y j A = 0. So this is not a new KMS relation and we should mod out by such identities which we can enumerate to be n (2 (n−1)−2 −1)(n−1)! = (2 n−3 −1)n!. Subtracting this from the previous count, we therefore see that the number of independent relations after accounting for the sJacobi identities is n!. This is however too many, since there are only n! − (n − 1)! independent Wightman functions after accounting for the KMS relations. Therefore among the many relations we derive by the above logic, there ought to be a further (n − 1)! relations among the identities obtained via use of the n = 2 FD theorem and the sJacobis. We shall call these relations among relations, the iterated KMS relations.
While the general construction of these (n − 1)! iterated KMS relations is clear, they are best understood by working them out explicitly at low orders of n. To write these down explicitly, we remind the reader of the notation (3.6) which emphasizes the primacy of the n = 2 FDT, which we therefore encapsulate by the new symbol: (3.32) We can now state in terms of this master relation, the iterated KMS relations quite succinctly.

Two-point iterated KMS relations:
The well-known 2-point FDT in our new notation reads as Three-point iterated KMS relations: We claim that the n = 3 relations are: Writing this out explicitly, we have (3.35) These relations can be proven by showing that any one of them is equal to

(3.37)
Four-point iterated KMS relations: We find the following set of relations for 4-point functions

Illustration: Harmonic Oscillator
We now illustrate the general relations obtained hitherto with the simplest possible model one can consider -a harmonic oscillator. Let X(t) be the coordinate of the particle moving in a harmonic oscillator of frequency µ. We will consider correlation functions of composite operators of the form X a (t) for a ∈ Z + . We will exhibit various KMS relations primarily for three point function, by computing all the Wightman functions. In passing, we will also show how one can also obtain these via an analytic continuation from Euclidean correlators, thus verifying explicitly some of the abstract statements encountered herein and in the literature.

Thermal Expectation values of Wightman correlators
We will be interested in a thermal correlation function where |n are the energy eigenstates of energy E n = µ (n + 1 2 ) and the partition sum is .
The position operator can be expanded in terms of creation and annihilation operators as (setting mass m = 1 for simplicity) 2) The operators act on the Hilbert space as and satisfy [a, a † ] = 1. Using this, we can compute thermal expectation values and explicitly verify various identities derived in previous sections. We provide further computational details in Appendix D.
Two-point functions: Using the above definitions, it is easy to show that (with t ij = t i −t j henceforth to account for time translational invariance) Summing over n with weights as in (4.1), we can compute the thermal propagator and explicitly verify the KMS condition: 15 2 µ G M (t 12 ) ≡ 2 µ X(t 1 ) X(t 2 ) β = e iµ t 12 + e βµ e −iµ t 12 e βµ − 1 = 2 µ X(t 2 ) X(t 1 + iβ) β . (4.5) In frequency space this gives the familiar form of the spectral function, for .

(4.6)
Three-point function: We can compute a three-point function using similar methods. Because we are dealing with a free theory, we can express the result in terms of the propagator: Similarly, one can evaluate other orderings. For example, (4.8) One can now explicitly check that the KMS conditions given in (2.11) hold. For instance, It is instructive to pass onto the nested correlators which (as expected) turn out to be much more simple than the Wightman correlators. For example (4.10) Four-point function: Finally, let us also record a couple of four point functions (4.11)

Euclidean correlators
The propagator: The Euclidean two-point function can be derived by summing the free propagator over the Matsubara modes with frequency ω n = 2π n β , n ∈ Z leading to the Green's function [32] (see Appendix D for details): e βµ e −µ|τ | + e µ|τ | e βµ − 1 (4.12) The analytic continuation to real time can again be performed to obtain (4.5) for Likewise the time ordered correlator evaluates to 2µ T X(t 1 )X(t 2 ) β = e βµ e −iµ |t 12 | + e iµ |t 12 | e βµ − 1 (4.14) Three-point function: The Euclidean 3-point function of normal ordered operators then evaluates to The real-time Wightman correlators obtained earlier can then be extracted by different analytic continuations. For instance, (4.7) follows from while (4.8) results from the analytic continuation: Chaos correlator: Finally, let us draw a connection between the OTO 4-point function that provides a diagnostic for chaos [20,30] and the Wightman function evaluated at complex times. To do this note that for any operator We take V ≡ X and W ≡ X 2 to find: where f = 1 e βµ −1 . The analytic continuation in t of this function decays at large µ only in the strip | (t)| ≤ β 2 . This means that in a field theoretical setting with an infinite number of harmonic oscillators (with arbitrarily high frequencies) the above function exists only within this strip. We give further expressions for other 4-point functions and some tremelo 6-point functions in Appendix D.

OTO classification of thermal correlators
Since generic Wightman functions do not respect time-ordering, their computation necessarily involves timefolds or out-of-time order (OTO) path integral contours. While this point has been well appreciated for a long time now, only recently has there been a systematic attempt to classify the OTO correlators [18]. As discussed there, the n! Wightman functions can be computed by contours with at most n+1 2 timefolds. Moreover, given a particular time ordering there is a canonical proper q-OTO number with 1 ≤ q ≤ n+1 2 which provides the simplest representation of the correlator (i.e., the minimum number of timefolds required to represent the correlator on a contour).
While the discussion of [18] was for generic initial density matrices, as we saw above, various Wightman correlators are related by KMS conditions. These relations as we have described involved cyclic permutations of the operators, which do not respect the timeordering. 16 In fact, one can already see from the simplest example of 3-point functions in (2.11) that the KMS conditions do not respect OTO number: while 123 is time-ordered and can be computed from a 1-OTO contour, 231 necessitates a 2-OTO contour; these two correlators (with some arguments shifted in imaginary time) are related by a KMS relation. More generally, the proper q-OTO number for a thermal n-point Wightman function lies in the range 1 ≤ q ≤ n 2 . As the KMS relations map all cyclic permutations of operator insertions in a given correlator to each other, for an n-point function this ends up equating correlators of different proper-OTO number. The question we propose to answer is to understand the OTO classification of the correlators in any cyclic orbit of the KMS action. Below we add some additional structure to the classification scheme described in [18] to allow for a cleaner discussion of cyclic features inherent in the KMS relations, and show how to decompose the space of n! Wightman functions into (n − 1)! independent spectral functions, using OTO-contours and KMS relations. Physically, it is interesting to understand these relations to ascertain which of the OTO correlators carry non-trivial and independent information in the thermal state.

Review of OTO classification
Before proceeding it will be helpful to have a quick overview of the nomenclature and results of [18], which we invoke extensively below. A k-OTO functional integral computes the generating function for thermal Wightman correlators and is defined to be: A pictorial representation of the functional integral contour in the complex time plane is given in Fig. 2.
When computing correlators, we essentially insert operators along this contour at the appropriate time, with a contour ordering prescription that relates directly to the permutation σ ∈ S n corresponding to the Wightman correlator of interest. The contour should be viewed as an abacus, with operators free to slide about from one level to another, as long at they are unobstructed by other operators (or the density matrix). For instance: Given a particular σ ∈ S n , there is a primitive contour of minimal number of switchbacks, which is the proper-OTO contour for the said correlator. A proper q-OTO Wightman function is one that can be minimally represented on a contour with q timefolds (and cannot be computed using one with less than q timefolds). For a given n, one can show that proper q-OTO contours with q = 1, 2, · · · , n+1 2 are necessary. The upper bound can be easily seen by considering the completely oscillating or tremelo permutation, which is a sequence where insertion times alternately increase/decrease along the contour.
The proper q-OTO contours provide a partitioning of the set of n-point Wightman functions. In particular, a q-OTO contour computes g n,q of the n! time-ordering correlators, which have proper-OTO characteristic q. One finds: Furthermore, each such proper q-OTO correlator can be computed from the generating functional (5.1) in a number of h (q) n,k equivalent ways (akin to sliding relations as illustrated in the above picture). It is a non-trivial fact that this number only depends on (n, k, q), but not on the specific correlator and order of insertion points. We therefore have a decomposition of the (2k) n contour n-point functions computed by (5.1) of the form .

(5.3)
It is also useful in the sequel to know that for even n, there is a maximally out-of-timeorder correlator corresponding to a tremelo permutation. 17 Equivalently, the run structure of the permutation alternates between ascents and descents in a sawtooth pattern.
In the following it will be helpful to know of turning-points and turning-point operators along the OTO contour. A future turning-point is the turning point at the right end of Fig. 2, i.e., the junctions between There are q-future and (q − 1)-past turning-points. Finally, an operator inserted just before a turning-point will be referred to as a turning-point operator (an example is O 2 in the pictorial representation of G(t 4 , t 1 , t 3 , t 2 ) depicted above).

Generalities
Let us now turn to the interplay between the OTO classification and KMS relations. Consider any Wightman correlator labeled by a permutation σ ∈ S n . Having made our choice to let t 1 > t 2 > · · · t n we can let σ be the permutation which takes a completely time ordered correlator to the correlator in question; eg., 123 . . . n is mapped onto 2134 . . . n by σ = (12).
While the physical picture should be clear from the discussion below, it is worthwhile describing our results at the outset in some formal terms. In standard combinatorics, permutations of n elements can be classified based on their run structure. Assuming the elements which are being permuted are ordered, runs are usually subsequences of ascent or descent as described in §3.4. An ascent of permutation σ ∈ S n acting on the ordered set {1, 2, · · · , n} is a sub-sequence of σ(1) σ(2) · · · σ(n) where σ(i) < σ(i + 1), while a descent requires σ(i) > σ(i + 1). It should be transparent that the OTO classification of Wightman correlators is isomorphic to the run structure of the permutation; each ascent corresponds to a forward evolution, while a descent to the backward evolution, and the entire sequence can be decomposed into a set of ascents and descents. The discussion of [18] should be then interpreted as an analysis of linear run structure, since one is only interested in the absolute ordering of the operator insertions.
The KMS conditions on the other hand, involve cyclic relations, which therefore implies that we should look to classify the cyclic run structure of permutations. This can be understood as follows: in the OTO contour of [18] (see also [19]), it was noted that the future and past turning-points of the contour should be interpreted as the insertion of an identity matrix, since the gluing conditions on the contour require that the future-most end of a forward directed segment of the contour, gets mapped back without change to the future-most end of a past directed segment. For generic density matrices, once the forward/backward evolution is complete, one takes the trace.
However, in the case of thermal density matrices, one can view the past-most point on the initial and final legs of the contour (denoted 1R and 1L in Fig. 2) as also being joined together through an imaginary time evolution. In other, words a thermal k-OTO contour is a thermal circle, with forward/backward excursions in real time. 18 Operators inserted on the contour can then be slid around, as long as they do not encounter another operator blocking their way. Sliding an operator through the density matrix (i.e., along the dashed line of Fig. 2) will result in a shift of the operator argument by −i β. In keeping with our conventions, we will only allow sliding counter-clockwise through the density matrix.
In the thermal state in addition to the future/past turning-points, we now also have the density matrix turning-point. Note that in the infinite temperature limitρ T reduces to the maximally mixed state given by the identity operator, which results in the density matrix turning-point becoming isomorphic (modulo normalization) to a past turning-point. The cyclic run structure is now transparent, since we are required to decompose n! ordering into sequences that can be nicely arranged on the thermal circle, with real-time forward/backward evolutions.
The KMS relations can then be simply understood in terms of sliding operators on the closed thermal OTO contour. Algorithmically we proceed as follows: • Begin with a proper q-OTO representation of the n-point Wightman function associated with the permutation σ.
• Embed the q-OTO contour into a redundant (q + 1)-OTO contour by appending one further timefolds at the bottom of the contour, i.e., inserted just before the density matrix. 19 • Starting at the very past of the 1R contour, take operators sequentially through the density matrix, traversing counter-clockwise. The operators should be brought to the first empty contour having the same direction of time evolution as the original contour from which they were taken.
• Each operator transition through the density matrix gives a KMS related correlator. At the end of the sliding, we canonicalize the OTO by erasing the redundant legs (adjacent forward/backward evolution).
We illustrate the above algorithm pictorially in Figure 3. An immediate consequence of the construction is that it makes manifest the fact that starting with a Wightman correlator with proper OTO number being q, in the cyclic orbit of σ, we are guaranteed to find proper OTO numbers being either q or q + 1 and none other.

OTO sliding and KMS relations
Let us now turn to giving a more concrete picture which will result in a decomposition of n-point Wightman functions into equivalence classes under KMS relations. To this end we will introduce some basic objects, q-and δ-lists which will be suitable sequences of integers that capture the OTO structure of a permutation and its cyclic cousins. The analysis below will complement the general picture above and result in an explicit formula for the count of KMS relations.

δ and q lists
Given a permutation σ ∈ S n , we associate with it an ordered set of proper OTO numbers {q 1 (σ), q 2 (σ), . . . , q n (σ)} which we will call the q-list. These are the proper OTO numbers of Wightman correlators associated with σ and its cyclic cousins. This list suffices to determine the OTO structure of KMS relations, for q j (σ) is the proper-OTO number when density matrix is sandwiched between the j th and (j − 1) st operators. It is actually more efficient to consider a more primitive object which characterizes the OTO structure of the Wightman correlator. For a given permutation σ, we begin by defining a binary sequence of n numbers which we refer to as the δ-list: We take δ i (σ) = 1 if the i th operator in the σ Wightman correlator is a past turning-point operator. 20 Otherwise, we take δ i (σ) = 0. Note that for i = 1 and i = n, we check the past turning-point by cyclically moving them away from the edges (since by definition, edge operators can never be a past turning-point operator).
The δ-list determines the q-list, since the total number of (non-edge) past turning-point operators in a correlator should be one less than its OTO number. Hence, where we take δ 0 (σ) = δ n (σ) and have defined q 0 ≡ 1 + n i=1 δ i (σ). The relation (5.4) can then be used to compute the q-list defined as the sequence {q 1 (σ), q 2 (σ), . . . , q n (σ)} for a given permutation σ. This analysis gives a simple way to obtain the proper OTO numbers associated with all correlators related to a given permutation σ by KMS condition. We illustrate this construction with some low point examples n ≤ 4 in Tables 1, 2, and 3 respectively. Figure 3 provides another detailed example. 20 We remind the reader that the past turning-point operator lies to the past of both its neighbours [18].
{1, 1} 21 Table 1. q-lists and δ-lists for two-point functions {2, 1, 1} 231 , 132 Table 2. q-lists and δ-lists for three-point functions We are interested in enumerating how many KMS relations can be represented on a proper q-OTO contour. For example for the three-point functions, KMS conditions relate the first entry in the third columns of Table 2 with each other, i.e., they relate the two 1-OTOs (e.g., 312 , 123 ) with a 2-OTO ( 231 ). Similarly KMS conditions also relate the second entry in the third columns of the tables with each other, viz., they relate the two 1-OTOs (i.e., 321 , 213 ) with a 2-OTO ( 132 ). Thus by knowing the q-list we will have identified how the KMS relations cycle through q(σ)-OTOs and (q(σ) + 1)-OTOs.

Properties of q and δ-necklaces.
Given the q-and δ-lists, we can now define the notion of q-and δ-necklaces. All the q (or δ) lists obtained by a cyclic permutation of the operators in the correlator is collectively termed as q (or δ)-necklace. As noted these are equivalence classes of lists once we account for the KMS relations in thermal states. We now proceed to give a concise way of thinking about these objects and also note some interesting and useful properties of these necklaces.
Any δ-necklace must obey the property that two successive entries can never both be 1 (because two successive operators cannot both be in the past of each other). In fact, if the total number of 1's in a δ-necklace is p, then by performing cyclic shifts, one can always bring  Table 4.
q-necklaces and KMS relations: These abstract considerations can now be put to use to understand some features of the KMS relations. A couple of obvious properties of q-necklaces are useful in this regard: • Since KMS conditions relate correlators within the q-necklace, it is now apparent that they relate a proper p-OTO correlator either to p-OTO or to (p + 1)-OTO correlators as promised.
• For n = 2p, there are special necklaces where all correlators related through the KMS conditions are proper p-OTO. These correspond to cyclic orbits of tremelo permutations where every alternate operator is a past turning-point operator, cf., last row of Table 3.
These considerations make clear that the δ-and q-necklaces of a permutation are sufficient to encode the KMS relations of interest.

Necklace degeneracies
Having identified the utility of working with the necklaces, we now turn to computing the number of necklaces we can form for a given n. The count can be efficiently organized in terms of ascertaining the number of potential necklaces that are allowed for a given n, and fixed proper-OTO number p, and thence determining the degeneracy with which each is encountered. We now give the essential ideas behind these counts, noting that the general picture should be such that at the end of the day we get a partitioning of (n − 1)! (the number of independent thermal Wightman functions) into necklace sets.
Enumerating the necklace types: In order to enumerate δ-necklaces, it suffices to count configurations with p 1's. This is a standard problem in Polya theory. We need the count P(n, p), corresponding to the partition of n into p objects, with each object bigger than 1, but ordered up to cyclic permutations. This counting also gives the possible independent q-necklaces since the latter is inferred from the δ-necklace directly.
Necklace degeneracies: Given a δ-necklace {m 1 , m 2 , . . . , m p } we will assume that we have somehow chosen some canonical representative for each {m 1 , m 2 , . . . , m p } which we denote as {m 1 , m 2 , . . . , m p } c . We now enumerate the number of such necklaces, which will give a complete decomposition of the n! Wightman functions in terms of KMS cognizant parts.
Consider then one of the elements among the d n,p These degeneracy factors encode the basic information about how many times (modulo symmetry factors, see below) a canonical δ-necklace shows up within a KMS orbit. In summary, we have provided an explicit decomposition of the n! Wightman functions into equivalence classes under KMS relations. This was aided by encoding the OTO structure into the δ-and q-necklaces. The key point to note is the degeneracies with which these necklaces appear, and understanding them gives a complete decomposition of the Wightman functions into KMS equivalence classes, as has been explicitly enumerated for some low-point functions in Table 4. We see explicitly from there that for each n-point function, there are q-necklaces of all proper p-OTO types with p ≤ n+1 2 . Further, each necklace contains a set of OTO numbers differing by at most one unit. However, for odd n, the highest proper-OTO number p = n+1 2 lies in the KMS orbit of p = n−1 2 , thereby informing us that the KMS relations are making the higher OTO correlator redundant in the thermal state.

Discussion
We have primarily focused on synthesizing known features of thermal correlation functions and arguing that they are best understood in the space of out-of-time-order (OTO) observables. While traditional presentations avoid the OTO correlators by explicitly convolving the KMS condition with time reversal or a CPT transformation to restore operator ordering, the generalization as discussed herein allows for a simpler interpretation. The KMS condition acts by cyclically permuting operators in a given correlator. Among other things we have shown that the proper-OTO numbers of all correlators in the KMS orbit of a given Wightman function will at most differ by unity. We also 'solved' the KMS relations explicitly, by constructing a causal basis for the independent thermal correlation functions in terms of the fully nested commutators (with one operator held in fixed position). We have given explicit formulae that express various complete classes of correlation functions in terms of this causal basis; these can be interpreted as the full set of generalized n-point fluctuation-dissipation relations.
There are several interesting avenues that are ripe for further exploration. For instance, one could attempt to understand the structures herein in kinetic theory and in explicit models. However, there are a set of conceptual questions that we have not addressed herein which are interesting to contemplate further:  Table 4. The list of allowed q-and δ-necklaces and their degeneracies (accounting for the symmetry factors). One can check that the sum of the entries in the last column for a given n is (n − 1)!, which confirms that the necklaces provide an OTO classification of thermal n-point Wightman functions according to KMS orbits.
• Thermalization of chaotic quantum systems is a subtle process. One relevant observable in the context of quantum chaos is the 4-point tremolo correlator (see below). We can ask: are there finer-grained or perhaps novel features of thermalization and chaotic behaviour, that are captured excluvively by higher 2k-point tremelo correlators? If so, are these features captured more simply in terms of the causal basis of nested correlators?
• How does the notion of thermal equivariance introduced in the context of hydrodynamic effective field theories [33], which are constructed in the Schwinger-Keldysh context, generalize to higher-point functions?
• What is the implication of 2-OTO thermal correlators which are related to Schwinger-Keldysh correlators in terms of transport?
The fact that KMS relations lead to relations between observables of different proper-OTO number can be interpreted to imply that non-trivial information in thermal correlators only appears at even-point functions. For instance, while there are 2-OTO 3-point functions, they lie in the KMS orbit of a 1-OTO correlator and their physics should thus not belong to the class of 2-OTO observables. The first place where we encounter non-trivial 2-OTO observables is in a 4-point function, where one can canonically choose it to be the chaos correlator [20,29,30]. The natural observable to pick, is the tremelo correlator, where every alternate operator is a past turning point operator. Equivalently the δ-list would comprise of n 2 alternating 1s, e.g., the last two roles of Table 3. These are distinguished by the fact that every correlator in their KMS orbit has maximal allowed OTO-number for the given number of operators.
Given that this observable is cognizant of the detailed dynamics of thermalization, one may wonder if higher even-point functions provide further detailed signatures. In particular, could one view the higher-point functions as higher moments of some distributions that captures equilibration? One useful avenue to examine is to explore how operator scrambling interplays with higher-OTOs. Similar considerations have inspired explorations of k-designs in chaotic quantum channels [22,34].
It is also worthwhile to explore in more detail the physical implication of 2-OTO thermal correlators which lie in the KMS orbit of a Schwinger-Keldysh correlator. Since a natural subset of the latter forms the basis of transport theory, it would be interesting to understand the OTO fluctuations in terms of more natural response coefficients. These considerations could provide useful generalizations of the Jarzynski relation [27], which we hope to explore elsewhere.
Any perturbation theoretic/path-integral approach to these correlators needs a way of effectively re-packaging the contour ordered correlators that arise from generalized Schwinger-Keldysh contours.
This would mean extending the existing results on Schwinger-Keldysh contour thermal correlators [5,6] whereby these correlators have been written as outer-products of certain 2-component column vectors thus bringing out their KMS structure. Like their Schwinger-Keldysh analogues, this form could be very useful for summering the non-linear fluctuation-dissipation relations, for establishing spectral representations and sum rules, for efficient hard thermal loop (HTL) approximations and for deriving generalized Kubo formulae.
Another interesting extension of the work here is to set up the analogue of Brownian motion which keeps track of higher OTO correlations of the bath. This is a natural higher OTO extension of the famous work by Feynman-Vernon [35] and Caldeira-Leggett [36].
Finally, while we have focused exclusively on thermal density matrices, it should be possible to generalize the above discussion to a general initial stateρ initial . Working with the modular Hamiltonian, K = − log(ρ initial ) we could effectively run the same arguments for the modular evolved operators (see also the discussion section of [19]). Sliding an operator through the density matrix as here would now result in the modular evolution of the operator. Various authors have noted the similarity between thermal and modular evolution, and in the holographic context, the modular evolution plays a role in construction of local bulk observables in the entanglement wedge [37][38][39]. It would be interesting to examine whether the modular KMS relations which one can derive by sliding operators through the density matrix have any useful information to impart for these considerations.

A Thermal Schwinger-Keldysh correlators
Motivated by developing a nonlinear generalization of the Fluctuation-Dissipation (FD) theorem for higher-point response functions, [7] used the Schwinger-Keldysh formalism as the appropriate framework for studying nonlinear response theory near thermal equilibrium. They worked with Schwinger-Keldysh contour correlators and used the KMS relations to derive the generalized FD relations for real time (1-OTO) thermal Green functions. The FD relations furnish a spectral representation for retarded and advanced thermal Green functions.
In this appendix we will compare the results of our paper with those of [7]. In the process, we will derive expressions for 3 and 4 point (1-OTO) thermal contour correlators in terms of our commutator basis. 21 The thermal n-point Green function in [7] is defined as G rrara.... ≡ (−i) n−1 2 nr−1 T SK 1 r 2 r 3 a 4 r 5 a .......
where i r stands for (A iR (t i ) + A iL (t i ))/2 (average basis), i a stands for A iR (t i ) − A iL (t i ) (difference basis), and T SK denotes Schwinger-Keldysh (i.e., 1-OTO) contour ordering. 22 Three-point spectral function relations: Consider 1-OTO 3-point contour correlators: We can use Keldysh rules to write the contour correlators as nested correlators: where we have introduced the time-ordering step-function: Θ ijk ≡ Θ(t i > t j > t k ). Passing onto our causal basis for nested correlators we can write expressions for these contour correlators in terms of nested commutators: This gives, for a thermal density matrix, an expression for 1-OTO 3-point contour correlators in terms of two spectral functions. This expression, in a different notation, was given in Eq. (31) of [5]. Wang-Heinz [7] write the FD relations for 3-point 1-OTO contour correlators in the Av-Dif basis as As above, we can evaluate all these contour correlators in our commutator basis (noting that complex conjugation reverses the time ordering in Theta functions), to find: (A.6) 21 Similar expressions are derived in [5,40]. 22 The factors of 2 are chosen so that they cancel on using the Keldysh rules to move to the nested commutator basis: nr is the number of average operators in the correlator . With all these expressions in hand we can compare both sides of (A.5) and confirm that they indeed hold.
Four-point spectral function relations: We now move to the 4-point function case. Wang-Heinz [7] write the FD relations for 4-point 1-OTO contour correlators in the Av-Dif basis in the following form: We use the Keldysh rules to get the following contour correlators in nested form: where the sums run over permutations of labels as indicated.
Using (3.15) (or more generally (3.26)) we can change all nestings to commutators and go to our causal basis: These equations (together with the Jacobi relations of the kind in (3.16)) give the SK correlators G rrrr , G rrra , G rraa , G raaa in our causal commutator basis [ [1σ(234)] ] . Note that for any fixed time ordering the above expressions simplify significantly. We can similarly get expressions for all the G's in (A.7) (the complex conjugation again reverses the time ordering in Theta functions) and explicitly check that it is true.

B Nested Thermal brackets and tJacobi relations
In this appendix we illustrate the use of the thermal nested brackets defined in (3.8) for deriving FD relations for nested correlators. We will consider nested thermal (anti-)commutators and present various tJacobi operator relations between them. These are analogous to the sJacobi relations discussed in [18] and in fact comprise a one-parameter deformation of those identities. An important tool for the following considerations is Lemma 4: Given the graded commutator (3.7) and the thermally deformed graded commutator (3.8), the following thermal involution identity holds: where 1 is the identity operator.
The terminology thermal involution refers to the fact that the right hand side of (B.1) is a nested thermally deformed correlator, where the nesting is inverted compared to the left hand side.
Proof: We will prove Lemma 4 inductively, illustrating the basic ideas and using the KMS relations in appropriate steps. Firstly, we have where we have used KMS relations in the first equality above to bring X to the left. These manipulations can be recursively applied to pass to higher orders. For instance replacing X by [X, C] ε C in the above equation we get: , we thus obtain and so on recursively. This process makes clear that the relations in (B.1) hold at arbitrary level.
Having established (B.1), let us do a counting of the number of correlators contained therein, and of the number of relations between those. There are 2 n n! thermal nested npoint structures of the form (B.5) These correlators should be expanded in a Wightman basis of n! structures. Thus there ought to be (2 n − 1)n! relations between them. We refer to these as tJacobi relations. Of these, (2 n−1 − 1)n! encode the information of (n + 1)-point sJacobi relations, which are valid for any density matrix [18]. The remaining (improper) ones are more directly related to a basic fluctuation-dissipation theorem. More precisely, we can understand all tJacobi relations according to the following two complementary observations.
• Improper tJacobi relations: These are the 2 n−1 n! tJacobi relations, which can be described as a simple nesting of the basic n = 1 "seed" of fluctuation-dissipation relations inside thermally deformed (anti-)commutators. The n = 1 FD seed is the trivial identity This relation captures the basic FD theorem, since we can multiply by B from the left and upon taking expectation value infer using (B.1) that the FD theorem (3.9) holds (see also below for illustration). Nesting this relation with (n − 1) further operators inside (n − 1) thermal brackets and considering all operator permutations gives 2 n−1 n! of the tJacobi relations.
• Relations between n-point tJacobis and (n + 1)-point sJacobis: To understand the remaining (2 n−1 − 1)n! tJacobi relations, we start with the observation that (2 n−1 − 1)n! is precisely the number of (n+1)-point sJacobi relations, where the (n+1) st operator is being given a fixed position (i.e., only considering permutations of the remaining n operators).
This is not coincidental: the thermal involution relation (B.1) allows us to relate any given n-point tJacobi relation to a standard sJacobi-type relation for (n + 1)-point functions. To see this, let us assume we have a tJacobi relation of the form σ,εσ X σ,εσ = 0 with X σ,εσ as in (B.5). We can now multiply with the (n + 1) st operator and find, using (B.1): The right hand side is of the form of an (n + 1)-point sJacobi relation.
Note that one the one hand it is justified to call the relations (B.7) as sJacobi identities: all thermal factors occurring on the left hand side have been removed on the right hand side and the relations thus do not rely on thermality! On the other hand, we should note that the right hand side of (B.7) is not written in the standard form of sJacobi relations explored in [18], since a particular operator A n+1 is singled out to be innermost in all terms. We conjecture that the mapping described above is nevertheless one-to-one.
Conjecture: Every n-point tJacobi relation is of one of two types: either it is one of the 2 n−1 n! improper ones (as described above), or it can be derived as descending from an (n+1)point sJacobi relation via (B.7) (now reading that equation from right to left).
Let us now illustrate these general constructions for small values of n. n = 1 tJacobi: This is the simplest relation since we have a single operator. From (B.1) we find: This relation captures the basic FD theorem, since we can multiply by B on the left and upon taking expectation value infer using (B.1) that (3.9) holds. One can similarly write down the tJacobi relations for higher values of n; explicit expressions for n = 3 tJacobi relations are given below. n = 3 tJacobi: The n = 3 tJacobi relations following from KMS relations are: The choice of signs give 2 n−1 n! = 24 of the tJacobi relations at n = 3. One can check explicitly that these relations are all linearly independent. (B.14) and their six permutations. We believe that the relations in (B.14) are linearly independent, and together with (B.12) and (B.13) give all the 42 tJacobi relations which whittle down the space of 4-point functions to 6.

C Nested correlator bases: derivations and proofs
In this appendix, we first derive a basis of n! nested correlators, which span the space of all n-point functions for any initial density matrix. In the second part, we prove that a canonical choice of (n − 1)! of these serve as a basis, if in addition we assume a thermal initial state with KMS condition (Lemma 2).

C.1 A nested basis for n-point functions (in generic states)
In the formulation of k-oto n-point correlators in terms of nested commutators and anticommutators, there are 2 n−2 n! basic correlators where S + n denotes the group of even permutations of n objects. We know that there are only n! independent Wightman functions, which means that the representation (C.1) is highly redundant. We will now describe a canonical choice of n! basis elements of the above form. All other nested correlators (and hence all Wightman functions) can then be expressed in terms of these using sJacobi identities described in [18]. We do not assume the KMS condition in this subsection.
We will construct this basis recursively as follows, to allow for ease of visualization of the process. At the end of the day we will prove that (3.22) is a suitable basis for our considerations in thermal field theories.
• n = 2: For 2-point function the basis is trivial to state: (C.2) • n = 3: For 3-point functions, we nest the 2-point function basis from above inside a commutator or anti-commutator and consider all 3 choices for the outermost operator:

(C.3)
Here we picked an anti-commutator for the outermost nesting only in the case where the operator with smallest index O 1 is outermost. For later convenience, we think of each row as a block labeled by the index of the outermost operator ∈ {1, 2, 3}. Each such block has exactly (n − 1)! elements induced by the 2-point function basis.
Then the remaining nested operators are given by the linear combinations 3 .

(C.4)
Let us now consider arbitrary values of n. Let B n−1 denote the basis of (n − 1)-point nested correlators. Then the basis B n is constructed as follows. We simply nest the (n − 1)! objects of the B n−1 basis inside a commutator or anti-commutator and consider all n possibilities for the outermost operator (which then accounts for n × (n − 1)! = n! choices).
We can choose to always pick a commutator for this nesting, except for one case, say, when the outermost operator is O 1 whence we pick an anti-commutator. 23 Explicitly, let us denote the n! elements of B n as {b  Proof: To prove that the construction (C.5) gives a basis of n! n-point functions, we use induction. For small values of n, we have given explicit constructions in the main text, so it remains to show that assuming B n−1 forms a basis of (n − 1)-point functions, then also B n forms a basis of n-point functions. We show this by explicitly constructing any given element of the Wightman basis in terms of B n .

(C.7)
It is straightforward to check this equation. We will now argue that every line of this equation is a combination of basis elements of B n . This involves two steps. First, it follows by induction that any string of n−1 operators appearing as the first entry of any of the (anti-)commutators above is a combination of the basis elements appearing in the basis B n−1 . This requires a different labeling of operators for each line, but the important part is that this can be done, simply because the first entry of each (anti-)commutator is an (n − 1)-point Wightman function. Secondly, we realize that these (relabeled permutations of) the basis elements of 23 There is some freedom here. Another possible choice would be to choose all the outermost brackets to be anti-commutators except for the case when the outermost operator is O1 and n is odd. Various other prescriptions are possible, but we choose a particularly simple one here.
B n−1 are precisely the objects that we nest inside (anti-)commutators in order to construct B n . Further, according to the construction (C.5), we always nest inside commutators, except when the outermost operator is 1, in which case we use an anti-commutator. This is exactly the same structure as that appearing in the linear combination (C.7). More explicitly, the first line of (C.7) is a linear combination of the (σ n )-th block of B n , the second line is a linear combination of the elements of the (σ n−1 )-th block, and so on. This completes the proof. As a corollary, we can immediately see from the explicit construction (C.7) that the linear combination that expresses any Wightman correlator in terms of the basis B n contains every single one of the n! elements of B n and all coefficients have the same absolute value. That is, we can write any Wightman n-point function as where σ ∈ S n and s i = ±1. In this sense the basis B n is very democratic: no particular Wightman function is any simpler or any more complicated than any other one, when expressed in this basis.

C.2 A causal basis for thermal correlators
In the previous subsection we constructed a basis B n of n! nested correlators, which form a basis of n-point functions in generic states. Let us now assume that the initial state is thermal, so we can use the KMS condition to further reduce the basis to (n − 1)! elements. One can easily see from the previous subsection that B n contains precisely (n−1)! nested commutators (i.e., nested correlators which involve no anti-commutators). These are all of the form [ [1σ(2) · · · σ(n)] ], i.e., the operator 1 is innermost and all others occur in all possible orders. We need to prove that all remaining elements of B n can be expressed in terms of the "causal basis" [ [1σ(2) · · · σ(n)] ] using the KMS condition. Alternatively, we can proceed by showing that every Wightman function can be expressed in terms of the causal basis. Since the KMS condition acts cyclically on Wightman functions (as in (2.13)), we can restrict to Wightman functions of the form 1σ(2) · · · σ(n) and show that they are in one-to-one correspondence with the causal correlators [ [1σ(2) · · · σ(n)] ] .
We will proceed in two steps: first, we give a general expansion for the causal basis correlators in terms of Wightman functions. While it is obvious that this exists, it will be useful to write it out more explicitly. Secondly, we will show the converse: every Wightman function of the form 1σ(2) · · · σ(n) can be written in terms of [ [1σ(2) · · · σ(n)] ] . Finally, we will also prove Lemma 3 to express generic nested correlators in terms of the causal basis, thus completing the arc reaching from the nested correlators to the basis B n , and its reduction to the causal basis in thermal states.
1. Causal basis in terms of Wightman correlators: Any element of the causal basis can be expressed in terms of Wightman functions of the form 1σ 2 · · · σ n (where σ i ≡ σ(i)).

Wightman correlators in terms of causal basis (proof of Lemma 2):
We now need to show that the relation (C.9) can be inverted to express a given Wightman function 1σ 2 · · · σ n in terms of causal nested correlators [ [1ρ 2 · · · ρ n ] ] . The result reads 1 σ 2 · · · σ n = ρ∈S n−1 (1 + f 1 ) (C.12) which was given for σ = id in Lemma 2 (we refer to §3.4 for more explanation). In the following, we also restrict to σ = id to simplify notation.
Proof: To prove (C.12), we need to show that it inverts (C.9). To this end, let us define the following thermal factor associated with two permutations ρ and σ of {2, . . . , n}: A slight variation of this object was encountered in Lemma 3. We claim that T (ρ,σ) {−} = δ ρ,id , i.e., it is unity if ρ = id, and vanishes otherwise. We will prove this statement below. First, however, let us note that (C.12) immediately follows once we show this. To see this, let us plug the equation T where we redefined ρ → π{s i } • ρ in the last step. The statement (C.12) now follows immediately from comparison with the already established formula (C.9).
(C. 19) This shows that there is a vanishing factor in (C.17), so T  with thermal factors given in (3.27). We start with an immediate generalization of (C.9) to the case where we allow for both commutators and anti-commutators (and we fix σ = id in (C.9) to declutter notation): [ [1 2 ε 2 · · · n εn ] ] = (1 + ε n e βωn ) where we exchanged the order of summation in the second step by relabelling ρ → π −1 {s i } • ρ. Formula (3.27) for the thermal factors now follows from the observation (C.16).

D Details on the harmonic oscillator
We provide here some details on the harmonic oscillator example that we skipped in §4.