Renormalisation group running effects in pp → t ¯ th in the Standard Model Effective Field Theory

We study the effects of renormalisation group running of the Wilson coefficients in Standard Model Effective Field Theory, using the process pp → t ¯ th as a showcase. We consider both strong and top Yukawa running effects, since the latter can be relevant in presence of large Wilson coefficients. We study the difference between the use of a dynamical and fixed renormalisation scale by exploring different scenarios for the higher-dimensional operators at the high energy scale of new physics, assumed to be at the TeV scale.


I. INTRODUCTION
After the discovery of a scalar particle compatible with the properties of the Standard Model (SM) Higgs boson in 2012 [1,2], all the particles predicted by the SM have been observed.However, various reasons suggest that the SM should be extended, as for instance massive neutrinos, the matter anti-matter asymmetry of the universe or the missing dark matter candidate.
Given that, so far, no striking new physics signal has been observed the new physics energy scale might lie much above the electroweak scale.This motivates to use an effective field theory approach to describe new physics in a model-independent way.In this context, the Standard Model Effective Field Theory (SMEFT) [3,4] represents a powerful tool to approach the search for heavy new physics.
Let Λ ≳ v = 246 GeV be the typical mass scale of the new degrees of freedom, lying beyond the experimental reach.The indirect effects of such heavy particles can be captured and parametrised with a tower of higher dimensional operators: where O i is an operator consisting of SM fields (and derivatives) invariant under the SM unbroken gauge group SU(3) C ⊗ SU(2) L ⊗ U(1) Y .Assuming lepton and baryon number conservation, the first terms in Eq. ( 1) arise at dimension-six level and are expected to give the dominant contributions in collider physics.A complete and non-redundant basis is defined in Ref. [4] (Warsaw basis) and in Refs.[5][6][7] (SILH basis).
Loop computations typically give an infinite result: to get rid of such divergences, the theory must be renormalised. 1An important consequence is that the parameters in the theory acquire a dependence on the energy scale, encoded in a set of differential equations known as Renormalisation Group Equations (RGEs).For what concerns the SMEFT at dimension-six and one-loop level, the RGEs have been computed in Refs.[9][10][11].
Typically, the SMEFT is used as the low-energy limit of some UV model (top-down approach).At some energy scale µ ∼ Λ, where Λ may be the typical mass scale of the new degrees of freedom, a matching between the UV model and the SMEFT is performed.The treelevel matching between the SMEFT and generic UV theories has been computed and automatised and one-loop matching is in development (see Refs. [12][13][14][15][16]).After the matching, the RGE running can be computed in the SMEFT, making it possible to address the impact of the higher dimensional operators on low-energy observables such as cross sections and differential distributions.This program offers a universal and pragmatic approach to study the phenomenological consequences of new physics including running effects that would otherwise need an ad hoc computation for each UV model.In the bottom-up approach instead one is agnostic on the Wilson coefficients of the higher-dimensional operators and uses them to parametrise deviations with respect to the SM in a model-independent way.In this way, the bounds on the 1 Formally, the Lagrangian in Eq. ( 1) is not renormalisable since it contains terms with mass dimension greater than four.This implies that the absorption of divergences requires an infinite number of counterterms.However, is it possible to renormalise the theory order by order in the expansion in 1/Λ, using a finite number of counterterms, see e.g.Ref. [8] for a more detailed discussion.
coefficients can be used to restrict the parameter space of any new physics scenario.
The goal of this paper is to test the effects of RGE running in the prediction of differential cross sections when adopting a dynamical scale choice.This means that rather than evolving the SMEFT RGEs down to some fixed electroweak scale we evolve them to a dynamical scale varying event by event.We will showcase this effect by considering the example of the process pp → t th.The running effects proportional to α s have been recently studied in Refs.[17][18][19][20].In this work, we take a step further and consider the fully coupled system of the SMEFT RGEs.The motivation is that, in particular for operators in the top sector, the running is often proportional to the top-quark Yukawa coupling.Such operators tend to be less constrained than operators whose running is proportional to α s (see Ref. [21]).For this reason, these effects can be hence of similar size or even bigger as the α s ones.In addition, we show that the first leading-log approximation, in which the RGEs are solved assuming that the anomalous dimension matrix is independent of the scale, is not sufficient to describe accurately the RGE running effects in the presence of large Wilson coefficients.This motivates to solve numerically the fully coupled system of RGEs in collider physics studies.
The paper is organized as follows: in Sec.II we describe the relevant set of operators and our computational strategy.In Sec.III we show the numerical results of our analysis in three different scenarios for the Wilson coefficients.Finally, in Sec.IV we give our conclusions.

