Emergent Supersymmetry in Local Equilibrium Systems

Many physical processes we observe in nature involve variations of macroscopic quantities over spatial and temporal scales much larger than microscopic molecular collision scales and can be considered as in local thermal equilibrium. In this paper we show that any classical statistical system in local thermal equilibrium has an emergent supersymmetry at low energies. We use the framework of non-equilibrium effective field theory for quantum many-body systems defined on a closed time path contour and consider its classical limit. Unitarity of time evolution requires introducing anti-commuting degrees of freedom and BRST symmetry which survive in the classical limit. The local equilibrium is realized through a $Z_2$ dynamical KMS symmetry. We show that supersymmetry is equivalent to the combination of BRST and a specific consequence of the dynamical KMS symmetry, to which we refer as the special dynamical KMS condition. In particular, we prove a theorem stating that a system satisfying the special dynamical KMS condition is always supersymmetrizable. We discuss a number of examples explicitly, including model A for dynamical critical phenomena, a hydrodynamic theory of nonlinear diffusion, and fluctuating hydrodynamics for relativistic charged fluids.


Introduction
The goal of many-body physics is to explain and predict macroscopic phenomena. Except for some very simple systems, however, it is rarely possible to compute macroscopic behavior of a system directly from its microscopic description. For static properties of equilibrium systems, we have the extremely successful Laudau-Ginsburg-Wilson paradigm, which provides an effective field theory (EFT) description of long-distance (IR) physics (1.1)

JHEP01(2018)040
In (1.1), φ denotes collectively external sources, χ denotes collectively gapless modes, and S eff is the low energy effective action of the gapless modes obtained by integrating out gapped degrees of freedom. While in practice such direct integrations are almost always impossible, one can deduce the general form of S eff on physical ground. Two elements are needed for this purpose: (i) choice of the IR dynamical variables χ which best capture gapless (collective) degrees of freedom; (ii) symmetries of S eff . One can then write down S eff as the most general local field theory consistent with the symmetries. For a non-equilibrium system or dynamical quantities of an equilibrium system, partition function is inadequate. A large class of non-equilibrium observables can be extracted from the generating functional defined on a closed time path (CTP) [1][2][3][4][5][6] where ρ 0 denotes the state (density matrix) of the system, and U (t 2 , t 1 ; φ) is the evolution operator of the system from t 1 to t 2 in the presence of external sources denoted by φ. The sources are taken to be slowly varying functions and there are two copies of them, one for each leg of the CTP contour. The second line (1.3) is the "microscopic" path integral description, with ψ 1,2 denoting microscopic dynamical variables for the two copies of spacetime of the CTP and S 0 [ψ; φ] the microscopic action. Whereas in (1.1) derivatives with respect to sources φ give thermodynamic quantities of a system, in (1.2) derivatives with respect to φ 1,2 give dynamical properties of a (non)equilibrium system such as (nonlinear) response and fluctuating functions. As in (1.1) we can consider integrating out short-lived degrees of freedom in (1.3) to obtain a non-equilibrium EFT I eff for slow modes (denoting them collectively by χ 1,2 and there are now two copies of them) e W [φ 1 ,φ 2 ] = Dχ 1 Dχ 2 e iI eff [χ 1 ,φ 1 ;χ 2 ,φ 2 ;ρ 0 ] . (1.4) Again to write down the general form of I eff one needs to specify appropriate dynamical variables χ 1,2 and the symmetries satisfied by the low energy effective action I eff , although these tasks normally become significantly more challenging in non-equilibrium situations.
In (1.4) ρ 0 is also encoded in the couplings of I EFT (below for notational simplicity we will suppress ρ 0 in I eff ). In general I eff does not have the factorized form of (1.3), and is complex. It is often convenient to introduce the so-called r − a variables [7,8] where as usual χ r correspond to physical quantities while χ a can be interpreted as noises.
The functional integral (1.4) defines a "bare" theory at some short distance (time) cutoff scale. 1 Physics at larger distance and time scales is obtained by further applying JHEP01(2018)040 renormalization group procedure. While I eff in principle contains an infinite number of terms with increasingly higher number of derivatives, in practice to describe macroscopic phenomena one only needs to keep track of a finite number of relevant interactions.
Non-equilibrium EFTs provide powerful tools for dealing with dynamical questions and non-equilibrium systems. The effective action I EFT incorporates dissipations and retardation effects from the bath of short-lived degrees of freedom (which have been integrated out) in a medium. Its general structure has recently been used to derive from first principle the local second law of thermodynamics [9], and a new formulation of fluctuating hydrodynamics has been proposed in terms of such an EFT [10,11] (see also [12][13][14][15][16][17]). See also [18] for a review of applications to driven open systems. When an I EFT is truncated to quadratic order in noises (i.e. a-variables) the path integral (1.4) reduces to a so-called Martin-Siggia-Rose-De Dominicis-Janssen [19][20][21] functional integral which is in turn equivalent to a stochastic Langevin equations (for a review see [22]).
Compared to EFTs for equilibrium systems, there are new elements in identifying both dynamical variables and symmetries for a non-equilibrium EFT (1.4). The unitarity of time evolution in (1.2) implies that the action should in addition satisfy the following conditions (see e.g. [9,10] for more details) Im I eff ≥ 0 (1.7) I eff [χ, φ; χ, φ] = 0, or I eff [χ r = χ, φ r = φ; χ a = 0, φ a = 0] = 0 , (1.8) where for definiteness we have taken χ 1,2 and sources φ 1,2 to be real. These conditions are, however, enough only for performing the functional integrals of (1.4) at tree level. With loops included one also has to worry about defining the integration measure Dχ 1 Dχ 2 precisely. To see this, in (1.2) taking φ 1 = φ 2 = φ, we then find that Tr U (+∞, −∞; φ)ρ 0 U † (+∞, −∞; φ) = Tr(ρ 0 ) = 1 ⇒ W [φ, φ] = 0 . (1.9) While equation (1.8) leads to (1.9) at tree-level, this is no longer so when including loops and one has to include an additional integration measure factor. See e.g. section I E of [10] for an explicit discussion. To ensure (1.9) at loop level can use the standard trick of parameterizing integration measures by introducing an anti-commuting partner for each bosonic variable, i.e. c r , c a for χ r , χ a respectively, and requiring the action to be invariant under the following BRST-type fermionic transformation [23] δχ r = ǫc r , δc a = ǫχ a . (1.10) Here ǫ is an anti-commuting constant. To show that (1.10) is enough to ensure (1.9) at loop level is quite simple and is reproduced in appendix A for completeness. In particular, the BRST invariance automatically leads to (1.8) for the bosonic part of the action. c r,a are anti-commuting but transform the same as their bosonic partners under spacetime rotations. They will be subsequently referred to as ghost variables following standard JHEP01(2018)040 terminology. The need for ghosts and BRST symmetry can also be anticipated from results on the functional integral forms of stochastic equations [23], and has been emphasized recently [10,16] in the context of fluctuating hydrodynamics. 2 There are three different regimes for (1.4). The first is the full quantum regime where path integrations describe both quantum and classical statistical fluctuations. The second is the classical regime with → 0. In the → 0 limit the path integrals survive and describe classical statistical fluctuations. The third is the level of equations of motion which corresponds to the thermodynamic limit with all classical and quantum fluctuations neglected. Since the constraints (1.6)-(1.8) concern only with the general structure of the action, they remain in the classical limit. Similarly the requirement W [φ, φ] = 0 also survives the classical limit, and so do ghost variables and the corresponding BRST symmetry. It is striking that a classical statistical system is significantly constrained by these remnants from quantum unitarity.
For many physical processes in nature, macroscopic physical quantities of interests typically vary over spatial and temporal scales much larger than microscopic molecular collision scales (or any microscopic interaction scales). Such a system is considered as in local equilibrium, for which an additional Z 2 symmetry should be imposed on I eff [9,11].
A subclass of local equilibrium systems correspond to thermal systems perturbed by slowly varying external sources, and in this case the need for this Z 2 symmetry can be readily understood as follows. For ρ 0 = 1 Z e −β 0 H , the generating functional (1.2) satisfies an additional constraint coming from combining the Kubo-Martin-Schwinger (KMS) condition [25][26][27] with time reversal invariance, for arbitrary θ ∈ [0, β 0 ]. Below we will simply refer to (1.11) as the KMS condition, but it should be kept in mind it also encodes consequences of microscopic time-reversal symmetry. 3 Additional condition(s) then need to be imposed on I eff for (1.4) to satisfy (1.11). For variables χ r,a associated with non-conserved quantities, the required symmetry is well known, probably since 70's [21,28,29]. In this case the couplings between χ r,a and external sources are the standard ones See also [24]. Refs. [16,24] appear to require two BRST generators while [23] and [10] require only one. 3 The KMS condition itself only relates W to a time-reversed one, and so does the time reversal symmetry of the microscopic theory. Only the combination of them leads to a nontrivial constraint on W itself [7,8]. See [10] for a detailed discussion. Depending on circumstances one could combine KMS with T or PT or CPT . For definiteness here we follow [10] to combine it with PT . It is simple to adapt (1.12) for a system with only T invariance by simply removing the minus signs before x.

