Source author record

M. Centelles

M. Centelles appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2022arXiv

The finite range simple effective interaction including tensor terms

The prediction of single particle level crossing phenomenon between $2p_{3/2}$ and $1f_{5/2}$ orbitals in $Ni$- and $Cu$-isotopic chains by the finite range simple effective interaction without requiring the tensor part is discussed. In this case the experimentally observed crossing could be studied as a function of nuclear matter incompressibility, $K(ρ_0)$. The estimated crossing for the neutron number $N$=46 could be reproduced by the equation of state corresponding to $K(ρ_0)$=240 MeV. However, the observed proton gaps between the $1h_{11/2}$ and $1g_{7/2}$ shells in $Sn$ and $Sb$ isotopic chain, and the neutron gaps between the $1i_{13/2}$ and $1h_{9/2}$ shells in $N$=82 isotones, as well as the shell closure properties at $N$=28 require explicit consideration of a tensor part as the central contribution is not enough to initiate the required level splittings.

preprint2020arXiv

Comment on the paper Eur. Phys. J. A (2019) 55:150

The conclusions of the study published as Eur. Phys. J. A (2019) 55:150 questioning the adequacy of the recently proposed Gogny D1M* interaction for finite-nuclei calculations using harmonic oscillator (HO) basis are revised. The existence of an instability in finite nuclei for D1M* when coordinate-space methods are used to solve the HF equations (as shown in Eur. Phys. J. A(2019) 55:150) is independently confirmed using a computer code based on a quasi-local approximation (QLA) to the HF energy density with finite-range forces. We confirm that the most affected quantity in the coordinate-space calculation is the spatial density at the origin. Our study reveals that some findings concerning these instabilities are not easy to reconcile with the arguments used in Eur. Phys. J. A (2019) 55:150. For instance, some nuclei such as $^4$He and $^{40}$Ca, which diverge in HF mesh-point calculations performed with D1M*, become perfectly stable when Coulomb force is switched off. We have also found instabilities in some nuclei when the D1M interaction is used. Finally, a connection between the occupancy of $s$-orbitals near the Fermi level and the appearance of instabilities is observed. Several convergence and stability studies are performed with HO basis of different sizes and oscillator parameters to demonstrate the robustness of the D1M* results for finite nuclei when the HO basis is used.

preprint2020arXiv

Constraints from GW170817 on the bulk viscosity of neutron star matter and the r-mode instability

We perform a systematic study of the dependence of the r-mode phenomenology in normal fluid pulsar neutron stars on the symmetry energy slope parameter $L$. An essential ingredient in this study is the bulk viscosity, which is evaluated consistently for several equations of state corresponding to different values of the slope parameter $L$. Direct Urca processes, which are allowed from a critical $L$-value onwards, enhance the bulk viscosity and have large influence on the $r$-mode features, such as the instability boundary and spin-down properties of newborn neutron stars. The magnitude of the changes in the $r$-mode properties induced by the direct Urca processes are driven by the $L$-value of the equation of state and the mass of the pulsar. The study has been done by using a family of equations of state of $β$-equilibrated neutron star matter obtained with the finite range simple effective interaction, which provides realistic results for nuclear matter and finite nuclei properties. These equations of state predict the same properties in symmetric nuclear matter and have the same value of the symmetry energy parameter, $E_s(ρ_0)$, but differ in the slope parameter $L$. The range chosen for the variation of $L$ is decided from the tidal deformability data extracted from the GW170817 event and the maximum mass constraint.

preprint2020arXiv

GW170817 constraints analyzed with Gogny forces and momentum-dependent interactions

A set of equations of state obtained from finite-range Gogny forces and momentum-dependent interactions is used to investigate the recent observation of gravitational waves from the binary neutron star merger GW170817 event. For this set of interactions, we have calculated the neutron star tidal deformabilities (related to the second Love number), the mass-radius diagram, and the moment of inertia~($I$). The $I$-Love relation has been verified. We also have found strong correlations among the tidal deformability of the canonical neutron star, its radius, and the derivatives of the nuclear symmetry energy at the saturation density. Most of the obtained results are located within the constraints of the tidal deformabilities extracted from the GW170817 detection.

preprint2020arXiv

Microscopic-Macroscopic Approach for Ground-State Energies Based on the Gogny Force with the Wigner-Kirkwood Averaging Scheme

In the previous paper I \cite{bhagwat20} we have shown that self-consistent Extended Thomas-Fermi (ETF) potentials and densities associated with a given finite-range interaction can be parametrized by generalized Fermi distributions. As a next step, a comprehensive calculation of ground-state properties of a large number of spherical and deformed even-even nuclei is carried out in the present work using the Gogny D1S force within the ETF scheme. The parametrized ETF potentials and densities of paper I are used to calculate the smooth part of the energy and the shell corrections within the Wigner-Kirkwood semiclassical averaging scheme. It is shown that the shell corrections thus obtained, along with a simple liquid drop prescription, yield a good description of ground-state masses and potential energy surfaces for nuclei spanning the entire periodic table.

preprint2020arXiv