II. DESCRIPTION OF THE COMPUTATION
The hadronic process pp → t th is generated by two different partonic channels: gg → t th and qq → t th, being q = u, d, c, s, b.In the SM this process arises at tree-level via diagrams like the ones in Fig. 1.
FIG. 1: Tree-level diagrams for pp → t th in the SM.
Considering in addition SMEFT operators gives rise to new interaction vertices, as well as to modifications of the existent ones.The complete set of Feynman rules for the SMEFT (in the broken phase), which we follow here, is given in Ref. [22].
The Lagrangian considered in this work is: The SM Lagrangian in Eq. ( 2) is given by where Q L and t R refer to the third family isospin doublet (left-handed) and isospin singlet (right-handed, uptype), respectively, φj = ε jk ϕ †k (with ε 12 = +1) and T A are the SU(3) C generators.When spontaneous symmetry breaking occurs, the physical Higgs boson acquires a vacuum expectation value v = 246 GeV leading to a mass given by m 2 H = 2λv 2 (in the SM).The part of the SMEFT Langrangian containing only the Higgs doublet and its gauge interactions is given by The interactions with the physical (neutral) Higgs boson h in the SMEFT in the unitary gauge are obtained setting ϕ = 1 √ 2 (0, v T + h[1 + c kin ]) T in Eq. (2).Higher dimensional operators produce O v 2 /Λ 2 corrections to the vacuum expectation value and give rise to a noncanonically normalised Higgs kinetic term.The latter can be put back in the canonical normalisation by means of a redefinition of the Higgs field, yielding a universal shifts of all the single Higgs couplings (see Ref. [23] for more details): ( The terms in Eq. ( 2) involving bosons and zero or two quarks are collected in While we follow Ref. [4] in the notation of the Wilson coefficients, it should be noticed that we are performing a slight change of notation: since we are focusing on the third generation, we fix C tG = C 33 uG (and in analogy for the other Wilson coefficients).
In the following, we assume the Wilson coefficients to be real, neglecting the possibility of CP violation.This choice is justified by the severe bounds on CP violating SMEFT operators from electric dipole moment, see Refs.[24,25]. 2   Finally, L 4F contains all the four-quark operators in the Warsaw basis having a non-vanishing matrix element for qq tt.Clearly these terms affect only the qq → t th channel (at tree level).
In the above expression, p, r, s, t are the generation indices, j, k are the SU(2) L indices in the fundamental representation and τ I are the Pauli matrices.It holds Obviously, we are interested in the combinations giving a non-zero matrix element between light quarks and top quarks (e.g.O 1133 uu ).A sample of the tree-level diagrams arising due to higher dimensional operators is given in Fig. 2.
In this work, we do not consider contributions where the Higgs boson couples to light quarks and gluons because such operators are not expected to show large running effects, being suppressed by light Yukawa couplings.
In Eqs. ( 4),( 6), (7) we include only operators which enter at tree-level in pp → t th.This choice does not take into account the fact that some operators (such as the chromomagnetic operator) cannot be generated at treelevel, assuming that the UV completion of the SM is weakly coupled and renormalisable (see Ref. [29,30]).In other words, Eqs. ( 4),( 6), (7) follow an expansion based only on the canonical dimension of the operators.
Four-top operators enter pp → t th at one-loop level.Schematically, their (bare) contribution goes as M ∼ 2 While some combinations of CP violation in different Higgs couplings can survive the bounds, LHC searches (see Refs. [26][27][28]) also in t th can provide complementary probes.We do not consider this possibility though in the following. FIG.
2: Tree-level diagrams for pp → t th in the SMEFT.The purple square denotes an insertion of O ϕG , the orange square denotes an insertion of O tG and the blue square denotes an insertion of L 4F .
Since the poles determine the anomalous dimension of the theory, including running effects in operators which enter at treelevel in the process means resumming the logarithmic terms.The finite terms represented by A can be phenomenologically relevant, as described in Ref. [31].However, we do not include such terms in our analysis since they do not depend on the choice of the renormalisation scale µ R and do not contribute to the difference of the result obtained with a fixed or dynamical renormalisation scale.However, for a precise phenomenological analysis, they should be included.In the same spirit, we do not include the SM NLO QCD and electroweak corrections but note that they moderately increase the cross section, see Refs.[32][33][34][35].
As it has been pointed out recently in Ref. [36], there is a non-trivial interplay between four-fermion operators and other SMEFT operators.Depending on of the continuation scheme chosen for the γ 5 in D ̸ = 4 space-time dimensions a four-fermion operator generated at treelevel, e.g.O 3333 qu (1) , automatically introduces a one-loop coefficient to e.g.O tG .We will for the moment refrain from considering this interplay but come back to this in Scenario 2 analysed in Sec.III.We note for the moment that the coefficients of the operators considered in this work can be considered independent when using the Breitenlohner-Maison-'t Hooft-Veltman (BMHV) scheme for γ 5 (see Refs. [37,38]).
We compute matrix elements at O 1 Λ 2 and we consider cross sections at O 1 Λ 4 . 3In other words, we consider only diagrams with a single insertion of a dimension-six operators but we do not include interference terms between the SM and dimension-eight operators.We generate the diagrams with qgraf-3.6.5 [39] and perform the algebra with FeynCalc [40][41][42].
While the squaring of the matrix element for the quarkinduced channels is straightforward, the sum over external polarisation states of the gluons is computationally expensive (being n a generic unit vector which must drop from the final result).We used a shortcut, based on the optical theorem (see Ref. [43] for the details).It is possible to use the simple expression we subtract incoherently the squared matrix element where the gluons are replaced by the ghosts (η).Let M gg→t th = ϵ µ1 (q 1 )ϵ µ2 (q 2 )A µ1µ2 : the unpolarized squared matrix element is given by We explicitely verified the correctness of this approach using the SM amplitude for gg → t th.
For the numerical evaluation of the differential cross section we rely on our own code written in Fortran, C++ and python.The cross section of the i-th event is given by , where m t is the top quark mass and m h the Higgs mass.The phase space PS i is generated by rambo [44] and X = gg, ūu, dd... denote the different partonic channels.The average over the initial states yelds S gg = 1 256 and S qq = 1 36 and we set the collider energy to E Collider = 13.6 TeV throughout the paper.The weight of the parton distribution function (PDF) is given by being µ F the factorisation scale set to µ F = H T /2, with H T ≡ p T,t + p T, t + p T,h .
We access the PDFs via LHAPDF-6.5.3 [45], employing the PDF set NNPDF40 lo as 01180 [46].require a positive cross section in the whole phase space.We note though that, if the cross section beecomes negative when considering O 1 Λ 2 terms only, the EFT validity might be questionable.
The total cross section is given by where N E denotes the number of events.The errors are computed as the square root of the bin sum of the σ i .
Using a dynamical scale implies that the Wilson coefficients of the SMEFT operators should be evaluated at this scale.Given the definition of the set of non-vanishing SMEFT operators at the scale µ = Λ, {C i (Λ)}, we have to solve the RGEs from Λ to µ R in order to obtain {C i (µ R )} and use them to compute the matrix element squared.We stress that this must be done for every event if one chooses a dynamical scale: high time efficiency is crucial in this context.To reduce the computational time is it possible to rely on some approximations, for instance like solving the RGEs analytically only for the dominant contribution (typically α s , see Ref. [17][18][19][20]) or employing the first leading-log approximation, with µ The first leading-log approximation is expected to provide an accurate result only if Γ ij (Λ)C j (Λ)log (µ R /Λ)/(16π 2 ) is small enough to remain inside the perturbative regime.We will come back on this topic in Scenario 3 in Sec.III.
We solve the RGEs numerically using RGESolver version v1.0.2, presented in Ref. [47]. 4 The numerical approach allows us to consider not only O (α s ) effects but also O (Y t ) effects and address their importance.
Differential cross sections are obtained by respective binning of the events.The initial conditions for the SM parameters are obtained evolving the pure SM RGEs up to Λ using the default input provided by RGESolver (meaning Table 4 of Ref. [47]).In the following, we distinguish between the three set-ups.In all these cases, we evolve from µ = Λ to µ = µ R .
• SM: same as before but setting at µ = Λ all the SMEFT coefficients to 0. This choice represents the benchmark with respect to which we compare the SMEFT effects, meaning that we use for the SM input parameters and their running exactly the same input as for the SMEFT case.
• SMEFT (Fixed Scale): a fixed scale µ R = m t is used for all the parameters of the theory.To isolate the contribution coming from pure SMEFT running, the strong coupling α s is still evaluated at the dynamical scale µ F = H T /2.Further discussion regarding this aspect can be found in the next section for the operator O (1) Qt .
The masses are renormalised on-shell: we use m t = 173 GeV, m h = 125 GeV.The effective coupling of Higgs to top quarks is given by: The vacuum expectation value is not directly measured.
It is derived from G F , which is measured in muon decay.When SMEFT operators are considered, this input value is modified by the operators O 11,22  Hl(3) and O 1221 ll (see Ref. [50]).We neglect this effect since such operators do not play any role in our analysis, see Refs.[51][52][53][54] for detailed analyses.

III. RESULTS
We show the dependence of the differential cross section on the running of the Wilson coefficients in three different scenarios.All the plots displayed in this section have been obtained simulating 10 5 events.The whole computation (event generation, running and matrix element evaluation) takes O (20 min) on a normal laptop.
We show in many plots the percentual difference between the fixed and dynamical scale in each bin, namely a. Scenario 1 : We turn on one Wilson coefficient at the time, using Λ = 2 TeV.The idea of this scenario is to demonstrate the importance of Y t running with respect to pure α s running.To do so, we compare two cases: we turn on individually the four-top operators O (1,8) Qt

≡ O 3333
qu (1,8) , choosing as initial conditions (i) C Qt = 20 /TeV 2 .This choice, consistent with the individual bounds at O 1/Λ 2 presented in Ref. [21], is motivated by the fact that these operators enter the β-functions of C tϕ with pure Yukawainduced terms in the combination C (1) Qt , being C F = 4/3 for SU(3) C (see Refs. [9,10]).This concretely means that both choices (i) and (ii) generate the same contribution to the running of C tϕ at leading-logarithm, hence the same Yukawa-induced running in both cases.However, O Qt contributes to the α s running of other operators entering at tree-level in pp → t th via, for example, penguin diagrams (e.g. it generates C 1133 uu (µ R ) ̸ = 0).This is not the case for O Qt .For this reason, a comparison between the two operators provides a good indicator of the importance of purely Yukawa-induced RGE running effects with respect to the ones proportional to the strong coupling for this process.
We report the differential distribution with respect to the Higgs transverse momentum p T,h for C (1) Qt in Fig. 3 and for C (8) Qt in Fig. 5.We observe in both cases differences in the distribution obtained with the dynamical scale with respect to the one obtained with the fixed scale below 10%.The two choices of µ R agree rather well in the [50,100] GeV bin.In both cases, we see a small enhancement at low p T,h , which can be qualitatively understood looking at the energy-scale dependence of g ht t, displayed in Fig. 4 (Fig. 6) for C Qt ).The small difference between Fig. 4 and Fig. 6 can be attributed to the non leading-log contributions in the running of C tϕ .