JHEP01(2018)040
which then immediately implies that for (1.4) to satisfy (1.11), the action should satisfy 4 (1.14) Following [9,11] we will refer to (1.14) as the dynamical KMS condition. In the absence of external sources it becomes a Z 2 symmetry (dynamical KMS symmetry) of the action. We will refer to transformations (1.15) on dynamical variables as dynamical KMS transformations.
The story for variables associated with conserved quantities (hydrodynamical variables) is more complicated since the couplings to external sources are more intricate making it more difficult to deduce the needed transformations on dynamical variables. In [10] a shortcut was proposed which imposes (1.11) on a contact-term action which in turn constrains the action for dynamical variables through the special structure of the couplings between dynamical variables and external sources. It was termed as the local KMS condition. Only very recently were dynamical KMS transformations on hydrodynamical variables finally found in [11]. We emphasize that a system in local equilibrium is not restricted to a thermal density matrix in the presence of slowly varying external sources. Such a system can, for example, be in a pure state. In these more general cases while there is no such requirement as (1.14), invariance of an action under dynamical KMS transformations ensures the system is in a local equilibrium. 5 For example, in the classical limit dynamical KMS condition ensures that first law and second law of thermodynamics, as well as fluctuation-dissipation and Onsager relations are all satisfied locally [9,11].
With dynamical KMS transformations for bosonic variables understood, in this paper we consider the extensions to ghost variables, which are needed to have a complete formulation of a non-equilibrium EFT. For example, to ensure (1.11) at loop level we need dynamical KMS transformations on all variables.
Furthermore, it has been long known in the context of functional representation for linear stochastic systems that there is an emergent supersymmetry as a consequence of fluctuation-dissipation relations [23,[30][31][32][33]. More recently, it was found in [10] that after imposing the local KMS condition and BRST symmetry there is also an emergent supersymmetry for a hydrodynamic theory of nonlinear diffusion. We would like to understand the precise origin and the full extent of this emergent supersymmetry. In particular, we would like to extend the discussion to a general non-equilibrium EFT including full fluctuating hydrodynamics. 6 We will restrict our discussion to the classical level with → 0. At the classical level, the dynamical KMS transformations dramatically simplify. For example, equations (1.12)

JHEP01(2018)040
and (1.15) becomẽ (1.17) which are local transformations combined with a spacetime reflection. We stress that these are finite Z 2 transformations. The dynamical KMS transformations for hydrodynamical variable, although more involved, have a similar structure (see section 5). The quantum regime has a number of additional complications and will not be pursued here (see section 6 for a brief discussion). We will show that any system in local equilibrium has an emergent supersymmetry at low energies. With increasing complications and generality we consider three classes of systems depending on whether or not a system has conserved quantities or dynamical temperature: (i) no conserved quantities with a fixed background temperature; (ii) with conserved quantities and a fixed background temperature; (iii) with conserved quantities and dynamical temperature. Clearly the third class includes all systems. As an example for class (i) we consider model A of critical dynamics [36], for class (ii) a theory of nonlinear diffusion, and for class (iii) a fluctuating hydrodynamics for charged fluids proposed in [10,11]. It turns out when expressed in terms of the right sets of variables, all three classes have essentially the same structure. Here is a summary of the main results: 1. We show that there is essentially a unique extension of dynamical KMS transformations to ghost variables which is self-consistent. The dynamical KMS transformation on ghost variables turn out to be a Z 4 operation, but is still a Z 2 operation of the action.
2. For any action the combination of BRST symmetry and dynamical KMS symmetry leads to an emergent fermionic symmetry which together with the BRST symmetry forms a supersymmetric algebra.
3. Starting with a supersymmetric action one can always construct an action which is both BRST and dynamical KMS invariant.
4. Supersymmetry does not impose the full dynamical KMS invariance, only a particular consequence of the dynamical KMS symmetry, to which we refer as the special dynamical KMS condition. Conversely we prove a theorem stating that any bosonic action satisfying the special dynamical KMS condition is always supersymmetrizable.
5. For a system for which temperature is non-dynamical (i.e. with a fixed constant temperature), one finds a global supersymmetry. For a system for which temperature is dynamical, such as a fluctuating hydrodynamics, one finds a local supersymmetry.
6. Supplementing the bosonic story of fluctuating hydrodynamics proposed in [10,11] with dynamics of ghosts, this paper finally gives a complete formulation of fluctuating hydrodynamics in the classical regime.

JHEP01(2018)040
The plan of the paper is as follows. In next section we present a general discussion of emergence of supersymmetry from BRST and dynamical KMS symmetries. In section 3-5 we discuss three classes of examples. We conclude in section 6 with a discussion of future directions. In appendix A we give further argument for the need of BRST symmetry. Appendix B contains details of a proof for a supersymmetrizability theorem.
While this paper is in preparation we learned that overlapping results have been obtained by Kristan Jensen, Natalia Pinzani-Fokeeva, and Amos Yarom [35].

Emergent supersymmetry: general structure
In this section we present a general discussion of emergence of supersymmetry from BRST and dynamical KMS symmetries.

General case
Consider an action I[F i ] with F i = (b i , f i ) which b i denotes collectively bosonic source and dynamical fields, and f i denotes collectively anti-commuting source and dynamical fields (ghost variables here). To make our equations compact we will use index i to denote both field species and spacetime points. We assume that the action is invariant under a BRST-type fermionic symmetry, i.e.
where ǫ is an anti-commuting constant, and Q is an anti-commuting operator satisfying Now let us suppose that I is invariant under another bosonic symmetry where K α is an invertible bosonic operator (i.e. maps bosons to bosons and ghosts to ghosts) and index α denotes different elements of the symmetry group. Note that while Q acts as a derivation, K α acts as a finite transformation. Acting on a product, K α transforms all factors at the same time.
Clearly the action is also invariant under the combined operations Q α = K α QK −1 α , where the notation on the right hand side means after evaluating (QF j ) ∂K −1

JHEP01(2018)040
Thus we find that for each symmetry transformation K α there is an emergent fermionic symmetry Q α . Note that the collection {Q α } also includes the original Q as {K α } includes the identity element. Note that Suppose we have an action I 0 [F i ] which is not invariant under a K α -transformation. Then it follows immediately from our definition that (2.9)

A special case
Now let us specialize to a situation which will be relevant for the rest of this paper, with {K α } = 1, K, K 2 , K −1 being a set of Z 4 transformations satisfying In this case for any action I 0 (which is not necessarily invariant under K) we have as an action is always even in the number of ghosts variables. Also note that Q K 2 = −Q and Q K −1 = −Q K , and thus the independent {Q α }'s are Q andQ ≡ Q K . From (2.9) and (2.11) we then have for any action I 0 Now suppose I 0 is BRST invariant, i.e. QI 0 = 0. We can construct a K-invariant action as But this action is in general not BRST invariant as QĨ 0 does not have to be zero. From (2.12) we conclude that for I to be both BRST and K-invariant, the sufficient and necessary condition is that I 0 should in addition be invariant underQ.

Strategy for extending dynamical KMS transformations to ghosts
Since ghost variables are introduced to give the correct integration measure and do not directly couple to external sources, there is no obvious principle to determine how they should transform under dynamical KMS symmetry. Our strategy is based on the following non-trivial self-consistency requirement: BRST and dynamical KMS invariance of the

JHEP01(2018)040
full action does not put further constraint on the pure bosonic part of the action. More explicitly, with the full action written in a form then the pure bosonic part I b [b i ] should coincide with the most general action one can construct based (1.6)-(1.8) and the bosonic dynamical KMS invariance. This requirement is due to that the bosonic action I b already provides a complete formulation of tree-level physics, thus extension of dynamical KMS symmetry to the ghost sector should not change that physics. The requirement is highly nontrivial mathematically as dynamical KMS invariance constrains the ghost part of the action I f which in turn constrains the bosonic part I b via BRST symmetry. Our discussion contains the following elements: 1. Applying the consistency requirement at quadratic level in dynamical variables uniquely determines the dynamical KMS transformation for ghost variables at linear level.
2. As a simplest possibility we postulate the linear transformation deduced from the quadratic action is the full transformation. Including both bosons and ghosts, the dynamical KMS transformations have the structure discussed around (2.10). We then constructQ explicitly from Q and K using (2.4). One finds that Q,Q form a supersymmetric algebra. In other words, a non-equilibrium EFT must be supersymmetric invariant.
3. We provide a strong support for the postulate of item 2 by proving that the selfconsistency requirement is indeed satisfied for the full nonlinear action.
From the above discussion and that around (2.14) we conclude that one can obtain a BRST and dynamical KMS invariant action by first writing down a most general supersymmetric action and then impose (2.14).

Special dynamical KMS symmetry and a theorem on supersymmtrizability
In this subsection we elaborate a bit further on the self-consistency requirement of the previous subsection. Consider a most general action I b [b i ] of bosonic variables which satisfies (1.6)-(1.8) and the dynamical KMS condition (1.14). Now adding a ghost partner f i for each bosonic variable to obtain a full action I[b i , f i ] which is BRST and dynamical KMS invariant. From our discussion above we learned that this full action must be supersymmetric. The self-consistency requirement requires that the bosonic part of I[b i , f i ] should coincide with the original I b . This in turn requires that I b be supersymmetrizable. Conversely, if I b is supersymmetrizable, then we can construct a BRST and dynamical KMS invariant action with the same bosonic part by first constructing a supersymmetric extension of I b and then using (2.14). Thus for I b to be supersymmetriable is both sufficient and necessary

JHEP01(2018)040
for constructing a full BRST and dynamical KMS invariant action. So the self-consistent condition boils down to the statement: a bosonic action which satisfies (1.6)-(1.8) and the bosonic dynamical KMS condition (1.14) should be supersymmetrizable under the supersymmetry generated by Q andQ.
We will be able to prove that this is indeed the case. In fact we will be able to prove a stronger statement which was first observed in [10] for a theory of nonlinear diffusion at cubic level. To describe the statement, we need to be a bit more specific on the general structure of dynamical KMS condition. Equations (1.6) and (1.8) imply that the bosonic Lagrangian density L b can be expanded in a-fields as where we use Λ r , Φ a to denote collectively r-and a-fields respectively. Note that the sum starts with n = 1 as n = 0 term is not allowed by (1.8). In the classical limit → 0, the dynamical KMS transformation on bosonic variables is a Z 2 transformation which can be schematically written as 7 where Φ r is a product of bosonic r-variables with altogether one derivative. The dynamical KMS condition (1.14) can then be written as whereL b is obtained by plugging (2.17) into (2.16) and taking x → −x, i.e.
is obtained from f (n) by flipping the signs of all derivatives and L b k denotes terms inL b with k factors of Φ a . Note that the k-sum starts with zero. Equating equation (2.18) order by order in the expansion of Φ a we then find an infinite number of conditions and where V µ k denotes terms containing k factors of Φ a . Alternatively we can also impose the dynamical KMS condition as follows. Take a Lagrangian density L 0 of the form (2.16). Due to Z 2 nature of the transformation, then  (2.20) as the special dynamical KMS condition. From (2.12)-(2.14) and that BRST symmetry implies (1.8) for the bosonic part, it then follows that supersymmetry of L 0 ensures the special dynamical KMS condition for its bosonic part. Conversely, we will be able to prove the following supersymmetrizability theorem: Any local bosonic Lagrangian which satisfies (1.8) and the special dynamical KMS condition (2.20) is supersymmetrizable.
Comparing with the discussion around (2.14), we see the procedure of imposing (2.22) commutes with supersymmetrization. One could either do it before or after.
3 With no conserved quantities at a fixed temperature: model A In this and next two sections we consider the extension of dynamical KMS transformations to ghosts and the associated supersymmetry for some explicit examples of non-equilibrium EFTs. In this section we will consider systems with no conserved quantities at a fixed background temperature, i.e. temperature is not a dynamical variable. In this case the story is technically much simpler, but captures all the essential elements.
As an illustration of a system with no conservation laws, we consider the critical dynamics of a n-component real order parameter χ i , i = 1, · · · , N at a fixed inverse temperature β 0 (i.e. model A [36,37]). The dynamical variables in (1.4) are then {χ ri , χ ai } and the action should be invariant under an SO(N ) symmetry which rotates χ ri , χ ai simultaneously. 8 The dynamical KMS transformation for bosonic variables is the same as (1.17) In this case the couplings (1.13) to external sources are rather trivial, so we will suppress the sources below. At quadratic order in χ a,r (but to all orders in derivatives) the bosonic part of the Lagrangian can be written as where G ra are G aa are some differential operators. Note that by definition G aa satisfies G aa = G * aa where G * aa denotes the operator obtained from G aa by taking all ∂ µ to −∂ µ . Imposing the dynamical KMS condition leads to the condition As discussed in the Introduction we should also introduce anti-commuting partners c ri , c ai for χ ri , χ ai respectively, and require the action to be invariant under the following BRST transformations δχ ri ≡ ǫQχ ri = ǫc ri , δc ai ≡ ǫQc ai = ǫχ ai , Qχ ai = Qc ri = 0 . (3.4)

JHEP01(2018)040
At quadratic level the most general Lagrangian invariant under (3.4) can be written as 9 where G rr is an arbitrary differential operator satisfying G rr = −G * rr from anti-commuting nature of c ri . At quadratic level the dynamical transformation on c r,a must be linear and requiring no further constraints on G ra we find that the only possibility is to require Lagrangian be invariant under 10 which in turn requires G rr = 0. We thus propose (3.6) as the dynamical KMS transformation for ghosts. Combining (3.1) and (3.6) we find the Z 4 structure discussed around equation (2.10) with Now applying (2.4) to (3.7)-(3.8) and (3.4), we find that It can also be readily checked that i.e. Q,Q form a supersymmetric algebra. It can also be checked explicitly that invariance of (3.5) under Q andQ indeed leads to (3.3).

Superspace
To impose supersymmetry, it is convenient to use superspace formalism [38]. We introduce two Grassmannian coordinates θ,θ and the superfield Q,Q can then be written in terms of the following differential operators with (3.4) and (3.9) given by

JHEP01(2018)040
Note that as usual acting on superfields with an opposite sign from (3.10). The corresponding covariant derivatives arē which satisfy Note that A general Lagrangian which is invariant under (3.4) and (3.9) can then be written as where F is a local expression constructed out of Ψ i , and their covariant and ordinary derivatives.

Proof of the supersymmetrizability theorem
We now present a proof of the supersymmetrizability theorem stated at the end of section 2.4. Here we will discuss the main steps. There is a key step whose proof is rather contrived, which we will leave to appendix B. Consider a general bosonic action I b [χ ri , χ ai ], which satisfies (1.6), (1.8), and the special dynamical KMS condition (2.20). Since the dynamical KMS transformation (3.1) is linear in fields, it does not change the total number of fields in a given term. In other words, suppose we expand L b in terms of the power of dynamical variables where L n contains altogether n factors of χ r,a , then different L n 's do not mix under dynamical KMS transformations. It is then enough to prove the theorem for a general L n . L n can be written schematically in a form where each term should be understood as

JHEP01(2018)040
and the indices I, J include both species indices and indices for all possible derivatives on them. It is easy to write them in momentum space, for example, The second term resulted from the replacement can be regrouped into terms with m − 1 χ a 's. Continuing this procedure we will then be left with terms with one factor χ a which we will denote as where g is symmetric in J 2 , · · · J n indices.
Note that under a dynamical KMS transformation, a term of the form (3.24) will always contain at least one factor of χ a due to theχ a factor there. Thus the special dynamical KMS condition (2.20) will only involve (3.26) which can be written as and in momentum space Recall that index J 1 · · · include both species indices and momenta, i.e. χ J k r ≡ χ ri (ω k , k k ) with k = 1, · · · n, and momentum conservation implies that n k=1 ω k = 0. From properties of symmetric polynomials one can show that (3.27) implies that I b can be written as (the proof of which is a bit involved and we leave it to appendix B) can be written in a form Using (3.17)-(3.18) such a term can be supersymmetrized as We thus have shown all terms in L n can be supersymmetrized, which concludes the proof.
To conclude this subsection let us note that the proof does not depend on the nature of the species index i, which can be generalized to any kinds of indices including spacetime indices, say for a tensor field. In particular we will see the proof applies also to the examples of next two sections.

Full formulation
To complete the formulation of the model A EFT, now let us consider the generalization of (1.6) to the full action. A natural generalization is with η r,a = ±1. It can be readily checked that only the choice η r = 1 and η a = −1 is compatible with the BRST symmetry, and thus we should have To see that (3.34) is also compatible withQ and K operations, let us define an operation S asŜ where F r,a denote respectively any r and a-type variables (including sources, bosonic and ghost dynamical variables). Equation and thus supersymmetry is preserved byŜ. Also note that [K,Ŝ] = 0 acting on bosonic fields and {K,Ŝ} = 0 acting on ghost fields, and thus acting on action K commutes with S due to the fact that an action must contain even number of ghost fields. This shows (3.34) is also compatible with dynamical KMS condition.
We also need to check the self-consistency: (3.34) should not put further constraints on the bosonic action. This amounts to showing the full action I obtained by supersymmetrizing a bosonic action satisfying (1.6) satisfiesŜI = −I. Note that from the discussion of last section any term in the bosonic action can be supersymmetrized to a single term in terms JHEP01(2018)040 of superfields. We thus only need to show that any term in the superspace has a definite eigenvalue underŜ (then this eigenvalue must agree with that of the bosonic part). To see this note that the superfield (3.11) has the following structure under transformation ofŜ: Ψ ∼ X + iθX,DΨ ∼ iX, DΨ ∼ X + iθX,DDΨ ∼ iX and DDΨ ∼ iX +θX where X denotes the type of fields with eigenvalue 1 underŜ. Note that X does not containθ. Thus any term consisting of products of such superfields will be the proportional to X n (X +iθX) m ∼ iθX n+m , where only oneθ survives. It has a definite eigenvalue ofŜ. Finally given that [K,Ŝ] = 0 acting on action, the step (2.14) does not change the eigenvalue ofŜ.
We can now present the full procedure for constructing the EFT for model A using supersymmetry: 1. Construct a most general supersymmetric action I 0 , which satisfies (3.34) (and of course whatever other symmetries of the system).
3. The bosonic part of the action should further be constrained by (1.7).
Instead of using supersymmetry one can of course directly impose BRST symmetry and the special dynamical KMS condition. With the powerful formalism of superspace, supersymmetry should in general be a faster route. Finally we note that in a most general supersymmetric action there can be terms which are not related to the pure bosonic action, i.e. terms involving ghosts transform among themselves under supersymmetric transformations. Whether one should include such terms requires further consideration.

With conserved quantities at a fixed temperature: nonlinear diffusion
In this section as an example of systems with conserved quantities at a fixed temperature we consider the hydrodynamic theory for nonlinear diffusion developed in [10]. For slow variables associated with conserved quantities, couplings to external sources play an important role in the formulation of the theory. So in this section we will turn on external sources from the beginning. We will see that the same structure as that of model A emerges. The discussion here generalizes and systemizes some previous observations in [10] regarding BRST invariance, KMS conditions and supersymmetry.
We consider the theory of diffusion mode associated with a U(1) conserved current at a fixed inverse temperature β 0 , ignoring possible couplings between the diffusion mode and other hydrodynamical modes. The dynamical variables are ϕ r , ϕ a with ϕ r interpreted as the diffusion mode and ϕ a the corresponding noise variable. The background sources are A rµ and A aµ which couple to conserved currents J µ a and J µ r respectively. The bosonic action satisfy the following conditions: 1. ϕ r , ϕ a must always be acted on by at least one derivatives. We will thus count ∂ µ ϕ r,a as having zeroth derivative. In the presence of background fields A rµ , A aµ , the action should depend only on the combinations The local chemical potential is given by µ = B r0 = A r0 + µ d with µ d = ∂ 0 ϕ r giving the dynamical part, and it is often convenient to useμ = β 0 µ andμ d = β 0 µ d .

The action is invariant under
The dynamical KMS transformation on bosonic variables arẽ and when including background fields We now introduce ghost partners c r,a for ϕ r,a respectively, and require the action (in the absence of background fields) to be invariant under transformation δϕ r = ǫc r , δc a = ǫϕ a . (4.6) In the presence of external sources it is convenient to introduce ghost partners η rµ , η aµ for A rµ , A aµ respectively and the action should be now be invariant under the combinations of (4.6) and 11 δA rµ = ǫη rµ , δη aµ = ǫA aµ . Introducing H rµ = η rµ + ∂ µ c r , H aµ = η aµ + ∂ µ c a . Extending the dynamical KMS transformation (4.5) to ghost fields proceeds in an identical manner as the example of section 3 and we find or in terms of source and dynamical fields separatelỹ (4.11) Similarly one findsQ transformation is given bȳ δB rµ = −ǭH aµ ,δB aµ =ǭiβ 0 ∂ 0 H aµ ,δH rµ =ǭ(B aµ + iβ 0 ∂ 0 B rµ ),δH aµ = 0 . (4.12) 11 See appendix A for motivation for introducing ghost partners and BRST transformations for external sources.

JHEP01(2018)040
Note that the above transformations for B rµ , B aµ , H rµ , H aµ are identical to those of section 3 for χ ri , χ ai , c ri , c ai . Thus all the results there can be directly carried over with simple change of notations. For example, the superfield now has the form Σ µ (x, θ,θ) = B rµ + θH rµ + H aµθ + θθB aµ . (4.13) The proof of the supersymmetrizability theorem also carries over as in the proof the nature of species indice i, j did not play any role. Here they are replaced by µ, ν of vector indices.

BRST and emergent supersymmetry for fluctuating hydrodynamics
As an example with both conserved quantities and dynamical temperature, in this section we consider the full fluctuating hydrodynamics for a relativistic charged fluid in the classical limit [10,11], which is the low energy effective theory for slow modes associated with stress tensor and a conserved U(1) current.
Here the story is much more complicated than those of the previous two examples. Remarkably, we will see in the end an almost identical structure to that of previous examples emerges when the theory is expressed in terms of an appropriate set of variables. There is also an important difference. In previous two examples with a fixed temperature we saw that the background temperature plays an important role in the supersymmetric algebra (3.10). Now with a dynamical temperature we will see that supersymmetry becomes local.
Supplementing the bosonic story of [10,11] with the ghost sector, the discussion here completes the formulation of fluctuating hydrodynamics in the classical regime.

Bosonic sector
The dynamical variables are given by χ r = (X µ (σ), ϕ r (σ), β(σ)) and χ a = (X µ a (σ), ϕ a (σ)) with µ the spacetime index. Here σ α = (σ 0 , σ i ) with i = 1, · · · d − 1 are coordinates of a "fluid spacetime" labelling fluid elements and their internal clocks. X µ (σ α ) gives "physical" spacetime coordinates X µ of a fluid element labelled by σ α as in the standard Lagrange description of fluid flows, with X µ a describing the corresponding noises. As in the diffusion example of section 4 ϕ r,a are the charge diffusion mode associated with U(1) and the corresponding noise. β(σ) is the local inverse temperature. We can write where T 0 = 1 β 0 is the temperature at infinities where we take all external sources and dynamical fields to vanish. β 0 is the parameter appearing in the KMS condition (1.11).
The external sources are g µν (x), A µ (x) and g aµν (x), A aµ (x), with g µν (x) the spacetime metric. They are defined in physical spacetime. Here x ≡ x µ denotes physical spacetime coordinates and should be distinguished from dynamical variables X µ (σ).
The theory can be formulated either in fluid spacetime as the above variables indicate or in physical spacetime by inverting X µ (σ α ). Below we will work in the fluid spacetime as it is more convenient for introducing ghost partners and writing down a supersymmetric action.

JHEP01(2018)040
Quantities can be pulled-back or pushed-forward between the fluid and physical spacetimes through X µ (σ α ). For example, from ϕ a (σ α ) we can obtain ϕ ′ a (x) = ϕ a (σ(x)) where σ α (x µ ) is the inverse function of X µ (σ). Unless otherwise specified below for notational simplicity we will always use the same notation for a quantity and its push-forward (or pull-back), i.e. write ϕ ′ a (x) simply as ϕ a (x) and distinguish ϕ a (x) from ϕ a (σ) either by its argument or from context.
In the classical limit the action depends on the dynamical variables and external sources only through certain combinations, more explicitly (for more details see [11]) with L Xa denotes the Lie derivative along the vector X µ a (x) ≡ X µ a (σ(x)). In addition to (1.6)-(1.8), the action should also be invariant under separate spatial and time diffeo- as well as (4.3). The local velocity and chemical potential are defined as We also introduce a local temperature vector in fluid spacetime and it push-forward in physical spacetime The dynamical KMS transformation on bosonic variables can be written as [11] The quantities in (5.3)-(5.4) then transform as where L β denotes the Lie derivative along the vector β α .

JHEP01(2018)040
For our discussion of BRST symmetry below, we will need to use vielbeins 12 for various quantities in (5.3) Under dynamical KMS transformation we havẽ The action can thus also be considered as Note that in writing down explicit terms we will often need to use the inverse f α A of f A α , whose transformation can be worked out from the above.

(5.21)
We should stress that the above transformations for e A µ , A µ , η aµ and m A aµ are defined in physical spacetime. When pulled back to fluid spacetime we then have, for example, We do not introduce any ghost partner for β(σ) and require the BRST transformation of β be such that δβ α = 0: Note that the transformation of β is complicated as h 00 consists of that of X µ and e A µ . From now on we can just simply regard β α as a BRST invariant field. As a result we have 12 Below expressions come from → 0 limit of g1µν = ηABe A 1µ e B 1ν and g2µν = ηABe A 2µ e B 2ν with g1,2µν = gµν ± 2 gaµν and e A 1,2µ = e A µ ± 2 e A aµ .

JHEP01(2018)040
Since in the action (5.19) only f A α , f A aα , B α , B aα can appear, let us now construct objects which contain them and at the same time have good transformation properties under (5.20)- (5.21). This process is facilitated by considering BRST superfields defined as follows: Motivated from (5.4) and (5.16) now consider and similarly with ∂ α X µ L Xa E A µ (X ). Various quantities in (5.28)-(5.30) are given by In the above expressions L γ , L γa are defined as usual, while which is covariant explicitly. We thus find that the BRST extensions of f A aα , B aα are respectively F A aα , B aα and Let us now consider the dynamical KMS transformation on ghost variables. As in previous examples, by examining the quadratic action we propose that Again the transformations of sources are given in physical spacetime and when pulled back to the fluid spacetime the arguments should also transform, e.g.

JHEP01(2018)040
where we have used (5.9). Applying (5.9)-(5.10) and (5.38)-(5.39) to (5.28)-(5.30) we find thatf and the action can be constructed using the following two superfields and their (super)-derivatives. In contrast to the examples of last two sections here the right hand side of the supersymmetric algebra (5.45) depends on dynamical fields. In particular, β α contains some complicated dependence on f A α (recall (5.7)). Thus a supersymmetric transformation no longer preserves the total power of fields, so the proof of supersymmetrization theorem of section 3.2 and appendix B cannot be immediately applied. This potential problem can be avoided as follows. By using the freedom of time reparameterization (5.5) we can set β α = β 0 δ α 0 by choosing Then the supersymmetry becomes global, i.e.
{Q,Q} = iβ 0 ∂ 0 (5.48) and the proof can be applied at least when expanding the action around an equilibrium configuration. Note that in this gauge local temperature is expressed through f A 0 via (5.47).

Conclusions and discussions
We first summarize the main results of the paper. Formulating a consistent non-equilibrium effective field theory for a system in local equilibrium requires imposing BRST and dynamical KMS symmetries. We showed that BRST and dynamical KMS symmetries always lead to an emergent supersymmetry. Conversely, supersymmetry provides a convenient way to impose BRST symmetry and the special dynamical KMS condition. Starting from JHEP01(2018)040 a supersymmetric action one can then construct a BRST and dynamical KMS invariant theory through a simple procedure (2.14). We have discussed a few explicit examples in detail, in particular completing the formulation of fluctuating hydrodynamics of [10,11] by understanding how to introduce ghosts and implement various symmetries in the ghost sector. Our results have important implications for studying non-equilibrium questions. Ghosts run in loops and thus their dynamics plays an important role in understanding physical effects of statistical fluctuations on a physical process or physics observables. For example, it would be interesting to re-examine various dynamical critical phenomena in this light. With a full action for fluctuating hydrodynamics, one could explore systematically many questions related to effects of fluctuations on transport coefficients, hydrodynamic instabilities, and so on, in particular in far-from-equilibrium situations. At a technical level supersymmetry may also help find nontrivial fixed points for such non-equilibrium effective field theories. There are still a number of conceptual and technical challenges to overcome in order to generalize the current discussion to the quantum regime. At conceptual level, quantum fluctuations operate at the scale of T , 13 which is essentially the cutoff scale for an EFT, and thus makes the theory intrinsically nonlocal. This can already be seen from (1.15) which involves a translation of order T , therefore requiring that the theory should be able to resolve such a scale. Technically the transformation involves an infinite number of derivatives and thus invalidates derivative expansion. While one could take into quantum effects perturbatively by developing an expansion in , it is an interesting question whether it is possible to capture full quantum effects at low energies by relaxing locality. Indeed at quadratic order around thermal equilibrium it was found in [10] that one could write down an "effective field theory" which includes an infinite number of derivatives. In particular, in this theory BRST and dynamical KMS symmetries lead to a quantum deformed supersymmetric algebra There are immediate technical difficulties in generalizing such an algebra to nonlinear level, as acting on a finite product of local fields the left hand side is a derivation while the right hand side is not. Finally there is a potentail ambiguity in the dynamical KMS transformations for hydrodynamical variables at quantum level [11]. We hope to return to these issues in the future.

JHEP01(2018)040
Consider an action I b [χ r , φ r ; χ a , φ a ] where χ r,a denote collectively the dynamical variables and φ r,a external sources. Now let us introduce BRST partners c r,a for dynamical fields and require the full action in the presence of φ r (with φ a = 0) to be invariant under We can write this as 14 for some F where I denotes the total action including ghosts. From (A.2) under variation of φ r we then have for some operator V and under a variation of φ r the full generating functional is where in the second equality we have used that I is BRST invariant and in the third equality we have used that Q d as defined in (A.2) is a total derivative under the path integration. We have thus shown that W is independent of φ r and can then be set to zero by choosing a normalization constant.
To ensure the action is BRST invariant for any φ r , it is convenient to introduce also BRST partners η r,a for external sources φ r,a and require the full action (with all external fields turned on) to be invariant under simultaneous transformation of (A.1) and the corresponding transformations on sources δφ r = ǫη r , δη a = ǫφ a . (A.5) We thus have where Q s denotes the source part of the BRST operator. Now setting η r,a and φ a to zero we then obtain an action which is automatically invariant under Q d for all φ r .
For the examples discussed in the main text, equation (A.4) clearly applies to model A and nonlinear diffusion with conserved quantities. For fluctuating hydrodynamics, there is an interesting subtlety. To guarantee physical spacetime diffeomorphisms X → X ′ (X), the path integral measure for X's should come in the form DX G(X), in which case there is then a nontrivial variation δ G(X)/δX inside the path integral and (A.4) appears to break down. But for the effective action to be physical spacetime diffeomorphism invariant, X µ 's fermionic partner γ µ introduced in section 5.2 should also transform as a physical spacetime vector, i.e. γ µ (X) ≡ γ µ (σ(X)) → γ ′ν (X ′ ) = ∂ µ X ′ν (X)γ µ (X). Since the Jacobian for Grassmanian field γ µ is the inverse of that of X µ , the measure DXDγ is invariant under physical spacetime diffeomorphisms. The same cancellation of Jacobian exists between X µ a and its partner γ µ a . Thus there is no need for including G(X) in the integration measure and (A.4) applies.

JHEP01(2018)040 B Proof of special KMS condition implying supersymmetrization
In this appendix we prove that the most general solution to equation (3.28) has the form (3.29). The general proof is a bit involved and notation heavy. We will start by considering some simpler cases, which already captures the essence of the proof.

B.1 Simplest case
Let us first consider a simplest case for which there is only one species and all fields only depend on time (i.e. a quantum mechanical example). Equation (3.28) can now be written explicitly as where f (ω 1 ; ω 2 , · · · , ω n ) is a polynomial symmetric in last n − 1 variables. The above equation can be written more explicitly as Expanding f in polynomials of ω i , equation (B.2) is valid degree by degree. Thus without loss of generality we can take f to have degree h and write it as where f k is a symmetric homogeneous polynomial of ω 2 , · · · , ω n of degree k.
The basic idea of the proof is to use symmetric polynomials to write f in a suitable form which then enables us to solve (B.2) explicitly. Now introducing the power sum basis for symmetric homogeneous polynomials of ω 1 , · · · , ω n q k = n s=1 ω k s , k = 0, · · · , n (B.4) with some manipulations of symmetric polynomials we can further write f as where g m is a symmetric homogeneous polynomial of ω 1 , · · · , ω n of degree m and can be expanded in basis {q k } as already indicated in (B.5). We stress that in contrast to (B. q n 2 · · · q n k + a (P ) 2 q n 1 ω n 2 −1 1 · · · q n k + · · · + a (P ) k q n 1 · · · q n k−1 ω n k −1 where P = {n 1 , · · · , n k } denotes a (non-ordered) partition of h + 1 (i.e. k i=1 n i = h + 1) and the sum is over all partitions (with all possible k and {n i }). Now consider a partition P with at least one n j = 1, for which the corresponding terms in (B.6) are of two types:

JHEP01(2018)040
one is q n 1 · · · q n j−1 q n j+1 · · · q n k and the other is proportional to q 1 . The former is fully symmetric in all the ω s 's and thus belongs to the first term in (3.29). The latter vanishes by momentum conservation. We thus find that where f (s) is fully symmetric in all ω s 's and P ′ denotes those partitions with n i > 1 for all i. Substituting (B.7) into (B.2) we find where the term containing f (s) vanishes by momentum conservation. Since q n 1 q n 2 · · · q n k for different partitions are independent we thus have − a for all partitions P ′ . Substituting the above equation back into f (a) we find that where and below slash means omitting that term. Now putting (B.10) back into and using permutations among χ r 's we can write a (P ′ ) j q n 1 · · · ✚ ✚ q n j · · · q n k−1 (ω which is the form of (3.31).

B.2 The next simpler case
Now we still consider a single field species, but with full momentum k µ = (ω, k i ), i = 1, · · · , d − 1 dependence in d spacetime dimension, i.e. f in (B.1) should be understood as f (k µ 1 ; k µ 2 , · · · , k µ n ) which is a polynomial of nd variables k µ s and is symmetric under permutations of k µ s for s = 2, · · · , n. Equation (B.2) now becomes The idea for solving (B.13) is exactly the same as in last subsection with the only difference being that we now need to use multi-symmetric polynomials that are natural generalizations of symmetric polynomials [39,40].
With some manipulations we can rewrite (B.14) as (see appendix B.4 for details) ) which is the counterpart of (B.5). In (B.15) we should sum over all powers of k 1 with multi-degree η with |η| = n − 1 and h − η ∈ N d . g γ 's are homogeneous multi-symmetric polynomials of all k µ s with multi-degree γ and can be expanded in the multi-power sum basis q γ 2 · · · q γ k + a (P ) 2 q γ 1 k γ 2 −e 0 1 · · · q γ k + · · · + a (P ) where P denotes a partition of h + e 0 (i.e. k i=1 γ i = h + e 0 ) and we sum over all partitions. Note that by definition γ i − e 0 ∈ N d for all i. We can again separate the summation over P into those which have at least one |γ j | = 1 and those with none (which we denote as P ′ ). The former is either zero by momentum conservation or fully symmetric in all the k µ s (if some γ j = e 0 ), and thus where f (s) is fully symmetric in all k µ s and in P ′ all |γ j | ≥ 2. Now (B.13) gives

JHEP01(2018)040
and then − a for all partitions P ′ . Substituting the above equation back into f (a) we again have and using permutations among χ r 's we then replace k again has the structure of (3.31).

B.3 General case
We now consider the general case where χ can have arbitrary species indices. The idea is exactly the same as before except that now we view (3.26) as a multi-variable polynomial of both momenta and fields. More explicitly, definē as d + N dimensional vectors. Then I b of (3.26) (when we expand all the coefficients in momenta) can be viewed as a polynomial I b (k 1a ;k 2 , · · · ,k n ) which is symmetric in the last n − 1 variables, and equation (3.28) (after symmetrization) becomes ω 1 I b (k 1 ;k 2 , · · · ,k n ) + ω 2 I b (k 2 ;k 1 , · · · ,k n ) + · · · + ω n I b (k n ;k 1 , · · · ,k n−1 ) = 0 . (B.25) The story is similar to before and we can solve (B.25) for I b with a given multi-degree h which is now a d + N -dimensional vector. There are two important differences: 1. Momentum conservation only applies to the first d components ofk. We thus separate the multi-agree fork as γ = (γ,γ) whereγ is d-dimensional andγ is N -dimensional, i.e.k γ = kγξγ. Similarly the total multi-degree h = (ĥ,h) and the basis areê a with a = 0, 1, · · · d − 1 andẽ i with i = 1, · · · N .
2. We only need to consider polynomials of the form (3.26), i.e. |γ s | for eachk s can only be 1 or 0 and |h| = n.
We can now expand I b (k 1 ;k 2 , · · ·k n ) in the form of (B.17) with q α now defined in terms ofk, k 1 replaced byk 1 , and e 0 replaced byê 0 . In particular for any i we should havê γ i −ê 0 ∈ N d . We still have and replacingk 1 byk 1a we find that where we should also replace all thek 1 's in q γ 's byk 1a . From item 2 discussed earlier the terms relevant for I b must have |γ i | = 1 or 0 for all i's and k i=1 |γ i | = n. For a γ i with |γ i | = 1 we will then haveγ i =ẽ j i for some j i . Those with |γ i | = 0 can be written as q γ i = qγ i ≡ n s=1 kγ i s . Since for all γ i −ê 0 ∈ N d , by exchanging indices in (B.27) we can always extract an overall factor ω and the remaining factor in (B.27) is then anti-symmetric betweenk 1a and k 2 , thus giving (3.31). To see this explicitly we must keep in mind that we should only select those terms in (B.27) which has the structure of (3.26). 15 There are three types of terms in (B.27): 1. |γ i | = |γ k | = 1 for which we havē Now the two terms on the right hand side must select different factors ξ s from the product of q j 's in (B.27) so that they will have the structure of (3.26). Take some will result in terms (B.31) 15 All the rest terms must cancel themselves by definition.

JHEP01(2018)040
Multiplying it with products of q i 's in (B.27) and selecting the appropriate terms for (3.26) we always find terms of the form (3.31), due to full permutation symmetry of q i 's.
This then concludes the full proof.

B.4 Justification of expansions using symmetric polynomials
We now show that (B.5) follows from (B.3). For this purpose we can expand f m in (B.3) in power sum basis (the fundamental theorem of symmetric polynomials) of ω 2 , · · · ω n p k = n s=2 ω k s , (k = 0, · · · , n − 1) (B.32) From (B.4) we have q k = ω k 1 + p k . Plugin this relation into (B.3) and we can rewrite f in terms of the following expansion: 1g 0 (q i ) + ω h−1 1g 1 (q i ) + · · · + ω 1gh−1 (q i ) +g h (q i ) (B.33) whereg i 's are polynomials of degree i and of q k . Note this expansion in terms of q k and ω 1 is not unique although the expansion (B.3) is. Define the symmetric polynomial basis of ω 1 to ω n as σ i = 1≤k 1 <···<k i ≤n ω k 1 · · · ω k i , σ 0 = 1 (B.34) for (i = 0, · · · , n) and that of ω 2 to ω n as s i for (i = 0, · · · , n − 1). Using the following relations This expansion is equivalent to (B.33) because σ i and q i can be uniquely expanded by each other. Following this process, one can always reduce h in (B.33) to be less than n and (B.5) results. We now justify (B.15). Like symmetric polynomials, any multisymmetric polynomial can be uniquely expanded by elementary multisymmetric polynomials σ α = {i j,k } 1≤i j,k ≤n, all different ω i 0,1 · · · ω i 0,α 0 k 1 i 1,1 · · · k 1 for all posible choices of α with 0 ≤ |α| ≤ n. We define σ α for k 1 to k n and s α for k 2 to k n . For example, σ (2,1,1) = i 1 <i 2 1≤i j ≤n, all different

JHEP01(2018)040
We define q α for k 1 to k n and p α for k 2 to k n . For example, These two basis can be expanded via each other uniquely. As a necessary condition, by the notation of d-vector α, we see the number of these two basis are the same. Expanding f γ in (B.14) in terms of p α and using the relations we can rewrite (B.14) in the form of 1g e j (q α ) + for any |α| = n and |β| ≥ 0. Using this method, we can reduce (B.42) to an expansion with the length of maximal exponent of k 1 to be less than n and gives (B.15).
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.