## Pairing phase transition in the odd-A nuclei : identification and classification

Yumeng Wang ( 王 宇 萌 ), 1 Yuhang Gao ( 高 宇 航 ), 1 and Lang Liu ( 刘 朗 ) 1, ∗ 1 School of Science, Jiangnan University, Wuxi 214122, China.

## Abstract

The investigation into the pairing phase transition in the odd-A nucleus 161 Dy utilizes a sophisticated blend of covariant density functional theory and the shell-model-like approach. It is discerned that variations in thermodynamic quantities at the critical temperature do not exclusively align with pairing phase transitions. The presence of an S-shaped heat capacity curve, often interpreted as indication of such transitions, does not offer a definitive confirmation. Additional factors, including the blocking effect, can modify the heat capacity curve and impede the transition process. The pairing phase transition in 161 Dy, occurring around 0.7 to 1.0 MeV, is unequivocally characterized as a first-order transition. Furthermore, an exploration into the impact of varying strengths of pairing correlations on these transitions reveals a nonlinear relationship, adding complexity to the transition dynamics.

∗ liulang@jiangnan.edu.cn

## I. INTRODUCTION

The pairing phase transition is widely recognized as essential for comprehending the thermodynamic properties of atomic nuclei [1]. Due to the precise measurement of level densities at low excitation energies, thermodynamic quantities have been deduced within the canonical ensemble framework, thereby confirming the S-shaped curve of heat capacity as a function of temperature [2-5]. The significance of this S-shaped curve has been extensively explored across various theoretical models, including the nuclear shell model [6-8], the mean field model [9-13], and others [14, 15]. It is considered the hallmark of the pairing phase transition from the superfluid phase to the normal phase. However, in the superfluid phase, particle number conservation is violated by the BCS theory or the Bogoliubov transformation within the mean field model. The importance of conservation laws in studying phase transitions of finite systems has been elucidated in Refs. [16-18]. To address these drawbacks, we employ the shell-model-like approach, which ensures particle number conservation and accurately deals with the blocking effect [19-21]. Additionally, covariant density functional theory (CDFT) has effectively described numerous nuclear phenomena [22-25]. The shell-model-like approach has been integrated within the CDFT framework [21]. Previous studies have examined the thermodynamic properties of even-even nuclei and odd-A nuclei using CDFT and the shell-model-like approach [26, 27].

Although the pairing and pairng phase transition has been extensively investigated in many aspects [28-31], few studies have focused on the order of these transitions, and it remains to be more clearly defined. The initial classification of general types of phase transitions, introduced by Paul Ehrenfest in 1933, has played a pivotal role in the thermodynamic study of critical phenomena [32]. Subsequently, Lee and Yang proposed a theorem on the distribution of roots of the grand partition function, predicting its broad applicability even in small systems [33]. In Refs. [34, 35], a classification scheme for phase transitions in finite systems was proposed, based on the distribution of zeros (DOZ) of the canonical partition function in complex temperature. This method has been successfully applied to determine the order of phase transitions in 162 Dy [36]. To our knowledge, it has not yet been utilized to study the order of the pairing phase transition in odd-A nuclei.

This work aims to further investigate the microscopic mechanism underlying the S-shaped heat capacity curve, examining whether its appearance entirely corresponds to the occur-

rence of pairing phase transitions. If a pairing phase transition does occur, the order of the transition in odd-A nuclei remains unclear, and the influence of the blocking effect on the phase transition needs to be studied. Additionally, this research will explore whether and how the strength of the pairing correlation affects the pairing phase transition, an aspect that has not been addressed in previous studies.

In this work, we use 161 Dy as an example. By employing the shell-model-like approach and covariant density functional theory, combined with the effective interaction PC-PK1 [37], we investigate the thermodynamic properties and the influence of the blocking effect on the pairing phase transition. The order of the pairing phase transition is determined according to a classification scheme based on the distribution of zeros of the pairing partition function in the complex temperature plane. Furthermore, we construct an ideal odd system and an ideal even system to explore the impact of pairing correlation strength on the pairing phase transition, as evidenced by changes in the order of the transition.

This paper is organized as follows: In Sec.II, we briefly introduce the theoretical framework. Sec.III presents the results and discussion. Finally, Sec.IV contains the summary.

## II. THEORETICAL FRAMEWORK

