Isoscalar giant monopole resonance for drip-line and super heavy nuclei in the framework of a relativistic mean field formalism with scaling calculation

Isoscalar giant monopole resonance for drip-line and super heavy nuclei in the framework of a relativistic mean field formalism with scaling calculation

S. K. Biswal S. K. Patra
September 12, 2019

We study the isoscalar giant monopole resonance for drip-lines and super heavy nuclei in the frame work of a relativistic mean field theory with scaling approach. The well known extended Thomas-Fermi approximation in the non-linear - model is used to estimate the giant monopole excitation energy for some selected light spherical nuclei starting from the region of proton to neutron drip-lines. The application is extended to super heavy region for Z=114 and 120, which are predicted by several models as the next proton magic number beyond Z=82. We compared the excitation energy obtained by four successful force parameters NL1, NL3, NL3 and FSUGold. The monopole energy decreases toward the proton and neutron drip-lines in an isotopic chain for lighter mass nuclei contrary to a monotonous decrease for super heavy isotopes. The maximum and minimum monopole excitation energies are obtained for nuclei with minimum and maximum isospin, respectively in an isotopic chain.

Relativistic Thomas-Fermi approximation, Monopole excitation energy, Compressibility modulus
21.10.Dr, 21.65.Cd, 21.60.-n, 21.10.-k

Editorial 11institutetext: Institute Of Physics,
Sachivalya Marg, 751 005, India.

I Introduction

The study of nuclei far away from the drip-lines has a current research interest due to their very different properties than the nuclei at the stability valley. New properties of these nuclei like the soft giant resonance, the change of magic number, the halo and skin structures and the new decay modes stimulate strongly the research using radioactive ion beam (RIB) tan85 ; han89 . On the other hand, the super heavy nuclei which are on the stability line, but extremely unstable due to the excessive Coulomb repulsion attract much theoretical attention for its resemblance to the highly asymmetry nuclear matter limit og ; kumar89 . These nuclei possess a large amount of collective excitation and their study along an isotopic chain is more informative for the structural evaluation of astrophysical objects like neutron star aru04 . Also, the nuclear symmetry energy, and consequently the proton to neutron ratio, are crucial factors in constructing the equation of state (EOS) for asymmetry nuclear matter.

The compressibility of a nuclear system depends on its neutron-proton asymmetry. Also it is well known that the EOS of an asymmetry accerating object like neutron star substantially influence by it‘s compressibility. Although, the compressibility at various asymmetry is an important quantities, it is not a direct experimental observable. Thus, one has to determine the from the linked experimental quantity (which is directly or indirectly related to ) like isoscalar giant monopole resonance (ISGMR) young99 ; bohigas79 . The ISGMR is a well defined experimental observable, which can be measured precisely through various experimental techniques. The drip-lines and super heavy nuclei are vulnerable and unstable in nature because of the presence of excess neutron and large number of protons, respectively. Thus it is instructive to know the giant monopole resonance, compressibility modulus and other related quantities for both drip-lines and super heavy nuclei. In this context, our aim is to study the giant monopole excitation energy and the compressional modulus of finite nuclei near the drip-line han89 as well as for recently discussed super heavy nuclei with proton numbers Z=114 and 120, which are predicted to be next magic numbers beyond Z=82 with various models sil04 ; rutz95 . In addition, the calculations of Refs. sobi ; patra suggest that these nuclei possess spherical ground state or low-lying spherical excited solutions. More specifically, we aimed to study the following within the frame-work of an extended relativistic Thomas-Fermi approximation:

  • How the isoscalar excitation energy and the finite nuclear compressibility varies in an isotopic chain in drip-lines and super heavy nuclei within a well tested model like relativistic extended Thomas-Fermi frame work using scaling and constrained approaches which is developed by some of us recently patra01 ; patra02 .

  • A comparative study of ISGMR obtained with various parameter sets such as NL1, NL3, NL3 and FSUGold for the same drip-lines and super heavy nuclei.

  • The resonance width , which is mostly the difference between the scaling and constraint excitation energies are analyzed in the isotopic chains of light and super heavy nuclei.

  • Finally, the relation between the finite nuclear compressibility with the infinite nuclear matter values in various force parameter sets are looked for.

In relativistic mean field (RMF) formalism the NL1 parameter set rufa86 is considered to be one of the best interaction for a long time to predict the experimental observables. The excessive large value of asymmetry coefficient MeV questions about the accuracy for the prediction of neutron radius near the drip-line. As a result, the discovery of NL3 parameter set lala97 complement the limitation of NL1 force and evaluates the ground state properties of finite nuclei in an excellent agreement with experiment boguta87 ; lala97 ; patra91 ; gamb90 ; sumi93 ; toki94 . It reproduces the proton or charge radius precisely along with the ground state binding energy. Unfortunately, the experimental data for neutron radius has a large error bar batty , which covers most of the prediction of all relativistic and non-relativistic models brow00 . The FSUGold parameter set todd05 ; fattoyev10 reproduces the ISGMR pretty well with the experimental data for Zr and Pb. There is also a possibility to solve the uncertainty in neutron radius problem roca11 using this interaction. The NL3 force parametrization lala09 claimed to be an improved version of NL3 to reproduce the experimental observables. We used all these forces and made a comparison of their predictive power for various experimental data. Then we selected NL3 as a suitable parameter set for our further investigations for ISGMR and related quantities.