Since O (1)
Qt does not mix with operators contributing at tree-level to pp → t th via QCD interactions, the fact that we find a (small) deviation between fixed and dynamical scale means that purely top Yukawa-induced running has an effect on the distribution.For this reason, it should be included in presence of large Wilson coefficients.Indeed, we find in this case that both the operators with Wilson coefficients O (1)

Qt and O (8)
Qt show deviations between fixed and dynamical scale of the same order, which lets us conclude that the effects stemming from top Yukawa-induced running can be as important as the ones propotional to α s .
It is also interesting to look at the correlation between p T,h (which is a physical observable) and µ R (which is an arbitrary energy scale).We quantify this correlation looking at d 2 σ dµR dp T ,h , reporting the result for C (1) Qt in Fig. 7.We expect roughly the same result for the bin p T,h ∈ [50, 100] × GeV, since most of the events falling in such interval produce µ R ∼ m t .This is precisely what can be inferred from Fig. 7.We also note that the largest part of Lastly, we comment about the choice of using a different renormalisation scale for the strong coupling (µ R = µ F = H T /2) and for the SMEFT coefficients (µ R = m t ) in the fixed-scale case.This choice allows to highlight the importance of the running of the SMEFT coefficients, separating it from pure SM running of α s .For the sake of clarity, we also report in Fig. 8 the results for C (1) Qt in which the renormalisation scale in the fixed-scale scenario is set to µ R = m t also for the strong coupling.With this choice, we can see a larger difference between the fixed scale and the dynamical scale distribution.However, part of this difference is due to pure SM strong running, as can be inferred from the lower panel of the same figure.SM a new scalar Φ ∼ (8, 2) 1 2 with a mass M Φ ≫ v: The tree-level matching, performed at µ = Λ, gives rise to (see Ref. [55]): At one-loop level the chromomagnetic operator can be generated. 5The value of its coefficient depends on the continuation scheme chosen for γ 5 : The index NDR stands for naïve dimensional regularisation.We refer to its exact definition to Ref. [36] which also discusses exhaustively the scheme-independence of the final result.Relevant for the following discussion is the fact that in the NDR scheme the four-top operators contribute at two-loop order to the running of C ϕG .Its β-function is thus given by (19) We note that the two terms are indeed of the same order, since O tG is generated at one-loop level in the matching.If one would adopt a loop counting of the operators this would be directly manifest within the RGEs.
Within this scenario we want to compare the importance of the latter term in comparison to the first one which has been included in previous analyses on running effects, see Refs.[19,56].We choose Y 2 Φ /M 2 Φ = 100/TeV 2 , which provides SMEFT coefficients which are within the marginalised bounds at O 1/Λ 2 presented in Ref. [21].We use the first leading-log running (see Eq. ( 13)) for all the operators and for the two-loop contributions as the full two-loop RGEs in SMEFT are yet unknown, with some partial results presented in Refs.[57,58].We report in Fig. 9a the results.The two set-ups do not exhibit a large difference (below 10%).The reason is mainly due to the fact that both terms in Eq. ( 19) are effectively two-loop order terms (while the SM process arises at tree-level), so their impact is smaller with respect to the one-loop running effects shown in Figs. 3 and 5.This can also be seen in Fig. 9b.
If the two-loop contribution is not included, C ϕG increases as the renormalisation scale approaches 0. When, instead the two-loop contribution of the four-top operators is included, C ϕG = 0 at all scales in the leadinglogarithmic approximation.We stress that this contribution arises only when NDR is used, being absent in  BMHV and that the BMHV and NDR results correspond to each other only if the full Eq.( 19) is considered and not just the piece proportional to C tG .
c. Scenario 3 : In this scenario we want to maximalise the effects of the RG running.Hence we use the maximal values allowed for the Wilson coefficients by global fits, i.e. adopting the limits from Ref. [21], at the high energy scale Λ.We note though that we do not take into account the correlation in the fit, hence our scenario is not compatible with it.However, the scope of this study is not to provide a phenomenologically viable scenario, but it aims to address the importance of running effects.
We first rely on the most conservative bounds, namely the marginalised ones obtained at O 1/Λ 4 .These bounds provide smaller Wilson coefficients for the operators O (1,8) Qt with respect to the previous scenarios.We consider all the four-quark operators (both four-top and two-top-two-light), as well as O ϕ□ , O ϕD , O tϕ , O tG , O ϕG (translating the basis used in Ref. [21] to ours).This choice of coefficients is referred to as L4.We set all the operators to one of the extremes of each interval.The result is shown in Fig. 10.We observe a difference between the dynamical and the fixed scale of ∼ 25% in the kinematic tails.
We now relax the assumptions and consider the O 1/Λ 2 bounds.We set the four-top operators to 100/Λ 2 and the two-top-two-light to 1/Λ 2 , with Λ = 2 TeV. 6We set the other operators to one of the extremes of the bound interval.This choice of coefficients is referred to as L2.The result is shown in Fig. 11.In this case, the effect is larger, up to 70% in the kinematic tails.
Finally, we present a comparison between two resummation strategies: numeric solution of the RGEs and first leading-logarithm approximation in Eq. ( 13).The former is more precise, but computationally more expensive (within RGESolver, this method is ∼ 50 times slower than the approximated solution).In order to isolate the contribution from pure SMEFT running, also in this case we use the one-loop SM running for the strong coupling constant g s .The percentual difference in this case is defined as The results are reported in Fig. 12 (Fig. 13) for the bounds in scenario L4 (scenario L2).In the first case, we do not observe a large deviation between the two methods, oscillating in [−5%, 5%].In the latter case the Wilson coefficients are less constrained and we can see a larger difference, up to 15%.We thus conclude that, in presence of large Wilson coefficients, the first leading-logarithm solution is not a valid method and to obtain precise results and a numeric solution should be employed.

