Magnetic dipole moments as a strong signature for $\alpha$-clustering in even-even self-conjugate nuclei

We investigate the magnetic dipole moments in even-even self-conjugate nuclei from ${}^{12}\mathrm{C}$ to ${}^{44}\mathrm{Ti}$. For the latter, the measured gyromagnetic factors of excited states turn out to assume the same value of $g \approx + 0.5$ within statistical errors. This peculiar feature can be interpreted on the basis of collective excitations of $\alpha$-clusters. Analogously, the behaviour of the same observable is studied for all isotopes obtained by adding one or two neutrons to the considered self-conjugate nuclei. It is found that for the $N = Z + 1$ isotopes the $\alpha$-cluster structure hardly contributes to the observed negative g- factor value, corroborating molecular $\alpha$-cluster models. The addition of a further neutron, however, restores the original $\alpha$-cluster g-factors, except for the semi-magic isotopes, in which the deviations from $g \approx + 0.5$ can be associated with the relevant shell closures. Secondly, we analyze the same observable in the framework of a macroscopic $\alpha$-cluster model on a finite lattice of side length $L$. We focus on the discretization effects induced in the magnetic dipole moments of the $2_1^+$ and the $3_1^-$ states of ${}^{12}\mathrm{C}$ at different values of the lattice spacing $a$.


Preamble
The tendency of nucleons to congregate into clusters is known since the early days of nuclear physics [1]. The intrinsic stability of the 4 He nucleus alongside with the energy gap of 20.21 MeV with respect to the lowest single-particle excited state, makes the α-particle a candidate for a conglomerate of nucleons capable of surviving relatively unperturbed within the nuclear mean field, for a recent review see [2]. It is in fact well established that nuclear spectra are a result of the interplay between pairing and cluster correlations and the mean field generated by the individual nucleons [3]. In absence of correlations between the constituents, nucleons would move as independent particles and a well-pronounced shell structure (the Fermi surface) would appear [3]. As soon as the residual attractive interactions are turned on, nucleons orbiting in the same subshell but with opposite projection of the angular momentum along the symmetry axis of the nucleus (i.e. connected by the antiunitary time-reversal operation, see refs. [4,5]) form Cooper pairs [6,7] or give rise to larger subsystems. Among the latter, the α-particle dominates in nuclei with even and equal numbers of protons and neutrons, whose higher binding energy per nucleon reaches a local maximum in the Segrè chart. Cluster correlations are found to reduce the neutron skin thickness with respect to a mean field value [8].
Additionally, in the pioneering work in ref. [9] it was noticed that the binding energies of light α-conjugate nuclei scale linearly with the number of bonds among the α-clusters, sitting at the vertices of regular polyhedra (cf. fig. 1 in ref. [10] and refs. [11,12,13]). In particular, the adoption of the equilateral triangle (D 3h ) [14,15,16], the tetrahedron (T d ) and the triangular bipyramid (D 3h ) as basis-configuration for the inspection of the structure properties of the 12 C [17] and 16 O [18,19] within the Algebraic Cluster Model (ACM) [12,20] and of the 20 Ne [21] within a geometric macroscopic α-cluster model, respectively, has revived the interest for the subject in recent times (cf. refs. [22,23]).
Furthermore, α-clustering is reinforced as the excitation energy grows. As an example, in 12 C at around 10 MeV the nuclear density reduces to less than one third with respect to the initial value (cf. refs. [23]) and the inter-α separation increases accordingly. An example is represented by the Hoyle state at 7.65 MeV, credited by miscroscopic studies as a triangular bent-arm arrangement of α-particles [24,25] and earlier as a linear chain configuration [26]. On the other hand, the triangular ground state keeps a significant overlap with the p 3/2 closed subshell configuration (cit. refs. [3]), as the independent-particle in a mean field feature gains ground. In the opposite limit, the nucleus dissociates completely into the twelve constituents at around 10 2 MeV [3].
The idea according to which the α-particle structure is not manifest in the ground state but emerges gradually with the increase of the internal energy of the system is at the origin of the diagrammatic representation proposed in refs. [27,28] for the 12 C, 16 O, 20 Ne, 24  Therefore, the adoption of pictures that neglect the inner degrees of freedom of the 4 He clusters for the analysis of excited state properties, such as intrinsic magnetic dipole moments, of this class of nuclei appears justified. Motivated by the long standing experimental interest in magnetic moments for 12 C [29], 16 O [30], 20 Ne [31], 24 Mg [32], 28 Si [33], 32 S [34], 36 Ar [33], 40 Ca [35,36] and 44 Ti [37], this paper is dedicated to the theoretical interpretation of these observables throughout light and medium-mass α-conjugate nuclei. Moreover, we point out that in part of the aforementioned bibliography, e.g. the measurement of the g-factor of the 3 − 1 state of 16 O in ref. [30], the α-cluster nature that characterizes this nucleus was marginalized if not ignored during the discussion of the observed data.
The primary target of the present article is to fill this gap in the literature, showing that experimental g-factor values lie in a reasonable neighbourhood of g = +0.5 and constitute a new unambiguous evidence for an α-cluster structure in even-even self-conjugate nuclei. In particular, the treatment begins with the presentation of the values of the g-factors for low-lying energy levels of light and medium mass nuclei computed in the framework of the nuclear shell model (cf. ref. [38]). In parallel, we extract for comparison the same quantities by means of the Schmidt estimates [39,40]. Subsequently, we repeat the calculation by adding one or two neutrons to the original nuclei, highlighting how this operation alters the g-factor values. We further contextualize the outcoming nuclides in the realm of molecular nuclei [3,41,42], where the added neutron(s) either assume the role of valence particles (e.g. 10 Be) [3] or become an integral part of the cluster themselves (e.g. 12 Be) [3].
Secondly, the recent experimental advances as well as the surge of interest for α-conjugate nuclei calls for a study of magnetic dipole moments, µ, and g-factors in the framework of purely αcluster models as the one in refs. [43,44] on a cubic lattice. In particular, we display the behaviour of the average value of µ computed on the 2 + 1 and 3 − 1 states of 12 C with maximum angular momentum projection as a function of the lattice spacing. Moreover, in sec. 3 we demonstrate, how multiplet averaged [44] and isotropically averaged [45] values of the latter observable succeed in reducing the artifacts associated with finite lattice spacing in the two energy levels of the same nucleus recently analysed in ref. [46].