The paper is presented as follows: In section II, we outline in brief the formalism used in the present work. In section III, we discuss our results for the ground state isoscalar giant monopole resonance (ISGMR) for drip-lines and super heavy nuclei. The isoscalar monopole excitation energy and the compressibility modulus of finite nuclei are also analyzed. We give the summary and concluding remarks in section IV.

Ii The Formalism

In this paper we shall make use of the principle of scale invariance to obtain the virial theorem for the relativistic mean field serot86 theory by working in the relativistic Thomas–Fermi (RTF) and relativistic extended Thomas-Fermi (RETF) approximations patra01 ; mario93 ; mario92 ; spei98 ; mario98 ; mario93a . Although, the scaling and constrained calculations are not new, the present technique is developed first time by Patra et al patra01 ; mario10 and not much have been explored for various regions of the periodic chart. Thus, it is interesting to apply the model specially for drip-lines and super heavy nuclei. The calculations will be explored to the region ranging from Z=8 to Z=114, 120, where we can simulate the properties of neutron matter from the neutron-rich finite nuclei. For this purpose we compute moments and average centroid energies of the isoscalar giant monopole resonance (ISGMR) through scaling and constrained self-consistent calculations for ground state.

The detail formalisms of the scaling method are given in Refs. patra01 ; patra02 . For completeness, we have outlined briefly some of the essential expressions which are needed for the present purpose. We have worked with the non-linear Lagrangian of Boguta and Bodmer boguta77 to include the many-body correlation arises from the non-linear terms of the meson self-interaction schiff51 for nuclear many-body system. The nuclear matter compressibility modulus also reduces dramatically by the introduction of these terms, which motivates to work with this non-liner Lagrangian. The relativistic mean field Hamiltonian for a nucleon-meson interacting system is written by serot86 ; patra01 :


Here , , and are the masses for the nucleon (with being the effective mass of the nucleon), , and mesons, respectively and is the Dirac spinor. The field for the -meson is denoted by , for -meson by , for -meson by ( as the component of the isospin) and for photon by . , , and =1/137 are the coupling constants for the , , -mesons and photon respectively. and are the non-linear coupling constants for mesons. By using the classical variational principle we obtain the field equations for the nucleon and mesons. In semi-classical approximation we can write the above Hamiltonian in term of density as:




and is the free part of the Hamiltonian. The total density is the sum of proton and neutron densities. The semi-classical ground-state meson fields are obtained by solving the Euler–Lagrange equations ().


The above field equations are solved self-consistently in an iterative method. The pairing correlation is not included in the evaluation of the equilibrium property as well as monopole excitation energy. The Thomas-Fermi approach is a semi-classical approximation and pairing correlation has a minor role in giant resonance. It is shown in xavier11 ; baldo13 that only for open-shell nuclei, it has a marginal effect on ISGMR energy. As far as pairing correlation is concerned, it is a quantal effect and can be included in a semi classical calculation as an average as it is adopted in semi-empirical mass formula. In Ref. baldo13 , perturbative calculation on top of a semi classical approach is done, and it suggests that pairing correlation is unimportant in such approach like relativistic Thomas-Fermi (RTF) or relativistic extended Thomas-Fermi (REFT) approximations. In our present calculations, the scalar density () and energy density () are calculated using RTF and RETF formalisms. The RETF is the correction to the RTF, where the gradient of density is taken care. This term of the density takes care of the variation of the density and involves more in the surface properties.




In order to study the monopole vibration of the nucleus we have scaled the baryon density patra01 . The normalized form of the baryon density is given by


is the collective co-ordinate associated with the monopole vibration. As Fermi momentum and density are inter-related, the scaled Fermi momentum is given by


Similarly , , and Coulomb fields are scaled due to self-consistence eqs. (7-10). But the field can not be scaled simply like the density and momentum, because the source term of field contains the field itself. In semi-classical formalism, the energy and density are scaled like


The symbol shows an implicit dependence of . With all these scaled variables, we can write the Hamiltonian as:


Here we are interested to calculate the monopole excitation energy which is defined as with is the restoring force and is the mass parameter. In our calculations, is obtained from the double derivative of the scaled energy with respect to the scaled co-ordinate at and is defined as patra01 :


and the mass parameter of the monopole vibration can be expressed as the double derivative of the scaled energy with the collective velocity as


where is the displacement field, which can be determined from the relation between collective velocity and velocity of the moving frame,


with is the transition density defined as


