text:
# I Introduction

## I Introduction

Wilsonian renormalisation group techniques  like the exact renormalisation group (ERG)  are important tools for addressing non-perturbative problems within quantum field theory (for recent reviews, see Refs. ). They are similar in spirit to the block-spin action invented in condensed matter physics, and their particular strength is their flexibility, allowing for systematic approximations without being tied to the small coupling region. The ERG is based on an infrared (IR) regularisation with momentum scale parameter $`k`$ of the full propagator, which turns the corresponding effective action into a scale dependent functional $`\mathrm{\Gamma }_k`$. The ERG flow describes the change of the effective action under an infinitesimal variation of the IR scale $`k`$. It thereby interpolates between the initial UV action $`\mathrm{\Gamma }_{k=\mathrm{\Lambda }}`$ and the full quantum effective action $`\mathrm{\Gamma }\mathrm{\Gamma }_{k=0}`$. Although the flow depends explicitly on the specific infrared regulator chosen, the endpoint of the integrated full flow does not.

An explicit computation of the IR effective theory based on the ERG flow requires the specification of the field content, the initial condition $`\mathrm{\Gamma }_\mathrm{\Lambda }`$ and the choice of a particular IR regulator. The UV initial condition is typically given by the classical action. Hence, the main physical information is contained in the ERG flow itself. Most problems of physical interest are too complex to be solved exactly and an application of this formalism – as of any other method – is bound to certain approximations. Furthermore, the flow equation is equivalent to infinitely many coupled partial differential equations, which would seem very difficult to solve exactly. Therefore, one has to resort to some approximations or truncations which allow, at least in principle, for a systematic computation of the full quantum effective action. In order to provide reliable physical predictions, like a high precision computation of universal critical exponents, it is mandatory to provide a good control for approximated ERG flows.

A number of systematic expansion schemes for flow equations are known, including standard perturbation theory. Non-perturbative expansions of the effective action, not bound to the weakly coupled regime, are the derivative expansion, expansions in powers of the fields, or combinations thereof. For example, the leading order of the derivative expansion retains only an effective potential and a standard kinetic term, and contains non-perturbative information as it corresponds to the resummation of infinitely many perturbative loop diagrams. The study of approximate quantum effective actions along these lines is a sensible procedure since the underlying expansions admit a systematic improvement to higher order.

Solutions to truncated flow equations display a spurious dependence on the IR regulator . This is similar to the scheme dependence of physical observables observed within perturbative QCD , or the truncation dependence of solutions to Schwinger-Dyson equations. Its origin is the following. The IR regulator couples, through the flow equation, to all vertex functions of the theory. The flow trajectory of the functional $`\mathrm{\Gamma }_k`$ in the space of all effective action functionals depends on the regulator. Hence, the regulator – while regulating the flow – also modifies the effective interactions at intermediate scales $`k0`$. In other words, the effective action at intermediate scales still has some memory of the details of how the integrating-out of degrees of freedom has been performed. This regulator dependence is of no relevance for the full flow. Eventually, the convergence towards the full quantum effective action for any regulator ensures that all regulator-induced interactions cancel out in the physical limit. Approximations imply that certain vertex functions are neglected. Then, not all regulator-induced interactions cancel out for $`k0`$: the missing back coupling of the neglected vertex functions is responsible for regulator-dependent terms in the physical limit. In consequence, approximations to the full quantum effective action depend spuriously on the scheme.

Recently, a new line of reasoning has been put forward which essentially turns this observation around : given that the solution of a truncated flow depends on the regulator, it should be possible to identify specific ones which “optimise” the physical content of a given approximation. Optimised regulators stabilise the flow and lead to a faster convergence of expansions, such that the main physical information is almost exclusively contained within a few leading terms, and higher order contributions remain small .

In Ref. , we have derived a simple and generic optimisation criterion for ERG flows, based only on the full inverse propagator at vanishing field. Given the set of possible IR regulators, the criterion allows to distinguish the quality of regulators in the sense outlined above. In the present paper we study a specific “optimised” regulator for both bosonic and fermionic degrees of freedom. To be more explicit, we introduce the ERG flow for the effective action . For bosonic fields $`\varphi `$, it is given by

$$_t\mathrm{\Gamma }_k[\varphi ]=\frac{1}{2}\mathrm{Tr}\left(\frac{\delta ^2\mathrm{\Gamma }_k}{\delta \varphi (q)\delta \varphi (q)}+R_k\right)^1_tR_k.$$
(1)
Here, the trace denotes a sum over all loop momenta and indices, and $`t=\mathrm{ln}k`$ is the logarithmic scale parameter. The flow has a simple one-loop structure. The Wilsonian “integrating-out” is achieved by the infrared regulator $`R_k`$. It regulates the propagator for small momenta, while the insertion $`_tR_k`$ cuts off the large-momentum contributions. In total, only a small momentum window about $`q^2k^2`$ contributes to the flow. Apart from a few constraints displayed later, the function $`R_k`$ can be chosen at will. A “good” choice for the regulator function is at the root for reliable physical predictions, and we consider, for the bosonic fields, the optimised regulator

$$R_k^{\mathrm{opt}}(q^2)=Z_k(k^2q^2)\mathrm{\Theta }(k^2q^2),$$
(2)
where $`Z_k`$ is an appropriately defined wave function renormalisation. This regulator is particularly simple: for loop momenta $`q^2>k^2`$ it vanishes identically and the effective propagator appearing in the flow equation is not modified; for loop-momenta $`q^2<k^2`$ it acts like a momentum-dependent mass term in such a way that the inverse effective propagator $`q^2+R_k(q^2)`$ becomes a momentum independent constant. In consequence, the effective infra-red propagator does no longer distinguish between the different modes with $`q^2<k^2`$.

Optimised flows based on Eq. (2) derive from a generic optimisation criterion , and have a number of remarkable properties. The optimised flow leads to the fastest decoupling of heavy modes, in accordance with the decoupling theorem . In the limit $`k0`$, optimised flows smoothly approach a convex effective action, owing to a simple analytic pole of the flow . At non-vanishing temperature, the optimised flow factorises: the contributions from thermal and quantum fluctuations are disentangled, unlike for generic flows. A similar factorisation of the flow holds to second order in the derivative expansion for field-independent wave function renormalisations, and a partial factorisation is found for the general case. Finally, the optimised flow has a very simple analytic structure. This facilitates their study and is helpful for both analytical or numerical considerations. All these properties lead to a stabilisation of the flow and an improved convergence towards the physical theory. Analogous results for fermionic flows are discussed as well.

We also study this question within an RG formalism based on a proper-time regularisation of the operator trace for the one-loop effective action , which we call “proper-time renormalisation group” (PTRG) for short. In contrast to the ERG, the PTRG has no path integral derivation, which makes the conceptual reasoning more difficult . Still, owing to the close similarity to the ERG at leading order in the derivative expansion, it is possible to identify the analogue of Eq. (2) for the PTRG.

The format of the paper is as follows. We introduce the physical ideas behind the generic optimisation condition. Explicit realisations for bosonic and fermionic degrees of freedom are introduced as well (Section II). The main characteristics of optimised flows are discussed to leading order in the derivative expansion, and contrasted with those of generic flows (Section III). We then turn to the discussion of quantum field theories at finite temperature. We show that optimised thermal flows factorise on the level of the flow equation, unlike generic flows. A simple physical explanation for the factorisation is provided (Section IV). Next, we consider the extension to higher orders in the derivative expansion. The cases of field dependent or independent wave-function renormalisations are both discussed, and a similar factorisation for optimised flows is established (Section V). Finally, we provide the corresponding optimised proper-time cut-off for the PTRG (Section VI). Due to the qualitative difference of the topics studied, we discuss our findings separately at the end of the corresponding sections. We close with a summary and an outlook (Section VII). Three Appendices contain technical details and explicit expressions for optimised flows.

## II Optimisation

In this section, we discuss a generic optimisation criterion for ERG flows for Euclidean quantum field theories. In particular, we provide a simple and explicit optimised regulator for both bosonic and fermionic flows. Prior to this, we have to review a few basic properties of IR regulator functions, which are at the root of the subsequent considerations.

### A Regulators

The flow equation (1) is defined through the infrared regulator functions $`R_k(q^2)`$ and $`R_{F,k}(q^2)`$, respectively . These operators depend on an infrared scale $`k`$, which induces a scale dependence. When written in terms of the scale-dependent effective action $`\mathrm{\Gamma }_k`$, the scale dependence is given precisely by the flow equation (1). The right-hand side of Eq. (1) contains the full inverse propagators and the trace denotes a sum over all indices and integration over all momenta.

The regulator scheme (RS) functions can be chosen at will, however, within some basic restrictions. These restrictions ensure that the flow equation is well-defined, thereby interpolating between an initial action in the UV and the full quantum effective action in the IR. More specifically, it is required that

$$\underset{q^2/k^20}{lim}R_k(q^2)>0.$$
(3)
This ensures that the effective propagator at vanishing field remains finite in the infrared limit $`q^20`$, and no infrared divergences are encountered in the presence of massless modes. This property makes $`R_k`$ an infrared regulator. If the limit (3) is finite, we call the corresponding regulator masslike. The second requirement is the vanishing of $`R_k`$ in the infrared,

$$\underset{k^2/q^20}{lim}R_k(q^2)0.$$
(4)
This guarantees that the regulator function is removed in the physical limit, where the scale-dependent effective action $`\mathrm{\Gamma }_k`$ reduces to the quantum effective action $`\mathrm{\Gamma }=lim_{k0}\mathrm{\Gamma }_k`$. The third condition to be met is

$$\underset{k\mathrm{\Lambda }}{lim}R_k(q^2)\mathrm{}.$$
(5)
This way it is ensured that $`\mathrm{\Gamma }_k`$ approaches the microscopic action $`S=lim_{k\mathrm{\Lambda }}\mathrm{\Gamma }_k`$ in the UV limit $`k\mathrm{\Lambda }`$. In the rest of the paper, we put $`\mathrm{\Lambda }=\mathrm{}`$ for the UV scale, although our main line of reasoning can be applied for finite $`\mathrm{\Lambda }`$ as well. With this choice, the regulator function depends only on $`q^2`$ and $`k^2`$, and it is convenient to introduce a dimensionless function $`r(q^2/k^2)`$ as

$$R_k(q^2)=Z_kq^2r(q^2/k^2).$$
(6)
with $`Z_k`$ an appropriate wave function renormalisation (cf. Section V); $`Z_k1`$ to leading order in the derivative expansion. Owing to the general conditions imposed on the regulator, the function $`r(y)`$ ranges between $`0r(y)\mathrm{}`$.

Another condition concerns the proper normalisation of the regulator. The normalisation fixes the scale at which the IR regulator becomes effective. Let us employ the condition

$$R_k(q^2=c_Bk^2)=Z_kc_Bk^2$$
(7)
for bosons (a similar condition holds for fermions, see below) and $`c_B>0`$.In Ref.  the convention $`c_B=1`$ has been used. The normalisation translates into the condition $`r(c_B)=1`$. Two different choices for $`c_B`$ can always be mapped onto each other through a rescaling of the IR scale $`k`$. Hence, a proper normalisation is only of relevance for a comparison of different regulators (as done in Ref. ), or for theories containing different bosonic and/or fermionic degrees of freedom, where the relative normalisation of the regulators can become important.