Magnetic dipole moments
Magnetic dipole moments are known to contain information on the microscopic structure of a nuclear system. These observables are indeed very sensitive to the occupancy of quasiparticle orbits of valence nucleons and serve as a testing ground for A-body wavefunctions and theoretical models in general. Besides, the measurement of the magnetic dipole moments represents an immediate verification for the assignment of a given angular momentum to nuclear ground and excited states in the shell model. In the many-body definition (cf. ref. [47,48]), the former are given by the superposition of the single-nucleon orbital angular momentum and spin operators [48], weighted by the spin and orbital gyromagnetic factors indexed by π (ν) for protons (neutrons).
Here, e denotes the unit electromagnetic charge and m n and m p the mass of the neutron and the proton, in order. For free nucleons in ref. [49] the g-factors assume the values g p = 1 g p s = 5.5856946893(16) , g n = 0 g n s = −3.82608545 (90) . ( In more refined calculations, g p s and g n s are quenched by a factor 0.60 ÷ 0.75 [50,51] in order to account for the polarization of the core of the nucleus and for the meson-exchange current (MEC) corrections [52]. Nonetheless, some authors [53] prefer to neglect these corrections and use the free values for charges and g-factors of nucleons in the microscopic calculations of electromagnetic moments and transition probabilities. Ab-initio estimates for the corrections to the g-factors can be extracted from the average values of the MECs in the zero-momentum-transfer limit in eq. (125) of ref. [54], derived in the framework of chiral effective field theory (χEFT). Note that the currents used by these authors have to be taken with a grain of salt [55]. For the Schmidt estimates (SE) presented in this section, we reduce the nucleon g-factors to 70 % of the free values and assume Z to be always even.

N = Z nuclei
We target the average values of the operator in eq. (1) in low-energy excited states of αconjugate nuclei. Motivated by the high stability of the 4 He nucleus, we suppress the single-nucleon degrees of freedom in spite of the former accomplishment and interpret the nuclear excitations in terms of the rotational and vibrational motion of the M ≡ Z/2 α-particles. These bosonic groupings are characterized by zero spin and isospin as well as vanishing magnetic dipole moment in the 0 + ground state. In such a macroscopic framework, the spin contribution to the magnetic dipole moment remains equal to zero in any excited state. Hence, we are allowed to set g α s = 0 and rewrite eq. (1) as where m α is the α-particle mass, the summation is performed over the M α-clusters, g α = 2 is the orbital component of the gyromagnetic factor for the α-particle and coincides with the nuclear g-factor. Its small deviation with respect to 0.5 is essentially due to the mass excess of the 4 He nucleus.
The available experimental values for g (α) L in light and medium-mass N = Z even-even selfconjugate nuclei are all compatible within the statistical errors with the outcome of eq. (3) for macroscopic models. This remarkable agreement emerges dramatically in fig. 1 and constitutes the cornerstone of the present interpretation. The validity of this assertion is not limited to the lowest 2 + or 3 − states of α-conjugate nuclei, but remains intact at growing excitation energies.

Nucleus Level Energy [MeV]
g (SM ) g (α) g (exp) Z/A Exemplary, in this respect, is the case of the 2 + 1 and the 4 + 1 levels of 20 Ne and 32 S as well as the 2 + 1 , 4 + 1 , 2 + 2 and 4 + 2 states of 24 Mg, where the measured g-factors do not display systematic deviations from the value in eq. (3) within their sizable errors (cf. tab. 1). In the framework of the shell-model (SM), g-factor values in reasonable neighbourhood of the macroscopic α-cluster value of ≈ +0.5034 are found. For 10 ≤ Z ≤ 16 nuclei the theoretical SM data are obtained through the phenomenological USDB interactions [65], taylored for sd -shell nuclei [66,67]. Conversely, for 36 Ar the WBT interaction in ref. [68] is used in the same model space, whereas for 40 Ca the calculations [69] assume admixtures of random phase approximation 1p-1h states and deformed 3p-3h states [75]. However, the SM configuration d −1 3/2 f 7/2 alone yields in 40 Ca g-factors in fair agreement with the experimental and large-scale theoretical estimates in ref. [69]. Eventually, for 44 Ti the FPD6 [70] nucleon-nucleon interaction is implemented in the full pf -shell model space, assuming a self-conjugate 40 Ca core [37].
Furthermore, for states of zero total isospin, as the bosonic ones, generated by excitations of α-particles (cf. tab. 1), the g-factor is predicted to correspond to its isoscalar component, where µ N is the nuclear magneton and µ J (T z = ±T ) are the magnetic dipole moments of two states with total angular momentum J and opposite z-component of the total isospin, i.e. mirror nuclei. In eq. (5), the charge dependence of nuclear force and the nucleon masses has been also ignored, as in eq. (3.2) of ref. [71]. The average value of S z in the former equation is computed among the states with maximum total angular momentum and isospin projection along the z-axis, Moreover, for purely collective excitations constructed on the J = 0 ground state, the total angular momentum coincides with L and the average value of the third component of the spin in eq. (3.2) of ref. [71] vanishes. As a consequence, the isoscalar part of the gyromagnetic factor becomes independent on J and reduces to 1/2 [72] (cf. eq. (2)). This result agrees with eq. (3), the deviations owing to the neglected charge dependence of nuclear forces and nucleon masses.
Although the original measurements of the gyromagnetic factors of the 3 − 1 state of 16 O and the 2 + 1 state of 24 Mg yielded g = +0.55(3) [73] and +0.51(2) [74], respectively, the more recent and more precise observations in refs. [30] and [60] highlight a 10 % deviation from the macroscopic α-cluster value in eq. (3). This bias is only marginally covered by the shell-model calculations in ref. [76], where configuration mixing among one-body states lying within major harmonic oscillator shells has been considered.
Nevertheless, in ref. [56] the rest of the discrepancy is successfully attributed to isospin mixing exerted by the T = 1 (isovector) 3 − 4 state at 13.26 MeV, whose reduced magnetic dipole transition probability to the 3 + 1 state is large. Interpreting the two 3 − states as an admixture of 1p-1h neutron and proton states of an α-particle such as p −1 1/2 d 5/2 and indentifying the mixing source with the Coulomb force V C , an isovector contribution with the expected magnitude is obtained, where |∆E 3 − | denotes the energy separation between the two levels and the compact notation |J, T ≡ |JT J z T z for the shell-model states is understood. In summary, the full gyromagnetic factor of the 3 − 1 state, g (0) J , becomes +0.555 [56]. The same considerations, applied to the 2 + 1 (T = 0) state at 4.44 MeV and the 2 + 4 (T = 1) state at 16.11 MeV of 12 C, result in a smaller isovector correction to the g-factor of about 1.5 % with respect to the shell-model value (cf. tab. 1).
Altogether, the main nuclear structure lends weight to the consistency of α-clustering.
Concerning the 2 + 1 state of 24 Mg, the gyromagnetic factor from the USDB Hamlitonian in ref. [65], together with the isovector correction (cf. ref. [77]) yields +0.520, is consistent with the most recent experimental counterpart [32]. The overlap improves when also the effects of the MECs on the bare nucleon g-factors in eq. (2) are included [78], leading to a final value of +0.544 [32]. For the other values of the gyromagnetic factors reported in tab. 1, times are not ripe for singling out possible deviations from the macroscopic α-cluster value in eq. (3), since the statistical errors affecting the measurements overcome on average the 10 % of the g-factor value itself.
Despite the sizable experimental errors affecting the available estimates and the small deviations from g (0) J , our α-cluster interpretation of the g-factors for this class of nuclei remains solid. Underpinning this construction is the gradual emergence of α-clustering at increasing excitation energy and, remarkably, close to α-decay thresholds [27]. Although yet unmeasured, the magnetic moments of the 3 − 1 and the 2 + 2 states of 12 C as well as the 3 − 5 and 5 + 1 states of 16 O in the close vicinity of the 3α-and 4α-decay thresholds, respectively, are expected to adhere even more faithfully with the predictions of the macroscopic α-cluster picture. Besides, as shown in sec. 2.2, the observed independence of the g-factors on the angular momentum and the excitation energy of the states to which they refer, is a peculiar feature, which the N = Z + 1 and Z + 2 nuclei with even Z do not preserve.

N = Z + 1 nuclei
In the following we consider the case in which a single neutron is added to an α-conjugate nucleus. The major question that naturally arises is whether the additional neutron enhances or weakens α-clustering. Although no general answer exists [3], clustered states are expected to appear in correspondence of the M α + n decay thresholds. For 9 Be, different approaches such as the linear combination of atomic orbitals (LCAO) [79] and the cluster shell model (CSM) [80] highlight an underlying dumbbell configuration with the added neutron moving in the cluster field with Z 2 symmetry [80].
Concerning 13 C, the states with excitation energy lying between 9.90 and 16.95 MeV are well reproduced by a molecular α-cluster model [81], in which the odd neutron lies between the 4 He nuclei arranged in a linear chain. Besides, for the (1/2) − 2 and (1/2) + 1 states of the same nucleus, the triangular α-cluster configuration with one neutron located in the same plane of the 4 He nuclei (π-orbit) and along the axis perpendicular to the 3α-plane (σ-orbit), respectively, has been proposed [81]. Additionally, the large isoscalar reduced electric monopole transition probability between the (1/2) − 2 and (1/2) − 3 states and the ground state [82] is also regarded as a signature of α-clustering in this nucleus [41]. Indeed, states lying below the 9 Be+α decay threshold are well predicted by the shell-model for p-shell nuclei [83].  Concerning 17 O, attempts to reproduce the spectrum of the nucleus above the 13 C+α decay threshold in terms of the 4α + n configuration have been made [84,85] in the framework of the generator coordinate method (GCM) [86]. Nevertheless, for the lowest-energy states, SMbased approaches as the cluster-orbital shell model (COSM) [87] seem to provide a rather faithful description [88,89]. Moreover, a recent CSM analysis of the whole low-energy spectrum and the electromagnetic multipole transition probabilities of 21 Ne [42] suggests, that particle and hole neutron states can coexist with the underlying 20 Ne core, thus exciting the internal degrees of freedom of one of the α-particles. The comparison with the experimental energies turns out to be favourable for the model, except for the missing K P = 1/2 + band starting at 17.34 MeV [42].
Besides, the g-factor of the (3/2) + 1 state, inferred from the magnetic moment [42] with the free value for g n s , is equal to −0.287, i.e. smaller in magntitude than the experimental one in tab. 2. The value is obtained by addition of the g-factor of the 5α bi-piramidal core (cf. eq. (3)) with the contribution of the unpaired neutron, −0.85 µ N , times suitable Clebsch-Gordan coefficients (cf. eq. (20) in ref. [42]).
Switching to heavier nuclei, the cross section of elastic and inelastic scattering of α-particles suggests that α-particles at least at the surface level are present and interact with the incoming α-particle [90]. Moreover, in 41 Ca and 45 Ti, α-clusters can be induced by the pairing interaction alone at large values of the strength parameter G or at physical values of the latter, provided a four-nucleon force is introduced in the Hamiltonian [90]. The experimental data of the gyromagnetic factors for this class of nuclei (cf. tab. 2) highlight large deviations from the α-cluster value in eq. (4), as a result of the significant spin contribution from the unpaired neutron. The fact itself that the magnetic moment is not parallel to the total angular momentum leads to strongly J-dependent values of the g-factor in excited states of the same nucleus. Exemplary is the case of the (1/2) − 1 and the (5/2) + 1 states of 13 C (cf. fig. 2), whose g-factors differ by almost two units and the sign.
As a consequence, we replace the macroscopic α-cluster result in eq. (3) by the Schmidt formula [39], that attaches the magnetic moment to the unpaired nucleon, where u = n or p. Due to the inequality between the orbital and the spin gyromagnetic factors, µ precedes about the angular momentum j = + s of the single-particle ground or excited state.
The reported values of the nuclear g-factors are extracted from the magnetic moments by taking the projection of µ on the total angular momentum operator of the nucleus, J, where J is assumed to coincide with the one of the single-particle state j and the inner products are evaluated via the Landé formulae, The description of µ, merely in terms of the unpaired particle, proved to be quite effective for nuclei in the vicinity of a closed shell [110]. In fact, the Schmidt estimates for the sd -shell nuclei improve and quite accurate predictions are obtained for the (5/2) + 1 state of 25 Mg, as well as for the (3/2) + 1 one of 33 S and the (3/2) + 1 and (7/2) − 1 of 37 Ar. Conversely, the values of g (SE) associated with the (7/2) − 1 ground state and the (5/2) − 1 excited state of 45 Ti, sharply disagree with the experimental data. It follows that a larger model space is required for this nucleus, lying in the pf -shell, in order to obtain order-ofmagnitude agreement with the observed g-factors. The small excitation energy of the (5/2) − 1 state itself is an indicator of collective motion.
Concerning the shell-model estimates for 9 Be, calculations with phenomenological interactions encompassing the active nucleons in the p-shell (0 ω space) as well as the excitations to the sd -shell (1 ω space) produce a g-factor for the (3/2) − 1 state differing by less than 10% from the experimental counterpart [91]. The deviation reduces to nearly 1 % in the microscopic 2α + n model in ref. [111], where g ≈ −0.779, thus lending further weight to the molecular treatment of the unpaired nucleon moving around the α-clusters.
For 13 C, the former shell-model calculation with phenomenological interaction yields a similar agreement with the experimental g-factor of the ground state, whereas a more recent one based on Warburton-Brown interactions underestimates the same g-factor by 10 % [112]. The closest available estimate to the measured g-factor in tab. 2 yields +0.707 [113] and is based on the 4 He core in the p-sd model space and with 2 − 3 ω excitations, adopting the bare values of the proton and neutron g-factors. The result for the (5/2) + 1 excited state turns out to be even closer to the experimental value, which is indeed well reproduced (g ≈ −0.62) by the minimal coupling model in ref. [114], that assumes a self-conjugate 12 C core.
On the other hand, the g-factors of the sd -shell nuclei up to 29 Si are obtained from USDB interaction [65], and some disagreement with the experimental data is detected only for the (3/2) + 1 and (5/2) + 1 states of 21 Ne [57]. The discrepancy is slightly filled by the ab-initio in-medium similarity renormalization group (IM-SRG) approach, yielding g = −0.443 and −0.140 for the two states, respectively.
Concerning the (3/2) + 1 of 33 S, the sd -shell model calculation in ref. [100] with the empirical Hamiltonian of ref. [66] reproduces the experimental g-factor within less than 10 % deviation.
Moreover, the gyromagnetic factor of the Schmidt estimate approaches the measured value, thus suggesting that the largest contribution comes from the odd neutron in the d 3/2 shell [100].
Assuming a 28 Si core with configuration mixing and effective g p s and g n s for the 1s 1/2 and 0d 3/2 shells, the g-factor of the analogous (3/2) + 1 state of 37 Ar has been computed in ref. [102], when the experimental counterpart was unavailable. From the different shell-model estimates for the magnetic dipole moments of tab. 1 in ref. [102] follows, that a better agreement with the observed value of this state is obtained when the g-factors of the free nucleons are adopted [102].
For the (7/2) + 1 state of the same nucleus, in tab. 2 the shell-model state with negative parity is constructed from a 36 Ar core with a single neutron in the pf -shell, exploiting the free nucleon g-factors [104]. Core polarization effects are included in the calculation, by considering both the 0 + 1 and the 1 + levels of the α-conjugate core. The resulting g-factor is compatible within four standard deviations with the experimental value.
In the 41 Ca case, the reference result arises from a large scale sd-pf shell-model calculation with 28 Si core and USD [67], and modified Kuo-Brown [115] interactions for the sd and pf shells, respectively, and LKS [116] interaction for the mixing of the two major shells. The g-factor in tab. 2 is found to reproduce with less than 10% deviation the experimental value, even if the free values for g p s and g n s are employed. Finally, the gyromagnetic factors of the (7/2) − 1 and (5/2) − 1 states of 45 Ti are drawn from the shell model calculation in the pf model space with phenomenological interactions in ref. [107], in which g p s and g n s in the pf shell are fitted to reproduce the experimental magnetic moment of the ground state of 41 Ca (and 41 Sc). The g-factors of the two levels of 45 Ti in tab. 2 turn out to be overestimated in magnitude, although less than the Schmidt estimates, but possess the correct sign. Nevertheless, the application of the Nilsson model in ref. [117] with oblate deformation yields a g-factor of −0.003 for the 7/2 − ground state [117], underestimated but sensibly closer to the measured counterpart.

N = Z + 2 nuclei
In the following, we add another neutron to the original even-even self-conjugate nuclei and investigate, how the two extra neutrons affect the inner α-cluster structure.
As in 9 Be case, for the 10 Be nucleus, a number of studies examining the low-lying spectrum and the reduced electric and magnetic transition probabilities exists in microscopic or semi-microscopic pictures. Among these, the LCAO [118] and the antisymmetrized molecular dynamics (AMD) [119] have established that α-clustering in 10 Be emerges already in the ground and excited states of the K P = 0 + band, where the two valence neutrons form structures equivalent to π bonds.
The separation between the 4 He clusters grows in the second K P = 0 + and in the negative parity band K P = 1 − , where the two neutrons give rise to σ-like molecular bonds [119].  Table 3: Nuclear g-factors for excited states of N = Z + 2 nuclei, including 10 C. As noticed in tab. 2, the ratios Z/A do not follow with the measured g-factors.
Moreover, the application of the AMD in combination of the GCM in ref. [120] unveiled that the 0 + 1 , the 2 + 1 and the 3 − 1 states of 14 C possess intermediate features between a triaxial α-clustered configuration and a SM state, in which the two neutrons lie in the closed p-shell, while the protons fill the p 3/2 level. On the other hand, the 0 + 2 , 2 + 2 and 4 + 1 have a significant overlap with a triangular α-cluster configuration, in which the two neutrons stand in the same plane of the 4 He clusters and fill two sd -like orbitals [120]. Linear-chain α-cluster states are, indeed, found to characterize the 0 + 5 , 2 + 6 and 4 + 6 states located in the vicinity of the α + 10 Be decay threshold [120]. Concerning 18 O, the levels 0 + 1 , 2 + 1 and 4 + 1 below the 14 C+α threshold in the GCM model in ref. [121] have a strong overlap with the SM 0 ω prolate or spherical configuration, whereas the states belonging to the second K P = 0 + band support the formation of the inner α+ 14 C structure [121].
A similar study [122] conducted on 22 Ne predicts the existence of a K P = 0 + α-clustered band as well as negative parity doublets located above the 18 O+α threshold. Although less abundant, some theoretical cluster studies for the heavier α-conjugate nuclei are available, such as the ones in the framework of the weak coupling approach of the α-particle (or hole) with respect to a 40 Ca core [123], or in the semi-microscopic algebraic cluster model (SACM) [124,125]. The existence of parity doublet bands with K P = 0 + and 0 − in the neighbourhood of the lowest α-decay threshold is considered as a signature for the decomposition of the nucleus into a core plus an α-particle [123].
Regarding the g-factors of excited energy levels, their values are evidently correlated to the shell-closures. In particular, from the behaviour of the double lines in fig. 2, we infer that in the semi-magic nuclei at the N = 8 (resp. 20) shell closures, namely 14 C and 18 O (resp. 38 Ar and 42 Ca) the effect of the two extra neutrons is maximum (cf. tab. 3). For these nuclei, single-particle estimates are expected to reproduce quite accurately the g-factors, even when the free values of g p s and g n s are adopted. For the other nuclei, located before, between or beyond the two shell closures (cf. tab. 3), no significant discrepancy is found with respect to the macroscopic α-cluster prediction, with the sole exception of the g-factor of the 2 + 1 state of 22 Ne. For all other nuclei, the experimental g-factors agree within one standard deviation with ≈ +0.5034. Additionally, for the 4 + 1 state of 22 Ne, the g-factor turns out to be again compatible with the prediction of eq. (4). Since the pairing force couples equal nucleons with opposite total angular momentum projection, this interaction is not responsible for the suppression of the isovector contribution to the g-factor. On the other hand, four-body correlations between pairs of equal nucleons can quench g (1) J and reduce the g (0) J to 1/2, due to the vanishing total spin.
The application of AMD + GCM in ref. [120] to the low-lying states of 14 C, suggests that the considered 3 − 1 level in tab. 2 has a spatial distribution, intermediate between a SM state with the protons in the p 3/2 subshell and the neutrons in the1p-1h excited configuration and a triangular α-cluster state. In fact, the single-particle Gaussian wave-packets of the two neutrons gather still around the origin [120], despite a slight overall triaxial deformation (β ≈ 0.22, γ ≈ 24 • ). Also, the GCM analysis in ref. [121] predicts a strong overlap with a SM state for the 2 + 1 state of 18 O in tab. 3. We can deduce that, when the wavefunctions of the two extra neutrons overlap significantly in the space, clustering is hindered [134] and the g-factor displays significant deviations from eq. (4). Otherwise, if the two neutrons are spatially more separated, e.g. when sitting between two different pairs of α-particles in the molecular picture [134], the g-factor approaches the value of ≈ 0.5034, as in α-conjugate nuclei.
Therefore, we can predict that at higher-lying energy levels of 14 C and 18 O such as the 2 + 2 state, where clustering is expected to be more pronounced [121,120], the g-factor reaches the one of the other levels of the open-shell even-even N = Z + 2 nuclei in tab. 3. The shell closures at N = 8 and 20 (cf. fig. 3), in fact, seem to favour a more compact spatial distribution of the two neutrons, that undermines α-clustering in the low-lying energy states of the nearby semi-magic nuclei, (cf. 14 C, 18 O, 38 Ar and 42 Ca in tab. 3).
Although of a milder type, the phenomenon takes place also in the two considered states of 22 Ne, where the measured g-factor of the 2 + 1 state differs by 35 % with respect to the macroscopic α-cluster value.
As in the previous section, we give the shell-model g-factors in tab. 2 and the Schmidt estimates, that we obtain through the quenched values of g p s and g n s . The g-factor of the excited states of this class of nuclei is now attributed to two neutrons lying on two distinct single-particle states with angular momentum j 1 and j 2 , which sums up to J.
Next, the g-factors corresponding to the two single-particle levels are combined together, giving the overall g J factor through the formula Applying eq. Concerning the SM estimates, the gyromagnetic factor of the 2 + 1 state of 10 Be, calculated via the Monte Carlo shell-model (MCSM) with nucleon-nucleon interactions drawn from Chiral Effective Field Theory in the unitary correlation operator method (UCOM) [135,136], turns out to reproduce faithfully the value of eq. (4), thus suggesting a well-developed α-cluster structure, with the two neutrons acting as valence particles in a π orbital of a binary molecule (cf. ref. [119]).
Although 10 Be is a rather long-lived isotope, the magnetic moments associated with its excited states have not yet been measured.
For 14 C, the SM calculation, reported in ref. [128], is based on a 16 O core with three holes in the p-shell and one particle in the sd -shell. Besides, for the nucleons in the p-shell, the Cohen-Kurath two-body matrix elements [83] are adopted, whereas the Gillet interaction [137] is exploited between the p and sd -shell nucleons. The result is closer in magnitude with respect to the more recent estimate in ref. [112], obtained with the phenomenological V M U , SFO and SDPF-M interactions.
The data obtained from the phenomenological USDB interaction [65] for the sd -shell (cf. tab. I in ref. [57]) permit to reproduce quite accurately the g-factors of the 2 + 1 state of 18 O as well as the 2 + 1 and 4 + 1 states of 22 Ne, assuming the free value of g p s (g n s ) and an effective charge of e p = 1.5e (e n = 0.5e) for the protons (neutrons). In contrast, for the 2 + 1 state of 26 Mg and 30 Si, less agreement with the central values of the experimental data is attained [57]. In both cases, a more satisfactory answer is provided by the IM-SRG approach, that yields g = +0.512 and g = +0.420 for the two nuclides, respectively [57], both more compatible with the experimental values.
Regarding 34 S, the g-factor obtained in the large-scale SM calculation with sd -pf configuration mixing in ref. [76] is in excellent agreement with the observed result, whose central value is the closest to the macroscopic α-cluster value and its statistical error is smaller than for the neighbouring nuclei in tab. 3. The mixing between the sd and pf major shells has not been considered in the 38 Ar case, where the g-factor of the 2 + 1 is obtained entirely from the sd model space with the Warburton-Brown interactions [68], but with the same convention for the nucleon charges as well as for g n s and g p s . The result is still compatible with the experimental g-factor (cf. tab. 3).
The large-scale sd -pf SM calculation in ref. [132] produces a compatible value for the g-factor of the 2 + 1 state of 42 Ca, which the Schmidt estimate in tab. 3 reproduces with the same magnitude but opposite sign. Underpinning this calculation is the enlarged model space, that assumes the α-conjugate nucleus 28 Si as a core and the USD interactions [67] in the sd -shell, the modified Kuo-Brown ones [115] in the pf-shell and the LKS potentials [116] between the two major shells.
The 40 Ca core excitations are found to provide substantial improvement with respect to the pfshell estimates reported in ref. [132], thus indicating a deviation from sphericity in the 40 Ca core.
Eventually, the pf -shell model (FSM) calculation in ref. [133] with modified Kuo-Brown (KB3) interactions [115] sensibly underestimates the g-factor of the lowest 2 + state of 46 Ti, as a result of the limited size of the model space. In 45 Ti ( 46 Ti), the smallness of the excitation energy of the 7/2 − 1 (2 + 1 ) state suggests an underlying collective nature. Hence, a larger number of valence nucleons in the SM calculation or an α-clustered 44 Ti core would be recommendable.

The 12 C nucleus
Here, we present the theoretical predictions for the magnetic dipole moment of two energy levels of the 12 C nucleus within a macroscopic α-cluster model. For the latter, we select the Hamiltonian with the α-α interaction given by the isotropic Ali-Bodmer [138] potential in ref. [43] with the same parameters as in sec. II A of ref. [44]. On top of the latter, reproducing the short-range repulsive and long-range attractive effects of the strong force (cf. eqs. (2) and (3) of ref. [46]), we add the Coulomb interaction, accounting for the spherical charge distribution of the 4 He nucleus with charge radius R α = 1.44 fm (cf. eq. (4) of ref. [46]), as well as the Gaussian 3α force (cf. eq. (5) of ref. [46]), whose strength and range parameters were originally fitted to reproduce the binding energy of 12 C and the 2 + 1 -0 + 2 energy difference respectively (cf. refs. [139,43]). However, due to the adopted isotropic Ali-Bodmer potential with the parameters given in ref. [44], the energy of the ground state of this nucleus coincides with the opposite of the Hoyle state gap in the continuum and infinite-volume limit.
As in refs. [44,46], we transpose and apply the associated projectors in the numerical procedure for the extraction of the lattice eigenstates (cf. sec. 6 of ref. [46]). As in refs. [44,46], we introduce the lattice counterpart of SO(2), i.e. the cyclic group of order four, C 4 , associated to the counterclockwise rotation of π/2 about the z axis. Labeling the irreps of the former, I z , by positive integers modulo 4 [44], and identifying the irreps of the cubic group with Γ [46], we represent each lattice eigenstate with a round bracket (cf. ref. [45]) as |N , Γ, I z , S, P), where N is a positive nonzero integer that denotes the order with which the lattice eigenstates with the same transformation properties appear in the spectrum, thus fulfilling the role of the principal quantum number. When the SO(3) multiplet to which each state belongs in the continuum and infinite volume is identified, the angular momentum quantum number, , is added to the lattice states. Starting from states with well-defined I z (cf. eq. (13)) and part of the same SO (3) multiplet, it is possible to obtain states with good angular momentum projection along the z-axis in the continuum and infinite volume limit. The replacement of the C 4 by the SO (2) where the label Γ has been dropped for brevity.
Exactly the states in eq. (14) constitute the main ingredient for the computation of the average values of the magnetic dipole moment operator (cf. [110]), defined in the relative reference frame where g (α) L is given by eq. (4) and (L iM ) z is the z-component of the angular momentum operator on the lattice. Its expression, consistently with the formula of the squared total angular momentum derived in sec. 3.3 of ref. [46], is defined in the relative frame by where a † iM (n iM ) and a iM (n iM ) are the creation and annihilation operators of the α-particle i relative to the α-particle M at the site (n iM ) of the lattice, whereas the coefficients C (1,K) k , that scale as the inverse of the lattice spacing, are defined in refs. [46,140]. Considering the fact that the center-of-mass degrees of freedom have been dropped from the Hamiltonian, the l.h.s. of eq. (15) is equivalent to the magnetic dipole moment in the absolute reference frame, where the coordinate of the M th particle is extracted from the relative ones of the other particles, subject to the constraint that r CM = 0.
Furthermore, the computation of the magnetic dipole moment associated with the two energy levels, provides the opportunity to continue the analysis of the artifacts induced by the lattice environment in average values of spherical tensor operators that started in ref. [45]. In this respect, we here focus on the behaviour of the average values of the former observable with the lattice spacing and, aiming at reducing the discretization errors, we compute the isotropic averages (ref. [45]) of the magnetic dipole moments. Extending the definition in eq. (21) of ref. [45] to the magnetic dipole moment in the relative frame in eq. (15), the isotropically averaged expression of µ(N , , P) (r) is obtained, where the lattice angular momentum operator is expressed in the spherical basis and the Clebsch-Gordan coefficents of SO (3), ( |mm m ), have been introduced in the notation of ref. [141].
The counterpart of the last equation in the absolute reference frame can be obtained by applying the same definition (cf. ref. [45]) to the matrix element (N S P|μ Applying the theorem, we define the reduced matrix elements ofμ between lattice eigenstates, in the relative and in the absolute reference frames, ρ ≡ a or r. The definition in eq. (19) coincides with the one given in eqs. (15) and (19) of ref. [45], up to a factor equal to the Clebsch-Gordan coefficient ( 1 |000) multiplied by the constant 3/4π.
As in the latter study, due to the breaking of rotational symmetry, the SO(2) labels have not been dropped in the definition of the reduced brackets in eq. (18). Additionally, we assume that the Clebsch-Gordan coefficients on the r.h.s. of eq. (19) do not vanish, a fact that is guaranteed by the triangular inequality between the angular momenta, | − 1| ≤ ≤ + 1 and the conservation of the third component of the angular momentum m +m = m. Besides, the latter three projections are also required to be nonzero, since the sum of the three angular momenta in eq. (19), 2 + 1, is odd [141].
Reduced matrix elements can be immediately extended to the isotropic averages (cf. eq. (18)), with the same convention on the ρ index. Both the quantities in eqs. (19) and (20) are expected to reproduce asymptotically the exact reduced matrix elements in the continuum and infinite-volume limit, , (21) where N is the counterpart of N in the R 3M −3 configuration space.
Let us now begin with the behaviour of the magnetic dipole moment of the 2 + 1 state at 4.4398(2) MeV of 12 C as a function of the lattice spacing, displayed in fig. 4. As a result of the tuning of the parameters of the phenomenological 2α and 3α potentials of the present model [44,46], in the continuum and infinite-volume limit the energy eigenvalue of this state converges to ≈ 3.35 MeV, as shown by the curves associated with the E and T 2 multiplets of the cubic group and by the corresponding multiplet-averaged (cf. eq. (50) of ref. [46]) solid line in figs. 28 and 32 in ref. [46]. By setting the side of the cubic lattice to L ≥ 19 fm as in refs. [44,46] we have reduced the finite-volume errors associated with the energy eigenvalues to ≈ 10 −2 MeV, in order to remove the latter artifacts in the present analysis of discretization effects. Figure 4: Average value of the magnetic dipole moment for the 2 + 1 energy level of 12 C as a function of the lattice spacing. The broken line with circles denotes the estimated magnetic dipole moments obtained from the multiplet-averaged values of the squared total angular momentum operator, L 2 [46,140], over the five wavefunctions composing the 2 + 1 level, degenerate in the continuum and infinite volume limit. The broken line with triangles denotes the matrix elements of the magnetic dipole moment operator (cf. eqs. (15) or (17)) over the lattice states with maximum projection along the z axis. The isotropically-averaged estimates of the same observable (cf. eq. (18)), extracted from the matrix element ofμ with lattice wavefunctions with a different third component of the total angular momentum operator, are displayed by the broken lines with squares. The theoretical value of ≈ +1.0068 µ N is marked with a dotted line and reproduced with 10 −3 precision by the values of 22|μ 0 |22 and (22|μ 0 |22) • at a ≈ 0.65 fm, both equal to ≈ 1.0018 µ N . Finite-volume effects are suppressed by the constraint N a ≥ 19 fm.
In its generality, the expression of the magnetic dipole moment operator for α conjugate nuclei in eq. (3) permits to obtain average values forμ among states with maximum angular momentum projection, starting from estimates of the maximumL z eigenvalue. In the present case, the latter can be extracted from the average values of the squared total angular momentum operator as a function of the lattice spacing analyzed in fig. 39 of ref. [46]. Exploiting eq. (3) and the expression of the eigenvalues of L 2 (cf. eqs. (24)-(27) of ref. [46]) in the continuum and infinite volume, ( + 1) 2 , we have obtained the mulitplet-averaged estimates of the magnetic dipole moments for the 2 + 1 state, that we denote with (μ 0 ) L 2 . In fig. 4, the open circles at different values of the lattice spacing indeed describe a curve that lies always above the expectation value of the magnetic dipole moment at +1.0068 µ N , marked by a dotted line in fig. 4.
In particular, the minimum in (μ 0 ) L 2 at a ≈ 2.80 fm corresponds to a deep minimum of the average value of L 2 in the E multiplet of O, constituting the 2 + 1 energy state in the continuum and infinite-volume limit (cf. fig. (39) of ref. [46]). In contrast to energy eigenvalues and the average values of the α-α separation (cf. E r and R in figs. 32 and 33 of ref. [46]), the minima of the total squared angular momentum operator and hence (μ 0 ) L 2 cannot be unambigously mapped to the local maxima of the probability density function associated to the lattice eigenstates they refer to.
Concerning the isotropic average of the magnetic dipole moment (cf. eq. (18)), its behaviour with the lattice spacing shown by the curve with solid circles in fig. 4, follows quite faithfully the one of the magnetic dipole moment (open triangles in fig. 4) computed from the lattice 2 + 1 eigenstates with maximum angular momentum projection.
More precisely, (22|μ 0 |22) • slightly quenches the discretization artifacts in the very large lattice spacing region (a 2.75 fm) as well as in a small region around a ≈ 2.37 fm. However, for a 2.20 fm no significant improvement in the estimate of the magnetic moment from the isotropic average is observed, whereas in the peak region around (a ≈ 2.60 fm) the deviation from the asymptotic value of µ(1, 2, +) becomes even larger than for (22|μ 0 |22).
In summary, the isotropic average slightly reduces the discrepancies in the magnetic dipole moment between the extremal points. The reasons underlying this behaviour can be better understood by considering the individual reduced matrix elements (cf. eq. (19)) contributing to the reduced isotropic average in eq. (20). As it can be proven by rewriting the angular momentum operator in spherical components, the matrix elements contributing to the isotropic average in eq. (18) are overall 12 and fulfill the following symmetry relations: that vanishes as the associated Clebsch-Gordan coefficient, and Finally, from comparison between the magnetic dipole moment in the 2 + 1 of 12 C in fig. 5 and the quadrupole and hexadecupole electric moments in the 2 + 1 state of 8 Be in figs. 2 and 3 of ref. [45], we can conclude that the former observable is significantly less affected by discretization errors, in the same range of lattice spacing. More precisely, while all reduced matrix elements (2m ||μ m ||2m) at a ≈ 1.75 fm differ by less than 10% to the continuum and infinite volume counterpart, at the same lattice spacing some of the contributions to the isotropically-averaged electric qudrupole and hexadecupole moment differ by more than 50 % from their respective asymptotic values.
Second, we consider the magnetic dipole moment of the 3 − 1 state at 9.641(5) MeV of 12 C. In the adopted macroscopic α-cluster approach, the energy eigenvalue of the lattice counterpart of the 3 − 1 state coverges to ≈ 5.85 MeV above the g.s. in the continuum and infinite-volume limit (cf. fig. 36 of ref. [46]).  fig. 40 of ref. [46]).
Concerning the magnetic dipole moment computed from the state with maximum angular momentum projection along the z axis, (33|µ 0 |33), the curve with open triangles in fig. 6 reaches the theoretical value of µ(1, 3, −) in the continuum limit from below and displays a minimum at ≈ 2.29 fm, as the curve for the 2 + 1 level of 12 C (cf. fig. 4). However, the decrease of (33|µ 0 |33) in the region 1.75 a 2.25 fm is steeper than the one of (22|µ 0 |22) in the same interval (cf. fig. 4).
Nonetheless, the calculation of the isotropicaly-averaged value of (33|µ 0 |33) (cf. eq. (18)) yields a sizable suppression of the discretization errors in the region 1.80 a 2.80 fm, with a peak of ∼ 15% around the local minimum at a ≈ 2.35 fm. The curve for (33|µ 0 |33) • , indeed, lies above the one for (33|µ 0 |33) in the whole interval 0.65 a 3.10 fm (cf. fig. 6), unlike the 2 + 1 case (cf. fig. 4). The origin of this behaviour can be better investigated by considering the individual contributions (cf. eq. (19)) to the reduced isotropic average in eq. (20). Specifically, the non-vanishing matrix elements on the r.h.s. of eq. (18) are now 18 and fulfill the following symmetry identities: Due to the latter relations, the independent contributions to (33||μ 1 ||33) • reduce to six in total (cf. fig. 7). Differently from the 2 + 1 multiplet in fig. 5, the behaviour of these reduced brackets as a function of the lattice spacing is more multifaceted. In particular, the diagonal matrix element with maximum angular momentum projection (cf. the curve with diagonal crosses in fig. 7) follows a path that almost overlaps the one of the element (30||μ −1 ||31) (pentagons) throughout the considered interval of lattice spacing, displaying a local minimum at a ≈ 2.33 fm, followed by a local maximum at a ≈ 2.67 fm and a plateau in the large lattice-spacing region.
In the opposite direction, the two matrix elements undergo a steep increase with an inflection point at a ≈ 1.90 fm, that eventualy results in a plateau in the continuum limit. Both which indeed lie closest to the isotropic average (solid curve with full circles), except for intersection regions at 1.85 and 3.18 fm and in the stretch between the two local extrema. An inflection point is also detected around a ≈ 3.20 fm, as for the (32||μ 0 ||32) curve (cf. fig. 7).
Concerning the diagonal matrix element with projection 1 (vertical crosses), its curve displays a second local maximum at a ≈ 1.97 fm, followed by another minimum at a ≈ 1.88 fm, in contrast with all other reduced matrix elements. Additionally, its convergence to the continuum and infinite-volume value of (31||μ 0 ||31) is the slowest among the brackets of the magnetic dipole moment operator, since 1% precision is reached only at a ≈ 0.85 fm.
However, the two minima of the curve of the (31||μ 0 ||31) bracket admit an interpretation on the basis of the range parameters η −1 1 and η −1 0 of the Ali-Bodmer potential [43], equal to 1.89 and 2.29 fm respectively. The probability density function associated with the lattice state |31), in fact, possesses local maxima that can be included in the lattice points when a ≈ η −1 1 (cf. sec. 2 of ref. [46]).