In this section, we introduce briefly the basic formalism of the CDFT theory, which has been summarized in Ref.[37]. CDFT begins with a Lagrangian. From this Lagrangian, the Dirac equation for nucleons with local scalar S ( r ) and vector V µ ( r ) potentials can be deduced as

$$\left [ \gamma _ { \mu } ( i \partial ^ { \mu } - V ^ { \mu } ) - ( m + S ) \right ] \psi _ { \xi } = 0,$$

where

$$S ( r ) = \Sigma _ { S }, V ( r ) = \Sigma ^ { \mu } + \vec { \tau } \cdot \vec { \Sigma } ^ { \mu } _ { T V },$$

and ψ ξ is Dirac spinor.

In this study, SLAP is applied within the CDFT framework to deal with the pairing correlations by diagonalizing the Hamiltonian in the multi-particle configurations (MPCs) space, which is constructed by the single-particle levels obtained through the CDFT method.

The total Hamiltonian reads

$$H = H _ { \text{s.p.} } + H _ { \text{pair} }$$

/negationslash

$$= \sum _ { i } \varepsilon _ { i } a _ { i } ^ { + } a _ { i } - G \sum _ { i, j > 0 } ^ { i \neq j } a _ { i } ^ { + } a _ { i } ^ { + } a _ { j } a _ { j },$$

where ε i is the single-particle energy gained from the Dirac equation (1), ¯ is the time-reversal i state of i , and G represents constant pairing strength.

For a system with an odd particle number N = 2 n +1, it is necessary to block the state occupied by the odd nucleons. The MPCs for s = 0 states can be expressed as,

$$| c _ { b } c _ { 1 } \bar { c } _ { 1 } \cdots c _ { n } \bar { c } _ { n } \rangle = a _ { b } ^ { + } ( a _ { c _ { 1 } } ^ { + } a _ { \bar { c } _ { 1 } } ^ { + } \cdots a _ { c _ { n } } ^ { + } a _ { \bar { c } _ { n } } ^ { + } ) | 0 \rangle,$$

in which b denotes the single particle level blocked by the odd nucleon.

During the calculation process, a cutoff energy E c is needed to truncate the MPCs space. The configurations with energies E m -E 0 ≤ E c are used to diagonalize the Hamiltonian (3), in which E m and E 0 are the energies of the m th configuration and the ground-state configuration, respectively.

After the diagonalization of the Hamiltonian (3), the nuclear many-body wave function is shown as

$$\mathbf n \, \text{as} \\ | \psi _ { \beta } \rangle & = \sum _ { c _ { 1 } \cdots c _ { n } } v _ { \beta, c _ { 1 } \cdots c _ { n } } | c _ { 1 } \bar { c } _ { 1 } \cdots c _ { n } \bar { c } _ { n } \rangle \\ & \quad + \sum _ { i, j } \sum _ { c _ { 1 } \cdots c _ { n - 1 } } v _ { \beta ( i j ), c _ { 1 } \cdots c _ { n - 1 } } | i \bar { j } c _ { 1 } \bar { c } _ { 1 } \cdots c _ { n - 1 } \bar { c } _ { n - 1 } \rangle \\ & \quad + \cdots,$$

where β = 0 corresponds the ground state, β = 1 2 3 , , , . . . correspond the excited states with the excitation energy E β , and v β means the coefficient after diagonalization.

Moreover, the pairing gap can be calculated by Refs. [38-40]

$$\lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda
 \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \text{$\mu$} \quad,$$

The thermodynamic quantities of the pairing interaction are calculated in the canonical ensemble [41]. The partition function Z can be calculated by

$$Z = \sum _ { \beta = 0 } ^ { \infty } \eta ( E _ { \beta } ) \, e ^ { - E _ { \beta } / T }.$$

where E β is the excitation energy which could be acquired through the CDFT + SLAP method, and the corresponding level density η E ( β ) is taken as 2 , s i.e., the degeneracy of each state.

The classification scheme for phase transitions relies completely on Ref. [42]. Only the main formulas of the classification scheme will be given below. The scheme relies on the DOZ of the canonical partition function in the complex temperature plane.

The inverse complex temperature is defined as:

$$\mathcal { B } = \beta + i \tau,$$

where β = 1 /T and τ corresponds the imaginary part of the inverse complex temperature which is measured in MeV -1 . The zeros of the canonical partition function typically line up on curves in the complex temperature plane. The zeros are denoted by ( β , τ j j ). If there are n zero points, j = 1 . . . n and j increases with the increasing distance from the real axis.

