Computing Transfer Amplitude between BCS-Pair Condensates

Z. Y. Xia

Department of Physics, University of Shanghai for Science and Technology, Shanghai, China.

DOI: 10.4236/jamp.2022.107145                

Abstract

I propose a new algorithm that uses recursive relations to compute spectroscopic factor, pair transfer amplitude and cluster transfer amplitude. I demonstrate the algorithm that it can be calculated very quickly and stored within small computer memory consumption. In BCS case, the particle number is always conserved and the time-consuming projection is avoided. I drive analytical expressions for the pair transfer amplitude and the cluster transfer amplitude expressed by asymmetry many-pair density matrix. This algorithm practically could be used in all of the nuclear double beta decay fields, the heavy cluster emission fields and the single-nucleon transfer of the odd-mass isotopes.

Keywords

BCS, Spectroscopic Factor, Pair-Transfer Amplitude, Cluster-Transfer Amplitude

1. Introduction

The Bardeen-Cooper-Schrieffer (BCS) theory was first proposed as a superconducting microscopic theory [1] [2]. Later, it was used to treat pairing correlations in nuclear physics [3] [4]. It is still one of the “standard” treatments after fifty years [5], owing to its simplicity and the ease with which higher-order correlations could be added [for example, via quasiparticle random-phase approximation (QRPA)]. However, as compared to macroscopic quantum systems, the theory applied to finite nuclei has two major drawbacks. First, it breaks particle-number conservation by introducing quasiparticles. Second, the BCS theory needs a minimum pairing strength for nuclear systems with finite level spacing. It only yields trivial (vanishing) solutions below that strength, whereas pairing always has an impact in reality.

The “pair condensate” [Equation (3), with determinate particle number] is a frequent improvement over the BCS “quasiparticle vacuum” as the variational ground state [6]. A method was recently proposed under the pair condensate to calculate pair-hopping amplitudes by computing the many-pair density matrices [7]. In the example of the spherical nuclear shell model, it is demonstrated that the many-pair density matrices can be calculated fast using recursive relations and stored easily in computer memory. The method is demonstrated in semi-magic nuclei ^46,48,50Ca, ¹¹⁶Sn, and ¹⁸²Pb.

However, this method can only compute the pair-hopping amplitude of one same nucleus with fixed particle number. Nuclear reactions such as the nuclear double beta decay, single-nucleon transfer and cluster decay of which the amplitudes are of two different nuclei play a crucial role in nuclear physics and particle physics.

One of the pair transfer nuclear reactions is the nuclear double beta decay. Mayer was the first to explore the nuclear double beta decay mode with obvious lepton number conservation in 1935 [8]. On the particle physics side, the inverse half-life of double beta decay is expressed as a product of a phase-space factor and the relevant double beta decay nuclear matrix element, which is free of unknown parameters. Thus, the value of the double beta decay nuclear matrix elements is directly determined by the measured experimental half lives of double beta decays. As a result, double beta decay provides a rigorous test of nuclear structure calculations. The double beta decay is already well established experimentally for a couple of isotopes. Because of the enormous energy release, the transition from the ground state 0⁺ of the initial to the ground state 0⁺ of the final nuclei is the most favored for experimental study of this unusual phenomenon. Transitions to the 2⁺ and 0⁺ excited states of the final nucleus have lately received more study [9] - [20].

Because many observable quantities are acquired from one-body operators, coefficients of fractional parentage (CFP) are quite useful. In addition, there is a resurgence of interest in experiments involving single nucleon transfer studies in the 1p shell region [21]. As an example, the odd-mass isotopes ^5,7He are particle-unbound [22]. In the framework of an R-matrix analysis, the data provide the neutron spectroscopic factor, which describes to what extent the total angular momentum quantum number $J^{\pi }=3/2^{-}$ ground state of ⁷He can be regarded as a ${}^6\text{He}\otimes 1p_{3/2}$ configuration, where the symbol $\otimes$ means tensor product of the single-particle Hilbert space. The spectral factor is a fundamental quantity that characterizes the single-particle nature of nuclear excitation and is hence consequently an important test of wave functions derived using newly established methods [23].

Heavy cluster emission and super-asymmetric fission have been studied theoretically since the late 1970s [24]. Cluster radioactivity is the phenomena of radioactive nuclei emitting particles heavier than the alpha particle spontaneously [25]. This process can be thought of as a case of strong asymmetric fission [26] or as a decay process involving cluster formation and tunneling through the barrier, analogous to alpha decay [27]. The penetrability of the pre-scission component of the barrier was proven to be similar to the super asymmetric fission model by interpreting the cluster preformation probability within a super asymmetric fission model as the penetrability of the pre-scission part of the barrier [28] [29] [30] [31]. Rose and Jones [32] validated this occurrence in the radioactive decay of ¹⁴C from ²²³Ra in 1984, and Alexandrov et al. [33] confirmed it a few months later. Since then, the ¹⁴C decay of many isotopes of Ra nuclei and many other heavy cluster decays have been observed [34].