IV. CONCLUSIONS
We have computed the Higgs transverse momentum distribution for pp → t th in the SMEFT.We have compared different choices for the renormalisation and factorisation scale: a dynamical scale µ R = µ F = H T /2 (varying event by event) and a fixed scale µ R = m t .Understanding if the latter choice is a valid approximation is crucial, due to the large computational cost of the former.We computed the renormalisation group effects at one-loop level numerically using RGESolver.This made it possible to include not only the running of the strong coupling constant as for instance done in Refs.[17][18][19][20] but also the contributions proportional to the top quark Yukawa coupling.
We studied three different scenarios for the Wilson coefficients.In Scenario 1 we analysed individually the effect of the four-top operators C (1,8) Qt .We noted a small difference between the two choices of renormalisation scale, with substantial agreement around p T,h ∼ 100 GeV.Moreover, we argued that for these operators the running effects proportional to the top Yukawa coupling can be important.
In Scenario 2 we studied the interplay of the same operators with the chromomagnetic operator, in the convenient framework of an UV model involving a heavy scalar Φ ∼ (8, 2) 1/2 .In particular, we studied the impact of the two-loop contributions from four-top operators to the Higgs-gluon coupling.Remarkably, these contributions depend on the γ 5 scheme, as detailed in Ref. [36].In this work, we have shown that their effect is rather small for pp → t th, being a two-loop effect, while the SM matrix element arises at tree-level.However, we stress that this effect can be more relevant for other processes.
Finally, in Scenario 3 we studied the case where the Wilson coefficients assume the maximal values allowed by the bounds presented in Ref. [21].While this scenario is somewhat not very realistic, at the same time it gives an estimate on how large the running effects can be.We found an effect up to 25% if we use the L4 scenario where we have set the Wilson coefficient to one of the bounds quoted at O 1/Λ 4 in Ref. [21].If we use the scenario L2 where we set the four-fermion operator coefficients to 100/TeV 2 we observed even bigger effects, up to 70 % in the high transverse momentum bins.The importance of renormalisation group effects depends crucially on the size of the Wilson coefficients.In many scenarios where the new physics couples dominantly to top quarks, leading to large coefficients for the four-top operators, such effects should be taken into account.
Lastly, we compared the numeric solution to the first leading-log approximation.When employing the L4 setting of the Wilson coefficients the difference between the two methods is small, below 5%.When instead the large Wilson coefficients of L2 are used, the two methods show a sizeable deviation from each other, up to 15%.We thus conclude that, when Wilson coefficients are allowed to be large, the approximate solution is not a reliable strategy, calling for a numeric solution of the renormalisation group equations.

