On the use of quantality in nuclei and many-body systems

Abstract

The use of quantality is discussed in the case of nuclei and other many-body systems such as atomic electrons. This dimensionless quantity is known to indicate when a many-body system behaves like a crystal or a quantum liquid. Its role is further analyzed, showing the emergence of a fundamental lengthscale, the limit radius, which corresponds to the hard-core of the nucleon-nucleon interaction in the case of nucleons, and to a value close to the Bohr radius in the case of atomic electrons. The occurrence of a cluster phase in nuclei is analyzed using the quantality through its relation to the localization parameter, allowing for the identification of both the number of nucleons and the density as control parameters for the occurrence of this phase. The relation of the quantality to the magnitude of the interaction also exhibits a third dimensionless parameter, monitoring the magnitude of the spin-orbit effect in finite systems, through the realization of the pseudo-spin symmetry. The impact of quantality on the spin-orbit effect is compared in various many-body systems. The role of quantality in the relative effect of the binding energy and the shell one is also analyzed in nuclei. Finally, additional dimensionless quantities are proposed from the generalization of the quantality. Nuclei are found to be exceptional systems because all their dimensionless quantities are close to the order of unity, at variance with other many-body systems.

I Introduction

Quantality, a powerful dimensionless ratio first introduced by de Boer [1], was used by Mottelson to indicate when a system behaves like a quantum liquid (QL) rather than a crystal. It is larger for quantum liquid states than for crystal ones [2]. The former have constituents with delocalized wave-functions, leading to a homogeneous density, whereas the latter have localized constituents at given nodes. Superfluid Helium systems, nuclei, and atomic electrons are of QL type [2, 3, 4]. More precisely, nuclei and atomic electrons are found to be the two most QL-like systems, having larger values of the quantality than superfluid Helium systems.

An intermediate phase could appear between the crystal and the quantum liquid phase, the clusterized one [5, 6, 7]. The introduction of another dimensionless quantity, the localization parameter [8], could help, for instance, to understand why light nuclei are more clustered than heavy ones [9]. This parameter depends on the spatial dispersion of the nucleonic wavefunction which is determined by the equation of motion. Therefore, the localization parameter is richer than the quantality because it takes into account finite-size effects. It is analogous to the Wigner or Brückner parameter introduced in plasma physics or in condensed matter physics [10]. It allows for a unified understanding of cluster phases as transitional states between crystal and quantum liquid phases, also in nuclei [11, 12, 13, 6]. The relation between the quantality and the localization parameter was studied in [14].

Understanding $\alpha$ clusterization in nuclei as a different phase, compared to the QL one is relevant. The use of the localization parameter calculated both analytically with the Harmonic Oscillator (HO) approximation, and with covariant Energy Density Functionnals (EDF), allowed to understand nuclear clusterization as a nuclear transitional phase, compared to the usual nuclear QL one [8, 14]. A transition from a gas of $\alpha$ particles to nuclear liquid was also related to the properties of the nucleon-nucleon interaction with the lattice calculations [15, 9]. More recently, Effective Field Theories (EFTs) also tackle the occurrence of cluster phase in nuclei [16]. Recent ab initio PGCM calculations also confirm clusteristation traces in light nuclei, especially in ¹⁶O [17].

Control parameters allowing for the appearance of the cluster phase in nuclei have been discussed through these models, providing a broader understanding of the occurrence of clusterization than the only vicinity of the multi-alpha particle threshold provided by the seminal Ikeda conjecture [18]: the nucleon-nucleon interaction [15], the depth of the confining potential [8], the nucleon number [14, 13], the nuclear deformation [19, 13], the nuclear density [20, 21, 22, 23], or the temperature [24] could drive the transition from a nuclear quantum liquid state to a clusterized one. These studies allow to consider nuclear clusterization as an intermediate phase, as in other systems [5, 6, 7], which can occur in nuclear ground states. On the experimental side, recent signals using alpha knock-out reactions [25, 26] attracted a lot of attention as they allow to probe the degree of clusterization over the nuclear chart, from light nuclei (Be) to heavy ones such as the Tin isotopes.

Despite its powerfulness, the role of quantality has been very scarcely studied in nuclear physics [2, 3, 13, 4]. In this work, we further explore how quantality can be used to understand the main differences between atomic electrons and nucleons in their quantum liquid phases. It is then used to identify the control parameters for the cluster phase in nuclei, namely the role of the number of nucleons and the density. The relation to the low-density transition towards clusterization is analyzed. The use of quantality, as a dimensionless ratio, to understand its relation with the interaction coupling constants leads to a fundamental relation, allowing to study the magnitude of the spin-orbit effect in finite nuclei. The relative magnitude of the binding energy to the shell-structure one is also shown to depend on the quantality.

The present work aims to analyze the behavior of many-body systems with simple but powerful indicators, by progressively refining their quantum description. First, the role of the Zero Point Energy (ZPE) is discussed. Then, the inclusion of the dispersion of the wave-function of the constituents, from the solution of a e.g. Schrödinger equation is considered. As a third step, the relation with the spin-orbit effect is given, based on the non-relativistic reduction of the Dirac equation. It will be shown how quantality plays a pivotal role during all these three stages. In section II, the emergence of a fundamental lengthscale, the limit radius, is studied in relation to the ZPE and the quantality. The relation between the quantality and the localization parameter is discussed in section III, showing how the appearance of clusters in the low-density regime is recovered. Section IV relates the quantality to the dimensionless magnitude of the interaction, showing how the spin-orbit effect in finite systems can be deduced. Finally, a systematic study of possible dimensionless quantities is provided in Section V, generalizing the concept of quantality.