With the advent of heavy ion beams for inducing multi-particle transfer reactions, a wealth of data from nuclei with a large variety of nucleon numbers has become available. The no-recoil, zero-range distorted wave Born approximation (DWBA) algorithms that are typically used to analyze light ion single-particle transfer data have been demonstrated to be insufficient in assessing the new studies. There have been some questions about how to properly extract structure information from less restrictive exact finite range DWBA computations [35]. Studying transfer amplitude between BCS-pair condensates and a proper technique for computing it is needed.

In this work, I define the asymmetry many-pair density matrix and drive the recursive relations of the asymmetry many-pair density matrix. I demonstrate that it can be calculated very quickly and stored within small computer memory consumption. I also drive the pair transfer amplitude for double beta decay and the cluster transfer amplitude for cluster decay expressed by asymmetry many-pair density matrix. The manuscript is organized as following. In Section 2, I briefly introduce the pair condensate with zero generalized seniority of the 2N-particle system and the many-pair density matrix. Then in Section 3, I drive the relation between the cluster transfer amplitude and the asymmetry many-pair density matrix. Finally, Section 4 summarizes the work.

2. Basic Expressions

In this section I briefly review the many-pair density matrix and the simple relation between the Pauli-blocked normalizations and many-pair density matrix. The pair-creation operator

$P_{\alpha }^{†}=a_{\alpha }^{†}{a}_{\tilde{\alpha }}^{†}=P_{\tilde{\alpha }}^{†}$ (1)

creates a pair of particles on the single-particle level $|\alpha \rangle$ and its time-reversed partner $|\tilde{\alpha }\rangle$ ( $|\widetilde{\stackrel{˜}{\alpha }}\rangle =-|\alpha \rangle$ ). The coherent pair-creation operator

$P^{†}=\frac{1}{2}{\displaystyle \underset{\alpha }{\sum }}{v}_{\alpha }{a}_{\alpha }^{†}{a}_{\tilde{\alpha }}^{†}={\displaystyle \underset{\alpha }{\sum }}{v}_{\alpha }{P}_{\alpha }^{†}$ (2)

creates a pair of particles with the real structure coefficients $v_{\alpha }$ that are distributed coherently over the whole single-particle space. The summation index in Equation (2) is a pair index that only accounts for half of the single-particle space (for instance, only account for those single-particle levels with a positive magnetic quantum number m). The state with zero generalized seniority of the 2N-particle system in the presence of pairing correlations is

$|{\varphi }_N\rangle =\frac{1}{\sqrt{{\chi }_N}}{\left(P^{†}\right)}^N|0\rangle \mathrm{,}$ (3)

where

${\chi }_N=\langle 0|P^N{\left(P^{†}\right)}^N|0\rangle$ (4)

is the normalization factor. The pair creation and annihilation operators, $P_{\alpha }^{†}$,

and $P_{\alpha }$ along with ${\hat{N}}_{\alpha }=\frac{1}{2}\left(a_{\alpha }^{†}{a}_{\alpha }+a_{\tilde{\alpha }}^{†}{a}_{\tilde{\alpha }}\right)$ are the generators of a quasi-spin

$su\left(2\right)$ algebra,

$\left[P_{\alpha }\mathrm{,}{P}_{\beta }^{†}\right]={\delta }_{\alpha \beta }\left(1-2{\hat{N}}_{\alpha }\right)\mathrm{,}$ (5)

$\left[{\hat{N}}_{\alpha }\mathrm{,}{P}_{\beta }^{†}\right]={\delta }_{\alpha \beta }{P}_{\alpha }^{†}\mathrm{.}$ (6)

The many-pair density matrix is introduced as

$t_{{\alpha }_1{\alpha }_2\cdots {\alpha }_p\mathrm{<mo>\; ;\; </mo>}{\beta }_1{\beta }_2\cdots {\beta }_q}^N\equiv \langle 0|P^{N-p}{P}_{{\alpha }_1}{P}_{{\alpha }_2}\cdots P_{{\alpha }_p}\times P_{{\beta }_1}^{†}{P}_{{\beta }_2}^{†}\cdots P_{{\beta }_q}^{†}{\left(P^{†}\right)}^{N-q}|0\rangle \mathrm{,}$ (7)

where all the indices ${\alpha }_1{\alpha }_2\cdots {\alpha }_p$, ${\beta }_1{\beta }_2\cdots {\beta }_q$ are distinct, owing to the Pauli principle. Note that the $t_{{\alpha }_1{\alpha }_2\cdots {\alpha }_p\mathrm{<mo>\; ;\; </mo>}{\beta }_1{\beta }_2\cdots {\beta }_q}^N$ is real because of the real $v_{\alpha }$, so $t_{{\alpha }_1{\alpha }_2\cdots {\alpha }_p;{\beta }_1{\beta }_2\cdots {\beta }_q}^N=t_{{\beta }_1{\beta }_2\cdots {\beta }_q;{\alpha }_1{\alpha }_2\cdots {\alpha }_p}^N$. The normalization ${\chi }_N$ defined in Equation (4) is the special case of Equation (7) when there is no $\alpha$ and $\beta$ index,