Conclusion
With the present work, we have complemented the vast literature on magnetic dipole moments [50,71,109], by paralleling the average values of this observable for even-even self-conjugate nuclei from 12 C to 44 Ti with the partitioning of nuclear matter into α-clusters. The measured gyromagnetic factors of excited states of these nuclei turn out to assume, within their errors, the same value of g ≈ +0.50. It is shown that this specific feature can be explained on the basis of collective excitations of 4 He nuclei, neglecting the degrees of freedom of the single nucleons.
Such macroscopic approaches are best suited to low-lying excited states located at the vicinity of α-decay thresholds.
The stability of the α-clusters has been further investigated for isotopes of the above chain of nuclei by appending one or two neutrons. Relying on Schmidt estimates [39], we confirmed that, in proximity of the shell closures Z, N = 8 and 20 the extra one or two neutrons determine the experimental g-factor values. The original positive values of the corresponding self-conjugate nuclei undergo a significant reduction generated by the contribution of the negative g-factor of the additional neutrons. For the considered energy levels of the N = Z + 1 and 2 semi-magic nuclei in our chains, α-clustering is mitigated by the added neutrons. Outside the shell closures, for N = Z + 2 nuclei the added neutrons loose almost all their influence, and the predictions of eq. (4) are aligned with the experimental g-factors.
Spurred by the recent literature on the subject [44,45,46,25], we have calculated the magnetic dipole moment of the 2 + 1 state at 4.44 MeV and the 3 − 1 state at 9.64 MeV of 12 C in the framework of a macroscopic α-cluster model on the lattice in ref. [44]. Suppressing the finite-volume effects, we have analyzed the discretization artifacts through the behaviour of this observable as a function of the lattice spacing, ranging from a ≈ 0.65 to 3.50 fm.
Unlike the energy (cf. fig. 32 in ref. [46]) or the squared angular momentum (cf. fig. 39 in ref. [46] and the curve with open circles in fig. 4), the magnetic dipole moment appears rather insensitive to discretization effects for a 2.0 fm, where the artifacts amount to, at most, the 15% of the continuum and inifinite-volume value, governed by eq. (3).
In particular, for the 2 + 1 state, the isotropic average, introduced in ref. [45], turns out not to improve significantly the estimate of the magnetic dipole moment provided by matrix element with maximum angular momentum projection along the lattice z-axis alone (cf. the curve with solid circles in fig. 4). The latter conclusion does not hold in the 3 − 1 case, where the isotropic average turns out to improve sizably the esitmate of this observable in the region with a 3.05 fm (cf. fig. 6).
If its weak coupling to angular momentum symmetry breaking is confirmed in further nuclear systems on the lattice, the magnetic dipole moment evaluated at finite and sizable values of the lattice spacing (e.g. a ∼ 2.0 fm) can corroborate the classification of lattice eigenstates in terms of SO(3) quantum numbers, which is exact only in the continuum and infinite-volume limit.