If the first three zeros closest to the real axis have been determined, the average inverse distance between zeros can be calculated as

$$\Phi ( \widetilde { \tau } _ { j } ) = \frac { 1 } { d _ { j } }, \\ \div \sqrt { ( \beta _ { j + 1 } - \beta _ { j } ) ^ { 2 } + ( \tau _ { j + 1 } - \tau _ { j } ) ^ { 2 } }. \text{ The function $\Phi$ can then be}$$

where ˜ τ j = ( τ j + τ j +1 ) / 2, d j = √ ( β j +1 -β j ) 2 +( τ j +1 -τ j ) 2 . The function Φ can then be approximated in the vicinity of the real axis by a power law of τ j only,

$$\Phi ( \tau _ { j } ) \in \tau _ { j } ^ { \alpha },$$

then α can be calculated by means of

$$\alpha = \frac { \ln \Phi ( \tau _ { 3 } ) - \ln \Phi ( \tau _ { 2 } ) } { \ln \tau _ { 3 } - \ln \tau _ { 2 } }.$$

According to the definition in Ref. [42], when α &lt; 0, the system displays a first-order phase transition. When 0 &lt; α &lt; 1, the phase transition is of second order. When α &gt; 1, higherorder phase transition occurs.

## III. RESULTS AND DISCUSSION

A. 161 Dy

For neutrons and protons, 16 single-particle levels and 9 valence nucleons are selected to construct a multi-particle configuration space, including 5 levels on and below the Fermi

surface and 11 energy levels above the Fermi surface. PC-PK1 is chosen as the effective interaction. For neutrons, the Fermi surface is occupied by a single particle, and the maximum seniority number s is 8. The span of the upper and lower energy levels of the model space is 7.56 MeV. The corresponding pairing strengths of neutrons and protons are 0.29 MeV and 0.32 MeV, which can reproduce the experimental odd-even mass differences well. In this paper, the treatment of blocking is consistent with the method in Ref. [27]. 'Bk1' corresponds to the blocking of the Fermi surface, while 'bkn' corresponds to the blocking of the (n-1)th energy level above the Fermi surface.

FIG. 1. (Color online) The pairing gap (a), heat capacity (b), average energy (c), and partition function (d) of neutrons (black solid lines) and protons (red dashed line) for 161 Dy as functions of temperature.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-b43bfa9e-5e82-4f45-bb1e-a8a791c6084c.jpg)

Figure 1 illustrates the behavior of the pairing gap (a), heat capacity (b), average energy (c), and partition function (d) of neutrons (indicated by black solid lines) and protons (represented by red dashed lines) in 161 Dy as functions of temperature. Subfigure (a) reveals

obvious increase in the neutron curve at temperature below 0.3 MeV, contrasting with the almost constant trend observed in the proton curve within this temperature range. This behavior in the neutron curve indicates the pairing re-entrance phenomenon, wherein the pairing correlation in nuclei with an odd number of neutrons or protons is bolstered as temperature rises. This enhancement is primarily attributed to the influence of the singleparticle level structure near the level occupied by the odd nucleons, which weakens the blocking effect and reinforces the pairing correlation. Regarding heat capacity, both protons and neutrons exhibit an S-shaped curve. The initial turning point occurs at approximately 0.3 MeV for protons and 0.2 MeV for neutrons, aligning with the turning points of the pairing gap and average excitation energy. A second turning point emerges at around 1 MeV for protons and 0.7 MeV for neutrons, corresponding to the critical temperature of the phase transition. Additionally, the neutron curve demonstrates fluctuations in the lowtemperature region. In the average energy diagram (c), both curves are nearly zero at low temperatures (around 0.25 MeV for neutrons and 0.3 MeV for protons), gradually increasing in a near-linear way thereafter.

Figure 2 (a) displays the heat capacity of neutrons in 161 Dy, considering total and different blocking scenarios, as functions of temperature. The fluctuations observed at low temperatures primarily stem from the contribution of 'bk1', attributable to the single-particle level structure. In the case of 'bk2', corresponding to the blocking of the first level above the Fermi surface, Figure 2 (b) illustrates that no levels exist around this energy level. Consequently, an energy input is required to excite nucleons, resulting in near-zero heat capacity at low temperatures. Subsequent blockings involve energy levels occupied by odd nucleons, which have less impact on the shape of the heat capacity curve. Hence, accurate treatment of Fermi surface blocking and its neighboring levels is crucial for investigating the relevant properties of odd-A nuclei.