### B Optimisation criterion

Here, we discuss an optimisation criterion for ERG flows, which ensures that flows like Eq. (1) and approximations to it have good convergence and stability properties. Following Ref.  (see also Ref. ), we first provide the general criterion for optimised choices of RS functions. Then, more specifically, we apply this idea to bosonic and fermionic theories with standard kinetic terms.

The physical information of the flow equation (1) is contained in the full effective inverse propagator, which is given by

$$\frac{\delta ^2\mathrm{\Gamma }_k[\varphi ]}{\delta \varphi (q)\delta \varphi (q)}+R_k(q^2).$$
(8)
Notice that Eq. (8) depends both on the fields and on the RS function. The ERG flow is well-defined as long as the full inverse propagator displays a gap,

$$\underset{q^20}{\mathrm{min}}\left(\frac{\delta ^2\mathrm{\Gamma }_k[\varphi ]}{\delta \varphi (q)\delta \varphi (q)}|_{\varphi =\varphi _0}+R_k(q^2)\right)=Ck^2>0.$$
(9)
The functional derivative is evaluated at a properly chosen expansion point $`\varphi _0`$. The existence of the gap $`C>0`$ implies an IR regularisation. Furthermore, the gap is a prerequisite for the ERG formalism. Elsewise, Eq. (1) becomes singular at points where the full inverse effective propagator develops zero modes.The case $`C=0`$ indicates that a saddle point expansion about $`\varphi _0`$ is not applicable. Those points $`\varphi _0`$ in field space with $`C=0`$ correspond to an instability. The problem can be cured by chosing a more appropriate expansion point such that $`C>0`$. For related literature, see Ref. . The size of the gap $`C`$ in Eq. (9) depends both on the RS function and on dimensionless parameters like $`\varphi _0^2/k^2`$ or mass ratios, specific to the particular theory studied.

A natural optimisation criterion based on Eq. (9) consists of maximising the gap $`C`$ over the space of all possible RS functions. Optimised RS functions are those for which the maximum of $`C`$ is attained. The optimisation ensures that the momentum-dependent kernel of the ERG flow is the most regular. Therefore we expect that optimised flows are much more stable against approximations and show better convergence properties.

The optimisation condition as formulated above is, essentially, only sensitive to the momentum dependence of the full inverse propagator. Dropping momentum-independent terms on the left-hand side of Eq. (9) changes the number $`C`$ accordingly, but leaves the explicit dependence on $`R_k(q^2)`$ unchanged. Therefore, the optimisation leads to the same set of optimised RS functions as long as the implicit dependence of $`\mathrm{\Gamma }_k^{(2)}[\varphi ]\delta ^2\mathrm{\Gamma }_k[\varphi ]/\delta \varphi (q)\delta \varphi (q)`$ on the RS function remains negligible. For this reason, the optimisation condition of Refs.  is based only on the momentum-dependent terms of Eq. (8).

From now on, we concentrate on a standard kinetic term. The effect of a field-dependent wave function renormalisation can be taken into account as well (see Section V below). We expand the full inverse propagator as $`Z_k(q^2+Z_k^1R_k(q^2)+\mathrm{})`$ about the regularised kinetic term. Finally, dropping the momentum-independent terms transforms Eq. (9) into

$$\underset{q^20}{\mathrm{min}}\left(q^2+Z_k^1R_k(q^2)\right)=Ck^2>0.$$
(10)
A far reaching consequence of the infrared regulator in Eq. (10) is the presence of a gap for all $`k>0`$, which follows trivially from Eq. (3). The decisive difference between Eq. (9) and Eq. (10) is that the size of the gap $`C>0`$ in Eq. (10) depends only on the particular choice for the RS, but not on the specific theory. Rewriting Eq. (10) in dimensions of $`k`$ leads to

$$P^2(y)q^2/k^2+R_k(q^2)/(Z_kk^2)=y[1+r(y)],$$
(11)
where $`yq^2/k^2`$. Expressed in terms of Eq. (11), the size of the gap is given by

$$C=\underset{y0}{\mathrm{min}}P^2(y).$$
(12)
Any RS function is now characterised by the associated gap $`C`$. The size of the gap can be made arbitrarily small. Effectively, this corresponds to removing the IR regulator in the first place. However, for fixed normalisation $`c_B`$, it cannot be made arbitrarily large, $`C<\mathrm{}`$. Hence, the natural optimisation condition, which is the requirement to maximise the gap, becomes

$$C_{\mathrm{opt}}=\underset{(\mathrm{RS})}{\mathrm{max}}\left(\underset{y0}{\mathrm{min}}P^2(y)\right).$$
(13)
A few comments are in order. The maximum in Eq. (13) is taken over the (infinite-dimensional) space of all possible RS functions. The number $`C_{\mathrm{opt}}`$ is uniquely determined and reads $`C_{\mathrm{opt}}=2c_B`$, where $`c_B`$ is the normalisation of bosonic regulators. From now on, we refer to Eq. (13) as an “optimisation condition”, and all RS functions for which $`C=C_{\mathrm{opt}}`$ are called solutions to the optimisation condition. The space of solutions to the optimisation condition is infinite-dimensional. Notice also that the condition to minimise the gap is not an extremisation linked to the regulator, because it corresponds to removing the IR regularisation. In Ref. , a variety of different solutions have been found, and some examples are given in Fig. 1 below.

In order to obtain Eq. (10), we have assumed a standard kinetic term for the fields. Therefore, the resulting optimisation condition Eq. (13) is independent on the specific theory. Once the momentum-dependent part of $`\mathrm{\Gamma }_k^{(2)}`$ depends on the fields, the corresponding optimisation condition based on the momentum-dependent part of Eq. (9) is sensitive to the specific theory. Within a derivative expansion, this happens starting from the second order (cf. the discussion in Section V).

The optimisation condition has a number of interpretations in more physical terms (cf. Refs. ). It has been shown that the radius of convergence for amplitude expansions is given by $`C`$. Therefore the optimisation condition improves their convergence. Furthermore, it leads to a smooth approach towards a convex effective potential in the IR limit $`k0`$. It has also been shown that it improves the convergence of the derivative expansion . Finally, it is worth emphasizing that the optimisation criterion is a rather mild condition: all regulator functions are described by at most countably infinitely many parameters, because $`R_k`$ is at least square integrable. Of these, only one parameter is fixed by the optimisation criterion.

We now turn to the discussion of fermionic degrees of freedom $`\psi `$ and $`\overline{\psi }`$ . The flow equation is given by

$$_t\mathrm{\Gamma }_k[\psi ,\overline{\psi }]=\mathrm{Tr}\left(\frac{\delta ^2\mathrm{\Gamma }_k}{\delta \psi (q)\delta \overline{\psi }(q)}+R_{F,k}\right)^1_tR_{F,k}.$$
(14)
As usual, the trace sums over all loop momenta and indices. The constraints on the function $`R_{F,k}`$ are similar to those on $`R_k`$ . Following Ref. , we chose the regulator proportional to $`q\text{/}`$ and introduce

$$R_{F,k}(q)=Z_{F,k}q\text{/}r_F(q^2/k^2)$$
(15)
We chose the normalisation as

$$R_{F,k}^2(q^2=c_Fk^2)=c_Fk^2.$$
(16)
This translates into the condition $`r_F(c_F)=1`$. It has been shown that the fermionic analogue of the function Eq. (11) is given by

$$P_F^2(y)=y[1+r_F(y)]^2.$$
(17)
Therefore, we can define the fermionic gap as

$$C_F=\underset{y0}{\mathrm{min}}P_F^2(y),$$
(18)
and the corresponding optimisation condition reads

$$C_{F,\mathrm{opt}}=\underset{(\mathrm{RS})}{\mathrm{max}}\left(\underset{y0}{\mathrm{min}}P_F^2(y)\right).$$
(19)
The optimised fermionic gap is uniquely determined through the normalisation $`c_F`$ as $`C_{F\mathrm{opt}}=4c_F`$. Conceptually, the fermionic case is treated in the same way as the bosonic case. The sole difference stems from the fact that the bosonic kinetic term contains two derivatives, while the fermionic kinetic term contains only one. Therefore, the functions (11) and (17) entering the optimisation condition are different.

### C Derivation of optimised bosonic and fermionic regulators

A lot of effort has been spent in order to provide explicit regulators which lead to sufficiently simple and analytic ERG flows. For example, the sharp cut-off provides a simple explicit flow to leading order in the derivative expansion. For this reason, it is one of the most intensively studied flows in the field (cf. Refs. ). Other attempts have been made based on power-like regulators $`R_kq^2(k^2/q^2)^b`$ for $`b=1`$ and $`b=2`$ , or variants of a mass-term regulator $`R_kk^2\mathrm{\Theta }(k^2q^2)`$. These regulators are still sufficiently simple from an algebraic point of view, and lead to reasonably simple flows.<sup>§</sup><sup>§</sup>§Of these, only the power-like regulator with $`b=2`$ solves the optimisation condition Eq. (13). However, in absence of an underlying “guiding principle” it was not obvious how to make progress given the plethora of possible regulators, and in particular, how to distinguish the “quality” of the corresponding flows.

Here, in turn, we take full advantage of the existence of a guideline provided by the optimisation criterion. We propose a regulator which (i) solves the optimisation criterion, (ii) is based on an additional stability criterion for approximate flows, and (iii) leads to simple explicit expressions for the corresponding flows. The heuristic derivation runs as follows. The space of regulators which solve the optimisation criterion is still infinite dimensional. Let us seek for a “simple” solution to Eq. (13). The simplest one corresponds to an inverse propagator which is flat, i.e. momentum-independent, $`P^2C_{\mathrm{opt}}`$. Take $`C_{\mathrm{opt}}=1`$. This immediately implies, using Eqs. (6) and (11), that $`R_k(q^2)=k^2q^2`$. Our naive Ansatz is consistent with Eqs. (3) and (5), but not with the main requirement Eq. (4) for small $`k^2<q^2`$. In order to fulfil Eq. (4), the regulator has to be cut-off above some loop momenta. Therefore, a natural proposal for the bosonic case consists in taking

$$R_k^{\mathrm{opt}}(q^2)=(k^2q^2)\mathrm{\Theta }(k^2q^2).$$
(20)
The ultraviolet modes $`q^2>k^2`$ are not touched by this regulator because Eq. (20) vanishes identically for $`q^2>k^2`$. In turn, for all modes with $`q^2k^2`$ the regulator acts as a momentum-dependent mass term $`(k^2q^2)`$ with the infrared limit $`k^2`$ for vanishing momenta. It is a masslike regulator. By construction, the inverse propagator at vanishing field Eq. (8) becomes momentum independent for all $`q^2k^2`$ (see Fig. 1). It is this property which is responsible for the main characteristics of the regulator: all infrared momentum modes below the scale $`k`$ are treated in the same way since the effective inverse propagator does no longer distinguish between them.

