Search

We present the lattice simulation of the renormalization group flow in the $3$-dimensional $O(N)$ linear sigma model. This model possesses a nontrivial infrared fixed point, called Wilson--Fisher fixed point. Arguing that the parameter space of running coupling constants can be spanned by expectation values of operators evolved by the gradient flow, we exemplify a scaling behavior analysis based on the gradient flow in the large $N$ approximation at criticality. Then, we work out the numerical simulation of the theory with finite $N$. Depicting the renormalization group flow along the gradient flow, we confirm the existence of the Wilson--Fisher fixed point non-perturbatively., Comment: 11 pages, 6 figures, talk presented at the 41st International Symposium on Lattice Field Theory (Lattice2024), July 28th - August 3rd, 2024, The University of Liverpool

Employing the modified Villain lattice formulation of the axion quantum electrodynamics, we present an alternative and much simpler derivation of the conclusion of~Ref.~\cite{Honda:2024sdz} that the sweep of the axial $U(1)$ non-invertible symmetry operator over the (non-genuine) gauge invariant 't~Hooft line operator with an integer magnetic charge does not leave any effect. The point is that such a 't~Hooft line can be represented by a boundary of a (non-topological) defect that is invariant under the axial transformation on the axion field., Comment: 15 pages, 1 figure. Added a note and an appendix. The final version to appear in PTEP

We study how the symmetry operator of the axial $U(1)$ non-invertible symmetry acts on the 't~Hooft line operator in the $U(1)$ gauge theory by employing the modified Villain-type lattice formulation. We model the axial anomaly by a compact scalar boson, the ``QED axion''. For the gauge invariance, the simple 't~Hooft line operator, which is defined by a line integral of the dual $U(1)$ gauge potential, must be ``dressed'' by the scalar and $U(1)$ gauge fields. A careful consideration on the basis of the anomalous Ward--Takahashi identity containing the 't~Hooft operator with the dressing factor and a precise definition of the symmetry operator on the lattice shows that the symmetry operator leaves no effect when it sweeps out a 't~Hooft loop operator. This result appears inequivalent with the phenomenon concluded in the continuum theory. In an appendix, we demonstrate that the half-space gauging of the magnetic $\mathbb{Z}_N$ 1-form symmetry, when formulated in an appropriate lattice framework, leads to the same conclusion as above. A similar result is obtained for the axion string operator., Comment: 21 pages, 1 figure. The final version to appear in PTEP

Recently, lattice formulations of Abelian chiral gauge theory in two dimensions have been devised on the basis of the Abelian bosonization. A salient feature of these 2D lattice formulations is that the gauge invariance is \emph{exactly\/} preserved for anomaly-free theories and thus is completely free from the question of the gauge mode decoupling. In the present paper, we propose a yet another lattice formulation sharing this desired property. A particularly unique point in our formulation is that the vertex operator of the dual scalar field, which carries the vector charge of the fermion and the ``magnetic charge'' in the bosonization, is represented by a ``hole'' excised from the lattice; this is the excision method formulated recently by Abe et al. in a somewhat different context., Comment: 16 pages, 3 figures, the final version to appear in PTEP

In conventional ionization chambers (ICs) using P-10 (Ar+CH4) gas, as the atomic number (Z) of the ion beams increases in the energy region of 200-300 MeV/u, the Z resolution deteriorates rapidly when Z>70. This degradation is attributed to substantial energy loss straggling caused by charge state fluctuation when the beams traverse a gas medium. The energy loss straggling increases when the beams cannot attain charge state equilibrium in the IC gas. In this study, a xenon-based gas (Xe+CH4), exhibiting a sufficiently large charge state changing cross section, was used in the IC to reach charge state equilibrium. The responses of ICs with P-10 and the xenon-based gases were examined using 238U beams and cocktail radioactive isotope (RI) beams with Z=40-90 at the RI Beam Factory (RIBF). For 238U beams at 165-344 MeV/u, the P-10 gas IC yielded an energy resolution of 1.9-3.0% in full width at half maximum (FWHM), which proved inadequate for Z identification in the uranium region. In contrast, the xenon-based gas IC demonstrated a satisfactory energy resolution of 1.4-1.6%. When using cocktail RI beams, a Z resolution of 1.28 and 0.74 was achieved by the P-10 and the xenon-based gas ICs, respectively, for beams with Z=84-88 at 200 MeV/u. The contrast in Z resolutions between the P-10 and the xenon-based gas ICs was effectively elucidated by the energy loss straggling model, incorporating collisional straggling and straggling due to charge state changes in the IC gases. The xenon-based gas IC, with more than 3sigma Z separation across a broad Z range (Z=40-90), emerged as a practical solution for Z identification of heavy ion beams., Comment: 24 pages, 13 figures, 1 table