To gain a detailed microscopic understanding of the nuclear pairing transition, it is valuable to examine the number of Cooper pairs that become broken with increasing temperature in 161 Dy. In this analysis, we follow the definition of the seniority component as outlined in Ref. [26]

$$\chi _ { s } = Z ^ { - 1 } \sum _ { \beta \in \{ s \} } \eta ( E _ { \beta } ) e ^ { - E _ { \beta } / T }.$$

The seniority component illustrates the contribution of excited states with specific seniority

FIG. 2. (Color online) (a) The heat capacity of neutrons with the total and different levels blocked for 161 Dy as functions of temperature. (b) The single particle levels of neutrons for 161 Dy. The levels in the orange area are used to construct the SLAP MPCs. λ corresponds to the Fermi surface.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-f0f5fdaa-2be5-4051-bf60-ab7d5882606e.jpg)

numbers to the total average energy. Figure 3 presents the seniority components for different seniority numbers s = 0, 2, 4, 6, 8 of neutrons (a) and s = 0, 2, 4, 6, 8, 10 of protons (b) in 161 Dy as functions of temperature. It is observed that s = 0 states contribute nearly 100% below T ≈ 0.3 MeV for protons and 0.2 MeV for neutrons, aligning with the vanishing heat capacity illustrated in Fig. 1. As the temperature exceeds 0.3 MeV for protons and 0.2 MeV for neutrons, the contribution of s = 0 states begins to decrease. Simultaneously, the curve for s = 2 states starts to rise, corresponding to the first turning point in the heat capacity curve and the temperature where the pairing gap curve begins to decline and the average energy curve begins to rise. The curves for s =2peak at 0.8 MeV for neutrons and 1 MeV for protons, respectively, corresponding to the second turning point of the heat capacity curve,

FIG. 3. (Color online) Seniority components with different seniority numbers s = 0,2,4,6,8 of neutrons (a) and s = 0,2,4,6,8,10 of protons (b) for 161 Dy as functions of temperature.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-20a7cfea-5344-47ac-8ea8-ea25a008da5d.jpg)

indicating the pairing phase transition from the superfluid phase to the normal phase. This transition is responsible for the S-shape of the heat capacity curve. Furthermore, around T ≈ 0.5 MeV, the contribution of s = 4 states begins to increase.

The level density serves as a valuable tool for characterizing pairing transitions in hot nuclei. To further elucidate the microscopic mechanisms underlying the S-shaped heat capacity curve, we investigate the number of excited states of nuclei in terms of different seniority contributions. Figure 4 illustrates the level densities with different seniority numbers s = 0, 2, 4, 6, 8 for neutrons (a) and s = 0, 2, 4, 6, 8, 10 for protons (b) in 161 Dy as functions of excitation energies. From the curves, distinct protuberances are evident for s = 2, 4, 6, and 8 for neutrons, and s = 2 and 4 for protons. Notably, these protuberances are absent at higher seniority numbers. This observation suggests that the initial fluctuations around 4 MeV are primarily attributed to the contribution of one-pair-broken states. The number of excited states increases rapidly due to the emergence of one-pair-broken states with increasing excitation energies. However, when T ≥ 4 MeV, the growth rate of the curves diminishes because no additional one-pair-broken states can be formed. This phenomenon corresponds

FIG. 4. (Color online) Level densities with different seniority numbers s = 0,2,4,6,8 of neutrons (a) and s = 0,2,4,6,8,10 of protons (b) for 161 Dy as functions of the excitation energies.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-5fbf30c1-7ec7-4d79-947d-53c65b3bea47.jpg)

to the first S-shaped feature in the heat capacity curve. Subsequently, with the emergence of the second fluctuation, the appearance of two-pair-broken states leads to rapid increase, followed by a decline after 8 MeV excitation energy for neutrons and 7 MeV for protons. Other protuberances may indicate additional S-shaped features in the heat capacity curve at higher energies.