II The quantality, the zero point energy, and the limit radius

Let us consider a many-body system composed of constituents of mass m and interacting with an interaction of typical magnitude V₀. The inter-constituent distance at equilibrium is r₀. These 3 quantities characterize the main behavior of the system as they correspond to the input of its equation of motion. Instead of solving it as exactly as possible, a complementary way to predict the qualitative behavior of the system is to use indicators, incorporating step by step quantum mechanical effects. As a first step, the quantality is defined as the ratio of the ZPE T₀ to the magnitude of the potential V₀:

$$\Lambda\equiv\frac{2T_{0}}{V_{0}}=\frac{\hbar^{2}}{mr_{0}^{2}V_{0}}$$ (1)

Here, the ZPE T₀ is defined as the kinetic energy due to the equilibrium distance r₀, through the Heisenberg relation:

$$T_{0}=\frac{\hbar^{2}}{2mr_{0}^{2}}$$ (2)

Quantality was used by Mottelson to characterize QL vs crystal behavior of systems [2]. The typical value calculated with Eq. (1) is $\Lambda\simeq$ 10^-2,-3 in the case of crystals like atoms and molecules, and $\Lambda\simeq$ 0.1-1 (i.e. more than one order of magnitude larger) in the case of QL such as ⁴He, nuclei or electrons in atoms (see Table I). In the case of nuclei, $\Lambda\simeq$ 0.4, using typical values of the nucleon-nucleon interaction, namely V₀ $\simeq$ 100 MeV and r₀ $\simeq$ 1 fm [2]. The quantality indicates the importance to analyze the relation between the ZPE and the potential energy in finite systems. Moreover, among the quantality values of Table I, nuclei and atomic electrons have the largest one, indicating strong quantum effects. We shall discuss these 2 systems. It should be noted that i) quantality is valid for both fermions and bosons, ii) the quantum liquid nature of electron in atoms has been discussed in Ref. [2], and iii) results have to be considered within one order of magnitude, due to the semi-quantitative nature of the approach.

Constituent	m	r0 (nm)	V0 (eV)	T0 (eV)	$\alpha$	$\eta$	$\cal{A}$	$\Lambda$	State	ZPEv	QV	TI
20Ne atom	20	0.31	31 10-4	1.1 10-5	4.9 10-6	-6.1 1012	12.0	0.007	crystal	3.3 10-8	5.5 10-21	8.0 10-19
H2 molecule	2	0.33	32 10-4	9.5 10-5	5.3 10-6	-5.9 1011	4.1	0.06	crystal	3.2 10-7	5.3 10-19	9.0 10-18
4He atom	4	0.29	8.6 10-4	6.0 10-5	1.2 10-6	-4.4 1012	2.7	0.14	QL	1.9 10-7	4.7 10-20	2.7 10-19
3He atom	3	0.29	8.6 10-4	0.8 10-5	1.2 10-6	-3.2 1012	2.3	0.19	QL	2.6 10-7	7.7 10-20	3.8 10-19
Nucleon	1	10-6	100 106	25 106	0.51	1.3	1.6	0.41	QL	0.18	0.02	0.04
e- in atoms	5 10-4	0.05	10	16	2.5 10-3	-5 104	0.6	3.1	QL	8.1 10-3	1.6 10-7	5.0 10-8

Table 1: Constituent mass m, interaction lengthscale r₀, its magnitude V₀, zero point kinetic energy T₀, effective coupling constant $\alpha$, spin-orbit parameter $\eta$ [27], action $\cal{A}$ and quantality $\Lambda$ [2] for various many-body systems. m is given in units of a nucleon mass for convenience. The Zero Point Energy velocity (ZPEv), quadratic velocity (QV) and total interaction (TI) are introduced in section V.

The inter-constituent distance r₀ has a lower boundary r_l, so that r${}_{0}>$r_l : below this limit radius r_l, the ZPE is large enough to overcome the potential energy and make the system unbound. The corresponding limit condition T₀=V₀ leads to

$$r_{l}=\frac{\hbar}{\sqrt{2mV_{0}}}$$ (3)

Fig. 1 gives a view of the interplay of V₀, r₀, r_l and the ZPE considering a schematic interconstituent interaction, in the three following cases: r${}_{0}<$r_l, r${}_{0}=$r_l and r${}_{0}>$r_l, for which the system is unbound, at the boundary limit and bound, respectively.

Figure 1: Schematic energy position of a 2 constituents system in the case of a square well interconstituent interaction for a) r${}_{0}<$r_l, b) r${}_{0}=$r_l and c) r${}_{0}>$r_l.

The quantality (1) can then be written as:

$$\Lambda=2\left(\frac{r_{l}}{r_{0}}\right)^{2}$$ (4)

Close to the limit case, such as for electrons in atoms, r₀=r_l and the quantality takes its maximal value, namely $\Lambda\simeq$2. This is confirmed by the phenomenological value of quantality [2], reproduced in Table I. In other systems, $\Lambda\leq$2. For instance, in the case of nuclei, $\Lambda\simeq$1. This is a specific feature of nuclei: $\Lambda$ is of the order of unity, as well as their other dimensionless quantities. We define it as specific naturalness and its origin shall be discussed in section IV.1. For instance, although both atomic electrons and nuclei have a quantality value close to unity, it will not be the case for other dimensionless quantities in the case of atomic electrons, contrarily to nuclei, as discussed in section IV. Eq (1) also shows that $\Lambda$ monitors the delocalization effects with respect to the minimal radius r_l: the larger the equilibrium distance r₀, the smaller the quantality with respect to the maximal value $\Lambda\simeq$2.

