<<

ADP-15-33/T935

Flavor-singlet in the graded approach to partially quenched QCD

Jonathan M. M. Hall and Derek B. Leinweber Special Research Centre for the Subatomic Structure of (CSSM), Department of Physics, University of Adelaide, Adelaide, South Australia 5005, Australia

Progress in the calculation of the electromagnetic properties of excitations in lattice QCD presents new challenges in the determination of sea- loop contributions to matrix elements. A reliable estimation of the sea-quark loop contributions represents a pressing issue in the accurate comparison of lattice QCD results with experiment. In this article, an extension of the graded symmetry approach to partially quenched QCD is presented, which builds on previous theory by explicitly including flavor-singlet baryons in its construction. The formalism takes into account the interactions among both octet and singlet baryons, octet , and their ghost counterparts; the latter enables the isolation of the quark-flow disconnected sea-quark loop contributions. The introduction of flavor-singlet states enables systematic studies of the internal structure of Λ-baryon excitations in lattice QCD, including the topical Λ(1405).

PACS numbers: 12.39.Fe 12.38.Gc 14.20.-c 11.30.Rd 12.38.Aw

I. INTRODUCTION

A. Quark-loop contributions to matrix elements

The calculation of hadronic matrix elements in nonper- turbative QCD typically encounters non-trivial contributions from quark-flow disconnected diagrams in which the current interacts with the of interest, first through a virtual sea-quark loop. An example process for a baryon undergo- FIG. 1. The quark-flow disconnected diagram for the process B → ing a virtual -baryon transition, B → B Φ, involving B Φ. The cross on the disconnected quark-loop line indicates the in- a disconnected sea-quark loop, is depicted in Fig. 1. These sertion of the matrix-element current. For an incoming baryon, B, disconnected-loop diagrams are difficult to calculate in nu- each of the three valence entering on the left hand side may contribute to the internal meson state, Φ. In the graded symmetry merical simulations of QCD due to the space-time coordinate approach, the conventional u, d and s sea-quark loops are comple- dependence of both the source and sink of the quark propa- mented by their commuting ghost counterparts, u˜, d˜and s˜. gator at the current insertion point. As the loop is correlated with other quark degrees of freedom only via exchange, resolving a nontrivial signal requires high statistics through all cases, the disconnected loops of Fig. 1 will make important innovative methods. contributions to the matrix elements of these states. While there has been recent success in isolating Recent lattice QCD calculations of the electromagnetic disconnected-loop contributions in ground-state baryon ma- form factors of the Λ(1405) have reported evidence of an im- trix elements, even at non-zero momenta [1, 2], the challeng- portant KN molecular component in the Λ(1405) [13]. On ing nature of these calculations suggests that resolving similar the lattice, the contribution to the magnetic form contributions for hadronic excitations in lattice QCD [3–16] factor of the Λ(1405) approaches zero, signaling the forma- will be elusive. tion of a spin-zero in a relative S-wave about the nu- However, it is here that the meson-baryon dressings of cleon. However, the light u and d valence-quark sector of the excitations are particularly important. Consider, for Λ(1405) makes a non-trivial contribution, which can be re- arXiv:1509.08226v2 [hep-lat] 10 Nov 2016 example, the Λ(1405) resonance which is of particular con- lated to the and electromagnetic form factors. temporary interest. It is widely agreed that there is a two-pole To understand the physics behind this contribution and en- structure for this resonance [17–28], associated with attrac- able a comparison with experiment, one needs to understand tive interactions in the πΣ and KN channels. The two poles the role of disconnected sea-quark loop contributions to the lie near the πΣ and KN thresholds and both channels can matrix elements of the Λ(1405). These contributions have contribute to the flavour singlet components of the resonance, not been included in present-day lattice QCD calculations and studied in detail herein. This two-pole structure, as first an- present a formidable challenge to the lattice QCD community. alyzed in Ref. [17], has now been entered into the When disconnected sea-quark loop-contributions are un- Data Group tables [29]. This achievement reflects decades of known, one often seeks observables insensitive to this con- research using the successful chiral unitary approaches [30– tribution. For example, the isovector baryon combination is 42], which have elucidated the two-pole structure [17]. In insensitive to a sea-quark contribution as it is similar for both describing the Λ(1405) resonance at quark masses larger than isospin variants, and cancels. For baryon states that are not in the physical masses, the introduction of a bare state dressed an isovector formation, a reliable estimation of the sea-quark by meson-baryon interactions has also been explored [28]. In loop contributions remains a pressing issue in the accurate 2 comparison of lattice results with experiment. group SU(3|3)V , as in Ref. [43]. The anticommuting be- havior of the standard quark degrees of freedom enable our symmetric baryon field in the graded symmetry Lagrangian to B. Flavor-singlet symmetry represent the flavor-singlet baryon. Our aim here is to disentangle quark-flow connected and In lattice QCD calculations, the Λ(1405) is excited from the disconnected couplings between singlet baryons and octet vacuum using interpolating fields dominated by a flavor sin- mesons and baryons in modern 2 + 1-flavor lattice QCD sim- glet operator. Thus, a complete understanding of the Λ(1405) ulations. As such, there are no issues with the artifacts of 0 necessarily involves understanding the role of disconnected quenched QCD, such as a light double-pole η meson. More- sea-quark loop contributions in the flavor-singlet sector of ef- over, we will consider quark masses similar to those of nature, fective field theory. Away from the SU(3)-symmetric limit, with mu ' md ≤ ms. flavor-singlet and flavor-octet interpolators mix in isolating the Λ(1405) [4] and therefore this new information in the flavor-singlet sector will complement the flavor-octet sector. D. Outline Flavor symmetry is important in effective field theory as it provides guidance on the meson-baryon couplings to be The structure of the paper is as follows. Section II expands used in coupled-channel analyses. It also has the ability to upon the graded symmetry approach to include the couplings describe the relative strength between quark-flow connected of flavor-singlet baryons to octet mesons and baryons, and the and disconnected contributions to observables. This separa- Lagrangian is established. In Section III, particle labels for tion is well understood in the flavor-octet sector of effective baryons containing a single ghost quark are explicitly written. field theory, but is unknown in the flavor-singlet sector. This This partly follows the method described in Ref. [46], with a gap in understanding is addressed in this article. few differences in the notation, which are introduced for con- It is perhaps important to note that the physical Λ(1405) venience in identifying the quantum numbers and symmetry is a mix of flavor symmetries due to the explicit breaking of properties of the baryons considered. As a relevant test case, SU(3) symmetry by the quark masses. Our task is not to de- the flavor-singlet component of the Λ(1405) is described in termine this mixing, but rather to disclose the separation of Sec. IV. connected and disconnected contributions in the flavor-singlet sector for the first time. Thus, in referring to flavor-singlet baryons in the following, it is to be understood that we are II. GRADED SYMMETRY FORMALISM referring to the flavor-singlet component of a baryon. In this manner, our results are applicable to Λ-baryons in general. In this section, the details of the graded symmetry approach are briefly reviewed. In order to derive couplings to flavor- singlet baryons, the standard graded symmetry approach re- C. Graded symmetry approach quires a complementary augmentation in order to include in- teractions between octet and singlet baryons, octet mesons, The separation of connected and disconnected contribu- and their ghost counterparts. tions to the baryon couplings can be handled using the graded For 2+1 flavor dynamical- lattice QCD calcula- symmetry approach [43–49] or the diagrammatic approach tions, the relevant graded group is SU(6|3)V , composed [50, 51]. Although it has been demonstrated that both meth- of three valence-quark flavors, three sea-quark flavors and ods are consistent with each other [50, 51], in this work, the three ghost-quark flavors [46]. While the sea quark masses graded symmetry approach is used because of its ability to be are free to differ from the valence sector, the ghost sec- easily generalized to include flavor-singlet baryons. The dia- tor masses are matched to the valence sector such that grammatic approach becomes cumbersome when the valence their fermion-determinant contributions exactly cancel the content of the baryon comprises three unique quark flavors fermion-determinant contributions of the valence sector, thus (u d s) [50]. ensuring the valence sector is truly valence only. Once the coupling strength of disconnected sea-quark loop However, a common approach used in practice, when sep- contributions to hadron two-point functions is understood, one arating connected and disconnected contributions, is to work can use this information to estimate the size of sea-quark loop within the simpler graded group SU(3|3)V of quenched QCD. contributions to matrix elements [51–57]. Here the standard quark sector is complemented by a ghost The formalism presented herein is based on the original quark sector with matching quark masses to render the quark work of Labrenz and Sharpe [43], and builds upon the con- sector to a valence quark sector. Identification of the ghost structions of that work in a consistent manner. We create a sector contributions to quark-flow diagrams is then sufficient new, completely symmetric baryon field in the graded symme- to identify the sea-quark loop contribution of the standard try Lagrangian, and couple this field to octet mesons and octet quark sector that is being removed. This approach is used baryons. Through the admission of the commuting ‘ghost’ herein. counterparts, introduced to isolate the disconnected loop com- We note, however, that there are important differences be- ponent of each coupling, the collectively fall into rep- tween the two graded groups, particularly for the flavor-singlet 0 0 resentations of the vector subgroup of the graded symmetry η meson. In SU(3|3)V of quenched QCD, the η appears 3