Acknowledgements
We thank T. Otsuka for the shell-model calculation of the g-factors of the 2 + 1 states of 10 Be and 10 C. We gratefully acknowledge funding from the Deutsche Forschungsgemeinschaft   Applying the great orthogonality theorem, the decomposition of the 2 + 1-dimensional representations of the rotation group into the #Cl cubic group irreps, is obtained (cf. tab. 5 and tab. II of ref. [44]), where the multiplicity of the latter into the original irrep D is denoted by In the last equation, the order of O of the finite group is at the denominator, whereas χ ν i and χ i represent respectively the characters of the irreps of the cubic and the rotation group in the conjugacy class Cl i with |Cl i | elements. Moreover, the map between the basis states of SO (3, Z) and SO(3) irreps is provided by the projectors in eq. (40) of ref. [46].
As a consequence of this descent in symmetry, in the cubic lattice environment the maximum rank of any irreducible tensor operator runs from one to three. As shown in sec. 4 of ref. [46] for the energy eigenstates, the non-null components q of T (Γ,k) and |Γ , admixture of the q mod 4 components of their SO(3) counterparts, can be unambigously denoted with the I z label. The ensuing distribution of m components of a spin-irrep into the ( , Γ) irreps of SO(3, Z) is known as subduction [143].
Furthermore, when the multiplicity coefficient q Γ of the irrep Γ of O is larger than one, further linear transformations on the outcoming states (cf. eq. (26)) or cubic tensor components (cf. eq. (27)) should be performed, in order to block-diagonalize the relevant projector and disentangle the repeated multiplets of states [140].
Eventually, for the construction of the tables (cf. tab. ??-17), the basis eigenfunctions for the '0' and '2' irreps of the cyclic group of order four, C 4 , generated by an element of the conjugacy class 6C 4 (π/2) of the cubic group (e.g. the counterclockwise rotation of π/2 about the z-axis of the cube) are assumed to be real, whereas the following phase convention for the basis state Ψ Γ,Iz=1 (Ψ Γ,Iz=3 ) belonging to the '1' ('3') irrep of C 4 is understood.