1 Introduction

CIPANP2018-Leskovec July 19, 2019

A Lattice QCD study of the resonance

Luka Leskovec

in collaboration with:

Constantia Alexandrou, Stefan Meinel, John Negele, Srijit Paul, Marcus Petschlies, Andrew Pochinsky, Gumaro Rendon, and Sergey Syritsyn

Department of Physics

University of Arizona, AZ-85721, USA

We present a lattice QCD study of the resonance with clover fermions at a pion mass of approximately MeV and lattice size fm. We consider two processes involving the . The first process is elastic scattering of two pions in P-wave with isospin . Using the Lüscher method we determine the scattering phase shift, from which we obtain the resonance mass and decay width . The second process is the radiative transition , where we follow the Briceño-Hansen-Walker-Loud approach to determine the transition amplitude in the invariant mass region near the resonance and for both space- and time-like photon momentum. This allows us to determine the coupling between the , the pion and the photon, and the resulting radiative decay width.


Thirteenth Conference on the Intersections of Particle and Nuclear Physics

Palm Springs CA, USA, May 29 – June 3, 2018

1 Introduction

The spectrum of hadrons is a jungle of particles arising from the strong interactions between quarks and gluons. We differentiate between stable hadrons, i.e. those that cannot decay via the strong interactions, and unstable hadrons, i.e. those that can. The simplest example of an unstable hadron is the meson, which is an isotriplet with quantum numbers . While the couples to multiple decay modes, for example , , and , we focus only on the latter two in this report. This simplification can be made because we perform our lattice QCD calculation at light quark masses corresponding to MeV, where the and thresholds are above the energy region we are focusing on. The coupling of the to the channel was previously investigated with lattice QCD in Refs. [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14], and the coupling to in Refs. [15, 16, 17].

To determine the coupling of the resonance to two pions and thus establish its strong decay width, we make use of the Lüscher formalism, which relates the finite-volume spectrum in various moving frames and irreducible representations with the infinite-volume scattering matrix [18, 19]. The appropriate formalism to handle radiative transitions of multi-hadron states, which we employ here, is the generalization of Lellouch and Lüscher’s work in Ref. [20] by Briceño, Hansen and Walker-Loud (BHWL) [21, 22].

2 Lattice Setup

We use a single lattice gauge-field ensemble with light quark masses corresponding to MeV and a strange-quark mass consistent with its physical value. The number of lattice points is . The lattice spacing , determined from the splitting calculated with NRQCD, is equal to fm, leading to a physical spatial volume of approximately fm.

3 -wave scattering in

We calculated the lattice spectra in several moving frames and irreducible representations [13] and determined the scattering phase shifts using the Lüscher method. To describe the phase shifts we use the Breit-Wigner formula


where is the invariant mass, is the -wave scattering phase shift, is the resonance mass and is the decay width. We investigate two different parametrizations:

  • BW I: -wave decay width:


    where is the coupling between the resonance and the scattering channel and is the scattering momentum, , and

  • BW II: -wave decay width modified with Blatt-Weisskopf barrier factors [23]:


    where is the scattering momentum at and is the centrifugal barrier radius.

  Model BW I BW II   Model BW I BW II

Figure 2: -wave scattering phase shift in as a function of invariant mass.
Figure 1: -wave scattering phase shift in as a function of invariant mass.
Table 3: Parameters of the Breit-Wigner parametrizations BW I and BW II obtained from fitting the discrete phase shift points with statistical uncertainties only [13].
Figure 1: -wave scattering phase shift in as a function of invariant mass.
Table 2: Parameters of the Breit-Wigner parametrizations BW I and BW II obtained from fitting the discrete phase shift points with statistical uncertainties only [13].

The resulting fits for both parametrizations are shown in Fig. 2, where the blue line corresponds to BW I and the red line to BW II. The numerical results for the parameters , and are listed in Table 2. There are only minor and statistically not very significant differences between the two parametrizations that appear in the high-energy region , where BW II describes the data slightly better [13]. Overall, we find that both parametrizations describe our results for the elastic -wave scattering well.