ACKNOWLEDGMENTS
The Feynman diagrams shown in this work were drawn with TikZ-Feynman, see Ref. [59].SDN would like to thank Gabriele Levati, Laura Reina and Luca Silvestrini for useful discussions.This project has received funding from the European Union's Horizon Europe research and innovation programme under the Marie Sk lodowska-Curie Staff Exchange grant agreement No 101086085 -  SDN also thanks the Lawrence Berkeley National Laboratory and the Berkeley Center for Theoretical Physics for hospitality.
Higgs transverse momentum distribution (upper panel) and percentual difference for each bin between the one-and two-loop running computed as ∆ ≡ dσ dp T ,h 1L − dσ dp T ,h 2L / dσ dp T ,h 2L (lower panel).

FIG. 9 :
FIG.9: Comparison between the one-and two-loop running of C ϕG .

FIG. 10 :
FIG.10: Higgs transverse momentum distribution (upper panel) and percentual difference for each bin between the dynamical scale and the fixed scale (lower panel).The coefficients are set at µ = Λ following L4, see text for the details.

FIG. 11 :FIG. 12 :FIG. 13 :
FIG.11: Higgs transverse momentum distribution (upper panel) and percentual difference for each bin between the dynamical scale and the fixed scale (lower panel).The coefficients are set at µ = Λ following L2, see text for the details.