Search

We apply the gradient flow on a color-electric two-point function that encodes the heavy quark momentum diffusion coefficient. The simulations are done on fine isotropic lattices in the quenched approximation at $1.5\,T_c$. The continuum extrapolation is performed at fixed flow time followed by a second extrapolation to zero flow time. Perturbative calculations of this correlation function under Wilson flow are used to enhance the extrapolations of the non-perturbative lattice correlator. The final estimate for the continuum correlator at zero flow time largely agrees with one obtained from a previous study using the multi-level algorithm. We perform a spectral reconstruction based on perturbative model fits to estimate the heavy quark momentum diffusion coefficient. The approach we present here yields high-precision data for the correlator and is also applicable for actions with dynamical fermions., Published version, 14 pages, 15 figures

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.

Motivated by the recent work of Chigusa, Moroi, and Shoji, we propose a new simple gradient flow equation to derive the bounce solution which contributes to the decay of the false vacuum. Our discussion utilizes the discussion of Coleman, Glaser, and Martin and we solve a minimization problem of the kinetic energy while fixing the potential energy. The bounce solution is derived as a scale-transformed of the solution of this problem. We also show that the convergence of our method is robust against a choice of the initial configuration., 6 pages, 8 figures; the version published in PRD

We study the pressure anisotropy in anisotropic finite-size systems in SU(3) Yang-Mills theory at nonzero temperature. Lattice simulations are performed on lattices with anisotropic spatial volumes with periodic boundary conditions. The energy-momentum tensor defined through the gradient flow is used for the analysis of the stress tensor on the lattice. We find that a clear finite-size effect in the pressure anisotropy is observed only at a significantly shorter spatial extent compared with the free scalar theory, even when accounting for a rather large mass in the latter., 10 pages, 5 figures, one figure added in Fig.4

We calculate the topological susceptibility at $2.5\text{ }\text{ }{T}_{\mathrm{c}}$ and $4.1\text{ }\text{ }{T}_{\mathrm{c}}$ in SU(3) pure Yang-Mills theory. We define topology with the help of gradient flow and we largely overcome the problem of poor statistics at high temperatures by applying a reweighting technique in terms of the topological charge, measured after a specific small amount of gradient flow. This allows us to obtain a sample of configurations which compares topological sectors with good statistics, with enhanced tunneling between topologies. We quote continuum extrapolated results at these two temperatures and conclude that our method is viable and can be extended without new conceptual problems to the case of full QCD with fermions.

It is (or should be) well known that specification of a heat bath requires both a temperature and a 4-velocity, the rest frame of the heat bath. In static spacetimes there is a very natural and unique candidate for the 4-velocity of the heat bath, the normalized timelike Killing vector. However in stationary non-static spacetimes the situation is considerably more subtle, and several different "natural" 4-velocity fields suitable for characterizing the rest frame of a heat bath can be defined - thus Buchdahl's 1949 analysis for the Tolman temperature gradient in a stationary spacetime is only part of the story. In particular, the heat bath most suitable for describing the Hawking radiation from a rotating black hole is best described in terms of a gradient flow normal to the spacelike hypersurfaces, not in terms of Killing vectors., 7 pages

Euclidean two-point correlators of the energy-momentum tensor (EMT) in SU(3) gauge theory on the lattice are studied on the basis of the Yang-Mills gradient flow. The entropy density and the specific heat obtained from the two-point correlators are shown to be in good agreement with those from the one-point functions of EMT. These results constitute a first step toward the first principle simulations of the transport coefficients with the gradient flow., 5 pages, 12 figures

We propose a new approach to circumvent the sign problem in which the integration path is optimized to control the sign problem. We give a trial function specifying the integration path in the complex plane and tune it to optimize the cost function which represents the seriousness of the sign problem. We call it the path optimization method. In this method, we do not need to solve the gradient flow required in the Lefschetz-thimble method and then the construction of the integration-path contour arrives at the optimization problem where several efficient methods can be applied. In a simple model with a serious sign problem, the path optimization method is demonstrated to work well; the residual sign problem is resolved and precise results can be obtained even in the region where the global sign problem is serious., 4 pages, 6 figures

It has become customary to use a smoothing algorithm called "gradient flow" to fix the lattice spacing in a simulation, through a parameter called $t_0$. It is shown that in order to keep the length $t_0$ fixed with respect to mesonic or gluonic observables as the number of colors $N_c$ is varied, the fiducial point for the flow parameter must be scaled nearly linearly in $N_c$. In simulations with dynamical fermions, the dependence of $t_0$ on the pseudoscalar meson mass flattens as the number of colors rises, in a way which is consistent with large $N_c$ expectations., Comment: 13 pages, 5 figures, published version of paper