Figure 5 illustrates contour plots of the partition function of neutrons and protons in the complex temperature plane. Zero points at β ≈ 1 05 MeV . -1 for protons and 1.3 MeV -1 for neutrons effectively reflect the critical temperature of the phase transition. The corresponding critical temperature of protons and neutrons are 0.95 MeV and 0.77 MeV, and the experimental value for 161 Dy is 0.52 MeV [2]. The contour plots of the partition function of neutrons for 161 Dy with different blocking scenarios in the complex temperature plane are presented in Fig. 6. It is evident that the partition functions of 'bk1' and 'bk2' contribute the most to the overall partition function. Conversely, patterns of blocking higher energy levels are similar, contributing minimally to the overall partition function. This observa-

FIG. 5. (Color online) The contour plots of the partition function of neutrons and protons for 161 Dy in the complex temperature plane.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-24f7cb08-6123-4b39-a3fc-398861b89699.jpg)

tion aligns with earlier discussions regarding the effects of different blocking scenarios on heat capacity. Furthermore, in cases of blocking higher energy levels, the partition function predominantly occupies the zero-point domain. This suggests weakening of the nuclear superfluidity, with the nucleus primarily existing in the normal phase.

According to Fig. 5, the minimums are determined iteratively as the zero points of the partition function. The distributions of zero points of the partition function for neutrons and protons in the complex temperature plane are displayed in Fig. 7, with the zero points of the pairing phase transition indicated by the solid line. By calculation using Eq. 11, the α values corresponding to the neutrons and protons of 161 Dy are -2.04 and -5.66, respectively. Based on the phase transition classification method, this indicates a first-order phase transition.

Figure 8 illustrates the contour plots of the partition function for 161 Dy in the complex temperature plane with different seniority numbers s = 0, 2, 4, 6, 8, and 10, where the upper and lower lines correspond to neutrons and protons, respectively. In the first column of the seniority number s = 0 (fully paired states), zero points exist only in the region where β is very small. Additionally, there are no zero points in the region of 0 to 2 MeV ( β &gt; 0.5 MeV -1 ), indicating that the nucleus is entirely in the superfluid phase. As we move to the s = 2 (one pair broken) state, the zero point domain begins to appear with a clear

FIG. 6. (Color online) The contour plots of the partition function of neutrons for 161 Dy with different blocking in the complex temperature plane.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-96e3fbac-76ce-4a64-96f9-f248bf9c8192.jpg)

boundary, signifying the coexistence of the superfluid and normal phases. Subsequently, differences emerge between neutrons in odd systems and protons in even systems. The zero domain for neutrons is smaller than that for protons in the s = 2 state, primarily due to the blocking effect of the odd nucleon hindering the evolution of the pairing phase transition. With increasing seniority number, the zero point domain gradually expands, and the normal phase gradually dominates. The evolution of the pairing phase transition can still be effectively demonstrated. By the time we reach the seniority number s = 8 or 10 (completely broken) state, the image is nearly entirely covered by the zero point domain, indicating that the nucleus is entirely in the normal phase.

FIG. 7. The DOZ of the partition function of the neutrons (a) , protons (b) for 161 Dy in the complex temperature plane. The solid line and dashed line represent two kinds of phase transitions.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-bea998c7-216d-492c-9ba6-dea8ee9ca33e.jpg)

FIG. 8. (Color online) The contour plots of the partition function with different seniority numbers s = 0,2,4,6,8,10 of neutrons (top row) and protons (bottom row) for 161 Dy in the complex temperature plane.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-53e6e26b-be99-4cba-baff-55fa662ac3db.jpg)

## B. The Toy Model

The study above found that the structure of single particle levels affects the heat capacity curve. Under the SLAP framework, single particle levels are derived by Eq. 1. In fact, it may be affected by many factors, such as nuclear deformation, shell structure, etc. Since

this paper only focuses on the influence of pairing correlation, it is necessary to remove the influence of other factors, so an ideal model can be established. In this study, the ideal models with odd particles and even particles are constructed with controllable pairing strength, and thermodynamic quantities are investigated under the same framework. For the even system, a configuration space with 10 levels and 10 particles is selected, with 5 levels above the Fermi surface. The energy levels range from 1 to 10 MeV, and the pairing strength varies from 0 to 2 MeV in increments of 0.1 MeV. The cut-off energy is set to 45 MeV, resulting in a configuration space dimension of 10 . In the case of the odd system, a configuration space 4 with 10 levels and 11 particles is chosen, with 4 levels above the Fermi surface. The single particle levels, pairing strength, and cut-off energy are consistent with the even system. The configuration space dimension is 1 4 . × 10 . By employing these ideal models, the study aims 4 to investigate the thermodynamic quantities under the same framework, thus allowing for clearer understanding of the sole influence of pairing correlation on the system.