The dimensionless regulator function $`r(y)`$ is defined in Eq. (6). With the choice made in Eq. (20) it follows that

$$r_{\mathrm{opt}}(y)=\left(\frac{1}{y}1\right)\mathrm{\Theta }(1y).$$
(21)
The regulator function is normalised with $`c_B=\frac{1}{2}`$. Such a normalisation can always be achieved. Other normalisations are considered at the end of Section III.

In Fig. 1, we have displayed the effective inverse propagator at vanishing field for different optimised regulators. The full line represents Eq. (20), and the thin dashed line corresponds to $`R=0`$ (no regulator). The set of dashed lines correspond to the optimised regulators discussed in Fig. 3 of Ref. . Here, they have been given for comparison. Notice that all curves cross in the normalisation point $`r(c_B)=1`$. All optimised propagators display the same gap $`C_{\mathrm{opt}}=2c_B`$, but differ essentially in the curvature around their minima.

The fermionic analogue of Eq. (20) is derived in essentially the same way. Starting with Eq. (17), imposing $`P_F^21`$ for small momenta, and the general conditions (3), (4) and (5), we finally end up with

$$R_{F,k}^{\mathrm{opt}}(q)=q\text{/}\left(\sqrt{\frac{k^2}{q^2}}1\right)\mathrm{\Theta }(k^2q^2),$$
(22)
normalised with $`c_F=\frac{1}{4}`$. In terms of a dimensionless function $`r_F(q^2/k^2)`$, Eq. (22) becomes

$$r_{F,\mathrm{opt}}(y)=\left(\frac{1}{\sqrt{y}}1\right)\mathrm{\Theta }(1y).$$
(23)
and $`r_F(\frac{1}{4})=1`$. The non-analyticity of Eq. (23) is a direct consequence of $`R_{F,k}`$ having only one mass dimension. We shall see below that it is of no harm for the computation of fermionic flows because Eq. (23) enters only in specific combinations such that the non-analyticity disappears.

## III Derivative expansion

The flow equation (1) is a functional differential equation, which, from a technical point of view, is equivalent to infinitely many coupled partial differential equations for the couplings parametrising the effective action $`\mathrm{\Gamma }_k`$. A number of different systematic approximation procedures for flows are known. In this section, we consider flows to leading order in the derivative expansion, based on expanding the operators of the effective action according to the number of derivatives . This leads to a closed set of coupled partial differential equations for the coefficient functions. We discuss the main structure of optimised flows and contrast it with generic ones.

### A Specific flows

In order to make our subsequent reasoning more transparent, it is useful to have an explicit example at hand. To that end, we consider an $`O(N)`$-symmetric real scalar field theory in $`d`$ dimensions, the linear sigma model. To leading order in the derivative expansion we make the Ansatz

$$\mathrm{\Gamma }_k=d^dx\left(U_k(\overline{\rho })+\frac{1}{2}Z_k(\overline{\rho })_\mu \varphi ^a_\mu \varphi _a+\frac{1}{4}Y_k(\overline{\rho })_\mu \overline{\rho }_\mu \overline{\rho }+𝒪(^4)\right)$$
(24)
for the effective action, with $`\overline{\rho }=\frac{1}{2}\varphi ^a\varphi _a`$. For $`N1`$, there are two independent wave function factors $`Z_k`$ and $`Y_k`$ beyond leading order in this expansion (cf. Section V). To leading order in the derivative expansion, the flow equation (1) reduces to a flow for the effective potential, $`_tU_k`$. The main physical applications concern the non-trivial Wilson-Fisher fixed point in $`d=3`$ and the computation of related universal quantities.

Inserting the Ansatz (24) into the basic flow equation, and using $`ZY1`$, yields

$$_tU_k=\frac{1}{2}\frac{d^dq}{(2\pi )^d}\left(\frac{(N1)_tR_k(q^2)}{q^2+R_k(q^2)+U_k^{}(\overline{\rho })}+\frac{_tR_k(q^2)}{q^2+R_k(q^2)+U_k^{}(\overline{\rho })+2\overline{\rho }U_k^{\prime \prime }(\overline{\rho })}\right).$$
(25)
It is a second order non-linear partial differential equation. One easily recognizes the contributions from the $`N1`$ “Goldstone” modes and the “radial” mode. A similar flow equation has been obtained for the wave function renormalisations $`Z_k`$ and $`Y_k`$ . The momentum integration is regularised in the UV, owing to the regulator term $`_tR_k(q^2)`$ in the numerator, and in the IR due to $`R_k(q^2)`$ in the denominator.

### B Generic flows

For convenience we perform the angular part of the momentum integration and rewrite the right-hand side of Eq. (25) in terms of so-called threshold functions  as

$$_tU_k(\overline{\rho })=2v_d(N1)k^d\mathrm{}_0^d\left(\frac{U_k^{}(\overline{\rho })}{k^2}\right)+2v_dk^d\mathrm{}_0^d\left(\frac{U_k^{}(\overline{\rho })+2\overline{\rho }U_k^{\prime \prime }(\overline{\rho }))}{k^2}\right).$$
(26)
The constants $`v_d`$ are given by

$$v_d^1=2^{d+1}\pi ^{d/2}\mathrm{\Gamma }(\frac{d}{2}),$$
(27)
and the functions $`\mathrm{}_n^d(\omega )`$ are defined as

$$\mathrm{}_n^d(\omega )=\left(\delta _{n,0}+n\right)_0^{\mathrm{}}𝑑yy^{{\scriptscriptstyle \frac{d}{2}}1}\frac{y^2r^{}(y)}{(P^2(y)+w)^{n+1}}.$$
(28)
While the flow (26) is specific for the theory defined by Eq. (24), the functions (28) are not. These functions describe the generic structure of the flow to leading order in the derivative expansion. The flows for different indices $`n>0`$ are related by

$$_\omega \mathrm{}_n^d(\omega )=(n+\delta _{n,0})\mathrm{}_{n+1}^d(\omega ).$$
(29)
Therefore, it suffices to study the flows $`\mathrm{}_0^d(\omega )`$.

The fermionic analogue of the flow (28) is

$$\mathrm{}_{F,n}^d(\omega )=\left(\delta _{n,0}+n\right)_0^{\mathrm{}}𝑑yy^{d/2}\frac{2y(1+r_F)r_F^{}}{(P_F^2(y)+w)^{n+1}}.$$
(30)
and Eq. (29) holds equally for Eq. (30). Notice the additional factor $`1r_F`$ in the integrand, which arises due to the Dirac structure of Eq. (22). We have used the normalisation condition $`r_F(\frac{1}{4})=1`$.

It is evident that the characteristics of the flow, determined by the choice of $`R_k`$, are entirely encoded within the functions Eqs. (28) and (30) (or similar functions to higher order in the derivative expansion). For a generic regulator, these are complicated functions of the fields, which can be computed explicitly only for very specific choices for the regulator.

Two properties of generic flows given in terms of $`\mathrm{}_0^d(\omega )`$ and $`\mathrm{}_{F,0}^d(\omega )`$ are worth being mentioned. First of all, from their very definition and the constraints imposed on the regulator function, we conclude that any function (28) for $`n=0`$ decays at most as $`1/\omega `$ for $`\omega \mathrm{}`$ . Therefore, they describe the decoupling of “heavy” modes from the flow, which is a manifestation of the decoupling theorem . Secondly, all flows have a pole in $`C+\omega `$, where $`C`$ denotes the gap. Both the analytical structure and the strength of the pole depend on the regulator. From the general requirements for regulators, and the explicit form of (28), it follows that the pole for $`n=0`$ cannot be stronger than a simple analytical pole $`1/(C+\omega )`$. The pole of threshold functions has important physical implications. It determines the approach to a convex effective potential for theories within a phase of spontaneous symmetry breaking .

### C Optimised flows

Now we turn our attention to the optimised regulators introduced in Eqs. (20) and (22). The evaluation of Eq. (28) is particularly simple because the $`\mathrm{\Theta }`$-function cuts off the momentum integration. Using Eq. (21), Eq. (28) reduces to two terms,

$$\mathrm{}_0^d(\omega )=\frac{1}{1+w}_0^1𝑑yy^{{\scriptscriptstyle \frac{d}{2}}1}+_0^{\mathrm{}}𝑑yy^{{\scriptscriptstyle \frac{d}{2}}1}\frac{(y1)\delta (1y)}{1y\mathrm{\Theta }(1y)+w}.$$
(31)
In the first term, the momentum integration is cut-off above $`y1`$. The function $`P^2(y)`$ remains a constant in this momentum regime, which allowed to move the $`\omega `$-dependent term in front of the momentum integration. The integrand of the second term contains products of distributions. Since the integrand is proportional to $`(1y)\delta (1y)`$ the second term vanishes identically, independently of the specific implementation for the $`\mathrm{\Theta }`$-function. The remaining momentum integration of the first term becomes trivial and gives

$$\mathrm{}_0^d(\omega )=\frac{2}{d}\frac{1}{1+\omega }.$$
(32)
We used the normalisation $`r(\frac{1}{2})=1`$ and hence $`P^2=1`$ for $`y1`$.

For fermionic flows (30) and the regulator (22), we find

$$\mathrm{}_{F,0}^d(\omega )=\frac{1}{1+\omega }_0^1𝑑yy^{d/21}+2_0^{\mathrm{}}𝑑yy^{{\scriptscriptstyle \frac{d}{2}}+1}\frac{[1+(\frac{1}{\sqrt{y}}1)\mathrm{\Theta }(1y)](\frac{1}{\sqrt{y}}1)\delta (1y)}{[\sqrt{y}+(1\sqrt{y})\mathrm{\Theta }(1y)]^2+w}.$$
(33)
The first term has a restricted momentum integration due to the cut-off provided by the $`\mathrm{\Theta }`$-function. The second term is more involved, and the integrand even contains products of distributions. Notice, however, that it contains the factor $`(\frac{1}{\sqrt{y}}1)\delta (1y)`$ which is proportional to $`(y1)\delta (1y)`$. Therefore, the second term vanishes identically and independent of the parametrisation of the distributions and their products. The evaluation of the first term gives finally

$$\mathrm{}_{F,0}^d(\omega )=\frac{2}{d}\frac{1}{1+\omega },$$
(34)
and is identical to the bosonic flow.

### D Discussion

The flows described by the functions (32) and (34) have the simplest asymptotic structure for $`\omega \mathrm{}`$. This implies that heavy modes decouple “the fastest” from the flow for optimised regulators. For comparison, the sharp cut-off leads only to a logarithmic decoupling $`\mathrm{ln}\omega `$. Also, the decoupling does not depend on the particular theory studied (i.e. the dimension), unlike the case for polynomial regulators $`R_kq^2(k^2/q^2)^b`$. Furthermore, the flow described by the functions (32) and (34) have the simplest and strongest pole structure for $`C+\omega 0^+`$. The pole is a simple analytic one, which is not the case for generic regulator functions. An immediate implication of this structure is that the optimised flows (32) and (34) lead to a logarithmically smooth approach towards a convex effective potential. This is very different from the sharp cut-off case, where the approach is only exponentially . A detailed presentation of these results is given elsewhere.

