Sigma resonance parameters from a N_{f}=2 lattice QCD simulation

Sigma resonance parameters from a lattice QCD simulation

R. Molina Institute of Physics of the University of São Paulo, Rua do Matão, 1371, Butantã, São Paulo, 05508-090, Brazil    D. Guo The George Washington University, Washington, DC 20052, USA    A. Alexandru The George Washington University, Washington, DC 20052, USA Department of Physics, University of Maryland, College Park, MD 20742, USA Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, CH-3012 Bern, Switzerland    M. Mai The George Washington University, Washington, DC 20052, USA    M. Döring The George Washington University, Washington, DC 20052, USA Thomas Jefferson National Accelerator Facility, Newport News, VA 23606, USA

In this work we present the analysis of the energy spectrum from a recent two-flavor () lattice QCD calculation for pion-pion scattering in the scalar, isoscalar channel (the -meson). The lattice simulation was performed for two quark masses corresponding to a pion mass of and . The -meson parameters are extracted using various parametrizations of the scattering amplitude. The results obtained from a chiral unitary parametrization are extrapolated to the physical point and read , where the uncertainties in the parentheses denote the stochastic and systematic ones. The behavior of the -meson parameters with increasing pion mass is discussed as well.
XIV International Workshop on Hadron Physics,
18-23 March 2018,
Florianópolis, SC, Brazil


=.08em \lightrulewidth=.05em \cmidrulewidth=.03em \belowrulesep=.65ex \belowbottomsep=0pt \aboverulesep=.4ex \abovetopsep=0pt \cmidrulesep=\cmidrulekern=.5em \defaultaddspace=.5em

I Introduction

The is one of the resonances more discussed in the literature. Several phenomenological models have pointed out that it is relevant for nuclear forces. However, its properties have been subject to a long discussion due to a variety of models which tried to determine its mass and width from experimental data, often showing incompatible results. Even its existence was several times questioned and the Review of Particle Properties did not show it for twenty years. The particle is nowadays a well-established resonance whose pole properties have been determined more precisely from dispersion relations and coupled-channel approaches which rely upon chiral Lagrangians, analyticity, and unitarity relations Pelaez:2015qba ().

Recently, the GWU group has extracted the energy spectrum and phase shifts for pion-pion scattering in the scalar, isoscalar channel, in a lattice QCD simulation for two different pion masses, and MeV. In this work, we present the analysis of these data by using several parameterizations and estimate the uncertainties due to the use of one or other. First, we consider a general expansion in an energy-variable conformally mapping the energy plane to the unit disk, similar to the analysis of Refs. Yndurain:2007qm (); Caprini:2008fc (). Second, we employ a model based on the chiral unitary approach (UChPT), used, e.g., in Refs. Hu:2017wli (); Hu:2016shf (); Guo:2016zos (). Subsequently, the UChPT amplitude is extrapolated to the physical point. Our final result, based on all lattice data presented here with and without the isovector channel data Guo:2016zos (), reads  MeV and agrees with the result of the most recent analysis of experimental data Pelaez:2015qba () within the quoted region. Additionally, the study of the pion mass dependence of resonance mass and coupling to the channel is studied in a broader range of , as it was done in previous works Doring:2016bdr (); Hanhart:2008mx ().

Ii Analysis of the scattering amplitude

Figure 1: Results of the individual fits, for  MeV or  MeV, to the lattice data in terms of using the conformal parametrization in two variants ([cm1] and [cm2]), and the chiral unitary approach ([chm2]).

The discrete energy spectrum is obtained from the finite volume scattering amplitude,


which depends explicitly on the form of the -matrix Doring:2012eu (). For boxes with asymmetry in the direction and in the rest frame,


and , being the momentum in the rest frame, , and is two pion loop-function written in the conventional dimensional regularization form Guo:2018zss (). The positions of poles in Eq. (1) determine the discrete energy spectrum on the lattice. The form of the function is not fixed by unitarity. In this work we use four versions of two different types of the -matrix to gauge the systematic uncertainty tied to a particular choice. We analyze data from different pion masses and channels individually and also simultaneously. The parameterizations, channels and pion masses employed and number of parameters (pars.) in the fits are:
- Individual fits, or MeV,

  • Conformal parameterizations, (2 pars.)

  • Chiral parameterization, and (4 pars.)