FIG. 9. (Color online) The heat capacity with different pairing strengths for the ideal even system (a) and the ideal odd system (b) as functions of temperature.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-8c9ebdea-f184-4962-b581-7f0d5a724ab3.jpg)

Figure 9 illustrates the heat capacity with varying pairing strengths for both the ideal even system (a) and the ideal odd system (b) as functions of temperature. As the pairing strength G increases, the temperatures corresponding to the first and second turning points of the

curves, as well as the peak values of the curves, also increase. In Fig. 9 (b), fluctuations are observed in the low-temperature region, similar to the situation discussed for 161 Dy. These fluctuations are primarily attributed to the pairing re-entrance phenomenon in the odd-nucleon system.

FIG. 10. (Color online) The DOZ of the partition function with different pairing strengths for the ideal even system in the complex temperature plane.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-86e619e8-48d6-4d03-a52a-3c054bccf653.jpg)

According to the contour plots of the partition function with varying pairing strengths for the ideal even system in the complex temperature plane, the Distribution of Zeroes (DOZ) of the partition function is displayed in Fig. 10. When the pairing strength is 0.1, a zero point is observed outside the plotted range, as indicated in the insert figure. In contrast, when the pairing strength is 0 (indicating no pairing correlation), there is no zero point in the corresponding β region. However, the heat capacity curve for G = 0 still exhibits an S shape, suggesting that the appearance of an S-shaped heat capacity curve does not solely

TABLE I. The value of α with different pairing strengths for the ideal even system.

| G (MeV)   |   0.1 |   0.2 |   0.3 |   0.4 |   0.5 |    0.6 |   0.7 |   0.8 |   0.9 |   1.0 |
|-----------|-------|-------|-------|-------|-------|--------|-------|-------|-------|-------|
| α         | -3.98 | -4.73 | -5.45 | -3.86 | -5.46 | -10.49 | -39.3 |   7.1 |  2.69 |  0.94 |

TABLE II. The value of α with different pairing strengths for the ideal odd system.

| G (MeV)   |   0.1 |   0.2 |   0.3 |   0.4 |   0.5 |   0.6 |    0.7 |   0.8 |   0.9 |   1.0 |
|-----------|-------|-------|-------|-------|-------|-------|--------|-------|-------|-------|
| α         |  -3.4 | -3.66 | -3.98 | -4.66 | -4.61 | -8.37 | -18.59 | 35.01 |   7.2 |   3.2 |

correspond to the pairing phase transition. As the pairing strength G increases, the branch of zero points near the real axis undergoes significant changes. The zero point corresponding to the peak value of the heat capacity curves gradually shifts to the region of smaller β due to the increasing temperature.

The corresponding α values are calculated and presented in Table I. For pairing strengths G ≤ 0.7, a first-order phase transition is observed. When G is between 0.8 and 0.9, a higherorder phase transition is identified. Finally, for G = 1, a second-order phase transition is noted. This analysis reveals a nonlinear relationship between the pairing strength and the order of the pairing phase transition, consistent with findings in Ref. [41].

Figure 11 illustrates the Distribution of Zeroes (DOZ) of the partition function with varying pairing strength for the ideal odd system in the complex temperature plane. Notably, a branch is observed when the pairing strength is 0. Presence of this branch is primarily attributed to changes in the zero point distribution caused by the blocking effect, and it does not correspond to a real pairing phase transition. This behavior contrasts with the case of even systems.

Based on the zero distribution depicted in Fig. 11, the corresponding α values are calculated and presented in Table II. For pairing strengths G ≤ 0.7, a first-order phase transition is identified. However, for G between 0.8 and 1, a higher-order phase transition is observed. This analysis demonstrates that the pairing strength also influences the pairing phase transition in the odd system.

FIG. 11. (Color online) The DOZ of the partition function with different pairing strengths for the ideal odd system in the complex temperature plane.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00423v1-figure-a03e55c0-41df-425a-98e0-3522b4ea761a.jpg)

## IV. SUMMARY

In summary, this work investigates the pairing phase transition in the odd-A nucleus 161 Dy using the covariant density functional theory and the shell-model-like approach. Thermodynamic quantities are evaluated within the canonical ensemble theory, and the phase transition order is determined using a classification scheme based on the distribution of zeros of the partition function in the complex temperature plane. The impact of different pairing strengths on the pairing phase transition is also examined.