Whether the particle is at equilibrium at a distance r₀ equal to or larger than r_l may depend on several causes, such as delocalization effects. In nuclei, using the typical magnitude V₀ of the nucleon-nucleon interaction gives r_l=0.5 fm from Eq. (3). This corresponds to the size of the repulsive core of the interaction. It can be understood with the specific naturalness of nuclear systems: as an order of magnitude reasoning, V₀r${}_{l}\simeq\hbar$c.

Eq (4) can be rewritten as

$$r_{0}=r_{l}\sqrt{\frac{2}{\Lambda}}$$ (5)

The equilibrium distance r₀ remains larger than r_l because $\Lambda\leq$2. However, the smaller $\Lambda$, the larger interconstituent distance r₀, and the smaller the delocalization effect, because the constituents are further apart. Therefore, $\Lambda$ drives localization effect. However, it misses finite-size effects. This will be taken into account in section III.

In nuclei the limit radius, being of the order of the hard-core one, gives a straightforward reason why the hard-core is washed out: it corresponds to the minimal radius allowed by the ZPE, and in practice, the mean inter-nucleon distance is r${}_{0}\simeq$ 1.2 fm, which is significantly larger. Therefore, due to a delocalization effect monitored by the quantality, the saturation density is lower than the packing one, which would correspond to a nucleonic system with interconstituent distance r_l, the size of the repulsive core of the interaction. As discussed by Mottelson, this is an important question to understand [28, 29]. More precisely, using the relation between the interconstituent distance and the density:

$$r_{0,l}=\left(\frac{3}{4\pi\rho_{0,l}}\right)^{\frac{1}{3}},$$ (6)

Eq (5) leads to

$$\rho_{l}=\left(\frac{2}{\Lambda}\right)^{3/2}\rho_{0}\simeq 10\rho_{0}$$ (7)

where the quantality value $\Lambda$$\simeq$0.4 for nuclei has been used. Therefore, the packing density in nuclei is about one order magnitude larger than the saturation one ($\rho_{0}\simeq$0.16 fm^-3), which can be understood from the value of quantality: $\rho_{l}\simeq$1.5 fm^-3.

Turning to the case of electrons in atoms, and using the Coulomb energy for V₀:

$$V_{0}=\frac{e^{2}}{a_{B}},$$ (8)

where a${}_{B}\equiv\hbar^{2}$/me² is the Bohr radius, leads to r_l=a_B/$\sqrt{2}$, using Eq. (3). Since the equilibrium distance is r₀=a_B, this means that r₀ is close to its minimal value and that the ZPE is close to its maximal value: r₀=a_B=r_l$\sqrt{2}$. The present approach allows to realize that atomic electrons are an exceptional system because it is close to a limit case. Its origin can be found in the specific 1/r behavior of the Coulomb potential, interplaying with the 1/r² behavior of the ZPE: equilibrium position involves spatial derivatives which are related, as the derivative of the potential energy can be expressed from the kinetic energy [30].

Therefore, the limit radius r_l is a fundamental quantity, recovering both the hard-core size in nuclei as well as a value close to the Bohr radius in the case of atomic electrons.

III Quantality and the cluster phase in nuclei

A step further in including quantum effects in quantities monitoring the behavior of the system can be achieved by the use of the localization parameter [8, 14]. It is defined as

$$\alpha_{loc}(r_{0})\equiv\frac{\lambda}{r_{0}}$$ (9)

where $\lambda$ is the typical spatial extension of the constituent wavefunction and r₀ the interconstituent distance at equilibrium. In order to study cluster phases, this dimensionless quantity is $>$1 for homogeneous quantum liquid phase, $<$1 in the crystal case and $\sim$1 in the transitional cluster case. This quantity is related to the so-called Brückner one in condensed matter (see Appendix A). The interest of the localization parameter is that it takes into account finite-size effects, compared to quantality, by introducing $\lambda$ as an additional quantity to consider in addition to m, r₀ and V₀.

In the case of nuclei, the spatial extension of the nucleonic wave function $\lambda$ can be calculated microscopically from e.g. EDF approaches. It is well approximated by the confining length of a harmonic oscillator (HO) potential, which parameters are chosen to match the radius $R$ of a given nucleus [11]:

$$\lambda\simeq\sqrt{\frac{\hbar}{m\omega}}=\frac{\sqrt{\hbar R}}{(2mV_{0})^{\frac{1}{4}}}$$ (10)

with $\hbar\omega$ the typical energy of the HO, $m$ the nucleon mass and $V_{0}$ the depth of the confining potential. The typical confining potential V₀ has been taken as the one of the magnitude of the nucleon-nucleon interaction of Eq. (1). This is well justified in nuclei, as the empirical values are 80 MeV [8] and 100 MeV [2] for the depth of the confining potential, and the magnitude of the nucleon-nucleon interaction, respectively. A more detailed justification for this approximation is given in the appendix of Ref. [19].

Using Eqs (3) and (6), Eq. (10) can be rewritten as

$$\lambda\simeq\sqrt{r_{l}R}=\sqrt{r_{0}R}\left(\frac{\Lambda}{2}\right)^{1/4}$$ (11)

Eq (11) shows that the typical nucleonic dispersion can be interpreted as the geometrical average of the limit radius from the ZPE and the nuclei one. It also emphasizes the contribution of both the quantality and finite-size effects on the dispersion. It should be noted that this equation can also be applied to atomic electrons, the HO confining potential being a first-order approximation to any bound system. Since both r_l and R are close to the Bohr radius, as discussed in the previous section, the dispersion of the electron wave function $\lambda$ is also close to the Bohr radius.

The localization parameter can be related to the quantality. Especially in nuclei, where R=r₀A^1/3, and using Eqs. (9),(11):