Structure and composition of inner crust of neutron stars from Gogny interactions

The detailed knowledge of the inner crust properties of neutron stars might be important to explain different phenomena such as pulsar glitches or the possibility of an {\it r-process} site in neutron star mergers. It has been shown in the literature that quantal effects like shell correction or pairing may play a relevant role to determine the composition of the inner crust of the neutron star. In this paper we construct the equation of state of the inner crust using the finite-range Gogny interactions, where the mean field and the pairing field are calculated with same interaction. We have used the semiclassical Variational Wigner-Kirkwood method along with shell and pairing corrections calculated with the Strutinsky integral method and the BCS approximation, respectively. Our results are compared with those of some popular models from the literature. We report a unified equation of state of the inner crust and core computed with the D1M* Gogny force, which was specifically fabricated for astrophysical calculations.

preprint2016arXiv

Exact versus Taylor-expanded energy density in the study of the neutron star crust-core transition

The importance of the fourth and higher order terms in the Taylor series expansion of the energy of the isospin asymmetric nuclear matter in the study of the neutron star crust-core phase transition is investigated using the finite range simple effective interaction. Analytic expressions for the evaluation of the second and fourth order derivative terms in the Taylor series expansion for any general finite range interaction of Yukawa, exponential or Gaussian form have been obtained. The effect of the nuclear matter incompressibility, symmetry energy and slope parameters on the predictions for the crust-core transition density is examined. The crustal moment of inertia is calculated and the prediction for the radius of the Vela pulsar is analyzed using different equations of state.

preprint2016arXiv

Model dependence of the neutron-skin thickness on the symmetry energy

The model dependence in the correlations of the neutron-skin thickness in heavy nuclei with various symmetry energy parameters is analyzed by using several families of systematically varied microscopic mean field models. Such correlations show a varying degree of model dependence once the results for all the different families are combined. Some mean field models associated with similar values of the symmetry energy slope parameter at saturation density $L$, and pertaining to different families, yield a greater-than-expected spread in the neutron-skin thickness of the $^{208}$Pb nucleus. The effective value of the symmetry energy slope parameter $L_{\rm eff}$, determined by using the nucleon density profiles of the finite nucleus and the density derivative $S^\prime(ρ)$ of the symmetry energy starting from about saturation density up to low densities typical of the surface of nuclei, seems to account for the spread in the neutron-skin thickness for the models with similar $L$. The differences in the values of $L_{\rm eff}$ are mainly due to the small differences in the nucleon density distributions of heavy nuclei in the surface region and the behavior of the symmetry energy at subsaturation densities.

preprint2016arXiv

Relativistic mean field interaction with density dependent meson-nucleon vertices based on microscopical calculations

Although ab-initio calculations of relativistic Brueckner theory lead to large scalar isovector fields in nuclear matter, at present, successful versions of covariant density functional theory neglect the interactions in this channel. A new high precision density functional DD-ME$δ$ is presented which includes four mesons $σ$, $ω$, $δ$, and $ρ$ with density dependent meson-nucleon couplings. It is based to a large extent on microscopic ab-initio calculations in nuclear matter. Only four of its parameters are determined by adjusting to binding energies and charge radii of finite nuclei. The other parameters, in particular the density dependence of the meson-nucleon vertices, are adjusted to non-relativistic and relativistic Brueckner calculations of symmetric and asymmetric nuclear matter. The isovector effective mass $m_{p}^{\ast}-m_{n}^{\ast}$ derived from relativistic Brueckner theory is used to determine the coupling strength of the $δ$-meson and its density dependence.

preprint2015arXiv

Study of spin polarized nuclear matter and finite nuclei with finite range simple effective interaction

The properties of spin polarized pure neutron matter and symmetric nuclear matter are studied using the finite range simple effective interaction, upon its parametrization revisited. Out of the total twelve parameters involved, we now determine ten of them from nuclear matter, against the nine parameters in our earlier calculation, as required in order to have predictions in both spin polarized nuclear matter and finite nuclei in unique manner being free from uncertainty found using the earlier parametrization. The information on the effective mass splitting in polarized neutron matter of the microscopic calculations is used to constrain the one more parameter, that was earlier determined from finite nucleus, and in doing so the quality of the description of finite nuclei is not compromised. The interaction with the new set of parameters is used to study the possibilities of ferromagnetic and antiferromagnetic transitions in completely polarized symmetric nuclear matter. Emphasis is given to analyze the results analytically, as far as possible, to elucidate the role of the interaction parameters involved in the predictions.

preprint2015arXiv

The neutron skin thickness from the measured electric dipole polarizability in $^{68}$Ni, $^{120}$Sn, and $^{208}$Pb