${\chi }_N\equiv t_{\mathrm{<mo>\; ;\; </mo>}}^N\mathrm{.}$ (8)

For convenience, I introduce $\left[{\gamma }_1{\gamma }_2\cdots {\gamma }_r\right]$ to represent a subspace of the original single-particle space, by removing Kramers pairs ${\gamma }_1{\gamma }_2\cdots {\gamma }_r$ of single-particle levels from the latter. The Pauli-blocked many-pair density matrix is defined as

$\begin{array}{c}{t}_{{\alpha }_1{\alpha }_2\cdots {\alpha }_p\mathrm{<mo>\; ;\; </mo>}{\beta }_1{\beta }_2\cdots {\beta }_q}^{\left[{\gamma }_1{\gamma }_2\cdots {\gamma }_r\right],N}\equiv \langle 0|P^{N-p}{P}_{{\gamma }_1}{P}_{{\gamma }_2}\cdots P_{{\gamma }_r}\times P_{{\alpha }_1}{P}_{{\alpha }_2}\cdots P_{{\alpha }_p}{P}_{{\beta }_1}^{†}{P}_{{\beta }_2}^{†}\cdots P_{{\beta }_q}^{†}\\ \times P_{{\gamma }_1}^{†}{P}_{{\gamma }_2}^{†}\cdots P_{{\gamma }_r}^{†}{\left(P^{†}\right)}^{N-q}|0\rangle ,\end{array}$ (9)

and Pauli-blocked normalization

${\chi }_N^{\left[{\gamma }_1{\gamma }_2\cdots {\gamma }_r\right]}\equiv t_{;}^{\left[{\gamma }_1{\gamma }_2\cdots {\gamma }_r\right],N}=\langle 0|P^N{P}_{{\gamma }_1}{P}_{{\gamma }_2}\cdots P_{{\gamma }_r}{P}_{{\gamma }_1}^{†}{P}_{{\gamma }_2}^{†}\cdots P_{{\gamma }_r}^{†}{\left(P^{†}\right)}^N|0\rangle$ (10)

could be easily derived. Also, there is no duplicated P operator, and duplicated $P^{†}$ operator in Equations (9) and (10), owing to the Pauli principle. The many-pair density matrix (7) and the Pauli-blocked many-pair density matrix (9) could be expressed by normalizations, and normalizations could be computed by recursive relations [36], which means for a certain nucleus, if we know all $v_{\alpha }$ (or $n_{\alpha }=\langle {\varphi }_N|{\hat{n}}_{\alpha }|{\varphi }_N\rangle =\langle {\varphi }_N|a_{\alpha }^{†}{a}_{\alpha }|{\varphi }_N\rangle$ ) of single-particle orbit, we could rapidly and accurately compute many-pair density matrix and normalizations [37].

3. Transfer Amplitude

In the Section 2, I introduce the many-pair density matrix, of which the bra and the ket are of a same nucleus ( $v_{\alpha }$ is the same). In this section I drive the cluster transfer amplitude to solve the issue with the inner product between two different nuclei. Defining

$|{{\varphi }^{\prime }}_N\rangle =\frac{1}{\sqrt{{{\chi }^{\prime }}_N}}{\left({P^{\prime }}^{†}\right)}^N|0\rangle \mathrm{,}$ (11)

where

${{\chi }^{\prime }}_N=\langle 0|{\left(P^{\prime }\right)}^N{\left({P^{\prime }}^{†}\right)}^N|0\rangle \mathrm{,}$ (12)

the inner product between $|{\varphi }_N\rangle$ and $|{{\varphi }^{\prime }}_N\rangle$ is

$\begin{array}{c}\langle {{\varphi }^{\prime }}_N|{\varphi }_N\rangle =\langle {\varphi }_N|{{\varphi }^{\prime }}_N\rangle \\ =\frac{1}{\sqrt{{{\chi }^{\prime }}_N{\chi }_N}}\langle 0|{\left(P^{\prime }\right)}^N{\left(P^{†}\right)}^N|0\rangle \\ =\frac{1}{\sqrt{{{\chi }^{\prime }}_N{\chi }_N}}\langle 0|{\left(P\right)}^N{\left({P^{\prime }}^{†}\right)}^N|0\rangle \\ =\frac{{\Upsilon }_N}{\sqrt{{{\chi }^{\prime }}_N{\chi }_N}}\mathrm{,}\end{array}$ (13)

where

$P^{\prime }={\displaystyle \underset{\alpha }{\sum }}{v^{\prime }}_{\alpha }{P}_{\alpha },$ (14)

and

${\Upsilon }_N=\langle 0|{\left(P^{\prime }\right)}^N{\left(P^{†}\right)}^N|0\rangle =\langle 0|{\left(P\right)}^N{\left({P^{\prime }}^{†}\right)}^N|0\rangle .$ (15)

