-
The nuclear mass (binding energy) is one of the most important topics in nuclear physics and astrophysics. In the semi-empirical mass formula based on the liquid drop model, the binding energy
B(A,Z) of a nucleus is the sum of the volume, surface, Coulomb, symmetry, and pairing terms, which can be expressed as a function of the mass number A and charge number Z [1,2],B(A,Z)=avA−asA2/3−aCZ2A1/3−asym(N−Z)2A+apA−1/2,
(1) where
av ,as ,aC ,asym , andap are the volume, surface, Coulomb, symmetry, and pairing energy coefficients, respectively. The coefficients are obtained by fitting to experimentally measured binding energies of nuclei. In contrast, the nuclear mass is usually indicated by the mass excess (nuclide mass minus the mass number in atomic mass units). The mass excess of the isobaric nuclei can be expressed as a parabola, as a function of the atomic number Z with a high degree of precision, referred to as the Bohr-Wheeler parabola [3]. In Ref. [4], the nuclide mass equation is expressed asM(A,Z)=M(A,ZA)+12bA(Z−ZA)2+PA−S(N,Z),
(2) where M is the mass excess of nucleus (A,Z) in MeV, and Z and A are the atomic and mass numbers of the nucleus.
M(A,ZA) is described as the minimum of the mass excess parabola, whereZA is a parameter and not the atomic number of an existing nuclei.bA is the curvature parameter of the isobaric mass parabola,PA is the pairing energy due to the odd-even variation, andS(N,Z) is the shell correction term. In the textbook [5], the relation between the nuclear mass excess and the binding energy is written asB(A,Z)=931.4943×(0.008665A−0.00084Z)−M(A,Z),
(3) where 0.008665 is the mass excess of the neutron, and 0.00084 is the hydrogen-neutron mass difference in atomic mass units. One atomic mass unit is equivalent to 931.4943 MeV. The question is, hence, whether the Coulomb energy coefficient and the symmetry energy coefficient in the binding energy formula can be expressed by parabolic parameters
ZA andbA .There are two ways to obtain the parameters
ZA andbA . One is by fitting the experimentally measured mass excess of isobaric nuclei with Eq.(2). The other is by fitting by measured β-decay energies of isobaric nuclei with Eq.(4). The first method requires that four parametersZA ,bA ,M(A,ZA) , andPA are determined by a parabola fit, while the second requires that three parametersZA ,bA , andPA are determined by fitting theβ -decay energies. Isobaric analyses ofβ -decay energies were conducted by several authors [3,6,7], where the early work was based on scanty experimental data. In Ref. [3], Bohr and Wheeler employed the least-squares method to determine the parameters of the mass parabola,ZA ,bA , andPA , for 20 A values by a parabola fit of the isobaric mass excess. In Ref. [6], J. W. Dewdney analyzed the isobaricβ -decay energies and used the least-squares method to determine the parameters of the mass parabola for 157 A values. The experimental data were taken from Refs. [8,9]. Later, X.Y. Li and co-workers [7] updated the three parameters of the isobaric mass parabola for 234 A values (10⩽A⩽253 ) in the same manner as Dewdney, however they adopted a different mass table AME1977 [10] including about 1000 nuclides. In the recent decades, with the development of scientific instruments and the progress of technology, a large number of unstable nuclei could be produced, and their masses could be measured with high precision. The experimental information or recommended values for nuclear and decay properties of 3437 nuclides are compiled in the mass table AME2016 [11]. The number and precision of nuclear masses in AME2016 have been considerably increased in comparison with the results compiled in the AME1977 database [10]. The available data sets are presently significantly larger, and it is possible to analyze them statistically.Double
β -decay is a popular topic. It is a rare transition between two nuclei with the same mass number, which changes the nuclear charge number by two units,(A,Z)→(A,Z±2) , and has been observed in many nuclei [12,13]. The energy released in the doubleβ -decay exhibits even a greater regularity than the singleβ -decay, as shown in Ref. [14]. However, detailed studies of the mass parabola have not been performed using the total doubleβ -decay energies. In this study, we propose a simple empirical formula to obtain the parametersZA andbA . The relation between the Coulomb energy coefficient and the symmetry energy coefficient in terms ofZA andbA is then presented. We use more than 2400 total doubleβ -decay energies, which are complied in the AME2016 database, to analyze the properties of the Bohr-Wheeler isobaric mass parabola in all modes of doubleβ -decay employing a theoretical method. The advantage of this approach is that there are only two parameters in the expression for the doubleβ -decay energies, as the pairing energyPA and the mass excessM(A,ZA) are removed by the mass difference between(A,Z) and(A,Z±2) . Furthermore, the overall simplicity of the doubleβ -decay energy pattern may point to a convenient empirical mass formula. -
The Q value is defined as the total energy released in a given nuclear reaction. The Q value of
β -decay is calculated by the mass difference between the two nuclei(A,Z) and(A,Z+1) , while the expressions forβ -decay energies may be derived from the empirical mass Eq. (2), when the shell correction is neglected. According to Eq.(2), theβ -decay energy is directly proportional toZ−ZA . For odd A isobars, there is only one isobaric mass parabola, since the pairing energy term is equal to zero, and theβ -decay energies plotted versus Z lie on a single straight line. However, for even A isobars, due to the extra stability associated with pairs of similar nucleons, two isobaric mass parabolas with the same shape are obtained, but they are displaced one below the other. The lower parabola contains even-even nuclides, while odd-odd nuclides lie on the upper parabola. Consequently, theβ -decay energies plotted versus Z do not lie on a single straight line, but on a pair of parallel straight lines. For negativeβ -decay and positiveβ -decay (or electron capture), one can combine them into one universal expressionQβ=M(A,Z)−M(A,Z+1)=−bA(Z−ZA+12)+ΔE.
(4) The reaction will proceed only when the Q value is positive: "
Qβ>0 " forβ− -decay occurs when the mass of atomM(A,Z) is greater than the mass of atomM(A,Z+1) ; "Qβ<0 " forβ+ -decay occurs when the mass of atomM(A,Z+1) is greater than that ofM(A,Z) . For even-even nuclidesΔE=2PA , for odd-odd nuclidesΔE=−2PA , and for odd-A nuclides (i.e., even-odd and odd-even)ΔE=0 .Our aim in this study is to propose a very simple empirical formula that only depends on the basic mass parabolic parameters
ZA andbA . We employ the total doubleβ -decay energy to analyze the mass parabola rather than the total singleβ -decay energy. The parametersM(A,ZA) andΔE are removed by the mass difference between(A,Z) and(A,Z+2) . We obtain a universal expression for the total doubleβ -decay energy in the same manner as above,Q2β=M(A,Z)−M(A,Z+2)=−2bA(Z−ZA+1),
(5) where “
Q2β>0 ” forβ−β− , and “Q2β<0 ” forβ+β+ . The nucleus(Z,A) will be called the decay nucleus, no matter whether the decay actually proceeds from(A,Z) to(A,Z+2) or vice versa. The total doubleβ -decay energies lie on the same lines for odd-odd and even-even nuclei for even-A values, and for adjacent odd-A and even-A values, they lie approximately on the same lines. -
We employ the three Eqs. (2), (4), and (5) to analyze the corresponding experimental data, i.e., the mass excess, total
β -decay energies, and total doubleβ -decay energies. The parametersZA andbA in the isobaric mass parabola are obtained using the least-squares fitting procedure. The calculated results are compared, and they are found to be almost identical for three different methods of fitting the corresponding experimental data.The first equation Eq.(2) requires at least three parameters,
ZA ,bA , andM(A,ZA) , to fit the experimental mass excess for odd-A nuclei. Here,PA = 0 is assumed, based on the values ofPA<0.3 MeV for odd-A nuclei, presented in Ref. [6]. For even-A nuclei, an isobaric slice through the mass surface tends to yield two parabolas of the same shape, but they are displaced one below the other. Even nuclides lie on the lower parabola, while the odd nuclides lie on the upper parabola. Hence, it is necessary to unify the parametersZA andbA for even-even and odd-odd nuclei,bA=0.5(bo−oA+be−eA) andZA=0.5(Zo−oA+Ze−eA) , wherebe−eA andZe−eA are for even-even nuclei, andbo−oA andZo−oA for odd-odd nuclei. Finally, the value ofΔE is obtained from the differenceΔE=Mo−o(A,ZA)−Me−e(A,ZA) . For Eq. (4), the procedure is similar to that for Eq. (2), and the value ofΔE equals to half of the difference between the two intercepts of a pair of parallel lines.To construct the formula for
Q2β , we start by plotting the total doubleβ -decay energyQ2β versus the charge number Z for 11 odd-A nuclei7⩽A⩽257 in Fig. 1. The experimental data forQ2β are taken from AME2016 [11]. The straight line relations in this plot are clearly distinguishable, and the slopes2bA of these lines exhibit the general trend thatQ2β decreases slightly with increasing mass number A. A similar linear relation is also shown for even-A values in Fig. 2.Figure 1. (color online) The total double
β -decay energyQ2β as a function of charge number Z for 11 odd-A values (7⩽A⩽257 ). The experimental data are taken from AME2016 [11]Figure 2. (color online) Two examples of the experimental data and their treatment. Mass excess (a),
β -decay energy (b) as a function of the charge number Z, and the fitted curves forA=131 . (c) and (d) are the same as (a) and (b) forA=132 . Solid curves in Fig. 2 (a) and (c) denote the results of fitting the experimental mass excess with Eq. (2). Solid lines in Fig. 2 (b) and (d) denote the results of fitting the experimental dataQβ with Eq. (4), and dashed lines in Fig. 2 (b) and (d) denote the results of fitting the experimental dataQ2β with Eq. (5).Figure 2 shows two examples of the treatment of experimental data for
A=131 andA=132 as a function of the nuclear charge number Z. Figure 2 (a) is a plot of the experimental isobaric mass parabola for A = 131. The mass parabola is related to theβ and2β -decay energies shown in the plot below (b). The interceptZA on theβ -decay energy plot corresponds to the minimum of the mass parabola. The slopebA of theβ -decay energy plot is a measure of the steepness of the mass parabola. The red dashed lines in Fig. 2 (b) and (d) denote the results of fitting the experimental data forQ2β by Eq. (5). Figure 2 (c) and (d) are similar plots for A = 132. Figure 2 (c) shows that the mass excess lies alternately on a pair of parabolas of identical shape. The parabola for even masses is below the odd-mass parabola byΔE . This implies that theβ -decay points lie on a pair of parallel straight lines, while the doubleβ -decay points lie on a single straight line for both odd-A and even-A.Figure 3 presents the mass parabolic parameters
ZA (a),bA (b), andΔE (c) as a function of the nuclear mass number A. The values ofZA andbA obtained by the aforementioned procedure are shown in Fig. 3 (a) and (b) with blue solid curves using Eq. (5).ΔE is depicted with blue open circles in Fig. 3 (c) calculated with Eq.(4). The solid curveΔE=2PA is the result of a fit of theβ -decay energy, and the pairing energy termPA=11.50δ/A1/2 is obtained, whereδ equals to 0 for odd-A, +1 for even-even and -1 for odd-odd nuclides. The calculated values are compared with the early results of Bohr and Wheeler [3] (solid squares), Dewdney [6] (solid circles) and X.Y.Li [7] (solid triangles). There are currently 358 known β-decay stable nuclides [15], plotted in Fig. 3 (a) with open circles. We found a very good agreement between our results and those obtained previously in Refs. [6,7], except for a few A values (ZA parameter for A = 111, andbA parameter for A = 144). We presume that our results for these cases are more robust, because the discontinuities occur only for the two A values shown in Fig. 3. Our results indicate thatZA = 47.7566 andbA = 1.7955 for A = 111, andZA = 59.9768 andbA = 1.5102 for A = 144.Figure 3. (color online) The mass parabola parameters
ZA (a),bA (b), andΔE (c) as a function of the nuclear mass number A. The calculated values (blue solid curves) are compared with early results of Bohr and Wheeler [3] (solid squares), Dewdney [6] (solid circles) and X.Y. Li [7] (solid triangles). Black solid curve in (c) is the pairing energy termΔE=2PA , which is obtained by fitting theβ -decay energies. -
The Coulomb energy coefficient
aC can be determined from the mass relation for mirror nuclei. On the one hand, due to the charge independence of the nuclear force, the binding energies of mirror pairs differ only in their Coulomb energy. The Coulomb energy difference between a pair of mirror nuclei is proportional toY=N−Z [16,17],−ΔB=EC[A,12(A+Y)]−EC[A,12(A−Y)]=bcY,
(6) where
bc is the proportionality coefficient, which indeed depends on A.On the other hand, the mass relations for mirror nuclei can be obtained from Eq. (2). First, one needs to express the charge number Z with the mass number A, i.e.,
Z=12(A±Y) for a pair of mirror nuclei. Then, the difference of mass excess for a pair of mirror nuclei for both odd and even A with differentY=1,2,3,... readsΔM=M[A,12(A+Y)]−M[A,12(A−Y)]=12bA(A−2ZA)Y.
(7) Combining the three Eqs. (3), (6), and (7), and taking the hydrogen-neutron mass difference
Δ(1H−n)= −0.7825 MeV, we obtainΔMY=bc=12bA(A−2ZA)−Δ(1H−n).
(8) This result is consistent with those in Refs. [18,19]. Our results are presented in Fig. 4, where the mass difference of 95 pairs of mirror nuclei, scaled by the charge difference Y, are plotted against
A2/3 , and are seen to lie on a straight line. The value of Y ranges from 1 (32 cases) to 6 (1 case). The solid and dashed curves denote the results with and without the shell correction energies, respectively. The shell correction energies are taken from Ref. [20]. The thick solid curve shows some oscillations and fluctuations due to the shell effect. When the shell corrections are taken into account, the fluctuations of the extractedbc are effectively reduced. The thin solid line represents the straight line fit. We obtainbc=0.691A2/3−0.8724 MeV by fitting the left hand of Eq. (8) for 251 A values, with an rms deviation of 0.384 MeV, which is the error associated with this method of determiningaC and c , and it is smaller than in Ref. [18]. The intercept at 0.8724 is the contribution of the Coulomb exchange term and the nuclear surface diffuseness correction term. If we assume the Coulomb energy expressionEC=aCZ2A1/3(1−cZ−2/3) , in which the contribution of the Coulomb exchange term is taken into account, we obtainaC=0.691 MeV andc=1.1914 MeV. The mirror-nuclide method can be extended to include all values of A in Eq. (8), although mirror nuclides are observed forA⩽75 . The availability of the Coulomb energy coefficients for the complete range of A values should be helpful in studies addressing the variation of the Coulomb energy with the mass number.Figure 4. (color online) Scaled mass difference (solid squares) of 95 pairs of mirror nuclei in the region 11
⩽A⩽ 75 as a function ofA2/3 . Solid and dashed curves represent the left-hand side of Eq.(8) with and without the shell correction energies, respectively, which are taken from Ref. [20]. Thin solid line represents the straight line fit. -
The symmetry energy coefficient plays an important role in nuclear physics and astrophysics. In this study, the symmetry energy coefficient
asym is deduced and expressed by the two parametersZA andbA . First, we insert Eqs. (1) and (3) into the expression (5) and take the Q value of the doubleβ -decay of the mother nuclide (A,Z-1), to obtainQ2β(A,Z−1)=M(A,Z−1)−M(A,Z+1)=B(A,Z+1)−B(A,Z−1)+2Δ(1H−n)=−(16asymA+4acA1/3)Z+8asym+1.5649.
(9) The Q value of the double
β -decay of (A,Z-1) in Eq. (5) is written asQ2β(A,Z−1)=−2bA(Z−ZA),
(10) then, the symmetry energy coefficient
asym is obtained by solving the combination of Eqs. (9) and (10), that isasym=(2bAZA−1.5649)8≃bAZA4.
(11) There are two forms for describing the mass-dependent symmetry coefficient. One was proposed by Danielewicz and Lee,
asym(A)=S0/(1+κA−1/3) [21], whereS0 is the volume-symmetry energy coefficient of the nucleus, andκ is the ratio of the surface-symmetry coefficient to the volume-symmetry coefficient. The other form is frequently employed, and it is written asasym(A)= S0(1−κA−1/3) [22–27]. Fig. 5 shows the symmetry energy coefficientasym(A) as a function of A from Eq. (11). The open circles denote the results without the shell corrections. These exhibit some odd-even staggering and fluctuations, which could be caused by the nuclear residual pairing interaction and the nuclear surface diffuseness correction effect. The blue solid curve and the red dashed curve denote the results of two analytic expressions, in which the coefficients are determined by fitting open circles. By performing a two-parameter fit ofasym(A) , attained previously with the relation Eq. (11), we obtained the values ofS0 andκ . Within the 95% confidence interval,S0=26.575±0.271 andκ=0.99±0.035 , and the corresponding rms deviation is 389 keV, if the mass dependence of the symmetry energy coefficient isasym(A)= S0(1−κA−1/3) . By adoptingasym(A)=S0/(1+κA−1/3) , one obtainsS0=30.102±0.741 MeV andκ=2.091±0.166 with an rms deviation of 398 keV. This result is shown in Fig. 5 by the red dashed curve. The obtained values ofS0 andκ are in agreement withS0=(31.1±1.7 ) MeV andκ=(2.31±038 ), given by Min Liu [28].Figure 5. (color online) The symmetry energy coefficient for finite nuclei
asym(A) as a function of the mass number A from Eq. (11) (open circles). Shell corrections S(N,Z) from Ref. [20] are removed from the mass excess of nuclei. Blue solid curve and red dashed curve denote the results of two analytic expressions, in which the coefficients are determined by fitting to the open circles. -
In summary, we proposed a method to determine two important parameters,
ZA andbA , in the well-known Bohr-Wheeler mass parabola. The linear relation between the total doubleβ -decay energy and the charge number Z is deduced based on the isobaric mass parabola, which only includes two parameters,ZA andbA . The values of these two parameters were obtained by fitting 251 A values of the experimental total doubleβ -decay energies of nuclei compiled in the AME2016 database. The advantage of this approach is that we could remove the pairing energy termPA in the process, which causes odd-even variation, as well as the mass excessM(A,ZA) of the most stable nuclide. The Coulomb energy coefficientaC=0.691 MeV was determined from the mass difference relation for mirror nuclei0.5bA(A−2ZA)=Δ(1H−n)+ aC(A2/3−1.0583c) , andc=1.1914 MeV. The symmetry energy coefficients in the mass dependence given in the two forms,asym(A) ,asym(A)=30.102/(1+2.091A−1/3) , andasym(A)=26.575(1− 0.99A−1/3) , were determined by fitting the relationasym(A)=0.25bAZA . The obtained values ofS0 andκ are in agreement with the results in other references. This implies that the proposed method for determining the Coulomb energy coefficient and the symmetry energy coefficient is reliable. Further work is in progress, which includes a comparison of the experimental data and the results of several theoretical mass tables.
Nuclear mass parabola and its applications
- Received Date: 2019-08-18
- Available Online: 2019-12-01
Abstract: We propose a method for extracting the properties of the isobaric mass parabola based on the total double