The information on the symmetry energy and its density dependence is deduced by comparing the available data on the electric dipole polarizability $α_D$ of ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb with the predictions of the Random Phase Approximation, using a representative set of nuclear energy density functionals. The calculated values of $α_D$ are used to validate different correlations involving $α_D$, the symmetry energy at the saturation density $J$, the corresponding slope parameter $L$ and the neutron skin thickness $Δr_{\!np}$, as suggested by the Droplet Model. A subset of models that reproduce simultaneously the measured polarizabilities in ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb are employed to predict the values of the symmetry energy parameters at saturation density and $Δr_{\!np}$. The resulting intervals are: $J\!=\!30 \text{-}35$ MeV, $L\!=\!20 \text{-} 66$ MeV; and the values for $Δr_{\!np}$ in ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb are in the ranges: 0.15\text{-}0.19 fm, 0.12\text{-}0.16 fm, and 0.13\text{-}0.19 fm, respectively. The strong correlation between the electric dipole polarizabilities of two nuclei is instrumental to predict the values of electric dipole polarizabilities in other nuclei.

preprint2015arXiv

Unified equation of state for neutron stars on a microscopic basis

We derive a new equation of state (EoS) for neutron stars (NS) from the outer crust to the core based on modern microscopic Brueckner-Hartree-Fock (BHF) calculations using the Argonne $v_{18}$ potential plus three-body forces computed with the Urbana model. To deal with the inhomogeneous structures of matter in the NS crust, we use the recent Barcelona-Catania-Paris-Madrid (BCPM) nuclear energy density functional that is directly based on the same microscopic BHF calculations, and which is able to reproduce the ground-state properties of nuclei along the periodic table. The EoS of the outer crust requires the masses of neutron-rich nuclei, which are obtained through Hartree-Fock-Bogoliubov calculations with the BCPM functional when they are unknown experimentally. To compute the inner crust, Thomas-Fermi calculations in Wigner-Seitz cells are performed with the same functional. Existence of nuclear pasta is predicted in a range of average baryon densities between $\simeq$0.067 fm$^{-3}$ and $\simeq$0.0825 fm$^{-3}$, where the transition to the core takes place. The NS core is computed from the nuclear EoS of the BHF calculation assuming non-exotic constituents (core of $npeμ$ matter). In each region of the star, we discuss the comparison of the new EoS with previous EoSes for the complete NS structure, in particular, with the Lattimer-Swesty EoS and with the Shen et al. EoS widely used in astrophysical calculations. The new microscopically derived EoS fulfills at the same time a NS maximum mass of 2~$M_\odot$ with a radius of 10 km, and a 1.5~$M_\odot$ NS with a radius of 11.7 km.

preprint2014arXiv

From the crust to the core of Neutron Stars on a microscopic basis

Within a microscopic approach the structure of Neutron Stars is usually studied by modelling the homogeneous nuclear matter of the core by a suitable Equation of State, based on a many-body theory, and the crust by a functional based on a more phenomenological approach. We present the first calculation of Neutron Star overall structure by adopting for the core an Equation of State derived from the Brueckner-Hartree-Fock theory and for the crust, including the pasta phase, an Energy Density Functional based on the same Equation of State, and which is able to describe accurately the binding energy of nuclei throughout the mass table. Comparison with other approaches is discussed. The relevance of the crust Equation of state for the Neutron Star radius is particularly emphasised.

preprint2014arXiv

Influence of the single-particle structure on the nuclear surface and the neutron skin

We analyze the influence of the single-particle structure on the neutron density distribution and the neutron skin in Ca, Ni, Zr, Sn, and Pb isotopes. The nucleon density distributions are calculated in the Hartree-Fock+BCS approach with the SLy4 Skyrme force. A close correlation is found between the quantum numbers of the valence neutrons and the changes in the position and the diffuseness of the nuclear surface, which in turn affect the neutron skin thickness. Neutrons in the valence orbitals with low principal quantum number and high angular momentum mainly displace the position of the neutron surface outwards, while neutrons with high principal quantum number and low angular momentum basically increase the diffuseness of the neutron surface. The impact of the valence shell neutrons on the tail of the neutron density distribution is discussed.

preprint2013arXiv

Density dependence of the symmetry energy from neutron skin thickness in finite nuclei

The density dependence of the symmetry energy around saturation density, characterized by the slope parameter L, is studied using information provided by the neutron skin thickness in finite nuclei. An estimate for L is obtained from experimental data on neutron skins extracted from antiprotonic atoms. We also discuss the ability of parity-violating elastic electron scattering to obtain information on the neutron skin thickness in 208Pb and to constrain the density dependence of the nuclear symmetry energy. The size and shape of the neutron density distribution of 208Pb predicted by mean-field models is briefly addressed. We conclude with a comparative overview of the L values predicted by several existing determinations.

preprint2013arXiv

Electric Dipole Polarizability in ${}^{208}$Pb: insights from the Droplet Model

