{'fraction_non_alphanumeric': 0.04980697424213657, 'fraction_numerical': 0.020062021391051198, 'mean_word_length': 3.8427827152926755, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Decay constant of B s and B ∗ s mesons from N f = 2 lattice QCD

We report on a two-flavor lattice QCD estimate of the Bs and B ∗ s leptonic decays parameterized by the decay constants fBs and fB∗ s . In addition to their relevance for phenomenology, their extraction has allowed us to investigate whether the “step scaling in mass” strategy is suitable with WlilsonClover fermions to smoothly extrapolate quantities of the heavy-strange sector up to the bottom scale. From the central value of fDs quoted by FLAG at Nf = 2 and our ratio fBs fDs , we obtain fBs = 212(10)(2)( +2 −10) MeV and fB∗ s /fBs = 1.05(2)( +0 −1).


I. INTRODUCTION
In the very active research of new effects in high-energy particle physics, flavour physics does play a key role at the so-called intensity frontier. Indeed, rare events are sensitive probes of New Physics (NP) scenarios with the exchange of extra particles in quantum loops with respect to what is known from the Standard Model (SM), However, theoretical uncertainties on hadronic quantities, for instance hadron decay constants, that encode the dynamics of QCD at large distance, severely weaken the constraints that are derived through the analysis of experimental data. Those hadronic constants cannot be reliably estimated in perturbation theory. b-quark physics is a particularly interesting place to search for NP effects and it has recently regained even stronger attention after experimental signs of several anomalies in B and B c decays. More precisely several ratios R D = Γ(B→Dτ ντ ) Γ(B→Dℓν ℓ ) ℓ=e,µ , R D * = Γ(B→D * τ ντ ) Γ(B→D * ℓν ℓ ) ℓ=e,µ , R J/ψ = Γ(Bc→J/ψτ ντ ) Γ(Bc→J/ψℓν ℓ ) ℓ=e,µ , R K = Γ(B→Kµ + µ − ) Γ(B→Ke + e − ) and R K * = Γ(B→K * µ + µ − ) Γ(B→K * e + e − ) show some discrepancy with SM expectations [1] - [12]. The three former might bring stringent constraints onbc currents, for instance mediated by the exchange of leptoquarks [13]. The further ratios R D ( * ) s = Γ(B→D ( * ) s τ ντ ) Γ(B→D ( * ) s ℓν ℓ ) ℓ=e,µ , under investigation at LHCb, will provide even more informations once, on the theory side, the hadronic matrix elements associated to B s → D ( * ) s are under comparable control by means of lattice QCD. Simulating the B s meson on the lattice is delicate as far as cut-off effects are concerned. Several strategies have been followed in the literature, including simulations of relativistic b-quarks using an action tuned so as to minimize discretization errors [14] - [17], the use of Non Relativistic QCD [18], [19], performing computations in Heavy Quark Effective Theory (HQET) [20] and the extrapolation of simulation results obtained in the region between the charm quark mass m c and a mass ∼ 3m c to the physical b-quark mass [21], [22]. As we plan to employ the latter approach to study B s decays with O(a) improved Wilson-Clover fermions, an intermediate step is to extract f Bs and f B * s , in order to validate the method. The lattice QCD community has made a significant effort to compute f Bs with N f = 2 [22], [23], N f = 2 + 1 [14], [15], [18], [24], [25] and N f = 2 + 1 + 1 [19], [26] - [28]. Recently, the SU(3) symmetry breaking f Bs /f B has been extracted at the physical point [29]. Concerning the spin-symmetry breaking ratio f B * s /f Bs only 2 lattice results are available, both at N f = 2 + 1 + 1 [30], [31]. Ratios f B * /f B and f B * s /f Bs have been investigated with other methods than lattice simulations, i.e. constituent quark models [32], [33] and QCD sum rules [34] [36].
The paper is organized as follows: in section II we recall what is the "step scaling in mass" strategy, in section III we present the simulations details and our raw data, and in section IV we describe our analysis and comment the results. Finally we conclude in section V. on a third side. In the following we focus on the two latter in the framework of N f = 2 Wilson-Clover fermions. With a given pion mass, through κ sea , the valence strange quark mass κ s and lattice spacing a, we consider the ratio where λ = mB s mD s 1 K , i = 0, · · · , K − 1, and M P (a, κ sea , κ h0 , κ s ) ≡ m Ds , up to mistuning effects.
C stat A (µ 1 , µ 2 ) and C stat V (µ 1 , µ 2 ), that will appear later in the paper, are the matching coefficients between the QCD currents J QCD A(V ) ≡Qγ 0 γ 5 q(Qγ 1 q) and their HQET counterpart J HQET where µ Q is a scale related to the heavy quark mass m Q , for instance the heavy-light pseudoscalar meson mass M P (Q, q). C stat A and C stat V are known up to 3-loop of perturbation theory [37] - [40] 1 . r P is independent of µ f because it involves the renormalization group equation.
Thanks to scaling laws in HQET, it is expected that r P has a simple expansion in the inverse heavy quark mass defined in a specific renormalization scheme, for instance the pole mass [41]. But, in the case of Wilson-Clover fermions and contrary to the case of twisted-mass fermions, there is so far no straightforward relation between the pole quark mass and the bare quark mass, through the renormalization group invariant (RGI) quark mass, if the quark is significantly heavier than the charm. Indeed, using the series of RGI masses m RGI already determined [42] and m RGI c known after the tuning of κ c , we define the improved RGI mass m RGI h(i) ∝ (1 + b m am h(i) )m h(i) , where m h(i) is the heavy vector Ward Identity quark mass m h(i) = 1 2κ h(i) − 1 2κ critical . The problem is that we can have 1 + b m am h(i) < 0 because the improvement coefficient b m is negative [43]. Negative RGI masses are of course not physical. The issue would be solved by adding the O(a 2 ) term in the definition of the RGI mass, which is unfortunately unknown. That is why we have decided to consider the inverse of pseudoscalar heavy-strange pseudoscalar meson masses M P (h(i), s) as the parameter expansion of r P and we define the steps as sequential ratios MP (h(i+1),s) MP (h(i),s) that should be constant with the regulator a and the sea quark mass κ sea . They are the analogous of the ratios of RGI quark masses

III. LATTICE CALCULATION DETAILS
We have performed our analysis from the CLS ensembles made of N f = 2 nonperturbatively O(a)improved Wilson-Clover fermions [44,45] and the plaquette gauge action [46]. In Table I we collect the main informations about the simulations. Three lattice spacings a β=5.5 = 0.04831 (38) fm, a β=5.3 = 0.06531(60) fm, a β=5.2 = 0.07513(79) fm, determined from a fit in the chiral sector [47], are considered with pion masses in the range [190 , 440] MeV. With respect to the work reported in [48], we have taken the bare strange quark masses at β = 5.2 from [49] and we have tuned the charm quark mass on those ensembles by imposing am Ds (a, κ sea , κ c , κ s ) = am physical

Ds
. The values we find for κ c are close to what is quoted 1 In this work we have taken the N 2 LO formulae for the matching between QCD and HQET at the scale µ Q .   in [49] where the tuning was realised thanks to a constraint on cut-off effect magnitude for the ratio of PCAC masses m PCAC c /m PCAC s . We have used the same procedure as in [48] to compute the statistical error at finite a and in the continuum limit, to compute stochastic all-to-all propagators and to reduce the contamination by excited states on 2-pt correlators by solving a 4 × 4 Generalized Eigenvalue Problem (GEVP) with one local and 3 Gaussian smeared interpolating fields. In our application of the step scaling in mass strategy we have chosen K = 6 steps. Similarly to [48] we have extracted the relevant matrix elements from projected correlators along the fixed ground state generalized eigenvectors v take the values (4a, 3a) at β = 5.2, (4a, 3a) at β = 5.3 and (6a, 5a) at β = 5.5. We collect in Tables X -XXV of the Appendix the whole set of raw data we need in our analysis, i.e. ratios of pseudoscalar and vector heavy-strange mesons masses for 2 subsequent heavy bare quark masses, ratios of pseudoscalar and vector heavy-strange meson decay constants for 2 subsequent heavy bare quark masses and PCAC quark masses m hs defined by: where , C A0P and C P P are axial-pseudoscalar and pseudoscalar-pseudoscalar 2-pt correlation functions of the heavy-strange meson defined by the bare quark masses (κ h , κ s ) and c A is the improvement coefficient of the axial bilinear of Wilson-Clover fermions determined in [50]. We show in Figure 1 the effective mass for two heavy-strange mesons, one at the first step scaling in mass and the other at the next to last step. For very heavy quarks the signal deteriorates quickly. It explains why we fixed shorter interval ranges to extract the hadronic properties, as indicated in Tables X -XXV.

IV. ANALYSIS AND DISCUSSION
A. Extraction of fB s We have performed extrapolations to the physical point by doing a global fit analysis. However, in a preparatory stage, we restrict our analysis to a given step in heavy mass i and study the pion mass and the cut-off dependence of r ′ P (i). Here we ignore the mistuning effects because our goal is to determine how large are the discretisation effects. We show in Figure 2 the extrapolation in a 2 and m 2 π of r P (i). We observe very big cut-off effects for the 4 th ratio (11%) and the 5 th ratio (18%). Hence we are not quite confident in using the ensembles at a = 0.075 fm and, most probably, those at a = 0.065 fm as well, at the final stage of our analysis. That is why we prefer, in the combined fit analysis, to exclude the data at the fourth and the fifth heavy quark mass at the lattice spacings 0.075 fm and 0.065 fm, Then, we have used the following fit ansatz: with r mistune (i) = mH s (i) λ i mH s (0) . We collect in Table II the fit parameters and χ 2 /d.o.f.. We show in Figure  3 the dependence of r P on 1/m Hs found with the fit formula (3). We retrieve the parametrically large cut-off effects on r P , justifying our decision to exclude some data of our 2 coarsest lattices in the analysis. (5) -0.014(2) 0.030 (9) 1.4  Then, in the continuum, we interpolate r P at the 6 points m Ds λ i+1 to get a set of 6 ratios r ′ P (i): We recall that r ′ P (i) is independent of the renormalization scale µ f . We collect in Table III values of r ′ P at the reference points 1 λ i+1 mD s . The last step is straightforward: f Bs /f Ds is obtained by a series of products: We get The effect on the statistical error of the correlation among the different terms of the product of r ′ P (i) is taken into account by the mean of computing the errors described in [48].
To address the systematic error, we have performed two other fits: -fit(A), adding to (3) a "next to leading order" chiral contribution in m 2 π ln(m 2 π ) -fit(B), adding to (3) a contribution in 1/m 2 Hs to count for a higher order in the heavy quark expansion -fit(C): fit (3) but using the matching coefficient C stat Other fits with extra terms in (am Hs ) 2 or in a 3 give non reliable results. We collect the corresponding fit parameters and χ 2 /d.o.f. in Table IV    A fourth source of systematics can be included by propagating the uncertainty on raw data if we change t min → t min + 2a to extract plateaus. In this case we get the following result: We collect the corresponding fit parameters in Table V. Adding together the different sources of systematics, we obtain f Bs f Ds = 0.87(4)( +1 −4 ).
where the first error is statistical and the second error counts for the systematic error. Concerning the D s meson decay constant, an update of the analysis reported in [48], now that we have the additional coarsest ensembles A5 and B6, meaning a third lattice spacing at our disposal, gives where the first error is statistical and the second error comes from the uncertainty on the lattice spacings. As in the previous analysis, a next to leading order contribution to the chiral fir destabilises the fit with chiral fit parameters compatible with zero.
Then we get where the first error is the statistical error, the second one counts for the systermatic error on f Ds while the third error corresponds to the systematic error on f Bs /f Ds . FLAG has recently made an update collection of lattice estimates of f Bs [51]. Our estimate of f Bs using the step scaling in mass strategy is compatible with the value obtained by the ALPHA Collaboration f Bs = 224(14) MeV [23] by a computation, performed over almost the same CLS ensembles as in this paper, of hadronic matrix elements in the framework of HQET with a non-perturbative matching of the HQET parameters with QCD. It is 2σ lower than the result reported by the ETM Collaboration [22] with N f = 2 twisted-mass fermions defined at maximal twist. The fact that we cannot constrain the static limit of the ratio r P to be equal to 1, due to mistuning effects of the heavy quark mass, explains probably that discrepancy.
As the HQET anomalous dimension of the axial and vector static-light operator are the same, applying the renormalization group equation makes R ′ f * independent of the renormalization scale µ f . We can impose the static limit constraint lim MP →∞ R m * = lim MP →∞ R ′ f * = 1 because those ratios are free of heavy quark mistuning effect. We collect in Table VI  Other terms in the fit lead to unstable and reliable results. We collect the fit parameters in Tables VII and VIII and we obtain: . As for f Bs /f Ds , we have counted for the systematic error coming from the change t min → t min + 2a in the plateaus extraction. We get the following results: Fit parameters are collected in Table IX. Eventually, we quote where the first error is statistical and the second error counts for the systematic error estimated by the fits (A') and (B') and contamination from excited states. There is no reason why the ratio

C. Comment
We collect in Figure 5 the lattice QCD estimates of f Bs at N f = 2 [22], [23], with the corresponding FLAG average [51] and those of f B * s /f Bs [30]. [31]. Of course the fact that we get f B * s /f Bs > 1 while the 2 other lattice QCD results read f B * s /f Bs < 1 is puzzling. However, a computation performed by the ETM Collaboration with N f = 2 dynamical quarks indicated the hierarchy f B * /f B > 1, f B * /f B = 1.050(16) [53]. So, a plausible explanation for the observed tension is the effect of the quenching of the strange quark in the spin-breaking contribution of the heavy quark symmetry to the ratio f B * s /f Bs . It might be of the same order of magnitude as in f D * s f Ds but with a more important qualitative impact because we examine a region of paramaters closer to the symmetric point f H * /f H = 1. In that respect, studies of this ratio with N f = 2 + 1 ensembles are welcome.

V. CONCLUSION
In that paper we have reported on a lattice estimate of f Bs and f B * s /f Bs . The main puropose of the work was testing the step scaling in mass method with Wilson-Clover fermions for which the RGI heavy quark mass can not be used yet as a physical parameter of the heavy quark expansion. Indeed, severe negative O(am) cut-effects need to be balanced by still unknown O(am) 2 improvement terms to define safely the RGI mass. Instead, we have chosen the (inverse of) the heavy-strange meson mass as the expansion parameter. We obtain a quite low result for f Bs compared to other lattice QCD estimates at N f = 2, though it is compatible with the one got using the same set of gauge ensembles as here but with a complete different approach to simulate the heavy quark. We have found the hierarchy f B * s /f Bs > 1, indicating a positive correction to the symmetric point when the strange quark is quenched. A look at the literature leads to the conclusion that this correction becomes negative when the strange quark is taken into account in the sea. The next step of our program is the investigation of the form factors associated to B s → D