A New Look at the $Y$ Tetraquarks and $\Omega_c$ Baryons in the Diquark Model

We analyze the hidden charm $P$-wave tetraquarks in the diquark model, using an effective Hamiltonian incorporating the dominant spin-spin, spin-orbit and tensor interactions. We compare with other $P$-wave system such as $P$-wave charmonia and the newly discovered $\Omega_c$ baryons, analysed recently in this framework. Given the uncertain experimental situation on the $Y$ states, we allow for different spectra and discuss the related parameters in the diquark model. In addition to the presently observed ones, we expect many more states in the supermultiplet of $L=1$ diquarkonia, whose $J^{PC}$ quantum numbers and masses are worked out, using the parameters from the currently preferred $Y$-states pattern. The existence of these new resonances would be a decisive footprint of the underlying diquark dynamics.


I. INTRODUCTION
The experimental discovery of exotic, hidden charm or beauty, states has opened a new field in hadron spectroscopy. The exotic states, called X, Y, Z and P c , have been analysed in a number of theoretical models. They have been claimed to be hybrid charmonia, molecules, disguised charmonium states, or just coupled-channel or threshold effects, see [1][2][3][4][5][6] for recent reviews and references therein. We concentrate on the alternative diquark-antidiquark interpretation, tetraquark in brief, introduced in [7,8] following the light pentaquark picture discussed in [9], which has the potential to include all exotic hadrons seen thus far in a single scheme.
The objects of our interest in this paper are the so called Y states, with J P C = 1 −− , described as tetraquarks with orbital angular momentum L = 1 and L = 3. Y -states have also been interpreted as hadron molecules in [10], with Y (4008) not foreseen in that case. The state Y (4260) has also been advocated as an example of a (cc) 8 g hybrid [11]. However, evidence of two resonant structures in e + e − → π + π − h c in the Y (4260) region has been presented by the BESIII collaboration [12]. This would imply in the hybrid picture the existence of two states, nearby in mass, having different cc-spins. Another analysis of the Y states in the diquark approach can be found in [13]. The univocal prediction of tetraquarks is that here must be only four Y states arising from the orbital angular momentum L = 1 and no radial excitation, as shown in Table I. Parity requires L to be odd, and charge conjugation requires the spin structure of these states to be symmetric under the exchange diquark ↔ antidiquark, due to the factor (−1) L introduced by the exchange of coordinates. Besides those of Table I, there are two other spin structures possible which are antisymmetric and in P wave they would give J P C = 1 −+ . One obtains J P C = 1 −− also from L = 3, but this state is expected to be considerably heavier and we do not consider it.
In [14], the four basic L = 1 resonances with J P C = 1 −− of the diquark-antidiquark spectrum were identified with Y (4008), Y (4260), Y (4290) (a broad structure in the h c channel), or Y (4220) (a narrow structure in the same channel) and a possible ΛΛ resonance around 4630 MeV. The Y (4360) and the Y (4660), also known at that time, were tentatively considered to be n = 2 radial excitations of Y (4008) and Y (4260), respectively. Since that paper appeared, the experimental situation has evolved considerably. The status of the Y (4008) is no longer established, and the Y (4260) is now claimed by BESIII as a double humped structure [15], which is resolved into two resonances: a lower component, Y (4220), with observed decays into h c π + π − and χ c0 ω, and a higher component, Y (4330), which decays into J/ψ π + π − . On the other hand, it was also observed that Y (4630) and Y (4660) could be fitted as a unique resonance, mainly decaying into ΛΛ [16].
SI assumes Y (4008) to be a real resonance and it features Y (4260) as a single state. Belle [17] finds that data are better fit with two resonances, Y (4260) and Y (4008), and the width of the Y (4008) is found to be a factor 2 larger than that of the Y (4260). However, the Y (4008) has been seen so far by Belle only and current analysis of this resonance from BESIII is inconclusive [15]. On this basis, SII discards the Y (4008) and it features the two lines resolving the Y (4260), according to BESIII [15], Y (4220) and Y (4330), as lowest resonances. The Y (4360) and Y (4390) appearing in SI and SII respectively, are considered as the same resonance seen in different experiments. Similarly, in both SI and SII, one considers Y (4660) and the proposed Λ −Λ resonance at 4630 MeV to correspond to the same state.
The spectra of the two scenarios extend over 400 − 600 MeV, and one could wonder if this is consistent with these states belonging to a single, fine structure multiplet, given that P -wave charmonium states are all within an interval of about 100 MeV. However, the difference can be defended by the different composition of the tetraquarks w.r.t. the charmonia, as explained below.
First, in tetraquarks the total quark spin goes up to S = 2, which amplifies the range of the spin-orbit and tensor couplings, and the effect of the tensor force in tetraquarks has not been investigated so far. Secondly, the constituents of Y states are diquarks and antidiquarks with spin 0 and 1, the "good" and "bad" diquarks in Jaffe's terminology, see [9]. For S-wave tetraquarks, the mass difference between Z(4020) and Z(3900) results in a mass difference of "bad"to "good" [cq] diquarks of about 120 MeV, and the Y states contain from zero to two "bad" diquarks. This splits the masses of the different components of a tetraquark multiplet considerably more than in the conventional charmonia. Note that, the diquark mass difference is related to the spin-spin coupling between the charm and the light quark in the diquark, which comes out to be 3 − 4 times the c − q coupling in the charmed baryons [14]. In QCD, these couplings are proportional to the quark overlap probability, |ψ(0)| 2 , and the result is simply understood to indicate a closer packing of diquarks in the tetraquarks than in baryons.
The upper range of masses in the two scenarios goes into the region where radial excitations of the lowest Pwave tetraquarks are expected and one may wonder if the highest and, possibly, the next to highest Y states may be the radial excitations of the two lowest ones. For definiteness we assume this not to be the case.
It is still possible, however, that a better experimental resolution may substantiate the differences observed between the 4360 − 4390 and the 4630 − 4660 peaks so as to indicate the presence of one or more radial excitation in the region, as was assumed in [14]. We examine the issue of possible radial excitations in Section IV. A similar issue has been raised for the excited Ω c = css states, whose mass spectrum has been measured by the LHCb collaboration [18], and confirmed recently by Belle [19], except perhaps the Ω c (3119). The LHCb mass spectrum is discussed in a number of papers [20][21][22][23], following the analysis of [20], in which all observed five states are assumed to be P -wave c quark and ss diquark. Also here, the highest mass states overlap with the 2S, positive parity, radial excitations of the S-wave Ω c [24,25]. This issue will be clarified as and when the J P quantum numbers of the excited Ω c states are experimentally determined.
Within the two Scenarios given above, we work out the mass spectrum derived from the spin-orbit, spin-spin and tensor coupling interactions, the latter was not included in [14]. The principal aim is to investigate whether the tetraquark picture may provide a satisfactory description of the presently determined Y -states in the cc sector, eventually distinguishing between SI and SII. For instance, since parameters are obtained from the solution of a second order equation, the mass formulae could produce complex parameters. In the case that there is no real-parameter solution possible, the underlying theoretical framework, namely the diquark picture, can be excluded as a physical template for the Y states. We use the reality condition on the parameters to eliminate some alternative assignments, as discussed below.
On the positive side, one would expect the value of the chromomagnetic, spin-spin coupling inside the diquark, [κ cq ] P , to be close to the analogous parameter derived for the S-wave tetraquarks, which is [κ cq ] S 67 MeV as discussed in [14]. For tight diquarks, this parameter should not be too much affected by the addition of one unit of orbital angular momentum. In addition, comparison with the S-wave tetraquark masses can give the energy for the excitation of one unit of orbital momentum.
The diquark mass being very similar to the charm quark mass, one expects the excitation energy of the tightly bound Y states to be similar to the one obtained from the comparison of the P and S-wave charmonia and for the P -wave Ω c states. This expectation is indeed satisfied by one solution in each of the two scenarios.
The tetraquark scheme predicts several other negative parity states with different J, and another 1 −− resonance arising from L = 3. Ignoring the L = 3 state, which is presumably rather heavy, we comment in the end on the composition of the full L = 1 supermultiplet and give an estimate of the expected masses.
The plan of the paper is as follows. In Sect. II, we repeat the analysis of the five L = 1 charmed-baryons Ω c , whose mass spectrum has recently been measured by the LHCb collaboration [18], following the analysis of [20]. In this connection, we offer an alternative calculation of the tensor couplings in terms of the Wigner's 6-j symbols.
In Sect. III, we introduce the two scenarios compatible with the present data, derive the mass formulae for the Y states and obtain the parameters of the Hamiltonian from the mass spectra. The role of radial excitations is discussed in Sect. IV. Error analysis and correlations among the parameters are presented in Sect. V. Results are discussed in Sect. VI.
A picture of the full L = 1 multiplet is reported in Sect. VII and Conclusions are presented in Sect. VIII.
In Appendix A we derive the tensor couplings for the Ω c baryons and diquarkonium, using Wigner's 6-j symbols. Correlation matrices in the analysis of the data on the Y states are given in Appendix B.