We study the electric dipole polarizability $α_D$ in ${}^{208}$Pb based on the predictions of a large and representative set of relativistic and non-relativistic nuclear mean field models. We adopt the droplet model as a guide to better understand the correlations between $α_D$ and other isovector observables. Insights from the droplet model suggest that the product of $α_D$ and the nuclear symmetry energy at saturation density $J$ is much better correlated with the neutron skin thickness $Δr_{np}$ of ${}^{208}$Pb than the polarizability alone. Correlations of $α_D J$ with $Δr_{np}$ and with the symmetry energy slope parameter $L$ suggest that $α_D J$ is a strong isovector indicator. Hence, we explore the possibility of constraining the isovector sector of thenuclear energy density functional by comparing our theoretical predictions against measurements of both $α_D$ and the parity-violating asymmetry in ${}^{208}$Pb. We find that the recent experimental determination of $α_D$ in ${}^{208}$Pb in combination with the range for the symmetry energy at saturation density $J=[31\pm (2)_{\rm est.}]$\,MeV suggests $Δr_{np}({}^{208}{\rm Pb}) = 0.165 \pm (0.009)_{\rm exp.} \pm (0.013)_{\rm theo.} \pm (0.021)_{\rm est.} {\rm fm}$ and $L= 43 \pm(6)_{\rm exp.} \pm (8)_{\rm theo.}\pm(12)_{\rm est.}$ MeV.

preprint2013arXiv

Electron scattering in isotonic chains as a probe of the proton shell structure of unstable nuclei

Electron scattering on unstable nuclei is planned in future facilities of the GSI and RIKEN upgrades. Motivated by this fact, we study theoretical predictions for elastic electron scattering in the N=82, N=50, and N=14 isotonic chains from very proton-deficient to very proton-rich isotones. We compute the scattering observables by performing Dirac partial-wave calculations. The charge density of the nucleus is obtained with a covariant nuclear mean-field model that accounts for the low-energy electromagnetic structure of the nucleon. For the discussion of the dependence of scattering observables at low-momentum transfer on the gross properties of the charge density, we fit Helm model distributions to the self-consistent mean-field densities. We find that the changes shown by the electric charge form factor along each isotonic chain are strongly correlated with the underlying proton shell structure of the isotones. We conclude that elastic electron scattering experiments in isotones can provide valuable information about the filling order and occupation of the single-particle levels of protons.

preprint2013arXiv

Giant monopole energies from a constrained relativistic mean-field approach

Background: Average energies of nuclear collective modes may be efficiently and accurately computed using a non-relativistic constrained approach without reliance on a random phase approximation (RPA). Purpose: To extend the constrained approach to the relativistic domain and to establish its impact on the calibration of energy density functionals. Methods: Relativistic RPA calculations are compared against the predictions of the corresponding constrained approach using two accurately calibrated energy density functionals. Results: We find excellent agreement-at the 2% level or better-between the predictions of the relativistic RPA and the corresponding constrained approach for magic (or semi-magic) nuclei ranging from 16O to 208Pb. Conclusions: An efficient and accurate method is proposed for incorporating nuclear collective excitations into the calibration of future energy density functionals.

preprint2013arXiv

Nuclear Symmetry Energy: constraints from Giant Quadrupole Resonances and Parity Violating Electron Scattering

Experimental and theoretical efforts are being devoted to the study of observables that can shed light on the properties of the nuclear symmetry energy. We present our new results on the excitation energy [X. Roca-Maza et al., Phys. Rev. C 87, 034301 (2013)] and polarizability of the Isovector Giant Quadrupole Resonance (IVGQR), which has been the object of new experimental investigation [S. S. Henshaw et al., Phys. Rev. Lett. 107, 222501 (2011)]. We also present our theoretical analysis on the parity violating asymmetry at the kinematics of the Lead Radius Experiment [S. Abrahamyan et al. (PREx Collaboration), Phys. Rev. Lett. 108, 112502 (2012)] and highlight its relation with the density dependence of the symmetry energy [X. Roca-Maza et al., Phys. Rev. Lett. 106, 252501 (2011)].

preprint2013arXiv

Symmetry energy of warm nuclear systems

The temperature dependence of the symmetry energy and symmetry free energy coefficients of infinite nuclear matter and of finite nuclei is investigated. For infinite matter, both these coefficients are found to have a weaker dependence on temperature at densities close to saturation; at low but homogeneous densities, the temperature dependence becomes stronger. For finite systems, different definitions of symmetry energy coefficients are encountered in the literature yielding different values. A resolution to this problem is suggested from a global liquid-drop-inspired fit of the energies and free energies of a host of nuclei covering the entire periodic table. The hot nucleus is modeled in a subtracted finite-temperature-Thomas-Fermi framework, with dynamical surface phonon coupling to nucleonic motion plugged in. Contrary to infinite nuclear matter, a substantial change in the symmetry energy coefficients is observed for finite nuclei with temperature.

preprint2012arXiv

Microscopic-Macroscopic Approach for Binding Energies with the Wigner-Kirkwood Method - II