In $U(1)$ lattice gauge theory with compact $U(1)$ variables, we construct the symmetry operator, i.e.\ the topological defect, for the axial $U(1)$ noninvertible symmetry. This requires a lattice formulation of chiral gauge theory with an anomalous matter content and we employ the lattice formulation on the basis of the Ginsparg--Wilson relation. The invariance of the symmetry operator under the gauge transformation of the gauge field on the defect is realized, imitating the prescription by Karasik in continuum theory, by integrating the lattice Chern--Simons term on the defect over \emph{smooth\/} lattice gauge transformations. The projection operator for allowed magnetic fluxes on the defect then emerges with lattice regularization. The resulting symmetry operator is manifestly invariant under lattice gauge transformations. In an appendix, we give another way of constructing the symmetry operator on the basis of a 3D $\mathbb{Z}_N$ topological quantum field theory, the level-$N$ BF theory on the lattice., Comment: 47 pages, 5 figures, the final version to appear in PTEP

In this talk, we give the lattice regularized formulation of the mixed 't Hooft anomaly between the $\mathbb{Z}_N$ $1$-form symmetry and the $\theta$ periodicity for $4$d pure Yang-Mills theory, which was originally discussed by Gaiotto $\textit{et al.}$ in the continuum description. For this purpose, we define the topological charge of the lattice $SU(N)$ gauge theory coupled with the background $\mathbb{Z}_N$ $2$-form gauge fields $B_p$ by generalizing L\"uscher's construction of the $SU(N)$ topological charge. We show that this lattice topological charge enjoys the fractional $1/N$ shift completely characterized by the background gauge field $B_p$, and this rigorously proves the mixed 't Hooft anomaly with the finite lattice spacings. As a consequence, the Yang-Mills vacua at $\theta$ and $\theta+2\pi$ are distinct as the symmetry-protected topological states when the confinement is assumed., Comment: 8 pages, 2 figures, talk presented at the 40th International Symposium on Lattice Field Theory (Lattice2023), July 31st - August 4th, 2023, Fermi National Accelerator Laboratory

In lattice compact gauge theories, we must impose the admissibility condition to have well-defined topological sectors. The admissibility condition, however, usually forbids the presence of magnetic operators, and it is not so trivial if one can study the monopole physics depending on the topological term, such as the Witten effect, on the lattice. In this paper, we address this question in the case of 2D compact scalars as it would be one of the simplest examples having analogues of the monopole and the topological term. To define the magnetic operator, we propose the ``excision method,'' which consists of excising lattice links (or bonds) in an appropriate region containing the monopole and defining the dual lattice in a particular way. The size of the excised region is $O(1)$ in lattice units so that the monopole becomes point-like in the continuum limit. We give the lattice derivation of the 't~Hooft anomalies between the electric and magnetic symmetries and also derive the higher-group-like structure related to the Witten effect., Comment: 21 pages, 4 figures

The gradient flow exact renormalization group (GFERG) is an idea that incorporates gauge invariant gradient flows into the formalism of the exact renormalization group (ERG). GFERG introduces a Wilson action with a cutoff while keeping vector gauge invariance manifestly. The details of the formalism are still to be worked out. In this paper, we apply GFERG to construct the Wilson action of massless Dirac fermions under the background chiral gauge fields. By formulating the chiral anomaly as a ``composite operator,'' we make the scale invariance of the anomaly manifest. We argue that the same result extends to QCD., Comment: 25 pages

We extend the definition of L\"uscher's lattice topological charge to the case of $4$d $SU(N)$ gauge fields coupled with $\mathbb{Z}_N$ $2$-form gauge fields. This result is achieved while maintaining the locality, the $SU(N)$ gauge invariance, and $\mathbb{Z}_N$ $1$-form gauge invariance, and we find that the manifest $1$-form gauge invariance plays the central role in our construction. This result gives the lattice regularized derivation of the mixed 't Hooft anomaly in pure $SU(N)$ Yang-Mills theory between its $\mathbb{Z}_N$ $1$-form symmetry and the $\theta$ periodicity., Comment: 20 pages, 3 figures; (v2) minor improvements, refs updated