4 transition amplitude

Because the is not a QCD asymptotic state, but rather a resonance in -wave scattering with , the observables related to the resonance photoproduction processes are obtained from the more general process . This process is described by the transition amplitude , which is a function of both the photon four-momentum transfer and the invariant mass . The transition can however be defined by analytical continuation to the pole located at .

Figure 3: 3-D plot of the transition amplitude . Lattice QCD results are shown as the vertical bars, where the widths and depths correspond to the uncertainties in and , and the magenta sections at the tops cover the range from to . For clarity, the opacity is reduced in data points with larger uncertainties. The surface shows the central value of the chosen parametrization “BWII F1 K2”.

The transition amplitude is associated with the infinite volume matrix element by


where is the photon four-momentum transfer. The transition amplitude has a manifest pole at , and can be written with the help of Watson’s theorem as


The form factor is thus free of poles in within the energy region of interest, . Following a general Taylor-expansion approach we parametrize using


where is the pole in the -channel, and the two variables and are defined as [24, 25, 26]:


We consider three different families of truncations of the series in Eq. 6 leading to several parametrizations of the transition amplitude [17]. A 3-D representation of the transition amplitude is shown for the chosen parametrization “BW II F1 K2” in Fig. 3.

In practice, the resonant form factor is determined by evaluating at the pole:


It becomes equal to the photocoupling at zero momentum transfer, . The physical observable we consider is the radiative decay width determined by the photocoupling ,


The resonant transition form factor is shown in the left panel of Fig. 4, where the inner shaded region represents the statistical and systematical uncertainties determined on the lattice. The outer shaded region indicates the parametrization uncertainty. The photocouplings determined for each of the parametrizations are shown in the right panel of Fig. 4; we find the photocoupling at MeV to be

Figure 4: Left panel: Resonant transition form factor . The inner shaded region represents the statistical and systematical uncertainties, and the outer shaded region represents the parametrization uncertainty. Right panel: Value of the photocoupling for the parametrizations under consideration. The data-point errorbars are associated with lattice statistical and systematical uncertainties. The inner shaded region corresponds to the errorbar of the chosen model, “BW II F1 K2”, while the outer shaded region represents the parametrization uncertainty.

Due to the larger-than physical pion mass in our calculation the thresholds are much closer to the resonances than in nature. By assuming that the photocoupling is a pion-mass-independent quantity, we use the physical values of the pion masses and mass to determine the radiative decay width to be . The number in the first bracket is the combined statistical and systematic uncertainty and the number in the second bracket is the parametrization uncertainty.

5 Summary

We have presented results of two of our recent lattice QCD studies, the determination of the strong decay width of the resonance and the calculation of the resonance radiative decay width. While our calculation was performed at a non-physical value of the quark mass, the couplings and are already close to their physical values. Future studies are needed to perform chiral extrapolations to the physical point [27, 28, 29] where direct comparisons to experiment can be made.


We are grateful to Kostas Orginos for providing the gauge field ensemble, which was generated using resources provided by XSEDE (supported by National Science Foundation Grant No. ACI-1053575). We thank R. A. Briceño, M. Hansen, and C. B. Lang for valuable discussions. SM and GR were supported in part by National Science Foundation Grant No. PHY-1520996; SM, GR, and LL were also supported in part by the U.S. Department of Energy Office of High Energy Physics under Grant No. DE-SC0009913. SM and SS futher acknowledge support by the RHIC Physics Fellow Program of the RIKEN BNL Research Center. JN and AP were supported in part by the U.S. Department of Energy Office of Nuclear Physics under Grant Nos. DE-SC-0011090 and DE-FC02-06ER41444. We acknowledge funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 642069. SP is a Marie Sklodowska-Curie fellow supported by the HPC-LEAP joint doctorate program. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under Contract No. DE-AC02-05CH11231. The computations were performed using the Qlua software suite [30].


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