II. EFFECTIVE HAMILTONIAN FOR THE
Ωc BARYONS In the diquark-quark description, the Hamiltonian for the Ω c states can be written as [20]: In Eq. (1), m c and m s are the masses of the c and the s quarks, respectively, κ ss is the spin-spin coupling of the quarks in the diquark, and L is the orbital angular momentum of the diquark-quark system. The coefficients a 1 and a 2 are the strengths of the spin-orbit terms involving the spin of the diquark S [ss] and the charm-quark spin S c , respectively, c is the strength of the spin-spin interaction between the diquark and the charm quark, and S 12 /4 represents the tensor interaction, defined by where S 1 and S 2 are the spins of the diquark and the charm quark, respectively, and n = r/r is the unit vector along the radius vector of a particle. The scalar operator of Eq. (2) can be expressed as the product 3S i 1 S j 2 N ij , where the tensor operator is To compute the matrix elements of this operator between states with the same fixed value L = 1 one uses the identity from Landau and Lifshitz [26] (see also [27]): One finds which requires the matrix elements S , L; J|L·S 1,2 |S, L; J .
The latter can be computed by applying the operators L·S 1,2 to the products of three angular momentum states, see [20]. More effectively, one can use Wigner's 6-j symbols (now implemented in computer codes), as is customary for analogous cases in atomic and nuclear physics and is explained in Appendix A.
In either way, one obtains: J = 5/2 : After diagonalisation, we get the mass corrections arising from the Hamiltonian (1). We remind that in all the five states, there is the common mass term In order to determine the parameters a 1 , a 2 , b and c, Karliner and Rosner [20] have used the spin averaged mass and have worked with the mass differences of the five Ω c states. We reproduce their values, summarised in Table II  with J P C = 1 −− are obtained for L = 1, 3. Spin wave functions are given in Table I We extend the Hamiltonian of P -wave tetraquarks given in [14] by including the tensor coupling contribution S 12 is defined as in Eq.
(2), with S 1,2 representing the spins of the diquark and antidiquark, and where N 1 is the number of spin 1, "bad", diquarks. Comparing to (1), we see that in this case the coefficients a 1 and a 2 are a 1 = a 2 ≡ 2a Y due to charge conjugation invariance. The spin-spin interaction between diquark and antidiquark is neglected here since in P -wave the overlap probability is suppressed [14]. In the Ω c case, the spinspin interaction, represented by c, Table II, is similarly suppressed and the same happens in P -wave charmonia. The calculation of the matrix elements of the L · S X operator, with S X = S [cq] and S [cq] , is described in Appendix A, Eqs. (A23) and (A25). We note here that: • tensor couplings are non vanishing only for the states with S Q = SQ = 1; • the operator L · S Q is not invariant under charge conjugation and it does mix the states Y 3 and Y 4 , with a J P C = 1 −+ state with the composition: • Y (+) appears as an intermediate state in the products (L · S Q )(L · SQ) and (L · SQ)(L · S Q ), giving contribution to both diagonal and non diagonal terms; charge conjugation invariance is of course restored when making the sum of the two products, which is block diagonal in the basis (Y 3 , Y 4 ) and Y (+) .
In conclusion, we have to consider the full (3 × 3) matrix L · S [cq] . Using (A10) and (A19) we find: where S, S = 0, 1, 2 and the curly brackets denote Wigner's 6-j symbols [30]. After Eqs. (5) and (9), tensor couplings over the Y 3 − Y 4 states are represented by the non-diagonal matrix: The eigenvalues of the mass matrix of Y -states derived from Eqs. (9) and (13), are written as: and M 00 = 2m Q + B Q . We have made explicit that the states Y 1,2 in Table I are eigenstates of the mass matrix, while M 3,4 are the eigenvalues of the matrix with In the scenario SI, we take the four J P C = 1 −− Y states to be Y (4008), Y (4260), Y (4360) and Y (4660), with masses (all in MeV) Masses are taken from PDG [28], except for the Y (4008), which is from Belle [17].
In the scenario SII the masses of Y (4220), Y (4330), Y (4390) and Y (4660) are (all in MeV): i.e. the state with the mass M 4 is the same as in SI.
For S-waves [14], the spin-spin interaction gives a larger mass to S = 1 diquarks with respect to S = 0 ones and the same for antidiquarks. For this reason, it is natural to associate Y 1 and Y 2 with the two lightest particles in increasing mass order, M 1 < M 2 . This association is forced by the fact that if we exchange the role of Y 1 and Y 2 , we obtain only complex solutions for the parameters of the Hamiltonian.
In the case of SII, the association agrees with the fact that Y 1 has a sizeable component over the state with vanishing c −c spin, in agreement with the observed decays of Y (4220) into h c , while Y 2 has pure c-c spin equal one, in line with the observed decays of Y (4330) into J/ψ. The assignment allows also to describe the decay Y (4330) → X(3872) + γ as an allowed electric-dipole transition, given that Y 2 has the same spin structure as one attributes in the model to X(3872).
On the other hand, for Y 3 and Y 4 , both containing two spin 1 diquarks, we shall allow both possibilities: Before proceeding to the estimate of the values of the parameters M 00 , a Y , κ cq and b Y , we first note their possible interdependence on each other. From Eq. (16) for M 4 − M 3 follows that this mass difference is invariant under the simultaneous sign change (a Y , b Y ) → (−a Y , −b Y ). Hence, from this mass difference alone, we have two solutions: a Y < 0 and a Y > 0. We shall call them case 1 and case 2, respectively. In line with the analysis for the Ω c states, given in Table II, only a Y > 0 should be kept. This is also the choice suggested by the natural mass ordering, in which the J = 3 state should have a higher mass than the J = 1 states. So, the only physically acceptable solution is the one which has positive value of a Y irrespective of the sign of the value of b Y .
However, as the errors on some of the masses are large, we shall see below that, including the errors, solutions whose central values have a Y < 0, are also allowed. In addition to Eq. (16), the mass difference M 2 − M 1 provides a constraint on the parameters a Y and κ cq : Thus, in both the scenarios for the Y i masses, κ cq > a Y , with the two approaching each other as this mass difference decreases.
The central values of the parameters a Y , b Y , κ cq , and M 00 are determined from the masses given in Eq. (17) for SI and in Eq. (18) for SII and presented in Table III. In each scenario, we indicate with (c1) and (c2) the two solutions obtained from the secular equation for Y 3 and Y 4 .

