No Arabic abstract
A systematic analysis of the structure of single-baryon correlation functions calculated with lattice QCD is performed, with a particular focus on characterizing the structure of the noise associated with quantum fluctuations. The signal-to-noise problem in these correlation functions is shown, as long suspected, to result from a sign problem. The log-magnitude and complex phase are found to be approximately described by normal and wrapped normal distributions respectively. Properties of circular statistics are used to understand the emergence of a large time noise region where standard energy measurements are unreliable. Power-law tails in the distribution of baryon correlation functions, associated with stable distributions and Levy flights, are found to play a central role in their time evolution. A new method of analyzing correlation functions is considered for which the signal-to-noise ratio of energy measurements is constant, rather than exponentially degrading, with increasing source-sink separation time. This new method includes an additional systematic uncertainty that can be removed by performing an extrapolation, and the signal-to-noise problem re-emerges in the statistics of this extrapolation. It is demonstrated that this new method allows accurate results for the nucleon mass to be extracted from the large-time noise region inaccessible to standard methods. The observations presented here are expected to apply to quantum Monte Carlo calculations more generally. Similar methods to those introduced here may lead to practical improvements in analysis of noisier systems.
Single state saturation of the temporal correlation function is a key condition to extract physical observables such as energies and matrix elements of hadrons from lattice QCD simulations. A method commonly employed to check the saturation is to seek for a plateau of the observables for large Euclidean time. Identifying the plateau in the cases having nearby states, however, is non-trivial and one may even be misled by a fake plateau. Such a situation takes place typically for the system with two or more baryons. In this study, we demonstrate explicitly the danger from a possible fake plateau in the temporal correlation functions mainly for two baryons ($XiXi$ and $NN$), and three and four baryons ($^3{rm He}$ and $^4{rm He})$ as well, employing (2+1)-flavor lattice QCD at $m_{pi}=0.51$ GeV on four lattice volumes with $L=$ 2.9, 3.6, 4.3 and 5.8 fm. Caution is given for drawing conclusion on the bound $NN$, $3N$ and $4N$ systems only based on the temporal correlation functions.
Path integrals describing quantum many-body systems can be calculated with Monte Carlo sampling techniques, but average quantities are often subject to signal-to-noise ratios that degrade exponentially with time. A phase-reweighting technique inspired by recent observations of random walk statistics in correlation functions is proposed that allows energy levels to be extracted from late-time correlation functions with time-independent signal-to-noise ratios. Phase reweighting effectively includes dynamical refinement of source magnitudes but introduces a bias associated with the phase. This bias can be removed by performing an extrapolation, but at the expense of re-introducing a signal-to-noise problem. Lattice Quantum Chromodynamics calculations of the $rho$ and nucleon masses and of the $XiXi$ binding energy show consistency between standard results obtained using earlier-time correlation functions and phase-reweighted results using late-time correlation functions inaccessible to standard statistical analysis methods.
We present first results for two-baryon correlation functions, computed using $N_f=2$ flavours of O($a$) improved Wilson quarks, with the aim of explaining potential dibaryon bound states, specifically the H-dibaryon. In particular, we use a GEVP to isolate the groundstate using two-baryon (hyperon-hyperon) correlation functions $big(langle C_{XY}(t)C_{XY}(0) rangle$, where $XY=LambdaLambda, SigmaSigma, NXi, cdotsbig)$, each of which has an overlap with the H-dibaryon. We employ a `blocking algorithm to handle the large number of contractions, which may easily be extended to N-baryon correlation functions. We also comment on its application to the analysis of single baryon masses ($n$, $Lambda$, $Xi$, $cdots$). This study is performed on an isotropic lattice with $m_pi = 460$ MeV, $m_pi L = 4.7$ and $a = 0.063$ fm.
We investigate -- as an alternative to usual Monte Carlo Renormalization Group methods -- the feasibility of extracting QCD beta-functions directly from a lattice analysis of correlations between the action and Wilson loops. We test this correlation technique numerically in four dimensional SU(2) gauge theory, on a 16^4 lattice at beta = 2.5 and find very promising results.
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.