$$\alpha_{loc}\simeq\left(\frac{\Lambda}{2}\right)^{1/4}A^{1/6}\simeq 0.7A^{1/6},$$ (12)

showing the A dependence introduced by the localization parameter, compared to the quantality. This expression has been studied in previous works [13, 11], and allows to describe the more frequent cluster occurrence in light nuclei, compared to heavy ones. A value smaller or larger than 1 corresponds to localized or delocalized states in quantal systems, respectively. It should be noted that Eq (12) has also been confirmed by the use of microscopic calculations of the spatial dispersion of the wave function [14, 11].

Having identified the number of nucleons as a control parameter for the occurrence of a cluster phase, the role of the density in such phases can now be analyzed with the present approach. The inter-constituent distance is related to the density of the system. Let us consider a many-body system where the inter-constituent distance can change compared to its equilibrium value, due for instance, to external conditions such as compression or dilatation of the system, as it can occur in heavy-ion collisions. It can also be that the density of the system itself can be spatially dependent and hence have a lower density on the surface. For instance, finite nuclei display large fluctuations around $\rho_{\text{sat}}$, such that nucleons evolving in the low-density parts of finite nuclei could self-arrange into alpha-clusters [20, 21, 22, 23], as probed by recent (p,p$\alpha$) reactions [25].

As an in-medium monitor for cluster phase occurrence where finite-size effect can be important, such as in nuclei or for atomic electrons, the localization parameter (9) can be generalized to any interconstituent distance $\bar{r}$:

$$\alpha_{loc}=\frac{\lambda}{\bar{r}}$$ (13)

The mean interparticle-distance is related to the density of the system by

$$\bar{r}=\left(\frac{3}{4\pi\rho}\right)^{\frac{1}{3}}.$$ (14)

In the case of nuclei, the equilibrium distance $\bar{r}$=r₀ corresponds to the saturation density $\rho_{0}$.

Using R=$\bar{r}$A^1/3 owing to the short range of the nucleon-nucleon interaction and Eqs (10), (13), (14) lead to:

$$\alpha_{loc}=\sqrt{\frac{r_{l}}{\bar{r}}}A^{1/6}=\left(\frac{\rho}{\rho_{l}}A\right)^{1/6}$$ (15)

where $\rho_{l}$ is the density corresponding to the minimal radius. In nuclei, this is the packing density [28, 29] because r_l corresponds to the hard-core size, as discussed in section II.

Clusters shall occur when the nucleon spatial extension is of the order the inter-nucleon distance, namely $\alpha_{loc}\lesssim$1 [8, 14, 13]. In this case, Eq (15) brings a condition for cluster states:

$$\rho A\lesssim\rho_{l}=\left(\frac{2}{\Lambda}\right)^{3/2}\rho_{0}\simeq 10\rho_{0},$$ (16)

using Eq. (7). Hence, cluster states can occur in light nuclei (A$\simeq$10) at saturation density or at a lower density than the saturation one in heavy nuclei (A$>$10): either the density or the number of nucleons can be a control parameter for the occurrence of cluster phase from a nuclear quantum liquid one. Interestingly, the threshold number of A$\sim$ 10 for clusterization is driven by the quantality. Eq. (16) also shows that in ¹²⁰Sn, clusterization could occur below a density of about one-tenth of the saturation density, that is on its very surface. These results are in agreement with nuclear phenomenology, where lighter nuclei and/or smaller density (the so-called Mott transition [31, 32, 33, 34]) trigger cluster states. This is also in qualitative agreement with predictions of more microscopic approaches based e.g. on EDF, with or without explicit inclusion of alpha degrees of freedom [22, 21, 23], illustrating the relevance of the present discussion.

Another way to trigger the cluster phase is to dilute the whole nucleus [22, 23]. In the case of the transition from a quantum liquid to an alpha clusterised state by a decrease of the density, a more precise value of the density below which all nucleons are likely to combine into alpha-particles and condensate can be given. Because of Pauli blocking effects, alpha-particles in the nuclear medium shall dissolve as they start overlapping with each others. The critical Mott density at which this happens turns out to be given by the dimensionless parameter $\alpha_{loc}$ where the ratio is now between the size of an alpha-particle $R_{\alpha}\simeq r_{0}A_{\alpha}^{\frac{1}{3}}$ and the mean internucleon distance $\bar{r}$. The Mott density is then solution of $\alpha_{loc}=1$, that is, using Eq. (14)

$$\frac{\rho^{\alpha}_{\text{Mott}}}{\rho_{0}}=\frac{1}{A_{\alpha}}\sim 0.25$$ (17)

Therefore, in the case of a global dilution of the nucleus, below about one fourth of the saturation density, quartetting correlations are no more suppressed by Pauli blocking effects and a full alpha-clustering should correspond to the favoured arrangement of alpha-conjugate nuclei. This is in agreement with both covariant [23] and Gogny EDF [22] calculations in dilute nuclei.

IV Quantality for spin-orbit effect and shell structure

IV.1 The key formula

Dimensionless quantities are known to be a powerful tool in physics [35]. In addition to quantality, the dimensionless coupling constant is another important one:

$$\alpha\equiv\frac{r_{0}V_{0}}{\hbar c}$$ (18)

Eqs (1) and (18) yields the following relation between the quantality and the coupling constant:

$$\eta^{\prime}\Lambda\alpha^{2}=1$$ (19)

where

$$\eta^{\prime}\equiv\frac{mc^{2}}{V_{0}}$$ (20)

It should be noted that $\eta$’ is related to velocity effects in the system. For instance, using the T$\simeq$V₀ approximation for liquids, where T is the typical kinetic energy of the system, leads to, in the non-relativistic case