The binding energies of deformed even-even nuclei have been analysed within the framework of a recently proposed microscopic-macroscopic model. We have used the semiclassical Wigner - Kirkwood $\hbar$ expansion up to fourth - order, instead of the usual Strutinsky averaging scheme, to compute the shells corrections in a deformed Woods - Saxon potential including the spin-orbit contribution. For a large set of 561 even-even nuclei with $Z\ge 8$ and $N\ge 8$, we find an {\it rms} deviation from the experiment of 610 keV in binding energies, comparable to the one found for the same set of nuclei using the FRDM of Möller and Nix (656 keV). As applications of our model, we explore its predictive power near the proton and neutron drip lines as well as in the superheavy mass region. Next, we systematically explore the fourth - order Wigner - Kirkwood corrections to the smooth part of the energy. It is found that the ratio of the fourth - order to the second - order corrections behaves in a very regular manner as a function of the asymmetry parameter $I=(N-Z)/A$. This allows to absorb the fourth - order corrections into the second - order contributions to the binding energy, which enables to simplify and speed up the calculation of deformed nuclei.

preprint2012arXiv

Nuclear symmetry energy and neutron skin thickness

The relation between the slope of the nuclear symmetry energy at saturation density and the neutron skin thickness is investigated. Constraints on the slope of the symmetry energy are deduced from the neutron skin data obtained in experiments with antiprotonic atoms. Two types of neutron skin are distinguished: the "surface" and the "bulk". A combination of both types forms neutron skin in most of nuclei. A prescription to calculate neutron skin thickness and the slope of symmetry energy parameter $L$ from the parity violating asymmetry measured in the PREX experiment is proposed.

preprint2012arXiv

The effects of medium on nuclear properties in multifragmentation

In multifragmentation of hot nuclear matter, properties of fragments embedded in a soup of nucleonic gas and other fragments should be modified as compared with isolated nuclei. Such modifications are studied within a simple model where only nucleons and one kind of heavy nuclei are considered. The interaction between different species is described with a momentum-dependent two-body potential whose parameters are fitted to reproduce properties of cold isolated nuclei. The internal energy of heavy fragments is parametrized according to a liquid-drop model with density and temperature dependent parameters. Calculations are carried out for several subnuclear densities and moderate temperatures, for isospin-symmetric and asymmetric systems. We find that the fragments get stretched due to interactions with the medium and their binding energies decrease with increasing temperature and density of nuclear matter.

preprint2012arXiv

The pygmy dipole strength, the neutron radius of ${}^{208}$Pb and the symmetry energy

The accurate characterization of the nuclear symmetry energy and its density dependence is one of the outstanding open problems in nuclear physics. A promising nuclear observable in order to constrain the density dependence of the symmetry energy at saturation is the neutron skin thickness of medium and heavy nuclei. Recently, a low-energy peak in the isovector dipole response of neutron-rich nuclei has been discovered that may be correlated with the neutron skin thickness. The existence of this correlation is currently under debate due to our limited experimental knowledge on the microscopic structure of such a peak. We present a detailed analysis of Skyrme Hartree-Fock (HF) plus random phase approximation (RPA) predictions for the dipole response in several neutron-rich nuclei and try to elucidate whether models of common use in nuclear physics confirm or dismiss its possible connection with the neutron skin thickness. Finally, we briefly present theoretical results for parity violating electron scattering on ${}^{208}$Pb at the conditions of the PREx experiment and discuss the implications for the neutron skin thickness of ${}^{208}$Pb and the slope of the symmetry energy.

preprint2011arXiv

Influence of the nuclear symmetry energy on the structure and composition of the outer crust

We review and extend with nonrelativistic nuclear mean field calculations a previous study of the impact of the nuclear symmetry energy on the structure and composition of the outer crust of nonaccreting neutron stars. First, we develop a simple "toy model" to understand the most relevant quantities determining the structure and composition of the outer crust: the nuclear symmetry energy and the pressure of the electron gas. While the latter is a well determined quantity, the former ---specially its density dependence--- still lacks an accurate characterization. We thus focus on the influence of the nuclear symmetry energy on the crustal composition. For that, we employ different nuclear models that are accurate in the description of terrestrial nuclei. We show that those models with stiffer symmetry energies ---namely, those that generate thicker neutron skins in heavy nuclei and have smaller symmetry energies at subnormal nuclear densities--- generate more exotic isotopes in the stellar crust than their softer counterparts.

preprint2011arXiv

Neutron skin of 208Pb, nuclear symmetry energy, and the parity radius experiment

A precise determination of the neutron skin thickness of a heavy nucleus sets a basic constraint on the nuclear symmetry energy (the neutron skin thickness is the difference of the neutron and proton rms radii of the nucleus). The parity radius experiment (PREX) may achieve it by electroweak parity-violating electron scattering (PVES) on 208Pb. We investigate PVES in nuclear mean field approach to allow the accurate extraction of the neutron skin thickness of 208Pb from the parity-violating asymmetry probed in the experiment. We demonstrate a high linear correlation between the parity-violating asymmetry and the neutron skin thickness in successful mean field forces as the best means to constrain the neutron skin of 208Pb from PREX, without assumptions on the neutron density shape. Continuation of the experiment with higher precision in the parity-violating asymmetry is motivated since the present method can support it to constrain the density slope of the nuclear symmetry energy to new accuracy.

preprint2011arXiv

Study of the neutron skin thickness of ${}^{208}$Pb in mean field models