0 only as a double-hairpin diagram, whereas in SU(6|3)V the η η˜d and η˜s mesons along the diagonal of Eq. (4). This is use- appears in a partially-quenched generalization of the standard ful when considering the ghost sector of a flavor-octet baryon η0. In our case of interest, it is sufficient to hold the masses transition to a flavor-singlet baryon and an octet meson. of each each flavor in the valence, sea and ghost sectors equal, The group transformation properties of the meson field can such that all mesons and baryons participating in the diagrams be made compatible with the baryon fields via the methods of take their standard properties. χPT, where the exponential of the mesons is incorporated into In deriving a Lagrangian that couples flavor-singlet baryons an axial-vector field to octet baryons and mesons under the graded symmetry sub- ξ(x) = exp(iΦ(x)/f ), (5) group SU(3|3)V , the constructions introduced in Ref. [43] π will be used as much as possible, so that the extension to sin- i A = (ξ∂ ξ† − ξ†∂ ξ). (6) glet couplings can draw on the previously established formal- µ 2 µ µ ism. µ In the graded symmetry approach, the quark fields are ex- The elements of A transform under the graded-symmetry tended to include bosonic ghost counterparts for each flavor. subgroup SU(3|3)V as T The new field, Q = (u, d, s, u,˜ d,˜ s˜) , transforms under the µ µ † Aij → Uii0 Ai0j0 Uj0j. (7) vector subgroup of the graded symmetry group SU(3|3)V as a fundamental representation, analogous to three-flavor chiral perturbation theory (χPT) B. Octet and singlet baryons Q → U Q , Q → Q U † , (1) i ij j i j ji In the graded symmetry approach, the baryon field is gen- where the indices i, j run from 1 to 6. eralized to a rank-3 tensor, allowing the three constituent fla- vors of each baryon to be specified separately by each flavor index. This is convenient, as it readily allows one to identify A. Mesons the valence quark content of a baryon simply by inspection. The introduction of graded symmetry into the transformation matrices, U, in Eq. (1) requires the introduction of a grading The meson field is expanded to a 6 × 6 matrix, with factor, defined as four blocks of mesons corresponding to the ordinary mesons, ( π (qq), the entirely-ghost mesons, π˜ (˜qq˜), and the fermionic 1 i = 1 3 † if – composite mesons χ (˜qq) and χ (qq˜) in the off-diagonal ηi = , (8) 0 if i = 4 – 6 blocks

π χ† which takes into account the bosonic nature of the ghost Φ = . (2) χ π˜ quarks – swapping two quark fields incurs a minus sign only if both are anticommuting (Dirac-spinor field) quarks. The elements of the ordinary meson block follow the standard The transformation properties of the octet baryon field can convention be derived from the standard interpolating field

 1 0 1 + +  γ α,a β,b γ,c α,a γ,c β,b √ π + √ η π K Bijk ∼ [Qi Qj Qk − Qi Qj Qk ]abc(C γ5)αβ , (9) 1 2 6 η0  π− − √1 π0 + √1 η K0  π = √ 2 6 + √ I3. where a – c are color indices, α – γ are Dirac indices, and C is 2  0  6 K− K − √2 η 6 the standard charge conjugation matrix. This construction en- (3) sures the antisymmetry of the latter two Dirac indices (β, γ), With the inclusion of ghost quarks, the flavor of the ghost and thus eliminates the flavor-singlet baryon [43]. quark or antiquark may specified through either the usual The pure flavor-singlet baryon field, denoted BS, on the quark mixing, to form π˜0, η˜ and η˜0 particles, or through neu- other hand, is totally symmetric in quark flavor, Q. This can tral meson labels η˜u, η˜d and η˜s. Both representations have be encoded into the form of the baryon tensor by permuting their merit. When the full flavor symmetry of the meson dress- the color and Dirac indices while maintaining the flavor-label ing is present, the standard representation with π˜0, η˜ and η˜0 is ordering (ijk) as in Eq. (9), such that the established baryon of particular utility as the mass of the intermediate meson is transformation rules are maintained † manifest. Therefore, we consider the following form for χ S γ h α,a β,b γ,c β,b γ,c α,a γ,c α,a β,b Bijk ∼ Qi Qj Qk + Qi Qj Qk + Qi Qj Qk  √1 0 √1 + ˜ +  π˜ + η˜ π˜ K β,b α,a γ,c γ,c β,b α,a α,a γ,c β,bi 2 6 0 + Q Q Q + Q Q Q + Q Q Q † 1  π˜− − √1 π˜0 + √1 η˜ K˜ 0  η˜ i j k i j k i j k χ = √  2 6 +√ I3 . 2  0  6 × abc (C γ5)αβ. (10) K˜ − K˜ − √2 η˜ 6 (4) By transforming the individual quark operators (Q) within On the other hand, the flavor symmetry can be incomplete. In the baryon fields of Eqs. (9) and (10), and collecting the trans- this case, it is more transparent to work with the neutral η˜u, formation matrices (U) together using the commutation rule 4

[47] Similarly, the baryon tensor containing the flavor-singlet baryon follows from Eq. (10) for BS ηi(ηj +ηk) QiUjk = (−1) Ujk Qi , (11) S 1 Bijk = √ + ijk0 Bkk0 + jkk0 Bik0 + kik0 Bjk0 the resultant transformation rules for the baryon and an- 3 2  tibaryons (suppressing the Dirac indices) are − jik0 Bkk0 − kjk0 Bik0 − ikk0 Bjk0 . (17) √ η 0 (ηj +η 0 )+(η 0 +η 0 )(ηk+η 0 ) −1 Bijk → (−1) i j i j k Uii0 Ujj0 Ukk0 Bi0j0k0 , Again, the normalization factor of (3 2) ensures the (12) squares of the coefficients sum to 1. For the Λ0, one con- siders the six terms, BS , BS , BS , BS , BS and BS , ηi0 (ηj +ηj0 )+(ηi0 +ηj0 )(ηk+ηk0 ) † † † 123 132 213 231 312 √321 B → (−1) B 0 0 0 U U 0 U 0 , kji k j i k0k j j i i each of which provides a contribution with magnitude Λ0/ 6. (13) Similarly √ S S 0 for both B and B . ijk Bijk = 6 Λ , (i – k = 1 – 3) . (18) In matching the elements of the generalized baryon field to those of the standard baryon matrix Bkk0 , the indices are The symmetry relations that describe the relationships restricted from 1 to 3, and the elements are assigned so that among the elements of the baryon tensors are important when the result is antisymmetric under i ↔ j considering their general form in the graded symmetry group, i.e. where the indices i – k run from 1 to 6. For the octet

Bijk ∝ ijk0 Bkk0 , (i – k = 1 – 3), (14) field B, the following two symmetry relations are sufficient to specify the two group invariants of the symmetry, α and β where the convention used for the standard baryon field of [43], which are related to the standard D and F couplings χPT is ηj ηk+1 Bikj = (−1) Bijk , (19a)

 1 0 1 +  ηiηj ηiηj +ηj ηk+ηkηi √ Σ + √ ΛΣ p Bijk = (−1) Bjik + (−1) Bkji . (19b) 2 6 Λ0 B =  Σ− − √1 Σ0 + √1 Λ n  + √ .  2 6  I3 − 0 2 2 3 Similarly, the following symmetry relations exist for the field Ξ Ξ − √ Λ S 6 B (15) 0 S ηiηj S Λ denotes a positive- flavor-singlet baryon. Similarly, Bjik =(−1) Bijk , (20a) a matrix B∗ can be defined with the same form, representing BS =(−1)ηj ηk BS , (20b) the lowest-lying negative-parity baryons. The negative-parity ikj ijk 0∗ S ηiηk+ηiηj S flavor-singlet baryon is then denoted by Λ . Bjki =(−1) Bijk , (20c) The baryon construction of Eq. (14) contains both octet and BS =(−1)ηiηk+ηj ηk BS , (20d) singlet fields. However, it is known that interactions that in- kij ijk S ηiηj +ηiηk+ηj ηk S clude a singlet baryon have couplings that are not related via Bkji =(−1) Bijk . (20e) flavor symmetry to the interaction couplings between octet baryons only. It is therefore important to separate the baryon field into ‘octet-only’ and ‘singlet-only’ pieces, so that the C. The Lagrangian terms in the Lagrangian can be assigned independent coef- ficients. The singlet Lagrangian will couple the singlet-baryon field To extract the purely octet component of Bijk, Labrenz to a combination of octet-meson and octet-baryon fields, A · and Sharpe antisymmetrized over the j and k indices of Q B. Given the maximal symmetry for BS in Eq. (10), there in Eq. (9), so that the symmetric component (containing the are three different ways of contracting the octet-baryon flavor singlet) is removed. In terms of the standard baryon field, indices. In keeping with the baryon construction of Labrenz with indices i, j, and k restricted to 1 – 3 and fermion anti- and Sharpe [43], we consider a Lagrangian term similar to S commutation taken into account, this takes the following form their term with coefficient α, Bkji Akk0 Bijk0 , where the first two indices of B are contracted with BS; a term similar to S 1 their term with coefficient β, Bkji Aik0 Bk0jk, where the latter Bijk = √ (ijk0 Bkk0 + ikk0 Bjk0 ) . (16) two indices of B are contracted with BS; and a third candidate 6 S √ Bkji Akk0 Bik0j, where the outer indices of B are contracted The factor of 1/ 6 is chosen so that the baryon wave func- with BS. Noting that the octet field of Eq. (9) has a symmetry tions have the correct normalization. This becomes apparent constraint under the exchange of the latter two indices as in when assigning conventional baryon labels to the elements of Eq. (19a), this third candidate is related to the first candidate the B tensor (discussed in Sec III). For example, extracting and may be discarded. Similarly, the relation of Eq. (20a) S the proton states involves selecting√ √ the elements√B112, B121 for B is opposite to that of Eq. (19a) for B and therefore and B211, with weights of 1/ 6, 1/ 6 and −2/ 6, respec- the second candidate vanishes. Thus, there is only one group tively, whose squares sum to 1. invariant combination. 5

In the general case where all indices run from 1 – 6, the labels that provide information about the charge, quark com- commutation/anticommutation behavior of the elements of position and flavor symmetry. Thus, in defining the baryon la- the transformation matrices U depends on their location in bels, there are some differences in our notation with respect to the block matrix. As a result, grading factors must be in- Ref. [46]. Here, the choice in baryon labels is to manifest the cluded when swapping the order of the indices. Noting that symmetry and flavor properties of baryons with a single ghost S Eq. (20e) indicates that Bijk is antisymmetric under exchange quark. If one of the quarks is a (commuting) ghost particle, of the outer indices in both the normal quark and single ghost- it may be combined with the two remaining ordinary quarks quark sectors, the interaction Lagrangians take the following in a Σ-like symmetric or a Λ-like antisymmetric configuration forms of ordinary quarks. For example, for a proton with a ghost up ˜ + ˜ + √ S quark, one defines new particles Σp,u˜, and Λp,u˜, which cor- (ηi+ηj )(ηk+ηk0 ) LS = 2 2 gs (−1) Bkji Akk0 Bijk0 , (21) respond to symmetric and antisymmetric combinations of the √ S two remaining quark flavors, respectively. (ηi+ηj )(ηk+ηk0 ) S L = 6 g (−1) B A 0 B 0 . (22) SS ss kji kk ijk In addition, unphysical states such as doubly-charged Σ˜ ++ (˜u, u, u) We note these are special cases of the general irreducible rep- ‘’ p,u˜ and negatively charged ‘’ Σ˜ − (d,˜ d, d) are members of the single ghost-quark baryons. resentations of graded Lie groups classified by Balantekin and n,d˜ Bars [58, 59]. The leading coefficient in LS is selected to de- These states occur in both the graded-symmetry approach and 0 fine the coupling associated with the process Λ → π0Σ0 as the diagrammatic approach by considering all possible quark- gs. flow topologies. Note that these states do not contribute in Our investigation of the separation of quark-flow connected full QCD, where the cancellation of the connected and discon- and disconnected contributions to standard baryons constrains nected quark-flow diagrams associated with them enforces the the indices i – k to 1 – 3, whereas k0 = 1 – 6 incorporates the Pauli exclusion principle. ghost meson-baryon contributions, which reveal the discon- Following Ref. [46], the flavor-symmetric variants are en- nected contributions. coded in a rank-3 tensor sijk, which is symmetric in the last In summary, the first-order interaction Lagrangian for low- two indices reserved for ordinary quarks lying octet and singlet mesons, octet baryons and low-lying s = Σ˜ ++, s = Σ˜ + , s = Σ˜ + , negative-parity singlet baryons, including the Dirac structure, 411 p,u˜ 511 p,d˜ 611 Σ,s˜ √ √ √ takes the form ˜ + ˜ 0 ˜ 0 s412 = Σp,u˜/ 2, s512 = Σn,d˜/ 2, s612 = ΣΣ,s˜/ 2, (1) 0 − − L = Lα + Lβ + (LS + h.c.) + LSS , s = Σ˜ , s = Σ˜ , s = Σ˜ , BA 422 n,u˜ 522 n,d˜ 622 Σ,s˜ (ηi+ηj )(ηk+ηk0 ) µ √ √ √ = 2α (−1) Bkji γµγ5 Akk0 Bijk0 ˜ + ˜ 0 ˜ 0 s413 = ΣΣ,u˜/ 2, s513 = ΣΣ,d˜/ 2, s613 = ΣΞ,s˜/ 2, µ √ √ √ + 2β Bkji γµγ5 Aik0 Bk0jk ˜ 0 ˜ − ˜ − √ s423 = ΣΣ,u˜/ 2, s523 = Σ ˜/ 2, s623 = ΣΞ,s˜/ 2, (ηi+ηj )(ηk+η 0 ) Σ,d + 2 2 gs (−1) k 0 − − s433 = Σ˜ , s533 = Σ˜ , s633 = Σ˜ . n S∗ µ µ S∗ o Ξ,u˜ Ξ,d˜ Ω,s˜ × B γ A B 0 + B 0 γ A B kji µ kk0 ijk k ji µ k0k ijk (24) √ S∗ (ηi+ηj )(ηk+ηk0 ) µ S∗ + 6 gss (−1) Bkji γµγ5 Akk0 Bijk0 , The baryon subscripts reflect the combined number of strange (23) quarks and strange ghost quarks and the superscripts record the total electric charge. The flavor of the residing ghost quark S∗ where Bijk denotes the odd-parity generalized field associated is explicit in the subscript. We note that although Ω has been ∗ 0∗ with B and thus Λ via Eq. (17). The loss of γ5 in the terms used to denote the (˜s, s, s) combination, this baryon is actually coupling octet- and singlet-flavor baryons reflects the odd par- associated with octet as opposed to decuplet baryons. ity of the low-lying Λ0∗. Using a similar notation, the Λ-like particles are encoded in an antisymmetric tensor, denoted tij t = Λ˜ 0 , t = Λ˜ − , t = Λ˜ − , III. GHOST BARYON STATE IDENTIFICATION 41 Σ,u˜ 51 Σ,d˜ 61 Ξ,s˜ ˜ + ˜ 0 ˜ 0 t42 = Λ , t52 = Λ ˜, t62 = ΛΞ,s˜, In this section we review briefly the process of relating the Σ,u˜ Σ,d S ˜ + ˜ 0 ˜ 0 elements of B and B to baryons having one ghost quark and t43 = Λp,u˜, t53 = Λn,d˜, t63 = ΛΣ,s˜. (25) two ordinary quarks. This is relevant to our investigation of the ordinary Λ0∗ baryon, where one ghost quark can be intro- These baryons are distributed over the elements of B in the duced through the disconnected sea-quark loop. The ordinary following manner. The manifestly symmetric part of B has baryons (those without ghost quark content) adopt the label- the simple relation ing convention assigned in Eqs. (15) through (17). The case r2 for one ghost quark in B was first present by Chen and Savage Bijk = sijk, (i = 4 – 6, j – k = 1 – 3), (26) [46]. We will extend the formalism to include BS. 3 In identifying the intermediate states of the meson-baryon which places the ghost quark in the first index position of B. channels dressing a baryon, we find it is useful to assign state The symmetric states with ghost quarks mapped to the second 6 and third index positions of B are combined with the antisym- A. Flavor-octet baryon dressings of the singlet baryon metric baryons of Eq. (25) in Here we consider the specific case relevant to the low- 1 1 Bjik = tip pjk + √ sijk , (27a) lying odd-parity Λ(1405) baryon, known to have an important 2 6 flavor-singlet component [4]. We commence by considering 0∗ 0∗ Bjki = −Bjik , (27b) the process Λ → B Φ → Λ via the Lagrangian of Eq. (23). In this case, the u, d and s valence quarks that comprise where i = 4 – 6, j – k = 1 – 3 and p = 1 – 3. These two the Λ0∗ state may each contribute to the Φ meson loop, shown equations map the symmetric√ baryon states to two elements in Fig. 1. The ghost quark sector may enter through the dis- of B with weight 1/ 6 such that the weighting of p2/3 in connected loop, and therefore there are three choices for the Eq. (26) provides a sum of squares equal to one. Similarly, antiquark plus three choices for the ghost antiquark. Noting Eqs. (27) distribute the antisymmetric ghost baryon states of that the flavor-singlet η0 meson cannot be combined with an Eqs. (25) over four elements of B with weight 1/2 such that a octet baryon to form a flavor singlet, there are (3×6)−1 = 17 0∗ sum of squares again equals one. terms in the Lagrangian for the Λ Φ B coupling. At one loop, The complete flavor antisymmetry of BS for ordinary the baryon B contains at most one ghost quark for which the quarks restricts its first ghost sector to baryon states having an- meson contains the corresponding antighost quark. tisymmetry in the ordinary quarks, i.e. the baryons of Eq. (25) The SU(3|3)V -flavor composition of the interaction La- grangian of Eq. (23) contains the regular full-QCD meson- baryon couplings to the singlet Λ0∗ S 1 0 Bijk = √ tip pjk , (28a) 6 QCD 0∗ n 0 − L 0∗ = gs Λ K n + K p BS = +BS , (28b) S,Λ jik ijk + − 0 0 − + S S + π Σ + π Σ + π Σ Bjki = +Bijk . (28c) o + K+ Ξ− + K0 Ξ0 + η Λ , (29) where i = 4 – 6 and j – k = 1 – 3. Eqs. (28), now with a prime to denote the association with the flavor-singlet baryon and the ghost meson to ghost baryon couplings sector, maps each baryon of Eqs. (25) over six elements of BS. S∗ Analogous results are obtained for the odd-parity B using  r 0 ijk Ghost 0∗ 2  ˜ − ˜ + ˜ ˜ 0 0∗ L 0∗ = −gs Λ K Λ + K Λ ˜ tip. S,Λ 3 p,u˜ n,d In studies of the Λ0∗, the low-lying ghost states reside in the +π ˜+ Λ˜ − +π ˜− Λ˜ + octet meson and octet baryon sectors. The odd-parity ghost Σ,d˜ Σ,u˜ states of BS∗ are of relevance to other states transiting to Λ0∗  ijk + K˜ + Λ˜ − + K˜ 0 Λ˜ 0 and a meson. As an example, consider an octet baryon tran- Ξ,s˜ Ξ,s˜ sition to a singlet baryon and an octet meson. Such physics 1  0 0  + √ π˜0 Λ˜ − π˜0 Λ˜ is relevant to understanding the disconnected loop contribu- 3 Σ,u˜ Σ,d˜ 0∗ 0∗ tions in Σ → π Λ or Ξ → K Λ . These become particu- 1   larly important in lattice QCD, where the initial state can be a + η˜Λ˜ 0 +η ˜Λ˜ 0 − 2η ˜Λ˜ 0 , 3 Σ,u˜ Σ,d˜ Σ,s˜ finite-volume excitation associated with a resonance. (30) This defines the method by which baryon labels are as- signed to different elements of the three-index baryon tensors, S where the Dirac structure has been suppressed. In the full- B and B , and their odd-parity analogues. By writing out the QCD case, the relative coupling strengths of the πΣ, KN, full Lagrangian in Eq. (23), one can simply pick out the rele- KΞ and ηΛ channels are consistent with early work [60] as vant contributions to each baryon in each interaction channel expected. SU(3|3) available for the symmetry subgroup V . The ghost terms, however, present novel features. Only an- tisymmetric Λ˜ ghost baryons contribute. Thus, the coupling of a baryon with a single ghost quark to a flavor-singlet baryon IV. FLAVOR-SINGLET Λ-BARYON COUPLINGS is constrained to have its two ordinary quarks in an antisym- metric formation. The Σ˜ baryons do not contribute. In this section we examine the coupling strengths for tran- In addition, flavor-singlet symmetry is manifest. Each me- sitions involving the flavor-singlet component, Λ0∗. As de- son charge state participating in the quark-flow disconnected scribed at the beginning of Sec. II, the couplings of the ghost loop contributes in a uniform manner with a coupling strength 2 2 sector in SU(3|3)V provide the couplings of quark-flow dis- of 3 gs . connected sea-quark loop contributions in SU(6|3)V that we The ghost-quark subscript on the baryon field enables the are seeking. In writing out the relevant terms of the La- identification of the flavor of the quark participating in the grangian, we report the meson-baryon → baryon transition loop. Thus, matrix elements can be separated into their quark- terms, noting that the baryon → meson-baryon terms are given flow connected and disconnected parts for each quark flavor. by the Hermitian conjugate. Similarly, the baryon subscripts facilitate the incorporation of 7

SU(3)-flavor breaking effects through variation of the hadron quark loop. After summing over the participating flavors, each masses. full-QCD process proceeds with one third of the strength in The essential feature to be drawn from this analysis is that quark-flow connected diagrams and two-thirds of the strength the relative contribution of quark-flow connected and dis- in quark-flow disconnected diagrams. connected diagrams is now known. Consider the process 0∗ − 0∗ 2 Λ → K p → Λ proportional to gs . This process pro- ceeds through a combination of quark-flow connected dia- B. Flavor-singlet baryon dressings of octet baryons grams and the quark-flow disconnected diagram of Fig. 1. See for example, Ref. [51]. The ghost Lagrangian reports a term It is also of interest to examine processes where the flavor- singlet baryon participates as an intermediate state dressing of r an octet baryon. The full Lagrangian describing the processes 2 0∗ −g Λ K˜ − Λ˜ + , (31) B → Φ BS∗ → B, where B is an octet baryon, can likewise be s 3 p,u˜ decomposed into full-QCD and ghost meson to ghost baryon components, suppressing the Dirac structure indicating that this process proceeds through the quark-flow 2 2 disconnected diagram with strength 3 gs . That is, two thirds QCD n 0 0∗ + 0∗ of the full-QCD process Λ0∗ → K− p → Λ0∗ proceeds by the LS,Λ0∗ = gs n (K Λ ) + p (K Λ ) quark-flow disconnected diagram with a u quark flavor partic- − 0 + + Σ (π− Λ0∗) + Σ (π0 Λ0∗) + Σ (π+ Λ0∗) ipating in the loop. The fact that two-thirds of the strength lies − 0 0 in the loop is perhaps surprising. + Ξ (K− Λ0∗) + Ξ (K Λ0∗) The full-QCD processes involving the neutral π0 or η o + Λ(η Λ0∗) , mesons are only slightly more complicated, as one needs to (32) sum over more than one quark flavor participating in the sea-

     Ghost − ˜ +0∗ ˜ 00∗ ˜ 0 ˜ 00∗ ˜ +0∗ + ˜ 00∗ ˜ + ˜ 00∗ LS,Λ0∗ = gs n π˜ Λp,u˜ +η ˜d Λn,d˜ + K ΛΣ,s˜ + p η˜u Λp,u˜ +π ˜ Λn,d˜ + K ΛΣ,s˜

−   +   + Σ π˜− Λ˜ 00∗ +η ˜ Λ˜ −0∗ + K˜ 0 Λ˜ −0∗ + Σ η˜ Λ˜ +0∗ +π ˜+ Λ˜ 00∗ + K˜ + Λ˜ 00∗ Σ,u˜ d Σ,d˜ Ξ,s˜ u Σ,u˜ Σ,d˜ Ξ,s˜   1 0 − +0∗ 00∗ 00∗ + −0∗ 0 00∗ + −0∗ + √ Σ − π˜ Λ˜ +η ˜u Λ˜ − η˜d Λ˜ +π ˜ Λ˜ − K˜ Λ˜ + K˜ Λ˜ 2 Σ,u˜ Σ,u˜ Σ,d˜ Σ,d˜ Ξ,s˜ Ξ,s˜ −  0  0  0  + Ξ K˜ − Λ˜ 00∗ + K˜ Λ˜ −0∗ +η ˜ Λ˜ −0∗ + Ξ K˜ − Λ˜ +0∗ + K˜ Λ˜ 00∗ +η ˜ Λ˜ 00∗ Σ,u˜ Σ,d˜ s Ξ,s˜ Σ,u˜ Σ,d˜ s Ξ,s˜

 0   1 − +0∗ 00∗ 00∗ 00∗ + −0∗ − +0∗ 00∗ 0 00∗ + −0∗ + √ Λ π˜ Λ˜ +η ˜u Λ˜ +η ˜d Λ˜ − 2˜ηs Λ˜ +π ˜ Λ˜ − 2 K˜ Λ˜ − 2 K˜ Λ˜ + K˜ Λ˜ + K˜ Λ˜ . 6 Σ,u˜ Σ,u˜ Σ,d˜ Σ,s˜ Σ,d˜ p,u˜ n,d˜ Ξ,s˜ Ξ,s˜ (33)

In each channel of Eq. (33), there are terms which represent there is no ss contribution in any channel apart from the octet the disconnected loop contributions of the normal full-QCD Λ. In the Σ, where both light flavor ηq mesons appear, the interactions of Eq. (32), and there are additional terms that flavor symmetry of the is apparent. Similarly, in the Λ serve to cancel the unphysical quark-flow connected diagrams channel, the flavor symmetry of the ηq contributions identi- associated with each channel. For example, in the case of the fies the octet-η participating to cancel the full-QCD process ˜ 0 ˜ 00∗ 0∗ neutron channel n, the relevant ghost term K ΛΣ,s˜ comes Λ(η Λ ). In all other cases, the light-quark ηu and ηd mesons with the same magnitude as the normal process, gs, indicat- propagate with the mass of the pion. ing that there is no completely-connected contribution for this process. The physical process proceeds through a quark-flow disconnected diagram. The other two terms in this channel, π˜− Λ˜ +0∗ and η˜ Λ˜ 00∗ occur to cancel the aforementioned un- p,u˜ d n,d˜ C. Flavor-singlet baryon dressings of the singlet baryon physical connected diagrams. These contributions are also im- portant and must be included when evaluating the quark-flow Finally, we also consider the quark-flow disconnected loop disconnected contributions to baryon matrix elements. contributions to the process Λ0∗ → Φ BS∗ → Λ0∗. We find Note that, unlike the flavor-singlet Λ-baryon couplings of 0 QCD 0∗ 0 0∗ Eq. (30), where the structure of the π˜ and η˜ is manifest, here LSS,Λ0∗ = gss Λ η Λ , (34) 8 and is known to play an important role [4]. There only the quark- flow connected contributions are calculated.  r 00 Because a determination of the quark-flow disconnected Ghost gss 0∗ 2  ˜ −0 ˜ +0∗ ˜ ˜ 00∗ L 0∗ = √ Λ K Λ + K Λ contributions to baryon excited states presents formidable SS,Λ 3 p,u˜ n,d˜ 2 challenges to the lattice QCD community, it is essential to +π ˜+0 Λ˜ −0∗ +π ˜−0 Λ˜ +0∗ Σ,d˜ Σ,u˜ have an understanding of the physics missing in these contem- porary lattice QCD simulations. An understanding of the cou- ˜ +0 ˜ −0∗ ˜ 00 ˜ 00∗  + K ΛΞ,s˜ + K ΛΞ,s˜ plings of the quark-flow disconnected contributions to baryon matrix elements is vital to understanding QCD and ultimately 1  00 ˜ 00∗ 00 ˜ 00∗  + √ π˜ Λ − π˜ Λ ˜ comparing with experiment. 3 Σ,u˜ Σ,d  We have discovered that the full-QCD processes 1  0 ˜ 00∗ 0 ˜ 00∗ 0 ˜ 00∗  + η˜ ΛΣ,u˜ +η ˜ Λ ˜ − 2η ˜ ΛΣ,s˜ . 0 3 Σ,d Λ0∗ → K− p → Λ0∗ , Λ0∗ → K n → Λ0∗ , (35) Λ0∗ → π+ Σ− → Λ0∗ , Λ0∗ → π0 Σ0 → Λ0∗ , Λ0∗ → π− Σ+ → Λ0∗ , Λ0∗ → K+ Ξ− → Λ0∗ , Here we have added a prime to the meson labels to remind the reader that in full QCD these contributions enter with the Λ0∗ → K0 Ξ0 → Λ0∗ and Λ0∗ → η Λ → Λ0∗ , 0 meson propagating with the mass of the flavor-singlet η . The proceed with two thirds of the full-QCD strength in the quark- symbols serve to indicate the manner in which specific quark- flow disconnected diagram. Details of the quark flavors par- flow disconnected loops contribute to the full-QCD process. ticipating in the loops and their relative contributions are pro- Again, we observe the physical process proceeding through vided in Eq. (30) of Sec. IV A. a quark-flow disconnected diagram. Here, the final term, indi- We have also considered the separation of quark-flow con- cating the ghost-η0 contributions, provides a total contribution nected and disconnected contributions for the flavor-singlet with strength gss, precisely removing the full-QCD contribu- baryon dressings of octet baryons in Sec. IV B, and flavor- tion of Eq. (34). The remaining ghost terms enter to ensure singlet baryon dressings of the the flavor-singlet baryon in unphysical contributions from quark-flow connected diagrams Sec. IV C. do not contribute in full QCD. Their contributions are impor- When combined with previous results for flavor-octet tant and must be included when evaluating the quark-flow dis- baryons, this information enables a complete understanding connected contributions to baryon matrix elements. of the contributions of disconnected sea-quark loops in the matrix elements of baryons and their excitations, including the Λ(1405). This allows an estimation of quark-flow dis- connected corrections to lattice QCD matrix elements, or an V. CONCLUSION adjustment of other observables prior to making a comparison with the quark-flow connected results of lattice QCD.

An extension of the graded symmetry approach to partially quenched QCD has been presented. The extension builds on previous theory by explicitly including flavor-singlet baryons ACKNOWLEDGMENTS in its construction. The aim is to determine the strength of both the quark-flow connected and quark-flow disconnected We thank Tony Thomas for inspiring this study, Phiala diagrams encountered in full-QCD calculations of hadronic Shanahan for helpful discussions and Stephen Sharpe for in- matrix elements containing flavor-singlet baryon components. sightful comments in the preparation of this manuscript. This This is particularly important in light of recent progress research is supported by the Australian Research Council in the calculation of the electromagnetic properties of the through Grants No. DP120104627, No. DP140103067 and Λ(1405) in lattice QCD [13], where flavor-singlet symmetry No. DP150103164 (D.B.L.).

[1] J. Green, S. Meinel, M. Engelhardt, S. Krieg, J. Laeuchli, [4] B. J. Menadue, W. Kamleh, D. B. Leinweber, and M. S. Mah- J. Negele, K. Orginos, A. Pochinsky, and S. Syritsyn, Phys. bub, Phys. Rev. Lett. 108, 112001 (2012), arXiv:1109.6716 Rev. D92, 031501 (2015), arXiv:1505.01803 [hep-lat]. [hep-lat]. [2] R. S. Sufian, Y.-B. Yang, A. Alexandru, T. Draper, K.-F. Liu, [5] R. G. Edwards, J. J. Dudek, D. G. Richards, and S. J. Wallace, and J. Liang, (2016), arXiv:1606.07075 [hep-ph]. Phys. Rev. D84, 074508 (2011), arXiv:1104.5152 [hep-ph]. [3] M. S. Mahbub, W. Kamleh, D. B. Leinweber, P. J. Moran, and [6] M. S. Mahbub, W. Kamleh, D. B. Leinweber, P. J. Moran, A. G. Williams (CSSM Lattice), Phys. Lett. B707, 389 (2012), and A. G. Williams, Phys. Rev. D87, 011501 (2013), arXiv:1011.5724 [hep-lat]. arXiv:1209.0240 [hep-lat]. 9

[7] C. Lang and V. Verduci, Phys. Rev. D87, 054502 (2013), th/0109006 [nucl-th]. arXiv:1212.5055. [34] T. Hyodo, A. Hosaka, E. Oset, A. Ramos, and M. J. V. Vacas, [8] M. S. Mahbub, W. Kamleh, D. B. Leinweber, P. J. Moran, Phys. Rev. C68, 065203 (2003), arXiv:nucl-th/0307005 [nucl- and A. G. Williams, Phys. Rev. D87, 094506 (2013), th]. arXiv:1302.2987 [hep-lat]. [35] B. Borasoy, U. G. Meissner, and R. Nissler, Phys. Rev. C74, [9] C. Morningstar, J. Bulava, B. Fahy, J. Foley, Y. Jhang, et al., 055201 (2006), arXiv:hep-ph/0606108 [hep-ph]. Phys. Rev. D88, 014511 (2013), arXiv:1303.6816 [hep-lat]. [36] Y. Ikeda, T. Hyodo, and W. Weise, Phys. Lett. B706, 63 (2011), [10] C. Alexandrou, T. Leontiou, C. N. Papanicolas, and E. Stiliaris, arXiv:1109.3005 [nucl-th]. Phys. Rev. D91, 014506 (2015), arXiv:1411.6765 [hep-lat]. [37] T. Hyodo, D. Jido, and A. Hosaka, Phys. Rev. C85, 015201 [11] B. J. Owen, W. Kamleh, D. B. Leinweber, M. S. Mah- (2012), arXiv:1108.5524 [nucl-th]. bub, and B. J. Menadue, Phys. Rev. D92, 034513 (2015), [38] M. Dring and U. G. Meiner, Phys. Lett. B704, 663 (2011), arXiv:1505.02876 [hep-lat]. arXiv:1108.5912 [nucl-th]. [12] J. M. M. Hall, A. C. P. Hsu, D. B. Leinweber, A. W. [39] S. X. Nakamura and D. Jido, PTEP 2014, 023D01 (2014), Thomas, and R. D. Young, Phys. Rev. D87, 094510 (2013), arXiv:1310.5768 [nucl-th]. arXiv:1303.4157 [hep-lat]. [40] T. Sekihara and S. Kumano, Phys. Rev. C89, 025202 (2014), [13] J. M. M. Hall, W. Kamleh, D. B. Leinweber, B. J. Menadue, arXiv:1311.4637 [nucl-th]. B. J. Owen, A. W. Thomas, and R. D. Young, Phys. Rev. Lett. [41] M. Mai, V. Baru, E. Epelbaum, and A. Rusetsky, Phys. Rev. 114, 132002 (2015), arXiv:1411.3402 [hep-lat]. D91, 054016 (2015), arXiv:1411.4881 [nucl-th]. [14] Z.-W. Liu, W. Kamleh, D. B. Leinweber, F. M. Stokes, A. W. [42] M. Mai and U.-G. Meiner, Eur. Phys. J. A51, 30 (2015), Thomas, and J.-J. Wu, Phys. Rev. Lett. 116, 082004 (2016), arXiv:1411.7884 [hep-ph]. arXiv:1512.00140 [hep-lat]. [43] J. N. Labrenz and S. R. Sharpe, Phys. Rev. D54, 4595 (1996), [15] D. Leinweber, W. Kamleh, A. Kiratidis, Z.-W. Liu, S. Mahbub, arXiv:hep-lat/9605034 [hep-lat]. D. Roberts, F. Stokes, A. W. Thomas, and J. Wu, Proceed- [44] J. N. Labrenz and S. R. Sharpe, Nucl. Phys. Proc. Suppl. 34, ings, 10th International Workshop on the Physics of Excited 335 (1994), arXiv:hep-lat/9312067 [hep-lat]. (NSTAR 2015), JPS Conf. Proc. 10, 010011 (2016), [45] M. J. Savage, Nucl. Phys. A700, 359 (2002), arXiv:nucl- arXiv:1511.09146 [hep-lat]. th/0107038 [nucl-th]. [16] Z.-W. Liu, W. Kamleh, D. B. Leinweber, F. M. Stokes, A. W. [46] J.-W. Chen and M. J. Savage, Phys. Rev. D65, 094001 (2002), Thomas, and J.-J. Wu, (2016), arXiv:1607.04536 [nucl-th]. arXiv:hep-lat/0111050 [hep-lat]. [17] J. Oller and U.-G. Meißner, Phys. Lett. B 500, 263 (2001). [47] S. R. Beane and M. J. Savage, Nucl.Phys. A709, 319 (2002), [18] C. Garcia-Recio, J. Nieves, E. R. Arriola, and M. J. V. Vacas, arXiv:hep-lat/0203003 [hep-lat]. Phys. Rev. D67, 076009 (2003), arXiv:hep-ph/0210311 [hep- [48] W. Detmold, B. C. Tiburzi, and A. Walker-Loud, Phys. Rev. ph]. D73, 114505 (2006), arXiv:hep-lat/0603026 [hep-lat]. [19] D. Jido, A. Hosaka, J. C. Nacher, E. Oset, and A. Ramos, Phys. [49] P. E. Shanahan, A. W. Thomas, K. Tsushima, R. D. Rev. C66, 025203 (2002), arXiv:hep-ph/0203248 [hep-ph]. Young, and F. Myhrer, Phys. Rev.Lett. 110, 202001 (2013), [20] D. Jido, J. A. Oller, E. Oset, A. Ramos, and U. G. Meissner, arXiv:1302.6300 [nucl-th]. Nucl. Phys. A725, 181 (2003), arXiv:nucl-th/0303062 [nucl- [50] D. B. Leinweber, Phys. Rev. D69, 014005 (2004), arXiv:hep- th]. lat/0211017 [hep-lat]. [21] T. Hyodo, S.-i. Nam, D. Jido, and A. Hosaka, Prog. Theor. [51] J. M. M. Hall, D. B. Leinweber, and R. D. Young, Phys. Rev. Phys. 112, 73 (2004), arXiv:nucl-th/0305011 [nucl-th]. D89, 054511 (2014), arXiv:1312.5781 [hep-lat]. [22] V. K. Magas, E. Oset, and A. Ramos, Phys. Rev. Lett. 95, [52] R. D. Young, D. B. Leinweber, A. W. Thomas, and S. V. 052301 (2005), arXiv:hep-ph/0503043 [hep-ph]. Wright, Phys. Rev. D66, 094507 (2002), arXiv:hep-lat/0205017 [23] L. S. Geng, E. Oset, and M. Doring, Eur. Phys. J. A32, 201 [hep-lat]. (2007), arXiv:hep-ph/0702093 [hep-ph]. [53] D. B. Leinweber, A. W. Thomas, and R. D. Young, Phys. Rev. [24] Y. Ikeda, T. Hyodo, and W. Weise, Phys. Lett. B 706, 63 Lett. 92, 242002 (2004), arXiv:hep-lat/0302020 [hep-lat]. (2011). [54] D. B. Leinweber, S. Boinepalli, I. C. Cloet, A. W. Thomas, [25] Z.-H. Guo and J. A. Oller, Phys. Rev. C 87, 035202 (2013). A. G. Williams, R. D. Young, J. M. Zanotti, and J. B. Zhang, [26] M. Mai and U.-G. Meißner, Nucl. Phys. A 900, 51 (2013). Phys. Rev. Lett. 94, 212001 (2005), arXiv:hep-lat/0406002 [27] M. Doring, D. Jido, and E. Oset, Eur. Phys. J. A45, 319 (2010), [hep-lat]. arXiv:1002.3688 [nucl-th]. [55] D. B. Leinweber, S. Boinepalli, A. W. Thomas, P. Wang, A. G. [28] Z.-W. Liu, J. M. M. Hall, D. B. Leinweber, A. W. Thomas, and Williams, R. D. Young, J. M. Zanotti, and J. B. Zhang, Phys. J.-J. Wu, (2016), arXiv:1607.05856 [nucl-th]. Rev. Lett. 97, 022001 (2006), arXiv:hep-lat/0601025 [hep-lat]. [29] C. Patrignani, Chin. Phys. C40, 100001 (2016). [56] P. Wang, D. B. Leinweber, A. W. Thomas, and R. D. Young, [30] N. Kaiser, P. B. Siegel, and W. Weise, Nucl. Phys. A594, 325 Phys. Rev. C79, 065202 (2009), arXiv:0807.0944 [hep-ph]. (1995), arXiv:nucl-th/9505043 [nucl-th]. [57] P. Wang, D. B. Leinweber, and A. W. Thomas, Phys. Rev. D89, [31] N. Kaiser, T. Waas, and W. Weise, Nucl. Phys. A 612, 297 033008 (2014), arXiv:1312.3375 [hep-ph]. (1997). [58] A. Baha Balantekin and I. Bars, J. Math. Phys. 22, 1149 (1981). [32] E. Oset and A. Ramos, Nucl. Phys. A 635, 99 (1998). [59] A. Baha Balantekin and I. Bars, J. Math. Phys. 22, 1810 (1981). [33] E. Oset, A. Ramos, and C. Bennhold, Phys. Lett. B527, [60] E. Veit, B. K. Jennings, R. Barrett, and A. W. Thomas, Phys. 99 (2002), [Erratum: Phys. Lett.B530,260(2002)], arXiv:nucl- Lett. B137, 415 (1984).