Unbound states in 17C were investigated via one-neutron removal from a 18C beam at an energy of 245 MeV/nucleon on a carbon target. The energy spectrum of 17C, above the single-neutron decay threshold, was reconstructed using invariant mass spectroscopy from the measured momenta of the 16C fragment and neutron, and was found to exhibit resonances at Er=0.52(2), 0.77(2), 1.36(1), 1.91(1), 2.22(3) and 3.20(1) MeV. The resonance at Er=0.77(2) MeV [Ex=1.51(3) MeV] was provisionally assigned as the second 5/2+ state. The two resonances at Er=1.91(1) and 3.20(1) MeV [Ex=2.65(2) and 3.94(2) MeV] were identified, through comparison of the energies, cross sections and momentum distributions with shell-model and eikonal reaction calculations, as p-shell hole states with spin-parities 1/2- and 3/2-, respectively. A detailed comparison was made with the results obtained using a range of shell-model interactions. The YSOX shell-model Hamiltonian, the cross-shell part of which is based on the monopole-based universal interaction, was found to provide a very good description of the present results and those for the neighbouring odd-A carbon isotopes - in particular for the negative parity cross-shell states., Comment: Physics Letters B 836 (2023) 137629, published

We study phase structure and critical point of finite-temperature QCD in the heavy-quark region applying the hopping parameter expansion (HPE). We first study finite-size scaling on the critical point on $N_t=4$ lattices with large spatial volumes taking the leading order (LO) and the next-to-leading order (NLO) effects of the HPE, and find that the critical scaling of the Z(2) universality class expected around the critical point of two-flavor QCD is realized when the aspect ratio of the lattice is larger than about 9. This enables us to determine the critical point in the thermodynamic limit with high precisions. By a study of the convergence of the HPE, we confirm that the result of the critical point with the LO (NLO) approximation of the HPE is fairly accurate for $N_t=4$ (6), while we need to incorporate higher order effects for larger $N_t$. To extend the study to large $N_t$ lattices, we then develop a method to take the effects of higher-order terms of the HPE up to a sufficiently high order. We report on the status of our study on $N_t = 6$ lattice adopting the new method., Comment: 9 pages, 5 figures, Proceedings of the 39th International Symposium on Lattice Field Theory (LATTICE2022), 8th-13th August 2022, Bonn, Germany

We construct a non-trivial $U(1)/\mathbb{Z}_q$ principal bundle on~$T^4$ from the compact $U(1)$ lattice gauge field by generalizing L\"uscher's constriction so that the cocycle condition contains $\mathbb{Z}_q$ elements (the 't~Hooft flux). The construction requires an admissibility condition on lattice gauge field configurations. From the transition function so constructed, we have the fractional topological charge that is $\mathbb{Z}_q$ one-form gauge invariant and odd under the lattice time reversal transformation. Assuming a rescaling of the vacuum angle $\theta\to q\theta$ suggested from the Witten effect, our construction provides a lattice implementation of the mixed 't~Hooft anomaly between the $\mathbb{Z}_q$ one-form symmetry and the time reversal symmetry in the $U(1)$ gauge theory with matter fields of charge~$q\in2\mathbb{Z}$ when $\theta=\pi$, which was studied by Honda and Tanizaki [J. High Energy Phys. \textbf{12}, 154 (2020)] in the continuum framework., Comment: 18 pages. The final version to appear in PTEP

Modern theories of physics tell that the vacuum is not an empty space. Hidden in the vacuum is a structure of anti-quarks $\bar{q}$ and quarks $q$. The $\bar{q}$ and $q$ pair has the same quantum number as the vacuum and condensates in it since the strong interaction of the quantum chromodynamics (QCD) is too strong to leave it empty. The $\bar{q}q$ condensation breaks the chiral symmetry of the vacuum. The expectation value $<\bar{q}q>$ is an order parameter. For higher temperature or higher matter-density, $|<\bar{q}q>|$ decreases reflecting the restoration of the symmetry. In contrast to these clear-cut arguments, experimental evidence is so far limited. First of all, the $\bar{q}q$ is nothing but the vacuum itself. It is neither visible nor perceptible. In this article, we unravel this invisible existence by high precision measurement of pionic atoms, $\pi^-$-meson-nucleus bound systems. Using the $\pi^-$ as a probe, we demonstrate that $|<\bar{q}q>|$ is reduced in the nucleus at 58% of the normal nuclear density by a factor of 77 $\pm$ 2% compared with that in the vacuum. This reduction indicates that the chiral symmetry is partially restored due to the extremely high density of the nucleus. The present experimental result clearly exhibits the existence of the hidden structure, the chiral condensate, in the vacuum.

