1.China Nuclear Data Center, China Institute of Atomic Energy, P.O. Box 275(41), Beijing 102413, China 2.Guangxi Normal University, Guilin, Guangxi Province, 541004, China Received Date:2019-07-02 Available Online:2019-11-01 Abstract:The structural effect is believed to have no influence on the decay properties of medium and heavy-mass nuclei at excitation energies above the pairing gap. These properties can be described by statistical properties using so-called photon strength functions for different multipolarities, and directly related to the photoabsorption cross-section ($\sigma_{\rm abs}$). $\sigma_{\rm abs}$ is dominated by the electric giant dipole resonance at $\gamma$ energy $\epsilon_\gamma \leqslant 40$ MeV. In this study, we construct two kinds of systematic giant dipole resonance parameters by fitting the experimental photoabsorption cross-sections. One is based on the microscopic relativistic quasiparticle random phase approximation approach, whereas the other is estimated by the phenomenological models within the Lorentzian representation. Both of them are demonstrated ot efficiently describe the experimental photoabsorption cross-sections available for medium to heavy-mass nuclei, and they can obtain more reliable predictions for the unknown nuclear system.
HTML
--> --> -->
2.Theoretical formalismThe theoretical photoabsorption cross-section $ \sigma_{\rm abs}(\epsilon_\gamma) $ as a function of $ \gamma $-ray energy $ \epsilon_\gamma $ is taken as a sum of the terms
where the component $ \sigma_{\rm GDR}(\epsilon_\gamma) $ corresponds to GDR excitation, and $ \sigma_{\rm QD}(\epsilon_\gamma) $ is a quasideuteron contribution (a photoabsorption by a neutron–proton pair), which is taken in accordance with the model proposed by Chadwick et al. [25]. The GDR component $ \sigma_{\rm GDR}(\epsilon_\gamma) $ of the total photoabsorption cross-section is equal to the photoabsorption cross-section of electric dipole $ \gamma $ rays $ \sigma_{\rm E1}(\epsilon_\gamma) $, which is proportional to the strength function $ S_{\rm E1}(\epsilon_\gamma) $
where the isovector dipole strength function $ S_{\rm E1}(E_\gamma) $ can be calculated by a microscopic model, such as: RQRPA equations [13–15, 26]. In standard notation, these equations can be written in matrix form as
Then, the reduced electric transition probability from the ground to an excited state with angular momentum J and excitation energy $ \omega_\nu $ for a multipole operator $ \hat{Q}^J $ is given by
where $ \Gamma(\epsilon_\gamma) $ is the width of the Lorentzian distribution. Phenomenological models assume that the total photoabsorption cross-section has a Lorentzian shape with various expressions for the "width" $ \Gamma_r(\epsilon_\gamma) $. The standard Lorentzian (SLO) model [19, 20] is based on the Brink hypothesis. In the SLO model, the width $ \Gamma_r(\epsilon_\gamma) $ is taken as the constant $ \Gamma_r $:
where $ E_r $ and $ \sigma_r $ are the GDR energy (in $ {\rm MeV} $) and the peak cross-section (in $ {\rm mb} $), respectively. The SLO approach can be used to describe the photon absorption data for medium and heavy-mass nuclei [27, 28]. However, for $ \gamma $-emission, the SLO model underestimates the $ \gamma $-decay spectra at low energies [29]. In the low energy range, the generalized Fermi liquid (GFL) [23] and the enhanced generalized Lorentzian (EGLO) [21, 22] models are based on Kadmenskij-Markushev-Furman (KMF) model [10, 30] and give a better description of the $ {\rm E1} $ strengths at the energies $ \epsilon_\gamma $ close to zero. The photoaborption cross-section of the EGLO is given by the following expression:
and the parameter k reproduces the experimental $ {\rm E1} $ strength around a reference energy $ \epsilon_0 $; $ \chi(E-\gamma = E_r) = 1 $ and $ C_K(E_r) = C_{\rm KMF} $. The value of the factor k depends on the model adopted to describe the nuclear state density, and was obtained from the average resonance capture data, where $ \epsilon_0 = 4.5 $ MeV. If the Fermi gas model is used, k is given by
where the quantity K depicts the Landau parameters of the quasi-particle interaction in the isovector channel of Fermi system, and the energy-dependent width $ \Gamma_m(\epsilon_\gamma) $ is written as:
with $ C_{\rm coll} $ determined by normalizing the total width $ \Gamma_m(\epsilon_\gamma) $ at $ \epsilon_\gamma = E_r $ and $ T_f = 0 $ to the GDR width of a cold nucleus, i.e. $ \Gamma_m(\epsilon_\gamma = E_r) = \Gamma_r $. The component $ \Gamma_{dq} $ is assumed in the following form:
where $ C_{dq} = \sqrt{5\ln2/\pi} = 1.05 $; $ s_2 = E_2+\bar{\beta}_2^2\approx217.16/A^2 $ with $ E_{2+} $ as the energy of the first vibrational quadrupole state, and $ \bar{\beta}_2 $ being the effective deformation parameter. It can be noted that the SLO, EGLO, and GFL expressions are not consistent with the general relations between the imaginary part of the response function and photoabsorption cross-section. To avoid this shortcoming, the Modified Lorentzian approach (MLO) was proposed [11]. The photoabsorption cross-section within the MLO model has the following form:
where at zero excitation energy $ \Gamma_r = \Gamma_\gamma(\epsilon_\gamma = E_r) $. In the case where the excitation energy is not too high and $ \epsilon_\gamma $ ranges from zero up to the GDR energy, the MLO1 width $ \Gamma(\epsilon_\gamma) $ can be expressed as
where $ a = \Gamma_r/E_r = C_{\rm KMF}\cdot \epsilon_\gamma $ if the normalization condition $ \Gamma_r = \Gamma_\gamma(\epsilon_\gamma = E_r) $ is adopted for cold nuclei. This model is denoted as the simplified modified Lorentzian (SMLO) model [3]. The values of the parameter a are obtained by fitting the SMLO shape to the experimental photoabsorption cross-sections in spherical nuclei.
3.Calculations and discussionIn this section, we construct two kinds of systematic GDR parameters for medium to heavy nuclei. One is based on the microscopic RQRPA calculation. The other is achieved by phenomenological models. The shape of parameters is obtained by fitting the theoretical calculation to the experiment photoabsorption cross-sections from the EXFOR library [31]. In Table 1 and 2 we have listed 47 nuclei, which we use to fit the systematic GDR parameters. Therein are 25 spherical nuclei and 22 axially deformed nuclei.
Table2.Experimental values of axially deformed nuclei used to fit systematic GDR parameters.
23.1.Microscopic GDR parameters
-->
3.1.Microscopic GDR parameters
In this section, we construct the microscopic GDR parameters within the relativistic quasiparticle random phase approximation (RQRPA) calculation, which is obtained with the NL3 interaction [26]. RQRPA has been proved to be able to provide a reliable description of the GDR centroid and fraction of the energy-weighted sum rule exhausted by the E1 mode. However, this can only be used to solve the spherical nuclei, and if we want to reproduce the photoabsorption cross-section of the experimental data, it is necessary to find a good width of the Lorentzian distribution. In previous RQRPA calculations, the width of the Lorentzian distribution is usually treated as a constant [15, 16]. Here, we follow the method in Ref. [51] and define an energy-dependent width parameter as $ \Gamma(\epsilon_\gamma) = \Gamma\sqrt{\epsilon_\gamma/E_{\rm E1}} $ in Eq.(5), where $ E_{\rm E1} $ is the peak energy of the giant dipole resonance, and $ \Gamma $ is one of the parameters we use to fit the experimental data. To reproduce the photoabsorption cross-section of the experimental data, we add another parameter G in Eq.(2) to adjust the height of the curve of the photoaborption cross-section
With the help of the MINUIT package, we adjust G and $ \Gamma $ to improve the RQRPA calculation of the photoabsorption cross-section in Table 1. As shown in Table 3, the total $ \chi^2_{\rm tot} $ of the RQRPA with experimental photoabsorption cross-section has been improved more than 10 times.
$\Gamma$
G
$\chi_{\rm tot}^2$
RQRPA
2.0 (C)
1.0
270.0
RQRPA(fitted)
1.59 (ED)
0.65
21.0
Table3.Parameters and $\chi_{\rm tot}^2$ of microscopic GDR parameters within RQRPA. C depicts constant width of Lorentzian distribution., $ED$ denotes energy-dependent width.
In Fig. 1, we present the comparison of our microscopic GDR parameters (RQRPA) calculations with the experimental photoabsorption cross-sections and other GDR calculations such as: QRPA, MLO1, and SMLO for 52Cr, 90Zr, 120Sn, and 208Pb. The QRPA calculations are based on the Sly4 Skyrme force [8], and the results of MLO1* are based on the parameters from RIPL-3 [2]. The curves of the SMLO are calculated by phenomenological GDR parameters, which we discuss in the next section. Fig. 1shows that the centroid energy of RQRPA is closer to the experimental GDR energy. By adjusting G and $ \Gamma $ parameters, we can reproduce the experimental photoabsorption cross-section as well as the MLO1*. Our microscopic GDR parameter calculations can be used to estimate the spherical nuclei, and the cross-section for heavy-mass nuclei are lower than the ones obtained by experimental data. Figure1. (color online) Comparison of SMLO (black solid line) and RQRPA (red dashed line) calculations with experimental photoabsorption cross-sections and other systematic GDR parameters calculations QRPA [8] (blue dash-dot line) and MLO1* from RIPL-3 [2] (magenta dash-dot-dot line) for 52Cr, 90Zr, 120Sn and 208Pb.
23.2.Phenomenological GDR parameters -->
3.2.Phenomenological GDR parameters
In this section, systematic GDR parameters, based on phenomenological models, are systematically constructed. As introduced in Section 2, the phenomenological models assume that the total photoabsorption cross-section has a Lorentzian shape. In spherical nuclei, this depends on GDR parameters: the energy $ E_r $, the value $ \sigma_r $ of the photoabsorption cross-section, and width $ \Gamma_r $ at GDR energy. In axially deformed nuclei, this is taken as a sum of two components with two sets of the GDR parameters. There are several different models of the photoabsorption cross-section description. The major difference of various models is the expressions for the width of the curve. In the standard Lorentzian model (SLO) [19, 20], the width is adopted as a constant that is equal to the GDR width $ \Gamma_r $. In other models, such as the enhanced generalized Lorentzian model (EGLO) [21, 22], generalized Fermi liquid model (GFL) [23], hybrid model (GH) [24], and modified Lorentzian model (MLO) [11], the width is taken as dependent on the $ \gamma $-ray energy $ \Gamma(\epsilon_\gamma) $. Although in the different models, $ \Gamma(\epsilon_\gamma) $ has various formulas, $ \Gamma_r $ is the only parameter of the $ \gamma $-ray energy dependent width. Hence, in our phenomenological GDR parameters, first we have to define three basic values for the photoabsorption cross-section. They are assumed in the form of the mass number A and deformation $ \beta_2 $:
where $ a_i(i = 1,\cdots,8) $ are adjustable parameters, $ I = (N-Z)/A $, and $ \sigma_{\rm{{TRK}}} = 60NZ/A $ is the classical dipole Thomas-Reiche-Kuhn sum rule [52]. Comparing with the Eq.(137) in RIPL-3 Ref. [2], we added three extra parts with three more parameters, which correspond to the neutron-proton asymmetry I and deformation $ \beta_2 $ of the nuclei. Our phenomenological GDR parameters are not only for the spherical nuclei, but also can be used to calculate axially deformed nuclei. We adopt the following formula:
where indexes $ 1 $ and $ 2 $ correspond to the collective motion along and perpendicular to the axis of symmetry, respectively, with the relative semi-axes of a spheroid. We added eight parameters $ b_i (i = 1,\cdots,8) $ to change the value of width $ \Gamma_{r} $ and photoabsorption cross-section $ \sigma_r $ with deformation $ \beta_2 $.
Effective quadrupole deformation parameters $ \alpha_2 $ or $ \beta_2 $ were determined from ground state deformation parameters $ \beta_n\equiv\alpha_n/\sqrt{(2n+1)/4\pi} $, with the nuclear radius expansion expressed in spherical harmonics. $ \beta_n $ parameters were provided by Ref. [53]. In Table 4, we display all the phenomenological GDR parameters and total $ \chi^2 $-values of MLO1, MLO2, MLO3, EGLO, GFL, SLO, GH, and SMLO models. Symbols SLO* and MLO1* denote the parameters in RIPL-3 [2]. In Fig. 2, we display the $ \chi^2 $ of each nucleus with different models. The red and blue points are $ \chi^2 $ values for the spherical and deformed nuclei, respectively. Comparing with the global GDR parameters of RIPL-3, our systematic approach, Eq. (17) and (18), represents a significantly better fit of the experimental photoabsorption cross-section, especially for the deformed nuclei, since the total $ \chi^2 $ value of the SMLO model is the smallest, and the formula of the SMLO is very simple. Therefore, the SMLO model performs better than others.
MLO1
MLO2
MLO3
EGLO
GFL
SLO
GH
SMLO
SLO*
MLO1*
$ a_1 $
42.722
41.504
42.563
43.1769
39.2025
39.7783
38.4423
42.0152
27.47
28.69
$ a_2 $
16.0865
16.7204
16.1593
16.0834
18.8849
17.6336
18.6947
16.4032
22.063
21.731
$ a_3 $
0.0111
0.0124
0.0122
0.0048
0.0104
0.012
0.01
0.0062
0.0277
0.0285
$ a_4 $
1.0375
1.0157
1.0071
1.1963
1.0186
1.0095
1.0114
1.1682
1.9
1.9
$ a_5 $
0.0452
0.0477
0.0457
0.0309
0.0395
0.0459
0.0384
0.0364
1.222
1.267
$ a_6 $
-9.0294
-10.4079
-10.3748
-5.00
-6.6123
-11.7035
-10.4802
-10.5713
0
0
$ a_7 $
8.8811
9.6039
10.3457
4.1809
4.4196
9.7693
7.6594
10.3444
0
0
$ a_8 $
0.0039
0.0049
0.0043
0.0017
0.0023
0.0051
0.0039
0.0014
0
0
$ b_1 $
?0.8837
?0.903
?0.8877
?0.7608
?0.8142
?0.8966
?0.8669
?0.7962
0
0
$ b_2 $
0.2753
0.5083
0.395
0.0086
0.2809
0.5838
0.5052
0.0629
0.026
0.026
$ b_3 $
1.0789
0.8666
0.9495
2.2425
1.1868
0.7884
0.9247
1.5836
1.91
1.91
$ b_4 $
0.0855
0.0503
0.0758
0.1945
0.4745
0.0757
0.1551
0.1821
0.026
0.026
$ b_5 $
1.4629
1.5903
1.4809
1.4222
0.9907
1.4819
1.3278
1.2703
1.91
1.91
$ b_6 $
?0.184
?0.1789
?0.1604
?0.1539
?0.3373
?0.1666
?0.2366
?0.1882
0
0
$ b_7 $
0.0815
?0.0221
0.036
1.0276
0.3179
0.1194
0.2426
0.2623
0
0
$ b_8 $
0.2448
0.2307
0.2474
0.2914
0.2475
0.229
0.2303
0.304
0.333
0.333
$ \chi^2 $
9.8769
12.8526
10.7246
9.742
10.1275
15.9214
15.6991
9.5161
25.8048
19.0095
Table4.Parameters and $ \chi^2 $-values of systematics within MLO1, MLO2, MLO3, EGLO, GFL, SLO, GH, and SMLO models.
Figure2. (color online) The $\chi^2$ of each nuclei with different models, red and blue points depict $\chi^2$ values for spherical and deformed nuclei, respectively.
In Figs. 1 and 3, we displayed the comparison of the results of SMLO with the experimental photoabsorption cross-sections and other systematic calculations such as: QRPA and MLO1*. In Fig. 1, we compare cross-sections of the spherical nuclei, and in Fig. 3 the deformed nuclei including: 28Si, 80Se, 127I and 186W are taken into account. The QRPA calculations are based on the Sly4 Skyrme force [8] and MLO1* is the result of the GDR parameters in RIPL-3 [2]. As observed in Figs. 1 and 3, the SMLO model can provide a good description of the GDR peak energy for both spherical and deformed nuclei. Around the GDR region, the agreement of the phtotoabsorption cross-sections between SMLO and experimental data are rather satisfactory. Figure3. (color online) Comparison of SMLO (solid line) calculations with experimental photoabsorption cross-sections for 28Si, 80Se, 127I, and 186W. The QRPA (dash dot line) calculations are based on the Sly4 Skyrme force [8], and MLO1* (dash dot dot line) are the results of GDR parameters recommended by RIPL-3 [2].