We study whether the neutron skin thickness $Δr_{np}$ of ${}^{208}$Pb originates from the bulk or from the surface of the neutron and proton density distributions in mean field models. We find that the size of the bulk contribution to $Δr_{np}$ of ${}^{208}$Pb strongly depends on the slope of the nuclear symmetry energy, while the surface contribution does not. We note that most mean field models predict a neutron density for ${}^{208}$Pb between the halo and skin type limits. We investigate the dependence of parity- violating electron scattering at the kinematics of the PREX experiment on the shape of the nucleon densities predicted by the mean field models for ${}^{208}$Pb. We find an approximate formula for the parity-violating asymmetry in terms of the central radius and the surface diffuseness of the nucleon densities of ${}^{208}$Pb in these models.

preprint2010arXiv

Analysis of bulk and surface contributions in the neutron skin of nuclei

The neutron skin thickness of nuclei is a sensitive probe of the nuclear symmetry energy having multiple implications for nuclear and astrophysical studies. However, precision measurements of this observable are difficult. The analysis of the experimental data may imply some assumptions about the bulk or surface nature of the formation of the neutron skin. Here, we study the bulk or surface character of neutron skins of nuclei following from calculations with Gogny, Skyrme, and covariant nuclear mean-field interactions. These interactions are successful in describing nuclear charge radii and binding energies but predict different values for neutron skins. We perform the study by fitting two-parameter Fermi distributions to the calculated self-consistent neutron and proton densities. We note that the equivalent sharp radius is a more suitable reference quantity than the half-density radius parameter of the Fermi distributions to discern between the bulk and surface contributions in neutron skins. We present calculations for nuclei in the stability valley and for the isotopic chains of Sn and Pb.

preprint2010arXiv

Influence of the symmetry energy on the giant monopole resonance of neutron-rich nuclei

We analyze the influence of the density dependence of the symmetry energy on the average excitation energy of the isoscalar giant monopole resonance (GMR) in stable and exotic neutron-rich nuclei by applying the relativistic extended Thomas-Fermi method in scaling and constrained calculations. For the effective nuclear interaction, we employ the relativistic mean field model supplemented by an isoscalar-isovector meson coupling that allows one to modify the density dependence of the symmetry energy without compromising the success of the model for binding energies and charge radii. The semiclassical estimates of the average energy of the GMR are known to be in good agreement with the results obtained in full RPA calculations. The present analysis is performed along the Pb and Zr isotopic chains. In the scaling calculations, the excitation energy is larger when the symmetry energy is softer. The same happens in the constrained calculations for nuclei with small and moderate neutron excess. However, for nuclei of large isospin the constrained excitation energy becomes smaller in models having a soft symmetry energy. This effect is mainly due to the presence of loosely-bound outer neutrons in these isotopes. A sharp increase of the estimated width of the resonance is found in largely neutron-rich isotopes, even for heavy nuclei, which is enhanced when the symmetry energy of the model is soft. The results indicate that at large neutron numbers the structure of the low-energy region of the GMR strength distribution changes considerably with the density dependence of the nuclear symmetry energy, which may be worthy of further characterization in RPA calculations of the response function.

preprint2010arXiv

Origin of the neutron skin thickness of 208Pb in nuclear mean-field models

We study whether the neutron skin thickness (NST) of 208Pb originates from the bulk or from the surface of the nucleon density distributions, according to the mean-field models of nuclear structure, and find that it depends on the stiffness of the nuclear symmetry energy. The bulk contribution to NST arises from an extended sharp radius of neutrons, whereas the surface contribution arises from different widths of the neutron and proton surfaces. Nuclear models where the symmetry energy is stiff, as typical relativistic models, predict a bulk contribution in NST of 208Pb about twice as large as the surface contribution. In contrast, models with a soft symmetry energy like common nonrelativistic models predict that NST of 208Pb is divided similarly into bulk and surface parts. Indeed, if the symmetry energy is supersoft, the surface contribution becomes dominant. We note that the linear correlation of NST of 208Pb with the density derivative of the nuclear symmetry energy arises from the bulk part of NST. We also note that most models predict a mixed-type (between halo and skin) neutron distribution for 208Pb. Although the halo-type limit is actually found in the models with a supersoft symmetry energy, the skin-type limit is not supported by any mean-field model. Finally, we compute parity-violating electron scattering in the conditions of the 208Pb parity radius experiment (PREX) and obtain a pocket formula for the parity-violating asymmetry in terms of the parameters that characterize the shape of the 208Pb nucleon densities.

preprint2009arXiv

Microscopic-Macroscopic Approach for Binding Energies with Wigner-Kirkwood Method

The semi-classical Wigner-Kirkwood $\hbar$ expansion method is used to calculate shell corrections for spherical and deformed nuclei. The expansion is carried out up to fourth order in $\hbar$. A systematic study of Wigner-Kirkwood averaged energies is presented as a function of the deformation degrees of freedom. The shell corrections, along with the pairing energies obtained by using the Lipkin-Nogami scheme, are used in the microscopic-macroscopic approach to calculate binding energies. The macroscopic part is obtained from a liquid drop formula with six adjustable parameters. Considering a set of 367 spherical nuclei, the liquid drop parameters are adjusted to reproduce the experimental binding energies, which yields a {\it rms} deviation of 630 keV. It is shown that the proposed approach is indeed promising for the prediction of nuclear masses.