We define a one-particle irreducible (1PI) Wilson action in the gradient flow exact renormalization group (GFERG) formalism as the Legendre transform of a Wilson action. We consider quantum electrodynamics in particular, and show that the GFERG flow equation preserves the invariance of the 1PI Wilson action (excluding the gauge-fixing term) under the \emph{conventional\/} $U(1)$ gauge transformation. This is in contrast to the invariance of the original Wilson action under a modified $U(1)$ gauge transformation. The global chiral transformation also takes the \emph{conventional\/} form for the 1PI Wilson action. Despite the complexity of the GFERG flow equation, the conventional form of the gauge and global chiral transformations may allow us to introduce a non-perturbative Ansatz for gauge and chiral invariant 1PI Wilson actions., Comment: 21 pages, the final version to appear in PTEP

We formulate quantum electrodynamics on the basis of gauge (or BRST) covariant diffusion equations of fields. This is a particular example of the gradient flow exact renormalization group (GFERG). The resulting Wilson action fulfills a simple gauge Ward--Takahashi identity. We solve the GFERG equation around the Gaussian fixed point to the second order in gauge coupling and obtain the 1-loop beta function and anomalous dimensions. The anomalous dimension of the electron field coincides with that of the fermion field diffused by a gauge covariant flow equation of L\"uscher., Comment: 37 pages

We study the latent heat and the pressure gap between the hot and cold phases at the first-order transition temperature $T=T_c$ of SU(3) Yang-Mills theory, using the small flow-time expansion (SF$t$X) method based on the gradient flow. We first examine alternative procedures in the SFtX method -- the order of the continuum and vanishing flow-time extrapolations. We confirm that the final results adopting the two orders, as well as other alternatives in which the perturbative order of the matching coefficients and the renormalization scale of the flow scheme are varied, are all consistent with each other. We also confirm $\Delta p$ is consistent with zero, as expected from the dynamical balance of two phases at $T_c$. For the latent heat in the continuum limit, we find $\Delta \epsilon /T^4 = 1.117(40)$ for the spatial volume $L^3$ corresponding to the aspect ratio $N_s/N_t=T_cL=8$ and $1.349(38)$ for $N_s/N_t=6$. From hysteresis curves, we show that the entropy density in the hot phase is sensitive to the spatial volume, while that in the confined phase is insensitive., Comment: 9 pages, 6 figures, poster presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), July 26-30 2021, Zoom/Gather@MIT, USA