- Combined fits, for and MeV,

  • Chiral parameterization, (3 pars.)

  • Chiral parameterization, and (4 pars.).

For the conformal parameterizations, Eqs. (7) and (8), and (26) and (27) from Ref. Caprini:2008fc (), are used. These are named [cm1] and [cm2] respectively. The relation of in Ref. Caprini:2008fc () with the K-matrix in Eq. (1), is , with . For the chiral parameterizations, we employ the chiral unitary approach described in Ref. Oller:1998hw (). We call them [chm1] and [chm2] for fits, and analyses, respectively. The data corresponding to the -channel was analyzed previously in Ref.  Guo:2016zos (). The same set of data is used in this work. The above parameterizations and details of the fits are described in detail in Refs. Guo:2016zos (); Guo:2018zss (). Below, we summarize our main results.

Figure 2: Results of the fits to the lattice data using chiral parametrization as described in the main text. The fits are performed simultaneously to the data of both pion masses ( and MeV), including only the data in the channel [chm1], or the and channel [chm2].

Iii Results

The results of the fits are collected in Table I of Ref. Guo:2018zss (), which all pass the Pearson’s test with the total lying inside of the 80% confidence interval. We discuss first the results from individual fits. The phase shifts are depicted in Fig. 1. There is a clear overlap of the fits with the data. Only for the heavy pion mass and energies deep below threshold there is some discrepancy that could be significant. Note again that the phase-shifts are not fitted directly, but rather the energies extracted from lattice QCD. There is a good agreement between the phase shifts from both conformal parameterizations used in the channel [cm1] and [cm2], that we are interested here, and also with the [chm2] chiral individual fits. Note that we used the conformal parameterizations with two parameters. Since the chiral model for alone has three parameters, uncertainties are larger if we only perform individual fits. Indeed, the chiral parameterization contains LECs which appear in both, and channels. Thus, the fit of both channels, and , constrains better the LECs used.

We perform an analytical continuation to the complex energy plane. On the second Riemann sheet of this plane we determine the position () and residuum () of the -resonance pole. The results are collected in Table 1. In the three different parameterizations, and for individual fits at or MeV, there is a good agreement in the pole position, which lies above threshold, being a resonance in both cases.

For combined fits, the corresponding phase-shifts are depicted in Fig. 2. As shown in this picture, fitting the channel alone leads to larger uncertainties in the chiral unitary model, than in combination with the channel. Note that uncertainites are also larger in this channel for the raw lattice data. The best fit parameters are collected in Table I of Ref. Guo:2018zss (). For the low-energy constants we observe agreement between the individual and combined fits results. The extrapolated phase-shifts to the physical point are depicted in Fig. 3, showing a good agreement with the experimental data below energies of MeV. The results for the pole properties are collected in the rows 5-8 of Table 1. The pole position and residuum at the physical mass is compatible in both chiral parameterizations, [chm1] and [chm2]. There is also agreement in the pole properties for the other pion masses. Note that, in the case of the heavy pion mass, the pole position is pushed towards the threshold in the combined fit because of the presence of the light pion mass, becoming a virtual state. The uncertainites for the pole properties at this mass are larger, partly because the state is near the transition to a virtual state, and partly because of the larger uncertainites in the raw lattice data.