For completeness we quote the results for the flows (32) and (34) for arbitrary normalisation. While the normalisation is of no relevance for a theory containing only bosonic or fermionic degrees of freedom, their relative normalisation can become important for theories containing bosons and fermions. The normalisation conditions (7) and (16) correspond to $`r(c_B)=1`$ and $`r_F(c_F)=1`$, which can always be imposed because the functions $`r(y)`$ and $`r_F(y)`$ range between $`0r,r_F\mathrm{}`$. The optimised gaps are $`C_{\mathrm{opt}}=2c_B`$ and $`C_{F\mathrm{opt}}=4c_F`$. For arbitrary $`c_B`$ the flow Eq. (32) is obtained as

$$\mathrm{}_0^d(\omega )=\frac{2}{d}\frac{(2c_B)^{{\scriptscriptstyle \frac{d}{2}}+1}}{2c_B+\omega }$$
(35)
In the fermionic case we find

$$\mathrm{}_{F,0}^d(\omega )=\frac{2}{d}\frac{(4c_F)^{{\scriptscriptstyle \frac{d}{2}}+1}}{4c_F+\omega }$$
(36)
for the rescaled analogue of Eq. (34).

Finally, we note that flows $`\mathrm{}_0^d(\omega )\frac{1}{1+\omega }`$ have been used earlier in the literature , however, without the explicit knowledge of the corresponding regulator function. These trial functions are sufficiently simple to allow for analytical considerations. The motivation for their use was based on the observation that the generic threshold function Eq. (28) decays at most as $`\omega ^1`$ for large $`\omega `$. This suggested that a regulator may exist which leads to $`\mathrm{}_0^d(\omega )=A_d(C+\omega )^1`$. Let us show how the normalisation $`A_d`$ can be fixed from consistency arguments. We use the universal relation $`\mathrm{}_n^{2n}(0)=1`$ , which holds for $`d=2n`$ dimensions, to identify the prefactor as $`A_{2n}=\frac{1}{n}C^{n+1}`$. The analytic continuation to arbitrary dimensions leads finally to our results (35) and (36). This reasoning shows that the Ansatz $`\mathrm{}_0^d(\omega )=\frac{2}{d}C^{d/2+1}(C+\omega )^1`$ is self-consistent. However, we rush to add that these consistency arguments are necessary conditions, but not sufficient ones: only the explicit form of the regulator – as given by Eqs. (20) and (22) – finally justifies the few earlier computations. In addition, Eqs. (20) and (22) are explicitly required for the computation of the flow at finite temperature (see Section IV) or to higher order (see Section V).

## IV Thermal fluctuations

In this section we apply our reasoning in the context of a quantum field theory coupled to a heat bath at temperature $`T`$, and to leading order in the derivative expansion. We show that optimised flows, as opposed to generic ones, disentangle the different contributions related to thermal and quantum fluctuations, respectively. These properties are realised, on the level of the flow equation, in terms of an important factorisation. This leads to better convergence properties of the flow itself. Approximate solutions of the flow correspond to better approximations of the physical theory.

### A Imaginary time formalism

To be explicit, we consider a bosonic or fermionic field theory at thermal equilibrium at the temperature $`T`$ within the Matsubara formalism. This implies that periodic (resp. antiperiodic) boundary conditions for the bosonic (resp. fermionic) fields are employed. As a consequence, the $`q_0`$ integration in the flow equation (1) is replaced by a sum over Matsubara modes $`m=0,\pm 1,\pm 2,\mathrm{}`$. The trace in Eq. (1) contains a momentum integration, which is then substituted as

$$\frac{d^dq}{(2\pi )^d}T\underset{m}{}\frac{d^{d1}q}{(2\pi )^{d1}}.$$
(37)
In the integrand of Eq. (1) the $`q_0`$ variable is replaced by

$$q_02\pi c_mT$$
(38)
where

$`c_m`$ $`=`$ $`m\mathrm{for}\mathrm{bosons}`$ (39)
$`c_m`$ $`=`$ $`m+\frac{1}{2}\mathrm{for}\mathrm{fermions}.`$ (40)
It is also useful to introduce the variable

$$\tau =2\pi T/k$$
(41)
for the following considerations. The replacement (38) implies that functions $`\mathrm{}_n^d(\omega )`$ turn into temperature dependent functions $`\mathrm{}_n^d(\omega ,\tau )`$. We show that this function factorises for the regulators (20) and (22).

### B Dimensional reduction and fermion decoupling

Let us review a few basic facts known for generic flows at finite temperature within the imaginary time formalism .

Bosonic fields within the Matsubara formalism display the phenomenon of dimensional reduction at high temperature. This means that for $`T`$ large enough all non-vanishing Matsubara modes are suppressed due to effective masses $`mT`$ for the Matsubara modes with $`m0`$. Only the $`m=0`$ mode survives in this limit, leading to an effective theory in $`(d1)`$-dimensions. For a generic bosonic regulator, the finite temperature flow is given as

$$\mathrm{}_0^d(\omega ,\tau )=\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{m}{}_0^{\mathrm{}}𝑑yy^{{\scriptscriptstyle \frac{d1}{2}}1}\frac{(y+c_m^2\tau ^2)^2r^{}(y+c_m^2\tau ^2)}{P^2(y+c_m^2\tau ^2)+\omega }.$$
(42)
The function $`P^2`$ is defined in Eq. (11). The asymptotic regime where only the $`m=0`$ Matsubara mode contributes is reached for $`\tau \mathrm{}`$. From Eqs. (42) and (29), we deduce

$$\mathrm{}_n^d(\omega ,\tau \mathrm{})=\frac{v_{d1}}{v_d}\frac{T}{k}\mathrm{}_n^{d1}(\omega ).$$
(43)
On the other hand, the limit $`\tau 0`$ eventually switches on all higher order Matsubara modes. It is straightforward to verify that

$$\mathrm{}_n^d(\omega ,\tau 0)=\mathrm{}_n^d(\omega ).$$
(44)
The asymptotic limits for $`T\mathrm{}`$, Eq. (43), and $`T0`$, Eq. (44), display dimensional reduction for bosons as a function of temperature for generic regulator function .

Fermions at finite temperature within the Masubara formalism can be treated in essentially the same way. However, they happen to have no $`m=0`$ mode as antiperiodic boundary conditions have to be used on the $`q_0`$-integration. Hence, fermions do not display dimension reduction. Rather, they decouple completely from the RG flow once the smallest Matsubara mode is larger than the scale $`k`$. These properties can be read-off from the temperature-dependent flow. For a generic fermionic regulator, the flow $`\mathrm{}_{F,0}^d(\omega )`$ at finite temperature is defined as

$$\mathrm{}_{F,0}^d(\omega ,\tau )=\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{m}{}_0^{\mathrm{}}𝑑yy^{{\scriptscriptstyle \frac{d1}{2}}1}\frac{2(y+c_m^2\tau ^2)^2r_F^{}(y+c_m^2\tau ^2)[1+r_F(y+c_m^2\tau ^2)]}{P_F^2(y+c_m^2\tau ^2)+\omega }.$$
(45)
The function $`P_F^2`$ is given in Eq. (17). The asymptotic regime where the fermions decouple completely is reached for $`\tau \mathrm{}`$. From Eq. (45), we deduce that

$$\mathrm{}_{F,n}^d(\omega ,\tau \mathrm{})=0.$$
(46)
Again, the limit $`\tau 0`$ eventually switches on all higher order Matsubara modes such that

$$\mathrm{}_{F,n}^d(\omega ,\tau 0)=\mathrm{}_{F,n}^d(\omega ).$$
(47)
The asymptotic limits (46) and (47) describe the decoupling of fermions in the high temperature limit for arbitrary dimension and generic regulator function.

### C Optimised thermal flows and factorisation

We now turn to the optimised regulators (20) and (22). For this case, the flow (42) can be computed explicitly. Inserting Eq. (20) into Eq. (42), and following a reasoning analogous to the one after Eq. (31), we find

$$\mathrm{}_n^d(\omega ,\tau )=B_d(\tau )\mathrm{}_n^d(\omega )$$
(48)
with the temperature dependent function

$$B_d(\tau )=\frac{d}{d1}\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{m}{}\left(1c_m^2\tau ^2\right)^{{\scriptscriptstyle \frac{d1}{2}}}\mathrm{\Theta }\left(1c_m^2\tau ^2\right).$$
(49)
Notice that the temperature effects have factorised. This implies that temperature cuts off all amplitudes $`\omega `$ in the same manner. This is not the case for a generic regulator.

Let us discuss the thermal threshold factor $`B_d(\tau )`$. In Fig. 2 the thermal threshold factor $`B_d(\tau )`$ is displayed for $`d=2,3`$ and $`4`$ dimensions. Every single Matsubara mode contributes to Eq. (49) proportional to

$$\tau \left(1c_m^2\tau ^2\right)^{{\scriptscriptstyle \frac{d1}{2}}}\mathrm{\Theta }\left(1c_m^2\tau ^2\right).$$
(50)
The $`\mathrm{\Theta }`$-function is a remnant of the regulator (21) and cuts the $`m`$-th Matsubara mode off as soon as $`k<c_mT/2\pi `$. The factor $`\tau `$ stems from the $`q_0`$-integration and the factor $`(1c_m^2\tau ^2)^{(d1)/2}`$ from the $`d1`$ dimensional integration over spatial loop momenta $`|𝐪|`$. These functions vanish outside the interval $`0\tau 1/c_m`$. At the upper end they behave like $`(1/c_m\tau )^{(d1)/2}`$ and vanish linearly with $`\tau `$ at the lower end. This structure explains the spikes observed in Fig. 2, which are located precisely at the points $`\tau =1/c_m`$ and due to the decoupling of the $`\pm c_m`$-th Matsubara modes. Indeed, for $`\tau >1`$ only the $`m=0`$ Matsubara mode yields a contribution to $`B_d(\tau )`$ in Eq. (49). The asymptotic regime where only the $`m=0`$ Matsubara mode contributes is reached already for $`\tau >1`$ with $`B_d(\tau 1)=\tau dv_{d1}/2\pi v_d(d1)`$, or

$$\mathrm{}_n^d(\omega ,\tau 1)=\frac{v_{d1}}{v_d}\frac{T}{k}\mathrm{}_n^{d1}(\omega ).$$
(51)
Notice the difference to Eq. (43). Decreasing $`\tau `$ below $`\tau =1/c_m`$ eventually switches on the $`\pm c_m`$ Matsubara modes. For $`\tau `$ close to the points $`1/c_m`$, the term (50) increases as $`(1/c_m\tau )^{(d1)/2}`$ for decreasing $`\tau `$. This power law explains why the spikes are more pronounced in lower dimensions.To higher order in the derivative expansion, the spikes are smoothed-out for non-trivial wave function renormalisation, cf. Sect.V C. In the limit $`\tau 0`$ it is straightforward to verify that $`B_d(\tau 0)1`$ which implies

$$\mathrm{}_n^d(\omega ,\tau 0)=\mathrm{}_n^d(\omega ).$$
(52)
This asymptotic limit is the same as Eq. (44).