taking . Then the mass parameter can be written as . In non-relativistic limit, and the scaled energy is . The expressions for and can be found in bohigas79 . Along with the scaling calculation, the monopole vibration can also be studied with constrained approach bohigas79 ; maru89 ; boer91 ; stoi94 ; stoi94a . In this method, one has to solve the constrained functional equation:


Here the constrained is . The constrained energy can be expanded in a harmonic approximation as


The second order derivative in the expansion is related with the constrained compressibility modulus for finite nucleus as


and the constrained energy as


In the non-relativistic approach, the constrained energy is related by the sum rule weighted . Now the scaling and constrained excitation energies of the monopole vibration in terms of the non-relativistic sum rules will help us to calculate , i.e., the resonance width bohigas79 ; mari05 ,


Iii Results and Discussions

Nucleus Set BE (calc.) BE (Expt.) (calc.) (Expt.)
O NL1 127.2(118.7) 127.6 2.772(2.636) 2.699
NL3 128.7(120.8) 2.718(2.591)
NL3 128.1(119.5) 2.724(2.603)
FSUG 127.4(117.8) 2.674(2.572)
Ca NL1 342.3(344.783) 342.0 3.501(3.371) 3.478
NL3 341.6(346.2) 3.470(3.343)
NL3 341.5(344.2) 3.470(3.349)
FSUG 340.8(342.2) 3.429(3.327)
Ca NL1 412.7(419.5) 416.0 3.501(3.445) 3.477
NL3 414.6(422.6) 3.472(3.426)
NL3 413.5(420.3) 3.469(3.429)
FSUG 411.2(418.0) 3.456(3.418)
Zr NL1 784.3(801.1) 783.9 4.284(4.232) 4.269
NL3 781.4(801.7) 4.273(4.219)
NL3 781.6(798.7) 4.267(4.219)
FSUG 778.8(797.3) 4.257(4.214)
Sn NL1 989.5(1013.7) 988.7 4.625(4.583) 4.625
NL3 985.4(1014.6) 4.617(4.571)
NL3 986.4(1011.0) 4.609(4.569)
FSUG 984.4(1010.7) 4.611(4.569)
Pb NL1 1638.1(1653.7) 1636.4 5.536(5.564) 5.501
NL3 1636.9(1661.2) 5.522(5.541)
NL3 1636.5(1655.2) 5.512(5.538)
FSUG 1636.2(1661.4) 5.532(5.541)
Set NL1 NL3 NL3 FSUG emperical
0.154 0.150 0.148 0.148 0.17
16.43 16.31 16.30 16.30 15.68
211.7 271.76 258.27 230.0
43.6 38.68 37.4 32.597
0.57 0.594 0.60 0.61 0.6
Table 1: The calculated binding BE and charge radius obtained from relativistic extended Thomas-Fermi (RETF) approximation is compared with relativistic Hartree (with various parameter sets) and experimental results wang12 ; angeli13 . The RETF results are given in the parenthesis. The empirical values bethe71 ; blaizot of nuclear matter saturation density , binding energy per nucleon BE/A, compressibility modulus , asymmetry parameter and ratio of effective mass to the nucleon mass are given in the lower part of the table. The energy is in MeV and radius is in fm.

iii.1 Force parameter of relativistic mean field formalism

First of all, we examined the predictive power of various parameter sets. In this context we selected NL1 as a successful set of past and few recently used forces like NL3, NL3 and FSUGold. The ground state observables obtained by these forces are depicted in Table 1. Along with the relativistic extended Thomas-Fermi (RETF) results, the values with relativistic Hartree are also compared with the experimental data wang12 ; angeli13 . The calculated RMF results obtained by all the force parameters considered in the present paper are very close to the experimental data wang12 ; angeli13 . A detail analysis of the binding energy and charge radius clear that NL1 and FSUGold have a superior predictive power for O in RMF level. The advantage of FSUGold decreases with increase of mass number of the nucleus. Although, the predictive power of the pretty old NL1 set is very good for binding energy, it has a large asymmetry coefficient , which may mislead the prediction in unknown territory, like neutron drip-line or super heavy regions. The RETF prediction of binding energy and charge radius (numbers in the parenthesis) is very poor with the experimental data as compared to the RMF calculations. However, for relatively heavier masses, the ERTF results can be used within acceptable error. In general, taking into account the binding energy BE and root mean square charge radius , one may prefer to use either NL3 or NL3 parametrization.