$$\eta^{\prime}\simeq 2\left(\frac{c}{v}\right)^{2}$$ (21)

where $v$ is the velocity of the constituents. Therefore, Eq (19) highlights the relation between the quantality, the coupling constant, and the velocity of the system.

The spin-orbit effect is well-known for being related to the velocity of the constituents, in a quantum framework. Hence, Eq (19) could be related to the spin-orbit effect. In a covariant approach based on a non-relativistic reduction of the Dirac equation, it has been shown the spin-orbit effect in many-body systems is driven by the quantity [27]

$$\eta\equiv\frac{mc^{2}}{V-S}$$ (22)

where V and S are the vector and scalar mean-field potentials generated by mediators of the interaction. A large value of $\eta$ (in absolute value) gives a small spin-orbit effect compared to shell effects, whereas a value close to 1 gives a similar magnitude between spin-orbit and shell effects. The sign of $\eta$ gives the energy ordering of the large j and small j states, whose degeneracy has been raised by the spin-orbit effect [27]: $\eta<$0 corresponds to the state with larger j value at larger energy. Typical values of $\eta$ are given in Table I.

In the case of atomic electrons, S=0, and there is only the attractive vector potential V. For instance, V=-13.6 eV for an electron in the Hydrogen atom. In the case of nuclei, V$\simeq$320 MeV is repulsive, whereas S$\simeq$-400 MeV is attractive. In all the cases, the magnitude of the central potential is V₀=-(V+S).

Eqs (20) and (22) leads to

$$\eta=\varepsilon_{VS}\eta^{\prime}$$ (23)

with

$$\varepsilon_{VS}\equiv\frac{-(V+S)}{V-S}=\frac{V_{0}}{V-S}$$ (24)

When there is no scalar mediator, such as in the vast majority of many-body systems, then $\varepsilon_{VS}$=-1. A departure from this value indicates the presence of a scalar mediator in the interaction. For instance, in the case of nuclei, using the typical values of V and S leads to $\varepsilon_{VS}$=1/9$\simeq$0.1. $\varepsilon_{VS}$ is hereafter called the Pseudo Spin Symmetry (PSS) breaking coefficient because the $\varepsilon_{VS}$=0 case corresponds to the realization of the PSS [36], namely V=-S. Moreover, the vector and the scalar potentials cannot be of the same sign, because of the repulsive/attractive nature of the spin triplet/singlet dependence of the interactions (check), implying $\mid\varepsilon_{VS}\mid\leq$1. Therefore, the case $\mid\varepsilon_{VS}\mid$=1 corresponds to the maximal PSS violating case.

Using Eqs (19) and (23) leads to :

$$\eta\Lambda\alpha^{2}=\varepsilon_{VS}$$ (25)

The relations (19) and (25) have a strong physics meaning because they relate quantities involved in various aspects of many-body systems. The interpretation of the coupling constant is obvious. Quantality drives localization aspects (crystal, cluster and QL phases), as discussed above. Finally, $\eta$ can be interpreted, in finite systems such as nuclei or atomic electrons, as the parameter driving the spin-orbit effect [27].

Figure 2 provides a more detailed overview. For a given fixed central potential value V₀=-(V+S), letting V-S vary gives a line in the ($\eta$,$\varepsilon_{VS}$) plane, whose slope is $\Lambda\alpha^{2}$, using Eqs (23),(19). In the case of most systems, like electrons in atoms or graphene, there is no scalar mediator in the interaction (S=0). Hence, Eq. (24) gives $\mid\epsilon_{VS}\mid$=1. In the case of nuclei, there are both vector and scalar mediators, implying $\mid\epsilon_{VS}\mid<$1. Moreover, the value of $\eta$ is close to 1 (see Table I).

Considering the quantality, QL systems have larger slopes than crystal ones. Concerning the interaction effect, a system at exact pseudo-spin symmetry has V=-S, i.e. V₀=$\alpha$=0, and is on the X axis. On the other hand, graphene has an infinite value of $\Lambda$ (see Appendix B), and is on the Y axis.

Figure 2: Sketch of the relation between the PSS coefficient $\varepsilon_{VS}$ and the spin-orbit parameter $\eta$ in the case of 3 bound systems: nuclei, electrons in atoms, and graphene, for a fixed value of V+S$<$0.

Let us discuss the above-mentioned dimensionless quantities on nuclei and atomic electrons. Nuclear states are known to usually behave as quantum liquid states and therefore $\Lambda\lesssim$ 1, as displayed on Table I. Moreover, the low energy QCD interaction yields $\alpha\sim$ 1. Eq. (19) then provides $\eta\sim$ 1. This value of $\eta$ allows for strong spin-orbit effect in finite systems [27]. The nucleus is therefore a very specific system where $\eta\sim\alpha\sim\Lambda\sim$1 within about one order of magnitude. This is due to naturalness.

In order to understand the origin of this specificity of naturalness in nuclei, the coupling constant Eq. (18) can be rewritten with the mass of the mediator m₀c${}^{2}\equiv\hbar$c/$r_{0}$:

$$\alpha=\frac{1}{\Lambda}\frac{m_{0}}{m}$$ (26)

Therefore, naturalness could be interpreted by the close mass of the nucleon and the mesons, namely the 1 GeV scale, themselves related to the scale of chiral symmetry breaking and to the chiral condensate [37, 38]: in nuclei, the mass of the mediator is of the same order of magnitude than the constituent particle mass (m_p/m_π $\simeq$7 at most). The m${}_{\pi,\sigma}\lesssim$m fact, responsible for $\alpha\sim$ 1, is related to QCD effects which provide a typical energy scale of few hundreds of MeV. This allows for similar masses between the lightest baryons and mesons. Eq. (26) exhibits the relationship between QCD driven masses ( $\alpha\sim$ 1), the strong spin-orbit effect in nuclei ($\eta\sim$ 1) and their QL behavior ($\Lambda\sim$ 1).