IV. VARIATIONS ON THE THEME
We briefly comment on radial excitations, considering a Scenario III, proposed in [5], which envisages the confirmation of Y (4008) and the doubling of Y (4260): Given the masses of the first three states, we obtain the parameters of the Hamiltonian as functions of the mass of the fourth ground state, M 4 . The parameters are real if this mass is such as to make positive the radicand in Eq. (16). Numerically, this implies: This is consistent with the fourth ground state being Y (4660), with the parameters very similar to those of Scenario I, but not with Y (4390), which has to be the radial excitation of Y (4008). Yet another possibility is to take the first three states of the Scenario II and leave undetermined M 4 . The reality condition gives a result close to (20), however with a parameter κ cq a bit smaller than expected from the S-wave masses. Assuming the fourth state to be the Y (4660), we go back to the Scenario II with κ cq 54 MeV, which is in the acceptable range. In both alternatives considered, the full range of masses in Scenarios I and II is acceptable for P -wave ground states.

V. ERRORS AND CORRELATIONS
To work out errors and correlations among the parameters, we have used the method of least squares to determine the best-fit values and the covariance matrices. For this, the χ 2 -function is calculated. In general [28], where y = (y 1 , . . . , y N ) is the set of the experimentally measured values which are assumed to be independent and ∆y i are their variances. The quantities µ i ( θ) are dependent on the unknown parameters which are collected as the vector θ = (θ 1 , . . . , θ m ) where m ≤ N . For the problem at hand, we take the parameter-dependent functions from Eq. (14), µ i ( θ) = M i , where i = 1, · · · , 4, and The best-fit estimations of the parameters θ k , obtained after minimising the χ 2 -function, are presented in Table III, as central values. Note that each scenario results into two solutions which differ by the sign of the best-fit value of a Y , in line with the discussion above. The variances of the parameters are also shown in Table III the correlation matrices are collected in Appendix B. The parameters (22) are strongly correlated as all the correlation moments in the corresponding matrices (B1)-(B4) are close in magnitude to unity. To show this, we plot two-dimensional confidence level (C.L.) contours involving some of the coefficients. The correlations among the parameters a Y and b Y in terms of the 68.3% (χ 2 = χ 2 min + 2.3 for two degrees of freedom) and 95.4% (χ 2 = χ 2 min + 6.18) C.L. contours are presented in Fig. 1. Similar contours demonstrating correlations among a Y and κ cq are shown in Fig. 2. The first and the second rows in these figures correspond to the scenario I and II, respectively, and in each row, the left panels are plotted for the negative best-fit value of a Y (c1), while the right panels are for the positive best-fit value (c2). Our analysis shows that the scenario I (c1) is not tenable, as, within 95.4% C.L. and even higher, a Y remains negative. Thus, the requirement of positive a Y disfavours case 1 in the scenario I. In the scenario II (c1), small positive values of a Y are allowed with a relatively large probability. In (c2), large positive values of a Y are predicted for both the scenarios I and II.