Nucleus NL1 NL3 NL3 FSUG Expt.
23.31 21.75 27.83 25.97 26.86 25.20 26.97 25.17 21.130.49 0.97 0.11 0.86
20.61 19.77 24.01 23.16 23.32 22.48 22.98 22.30 19.200.40 0.69 0.43 1.03
19.51 18.67 22.69 21.73 22.01 21.11 21.72 20.88 19.900.20 0.68 0.29 0.97
16.91 16.41 19.53 19.03 18.97 18.50 18.60 18.21 0.56 0.37 0.97
15.97 15.50 18.42 17.94 17.90 17.44 17.52 17.13 0.52 0.38 0.90
15.87 15.39 18.29 17.81 17.78 17.32 17.42 17.02 16.10.10 0.51 0.36 0.86
15.76 15.28 18.16 17.67 17.65 17.18 17.31 16.90 15.90.10 0.51 0.34 0.85
15.63 15.19 18.02 17.52 17.51 17.04 17.19 16.77 15.800.10 0.51 0.32 0.83
15.51 15.03 17.87 17.36 17.37 16.89 17.07 16.63 15.60.10 0.50 0.30 0.80
15.38 14.90 17.72 17.20 17.22 16.73 16.94 16.49 15.40.20 0.50 0.28 0.78
15.24 14.76 17.56 17.03 17.07 16.57 16.81 16.34 15.00.20 0.49 0.24 0.77
15.11 14.61 17.40 16.85 16.91 16.40 16.67 16.19 14.800.20 0.48 0.24 0.72
12.69 12.11 14.58 13.91 14.18 13.55 14.04 13.44 14.170.28 0.40 0.14 0.54
11.32 10.60 13.00 12.14 12.64 11.83 12.55 11.79 0.36 0.09 0.45
11.05 10.31 12.68 11.80 12.33 11.50 12.29 11.53 0.35 0.04 0.37
11.28 10.53 12.96 12.07 12.60 11.76 12.48 11.69 0.36 0.12 0.27
11.04 10.28 12.67 11.77 12.33 11.47 12.25 11.47 0.34 0.08 0.42
Table 2: The results of isoscalar giant monopole resonance with various parameter sets for some known nuclei are compared with recent experimental data young13 . The calculations are done with relativistic extended Thomas-Fermi (RETF) approximation using both scaling and constrained schemes. The values of , and are obtained by subtracting the results of (NL3, NL3), (NL3, FSUG) and (NL3, FSUG), respectively. The monopole excitation energies with scaling and constrained are in MeV.

Before accepting NL3 or NL3 as the usuable parameter set for our further calculations, in Table 2, we have given the excitation energy of some selective nuclei both in light and super heavy regions with various parameter sets for some further verification. The isoscalar giant monopole energies and are evaluated using both scaling and constraint calculations, respectively. The forces like NL1, NL3, Nl3 and FSUG have a wide range of compressibility starting from 211.7 to 271.7 MeV (see Table 1). Because of the large variation in of these sets, we expect various values of and with different parametrization. From Table 2, it is noticed that the calculated results for O and Ca differ substantially from the data. Again this deviation of calculated result goes on decreasing with increase of mass number, irrespective of the parameter set. This may be due to the use of semi-classical approximation like Thomas-Fermi and extended Thomas-Fermi. In these approaches, quantal corrections are averaged out. When we are going from light to the heavy and then super heavy nuclei, the surface correction decreases appreciably. Consequently, the contribution to monopole excitation energy decreases with mass number A. In column 11, 12 and 13 of Table 2, the differences in obtained from various parameter sets are given, namely, is the difference in monopole excitation energy obtained by NL3 and NL3. Similarly, and are the ISGMR difference with (NL3, FSUG) and (NL3, FSUG), respectly. The values of , or goes on decreasing with increase of mass number of the nucleus without depending on the parameter used. In other word, we may reach to same conclusion in the super heavy region irrespective of the parameter set. However, it is always better to use a successful parameter set to explore an unknown territory. In this context, it is safer to choose NL3 force for our further exploration. The second observation is also apparent from the Table. It is commonly believe that, mostly the compressibility of the force parameter affect the excitation energy of ISGMR of the nucleus. That means force parameters having different have different excitation energy for the same nucleus. For example, has excitation energy 14.58 and 14.04 MeV with NL3 and FSUGold, respectively. Although, the ground state binding energy of Pb, either with Hartree (RMF) or REFT approximation matches well with NL3 and FSUGold parameter sets (see Table 1), their ISGMR differ by 0.54 MeV, which is quite substantial. This unequal prediction of may be due to the difference in nuclear matter compressibility of the force parametrizations.

Nucleus RMF FRDM Expt.
O 12 29 13 30 12 30 12 27 12 26 12 28#
Ca 34 69 33 71 34 71 34 66 30 73 35# 58 #
Ni 49 94 50 98 50 98 51 94 46 99 48 79
Sn 99 165 100 172 100 172 99 1 64 94 169 99# 138#
Pb 178 275 180 281 180 280 179 269 175 273 178 220#
114 267 375 271 392 274 390 271 376 269 339 285# 289#
120 285 376 288 414 288 410 289 396 287 339 - -
Table 3: The predicted proton and neutron drip-lines PDL and NDL for O, Ca, Ni, Sn, Pb, Z=114 and Z=120 in relativistic mean filed formalism (RMF) with various parameter sets are compared with experimental (where ever available) and Finite Range Droplet Model (FRDM) prediction.

iii.2 Proton and neutron drip-lines