In the case of electrons in atoms, $\Lambda$ is close to 1 as well (even being close to its maximal value, see section II), which also behave as a QL [2]. However, this is not a specific naturalness system since the coupling constant is the electromagnetic one: $\alpha$=1/137, yielding $\eta\sim$10⁴ using Eq.(19). Namely, values are far from unity for $\alpha$ and $\eta$. Indeed, the definition of the spin-orbit parameter (Eq. (22)) yields $\eta\sim$m_ec²/V${}_{0}\sim$ (10⁵ eV/10 eV) $\sim$10⁴, which is in agreement with the experimentally observed magnitude of the spin-orbit effect, compared to the shell one [27], i.e. the fine structure: 1/$\eta\sim$10^-4.

IV.2 Spin-orbit, quantality, and coupling constant analysis of many-body systems

Let us further investigate the relation (25) between the quantality, the interaction, and the spin-orbit effect, and show how quantality impacts the spin-orbit effect in finite systems. The spin-orbit rule [27] is derived from the non-relativistic reduction of the Dirac equation and, therefore, encompasses information on the equation of motion. It gives the magnitude $x$ of the ratio of the HO shell gap to the one induced by the spin-orbit effect, as a function of $\eta$ [27]. 1/$x$ monitors the importance of the spin-orbit effect with respect to the main shell effect: 1/$x$ $\sim$ 1 in nuclei, in agreement with nuclear structure phenomenology showing a spin-orbit degeneracy raising able to build magic gaps. On the contrary, 1/$x$ $\sim$ 10^-4 for electrons in atoms, in agreement with the hyperfine structure in atoms.

Using the key relation (25) into the spin-orbit rule [27], yields:

$$x(\Lambda)=\left|\frac{\epsilon_{VS}}{\Lambda\alpha^{2}}-1+\frac{\Lambda\alpha^{2}}{4\epsilon_{VS}}\right|$$ (27)

where $x$ is the ratio of the HO shell gap to the spin-orbit one.

Fig. 3 displays the $x$($\Lambda$) relation, for both the QL and crystal cases. In the specific case of a quantum liquid ($\Lambda\simeq$1) and in the case of the electromagnetic interaction ($\alpha\simeq$10^-2,$\epsilon_{VS}$=-1), Eq. (27) reduces to

$$x(\Lambda=1)\simeq\frac{1}{\alpha^{2}}\simeq 10^{4}$$ (28)

recovering the specific fact that the coupling constant is driving the magnitude of the spin-orbit effect of electron in atoms [39]. Eq. (27) therefore provides a general understanding of the impact of the coupling constant on the spin-orbit effect, in various systems. Figure 3 shows that in a QL ($\Lambda\simeq$1), the large spin-orbit coupling ($x\lesssim$1) can only be reached in strongly interacting systems whereas the small spin-orbit coupling ($x>>$1) can only be reached in electromagnetic interacting systems: a large spin-orbit effect requires both a coupling constant of the order of unity (such as in the strong interaction case) and a quantum liquid behavior, as seen on Eq. (27).

In the case of a crystal ($\Lambda\lesssim$ 1), the spin-orbit effect is reduced, by a factor corresponding to the drop of the quantality value from a QL to a crystal, typically 2-3 orders of magnitude (Table I). This result is in agreement with dedicated calculations of the spin-orbit effect in crystals [40].

Figure 3: The relative magnitude of the HO to spin-orbit gaps in a many-body system as a function of the quantality (Eq. (27)), for various dimensionless coupling constant: the electromagnetic (EM) case (the fine structure constant), the low-energy strong interaction (S) between nucleons, and the case of graphene ($\alpha$=2.5).

Another interesting case is graphene (see appendix B), where $\eta$=0, $\Lambda$ goes to infinity, and $\alpha$ remains finite at about 2.5.

IV.3 Binding energy vs shell effects dependence on quantality

In Fermionic many-body systems, it could be relevant to analyze the respective value of the binding energy to the typical energy scale at work, namely the major shell gap. A constituent-dominated system means that the binding energy per particle is larger than the typical shell gap, and the Fermi energy of the system doesn’t change drastically. On the other hand, a shell-structure-dominated system means that the Fermi energy can be small compared to the shell gap. Let us show how quantality also monitors this relative effect.

In nuclei, the binding energy per nucleon is B/A $\sim$ 8 MeV, a bit larger than the typical magic shell gap which is about 5 MeV between the last level of a shell and the first one of the next shell. As a general question, the competition in finite systems between the energy per constituent and the shell effect shall be considered. This feature explains why the addition of one nucleon to the system has a larger effect than magicity effects, as seen on the plots of the evolution of the separation energy as a function of the neutron or proton number (see e.g. Fig. 1 of Ref. [41], where the lines usually do not cross each other).

Therefore, the relevant quantity to analyse this behavior is the y ratio defined as

$$y\equiv\frac{B/A}{\hbar\omega}\simeq\frac{\epsilon_{F}}{\hbar\omega}$$ (29)

where $\hbar\omega$ is the typical HO shell gap. B/A is typically the Fermi energy of the system: Fig. 4 displays the quantities to be used with respect to the confining potential of the system. y $>$ 1 is a constituent-dominated system where the addition of one constituent can drastically impact the system: this is the case for nuclei, as known from its phenomenology. y $<$ 1 is a shell structure dominated system.