The heat capacity curve of 161 Dy exhibits an S shape, indicating a pairing phase transition from the superfluid phase to the normal phase. Extending the calculation to the complex temperature plane reveals that the pairing phase transition is of first order. Blocking the Fermi surface and the first level above it significantly contribute to the total partition func-

tion, highlighting the importance of considering the blocking of the nearest levels above the Fermi surface in the study of odd-A nuclei. The contour plots of the partition function with different seniority show that the blocking effect of the odd nucleon hinders the evolution of the pairing phase transition.

Furthermore, ideal models with odd and even particles are constructed to explore the influence of different pairing strengths on the pairing phase transition in the complex temperature plane. In the ideal even system, the absence of pairing correlation does not lead to a zero point corresponding to the critical temperature in the large β region, yet the heat capacity curve still exhibits an S shape. This suggests that the appearance of an S-shaped heat capacity curve does not exclusively indicate a pairing phase transition. Conversely, in the ideal odd system, the presence of a zero point branch in the partition function plot is due to the blocking effect. The analysis also demonstrates that the pairing strength has nonlinear influence on the pairing phase transition. In the ideal even system, the phase transition order changes from first order to higher order and then to second order as the pairing strength increases. For the ideal odd system, the phase transition order changes from first order to higher order with increasing pairing strength. In the near future research, we will focus more on the experimental spectroscopic properties and other properties to study the pairing.

## ACKNOWLEDGMENTS

This work was supported by National Natural Science Foundation of China under Grant No. 11775099 and the Jiangnan University Basic Research Program (JUSRP202406002).

- [1] A. Bohr and B. Mottelson, Nuclear Structure , Vol. I (World Scientific, 1998).
- [2] A. Schiller, A. Bjerve, M. Guttormsen, M. Hjorth-Jensen, F. Ingebretsen, E. Melby, S. Messelt, J. Rekstad, S. Siem, and S. W. Ødeg˚rd, Phys. Rev. C a 63 , 021306 (2001).
- [3] E. Melby, M. Guttormsen, J. Rekstad, A. Schiller, S. Siem, and A. Voinov, Phys. Rev. C 63 , 044309 (2001).
- [4] E. Melby, L. Bergholt, M. Guttormsen, M. Hjorth-Jensen, F. Ingebretsen, S. Messelt, J. Rekstad, A. Schiller, S. Siem, and S. W. Ødeg˚rd, Phys. Rev. Lett. a 83 , 3150 (1999).

