Isoscalar giant monopole resonance for dripline and super heavy nuclei in the framework of a relativistic mean field formalism with scaling calculation
Abstract
We study the isoscalar giant monopole resonance for driplines and super heavy nuclei in the frame work of a relativistic mean field theory with scaling approach. The well known extended ThomasFermi approximation in the nonlinear  model is used to estimate the giant monopole excitation energy for some selected light spherical nuclei starting from the region of proton to neutron driplines. 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 driplines 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.
pacs:
21.10.Dr, 21.65.Cd, 21.60.n, 21.10.kEditorial
^{1}^{1}institutetext:
Institute Of Physics,
Sachivalya Marg, 751 005, India.
I Introduction
The study of nuclei far away from the driplines 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 neutronproton 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 driplines 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 driplines 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 dripline 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 lowlying spherical excited solutions. More specifically, we aimed to study the following within the framework of an extended relativistic ThomasFermi approximation:

How the isoscalar excitation energy and the finite nuclear compressibility varies in an isotopic chain in driplines and super heavy nuclei within a well tested model like relativistic extended ThomasFermi 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 driplines 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 dripline. 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 nonrelativistic 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 driplines 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 ThomasFermi (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 driplines 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 neutronrich finite nuclei. For this purpose we compute moments and average centroid energies of the isoscalar giant monopole resonance (ISGMR) through scaling and constrained selfconsistent 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 nonlinear Lagrangian of Boguta and Bodmer boguta77 to include the manybody correlation arises from the nonlinear terms of the meson selfinteraction schiff51 for nuclear manybody system. The nuclear matter compressibility modulus also reduces dramatically by the introduction of these terms, which motivates to work with this nonliner Lagrangian. The relativistic mean field Hamiltonian for a nucleonmeson interacting system is written by serot86 ; patra01 :
(1)  
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 nonlinear coupling constants for mesons. By using the classical variational principle we obtain the field equations for the nucleon and mesons. In semiclassical approximation we can write the above Hamiltonian in term of density as:
(2) 
where
(3) 
(4) 
(5) 
(6) 
and is the free part of the Hamiltonian. The total density is the sum of proton and neutron densities. The semiclassical groundstate meson fields are obtained by solving the Euler–Lagrange equations ().
(7) 
(8) 
(9) 
(10) 
The above field equations are solved selfconsistently in an iterative method. The pairing correlation is not included in the evaluation of the equilibrium property as well as monopole excitation energy. The ThomasFermi approach is a semiclassical approximation and pairing correlation has a minor role in giant resonance. It is shown in xavier11 ; baldo13 that only for openshell 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 semiempirical 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 ThomasFermi (RTF) or relativistic extended ThomasFermi (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.
(11) 
with
(12) 
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
(13) 
is the collective coordinate associated with the monopole vibration. As Fermi momentum and density are interrelated, the scaled Fermi momentum is given by
(14) 
Similarly , , and Coulomb fields are scaled due to selfconsistence eqs. (710). But the field can not be scaled simply like the density and momentum, because the source term of field contains the field itself. In semiclassical formalism, the energy and density are scaled like
(15)  
(16) 
The symbol shows an implicit dependence of . With all these scaled variables, we can write the Hamiltonian as:
(17)  
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 coordinate at and is defined as patra01 :
(18)  
and the mass parameter of the monopole vibration can be expressed as the double derivative of the scaled energy with the collective velocity as
(19) 
where is the displacement field, which can be determined from the relation between collective velocity and velocity of the moving frame,
(20) 
with is the transition density defined as
(21) 
taking . Then the mass parameter can be written as . In nonrelativistic 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:
(22) 
Here the constrained is . The constrained energy can be expanded in a harmonic approximation as
(23) 
The second order derivative in the expansion is related with the constrained compressibility modulus for finite nucleus as
(24) 
and the constrained energy as
(25) 
In the nonrelativistic 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 nonrelativistic sum rules will help us to calculate , i.e., the resonance width bohigas79 ; mari05 ,
(26) 
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 
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 ThomasFermi (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 dripline 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 
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 semiclassical approximation like ThomasFermi and extended ThomasFermi. 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.  

NL1  NL3  NL3  FSUG  
PDL  NDL  PDL  NDL  PDL  NDL  PDL  NDL  PDL  NDL  PDL  NDL  
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     
iii.2 Proton and neutron driplines
In Table 3 we have shown the proton and neutron driplines (PDL and NDL) for various parameter sets. The neutron (or proton) dripline 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 driplines. If one compares the driplines of NL3 and NL3, then their predictions are almost identical, explicitly for lighter mass nuclei. Thus, the location of dripline 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 driplines for O isotopes are 12 and 13, and the neutron driplines are 29 and 30, respectively. The similar effects are noticed for other isotopes of the considered nuclei (see Table 3).
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 driplines. 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 driplines (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 driplines (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 dripline to neutron dripline. 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 dripline 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 dripline 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 dripline 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 .
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 nonrelativistic 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 driplines. As one proceed from proton to neutron dripline, 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 driplines 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
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 driplines, 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 driplines. This is because of the the instability originating from the repulsive part of the nuclear force, revealing a rich neutronproton 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 dripline 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 driplines. The expression for finite nucleus compressibility can be written as
(27) 
where . The coefficient is negative, so compressibility decreases with . For Ca chain, the compressibility obtained by scaling and constrained calculations are compared with the HartreeFock plus RPA method blaizot in Fig. 5. From Fig. 5, one can see that evaluated by semiclassical 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 semiclassical 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 dripline and decreases monotonically towards the neutron dripline. This is because, for high Zseries, the proton dripline appears at greater value of N in contrast to the lighter mass region. Again, the compressibility decreases with neutron number from proton to neutron driplines.
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 driplines. 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 driplines. 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 neutronrich matter, like neutron star as well as dripline nuclei are less compressible than the normal nuclei. In case of finite dripline 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.
References
 (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, ChengLi 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 XiJun 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. ToddRutel 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. RocaMaza, 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 RocaMaza, 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).