preprint2006arXiv

Average ground-state energy of finite Fermi systems

Semiclassical theories like the Thomas-Fermi and Wigner-Kirkwood methods give a good description of the smooth average part of the total energy of a Fermi gas in some external potential when the chemical potential is varied. However, in systems with a fixed number of particles N, these methods overbind the actual average of the quantum energy as N is varied. We describe a theory that accounts for this effect. Numerical illustrations are discussed for fermions trapped in a harmonic oscillator potential and in a hard wall cavity, and for self-consistent calculations of atomic nuclei. In the latter case, the influence of deformations on the average behavior of the energy is also considered.

preprint2006arXiv

Nuclear expansion with excitation

The expansion of an isolated hot spherical nucleus with excitation energy and its caloric curve are studied in a thermodynamic model with the SkM* force as the nuclear effective two-body interaction. The calculated results are shown to compare well with the recent experimental data from energetic nuclear collisions. The fluctuations in temperature and density are also studied. They are seen to build up very rapidly beyond an excitation energy of 9 MeV/u. Volume-conserving quadrupole deformation in addition to expansion indicates, however, nuclear disassembly above an excitation energy of 4 MeV/u

preprint2006arXiv

Thomas-Fermi theory for atomic nuclei revisited

The recently developed semiclassical variational Wigner-Kirkwood (VWK) approach is applied to finite nuclei using external potentials and self-consistent mean fields derived from Skyrme interactions and from relativistic mean field theory. VWK consists of the Thomas-Fermi part plus a pure, perturbative hbar^2 correction. In external potentials, VWK passes through the average of the quantal values of the accumulated level density and total energy as a function of the Fermi energy. However, there is a problem of overbinding when the energy per particle is displayed as a function of the particle number. The situation is analyzed comparing spherical and deformed harmonic oscillator potentials. In the self-consistent case, we show for Skyrme forces that VWK binding energies are very close to those obtained from extended Thomas-Fermi functionals of hbar^4 order, pointing to the rapid convergence of the VWK theory. This satisfying result, however, does not cure the overbinding problem, i.e., the semiclassical energies show more binding than they should. This feature is more pronounced in the case of Skyrme forces than with the relativistic mean field approach. However, even in the latter case the shell correction energy for e.g. 208Pb turns out to be only -6 MeV what is about a factor two or three off the generally accepted value. As an ad hoc remedy, increasing the kinetic energy by 2.5%, leads to shell correction energies well acceptable throughout the periodic table. The general importance of the present studies for other finite Fermi systems, self-bound or in external potentials, is pointed out.

preprint2005arXiv

Atomic Parity Non-Conservation, Neutron Radii, and Effective Field Theories of Nuclei

Accurately calibrated effective field theories are used to compute atomic parity non-conserving (APNC) observables. Although accurately calibrated, these effective field theories predict a large spread in the neutron skin of heavy nuclei. While the neutron skin is strongly correlated to a large number of physical observables, in this contribution we focus on its impact on new physics through APNC observables. The addition of an isoscalar-isovector coupling constant to the effective Lagrangian generates a wide range of values for the neutron skin of heavy nuclei without compromising the success of the model in reproducing well constrained nuclear observables. Earlier studies have suggested that the use of isotopic ratios of APNC observables may eliminate their sensitivity to atomic structure. This leaves nuclear structure uncertainties as the main impediment for identifying physics beyond the standard model. We establish that uncertainties in the neutron skin of heavy nuclei are at present too large to measure isotopic ratios to better than the 0.1% accuracy required to test the standard model. However, we argue that such uncertainties will be significantly reduced by the upcoming measurement of the neutron radius in 208Pb at the Jefferson Laboratory.

preprint2005arXiv

Sum Rule Approach to the Isoscalar Giant Monopole Resonance in Drip Line Nuclei

Using the density-dependent Hartree-Fock approximation and Skyrme forces together with the scaling method and constrained Hartree-Fock calculations, we obtain the average energies of the isoscalar giant monopole resonance. The calculations are done along several isotopic chains from the proton to the neutron drip lines. It is found that while approaching the neutron drip line, the scaled and the constrained energies decrease and the resonance width increases. Similar but smaller effects arise near the proton drip line, although only for the lighter isotopic chains. A qualitatively good agreement is found between our sum rule description and the presently existing random phase approximation results. The ability of the semiclassical approximations of the Thomas-Fermi type, which properly describe the average energy of the isoscalar giant monopole resonance for stable nuclei, to predict average properties for nuclei near the drip lines is also analyzed. We show that when hbar corrections are included, the semiclassical estimates reproduce, on average, the quantal excitation energies of the giant monopole resonance for nuclei with extreme isospin values.

preprint2004arXiv

Superheavy nuclei in relativistic effective Lagrangian model