Turning to the optimal fermionic regulator (23), the momentum integration in Eq. (45) can be performed explicitly to give

$$\mathrm{}_{F,n}^d(\omega ,\tau )=F_d(\tau )\mathrm{}_{F,n}^d(\omega ).$$
(53)
As in the bosonic case, the temperature effects factorise from the threshold effects. The fermionic thermal threshold factor $`F_d(\tau )`$ is given by

$$F_d(\tau )=\frac{d}{d1}\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{m}{}\left(1c_m^2\tau ^2\right)^{{\scriptscriptstyle \frac{d1}{2}}}\mathrm{\Theta }\left(1c_m^2\tau ^2\right)$$
(54)
Eq. (54) is identical to its bosonic counterpart (49) except for the Matsubara sum which runs over $`c_m=\pm \frac{1}{2},\pm \frac{3}{2},\mathrm{}`$ in Eq. (54). In Fig. 3 we have displayed the function $`F_d(\tau )`$ for $`d=2,3`$ and $`4`$ dimensions. Again, the spikes have the same origin as in the bosonic case and the same reasoning applies. The high temperature limit at which the fermions decouple completely, is already reached for $`k\pi T`$,

$$\mathrm{}_{F,n}^d(\omega ,\tau 2)=0.$$
(55)
Notice the important difference to Eq. (46), where the decoupling of fermions is only asymptotic. The limit $`\tau 0`$ is equivalent to Eq. (47).

### D Discussion

The optimised regulators (20) and (22) correctly describe dimensional reduction and fermion decoupling. In addition, they lead to a thermal factorisation of the flow as observed in Eqs. (48) and (53). From a physical point of view, this fact is easily understood. The imaginary time formalism compactifies the time direction and the temperature modifies the temporal momentum modes of the fields. The corresponding Matsubara mode, when compared to the infrared scale $`k`$, leads to a thermal decoupling. To leading order in the derivative expansion, the optimised regulator makes the temperature blind for the quantum fluctuations. It cannot distinguish between amplitudes of constant fields. In turn, the quantum fluctuations are sensitive to the field amplitudes, which are responsible for the mass decoupling, similar to the case for vanishing temperature. Clearly, these two effects have different physical origins. As a consequence, it is natural to employ a regulator which reflects this factorisation on the level of the flow equation.

For a generic regulator function, the flows (42) and (45) are complicated functions of both the temperature and the field amplitudes. They reflect dimensional reduction and fermion decoupling. Typically, however, they do not factorise. This simply means that a generic ERG flow entangles thermal and quantum fluctuations even to leading order in the derivative expansion. This is a direct consequence of the regulator term, whose coupling to the different operators in the effective action is leading to a field-dependent thermal decoupling of the different modes on the level of the flow equation. This entanglement is of no relevance if the flow can be solved exactly. In turn, for an approximate solution of the flow, the factorisation on the level of the flow equation is most helpful. It avoids a mixing of thermal and quantum fluctuations in a regime where they can be disentangled, thereby minimising possible artifacts due to the specific regulator function. As a consequence, the flow itself is stabilised, and expansions of the flow show much better convergence behaviour towards the physical theory. More generally, it is expected that this line of reasoning applies for generic optimised regulators.

Finally, the factorisation is very helpful for numerical solutions of flow equations. In the generic case, one two-parameter function has to be fitted in order to describe the flows (42) or (45). In turn, only two one-parameter functions are needed once they factorise as in Eqs. (48) and (53). This simplification is substantial, and even more so, because the functions (48) and (53) have a very simple analytical form.

## V Derivative expansion to second order

In this section, we apply the optimised regulator to higher order in the derivative expansion. We first discuss the general structure of the equations. Furthermore, we show that a simple factorisation of the flow takes place for field-independent wave function renormalisation. The physical origin of the factorisation is discussed, and its realisation on the level of the flow equation leads, as in the thermal case, to better convergence properties of the flow and the derivative expansion. For technical details on the computations, we defer to the Appendices .

### A Wave-function renormalisation

In the preceding sections we have restricted the discussion to the leading order in a derivative expansion. This implied the vanishing of the anomalous dimensions $`\eta _\varphi =\eta _\psi =0`$ or $`Z_\varphi Z_\psi 1`$. To higher order in the derivative expansion the multiplicative renormalisation of the fields has to be taken into account, according to

$`\varphi `$ $``$ $`Z_{\varphi ,k}^{1/2}\varphi `$ (56)
$`\psi `$ $``$ $`Z_{\psi ,k}^{1/2}\psi .`$ (57)
Once higher derivative terms are included in the Ansatz for the effective action, additional flow equations for the corresponding coefficient functions like $`Z_{\varphi ,k}(q^2,\overline{\rho })`$ and $`Z_{\psi ,k}(q^2,\overline{\rho })`$ have to be studied. The wave function renormalisations are functions of the scale parameter $`k`$ and can depend as well on momenta $`q^2`$ or on the mean fields $`\overline{\rho }`$. To second order in the derivative expansion the wave function renormalisation is evaluated at a particular momentum scale $`q^2=k_0^2`$ which fixes the renormalisation conditions. Typical choices for $`k_0`$ are $`k_0=0`$ and $`k_0=k`$.

The most important new ingredient at this order is the scale- and field-dependence of the wave function renormalisations. In the example defined through Eq. (24), these are given by the functions $`Z_k`$ and $`Y_k`$. Here, the function $`Z_k`$ is responsible for the renormalisation of the $`N1`$ “Goldstone” modes, which differs from $`\stackrel{~}{Z}_k=Z_k+\overline{\rho }Y_k`$ for the “radial” mode. The fact that different wave function renormalisations appear to second order in the derivative expansion (depending of the theory considered) is of no relevance for the following discussion of the flows. The parametric dependence of the flow on either $`Z_k`$ or $`\stackrel{~}{Z}_k`$ is the same.

Let us introduce an additional function $`z_k(\rho )`$ as

$$Z_k(\overline{\rho })=Z_k(\overline{\rho }_0)z_k(\overline{\rho }).$$
(58)
We have factored out a constant term $`Z_k(\overline{\rho }_0)`$ chosen at an arbitrary reference point. We have chosen the reference point $`\overline{\rho }=\overline{\rho }_0`$ which fixes the renormalisation of the fields for all momenta. Typical choices are either $`\rho _0=0`$, or $`\rho _0=`$ the minimum of the scale-dependent potential. The split (58) allows to separate the non-trivial field-dependence, contained in $`z_k(\overline{\rho })`$, from an overall renormalisation contained in $`Z_k(\overline{\rho }_0)`$. The factor $`z_k`$ is normalised as $`z_k(\overline{\rho }_0)=1`$.

In order to provide a simple form for the flow it is useful to introduce the field-independent factor $`Z_k(\overline{\rho }_0)`$ into the regulator function,

$`R_k`$ $``$ $`Z_{\varphi ,k}R_k`$ (59)
$`R_{F,k}`$ $``$ $`Z_{\psi ,k}R_{F,k}.`$ (60)
The flow equation, when written in terms of renormalised variables (56) and (57), receives additional contributions proportional to the anomalous dimensions

$`\eta _\varphi `$ $`=`$ $`_t\mathrm{ln}Z_{\varphi ,k},`$ (61)
$`\eta _\psi `$ $`=`$ $`_t\mathrm{ln}Z_{\psi ,k},`$ (62)
because the derivative $`_tR_k`$ in the flow equation now acts also on the explicit scale-dependence contained in $`Z_k`$.

### B Generic flows to second order

To second order in the derivative expansion the flow has turned into a function of the field amplitudes $`\omega `$, the anomalous dimension $`\eta `$ and the field dependent function $`z`$. The corresponding bosonic flow is defined as

$$\mathrm{}_0^d(\omega ,z,\eta )=_0^{\mathrm{}}𝑑yy^{{\scriptscriptstyle \frac{d}{2}}1}\frac{y^2r^{}(y)\frac{1}{2}\eta yr(y)}{y[z+r(y)]+\omega }.$$
(63)
Notice that the pole structure of the flow is changed, owing to the function $`z(\rho )`$. The effective inverse propagator becomes a function of the fields:

$$P^2(\rho ,y)=y[z(\rho )+r(y)].$$
(64)
The location of the pole of Eq. (63) at $`\omega =C(\rho )=\mathrm{min}_{y0}P^2(\rho ,y)`$ has turned into a function of the fields. For the optimised regulator (21) the pole is located at $`C=\mathrm{min}\{1,z\}`$. Compared to the leading order in the derivative expansion, the pole structure is modified once $`z1`$.

For the fermionic case, the flows are given as

$$\mathrm{}_{F,0}^d(\omega ,z_F,\eta _\psi )=_0^{\mathrm{}}𝑑yy^{d/2}\frac{(2yr_F^{}+\eta _\psi r_F)(z_F+r_F)}{y[z_F+r_F(y)]^2+\omega }.$$
(65)
The pole structure changed as well, as follows from

$$P_F^2(\rho ,y)=y[z_F(\rho )+r_F(y)]^2.$$
(66)
The location of the pole of Eq. (63) at $`\omega =C_F(\rho )=\mathrm{min}_{y0}P_F^2(\rho ,y)`$ has turned into a function of the fields. For the regulator (23) the pole is located at $`C_F=\mathrm{min}\{1,z_F\}`$. The pole structure is modified compared to the leading order in the derivative expansion once $`z_F1`$.

### C Optimised flows to second order

We now turn to the optimised flows and discuss their structure at second order in the derivative expansion. We refer to the Appendices for all technical details.

In the case of a generic wave-function renormalisation where $`z_k(\overline{\rho })`$ is a non-trivial function of the fields, the function $`\mathrm{}_0^d(\omega ,z,\eta )`$ as defined in Eq. (63) can be evaluated explicitly for the regulator (20). The structure of the flow is as follows. Consider the denominator of Eq. (63), given by $`y(z+r)+\omega `$. It can be rewritten as $`1+\omega +y(z1)+(y(1+r)1)`$. The last term in brackets vanishes for the optimised regulator (20) because the integration is restricted to $`y1`$. The remaining term can be written as the product $`(1+\omega )(1\frac{1z}{1+\omega }y)`$. Notice also that the numerator of Eq. (63) neither depends on $`z`$ nor on $`\omega `$. These observations lead to the following conclusions. First, and apart from an overall $`\omega `$-dependence $`(1+\omega )^1`$, the optimised flow depends on $`z`$ and $`\omega `$ only through the variable

$$\xi \frac{1z}{1+\omega }.$$
(67)
Second, the optimised flow factorises into a leading order term (32) and a remaining factor $`B_d(\xi ,\eta )`$,