Figure 4: The confining potential and relevant quantities

The binding energy per nucleon can be approximated by B/A $\sim\epsilon_{F}=V_{0}-E_{F}$ (Fig. 4). E_F shall be evaluated from a simple Fermi gas approximation which reads [42]:

$$E_{F}=\frac{\left(9\pi\right)^{2/3}}{8}\frac{\hbar^{2}}{mr_{0}^{2}}(1+\delta)^{2/3}$$ (30)

where $\delta$=(N-Z)/A is the neutron excess parameter. Over the nuclear chart one safely has 0$\lesssim\delta\lesssim$0.5.

The typical harmonic oscillator gap is [28]

$$\hbar\omega_{0}={\hbar\over R}\sqrt{2V_{0}\over m},$$ (31)

In both cases, approximated ways of solving the Schrödinger equation have been used (Fermi gas and HO).

For the numerical evaluation, two points should be considered: i) the Fermi gas model provides a quantitative value of B/A $\sim$ 8 MeV when V₀ $\simeq$ 45 MeV. This rather small value of V₀ is due to the approximations of the Fermi gas model. ii) injecting this value in (31) gives $\hbar\omega\simeq$10 MeV. The energy gap between the last subshell of a given major shell and the first one of the next shell (magic gap) is therefore rather of the order of $\hbar\omega$/2$\simeq$ 5 MeV.

This yields for the y ratio discussed above, using Eqs (30) and (31):

$$y\simeq\frac{V_{0}-E_{F}}{\hbar\omega_{0}/2}=\frac{\sqrt{2}R}{\hbar}\sqrt{mV_{0}}-\frac{\left[9\pi(1+\delta)\right]^{2/3}}{4}\frac{\hbar R}{r_{0}^{2}\sqrt{2mV_{0}}}$$ (32)

which can be rewritten, with the quantality (1):

$$y(\Lambda)\simeq A^{1/3}\left(\frac{2}{\sqrt{\Lambda}}-\frac{\left[9\pi(1+\delta)\right]^{2/3}}{4}\sqrt{\Lambda}\right)$$ (33)

using R=r₀A^1/3. Considering the quantality value for nucleons ($\Lambda\simeq$0.4, Table I) and the possible range of $\delta$ over the nuclear chart gives

$$y\simeq(1.2;1.7)A^{1/3}$$ (34)

Therefore, in nuclei y$>1$, meaning that the addition of a nucleon is larger on the binding energy than the main shell gap.

V Generalization of the quantality

The present study shows the pivotal role of the $\alpha$ (18), $\eta$ (22) and $\Lambda$ (1) ratios to characterize several fundamental properties of many-body states. More generally, it is possible to build in a systematic way dimensionless quantities from the 3 basic quantities of the system: V₀, r₀ and m. New dimensionless quantities could be produced in the spirit of the quantality using the Pi theorem [43]. In the following, in order to have only energy units, and for the sake of simplicity, the r₀ lengthscale can equivalently be described by the following mass:

$$m_{0}c^{2}=\frac{\hbar c}{r_{0}}$$ (35)

Setting the V₀, m₀, and m values imposes two independent ratios built from these. We choose the one driving the spin-orbit parameter (20):

$$\eta^{\prime}\equiv\frac{mc^{2}}{V_{0}}$$ (36)

and the dimensionless coupling constant:

$$\alpha=\frac{V_{0}}{m_{0}c^{2}}$$ (37)

In the case of up to maximum second order for any of the 3 basic quantities of the system (V₀, m₀ and m), four independent dimensionless quantities can be built. Table 1 summarizes the various dimensionless quantities discussed below for several many-body systems.

The first quantity

$$\Lambda\equiv\frac{1}{\eta^{\prime}\alpha^{2}}=\frac{m_{0}^{2}c^{4}}{mc^{2}V_{0}}$$ (38)

(38) is the quantality $\Lambda$ and the key formula (see Eq. (19)). It is related to the action of the system normalized to $\hbar$:

$${\cal A}\equiv\frac{r_{0}\sqrt{mV_{0}}}{\hbar}$$ (39)

by

$$\frac{1}{{\cal A}^{2}}=\Lambda$$ (40)

This shows that the action and the quantality of a system provide the same information. Indeed, it is well known that quantum effects in a system are large when its action is close to $\hbar$ [30]. This corresponds to $\cal{A}$$\gtrsim$ 1 and $\Lambda$ $\lesssim$ 1. This is the quantum liquid case. When quantum effects are smaller, such as in the crystal case, the action of the system is significantly larger than $\hbar$: $\cal{A}$$\gg$ 1 and therefore $\Lambda$ $\ll$ 1.

The second quantity which can be built from $\eta^{\prime}$ and $\alpha$ is

$$ZPEv\equiv\frac{1}{\eta^{\prime}\alpha}=\frac{m_{0}}{m}=\sqrt{\frac{2T_{0}}{mc^{2}}}\equiv\frac{v_{0}}{c}\equiv\beta_{0}$$ (41)

It corresponds to the velocity v₀ of the constituents, due to the ZPE. Its values in various systems are given in Table I. In the case of nuclei, it has a non-negligible value compared to one. Eq (41) shows that it comes from the naturalness of this specific system. Following this interpretation, the velocity of the nucleons is non-negligible compared to the light speed (i.e. at the upper limit of the non-relativistic approximation) because the spin-orbit effect is important ($\eta\simeq$1), as well as the magnitude of the interaction ($\alpha\simeq$1).

As the third quantity:

$$TI\equiv\frac{\alpha}{\eta^{\prime}}=\frac{V^{2}_{0}}{m_{0}c^{2}mc^{2}}$$ (42)