Isotopic and isotonic chains of superheavy nuclei are analyzed to search for spherical double shell closures beyond Z=82 and N=126 within the new effective field theory model of Furnstahl, Serot, and Tang for the relativistic nuclear many-body problem. We take into account several indicators to identify the occurrence of possible shell closures, such as two-nucleon separation energies, two-nucleon shell gaps, average pairing gaps, and the shell correction energy. The effective Lagrangian model predicts N=172 and Z=120 and N=258 and Z=120 as spherical doubly magic superheavy nuclei, whereas N=184 and Z=114 show some magic character depending on the parameter set. The magicity of a particular neutron (proton) number in the analyzed mass region is found to depend on the number of protons (neutrons) present in the nucleus.

preprint2004arXiv

Versatility of field theory motivated nuclear effective Lagrangian approach

We analyze the results for infinite nuclear and neutron matter using the standard relativistic mean field model and its recent effective field theory motivated generalization. For the first time, we show quantitatively that the inclusion in the effective theory of vector meson self-interactions and scalar-vector cross-interactions explains naturally the recent experimental observations of the softness of the nuclear equation of state, without losing the advantages of the standard relativistic model for finite nuclei.

preprint2003arXiv

Semiclassical evaluation of average nuclear one and two body matrix elements

Thomas-Fermi theory is developed to evaluate nuclear matrix elements averaged on the energy shell, on the basis of independent particle Hamiltonians. One- and two-body matrix elements are compared with the quantal results and it is demonstrated that the semiclassical matrix elements, as function of energy, well pass through the average of the scattered quantum values. For the one-body matrix elements it is shown how the Thomas-Fermi approach can be projected on good parity and also on good angular momentum. For the two-body case the pairing matrix elements are considered explicitly.

preprint2002arXiv

Surface Incompressibility from Semiclassical Relativistic Mean Field Calculations

By using the scaling method and the Thomas-Fermi and Extended Thomas-Fermi approaches to Relativistic Mean Field Theory the surface contribution to the leptodermous expansion of the finite nuclei incompressibility has been self-consistently computed. The validity of the simplest expansion, which contains volume, volume-symmetry, surface and Coulomb terms, is examined by comparing it with self-consistent results of the finite nuclei incompressibility for some currently used non-linear sigma-omega parameter sets. A numerical estimate of higher-order contributions to the leptodermous expansion, namely the curvature and surface-symmetry terms, is made.

preprint2001arXiv

Pairing Properties In Relativistic Mean Field Models Obtained From Effective Field Theory

We apply recently developed effective field theory nuclear models in mean field approximation (parameter sets G1 and G2) to describe ground-state properties of nuclei from the valley of $β$-stability up to the drip lines. For faster calculations of open-shell nuclei we employ a modified BCS approach which takes into account quasi-bound levels owing to their centrifugal barrier, with a constant pairing strength. We test this simple prescription by comparing with available Hartree-plus-Bogoliubov results. Using the new effective parameter sets we then compute separation energies, density distributions and spin--orbit potentials in isotopic (isotonic) chains of nuclei with magic neutron (proton) numbers. The new forces describe the experimental systematics similarly to conventional non-linear $σ-ω$ relativistic force parameters like NL3.

preprint2001arXiv

Scaling Calculation of Isoscalar Giant Resonances in Relativistic Thomas-Fermi Theory

We derive analytical expressions for the excitation energy of the isoscalar giant monopole and quadrupole resonances in finite nuclei, by using the scaling method and the extended Thomas-Fermi approach to relativistic mean field theory. We study the ability of several non-linear sigma-omega parameter sets of common use in reproducing the experimental data. For monopole oscillations the calculations agree better with experiment when the nuclear matter incompressibility of the relativistic interaction lies in the range 220-260 MeV. The breathing-mode energies of the scaling method compare satisfactorily with those obtained in relativistic RPA and time-dependent mean field calculations. For quadrupole oscillations all the analyzed non-linear parameter sets reproduce the empirical trends reasonably well.

preprint2000arXiv

Ground-state properties and spins of the odd {\boldmath$Z=N+1$} nuclei {\boldmath$^{61}$}Ga{\boldmath$-$}{\boldmath$^{97}$}In}

Binding energies, quadrupole deformation parameters, spins and parities of the neutron-deficient odd $Z=N+1$ nuclei in the $A\sim 80$ region are calculated in the relativistic mean field approximation. The ground-state and low-lying configurations of the recently observed $^{77}$Y, $^{79}$Zr and $^{83}$Mo nuclei are analyzed. The calculated results are compared with other theoretical predictions.

preprint1998arXiv

Nuclear surface properties in relativistic effective field theory

We perform Hartree calculations of symmetric and asymmetric semi-infinite nuclear matter in the framework of relativistic models based on effective hadronic field theories as recently proposed in the literature. In addition to the conventional cubic and quartic scalar self-interactions, the extended models incorporate a quartic vector self-interaction, scalar-vector non-linearities and tensor couplings of the vector mesons. We investigate the implications of these terms on nuclear surface properties such as the surface energy coefficient, surface thickness, surface stiffness coefficient, neutron skin thickness and the spin-orbit force.