$$\mathrm{}_0^d(\omega ,\xi ,\eta )=\mathrm{}_0^d(\omega )B_d(\xi ,\eta ),$$
(68)
because the denominator of Eq. (63) contains a momentum-independent factor $`(1+\omega )`$. Here we have introduced $`\mathrm{}_0^d(\omega ,\xi (\omega ,z),\eta )\mathrm{}_0^d(\omega ,z,\eta )`$. It it interesting to see that the structure of the optimised flow is still quite simple. An integral representation of $`B_d(\xi ,\eta )`$ is given in Eq. (B.107). (For all $`\xi <1`$, $`B_d(\xi ,\eta )`$ can be expressed in terms of hypergeometric functions, cf. Eqs. (B.108), (B.109) and (B.110); closed expressions of Eq. (68) for $`d=4,3`$ and $`2`$ dimensions are given in Eqs. (B.112), (B.114) and (B.115), respectively.) For $`|\xi |<1`$ the function $`B_d(\xi ,\eta )`$ can be Taylor-expanded in arbitrary dimensions, to wit

$$B_d(\xi ,\eta )=\underset{n=0}{\overset{\mathrm{}}{}}\frac{\xi ^n}{1+\frac{2n}{d}}\left(1\frac{\eta }{d+2+2n}\right)=1\frac{\eta }{d+2}+\frac{\xi }{1+\frac{2}{d}}+𝒪(\xi ^2,\xi \eta ).$$
(69)
The series representation (69) is best suited for the flow as long as $`|\xi |`$ remains small. This corresponds to either the limit of a field-independent wave function renormalisation $`z(\rho )1`$, or, for any $`z`$, to the limit of large amplitudes $`\omega `$. From Eq. (69), we obtain for Eq. (68) to zeroth order in $`\xi `$

$$\mathrm{}_0^d(\omega ,\xi ,\eta )=\mathrm{}_0^d(\omega )\left(1\frac{\eta }{d+2}\right)+𝒪(\frac{\xi }{1+\omega },\frac{\xi \eta }{1+\omega }).$$
(70)
We note that the $`\eta `$-dependent correction in Eq. (70) has the same functional dependence on the amplitude $`\omega `$ as Eq. (32). Stated differently, the optimised regulator leads to a simple factorisation in both the decoupling limit $`\omega 1`$ and for the case of a field-independent anomalous dimension $`\xi 0`$.

For completeness we give also the result for the bosonic flow at finite temperature. The corresponding flow $`\mathrm{}_0^d(\omega ,\xi ,\tau ,\eta )`$ still factorises as

$$\mathrm{}_0^d(\omega ,\xi ,\tau ,\eta )=B_d(\xi ,\tau ,\eta )\mathrm{}_0^d(\omega ).$$
(71)
This is the generalisation of Eq. (48) to second order in the derivative expansion. We only have to replace the function $`B_d(\xi ,\eta )`$ by its temperature-dependent counterpart $`B_d(\xi ,\tau ,\eta )`$ \[cf. Eq. (A.89)\]. It is straightforward, if tedious, to establish explicitly that the function $`B_d(\xi ,\tau ,\eta )`$ represents dimensional reduction in precisely the same way as $`B_d(\tau )B_d(\xi =0,\tau ,\eta =0)`$. Let us consider the most interesting case, which is the leading order in $`\xi 1`$. In this limit, Eq. (71) reads

$$B_d(\xi ,\tau ,\eta )=B_d(\tau ,\eta )+𝒪(\xi ,\xi \eta ).$$
(72)
The function $`B_d(\tau ,\eta )`$ can be expressed as

$$B_d(\tau ,\eta )=B_d\left(\tau \right)\eta \frac{d}{d^21}\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{m}{}\left(1c_m^2\tau ^2\right)^{{\scriptscriptstyle \frac{d+1}{2}}}\mathrm{\Theta }\left(1c_m^2\tau ^2\right).$$
(73)
The function $`B_d(\tau )`$, Eq. (49), has been discussed in Section IV. The new ingredient, beyond leading order, are given by the corrections $`\eta `$ in Eq. (73). Every single Matsubara mode contributes as

$$\tau \left(1c_m^2\tau ^2\right)^{{\scriptscriptstyle \frac{d+1}{2}}}\mathrm{\Theta }\left(1c_m^2\tau ^2\right).$$
(74)
Compared to the leading order contributions (50), we notice that Eq. (74) follows from Eq. (50) for $`dd+2`$. The reason is very simple. In the flow equation, the anomalous dimension is proportional to a term containing an additional factor $`q^2`$, which effectively increases the momentum measure by two dimensions. This has an immediate consequence. The thermal decoupling in Eq. (73) proportional to $`\eta `$ is much smoother than the leading order decoupling, simply because the spikes are less pronounced the higher the dimension. Therefore, the spikes observed in Fig. 2 are smoothed-out once $`\eta `$ (and $`\xi `$) are non-vanishing.

In the opposite regime where $`\frac{|1z|}{1+w}1`$, only a few leading terms of the series (69) have to be retained. This limit is of relevance close to the pole region of the flows $`\omega 1`$, or in the region of large $`z1`$. From the explicitly resummed expressions (B.112), (B.114) and (B.115), we conclude that a factorisation as

$$\mathrm{}_0^d(\omega ,z,\eta )=f_d(\omega ,z)\left(1\frac{\eta }{d}\right)$$
(75)
holds true, and $`f_4(z)=(z1)^1`$, $`f_3(z)=2f_4(z)`$ and $`f_2(\omega ,z)=f_4(z)\mathrm{ln}(\frac{z+\omega }{1+\omega })`$.

It is not surprising that a similar structure is found for fermionic flows. The correction term due to the substitution Eq. (59) simplifies Eq. (65) to

$$\mathrm{}_{F,0}^d(\omega ,z_F,\eta _\psi )=_0^1𝑑yy^{d/21}\frac{[1+\sqrt{y}(z_F1)][1\eta _\psi (1\sqrt{y})]}{[1+\sqrt{y}(z_F1)]^2+w}.$$
(76)
Eq. (76) factorises as

$$\mathrm{}_{F,0}^d(\omega ,z_F,\eta _\psi )=\mathrm{}_{F,0}^d(\omega )F_d(\omega ,z_F,\eta _\psi ).$$
(77)
The function $`F_d(\omega ,z,\eta )`$ can be expressed in terms of hypergeometric functions. At finite temperature, and for $`z=1`$ and $`\eta =0`$, it reduces to Eq. (54). Here, we are only interested in the structure of the flow for a nearly field-independent wave function renormalisation, $`z1`$, or for the decoupling limit. We find

$$F_d(\omega ,z_F,\eta _\psi )=1\frac{\eta _\psi }{d+1}+\frac{d(z_F1)}{d+1}\left(1\frac{2}{1+\omega }\right)\left(1\frac{\eta _\psi }{d+2}\right)+\mathrm{}$$
(78)
The two leading terms in Eq. (78) show that also fermionic flows factorise for field-independent wave function renormalisation.

### D Discussion

The structure of the flow has increased to second order in the derivative expansion. Let us discuss first the case of a field-independent wave function renormalisation $`z1`$. The corresponding flows (70) and (78) for the optimised regulators factorise, similar to the thermal case to leading order in the derivative expansion. Physically speaking, this structure can be made plausible as follows. The flow, when written in terms of the renormalised fields – and under the assumption that the renormalisation is momentum- and field-independent – depends, in addition to the fields, only on the anomalous dimension. The anomalous dimension is field independent, and, in consequence, unable to distinguish between fields of different amplitudes contained in $`\omega `$, which parametrise the quantum fluctuations. Therefore, it is natural that the flow factorises the contributions induced through $`\eta `$ from those induced by the amplitudes $`\omega `$. The disentanglement is realised by the optimised regulators.From the definition of Eqs. (63) and (65) it follows that all homogeneous regulators with $`r(y)yr^{}(y)`$ (or $`r_F(y)yr_F^{}(y)`$, respectively) factorise the anomalous dimension from the field-dependent part of the flow.

In turn, a generic flow does not reflect this factorisation. Rather, it leads to an entanglement between the renormalisation of the effective potential induced by the infrared regulator, and the renormalisation parametrised by a field-independent anomalous dimension. This is immediately evident from the observation that the $`\eta `$-dependent and the $`\eta `$-independent contributions to the flow of the effective potential have different functional forms as functions of the fields. At this level, the entanglement is due to the regulator, which modifies the coupling amongst all operators of the effective action. As mentioned in the thermal case, the entanglement is of no importance for the full solution to the flow. In turn, the factorisation is very useful for approximate solutions. It leads to more stable flows because irrelevant couplings, entirely due to the regulator, are removed. The same reasoning as given at the end of the previous section applies.

For the thermal bosonic flow (71), we notice that the dependence on the anomalous dimension enters the thermal factor $`B_d(\tau ,\eta )`$. In particular, the thermal corrections do not factorise from those due to a field-independent anomalous dimension. This structure can be understood as follows. The wave function renormalisation enters the momentum trace as a multiplicative renormalisation proportional to the kinetic term $`q^2`$. At finite temperature within the imaginary time formalism, the spatial and the temporal loop momenta are treated in an unequal way. Hence, thermal fluctuations couple in a non-trivial manner to the anomalous dimension of the fields. This implies that the temperature-dependent factor itself is modified due to the anomalous dimension, which provides the physical reason way no factorisation of the temperature effects from the anomalous dimension are expected in the first place.

For the case of a field-dependent wave-function renormalisation, a simple factorisation similar to Eq. (70) is not expected, simply because the wave function renormalisation is a function of the fields. Hence, the wave function renormalisation can distinguish different field amplitudes, in contrast to the field-independent case. However, two observations are still worth mentioning. First of all, we observe a partial factorisation, which is evident from Eqs. (68) and (71). This structure is based on the fact that the $`z`$-dependence enters only through the variable (67), as opposed to the generic case. For $`\frac{|1z|}{1+w}1`$, only a few leading terms have to be retained from the explicit series (69). It follows that each power of $`\frac{|1z|}{1+w}`$ is renormalised proportional to the anomalous dimension and an order-dependent numerical coefficient. Secondly, the limit for $`\frac{|1z|}{1+w}1`$ again allows for a simple factorisation, as follows from Eq. (75). Here, the wave function renormalisation can no longer distinguish field amplitudes, allowing for this simple structure.

A final comment concerns the numerical prefactors $`\eta `$ as found in Eqs. (70) and (78). We emphasise that the coupling of the anomalous dimensions to the effective potential is, apart from the field dependence, dimensionally suppressed – by factors $`\frac{1}{d+2}`$ for bosons and $`\frac{1}{d+1}`$ for fermions – as opposed to the leading order contributions. This additional suppression is noteworthy because the convergence of the derivative expansion is controlled by small anomalous dimensions of the fields. Here, we have just shown that an expansion performed with an optimised regulator leads to an additional dimensional suppression of the back-coupling of the anomalous dimension to the effective potential. A more detailed discussion of this observation will be given elsewhere.

## VI  Proper-time regularisation

In this section we leave aside the conceptual framework of the ERG based on a momentum-scale regularisation and address flows based on an operator cut-off regularisation. Our aim is to provide the analog of the optimised regulator (20) within the proper-time regularisation method. For a more detailed comparison with the exact renormalisation group, we refer the reader to Ref. .

A simple flow has been derived from a one-loop expression for the effective action which is UV and IR regularised using a Schwinger proper-time representation of the operator trace , amended by a regulator function $`f_k^{(d)}(\mathrm{\Lambda },s)`$ within the proper-time integral . The flow with respect to the infrared scale parameter $`k`$ follows from a 1-loop improvement as