VI. DISCUSSION
With the current uncertainty of the experimental scenarios and many parameters one cannot draw quantitative conclusions, except observing that the values of the parameters are qualitatively similar to those derived in the P -wave Ω c -states in three of the four solutions. One can, however, underline two criteria that lead to some preference for the scenario II.
The first is the value of the chromomagnetic coupling κ cq . We expect the fitted parameter to be close to the analogous parameter derived for the S-wave tetraquarks, which is [κ cq ] S 67 MeV as discussed in [14]. Indeed, there are no reasons to believe that the chromomagnetic coupling κ cq in the diquark should change with the addition of one unit of orbital angular momentum. At 95% C.L., the allowed value of κ cq from the Y states in the scenario II (c1) comes out somewhat smaller than anticipated, while it is somewhat larger in the scenario II (c2). (See, the lower two frames in Fig. 2). Thus, this criterion would favour the scenario II. A second expectation is for the Hamiltonian in Eq. (9) to describe both S and P -wave states, with the same value of the diquark mass. As commented in [14], Y 2 , which in SII corresponds to Y (4330), is in the same spin state as the X(3872) except that there is a gap in mass between the two, which here is fully accounted by B Q and by the spin-orbit interaction. If this is the case, one can derive the excitation energy of one unit of orbital momentum from the equation Using the input from Table III We may see what happens if we force the diquark spinspin couplings to be equal, [κ cq ] P = [κ cq ] S = 67 MeV, The χ 2 analysis is redone with four experimental input values and three unknown variables θ k = (M 00 , a Y , b Y ). We have one degree of freedom and can discriminate minima according to χ 2 min . The best-fit values and variances of the parameters M 00 , a Y , and b Y (all in MeV) corresponding to the minima with the lowest χ 2 min in each scenario are reported in Table IV. There are other minima in both the scenarios, but their χ 2 min are larger, and hence we don't discuss the resulting parameters.
With the parameters from Table IV we obtain: The value for B Q from the analysis of the Ω c -baryon resonances in the diquark model can be obtained from the expression for M 0 , given in Eq. (8). Using the input values of the quark masses and κ [ss] from [29] yields The values obtained for B Q in Eqs. (25) and (26) can be compared with the orbital angular momentum excitation energy in charmonium, given by the analogous formula (27) The combination of J/ψ-and η c -meson masses eliminates the contribution of the S-wave spin-spin interaction in the J/ψ-meson, absent in the h c , which has S cc = 0.
The similarity of the results for different P -wave systems with subcomponents in color 3 and3 is interesting and may suggest more precise calculations for the Y states, with the potential methods applied successfully to charmonia.
The states with masses M 3 and M 4 are linear combinations of Y 3 (S = 0) and Y 4 (S = 2). We note that in both SI and SII, the eigenvectors corresponding to M 3 and M 4 in c1 are close to S = 0 and S = 2, respectively, while in c2, it is the opposite, i.e., they are close to S = 2 and S = 0, respectively. Table V (column 1 and 2) gives the components of the eigenvector associated with M 4 , which is called v 4 , for different scenarios and solutions. The orthogonal vector v 3 is not shown. The eigenvectors carry interesting information; the projection of the eigenvector on the state with cc spin =1 is related to the probability of this state to decay into a J/ψ (S cc = 1) rather than in h c (S cc = 0). The fourth column gives the probability of finding S cc = 1 in v 4 . The table indicates that Y (4660) in solutions c2 should have a good probability to decay into h c while in c1 the J/ψ should dominate. This is quantified in the entries in Table V (third column).

