Hyperonnucleon potentials from lattice QCD
Abstract
We calculate potentials from the equaltime BetheSalpeter amplitude measured in the quenched QCD simulation with the spatial lattice volume, (4.4 fm). The standard Wilson gauge action with the gauge coupling on lattice together with the standard Wilson quark action are used. The hopping parameter is chosen for and quarks, which corresponds to GeV. The physical strange quark mass is used by taking the parameter which is deduced from the physical meson mass. The lattice spacing fm is determined by the physical meson mass. We find that the potential has strong spin dependence. Strong repulsive core is found in channel while the effective central potential in the channel has relatively weak repulsive core. The potentials also have weak attractive parts in the medium to long distance region ( fm fm) in both of the and channels.
Hyperonnucleon potentials from lattice QCD
Tetsuo Hatsuda
Department of Physics, University of Tokyo, Tokyo 1130033, Japan
Email: hatsuda@phys.s.utokyo.ac.jp
\abstract@cs
1 Introduction
Study of hypernuclei is one of the frontiers in nuclear physics. The strange degree of freedom gives a new dimension to the description of nuclear structure. A hyperon (or a strange quark) embedded in nuclei plays a characteristic role as an “impurity” or a probe in the many body system[1].
Modern nucleonnucleon () potentials give a successful description of the scattering data and have been used to make precise calculations of light nuclei[2, 3]. In contrast, hyperonnucleon () and hyperonhyperon () interactions have large uncertainties, because the scattering experiments are either difficult or impossible due to the short lifetime of hyperons. The full phase shift analysis at the same level of the case is not available yet. Although a lot of theoretical models which describe interaction together with the interaction have been published, the different model predicts different phase shifts and scattering parameters, e.g., of the potential in spite of the nice description of the sector[4, 5, 6, 7, 8, 9].
Recently, a lattice QCD study of the potential has been performed[10]. This approach may lead to a new paradigm to study the and interactions too, since the lattice QCD is an ab initio method of treating the fundamental theory of strong interaction. (See also the conventional approach to the phase shifts using the Lüscher’s finite volume formula[11].)
The purpose of the present report is to explore the and potentials from the lattice QCD simulation on the basis of the methodology developed in Refs. [10, 12]. The extension from the potential to and potentials is relatively straightforward. In the case of the potential, there are only two representations in the isospin channel, i.e., in , which correspond to isovector and isoscalar channels. Including the strange degrees of freedom extends the arithmetic into flavor , Here the isovector (isoscalar) channel of the sector is assigned to be a subset in the plet (plet) representation. The potentials for newly arising channels are hardly determined from the real experiment so far. The lattice QCD simulation with physical strange quark mass provides new numerical “data” for these strange channels.
In this report, we focus on the potential in the isovector () channel as a first step. We keep away, at present, from the isoscalar () channel of the potential, since it is not the lowest state of the isoscalar channel and strong decay mode may open below the . There are almost no experimental information on the interaction, although a few experimental data[13, 14, 15] suggest that the nucleus potential would be weakly attractive. Moreover, the nucleus interaction will be studied as a day one experiment in the near future at JPARC[16] through reaction on the nuclear target such as C.
2 Formulation
As we mentioned, the methodology to obtain the potential is along the lines of Refs. [12] and [10]. The latter describes, successfully for the first time, the potentials from lattice QCD simulation. We start from an effective Schrödinger equation for system at low energies:
(2.0) 
where and are the reduced mass of the system and the nonrelativistic energy in the centerofmass frame, respectively. The nonlocal potential can be represented by the derivative expansion as
(2.0) 
The general expression of the potential would be
(2.0)  
Here is the tensor operator with , symmetric () and antisymmetric () spin operators, the relative angular momentum operator, and () is isospin operator for (). We note that the antisymmetric spinorbit forces ( and ) newly come up because the constituents ( and ) are not identical.
According to the above expansion of the potential, the wave function should be classified by the total isospin , the total angular momentum and parity with . Particular spin (isospin) projection can be made in terms of (), e.g., for the isospin projection we have and . In this work, we focus only the isospin , wave component of the wave function, so as to obtain the (effective) central potential through
(2.0) 
The wave wave function is measured from the equaltime BetheSalpeter (BS) amplitude as
(2.0)  
(2.0)  
(2.0) 
where and denote the Dirac indices, , and the color indices, and the charge conjugation matrix. The summation over is taken for cubic transformation group to project onto the wave, and the summation over for zero total spatial momentum. and are the local field operators for the proton and . We take the upper components of the Dirac indices and to construct the spin singlet (triplet) channel by ().
The with is understood as the probability amplitude to find “nucleonlike” three quarks located at point and “like” three quarks located at point . The includes not only the elastic amplitude but also the inelastic amplitudes such as and , and so on. Note that, at low energies below the thresholds, however, the asymptotic behavior of is not affected by the inelastic contributions, since they decrease exponentially in the asymptotic region. (In the present calculation with the isospin , the channel is closed due to the isospin conservation.) On the other hand, and hence the potential may depend on the interpolating fields in the interaction region. Further study on this issue is found in Ref.[17] for the potential.
In the actual simulations, the BS amplitude is obtained through the fourpoint correlator,
(2.0)  
(2.0) 
Here is a source term located at . We utilize the wall source for in this work in order to enhance the lowest scattering state of the system. is the energy of the state, , and .
We also calculate the twopoint correlator, , for the octet baryons (), in order to check whether various two baryon (,, and ) thresholds are reproduced in the correct order. The interpolating fields for and , employed in this work, are given by
(2.0) 
(2.0) 
3 Numerical calculation
We use the standard Wilson gauge action at the gauge coupling on the lattice together with the standard Wilson quark action. See Ref.[17] for details. The hopping parameter of is chosen for the and quarks, which corresponds to GeV, GeV, and GeV. In order to determine the parameter for the strange quark mass (), we first measure the correlators of pseudo scalar and vector mesons by using the interpolating fields given by
(3.0)  
(3.0) 
with applying six set of hopping parameters; first three sets for ( and mesons) with taking a number from , and another three sets for ( and mesons) with taking different two numbers from . Assuming the following functional forms for the pseudo scalar meson mass squared and for the vector meson mass,
(3.0)  
(3.0) 
we obtained critical hopping parameter , and the physical parameter, , from . The lattice scale is determined so as to be fm from the physical meson mass. The parameter for the strange quark mass is determined as from the physical meson mass (494 MeV).
4 Results and Discussion
367(1)  811(4)  552.6(5)  882(2)  1164(7)  1379(6)  1263(6)  1312(6) 
Table 1 lists the hadron masses measured from the present lattice QCD simulation. 1283 gauge configurations are used to calculate the hadron masses and wave functions. (17 exceptional configurations are not used out of totally 1300 gauge configurations.) We note that the present results for the baryon masses provide the correct order of particular threshold energy of two baryon states in the strangeness sector; MeV (this channel is not allowed in the present case because of isospin conservation), MeV, MeV, and MeV. This warrants the desirable asymptotic behavior of the wave function.
Figure 1 shows the wave function obtained at the time slice . The () channel is plotted by circles (triangles), which are normalized at the spatial boundary . All the data are taken into account for fm, while only the data on the , , and axis and their nearest neighbors are used to plot for the outer region. As seen in the Figure, the wave functions are suppressed in the short distance region, and a slight enhancement is found in the medium range region for both and channels. There is sizable difference between the and , particularly of suppression in the short distance, suggesting that the channel has stronger repulsive core.
Figure 2 shows the (effective) central potentials for in the and channels. These results are still preliminary, since the potentials are obtained by assuming : The energy should be determined by fitting the asymptotic behavior of the wave function with the use of the Green’s function which is a solution of the Helmholtz equation on the lattice[12, 18]. Preliminary calculation suggests that would be small negative values for both and channels, similar to the case of the potential[10, 19].
In order to see the ground state saturation of the present results, we plot, in Fig. 3, the timeslice dependence of the potential in both of the (lefthandside) and (righthandside) channels at several radial distances; , , , , and fm. We can see that the saturation is achieved for within errors.
The present work is a first step toward the and potentials from the lattice QCD simulation. Systematic studies of the various channels such as , , , and so on are all interesting and important because they are intimately related not only to the structure of hypernuclei but also to the internal structure of neutron stars. We will present such studies in the near future.
Acknowledgments
Lattice QCD Monte Carlo calculation has been done with IBM Blue Gene/L computer at KEK. H. N. is supported by the Special Postdoctoral Researchers Program at RIKEN. This research was partly supported by GrantsinAid for Young Scientists (B) (No. 17740174) from the Japan Society for Promotion of Science (JSPS), and by the Ministry of Education, Science, Sports and Culture, GrantinAid (Nos. 13135204, 15540251, 15540254, 18540253, 19540261).
References
 [1] H. Tamura, Eur. Phys. J. A 13, 181 (2002).
 [2] A. Nogga, H. Kamada, W. Glöckle, B.R. Barrett, Phys. Rev. C 65, 054003 (2002).
 [3] S.C. Pieper, K. Varga, R.B. Wiringa, Phys. Rev. C 66, 044310 (2002).
 [4] Th.A. Rijken, Y. Yamamoto, Phys. Rev. C 73 044008 (2006).
 [5] Th.A. Rijken, V.G.J. Stoks, Y. Yamamoto, Phys. Rev. C 59, 21 (1999).
 [6] P.M.M. Maessen, Th.A. Rijken, J.J. de Swart, Phys. Rev. C 40, 2226 (1989).
 [7] Y. Fujiwara, C. Nakamoto, Y. Suzuki, Phys. Rev. C 54, 2180 (1996).
 [8] Y. Fujiwara, M. Kohno, C. Nakamoto, Y. Suzuki, Phys. Rev. C 64, 054001 (2001).
 [9] I. Arisaka, K. Nakagawa, M. Wada, S. Shinmura, Prog. Theor. Phys. 113, 1287 (2005).
 [10] N. Ishii, S. Aoki, T. Hatsuda, Phys. Rev. Lett. 99, 022001 (2007).
 [11] S. R. Beane et al. [NPLQCD Collab.], Nucl. Phys. A 794, 62 (2007).
 [12] S. Aoki, et al. [CPPACS Collab.], Phys. Rev. D71, 094504 (2005).
 [13] K. Nakazawa, Nucl. Phys. A 639, C345 (1998).
 [14] T. Fukuda, et al., Phys. Rev. C58, 1306 (1998).
 [15] P. Khaustov, et al., Phys. Rev. C 61, 054603 (2000).

[16]
See
for example,
http://jparc.jp/NuclPart/index_e.html
 [17] N. Ishii, S. Aoki, T. Hatsuda, in preparation.
 [18] M. Lüscher, Nucl. Phys. B 354, 531 (1991).
 [19] N. Ishii, S. Aoki, T. Hatsuda, in these proceedings.