Figure 3: Phase-shift extrapolated to the physical point as a result of the fits to the energy levels in the (red) and the (blue) channels, in comparison with the experimental data, see Ref. Guo:2018zss ().
Parametrization Fitted data                           
cm1   –   –   –            –   –   –
cm1   –   –   –   –   –   –         
cm2   –   –   –            –   –   –
cm2   –   –   –   –   –   –         
Ref. Pelaez:2015qba () experimental            –   –   –   –   –   –
Table 1: Pole positions ( in MeV) and corresponding couplings to the channel ( in GeV) from conformal mapping ([cm1] and [cm2]) and chiral unitary approach ([chm1] and [chm2]).
Figure 4: Left: -dependence of the pole position of the -resonance in the complex plane of the center of mass three-momentum from the [chm2] fit to the data from both pion masses. The dashed line represents the real -axis, which connects the first () and second () Riemann sheets, and the subscript denotes the positive/negative half-plane, respectively. The encircled numbers represent the pion mass in units of the physical one, while shows the result of the simultaneous fit to and at light (orange) and heavy (dark green) pion masses with corresponding error areas. Right: -dependence of the coupling of the -resonance to the channel in the same color coding as in the left panel.

With this in mind we make a prediction of the pole position and the corresponding coupling to the channel as a continuous function of the pion mass based on [chm2] fitted to both sets of lattice data, and , and both pion masses simultaneously. To remind the reader these sets are obtained from calculations at and . The result of the extrapolation is depicted in Fig. 4. With increasing pion mass, both poles for positive and negative travel in the 2on Riemann sheet towards the threshold, coupling more strongly to this channel. For , both poles meet at the real energy axis below threshold on the second Riemann sheet becoming virtual bound states. For higher pion mass, the poles evolve on the real axis towards and away from the threshold (), respectively. One of the pole reaches the two-pion threshold at , where the coupling vanishes, and becomes a bound state for heavier pion masses, then, the coupling to two pions increases in this region monotonically. This behavior is remarkably similar to the one of Ref. Hanhart:2008mx ().

Iv Conclusions

To extract the parameters of the -resonance, we have to use a parameterization that satisfies physical constraints, in particular unitarity, analyticity, and proper chiral behavior. To gauge the systematics associated with the choice of such parametrizations, we used two types of approaches (each in several variants): a generic one that makes no assumption about the underlying dynamics, and a chiral perturbation theory inspired one that allows to extrapolate the resonance parameters to different (i.e., physical) pion mass. The systematic errors associated with the choice of parametrization are about 10% for the pole position.

One of the strengths of the chiral parametrization is that it allows us to fit simultaneously both the and channel, for both pion masses. We find that the model describes the data well and that the results extracted from the simultaneous fit to both channels agree well with the -channel fit results. We use the combined fit to extrapolate to the physical point and, based on the position of the pole in the complex energy plane, we find that . Here the first error is the stochastic error and the second one is the combined systematic error discussed above.

The extrapolation to the physical point agrees with the experimental phase-shifts and the pole mass and width of the is compatible with the result of recent analyses based on experimental data.


  • (1) J. R. Pelaez, Phys. Rept. 658, 1 (2016)
  • (2) F. J. Yndurain, R. Garcia-Martin and J. R. Pelaez, Phys. Rev. D 76, 074034 (2007)
  • (3) I. Caprini, Phys. Rev. D 77, 114019 (2008)
  • (4) B. Hu, R. Molina, M. Döring, M. Mai and A. Alexandru, Phys. Rev. D 96, no. 3, 034520 (2017)
  • (5) B. Hu, R. Molina, M. Döring and A. Alexandru, Phys. Rev. Lett. 117, no. 12, 122001 (2016)
  • (6) D. Guo, A. Alexandru, R. Molina and M. Döring, Phys. Rev. D 94, no. 3, 034501 (2016)
  • (7) M. Döring, B. Hu and M. Mai, arXiv:1610.10070 [hep-lat].
  • (8) C. Hanhart, J. R. Pelaez and G. Rios, Phys. Rev. Lett. 100, 152001 (2008)
  • (9) M. Doring, U. G. Meissner, E. Oset and A. Rusetsky, Eur. Phys. J. A 48, 114 (2012)
  • (10) D. Guo, A. Alexandru, R. Molina, M. Mai and M. Döring, arXiv:1803.02897 [hep-lat].
  • (11) J. A. Oller, E. Oset and J. R. Pelaez, Phys. Rev. D 59, 074001 (1999) Erratum: [Phys. Rev. D 60, 099906 (1999)] Erratum: [Phys. Rev. D 75, 099903 (2007)]
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description