In this work we study the effects of gradient flow and open boundary condition in the temporal direction in 1+1 dimensional lattice $\phi^4$ theory. Simulations are performed with periodic (PBC) and open (OPEN) boundary conditions in the temporal direction. The Effects of gradient flow and open boundary on the field $\phi$ and the susceptibility are studied in detail along with the finite size scaling analysis. In both cases, at a given volume, the phase transition point is shifted towards a lower value of lattice coupling $\lambda_0$ for fixed $m_0^2$ in the case of OPEN as compared to PBC with this shift found to be diminishing as volume increases. We compare and contrast the extraction of the boson mass from the two point function (PBC) and the one point function (OPEN) as the coupling, starting from moderate values, approaches the critical value corresponding to the vanishing of the mass gap. In the critical region, boundary artifacts become dominant in the latter. Our studies point towards the need for a detailed finite volume (scaling) analysis of the effects of OPEN in the critical region.

Using a finite volume gradient flow renormalization scheme with Schr\"odinger Functional boundary conditions, we compute the nonperturbative running coupling in the range $2.2\ensuremath{\lesssim}{\overline{g}}_{\mathrm{GF}}^{2}(L)\ensuremath{\lesssim}13$. Careful continuum extrapolations turn out to be crucial to reach our high accuracy. The running of the coupling is always between one loop and two loop and very close to one loop in the region of $200\text{ }\text{ }\mathrm{MeV}\ensuremath{\lesssim}\ensuremath{\mu}=\phantom{\rule{0ex}{0ex}}1/L\ensuremath{\lesssim}4\text{ }\text{ }\mathrm{GeV}$. While there is no convincing contact to two-loop running, we match nonperturbatively to the Schr\"odinger functional coupling with background field. In this case, we know the $\ensuremath{\mu}$-dependence up to $\ensuremath{\sim}100\text{ }\text{ }\mathrm{GeV}$ and can thus connect to the $\mathrm{\ensuremath{\Lambda}}$-parameter.

We present results for the equation of state in ($2+1$)-flavor QCD using the highly improved staggered quark action and lattices with temporal extent ${N}_{\ensuremath{\tau}}=6$, 8, 10, and 12. We show that these data can be reliably extrapolated to the continuum limit and obtain a number of thermodynamic quantities and the speed of sound in the temperature range 130--400 MeV. We compare our results with previous calculations and provide an analytic parameterization of the pressure, from which other thermodynamic quantities can be calculated, for use in phenomenology. We show that the energy density in the crossover region, $145\text{ }\text{ }\mathrm{MeV}\ensuremath{\le}T\ensuremath{\le}163\text{ }\text{ }\mathrm{MeV}$, defined by the chiral transition, is ${\ensuremath{\epsilon}}_{c}=(0.18--0.5)\text{ }\text{ }\mathrm{GeV}/{\mathrm{fm}}^{3}$, i.e., $(1.2--3.1)\text{ }{\ensuremath{\epsilon}}_{\text{nuclear}}$. At high temperatures, we compare our results with resummed and dimensionally reduced perturbation theory calculations. As a byproduct of our analyses, we obtain the values of the scale parameters ${r}_{0}$ from the static quark potential and ${w}_{0}$ from the gradient flow.

We compute the topological charge and its susceptibility in finite temperature (2+1)-flavor QCD on the lattice applying a gradient flow method. With the Iwasaki gauge action and nonperturbatively $O(a)$-improved Wilson quarks, we perform simulations on a fine lattice with~$a\simeq0.07\,\mathrm{fm}$ at a heavy $u$, $d$ quark mass with $m_\pi/m_\rho\simeq0.63$ but approximately physical $s$ quark mass with $m_{\eta_{ss}}/m_\phi\simeq0.74$. In a temperature range from~$T\simeq174\,\mathrm{MeV}$ ($N_t=16$) to $697\,\mathrm{MeV}$ ($N_t=4$), we study two topics on the topological susceptibility. One is a comparison of gluonic and fermionic definitions of the topological susceptibility. Because the two definitions are related by chiral Ward-Takahashi identities, their equivalence is not trivial for lattice quarks which violate the chiral symmetry explicitly at finite lattice spacings. The gradient flow method enables us to compute them without being bothered by the chiral violation. We find a good agreement between the two definitions with Wilson quarks. The other is a comparison with a prediction of the dilute instanton gas approximation, which is relevant in a study of axions as a candidate of the dark matter in the evolution of the Universe. We find that the topological susceptibility shows a decrease in $T$ which is consistent with the predicted $\chi_\mathrm{t}(T) \propto (T/T_{\rm pc})^{-8}$ for three-flavor QCD even at low temperature $T_{\rm pc} < T\le1.5 T_{\rm pc}$., Comment: 18 pages, 4 figures, publised version