The asymmetry many-pair density matrix is defined as

$u_{{\alpha }_1{\alpha }_2\cdots {\alpha }_p;{\beta }_1{\beta }_2\cdots {\beta }_q}^N\equiv \langle 0|{\left(P^{\prime }\right)}^{N-p}{P}_{{\alpha }_1}{P}_{{\alpha }_2}\cdots P_{{\alpha }_p}\times P_{{\beta }_1}^{†}{P}_{{\beta }_2}^{†}\cdots P_{{\beta }_q}^{†}{\left(P^{†}\right)}^{N-q}|0\rangle .$ (16)

Now I drive the relation between $u_{{\alpha }_1{\alpha }_2\cdots {\alpha }_p;{\beta }_1{\beta }_2\cdots {\beta }_q}^N$ and ${\Upsilon }_N$. By substituting Equation (14) into ${\left(P^{\prime }\right)}^{N-p}$ and polynomially expanding, terms with $P_{{\alpha }_1},P_{{\alpha }_2},\cdots$ or $P_{{\alpha }_p}$ vanish due to the Pauli principle, which is equivalent to Pauli block the ${\alpha }_1{\alpha }_2\cdots {\alpha }_p$ indices from ${\left(P^{\prime }\right)}^{N-p}$. Thus, ${\left(P^{\prime }\right)}^{N-p}$ could be replaced by

${\left({P^{\prime }}_{\left[{\alpha }_1{\alpha }_2\cdots {\alpha }_p\right]}\right)}^{N-p}$, where

${P^{\prime }}_{\left[{\alpha }_1{\alpha }_2\cdots {\alpha }_p\right]}\equiv P^{\prime }-{v^{\prime }}_{{\alpha }_1}{P}_{{\alpha }_1}-{v^{\prime }}_{{\alpha }_2}{P}_{{\alpha }_2}-\cdots -{v^{\prime }}_{{\alpha }_p}{P}_{{\alpha }_p}\mathrm{.}$ (17)

Next, use

$\begin{array}{c}{\left({P^{\prime }}_{\left[{\alpha }_1{\alpha }_2\cdots {\alpha }_p\right]}\right)}^{N-p}={\left({P^{\prime }}_{\left[{\alpha }_1{\alpha }_2\cdots {\alpha }_p{\beta }_1{\beta }_2\cdots {\beta }_q\right]}-{v^{\prime }}_{{\beta }_1}{P}_{{\beta }_1}-{v^{\prime }}_{{\beta }_2}{P}_{{\beta }_2}-\cdots -{v^{\prime }}_{{\beta }_q}{P}_{{\beta }_q}\right)}^{N-p}\\ =P_{N-p}^q{v^{\prime }}_{{\beta }_1}{P}_{{\beta }_1}\cdots {v^{\prime }}_{{\beta }_q}{P}_{{\beta }_q}{\left({P^{\prime }}_{\left[{\alpha }_1\cdots {\beta }_q\right]}\right)}^{N-p-q}+\text{Others},\end{array}$ (18)

where “Others” do not contribute and each contributing term must have the factor $P_{{\beta }_1}{P}_{{\beta }_2}\cdots P_{{\beta }_q}$ to annihilate $P_{{\beta }_1}^{†}{P}_{{\beta }_2}^{†}\cdots P_{{\beta }_q}^{†}$ resulting in the permutations $P_{N-p}^q$ in the last step.

Treating ${\left(P^{†}\right)}^{N-q}$ similarly, Equation (16) becomes

$\begin{array}{c}{u}_{{\alpha }_1{\alpha }_2\cdots {\alpha }_p;{\beta }_1{\beta }_2\cdots {\beta }_q}^N=P_{N-p}^q{P}_{N-q}^p{v}_{{\alpha }_1}{v}_{{\alpha }_2}\cdots v_{{\alpha }_p}{v^{\prime }}_{{\beta }_1}{v^{\prime }}_{{\beta }_2}\cdots {v^{\prime }}_{{\beta }_q}\\ \times \langle 0|{\left({P^{\prime }}_{\left[{\alpha }_1\cdots {\beta }_q\right]}\right)}^{N-p-q}{P}_{{\beta }_1}\cdots P_{{\beta }_q}{P}_{{\alpha }_1}\cdots P_{{\alpha }_p}\\ \times P_{{\alpha }_1}^{†}\cdots P_{{\alpha }_p}^{†}{P}_{{\beta }_1}^{†}\cdots P_{{\beta }_q}^{†}{\left(P_{\left[{\alpha }_1\cdots {\beta }_q\right]}^{†}\right)}^{N-p-q}|0\rangle \\ =\frac{\left(N-p\right)!\left(N-q\right)!}{{\left[\left(N-p-q\right)!\right]}^2}\times v_{{\alpha }_1}{v}_{{\alpha }_2}\cdots v_{{\alpha }_p}{v^{\prime }}_{{\beta }_1}{v^{\prime }}_{{\beta }_2}\cdots {v^{\prime }}_{{\beta }_q}\\ \times {\Upsilon }_{N-p-q}^{\left[{\alpha }_1{\alpha }_2\cdots {\alpha }_p{\beta }_1{\beta }_2\cdots {\beta }_q\right]},\end{array}$ (19)