The gradient flow exact renormalization group (GFERG) is an exact renormalization group motivated by the Yang--Mills gradient flow and its salient feature is a manifest gauge invariance. We generalize this GFERG, originally formulated for the pure Yang--Mills theory, to vector-like gauge theories containing fermion fields, keeping the manifest gauge invariance. For the chiral symmetry we have two options: one possible formulation preserves the conventional form of the chiral symmetry and the other simpler formulation realizes the chiral symmetry in a modified form \`a la Ginsparg--Wilson. We work out a gauge-invariant local Wilson action in quantum electrodynamics to the lowest nontrivial order of perturbation theory. This Wilson action reproduces the correct axial anomaly in~$D=2$., Comment: 23 pages, the final version to appear in PTEP

We theoretically clarify the functional form to be used in $t \to 0$ extrapolation in the small flow time expansion (SF$t$X) method for the energy-momentum tensor (EMT), which facilitates lattice simulation of the EMT based on the gradient flow. We argue that in the $t \to 0$ extrapolation analysis, lattice data should be fitted by a power function in $g(\mu(t))$, the flow time dependent running coupling, where the power is determined by the perturbation order we consider. From actual lattice data, we confirm the validity of the extrapolation function. Using the new extrapolation function, we present updated lattice results for thermodynamics quantities in quenched QCD; our results are consistent with the previous study [arXiv:1812.06444] but we obtain smaller errors due to reduction of systematic errors., Comment: 46 pages, 8 figures, version to appear in PTEP

The gradient flow bears a close resemblance to the coarse graining, the guiding principle of the renormalization group (RG). In the case of scalar field theory, a precise connection has been made between the gradient flow and the RG flow of the Wilson action in the exact renormalization group (ERG) formalism. By imitating the structure of this connection, we propose an ERG differential equation that preserves manifest gauge invariance in Yang--Mills theory. Our construction in continuum theory can be extended to lattice gauge theory., Comment: 22 pages, the final version to appear in PTEP

We study latent heat and the pressure gap between the hot and cold phases at the first-order deconfining phase transition temperature of the SU(3) Yang-Mills theory. Performing simulations on lattices with various spatial volumes and lattice spacings, we calculate the gaps of the energy density and pressure using the small flow-time expansion (SFtX) method. We find that the latent heat $\Delta \epsilon$ in the continuum limit is $\Delta \epsilon /T^4 = 1.117 \pm 0.040$ for the aspect ratio $N_s/N_t=8$ and $1.349 \pm 0.038$ for $N_s/N_t=6$ at the transition temperature $T=T_c$. We also confirm that the pressure gap is consistent with zero, as expected from the dynamical balance of two phases at $T_c$. From hysteresis curves of the energy density near $T_c$, we show that the energy density in the (metastable) deconfined phase is sensitive to the spatial volume, while that in the confined phase is insensitive. Furthermore, we examine the effect of alternative procedures in the SFtX method - the order of the continuum and the vanishing flow-time extrapolations, and also the renormalization scale and higher-order corrections in the matching coefficients. We confirm that the final results are all very consistent with each other for these alternatives., Comment: 29 pages, 22 figures

To study the CP-violation using the $K_0-\bar{K}_0$ oscillation, we need the kaon bag parameter which represents QCD corrections in the leading Feynman diagrams. The lattice QCD provides us with the only way to evaluate the kaon bag parameter directly from the first principles of QCD. However, a calculation of relevant four quark operators with theoretically sound Wilson-type lattice quarks had to carry a numerically big burden of extra renormalizations and resolution of extra mixings due to the explicit chiral violation. Recently, the Small Flow-time eXpansion (SFtX) method was proposed as a general method based on the gradient flow to correctly calculate any renormalized observables on the lattice, irrespective of the explicit violations of related symmetries on the lattice. To apply the SFtX method, we need matching coefficients, which relate finite operators at small flow-times in the gradient flow scheme to renormalized observables in conventional renormalization schemes. In this paper, we calculate the matching coefficients for four quark operators and quark bi-linear operators, relevant to the kaon bag parameter., Comment: 23 pages, 3 figures. Published version: a discussion in the introduction supplemented, typos corrected, figures replaced by their vector (pdf) versions, results and conclusions not changed

We present additional observations to previous studies on the infrared (IR) renormalon in $SU(N)$ QCD(adj.), the $SU(N)$ gauge theory with $n_W$-flavor adjoint Weyl fermions on~$\mathbb{R}^3\times S^1$ with the $\mathbb{Z}_N$ twisted boundary condition. First, we show that, for arbitrary finite~$N$, a logarithmic factor in the vacuum polarization of the "photon" (the gauge boson associated with the Cartan generators of~$SU(N)$) disappears under the $S^1$~compactification. Since the IR renormalon is attributed to the presence of this logarithmic factor, it is concluded that there is no IR renormalon in this system with finite~$N$. This result generalizes the observation made by Anber and~Sulejmanpasic [J. High Energy Phys.\ \textbf{1501}, 139 (2015)] for $N=2$ and~$3$ to arbitrary finite~$N$. Next, we point out that, although renormalon ambiguities do not appear through the Borel procedure in this system, an ambiguity appears in an alternative resummation procedure in which a resummed quantity is given by a momentum integration where the inverse of the vacuum polarization is included as the integrand. Such an ambiguity is caused by a simple zero at non-zero momentum of the vacuum polarization. Under the decompactification~$R\to\infty$, where $R$ is the radius of the $S^1$, this ambiguity in the momentum integration smoothly reduces to the IR renormalon ambiguity in~$\mathbb{R}^4$. We term this ambiguity in the momentum integration "renormalon precursor". The emergence of the IR renormalon ambiguity in~$\mathbb{R}^4$ under the decompactification can be naturally understood with this notion., Comment: 30 pages, 3 figures, the final version to appear in PTEP

We study thermodynamic properties of Nf=2+1 QCD on the lattice adopting O(a)-improved Wilson quark action and Iwasaki gauge action. To cope with the problems due to explicit violation of the Poincare and chiral symmetries, we apply the Small Flow-time eXpansion (SFtX) method based on the gradient flow, which is a general method to correctly calculate any renormalized observables on the lattice. In this method, the matching coefficients in front of operators in the small flow-time expansion are calculated by perturbation theory. In a previous study using one-loop matching coefficients, we found that the SFtX method works well for the equation of state, chiral condensates and susceptibilities. In this paper, we study the effect of two-loop matching coefficients by Harlander et al. We also test the influence of the renormalization scale in the SFtX method. We find that, by adopting the mu_0 renormalization scale of Harlander et al. instead of the conventional mu_d=1/sqrt{8t} scale, the linear behavior at large t is improved so that we can perform the t -> 0 extrapolation of the SFtX method more confidently. In the calculation of the two-loop matching coefficients by Harlander et al., the equation of motion for quark fields was used. For the entropy density in which the equation of motion has no effects, we find that the results using the two-loop coefficients agree well with those using one-loop coefficients. On the other hand, for the trace anomaly which is affected by the equation of motion, we find discrepancies between the one- and two-loop results at high temperatures. By comparing the results of one-loop coefficients with and without using the equation of motion, the main origin of the discrepancies is suggested to be attributed to O((aT)^2)=O(1/N_t^2) discretization errors in the equation of motion at N_t =< 10., Comment: 45 pages, 13 figures. Erratum included to the published version. Errors in Eqs. (50) and (51) as well as in Figs. 2-7, 10, and 12 corrected. See PhysRevD.102.059903 for details. Analyses have been done adopting the correct formulas, and thus the results and conclusions are unchanged

By employing the $1/N$ expansion, we compute the vacuum energy~$E(\delta\epsilon)$ of the two-dimensional supersymmetric (SUSY) $\mathbb{C}P^{N-1}$ model on~$\mathbb{R}\times S^1$ with $\mathbb{Z}_N$ twisted boundary conditions to the second order in a SUSY-breaking parameter~$\delta\epsilon$. This quantity was vigorously studied recently by Fujimori et\ al.\ using a semi-classical approximation based on the bion, motivated by a possible semi-classical picture on the infrared renormalon. In our calculation, we find that the parameter~$\delta\epsilon$ receives renormalization and, after this renormalization, the vacuum energy becomes ultraviolet finite. To the next-to-leading order of the $1/N$ expansion, we find that the vacuum energy normalized by the radius of the~$S^1$, $R$, $RE(\delta\epsilon)$ behaves as inverse powers of~$\Lambda R$ for~$\Lambda R$ small, where $\Lambda$ is the dynamical scale. Since $\Lambda$ is related to the renormalized 't~Hooft coupling~$\lambda_R$ as~$\Lambda\sim e^{-2\pi/\lambda_R}$, to the order of the $1/N$ expansion we work out, the vacuum energy is a purely non-perturbative quantity and has no well-defined weak coupling expansion in~$\lambda_R$., Comment: 24 pages, 3 figures. Fig. 2 was replaced by a correct one and corresponding descriptions in the text were corrected. A note added. The final version to appear in PTEP

We calculate the PCAC mass for $(2+1)$ flavor full QCD with Wilson-type quarks. We adopt the Small Flow-time eXpansion (SFtX) method based on the gradient flow which provides us a general way to compute correctly renormalized observables even if the relevant symmetries for the observable are broken explicitly due to the lattice regularization, such as the Poinc\'{a}re and chiral symmetries. Our calculation is performed on heavy $u, d$ quarks mass ($m_{\pi}/m_{\rho}\simeq0.63$) and approximately physical $s$ quark mass with fine lattice $a \simeq 0.07$~fm. The results are compared with those computed with the Schr\"odinger functional method., Comment: 7 pages, 6 figures, Talk presented at the 37th International Symposium on Lattice Field Theory (Lattice2019), 16-22 June 2019, Wuhan, China