- [5] M. Guttormsen, R. Chankova, M. Hjorth-Jensen, J. Rekstad, S. Siem, A. Schiller, and D. J. Dean, Phys. Rev. C 68 , 034311 (2003).
- [6] S. Rombouts, K. Heyde, and N. Jachowicz, Phys. Rev. C 58 , 3295 (1998).
- [7] S. Liu and Y. Alhassid, Phys. Rev. Lett. 87 , 022501 (2001).
- [8] K. Langanke, D. J. Dean, and W. Nazarewicz, Nucl. Phys. A 757 , 360 (2005).
- [9] J. L. Egido, L. M. Robledo, and V. Martin, Phys. Rev. Lett. 85 , 26 (2000).
- [10] B. K. Agrawal, T. Sil, J. N. De, and S. K. Samaddar, Phys. Rev. C 62 , 044307 (2000).
- [11] N. Sandulescu, O. Civitarese, and R. J. Liotta, Phys. Rev. C 61 , 044317 (2000).
- [12] Y. F. Niu, Z. M. Niu, N. Paar, D. Vretenar, G. H. Wang, J. S. Bai, and J. Meng, Phys. Rev. C 88 , 034308 (2013).
- [13] J. J. Li, J. Margueron, W. H. Long, and N. Van Giai, Phys. Rev. C 92 , 014302 (2015).
- [14] M. Guttormsen, M. Hjorth-Jensen, E. Melby, J. Rekstad, A. Schiller, and S. Siem, Phys. Rev. C 63 , 044301 (2001).
- [15] M. Guttormsen, M. Hjorth-Jensen, E. Melby, J. Rekstad, A. Schiller, and S. Siem, Phys. Rev. C 64 , 034319 (2001).
- [16] K. Esashika, H. Nakada, and K. Tanabe, Phys. Rev. C 72 , 044303 (2005).
- [17] D. Gambacurta, D. Lacroix, and N. Sandulescu, Phys. Rev. C 88 , 034324 (2013).
- [18] C. Esebbag and J. L. Egido, Nucl. Phys. A 552 , 205 (1993).
- [19] J.-Y. Zeng and T.-S. Cheng, Nucl. Phys. A 405 , 1 (1983).
- [20] J. Zeng, T. Jin, and Z. Zhao, Phys. Rev. C 50 , 1388 (1994).
- [21] J. Meng, J. Y. Guo, L. Liu, and S. Q. Zhang, Front. Phys. China 1 , 38 (2006).
- [22] P. Ring, Prog. Part. Nucl. Phys. 37 , 193 (1996).
- [23] D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, and P. Ring, Phys. Rept. 409 , 101 (2005).
- [24] J. Meng, H. Toki, S. Zhou, S. Zhang, W. Long, and L. Geng, Prog. Part. Nucl. Phys. 57 , 470 (2006).
- [25] J. Meng, J. Peng, S. Zhang, and P. W. Zhao, Front. Phys. 8 , 55 (2013).
- [26] L. Liu, Z. H. Zhang, and P. W. Zhao, Phys. Rev. C 92 , 044304 (2015).
- [27] T. Yan, Y. L. Lin, and L. Liu, Phys. Rev. C 104 , 024303 (2021).
- [28] M. Horoi and V. Zelevinsky, Phys. Rev. C 75 , 4303 (2007).
- [29] Y.-F. Gao, B.-S. Cai, and C.-X. Yuan, Nucl. Sci. Tech. 34 , 9 (2023).
- [30] B. S. Cai and C. X. Yuan, Nucl. Sci. Tech. 34 , 271 (2023).

- [31] H. Tajima, T. Hatsuda, P. van Wyk, and Y. Ohashi, Sci. Rep. 9 , 18477 (2019).
- [32] G. Jaeger, Arch. Hist. Exact. Sci. 53 , 51 (1998).
- [33] T. D. Lee and C. N. Yang, Phys. Rev. 87 , 410 (1952).
- [34] P. Borrmann, O. M¨lken, u and J. Harting, Phys. Rev. Lett. 84 , 3511 (2000).
- [35] O. M¨ ulken, P. Borrmann, J. Harting, and H. Stamerjohanns, Phys. Rev. A 64 , 013611 (2001).
- [36] Y. H. Gao, Y. L. Lin, and L. Liu, Int. J. Mod. Phys. E 32 , 2350050 (2023).
- [37] P. W. Zhao, Z. P. Li, J. M. Yao, and J. Meng, Phys. Rev. C 82 , 054319 (2010).
- [38] L. F. Canto, P. Ring, and J. O. Rasmussen, Phys. Lett. B 161 , 21 (1985).
- [39] J. L. Egido, P. Ring, S. Iwasaki, and H. J. Mang, Phys. Lett. B 154 , 1 (1985).
- [40] Y. R. Shimizu, J. D. Garrett, R. A. Broglia, M. Gallardo, and E. Vigezzi, Rev. Mod. Phys. 61 , 131 (1989).
- [41] T. Sumaryada and A. Volya, Phys. Rev. C 76 , 024319 (2007).
- [42] A. Schiller, M. Guttormsen, M. Hjorth-Jensen, J. Rekstad, and S. Siem, Phys. Rev. C 66 , 024322 (2002).
Yumeng Wang Yuhang Gao Lang Liu
Pairing phase transition in the odd-A nuclei: identification and classification
The investigation into the pairing phase transition in the odd-A nucleus
161Dy utilizes a sophisticated blend of covariant density functional theory and
the shell-model-like approach. It is discerned that variations in thermodynamic
quantities at the critical temperature do not exclusively align with pairing
phase transitions. The presence of an S-shaped heat capacity curve, often
interpreted as an indication of such transitions, does not offer a definitive
confirmation. Additional factors, including the blocking effect, can modify the
heat capacity curve and impede the transition process. The pairing phase
transition in 161Dy, occurring around 0.7 to 1.0 MeV, is unequivocally
characterized as a first-order transition. Furthermore, an exploration into the
impact of varying strengths of pairing correlations on these transitions
reveals a nonlinear relationship, adding complexity to the transition dynamics.