where

$\begin{array}{c}{\Upsilon }_{N-p-q}^{\left[{\alpha }_1{\alpha }_2\cdots {\alpha }_p{\beta }_1{\beta }_2\cdots {\beta }_q\right]}=\langle 0|{\left({P^{\prime }}_{\left[{\alpha }_1\cdots {\beta }_q\right]}\right)}^{N-p-q}{P}_{{\beta }_1}\cdots P_{{\beta }_q}{P}_{{\alpha }_1}\cdots P_{{\alpha }_p}\\ \times P_{{\alpha }_1}^{†}\cdots P_{{\alpha }_p}^{†}{P}_{{\beta }_1}^{†}\cdots P_{{\beta }_q}^{†}{\left(P_{\left[{\alpha }_1\cdots {\beta }_q\right]}^{†}\right)}^{N-p-q}|0\rangle \\ =\langle 0|{\left(P^{\prime }\right)}^{N-p-q}{P}_{{\beta }_1}\cdots P_{{\beta }_q}{P}_{{\alpha }_1}\cdots P_{{\alpha }_p}\\ \times P_{{\alpha }_1}^{†}\cdots P_{{\alpha }_p}^{†}{P}_{{\beta }_1}^{†}\cdots P_{{\beta }_q}^{†}{\left(P^{†}\right)}^{N-p-q}|0\rangle .\end{array}$ (20)

From Equation (2), (5) and (6) it is easy to derive the identity

$P_{\alpha }{\left(P^{†}\right)}^{N+1}|0\rangle =v_{\alpha }\left(N+1\right){\left(P^{†}\right)}^N|0\rangle -{\left(v_{\alpha }\right)}^{2}N\left(N+1\right)P_{\alpha }^{†}{\left(P^{†}\right)}^{N-1}|0\rangle \mathrm{.}$ (21)

Now premultiplication with $\langle 0|{\left(P^{\prime }\right)}^N$ gives the result

$u_{\alpha ;}^{N+1}=v_{\alpha }\left(N+1\right){\Upsilon }_N-{\left(v_{\alpha }\right)}^{2}N\left(N+1\right)u_{;\alpha }^N.$ (22)

Based on Equations (14) and (15), ${\Upsilon }_N$ could be expressed as

${\Upsilon }_N={\displaystyle \underset{\alpha }{\sum }}{v^{\prime }}_{\alpha }\langle 0|{\left(P^{\prime }\right)}^{N-1}{P}_{\alpha }{\left(P^{†}\right)}^N|0\rangle ={\displaystyle \underset{\alpha }{\sum }}{v^{\prime }}_{\alpha }{u}_{\alpha ;}^N.$ (23)

Substituting $u_{\alpha ;}^{N+1}=v_{\alpha }\left(N+1\right){\Upsilon }_N^{\left[\alpha \right]}$ and $u_{;\alpha }^N={v^{\prime }}_{\alpha }N{\Upsilon }_{N-1}^{\left[\alpha \right]}$ implied from Equation (19) into (22) and (23), I get the recursive relations

${\Upsilon }_N=N{\displaystyle \underset{\alpha }{\sum }}{v}_{\alpha }{v^{\prime }}_{\alpha }{\Upsilon }_{N-1}^{\left[\alpha \right]},$ (24)

${\Upsilon }_N-{\Upsilon }_N^{\left[\alpha \right]}=N^2{v}_{\alpha }{v^{\prime }}_{\alpha }{\Upsilon }_{N-1}^{\left[\alpha \right]},$ (25)

with initial value ${\Upsilon }_0=\langle 0|0\rangle =1$ and ${\Upsilon }_0^{\left[\alpha \right]}=\langle 0|P_{\alpha }{P}_{\alpha }^{†}|0\rangle =1$. This is because by definition the norm of all the orthonormal basis for the Hilbert space is 1 in quantum mechanics, with no exceptions for vacuum. Knowing ${\Upsilon }_0^{\left[\alpha \right]}$, I could compute ${\Upsilon }_1$ by Equation (24), and then ${\Upsilon }_1^{\left[\alpha \right]}$ by Equation (25). Equations (24) and (25) are also valid in the blocked subspace $\left[{\gamma }_1{\gamma }_2\cdots {\gamma }_r\right]$. For example, in $\left[\beta \right]$ they read

${\Upsilon }_N^{\left[\beta \right]}=N{\displaystyle \underset{\alpha }{\sum }}{v}_{\alpha }{v^{\prime }}_{\alpha }{\Upsilon }_{N-1}^{\left[\alpha \beta \right]},$ (26)

