1. Introduction
This work explores new fermion features, related to order–disorder issues, as reflected by SU2 quantum models. More specifically, we discuss novel many-body thermodynamic traits that refer to the energy cost associated with changes in the system’s particle number. We envision a thermal quantum statistical scenario endowed with suitable fermion–fermion interactions.
As is well known, the notion of entropy is linked with the idea of disorder. In this work, we used a special order quantifier denominated disequilibrium
D (see [
1] and the references therein).
D is thought to be a Euclidean distance in probability space between the specific probability distribution of current interest and the uniform one. The latter can be regarded as reflecting maximal “disorder”, so that the bigger
D value should correspond to the larger possible order degree. This
D description of order yields valuable insights.
1.1. Our Present Fermion Models
Of course, atomic nuclei are paradigmatic examples of many-fermion systems [
2,
3,
4]. We used here two simple exactly solvable fermion models, whose treatment allows insight into the intricacies of the quantum many-body problem, without the need to resort to huge Hamiltonian matrices [
2,
3,
4]. We think that a good way to understand a subject is to first grasp specific, well-chosen cases and, only afterwards, worry about possible generalizations. In fact, it has been of immense usefulness in fermion-related theoretical research to appeal to the exactly solvable Lipkin Model (LM) [
2,
3,
4,
5,
6,
7], which has been fruitfully used to investigate manifold traits of the many-body problem [
8]. The two-energy-level LM (containing
N fermions) revolves about an SU2 algebra generated by operators called quasi-spin ones.
We appeal in this paper to Lipkin model variants that possess exact analytical solutions (instead, the LM needs numerical diagonalization). We can then compare our analytical exact solutions with other solutions that arise from many kinds of approximate N fermion manipulations.
A very important Casimir Operator (CO),
, pervades the workings of both the LM [
5] and our present variants of it.
N fermions are distributed amongst two N degenerated singe-particle levels. Our CO has affixed to it different multiplets. In particular, there exists an unperturbed ground state multiplet of great relevance [
5]. We emphasize that each of our two energy levels is degenerate and can accommodate by itself all our 2
possible fermion locations in the system. Each of these locations is called a site, which can be either occupied or empty. We call sister sites those in the same position in the upper and lower levels.
Cambiaggio and Plastino (CP) [
1,
9] advanced an SU2 × SU2 Lipkin extension so as to deal with the excited Lipkin multiplets. Such an extension led to the enactment, in quasi-spin language, of a BCS-like formalism, which allows mimicking nuclear superconductivity [
10,
11,
12], producing exact analytic results. In the CP-model, the BCS solution coincides with the exact solution. The CPM treatment contemplates a variable particle number.
1.2. Present Objective
Some physicists refer to level-crossings as phase transitions. This is what we also do in what follows. We used novel statistical quantifiers, in a Lipkin-related context, to compare the many-body workings of the pairing interaction with those of a spin–flip force [
3] that arises in a Plastino–Moskowski (PM) model of the last interaction [
13,
14]. These two fermion–fermion interactions individually display “phase transitions”. We wish to compare them. Our exactly solvable juxtaposition (CP and PM models) displays interesting new structural features that generate insight into fermion dynamics. Here, we focus our attention on the statistical characteristics that appear when one varies the number of fermions
N, and as a consequence, changes in energy and in the degree of order ensue.
1.3. Route to Be Followed
First, we establish a statistical indicator of order called the disequilibrium D.
Secondly, we establish how much free energy F it takes to change D. The D changes we are interested in ensue as a result of variations in the value of the coupling constants of the fermions’ Hamiltonian H.
The associated free energy costs of having D variations are quantified via a statistical quantifier, which we call (order variation energy cost).
Since we have two different interaction terms in H, with coupling constants named V and G, respectively, we need to deal with two different types of , which we call and .
Finally, we describe the behavior of and as a function of the number of fermions N. This behavior exhibits interesting features.
2. Description of Our Fermion System
We deal with an unperturbed Hamiltonian
plus two different interactions: a pairing one
(dimensionless coupling constant
G) and a spin–flip one [
3] (dimensionless coupling constant
V):
A canonical ensemble analysis is performed on H. One encounters two attractors (one for and the other for ). For the M attractor, each sister pair of sites is singly occupied. In the attractor, each sister pair is either doubly occupied or empty. This situation can be viewed as an “ordered” one. Instead, we view as “disordered” those instances in which the two energy levels possess distinct occupation numbers.
The Hamiltonian (
1) is of the sort that one encounters in many atomic nuclei: a long-range distance interaction (e.g., a quadrupole one) competes with a short-range one, the pairing interaction. Its pairing component displays the same mathematical traits as an actual nuclear pairing force.
The System’s Statistical Quantifiers S, D, and C
A very important one is the so-called disequilibrium
D (see, for instance, [
15,
16,
17,
18]). Consider two extreme circumstances: (i) perfect order or (ii) total randomness (no correlations) [
15]. Many degrees of correlation exist in intermediate circumstances. In [
15], such degrees were quantitatively characterized in terms of a probability distribution (PD). Maximum randomness (maximal disorder) is associated with a uniform PD (
). The remaining degrees of correlations are represented for different PDs via the Euclidean distances between these PDs and the
. These distances were baptized as disequilibria
D in [
15,
16]. A sort of hierarchy emerges according to the numerical values for
D, which would reflect the existence of privileged states among the accessible ones. By multiplying
D by the entropy quantifier
S, L. Ruiz, Mancini, and Calvet (LMC) [
15] gave birth to one of today’s most-employed notions of structural statistical indicator
C:
It is often stated that
C grasps the “Structural statistical indicator” in the very same fashion that entropy does with the concept of randomness [
15]. If one deals with
M accessible states [
15], one writes
Here,
are the individual normalized probabilities (
) [
15].
D attains its maximum for a perfectly ordered state and vanishes for
. We also have
The LMC’s scheme has enjoyed immense attention, with many applications for the canonical, micro-canonical, and grand canonical Gibbs’ ensembles ([
15,
16,
17,
18] are a tiny sample).
3. Special One-Body Operators and Our Couple of Two-Body Interactions
We face a spin–flip interaction [
3,
14] and a pairing superconducting one [
1], which mimics nuclear superconductivity [
10]. The symmetry background is that of the SU2 × SU2 group. It involves here
N fermions allocated to 2
-fold degenerate single-particle (sp) levels (
). There is an energy gap
(arbitrary units) between them. The system’s sp states are individualized by two quantum numbers:
, with
, and
p is (1) denominated a quasi-spin quantum number and (2) considered as a “site” [
5]. In this vein, one introduces SU2 quasi-spin operators
[
5] in terms of creation and destruction operators
:
Additional SU2 operators were added in [
9]. They are angular momentum-like “pairing” ones:
Obviously,
creates and
destroys a couple of particles, giving a null contribution to the
value. They are paired to
and, thus, do not contribute to the total
value. Any
J-operator will commute with all
Q operators, and vice versa (SU2 × SU2). In this scenario, the associated Hilbert orthonormal basis is characterized by the eigenvalues of
, with eigenstates
. One has [
9]
In the Lipkin model, and also here, we set 2
for a model containing
N fermions [
5]. The unperturbed ground state has
,
, and
[
9]. Going back
, it is written as
cast in terms of
J,
, and a coupling constant
V. Obviously, if
, then
. The eigenstates are
and
where the eigenvalues of
are called
M. It is clear that the
-eigenvalues are [
14]
The energy of the unperturbed (
) gs (ugs)
is called
The most relevant trait of
is that, as
V grows, the system starts going through
level crossings or phase transitions. Indeed, the ground state (at
) ceases being determined by
for
. It comes to be associated with increasingly bigger
values till one reaches
at
[
14]. One can assert that the ket
is an “attractor” for the state of the system as
V is augmented (at
).
We now enter the game the pairing interaction through the interaction expression
, i.e.,
adding to
the pairing bequest
:
displays its own phase transition at
[
9], where the system becomes a superconductor [
9].
, for large enough
G, exhibits a second
“attractor” II state. This superconducting state is associated with
with half of the p sites being doubly occupied.
5. New Quantifier Called the Order Variation Energy Cost
This quantifier tells us how much energy is required to modify the degree of order–disorder in the system as quantified by the disequilibrium D. Our main interest in this work is to ascertain how this quantifier depends on the number of fermions N in the system. In other words, tells us about many-body energy cost effects.
We deal here with two dimensionless coupling constants
and
. A variation from
X to
will result in a change in
D after performing some work
W on the system. This is the change to which we pay attention from now on. In the wake of [
20], we appeal to an efficiency degree
defined as
where
is the work performed on the system as a result of the
change. The work is measured in terms of the free energy
F.
becomes the decrease (growth) in
D that arises from each unit of work (measured in terms of
F) performed on the system. If
is negative, it is the system itself that does the work. We think of quasi-static processes.
We always select the temperature Kelvin. This is a very low value for nuclear physics. This fact is of the essence for our purposes, since we wish to study fermions’ ground state dynamics using statistical mechanics tools. The low T allows for such a desire without unduly disturbing the dynamics.
The quantity X is
V (and take )
G (and take )
We are interested in
involving the Helmholtz free energy
F. More explicitly,
Let us reiterate that the changes in
F can be linked to the work performed on the system
[
20]. Notice that, if
(alternatively
) decreases, it takes energy for this to happen.
6. Results for Our New Quantifier
First of all, note that our statistical quantifiers
S,
F,
D,
, etc., all depend on both
V and
G (dimensionless coupling constants). Thus, at finite temperature, a change
in the statistical quantifier
X that arises out of a variation, say
in
G, will be also connected with the actual
V value, that is
. Thus, in a statistical quantifier
X variation, one observes the emergence of a subtle interplay between
G and
V that is not imaginable by a mere inspection of the Hamiltonian form. Thus, we encounter a statistical correlation
in
entailing, we repeat, a purely statistical interconnection between
V and
G that is not given by the pertinent Hamiltonian. This is an artifact of working at finite temperatures, which is of the essence for our present purposes Importantly enough, such an interconnection will affect the energies’ costs, as we display in
Figure 1 and
Figure 2. This in turn will make it evident how different the two dynamics at play here (superconducting and spin–flip) are. We repeat that
is an order–disorder indicator. A more holistic interpretation of our scenario would posit that changing
X modifies the order–disorder arrangement of the system. Such a modification should then involve both
V and
G.
Figure 1 depicts the order–disorder indicators
(upper figure) and
(lower figure) versus
V at
and
for different
N values. In the first case, there is no energy cost, except at the special critical values
, where remnants of the spin–flip phase transitions are evident. We can call these remnants
singularities. The moral here is that
V critical values are detected by
. In the second instance (lower figure), we see clearly how these critical
V values affect
as well. How is this? This is by considerably lowering the associated energy cost of changing
G. Thus, we encounter a perhaps unexpected subtle interplay between
V and
G in the behavior of the energy cost. Finite temperature force
V and
G interact among themselves.
Let us pass now a quite similar order–disorder panorama, as exhibited by
Figure 2, whose protagonists are the variations of the energy cost with
G. We see that a similar picture emerges as that provided by
Figure 1. The
interrelation, a purely statistical interconnection between
V and
G that is not given by the pertinent Hamiltonian, is at work. First, at the critical
G values, the energy cost
grows. (There is some numerical noise in the graph.) Here, the salient feature is that the energy cost
descends to zero at the critical values
.
Finally, in
Figure 3 and
Figure 4, we depict (versus
N) the
N-dependent critical values of the coupling constants, respectively,
and
. As is well known [
1,
3], for fixed
N, the system displays a spin–flip phase transition at
and a superconducting one at
. Naturally, we speak of “critical values” for the coupling constants. At such values, as we saw above, our statistical quantifiers
and
display “singularities”. Thus, we may say that
Figure 3 and
Figure 4 display the
-singularities or critical points versus
N, as we state in the captions. The values of these critical points diminish as
N grows.
7. Conclusions
In the present work, we proceeded as follows.
First, we established a statistical indicator of order called the disequilibrium D.
Secondly, we ascertained how much free energy F it takes to change D. The particular changes we were interested are those that ensue as a result of variations in the value of the coupling constants of the fermions’ Hamiltonian H.
The associated free energy costs of having these D variations were quantified via a statistical quantifier, which we called (order variation energy cost).
Since we have two different interaction terms in H, with coupling constants named V and G, respectively, we had to deal with two different types of , which we called and .
Finally, we described the behavior of and as a function of the number of fermions N. This behavior exhibited interesting features.
Summing up, we here investigated the interconnection between statistical order and disorder. Our main tool to this effect was the new quantifier characterizing alterations in an order–disorder quantifier D. We focused our attention on changes in the total number of fermions N.
More specifically, we analyzed the behavior of when the coupling constants V and/or G underwent modifications and studied how the associated changes depend on N. We compared the statistical order–disorder traits of our two kinds of fermion–fermion interactions. They are indeed different, when analyzed through the lens provided by . In particular, they differ in relation to their response to variations in N.