VII. THE FULL L = 1 SUPERMULTIPLET
We expect many particles in the supermultiplet of L = 1 diquarkonia, analogous to the χ-states of char-  The total number of states coincides with the total number of quark spin and orbital momentum states, i.e., 2 4 × 3 = 48, as one verifies easily. Spin compositions are given in Tables VI and I, and tentative masses are presented inTable VI.
Indications exist for two 0 −+ states. However, in the same channel there should appear two conventional radially excited charmonia, η c (3S) and η c (4S), for a total of four states, with possible mixing and corresponding distortions of the spectrum.

VIII. CONCLUSIONS
We have analysed the masses of the four lightest Y states, using two experimental scenarios proposed in [1,2] and the effective Hamiltonian appropriate for L = 1, J P C = 1 −− tetraquarks, already introduced for the P -wave charmonia and for the excited Ω c states.
The current uncertainties on the spectrum of Y -states hinder us to reach a completely quantitative conclusion. However, we find (i) the coefficient of the spin orbit interaction to be positive, within errors, and comparable to that found for the Ω c states; (ii) the mass difference of the "bad" and "good" diquarks to be similar to what was found previously for the S-wave tetraquarks; (iii) the energy of the orbital excitation is found to be quite comparable to the values for charmonia and Ω c ; (iv) at variance with the latter cases, the coefficient of the tensor coupling turns out to be large and negative. The scenario with five Y states proposed in [5], including Y (4008) and the two components of the previous Y (4260), is also consistent if one assumes Y (4390) to be the radial excitation of Y (4008).
Features (i) to (iii) are coherent with our a priori expectations, while we have no particular objection or explanation of (iv). A slight preference for the Scenario II results, if one insists on enforcing the exact equality of VI: Spin composition, couplings and tentative masses (in MeV) of the particles in the P -wave supermultiplet, in addition to the states in Table I. Mass formulae are derived from Eq. (9), with M00, aY , bY taken from the second line of Table IV and κcq = 67 MeV. N1 is the number of spin-1 "bad" diquarks, defined in Eq. (10).
Mass(MeV) best fit Table IV 3 the mass difference of the "bad" and "good" diquarks in S and P wave states. Hopefully, some clarification on the composition of the Y states will be provided by BESIII, Belle II, and LHCb. With precise measurements, the parameters of the effective Hamiltonian can be determined more accurately, providing a quantitative test of the underlying diquark model. Tetraquarks require many more states in P wave other than the Y states reported in Table I Combining three angular momenta, j 1 , j 2 , j 3 to a given J, one may follow two paths, characterised by the values of the intermediate angular momentum to which the first two are combined, e.g. j 1 and j 2 to j 12 or j 2 and j 3 to j 23 , each path corresponding to different base vectors. In the formulae given below, these two bases are characterised as follows or, with a shorter notation where it is understood that j 1 , j 2 , j 3 and J are held fixed. Vectors in the two bases are, of course, related by a unitary transformation: Besides j 12 and j 23 , the Clebsch-Gordon coefficients C depend upon the angular momenta that are being held fixed, j 1 , j 2 , j 3 and J, that is the Cs depend on six angular momenta. To maximise the symmetry properties, one rewrites Eq. (A3) as [30]: Wigner's 6-j symbols are represented by the curly brackets. They appear in the calculation of the matrix elements of the spin-orbit Hamiltonian or the tensor coupling for two particles with spins S 1 and S 2 and different masses in the orbital angular momentum L. Examples are the P -wave Ω c baryons and the diquark-antidiquark tetraquarks in P -wave, considered in the present paper. In these cases, to classify states it is convenient to couple S 1 and S 2 to a total spin S and couple S to L to obtain the total J, that is: In this basis the matrix elements of the total spin-orbit operator are easily computed according to the formula: In the spin-orbit interaction and in the tensor coupling, however, one encounters the matrix elements of the operator L · S 1 = j 1 · j 2 , which would require a complicated calculation based on writing explicitly the states as products of three angular momentum states and applying the operator L · S 1 to them.
A more convenient way to proceed is to use Eq. (A4) and set j 1 = L, j 2 = S 1 , j 12 = L + S 1 , j 3 = S 2 . (A7) In this basis, is diagonal on the basis vectors.
Tables of 6-j symbols can be easily implemented in a computer code and they are already available, making use of the command SixJSymbol[{j 1 , j 2 , j 3 }, {j 4 , j 5 , j 6 }], in the symbolic computer algebra system Mathematica [31]. Therefore the result in Eq. (A10) can be obtained with a program of a few lines [32]. In the following, we give the explicit formulae for the cases considered in the paper. We will call j 12 = j LS [ss] and j 23 = S = 1/2, 3/2. We have to consider the matrix L · S [ss] in the two cases: J = 1/2 and J = 3/2. In the J = 1/2 case, j LS [ss] can take the values 0, 1 and Eq. (A10) reads: where S, S = 1/2, 3/2. This sum can be calculated easily when the required values of 6-j symbols are known: Therefore, in the basis of states ( 4 P 1/2 , 2 P 1/2 ) we have (the notation 2S+1 P J is the same as in [33]): In the same way with the 6-j symbol values: giving in the basis ( 4 P 3/2 , 2 P 3/2 ) Both results agree with [20].  Note that the state with S = L = 1 has positive charge conjugation, C, opposite to the value of C of the other two states and of the Y states.
The spin-orbit coupling must be even under C and, therefore, it is represented by which is diagonal on the states with S = 0, 2. However, the C-even combination of the spin-orbit couplings appearing in the tensor coupling is where S, S = 1, 2. For the matrices we obtain: (A30) In addition to the values of the 6-j symbols presented in (A24), for the derivation of (A27) and (A30) one needs to have two more:

Appendix B: Correlation Matrices
In this appendix the correlation matrices in the analysis of the data on Y states are collected. We label them in accordance with the notations used in Table III. SI (c1):