${\Upsilon }_N^{\left[\beta \right]}-{\Upsilon }_N^{\left[\alpha \beta \right]}=N^2{v}_{\alpha }{v^{\prime }}_{\alpha }{\Upsilon }_{N-1}^{\left[\alpha \beta \right]},$ (27)

with initial value ${\Upsilon }_0^{\left[\alpha \beta \right]}=\langle 0|P_{\alpha }{P}_{\beta }{P}_{\beta }^{†}{P}_{\alpha }^{†}|0\rangle =1$.

So by definition, the spectroscopic factor between an odd nucleus $a_{\alpha }^{†}|{{\varphi }^{\prime }}_{N-1}\rangle$ and an even nucleus $|{\varphi }_N\rangle$ is

$\begin{array}{c}{\kappa }_{\alpha {\alpha }_1}=\langle {{\varphi }^{\prime }}_{N-1}|a_{\alpha }{a}_{{\alpha }_1}|{\varphi }_N\rangle \\ =\frac{{\delta }_{\tilde{\alpha }{\alpha }_1}}{\sqrt{{{\chi }^{\prime }}_{N-1}{\chi }_N}}\langle 0|{\left(P^{\prime }\right)}^{N-1}{P}_{\alpha }{\left(P^{†}\right)}^N|0\rangle \\ =\frac{{\delta }_{\tilde{\alpha }{\alpha }_1}}{\sqrt{{{\chi }^{\prime }}_{N-1}{\chi }_N}}{u}_{\alpha ;}^N\\ ={\delta }_{\tilde{\alpha }{\alpha }_1}\frac{N{v}_{\alpha }{\Upsilon }_{N-1}^{\left[\alpha \right]}}{\sqrt{{{\chi }^{\prime }}_{N-1}{\chi }_N}},\end{array}$ (28)

the pair transfer amplitude is

${\kappa }_{\tilde{\alpha }\alpha }=\langle {{\varphi }^{\prime }}_{N-1}|P_{\alpha }|{\varphi }_N\rangle =\frac{N{v}_{\alpha }{\Upsilon }_{N-1}^{\left[\alpha \right]}}{\sqrt{{{\chi }^{\prime }}_{N-1}{\chi }_N}},$ (29)

and the cluster transfer amplitude is

$\begin{array}{c}{\kappa }_{\alpha {\alpha }_1{\tilde{\alpha }}_2{\alpha }_2\cdots {\tilde{\alpha }}_p{\alpha }_p}={\delta }_{\tilde{\alpha }{\alpha }_1}\langle {{\varphi }^{\prime }}_{N-1}|P_{{\alpha }_1}{P}_{{\alpha }_2}\cdots P_{{\alpha }_p}|{\varphi }_N\rangle \\ =\frac{{\delta }_{\tilde{\alpha }{\alpha }_1}}{\sqrt{{{\chi }^{\prime }}_{N-p}{\chi }_N}}\langle 0|{\left(P^{\prime }\right)}^{N-p}{P}_{{\alpha }_1}\cdots P_{{\alpha }_p}{\left(P^{†}\right)}^N|0\rangle \\ =\frac{{\delta }_{\tilde{\alpha }{\alpha }_1}}{\sqrt{{{\chi }^{\prime }}_{N-p}{\chi }_N}}{u}_{{\alpha }_1{\alpha }_2\cdots {\alpha }_p;}^N\\ ={\delta }_{\tilde{\alpha }{\alpha }_1}\frac{N!v_{{\alpha }_1}{v}_{{\alpha }_2}\cdots v_{{\alpha }_p}{\Upsilon }_{N-p}^{\left[{\alpha }_1{\alpha }_2\cdots {\alpha }_p\right]}}{\left(N-p\right)!\sqrt{{{\chi }^{\prime }}_{N-p}{\chi }_N}}.\end{array}$ (30)

The Equation (28), (29) and (30) can be used in computer program to compute any kinds of nuclear transfer amplitude between two given nuclei, since all parameters [the paired particle number N, the number of paired particle to transfer p, the structure coefficients $v_{\alpha }$ and ${v^{\prime }}_{\alpha }$, the normalization factors ${\chi }_N$ and ${{\chi }^{\prime }}_N$, and the transfer factors ${\Upsilon }_N$ which can be compute by recursive relations (24) and (25)] are known.

4. Summary

In summary, I propose a new algorithm that uses recursive relations to compute spectroscopic factor, pair transfer amplitude and cluster transfer amplitude. This work starts from the inner product between two different nuclei. In order to compute the product, I define the asymmetry many-pair density matrix and drive its recursive relations. The recursive relations can be calculated very quickly and stored within small computer memory consumption, which is the foundation of performing the next step. Due to the properties of the recursive relations, it could be easy to compute the transfer amplitude such as the spectroscopic factor, the pair transfer amplitude and the cluster transfer amplitude using the asymmetry many-pair density matrix.

Acknowledgements