$$_t\mathrm{\Gamma }_k=\frac{1}{2}_0^{\mathrm{}}\frac{ds}{s}\left(_tf_k^{(d)}(\mathrm{\Lambda },s)\right)\mathrm{Tr}\mathrm{exp}\left(s\mathrm{\Gamma }_k^{(2)}\right).$$
(79)
We refer to this flow as the “proper-time renormalisation group” (PTRG). It describes the partial resummation of perturbative diagrams. The proper-time regulator function plays the role of the momentum regulator $`R_k`$ within the ERG. The flow (79) is governed by the IR scale $`k`$. Following Ref. , we introduce a dimensionless function $`f(x)`$ as $`f_k^{(d)}(\mathrm{\Lambda },s)=f(\mathrm{\Lambda }^2s)f(k^2s)`$ and require $`f(x\mathrm{})=1`$ and $`f(x0)=0`$. This ensures that the usual Schwinger proper time representation is reached in the UV limit.

We are not aware of a simple and generic optimisation criterion, analogous to Eq. (13), which derives from within the PTRG formalism. Furthermore, the flow (79) has no path integral derivation, which makes a conceptual reasoning much more difficult. However, it is still possible to show that a function $`f_{\mathrm{opt}}(x)`$ exists which is equivalent to the optimised ERG regulator (20) to the leading order in the derivative expansion.

To that end, we apply Eq. (79) to a $`N`$-component real scalar theory in $`d`$ dimensions and to leading order in the derivative expansion. Using the Ansatz (24) the flow for the effective potential $`U_k(\overline{\rho })`$ with $`\overline{\rho }=\frac{1}{2}\varphi _a\varphi _a`$ becomes

$$_tU_k(\overline{\rho })=\frac{1}{2}(4\pi )^{{\scriptscriptstyle \frac{d}{2}}}_0^{\mathrm{}}\frac{ds}{s^{1+d/2}}_tf_k^{(d)}(\mathrm{\Lambda },s)\left[e^{s\left(U_k^{}(\overline{\rho })+2\overline{\rho }U_k^{\prime \prime }(\overline{\rho })\right)}+(N1)e^{sU_k^{}(\overline{\rho })}\right].$$
(80)
This flow is identical in form to the ERG flow (26), if we replace the ERG flow in Eq. (26) by the proper-time flow

$$\mathrm{}_0^d(\omega )=\frac{1}{2}\mathrm{\Gamma }(\frac{d}{2})_0^{\mathrm{}}𝑑xx^{1{\scriptscriptstyle \frac{d}{2}}}(_tf(x))\mathrm{exp}(x\omega ).$$
(81)
Here, the integration variable is $`x=k^2s`$ and stems from the proper-time integration, in contrast to Eq. (28), where $`y=q^2/k^2`$ stems from the momentum trace. Now, consider a specific class of proper-time regulator functions:

$$f(x)=\frac{\mathrm{\Gamma }(m,x)}{\mathrm{\Gamma }(m)},_tf(x)=\frac{2x^me^x}{\mathrm{\Gamma }(m)}$$
(82)
We have introduced a free parameter $`m`$ describing different regulators, and the incomplete $`\mathrm{\Gamma }`$-function $`\mathrm{\Gamma }(m,x)=_0^x𝑑yy^{m1}e^y`$. This yields the simple expression

$$\mathrm{}_0^d(\omega )=\frac{\mathrm{\Gamma }(m\frac{d}{2})\mathrm{\Gamma }(\frac{d}{2})}{\mathrm{\Gamma }(m)}(1+\omega )^{{\scriptscriptstyle \frac{d}{2}}m}$$
(83)
which agrees with Eq. (32) for $`m=1+\frac{d}{2}`$, or

$$f^{\mathrm{opt}}(x)=\frac{\mathrm{\Gamma }(\frac{d}{2}+1,x)}{\mathrm{\Gamma }(\frac{d}{2}+1)},_tf^{\mathrm{opt}}(x)=\frac{2x^{1+d/2}e^x}{\mathrm{\Gamma }(\frac{d}{2}+1)}.$$
(84)
The optimised proper-time regulator (84) corresponds to the optimised regulator (20) within the ERG approach. Hence, it is possible to identify an optimal regulator function for proper-time flows, owing to their close similarity to the ERG to leading order in the derivative expansion.

## VII Conclusions and outlook

This study was motivated by two observations. First, an application of the ERG to realistic physical problems is bound to certain approximations. Second, approximate solutions of flow equations depend spuriously on the infrared regulator. Combining these observations, it became obvious that an understanding of the spurious scheme dependence is mandatory in order to provide predictive power for approximate solutions. Previously, we showed that the gap of the full inverse propagator controls convergence properties of approximate solutions . It has also been shown, based on the computation of critical exponents for the Ising universality class, that the convergence of the derivative expansion is controlled by the gap . These observations lead to the conclusion that the freedom in the choice for the IR regulator can be used to maximise the physical information contained within a given approximation or truncation.

An interpretation of the interplay between the RS function on one side, and convergence of approximate flows on the other, is as follows. The IR regulator – by regulating the flow – modifies the interactions at intermediate scales $`k0`$ amongst all operators of the theory. Eventually, these cancel out for the integrated full flow, but not for approximated ones. Hence, changing the RS function for an approximated flow modifies some remaining RS dependent terms which cannot be cancelled due to the missing contributions from neglected operators. Therefore, a “fine-tuning” of the RS function allows to partly incorporate higher-order effects within the lower orders of a given approximation. This corresponds to an optimisation.

The present derivation of optimised ERG flows had two ingredients. First, we made use of a generic optimisation criterion for bosonic and fermionic fields , which states that the gap of the full inverse propagator, as a function of momenta, should be as large as possible. This way, the ERG flow is the least singular, and approximations to such flows are expected to be much more stable, leading to improved convergence of expansions. Second, we added the specific requirement that the effective inverse propagator be momentum-independent in the IR region (cf. Fig. 1).

We have studied specific optimised ERG flows for bosonic or fermionic theories up to second order in the derivative expansion and at vanishing and non-vanishing temperature. Their specific properties have been discussed in detail at the end of the corresponding sections. Here, let us only mention the perhaps most surprising property of optimised flows, which is the disentanglement of thermal and quantum fluctuations to leading order in the derivative expansion. A similar factorisation occurs for field-independent wave function renormalisations.

More generally, optimised ERG flows owe their main properties to the “flatness” of the effective inverse propagator, which extends over the entire momentum region $`q^2k^2`$ for the specific regulator studied here. Other regulators can lead to similar factorisation and convergence properties. Prime candidates are given by solutions to the optimisation criterion: as is evident from Fig. 1, they automatically lead to flat effective propagators – at least within a small region about the minimum of the effective inverse propagator. If this region extends over the domain where the flow equation receives its main contributions, we expect to find equally good flows.

An important conclusion is that the optimisation ideas discussed here should be useful for high-precision computations based on this formalism. Increasing the precision normally implies a full computation at the following order of the expansion. Here, we argued that the physical results can be improved already within a fixed order of the expansion. Immediate applications of optimised flows concern the computation of universal critical exponents for $`N`$-component scalar theories in three dimensions, or the study of convex effective potentials within a phase with spontaneous symmetry breaking. On the conceptual side, it is possible to show that the optimisation criterion can be interpreted as a natural minimum sensitivity condition, somewhat similar to the principle of minimum sensitivity as employed within perturbative QCD. We will leave a detailed discussion of these results to a future publication .

Our analysis can be extended in a number of directions. For gauge theories, modified Ward or BRST identities ensure the gauge invariance of physical Green functions , and the optimisation criterion is compatible with such additional constraints. This optimisation can also be implemented for field theories at finite temperature within the real-time formalism . While our present analysis is based on the derivative expansion, it seems worthwhile to study optimisations for other systematic expansions like expansions in powers of the fields. Finally, it would be interesting to see how these ideas apply to flows at finite density or to Hamiltonian flows .

Acknowledgement: This work has been supported by a Marie-Curie fellowship under EC contract no. HPMF-CT-1999-00404.

## Appendix A Flows to second order in the derivative expansion

In this Appendix, we derive explicit expressions for the optimised flow for the effective potentials to second order in the derivative expansion at both vanishing and non-vanishing temperature.

Our starting point is the flow for the effective potential to second order in the derivative expansion. We consider the bosonic flow at finite temperature. All information about the flow is parametrised by

$$\mathrm{}_0^d(\omega ,z,\tau ,\eta )=\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{n}{}_0^{1c_n^2\tau ^2}𝑑yy^{{\scriptscriptstyle \frac{d3}{2}}}\frac{[1\frac{1}{2}\eta (1yc_n^2\tau ^2)]\mathrm{\Theta }(1c_n^2\tau ^2)}{(z1)(y+c_n^2\tau ^2)+1+\omega }$$
(A.85)
Here, $`\omega (\rho )`$ is the field variable, $`z(\rho )`$ a field dependent wave function renormalisation, $`\eta `$ the anomalous dimension, $`\tau =2\pi T/k`$ the rescaled dimensionless temperature, and $`c_n=n`$ the Matsubara modes in the bosonic case, Eq. (39). The constants $`v_d`$ are defined in Eq. (27).

The leading order behaviour of Eq. (A.85) is given by the function

$$\mathrm{}_0^d(\omega )=\frac{2}{d}\frac{1}{1+\omega },$$
(A.86)
which follows from Eq. (A.85) for $`z=1`$, $`\tau =0`$ and $`\eta =0`$. Factorizing the main building block Eq. (A.86) from Eq. (A.85), we notice that the remaining factor depends on both $`\omega `$ and the variable $`z`$ only through the combination

$$\xi \frac{1z}{1+\omega }.$$
(A.87)
Therefore, it is most natural to make the variable transform Eq. (A.87) by writing $`\mathrm{}_0^d(\omega ,z,\tau ,\eta )\mathrm{}_0^d(\omega ,\xi (z,\omega ),\tau ,\eta )`$, and to rewrite the flow (A.85) as

$$\mathrm{}_0^d(\omega ,\xi ,\tau ,\eta )=\mathrm{}_0^d(\omega )B_d(\xi ,\tau ,\eta ),$$
(A.88)
where

$$B_d(\xi ,\tau ,\eta )=\frac{d}{2}\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{n}{}_0^{1c_n^2\tau ^2}𝑑yy^{{\scriptscriptstyle \frac{d}{2}}{\scriptscriptstyle \frac{3}{2}}}\frac{1\frac{\eta }{2}(1yc_n^2\tau ^2)}{1\xi (y+c_n^2\tau ^2)}\mathrm{\Theta }(1c_n^2\tau ^2)$$
(A.89)
Below, if not stated otherwise, we adopt a simplified notation: functions are evaluated at the points $`\xi =0`$, $`z=1`$, $`\tau =0`$ or $`\eta =0`$, if the corresponding arguments are not displayed. With these definitions at hand, we can face the explicit computation of Eq. (A.85).

