No Arabic abstract
The approximated partial wave decomposition method to the discrete data on a cubic lattice, developed by C. W. Misner, is applied to the calculation of $S$-wave hadron-hadron scatterings by the HAL QCD method in lattice QCD. We consider the Nambu-Bethe-Salpeter (NBS) wave function for the spin-singlet $Lambda_c N$ system calculated in the $(2+1)$-flavor QCD on a $(32a~mathrm{fm})^3$ lattice at the lattice spacing $asimeq0.0907$ fm and $m_pi simeq 700$ MeV. We find that the $l=0$ component can be successfully extracted by Misners method from the NBS wave function projected to $A_1^+$ representation of the cubic group, which contains small $lge 4$ components. Furthermore, while the higher partial wave components are enhanced so as to produce significant comb-like structures in the conventional HAL QCD potential if the Laplacian approximated by the usual second order difference is applied to the NBS wave function, such structures are found to be absent in the potential extracted by Misners method, where the Laplacian can be evaluated analytically for each partial wave component. Despite the difference in the potentials, two methods give almost identical results on the central values and on the magnitude of statistical errors for the fits of the potentials, and consequently on the scattering phase shifts. This indicates not only that Misners method works well in lattice QCD with the HAL QCD method but also that the contaminations from higher partial waves in the study of $S$-wave scatterings are well under control even in the conventional HAL QCD method. It will be of interest to study interactions in higher partial wave channels in the HAL QCD method with Misners decomposition, where the utility of this new technique may become clearer.
The reweighting method is widely used in numerical studies of QCD, in particular, for the cases in which the conventional Monte-Carlo method cannot be applied directly, e.g., finite density QCD. However, the application range of the reweighing method is restricted due to several problems. One of the most severe problems here is the overlap problem. To solve it, we examine a multipoint reweighting method in which simulations at several simulation points are combined in the data analyses. We systematically study the applicability and limitation of the multipoint reweighting method in two-flavor QCD at zero density. Measuring histograms of physical quantities at a series of simulation points, we apply the multipoint reweighting method to calculate the meson masses as continuous functions of the gauge coupling $beta$ and the hopping parameters $kappa$. We then determine lines of constant physics and beta functions, which are needed in a calculation of the equation of state at finite temperature.
In this article, we review the HAL QCD method to investigate baryon-baryon interactions such as nuclear forces in lattice QCD. We first explain our strategy in detail to investigate baryon-baryon interactions by defining potentials in field theories such as QCD. We introduce the Nambu-Bethe-Salpeter (NBS) wave functions in QCD for two baryons below the inelastic threshold. We then define the potential from NBS wave functions in terms of the derivative expansion, which is shown to reproduce the scattering phase shifts correctly below the inelastic threshold. Using this definition, we formulate a method to extract the potential in lattice QCD. Secondly, we discuss pros and cons of the HAL QCD method, by comparing it with the conventional method, where one directly extracts the scattering phase shifts from the finite volume energies through the Luschers formula. We give several theoretical and numerical evidences that the conventional method combined with the naive plateau fitting for the finite volume energies in the literature so far fails to work on baryon-baryon interactions due to contaminations of elastic excited states. On the other hand, we show that such a serious problem can be avoided in the HAL QCD method by defining the potential in an energy-independent way. We also discuss systematics of the HAL QCD method, in particular errors associated with a truncation of the derivative expansion. Thirdly, we present several results obtained from the HAL QCD method, which include (central) nuclear force, tensor force, spin-orbital force, and three nucleon force. We finally show the latest results calculated at the nearly physical pion mass, $m_pi simeq 146$ MeV, including hyperon forces which lead to form $OmegaOmega$ and $NOmega$ dibaryons.
We make a detailed comparison between the direct method and the HAL QCD potential method for the baryon-baryon interactions, taking the $XiXi$ system at $m_pi= 0.51$ GeV in 2+1 flavor QCD and using both smeared and wall quark sources. The energy shift $Delta E_mathrm{eff}(t)$ in the direct method shows the strong dependence on the choice of quark source operators, which means that the results with either (or both) source are false. The time-dependent HAL QCD method, on the other hand, gives the quark source independent $XiXi$ potential, thanks to the derivative expansion of the potential, which absorbs the source dependence to the next leading order correction. The HAL QCD potential predicts the absence of the bound state in the $XiXi$($^1$S$_0$) channel at $m_pi= 0.51$ GeV, which is also confirmed by the volume dependence of finite volume energy from the potential. We also demonstrate that the origin of the fake plateau in the effective energy shift $Delta E_mathrm{eff}(t)$ at $t sim 1$ fm can be clarified by a few low-lying eigenfunctions and eigenvalues on the finite volume derived from the HAL QCD potential, which implies that the ground state saturation of $XiXi$($^1$S$_0$) requires $t sim 10$ fm in the direct method for the smeared source on $(4.3 mathrm{fm})^3$ lattice, while the HAL QCD method does not suffer from such a problem.
A comparative study between the Luschers finite volume method and the time-dependent HAL QCD method is given for the $XiXi$($^1mathrm{S}_0$) interaction as an illustrative example. By employing the smeared source and the wall source for the interpolating operators, we show that the effective energy shifts $Delta E_{rm eff} (t)$ in Luschers method do not agree between different sources, yet both exhibit fake plateaux. On the other hand, the interaction kernels $V(vec{r})$ obtained from the two sources in the HAL QCD method agree with each other already for modest values of $t$. We show that the energy eigenvalues $Delta E(L)$ in finite lattice volumes ($L^3$) calculated by $V(vec{r})$ indicate that there is no bound state in the $XiXi(^1mathrm{S}_0)$ channel at $m_{pi}=0.51$ GeV in 2+1 flavor QCD.
In this paper, employing an all-to-all quark propagator technique, we investigate the kaon-nucleon interactions in lattice QCD. We calculate the S-wave kaon-nucleon potentials at the leading order in the derivative expansion in the time-dependent HAL QCD method, using (2+1)-flavor gauge configurations at the lattice spacing $a approx 0.09$ fm on $32^3 times 64$ lattices and the pion mass $m_{pi} approx 570$ MeV. We take the one-end trick for all-to-all propagators, which allows us to put the zero momentum hadron operators at both source and sink and to smear quark operators at the source. We find the stronger repulsive interaction in the $I=1$ channel than in the $I=0$. The phase shifts obtained by solving the Schr{o}dinger equations with the potentials qualitatively reproduce the energy dependence of the experimental phase shifts, and have the similar behavior to the previous results from lattice QCD without all-to-all propagators. Our study demonstrates that the all-to-all quark propagator technique with the one-end trick is useful to study interactions for meson-baryon systems in the HAL QCD method, so that we will apply it to meson-baryon systems which contain quark-antiquark creation/annihilation processes in our future studies.