The author gratefully acknowledges discussions with Prof. L. Y. Jia.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1]	Bardeen, J., Cooper, L.N. and Schrieffer, J.R. (1957) Microscopic Theory of Superconductivity. Physical Review Journals Archive, 106, 162-164. https://doi.org/10.1103/PhysRev.106.162
[2]	Bardeen, J., Cooper, L.N. and Schrieffer, J.R. (1957) Theory of Superconductivity. Physical Review Journals Archive, 108, 1175-1204. https://doi.org/10.1103/PhysRev.108.1175
[3]	Bohr, A., Mottelson, B.R. and Pines, D. (1958) Possible Analogy between the Excitation Spectra of Nuclei and Those of the Superconducting Metallic State. Physical Review Journals Archive, 110, 936-938. https://doi.org/10.1103/PhysRev.110.936
[4]	Belyaev, S.T. (1959) Effect of Pairing Correlations on Nuclear Properties. Det Kongelige Danske Videnskabernes Selskab. Matematisk-fysiske Meddelelser, 31, 1-55.
[5]	Broglia, R.A. and Zelevinsky, V. (2013) Fifty Years of Nuclear BCS: Pairing in Finite Systems. World Scientific, Singapore. https://doi.org/10.1142/8526
[6]	Dietrich, K., Mang, H.J. and Pradal, J.H. (1964) Conservation of Particle Number in the Nuclear Pairing Model. Physical Review Journals Archive, 135, B22-B34. https://doi.org/10.1103/PhysRev.135.B22
[7]	Jia, L.Y. (2017) Generalized Seniority on a Deformed Single-Particle Basis. Physical Review C, 96, Article ID: 034313. https://doi.org/10.1103/PhysRevC.96.034313
[8]	Goeppert-Mayer, M. (1935) Double Beta-Disintegration. Physical Review Journals Archive, 48, 512-516. https://doi.org/10.1103/PhysRev.48.512
[9]	Barabash, A.S. (1990) A Possibility for Experimentally Observing Two Neutrino Double Beta Decay. Soviet Journal of Experimental and Theoretical Physics Letters, 51, 207.
[10]	Barabash, A., Kopylov, A. and Cherehovsky, V. (1990) Search for Double β-Decay of 100Mo and 116Cd to the Excited States of 100Ru and 116Sn. Physics Letters B, 249, 186-190. https://doi.org/10.1016/0370-2693(90)91240-C
[11]	Blum, D., Busto, J., Campagne, J., Dassié, D., Hubert, F., Hubert, P., Isaac, M., Izac, C., Jullian, S., Kouts, B., Kropivyansky, B., Lalanne, D., Lamhamdi, T., Laplanche, F., Leccia, F., Linck, I., Longuemare, C., Mennrath, P., Natchez, F., Scheibling, F., Szklarz, G., Tretyak, V. and Zdesenko, Y. (1992) Search for γ-Rays Following ββ Decay of 100Mo to Excited States of 100Ru. Physics Letters B, 275, 506-511. https://doi.org/10.1016/0370-2693(92)91624-I
[12]	Piepke, A., Beck, M., Bockholt, J., Glatting, D., Heusser, G., Klapdor-Kleingrothaus, H., Maier, B., Petry, F., Schmidt-Rohr, U., Strecker, H., Völlinger, M., Barabash, A., Umatov, V., Müller, A. and Suhonen, J. (1994) Investigation of the ββ Decay of 116Cd into Excited States of 116Sn. Nuclear Physics A, 577, 493-510. https://doi.org/10.1016/0375-9474(94)90930-X
[13]	Barabash, A., Avignone, F., Collar, J., Guerard, C., Arthur, R., Brodzinski, R., Miley, H., Reeves, J., Meier, J., Ruddick, K. and Umatov, V. (1995) Two Neutrino Double-Beta Decay of 100Mo to the First Excited 0+ State in 100Ru. Physics Letters B, 345, 408-413. https://doi.org/10.1016/0370-2693(94)01657-X
[14]	Suhonen, J. and Civitarese, O. (1993) Two-Neutrino ββ Decay to Excited States. The 0+→2+ Decay of 136Xe. Physics Letters B, 308, 212-215. https://doi.org/10.1016/0370-2693(93)91273-P
[15]	Suhonen, J. and Civitarese, O. (1994) Quasiparticle Random Phase Approximation Analysis of the Double Beta Decay of 100Mo to the Ground State and Excited States of 100Ru. Physical Review C, 49, 3055-3060. https://doi.org/10.1103/PhysRevC.49.3055
[16]	Stoica, S. (1994) Half-Lives for Two Neutrino Double-Beta-Decay Transitions to First 2+ Excited States. Physical Review C, 49, 2240-2243. https://doi.org/10.1103/PhysRevC.49.2240
[17]	Dhiman, S.K. and Raina, P.K. (1994) Two-Neutrino Double-Beta Decay Matrix Elements for Ground and Excited States of 76Ge and 82Se Nuclei. Physical Review C, 50, R2660-R2663. https://doi.org/10.1103/PhysRevC.50.R2660
[18]	Schwieger, J., Simkovic, F., Faessler, A. and Kaminski, W.A. (1997) Double Beta Decay to Excited States in 76Ge within Renormalized QRPA. Journal of Physics G: Nuclear and Particle Physics, 23, 1647-1653. https://doi.org/10.1088/0954-3899/23/11/012
[19]	Suhonen, J., Toivanen, J., Barabash, A.S., Vanushin, I.A., Umatov, V.I., Gurriarán, R., Hubert, F. and Hubert, P. (1997) Renormalized Proton-Neutron QRPA and Double Beta Decay of 82Se to Excited States in 82Kr. Zeitschrift für Physik A Hadrons and Nuclei, 358, 297-301. https://doi.org/10.1007/s002180050333
[20]	Faessler, A. and Simkovic, F. (1998) Double Beta Decay. Journal of Physics G: Nuclear and Particle Physics, 24, 2139-2178. https://doi.org/10.1088/0954-3899/24/12/001
[21]	Cohen, S. and Kurath, D. (1967) Spectroscopic Factors for the 1p Shell. Nuclear Physics A, 101, 1-16. https://doi.org/10.1016/0375-9474(67)90285-0
[22]	Jonson, B. (2004) Light Dripline Nuclei. Physics Reports, 389, 1-59. https://doi.org/10.1016/j.physrep.2003.07.004
[23]	Beck, F., Frekers, D., von Neumann-Cosel, P., Richter, A., Ryezayeva, N. and Thompson, I. (2007) Spectroscopic Factor of the 7He Ground State. Physics Letters B, 645, 128-132. https://doi.org/10.1016/j.physletb.2006.12.036
[24]	Sandulescu, A., Poenaru, D.N. and Greiner, W. (1980) New Type of Decay of Heavy Nuclei Intermediate between Fission and Cap Alpha Decay. Soviet Journal of Particles and Nuclei, 11, 528-541.
[25]	Sandulescu, A., Florescu, A. and Greiner, W. (1989) Cold Fission as Emission of Fragments. Journal of Physics G: Nuclear and Particle Physics, 15, 1815. https://doi.org/10.1088/0954-3899/15/12/008
[26]	Poenaru, D.N., Ivascu, M. and Sandulescu, A. (1979) Alpha Decay as a Fission-Like Process. Journal of Physics G: Nuclear Physics, 5, L169-L173. https://doi.org/10.1088/0305-4616/5/10/005
[27]	Blendowske, R., Fliessbach, T. and Walliser, H. (1987) Microscopic Calculation of the 14C Decay of Ra Nuclei. Nuclear Physics A, 464, 75-89. https://doi.org/10.1016/0375-9474(87)90423-4
[28]	Poenaru, D.N. and Greiner, W. (1991) Cluster Preformation as Barrier Penetrability. Physica Scripta, 44, 427-429. https://doi.org/10.1088/0031-8949/44/5/004
[29]	Poenaru, D.N. and Greiner, W. (1991) Rare Decay Modes by Cluster Emission from Nuclei. Journal of Physics G: Nuclear and Particle Physics, 17, S443. https://doi.org/10.1088/0954-3899/17/S/045
[30]	Poenaru, D.N., Greiner, W. and Hourani, E. (1995) 12C Emission from 114Ba and Nuclear Properties. Physical Review C, 51, 594-600.
[31]	Santhosh, K.P. and Joseph, A. (2002) Exotic Decay: Transition from Cluster Mode to Fission Mode. Pramana, 59, 599-609. https://doi.org/10.1007/s12043-002-0071-y
[32]	Rose, H.J. and Jones, G.A. (1984) A New Kind of Natural Radioactivity. Nature, 307, 245-247. https://doi.org/10.1038/307245a0
[33]	Alexandrov, D., Belyatsky, A., Glukhov, Y.A., Novatsky, B., Oglobin, A., Stepanov, D., et al. (1984) Cluster Radioactivity in 127I. JETP Letters, 40, 909.
[34]	Santhosh, K.P. and Biju, R.K. (2008) Alpha Decay, Cluster Decay and Spontaneous Fission in 294-326122 Isotopes. Journal of Physics G: Nuclear and Particle Physics, 36, Article ID: 015107. https://doi.org/10.1088/0954-3899/36/1/015107
[35]	Draayer, J.P. (1975) Alpha-Particle Spectroscopic Amplitudes for SD Shell Nuclei. Nuclear Physics A, 237, 157-181. https://doi.org/10.1016/0375-9474(75)90470-4
[36]	Jia, L.Y. (2015) Practical Calculation Scheme for Generalized Seniority. Journal of Physics G: Nuclear and Particle Physics, 42, Article ID: 115105. https://doi.org/10.1088/0954-3899/42/11/115105
[37]	Jia, L.Y. (2019) Application of the Variational Principle to a Coherentpair Condensate: The BCS Case. Physical Review C, 99, Article ID: 014302. https://doi.org/10.1103/PhysRevC.99.014302