In Table 3 we have shown the proton and neutron drip-lines (PDL and NDL) for various parameter sets. The neutron (or proton) drip-line of an isotope is defined when the neutron (or proton) separation energy (or , where or with BE(N,Z) is the binding energy of a nucleus with N neutron and Z proton. From the table, it is seen that all the interactions predict almost similar proton and neutron drip-lines. If one compares the drip-lines of NL3 and NL3, then their predictions are almost identical, explicitly for lighter mass nuclei. Thus, the location of drip-line with various forces does not depend on its nuclear matter compressibility or asymmetry coefficient. For example, the asymmetry coefficient MeV and MeV for NL1 set and these are 38.68 and 271.76 MeV with NL3 parametrization. The corresponding proton drip-lines for O isotopes are 12 and 13, and the neutron drip-lines are 29 and 30, respectively. The similar effects are noticed for other isotopes of the considered nuclei (see Table 3).

Figure 1: The isoscalar giant monopole resonance (ISGMR) for O, Ca, Ni and Sn isotopes from proton to neutron drip-lines as a function of mass number.
Figure 2: The isoscalar giant monopole resonance (ISGMR) for Pb, Z=114 and Z=120 isotopes starting from proton to neutron drip-lines as a function of mass number.

iii.3 Isoscalar giant monopole resonance

It is well understood that the isoscalar giant monopole resonance has a direct relation with the compressibility of nuclear matter which decides the softness or stiffness of an equation of state blaizot . This EOS also estimates the structure of neutron stars, like mass and radius. Thus, the ISGMR is an intrinsic property of finite nuclei as well as nuclear equation of states and needed to be determined to gain some light for nuclear properties. The ISGMR for O, Ca, Ni, Sn, Pb, Z=114 and Z=120 isotopic series are given in Figs. 1 and 2. The results are calculated by using both constrained and scaling approaches in the isotopic chain, starting from proton to neutron drip-lines. We use the notation with the mass parameter . The figure shows that excitation energy obtained from scaling calculation is always greater than the constrained value. The difference between the monopole excitation of scaling and constrained calculations, generally gives the resonance width , with and in terms of the ratios of the integral moments of the RAP strength function mario92 . It is also equivalent to and from dielectric theorem, we have .

Now consider Fig. 1, where the excitation energy of giant isoscalar monopole resonance for lighter mass nuclei are plotted. For Z=8 the excitation energy decreases towards both proton (A=12, = 22.51 MeV) and neutron drip-lines (A=26, =21.22 MeV). This excitation energy has maximum value near N=Z (here it is a double closed isotope with Z=8, N=8, = 27.83 MeV). Similar trends are followed in isotopic chain of Ca with Z=20. We find maximum excitation energy at Ca ( MeV), whereas is found to be smaller both in proton (A=34, MeV) and neutron drip-lines (A=71, MeV). However, the trends are somewhat different for isotopic chain of higher Z like Z=50, 82, 114 and 120. In these series of nuclei, the excitation energy monotonically decreases starting from proton drip-line to neutron drip-line. For example, Pb and Pb are the proton and neutron drip nuclei having excitation energy and MeV, respectively. Fig. 2 shows clearly the monotonous decrease of excitation energy for super heavy nuclei. This discrepancy between super heavy and light nuclei may be due to Coulomb interaction and large value of isospin difference. For lighter value of Z, the proton drip-line occurs at a combination of proton and neutron where the neutron number is less than or nearer to the proton number. But for higher Z nuclei, the proton drip-line exhibits at a larger isospin. As the excitation energy of a nucleus is a collective property, it varies smoothly with its mass number, which also reflects in the figures. Consider the isotopic chain of Z=50, the drip-line nucleus (A=100) has excitation energy 18.84 MeV and the neutron drip nucleus A=171 has MeV. The difference in excitation energy of these two isotopes is 5.32 MeV. This difference in proton and neutron drip nuclei is 4.31 MeV for Z=82 and this is 2.37 MeV in Z=114. In summary, for higher Z nuclei, the variation of excitation energy in an isotopic chain is less than the lighter Z nucleus. Again, by comparing with the empirical formula of , our predictions show similar variation through out the isotopic chains. Empirically, the value of is found to be 80 bert76 , however if we select for lighter mass isotopes and for super heavy region, then it fits well with our results, which is slightly different than C=80 obtained by fitting the data for stable nuclei bert76 .

Figure 3: The difference between the monopole excitation energies of scaling and constrained calculations as a function of mass number A for O, Ca, Ni and Sn.
Figure 4: Same as Fig. 3, but for Pb, Z=114 and 120.

There is no direct way to calculate in the scaling or constrained method as random phase approximation (RPA). If we compare the excitation energy obtained from scaling calculation with the non-relativistic RPA result, then it is evident that the scaling gives the upper limit of the energy response function. On the other hand, the constrained calculation predicts the lower limit bohigas79 . As a result, the response width is obtained from the root square difference of and . We have plotted the for the light nuclei in Fig. 3 and for super heavy in Fig. 4. For lighter nuclei, is large both in proton and neutron drip-lines. As one proceed from proton to neutron drip-line, the value of decreases up to a zero isospin combination (N=Z or double close) and then increases. For example, = 10.92, 5.0 and 21.62 MeV for , and , respectively. Similar trends are also followed by Z=20, 28 and 50 isotopic chains. This conclusion can be drawn from the results of the excitation energy also (see Figs. 1 and 2), i.e., the difference between the scaling and constrained excitation energies are more in proton and neutron drip-lines as compared to the Z=N region. The value of in an isotopic chain depends very much on the proton number. It is clear from the isotopic chains of for O, Ca, Ni, Sn, Pb and Z=114, 120. All the considered series have their own behavior and show various trends. Generally, for lighter elements, it decreases initially to some extent and again increases monotonously. On the other hand for heavier nuclei like Pb, Z=114 and 120 this character of with mass number is somewhat different and can be seen in Fig. 4.

iii.4 Compressibility modulus for finite nuclei

Figure 5: The compressibility modulus obtained by both scaling and constrained approaches in the isotopic series of O, Ca, Ni and Sn.
Figure 6: Same as Fig. 5, but for Pb, Z=114 and 120.

The nuclear matter compressibility is a key quantity in the study of equation of state. It is the second derivative of the energy functional with respect to density at the saturation and is defined as , which has a fixed value for a particular force parametrization. It is well understood that a larger of a parameter set, gives stiff EOS and produce a massive neutron star aru04 . It has also a direct relation with the asymmetry energy coefficient of the parameter set estal01 . In the limit A approaches to infinitely large, the finite nucleus can be approximated to infinite nuclear matter system (N=Z for symmetry and for asymmetry matter). Thus, it is instructive to study the nature of compressibility of finite nucleus in the isotopic chain of finite nucleus. Here, we calculate the as a function of mass number for the light nuclei considered in the present study (O, Ca, Ni, Sn) and then extend the calculations to Pb, Z=114 and 120 in the super heavy region. Our calculated results are shown in Figs. 5 and 6. The compressibility of finite nuclei follows same trend as the excitation energy. For light nuclei, the compressibility has small value for proton and neutron drip-lines, whereas it is maximum in the neighborhood of double close combination.

It can be easily understood from Fig. 5 that, at a particular proton to neutron combination, the is high, i.e., at this combination of N and Z, the nucleus is more compressible. In other word, larger the compressibility of a nucleus, it will be more compressible. Here, it is worthy to mention that the nuclear system becomes less compressibility near both the neutron and proton drip-lines. This is because of the the instability originating from the repulsive part of the nuclear force, revealing a rich neutron-proton ratio, which progressively increases with the neutron/proton number in the isotope without much affecting to the density satpathy04 . Similar to the excitation energy, it is found that obtained by scaling method is always higher than the constrained calculation. The decrease in compressibility in the drip-line regions are prominent in constrained calculation than the scaling results. From leptodermous expansion blaizot , we can get some basic ideas about this decreases in the vicinity of drip-lines. The expression for finite nucleus compressibility can be written as


where . The coefficient is negative, so compressibility decreases with . For Ca chain, the compressibility obtained by scaling and constrained calculations are compared with the Hartree-Fock plus RPA method blaizot in Fig. 5. From Fig. 5, one can see that evaluated by semi-classical approximation deviates from RPA results for lighter isotopes contrary to the excellent matching with the heavier mass of Ca isotopes. This is because of the exclusion of the quantal correction in the semi-classical formalism. At higher mass nuclei, this correction becomes negligible and compares to the RPA predictions. This result is depicted in Fig. 6 for Pb and super heavy chain of nuclei. Here the results show completely different trends than the lighter series. The compressibility has higher value in the vicinity of proton drip-line and decreases monotonically towards the neutron drip-line. This is because, for high Z-series, the proton drip-line appears at greater value of N in contrast to the lighter mass region. Again, the compressibility decreases with neutron number from proton to neutron drip-lines.

Figure 7: Compressibility for finite nuclei obtained by scaling calculation versus nuclear matter compressibility .

Finally, we would like to see the trend of with nuclear matter compressibility for various force parameters and also with the size of a nucleus which can reach the infinite nuclear matter limit. For this we choose , and Ca as the selected candidates and shown in Fig. 7. Although, the super heavy nuclei approach the nuclear matter limit, we can not reproduce the from . This may be due to the asymmetry needed to form a bound nucleus, which is the reason for the deviation. That means, the asymmetry of and differs significantly (where ), which is the main source of deviation of from . Also, this deviation arises due to the surface contribution of the finite nuclei. For a quantitative estimation, we have calculated the for different force parameters having various at saturation. We find almost a linear variation of with for the considered nuclei as shown in Fig. 7. For Ca isotopes also we find a similar nature, but smaller than the super heavy nuclei.

Iv Summary and Conclusions

In summary, we have calculated the isoscalar giant monopole resonance for O, Ca, Ni, Sn, Pb, Z=114 and Z=120 isotopic series starting from the proton to neutron drip-lines. The recently developed scaling approach in a relativistic mean field theory is used. A simple, but accurate constrained approximation is also performed to evaluate the isoscalar giant monopole excitation energy. From the scaling and constrained ISGMR excitation energies, we have evaluated the resonance width for the whole isotopic series. This is obtained by taking the root square difference of and . The value of is always higher than the constrained result . In sum rule approach, the can be compared with the higher and as the lower limit of the resonance width. In general, we found an increasing trend of for both lighter and super heavy region near the proton and neutron drip-lines. The magnitude of is predicted to be minimum in the vicinity of N=Z or in the neighborhood of double close nucleus and it is maximum for highly asymmetry systems. In the present paper, we have also estimated the compressibility of finite nucleus. For some specific cases, the compressibility modulus is compared with the nuclear matter compressibility and found a linear variation among them. It is also concluded that the nucleus becomes less compressible with the increase of neutron or proton number in an isotopic chain. Thus the neutron-rich matter, like neutron star as well as drip-line nuclei are less compressible than the normal nuclei. In case of finite drip-line nuclei, the nucleus is incompressible, although it possess a normal density.

V Acknowledgment:

We thank Profs. X. Viñas and M. Centelles for a careful reading of the manuscript. We also thank Mr. S. K. Singh and Mr. M. Bhuyan for discussions.


  • (1) I. Tanihata, H. Hamagaki, O. Hashimoto, Y. Shida, N. Yoshikawa, K. Sugimoto, O. Yamakawa, T. Kobayashi, and N. Takahashi, Phys. Rev. Lett. 55, 2676 (1985); A. Ozawa, T. Kobayashi, T. Suzuki, K. Yoshida, and I. Tanihata Phys. Rev. Lett. 84, 5493 (2000); I. Tanihata, J. Phys. G: Nucl. Part. Phys. 22, 157 (1996).
  • (2) P. G. Hansen and B. Jonson, Europhys. Lett. 4, 409 (1989).
  • (3) Yu. Ts. Oganessian et al., Phys. Rev. Lett. 104, 142502 (2010); Yu. Ts. Oganessian et al., Phys. Rev. Lett. 83, 3154 (1999).
  • (4) K. Kumar, Superheavy Elements, (Adam Hilger, Bristol, 1989).
  • (5) P. Arumugam, B.K. Sharma, P.K. Sahu, S.K. Patra, Tapas Sil, M. Centelles and X. Viñas, Phys. Lett. B 601, 51 (2004).
  • (6) D. H. Youngblood, H.L. Clark and Y.-W. Lui, Phys. Rev. Lett. 82, 691 (1999).
  • (7) O. Bohigas, A. Lane and J. Martorell, Phys. Rep. 51, 267 (1979).
  • (8) T. Sil, S. K. Patra, B. K. Sharma, M. Centelles and X. Viñas, Phys. Rev. C 69, 054313 (2004); M. Bhuyan and S. K. Patra, Mod. Phys. Lett. A 27, 1250173 (2012).
  • (9) K. Rutz, J.A. Maruhn, P.-G. Reinhard and W. Greiner, Nucl. Phys. A 590, 680 (1995).
  • (10) P. Jachimowicz, M. Kowal and J. Skalski, Phys. Rev. C 83, 054302 (2011); A. Sobiczewski, F. A. Gareev and B. N. Kalinkin, Phys. Lett. 22, 500 (1966); A. Sobiczewski, Z. Patyk and S. C. Cwiok, Phys. Lett. 224, 1 (1989); A. Sobiczewski, Acta Phy. Pol. B 41, 157 (2010); Z. Ren, Phys. Rev. C 65, 051304(R) (2002).
  • (11) M. Bhuyan, S. K. Patra and Raj K. Gupta, Phys. Rev. C 84, 014317 (2011); S. K. Patra, M. Bhuyan, M. S. Mehta and Raj K. Gupta, Phys. Rev. C 80, 034312 (2009); S.K. Patra, Cheng-Li Wu, C.R. Praharaj, Raj K. Gupta, Nucl. Phys. A 651, 117 (1999).
  • (12) S. K. Patra, X. Viñas, M. Centelles and M. Del Estal, Nucl. Phys. A 703, 240 (2002); ibid Phys. Lett. B 523, 67 (2001); Chaoyuan Zhu and Xi-Jun Qiu, J. Phys. G 17, L11 (1991).
  • (13) S. K. Patra, M. Centelles, X. Viñas, and M. Del Estal, Phys. Rev. C 65, 044304 (2002).
  • (14) P. G. Reinhard et al. Z. Phys. A 323, 13 (1986).
  • (15) W. Pannert, P. Ring and J. Boguta, Phys. Rev. Lett. 59, 2420 (1987).
  • (16) G. A. Lalazissis, J. König and P. Ring, Phys. Rev. C 55, 540 (1997).
  • (17) S. K. Patra and C. R. Praharaj, Phys. Rev. C 44, 2552 (1991);
  • (18) A. R. Bodmer, Nucl. Phys. A 526, 703 (1991); A. R. Bodmer and C. E. Price, Nucl. Phys. A 505, 123 (1989); Y. K. Gambhir, P. Ring, and A. Thimet, Ann. Phys. (N.Y.) 198, 132 (1990).
  • (19) K. Sumiyashi, D. Hirata, H. Toki and H. Sagawa, Nucl. Phys. A 552 437 (1993).
  • (20) Y. Sugahara and H. Toki, Nucl. Phys. A 579, 557 (1994); S. Gmuca, J. Phys. G 17, 1115 (1991); Z. Phys. A 342, 387 (1992); Nucl. Phys. A 547, 447 (1992).
  • (21) C. J. Batty, E. Friedman, H. J. Gils and H. Rebel, Adv. Nucl. Phys. 19, 1 (1989).
  • (22) B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000).
  • (23) B. G. Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett. 95, 122501 (2005).
  • (24) F. J. Fattoyev, C. J. Horowitz, J. Piekarewicz and G. Shen, Phys. Rev. C82, 055803 (2010).
  • (25) X. Roca-Maza, M. Centelles, X. Viñas and M. Warda, Phys. Rev. Lett. 106, 252501 (2011).
  • (26) G. A. Lalazissis, S. Karatzikos, R. Fossion, D. Pena Arteaga, A. V. Afanasjev, P. Ring, Phys. Lett. B 671, 36 (2009).
  • (27) B. D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986).
  • (28) M. Centelles, X. Viñas, M. Barranco and P. Suhuck, Ann. Phys. (NY) 221, 165 (1993).
  • (29) M. Centeless, X. Viñas, M. Barranco, S. Marco and R. J. Lombard, Nucl. Phys. A 537, 486 (1992).
  • (30) C. Speichers, E. Engle and R. M. Dreizler, Nucl. Phys. A 562, 569 (1998).
  • (31) M. Centeless, M. Del Estal and X. Viñas, Nucl. Phys. A 635, 193 (1998).
  • (32) M. Centelles, X. Viñas, M. Barranco, N. Ohtsuka, A. Faessler, Dao T. Khoa and H. Müther, Phys. Rev. C 47, 1091 (1993).
  • (33) M. Centelles, S. K. Patra, X Roca-Maza, B. K. Sharma, P. D. Stevenson and X. Viñas, J. Phys. G: Nucl. Part. Phys. 37, 075107 (2010).
  • (34) J. Boguta, A.R. Bodmer, Nucl. Phys. 292, 413 (1977).
  • (35) L. I. Schiff, Phys. Rev. 80, 137 (1950); 83, 239 (1951); 84, 1 (1950).
  • (36) M. Baldo, L. M. Robledo, P. Schuck and X. Viñas, Phys. Rev. C87, 064305 (2013).
  • (37) X. Viñas, P. Schuck and M. Farine J. Phys. Conf. Ser. 321, 012024 (2011).
  • (38) T. Maruyama and T. Suzuki, Phys. Lett. B 219, 43 (1989).
  • (39) H. F. Boersma, R. Malfliet and O. Scholten, Phys. Lett. B 269, 1 (1991).
  • (40) M. V. Stoitov, P. Ring, and M. M. Sharma, Phys. Rev. C 50, 1445 (1994).
  • (41) M. V. Stoitsov, M. L. Cescato, P. Ring and M. M. Sharma, J. Phys. G: Nucl. Part. Phys. 20, LI 49 (1994).
  • (42) M. Centelles, X. Viñas, S. K. Patra, J. N. De and Tapas Sil, Phys. Rev. C 72, 014304 (2005).
  • (43) M. Wang, G. Audi, A. H. Wapstra, F. G. Kondev, M. MarCormick, X. Xu and B. Pfeiffer, Chinese Physics C 36, 1603 (2012).
  • (44) I. Angeli and K.P. Marinova, At. Data Nucl. Data Tables, 99, 69 (2013).
  • (45) H. A. Bethe, Ann. Rev. Nucl. Sci. 21, 93 (1971).
  • (46) J. P. Blaizot, Phys. Rep. 64, 171 (1980).
  • (47) D. H. Youngblood, Y.-W. Lui, Krishichayan, J. Button, M. R. Anders, M. L. Gorelik, M. H. Urin and S. Shlomo Phys. Rev. C 88, 021301 (2013).
  • (48) F. E. Bertrand, Ann. Rev. Nucl. Part. Sc. 26, 456 (1976); J. Speth, A. V. Woude, Rep. Porg. Phys. 44, 719 (1981); K. Goeke, J. Speth, Ann. Rev. Nucl. Part. Sc. 32, 65 (1982).
  • (49) M. Del Estal, M. Centelles, X. Vin˜as, and S. K. Patra, Phys. Rev. C 63, 024314 (2001).
  • (50) L. Satpathy and S. K. Patra, J. Phys. G: Nucl. Part. Phys. 30, 771 (2004); S K Patra, R K Choudhury and L Satpathy, J. Phys. G: Nucl. Part. Phys. 37, 085103 (2010).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description