This quantity is large when the interaction is large and/or the spin-orbit effect is large (small value of $\eta$). Therefore, it monitors the magnitude of the total interaction, namely the central+spin-orbit effects. Nuclei are expected to have the largest value of the total interaction, because of the both large spin-orbit and interaction magnitudes.

The last independent quantity is

$$QV\equiv\frac{1}{\eta^{\prime 2}\alpha}=\frac{m_{0}c^{2}V_{0}}{m^{2}c^{4}}=\frac{\beta_{0}}{\eta^{\prime}}$$ (43)

It is, therefore, the ratio of the zero point motion velocity to the spin-orbit parameter. A large value corresponds to a large zero-point motion velocity and a large spin-orbit effect. This could be interpreted as a quadratic velocity quantity, namely the velocity of the system due to both zero-point motion and spin-orbit effect: as an approximation, let us consider the kinetic energy of the constituent of the order of V₀, as it can be in liquids. In this case, $\eta$’$\simeq$2$\beta$, where $\beta$ is the velocity of the constituent over c. Hence,

$$QV\simeq 2\beta_{0}\beta$$ (44)

The typical values of the quadratic velocity for various systems are given in Table I. Once again, in nuclei, this value is closer to one and the largest.

VI Conclusions

The quantality is a powerful tool allowing to identify various phases in many-body systems. Its upper value corresponds to a minimal radius, which is close to the Bohr radius in the case of atomic electrons, and leads, in the case of nuclei, to the size of the hard-core of the nucleon-nucleon. To analyze cluster phases in nuclei, the localization parameter can be related to the quantality. The more general definition of the localization parameter allows to show that both the nucleon number and the low density are control parameters for the occurrence of the cluster phase.

When related to the dimensionless coupling constant, the quantality drives the spin-orbit effect in finite systems. The nucleus is a very specific system, due to its naturalness (i.e. its dimensionless constants of the order of the unity). The nuclear clusterization condition could even be considered as a supernaturalness state: all the dimensionless parameters are close to one, including the localization parameter. The impact of quantality on the spin-orbit effect is exhibited, using the spin-orbit rule obtained from a non-relativistic reduction of the Dirac equation. Also, the quantality is shown to drive the binding energy per particle compared to the shell gap.

Three additional dimensionless quantities have been discussed. It should be reminded that the quantities used in the present work, contain information from the solved Schrödinger equation (of course, in an approximated way), through e.g. the evaluation of the dispersion of the nucleonic wave function obtained by microscopic EDF or HO approximation, or the behavior of the spin-orbit effect through a non-relativistic reduction of the Dirac equation. The present results show how quantality can be useful to nuclear structure, beyond a mere evaluation of its value in various many-body systems.

Appendix A The Brückner parameter

The Brückner parameter, used in condensed matter, can be related to the localization parameter (13), in the case of nuclei. Behaving quantum-mechanically, the effective strength of the interaction between nucleons, i.e. the extent to which interactions impact the properties of nuclei and make the latter deviates from the ideal free case, can be measured by a dimensionless ratio between the mean potential $\braket{V}$ and the mean quantal kinetic $\braket{K}$ energies, or equivalently the ratio between the mean interparticle distance $\bar{r}$ and the Bohr radius $a_{B}$. This is the so-called Brückner (also known as the Wigner) parameter (see e.g. [44]) $r_{s}=\frac{\braket{V}}{\braket{K}}\sim\frac{\bar{r}}{a_{B}}$. In the nuclear case, we have

$$r_{s}=\frac{\hbar\omega}{E_{0}},$$ (45)

where the energy of the harmonic oscillator reads $\hbar\omega=\frac{\hbar}{R}\sqrt{\frac{2V_{0}}{m}}$ and $E_{0}=\frac{\hbar^{2}}{2m}\left(\frac{3\pi^{2}\rho}{2}\right)^{\frac{2}{3}}$ stands for the Fermi energy of nuclear matter at density $\rho$ [42]. One gets from Eqs. (15) and (45):

$$r_{s}\simeq\alpha_{loc}^{-2}$$ (46)

showing that the localisation parameter captures the relevant effects.

It should be noted that in Plasma physics, these dimensionless parameters are also used to explore the temperature dependence of various phases [44].

Appendix B Graphene

Graphene can be described by massless effective charge carriers using a 2D Dirac equation [45]. The effective fine structure constant in graphene is $\alpha\simeq$2.5 [46]. The interaction is of electrodynamics type. Therefore, $\varepsilon_{VS}$=1 because S=0 (see Eq. (24)). Let us investigate whether the present approach could describe the case of graphene.

In the case of massless constituent (m=0), only two quantities describe the system: V₀ and r₀ which are the magnitude and the range of the effective interaction, respectively. Hence, only one dimensionless parameter can be built from these two quantities. In the case of graphene, and following the present lines, it shall be the coupling constant $\alpha$ (Eq. (18)).

The quantality (1) behaves as $\Lambda\rightarrow\infty$ in the case of graphene because m=0: graphene could be considered as a “super”-QL, and the wave functions of the charge careers shall be much delocalized. This is in agreement with theoretical studies on systems obeying to 2D Dirac equations [47].

In the case of graphene, the spin-orbit parameter becomes $\eta$=0 (because m=0), meaning that the spin-orbit effect vanishes ( $\Lambda\rightarrow\infty$ in Eq. (27)). However, the key relation (19) still holds, i.e. Eq. (25) with $\varepsilon_{VS}$=1:

$$\eta\Lambda\alpha^{2}=1$$ (47)

Therefore, graphene case can be described by the present approach, as a limit case where $\eta$=0, $\Lambda$=$\infty$ and $\alpha$ remains finite.