Let us compute the function $`B_d(\xi ,\tau ,\eta )`$ more explicitly. Since the anomalous dimension enters only linear in Eq. (A.85), it is helpful to rewrite Eq. (A.89) as

$$B_d(\xi ,\tau ,\eta )=B_d(\xi ,\tau )\eta \overline{B}_d(\xi ,\tau ).$$
(A.90)
For $`\xi <1`$, the remaining integration over the momentum variable in Eq. (A.89) can be performed. This leads to

$`B_d(\xi ,\tau )`$ $`=`$ $`{\displaystyle \frac{d}{d1}}{\displaystyle \frac{v_{d1}}{v_d}}{\displaystyle \frac{\tau }{2\pi }}{\displaystyle \underset{n}{}}A_n(\xi ,\tau )(1c_n^2\tau ^2)^{{\scriptscriptstyle \frac{d3}{2}}}\mathrm{\Theta }(1c_n^2\tau ^2)\times `$ (A.92)
$`{}_{2}{}^{}F_{1}^{}(1,\frac{d1}{2};\frac{d+1}{2};\xi A_n(\xi ,\tau ))`$
The function

$${}_{2}{}^{}F_{1}^{}(a,b;c;z)=\frac{\mathrm{\Gamma }(c)}{\mathrm{\Gamma }(b)\mathrm{\Gamma }(cb)}_0^1𝑑tt^{b1}(1t)^{cb1}(1tz)^a$$
(A.93)
with $`{}_{2}{}^{}F_{1}^{}(a,b;c;z)={}_{2}{}^{}F_{1}^{}(b,a;c;z)`$ denotes the hypergeometric function . We also introduced the thermal amplitude function

$$A_n(\xi ,\tau )=\frac{1c_n^2\tau ^2}{1\xi c_n^2\tau ^2}.$$
(A.94)
The index $`n`$ corresponds to the Matsubara mode. The factors $`A_n`$ only appear in combination with the factor $`\mathrm{\Theta }(1c_n^2\tau ^2)`$. At the limits,

$`A_n(\xi ,0)=A_0(\xi ,\tau )`$ $`=`$ $`1`$ (A.95)
$`A_n(\xi ,c_n^1)`$ $`=`$ $`0.`$ (A.96)
In the same way, we find from Eq. (A.89) for the term $`\eta `$ in Eq. (A.88) the explicit expression

$`\overline{B}_d(\xi ,\tau )`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{d}{d1}}{\displaystyle \frac{v_{d1}}{v_d}}{\displaystyle \frac{\tau }{2\pi }}{\displaystyle \underset{n}{}}A_n(\xi ,\tau )(1c_n^2\tau ^2)^{{\scriptscriptstyle \frac{d1}{2}}}\mathrm{\Theta }(1c_n^2\tau ^2)\times `$ (A.98)
$`\left[{}_{2}{}^{}F_{1}^{}(1,\frac{d1}{2};\frac{d+1}{2};\xi A_n(\xi ,\tau ))\frac{d1}{d+1}{}_{2}{}^{}F_{1}^{}(1,\frac{d+1}{2};\frac{d+3}{2};\xi A_n(\xi ,\tau ))\right]`$
Combining Eq. (A.92) with Eq. (A.98) gives Eq. (A.90) and hence Eq. (A.88) explicitly. Eq. (A.88) is the most general expression for the factorisation at second order in the derivative expansion and at finite temperature.

The temperature dependence of the function $`B_d(\xi ,\tau ,\eta )`$ describes dimensional reduction. In particular, it obeys the limits

$`B_d(\xi ,\tau =0,\eta )`$ $``$ $`B_d(\xi ,\eta )`$ (A.99)
$`B_d(\xi ,\tau 1,\eta )`$ $`=`$ $`{\displaystyle \frac{v_{d1}}{v_d}}{\displaystyle \frac{d}{d1}}{\displaystyle \frac{\tau }{2\pi }}B_{d1}(\xi ,\eta )`$ (A.100)
The low- and high-temperature limits (A.99) and (A.100) are discussed in the following appendices.

Let us consider the case where $`\xi =0`$. It corresponds either to the case of a field-independent wave function renormalisation $`z1`$, and/or to the decoupling regime $`\omega 1`$. A few properties of Eq. (A.88) have been discussed in the main text. For $`\xi =0`$ the factor (A.90) reduces to

$$B_d(\tau ,\eta )=\frac{d}{d1}\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{n}{}\left(1c_n^2\tau ^2\right)^{{\scriptscriptstyle \frac{d1}{2}}}\left(1\eta \frac{1c_n^2\tau ^2}{d+1}\right)\mathrm{\Theta }\left(1c_n^2\tau ^2\right).$$
(A.101)
This corresponds to Eq. (49) discussed in Section IV. In addition, we notice that

$`B_d(\tau ,\eta 0)`$ $`=`$ $`B_d(\tau )`$ (A.102)
$`B_d(\tau 0,\eta )`$ $`=`$ $`1{\displaystyle \frac{\eta }{d+2}}`$ (A.103)
$`B_d(\tau 1,\eta )`$ $`=`$ $`{\displaystyle \frac{d}{d1}}{\displaystyle \frac{v_{d1}}{v_d}}{\displaystyle \frac{\tau }{2\pi }}\left(1{\displaystyle \frac{\eta }{d+1}}\right).`$ (A.104)
Eq. (A.102) corresponds to Eq. (32), Eq. (A.103) to the low-temperature limit (70), and Eq. (A.104) to the high-temperature limit (51).

Similar results are found for the fermionic case, though not discussed explicitly.

## Appendix B Low temperature limit

In the low temperature limit $`\tau 0`$, the flow (A.85) simplifies to

$$\mathrm{}_0^d(\omega ,z,\eta )=_0^1𝑑yy^{{\scriptscriptstyle \frac{d2}{2}}}\frac{1\frac{1}{2}\eta (1y)}{(z1)y+1+\omega }.$$
(B.105)
The remaining integration in Eq. (B.105) is solved to give

$$\mathrm{}_0^d(\omega ,\xi ,\eta )=\mathrm{}_0^d(\omega )B_d(\xi ,\eta )$$
(B.106)
with

$$B_d(\xi ,\eta )=\frac{d}{2}_0^1𝑑yy^{{\scriptscriptstyle \frac{d}{2}}1}\left(\frac{1}{1\xi y}\frac{\eta }{2}\frac{1y}{1\xi y}\right)$$
(B.107)
in arbitrary dimensions. For $`\xi \frac{1z}{1+\omega }<1`$, and hence $`z>\omega `$, the integration can be peformed analytically. We find

$$B_d(\xi ,\eta )=B_d(\xi )\eta \overline{B}_d(\xi )$$
(B.108)
where

$$B_d(\xi )={}_{2}{}^{}F_{1}^{}(1,\frac{d}{2};1+\frac{d}{2};\xi )$$
(B.109)
and

$`\overline{B}_d(\xi )`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left[{}_{2}{}^{}F_{1}^{}(1,\frac{d}{2};1+\frac{d}{2};\xi ){\displaystyle \frac{d}{d+2}}{}_{2}{}^{}F_{1}^{}(1,1+\frac{d}{2};2+\frac{d}{2};\xi )\right].`$ (B.110)
For $`|\xi |<1`$, Eq. (B.108) can be Taylor-expanded in $`\xi `$, leading to the first equation given in Eq. (69).

For applications, it will be useful to obtain explicit analytical expressions for the functions (B.105) for fixed dimensions. For $`d=4`$, we find

$`\mathrm{}_0^4(\omega ,z,\eta )`$ $`=`$ $`{\displaystyle \frac{1}{z1}}{\displaystyle \frac{1+\omega }{(z1)^2}}\mathrm{ln}{\displaystyle \frac{z+\omega }{1+\omega }}`$ (B.112)
$`\eta \left[{\displaystyle \frac{1+2\omega +z}{4(z1)^2}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{(1+\omega )(\omega +z)}{(z1)^3}}\mathrm{ln}{\displaystyle \frac{z+\omega }{1+\omega }}\right]`$
In $`d=3`$ dimensions, we find

$`\mathrm{}_0^3(\omega ,z,\eta )`$ $`=`$ $`{\displaystyle \frac{2}{(z1)}}{\displaystyle \frac{2\sqrt{1+\omega }}{(z1)^{3/2}}}\mathrm{arctan}\sqrt{{\displaystyle \frac{z1}{1+\omega }}}`$ (B.114)
$`\eta \left[{\displaystyle \frac{1+3\omega +2z}{3(z1)^2}}\sqrt{{\displaystyle \frac{1+\omega }{z1}}}{\displaystyle \frac{z+\omega }{(z1)^2}}\mathrm{arctan}\sqrt{{\displaystyle \frac{z1}{1+\omega }}}\right]`$
and the region for $`z<1`$ is obtained through analytical continuation. Finally, for $`d=2`$ we find

$$\mathrm{}_0^2(\omega ,z,\eta )=\frac{1}{z1}\mathrm{ln}\frac{z+\omega }{1+\omega }\eta \left[\frac{1}{2}\frac{1}{z1}+\frac{z+\omega }{2(z1)^2}\mathrm{ln}\frac{z+\omega }{1+\omega }\right]$$
(B.115)
Notice that the functions (B.112), (B.114) and (B.115) behave smoothly for $`z1`$, which follows either from Eq. (69) or by an explicit check.

For $`z=1`$ these expressions reduce to the result Eq. (70).

## Appendix C High temperature limit

The high-temperature limit is reached for $`\tau 1`$, or $`T2\pi k`$. Then, only the $`n=0`$ Matsubara mode contributes to the flow. Using Eqs. (A.92), (A.95) and (B.109), we find:

$$B_d(\xi ,\tau 1,\eta )=\frac{d}{d1}\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }B_{d1}(\xi ,\tau =0,\eta )$$
(C.116)
For $`|\xi |<1`$, Eq. (C.116) can be Taylor-expanded as

$$B_d(\xi ,\tau 1,\eta )=\frac{d}{d1}\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }\underset{n=0}{\overset{\mathrm{}}{}}\frac{(d1)\xi ^n}{2n+d1}(1\frac{\eta }{d+1+2n}).$$
(C.117)
Splitting $`B_d(\xi ,\tau 1,\eta )`$ as in Eq. (A.90), we have

$$B_d(\xi ,\tau 1)=\frac{d}{d1}\frac{v_{d1}}{v_d}\frac{\tau }{2\pi }{}_{2}{}^{}F_{1}^{}(1,\frac{d1}{2};\frac{d+1}{2};\xi ).$$
(C.118)
In turn, the flow proportional to the anomalous dimension reduces to

$`\overline{B}_d(\xi ,\tau 1)`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{d}{d1}}{\displaystyle \frac{v_{d1}}{v_d}}{\displaystyle \frac{\tau }{2\pi }}\times `$ (C.120)
$`\left[{}_{2}{}^{}F_{1}^{}(1,\frac{d1}{2},\frac{d+1}{2},\xi )\frac{d1}{d+1}{}_{2}{}^{}F_{1}^{}(1,\frac{d+1}{2},\frac{d+3}{2},\xi )\right].`$