Search

The equation of state of hot neutron star matter of n+p+e+$\mu$ composition in $\beta$-equilibrium is studied for both neutrino-free isothermal and neutrino-trapped isentropic conditions, using the formalism where the thermal evolution is built upon its zero-temperature predictions in a self-consistent manner. The accuracy of the parabolic approximation, often used in the finite temperature calculation of hot neutron star matter, is verified by comparing with the results obtained from the exact evaluation in the neutrino-free neutron star matter. The equation of state of neutrino-trapped isentropic matter at low entropic condition, relevant to the core-collapsing supernovae, is formulated. In the isentropic matter, the particle fractions and equation of state have marginal variance as entropy per particle varies between 1 to 3 (in the unit of k$_B$), but the temperature profile shows marked variation. The isentropes are found to be much less sensitive to the nuclear matter incompressibility, but have a large dependence on the slope parameter L. The bulk properties of the neutron stars predicted by the isentropic equation of state for different entropy are calculated. A model calculation for the early stage evolution of the protoneutron star to neutron star configuration is also given., Comment: 30 pages, 11 figures, submitted to Journal of Physics G

The correlation between the charge radii differences in mirror nuclei pairs and the neutron skin thickness has been studied with the so-called finite range simple effective interaction over a wide mass region. The so far precisely measured charge radii difference data within their experimental uncertainty ranges in the 34Ar-34S, 36Ca-36S, 38Ca-38Ar, and 54Ni-54Fe mirror pairs are used to ascertain an upper limit for the slope parameter of the nuclear symmetry energy L $\approx$ 100 MeV. This limiting value of L is found to be consistent with the upper bound of the NICER PSR J0740+6620 constraint at 1$\sigma$ level for the radius R$_{1.4}$ of 1.4 M$_\odot$ neutron stars. The lower bound of the NICER R$_{1.4}$ data constrains the lower limit of L to $\approx$ 70 MeV. Within the range for L = 70-100 MeV the tidal deformability $\Lambda^{1.4}$ constraint, which is extracted from the GW170817 event at 2$\sigma$ level, and the recent PREX-2 and CREX data on the neutron skin thickness are discussed., Comment: 10 pages, 5 figures, 2 tables

In this work, we study the elastic scattering of some light particles, such as $^2$H, $^3$H, $^3$He, and $^4$He, by heavy target nuclei with an extended Watanabe model, which uses as input the neutron-nucleus and proton-nucleus optical potentials and the ground-state wave functions of the projectile. The nucleon-nucleus optical potential used in this work was obtained within a semi-microscopic nuclear matter approach, whose real and imaginary parts are provided by the first and second-order terms, respectively, of the Taylor expansion of the Brueckner-Hartree-Fock mass operator obtained with the reaction G-matrix built up with the Gogny force \cite{lopez21}. The angular distributions of the scattering of $^2$H, $^3$H, $^3$He, and $^4$He from different target nuclei and at a different incident energy of the projectile computed with this model are analyzed. The reaction cross-sections corresponding to some of these scattering processes are also calculated. Our results are compared with the experimental values as well as with another Watanabe calculation where the nucleon-nucleus the optical potential is provided by the phenomenological K\"oning-Delaroche model. The limitations of the extended Watanabe model used in this work are also discussed.

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(\rho_0)$. The estimated crossing for the neutron number $N$=46 could be reproduced by the equation of state corresponding to $K(\rho_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.

Recent experiments performed in neutron-rich copper isotopes have revealed a crossing in the nucleus $^{75}$Cu between the $3/2^-$ and $5/2^-$ levels, which correspond to the ground-state and the first excited state in isotopes with mass number below $A =75$. Due to the strong single-particle character of these states, this scenario can be investigated through the analysis of the proton spectrum provided by mean-field models in nickel isotopes with neutron numbers between $N$=40 and $N$=50. In this work we show that the aforementioned crossing is mainly driven by the mean-field provided by the effective nucleon-nucleon and spin-orbit interactions. We also analyze the impact of the tensor interaction, and find that in some mean-field models it is essential to reproduce the crossing of the 2$p_{3/2}$ and 1$f_{5/2}$ proton single-particle levels, as in the case of the SAMi-T Skyrme force and the D1M Gogny interaction, whereas in other cases, as for example the SLy5 Skyrme force, a reasonable tensor force appears to be unable to modify the mean-field enough to reproduce this level crossing. Finally, in the calculations performed with the so-called simple effective interaction (SEI), it is shown that the experimental data in nickel and copper isotopes considered in this work can be explained satisfactorily without any explicit consideration of the tensor interaction., Comment: 7 pages, 3 figures, 1 table

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., Comment: 17 pages, 8 Figures and 2 Tables

Analytic average mean-field potentials of the Fermi-function (Woods-Saxon) type for the whole nuclear chart with space-dependent effective mass are deduced from the D1S Gogny force. Those ready-for-use potentials are advertised as an alternative to other existing phenomenological mean-field potentials., Comment: 14 pages, 5 Figures, 2 Tables

We study the superfluid dynamics of the outer core of neutron stars by means of a hydrodynamic model made of a neutronic superfluid and a protonic superconductor, coupled by both the dynamic entrainment and the Skyrme SLy4 nucleon-nucleon interactions. The resulting nonlinear equations of motion are probed in the search for dynamical instabilities triggered by the relative motion of the superfluids that could be related to observed timing anomalies in pulsars. Through linear analysis, the origin and expected growth of the instabilities is explored for varying nuclear-matter density. Differently from previous findings, the dispersion of linear excitations in our model shows rotonic structures below the pair-breaking energy threshold, which lies at the origin of the dynamical instabilities, and could eventually lead to emergent vorticity along with modulations of the superfluid density., Comment: 10 pages, 5 figures. v2: minor changes

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 $\beta$-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(\rho_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., Comment: 48 pages, 19 figures, 5 tables

Average, global pairing behaviours of electronic devices, like films, wires, and grains are studied with semiclassical methods, as Weyl and Thomas-Fermi approximations, in the ultrasmall, i.e., quantal regime, which exhibits strong quantum fluctuations and shell effects. This is the case also for superfluid nuclei. Results, mostly analytic, are elaborated for average size dependences, ready for easy use also in other circumstances. Comparisons with experimental results are given where possible.

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., Comment: 19 pages, 9 figures, Fig 9 updated and Tables are revised for Pressure

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., Comment: 7 pages, 4 figures. Accepted for publication in Phys. Lett. B

The properties of the core-crust transition in neutron stars are investigated using effective nuclear forces of finite-range. Special attention is paid to the so-called dynamical method for locating the transition point, which, apart from the stability of the uniform nuclear matter against clusterization, also considers contributions due to finite-size effects. In particular, contributions to the transition density and pressure from the direct and exchange energies are carefully analyzed. To this end, finite-range forces of Gogny, Modified Gogny Interaction (MDI) and Simple Effective Interaction (SEI) types are used in the numerical applications. The results from the dynamical approach are compared with those from the popular thermodynamical method that neglects the surface and Coulomb effects in the stability condition. The dependence of the core-crust transition on the stiffness of the symmetry energy of the finite-range models is also addressed. Finally, we analyze the impact of the transition point on the mass, thickness and fraction of the moment of inertia of the neutron star crust. Prominent differences in these crustal properties of the star are found between using the transition point obtained with the dynamical method or the thermodynamical method. It is concluded that the core-crust transition needs to be ascertained as precisely as possible in order to have realistic estimates of the observed phenomena where the crust plays a significant role.

The most successful Gogny interactions of the D1 family, namely D1S, D1N and D1M, suffer the common problem of a too soft neutron matter equation of state at high density, which prevents them from predicting a maximal mass of neutron stars of two solar masses, as required by recent astronomical observations. To cure this deficiency, we have proposed recently a reparametrization of the D1M force by fine tuning the slope of the symmetry energy in such a way that it preserves the ground-state properties of D1M in finite nuclei and also describes successfully the global properties of neutron stars, in particular its maximal mass, in consonance with the observational data. In this contribution we revisit this reparametrization by discussing two modified Gogny forces, dubbed D1M$^*$ and D1M$^{**}$., Comment: 8 pages, 3 figures, 2 tables Proceedings for 37 International Workshop on Nuclear Theory, Rila Mountains (Bulgaria) June 2018

The stacks of Bose-Einstein condensates coupled by long Josephson junctions present a rich phenomenology feasible to experimental realization and specially suitable for technological applications as the nonlinear-optics and superconducting analogues have already proved. Among this, we show that transverse Bloch waves excited in arrays of one-dimensional coupled condensates can carry tunneling superflows whose dynamical stability depends on the quasimomentum. Across the stacks with periodic boundary conditions, forming closed ring-shaped systems, such Bloch states yield transverse Josephson vortices with a generic non-integer circulation in units of $h/m$. Additionally, the superpositions of degenerate linear Bloch waves can suppress the supercurrents and give rise to families of nonlinear standing-wave states with strong (transverse) spatial localization. Stable states of this type can also be found in finite size systems., Comment: 14 pages, 8 figures

The conclusions of the manuscript 1806.02080v1 questioning the adequacy of the recently proposed Gogny D1M* interaction for finite nuclei calculations using harmonic oscillator (HO) basis are revised. Several convergence and stability studies are performed with HO basis of different sizes and oscillator parameters and the results show the robustness of the D1M* results for finite nuclei. This analysis is also extended to beyond mean-field calculations of generator-coordinate-method type with D1M*. On the other hand, the existence of a finite-size instability in finite nuclei when coordinate space methods are used to solve the HF equations (as shown in 1806.02080v1) is independently confirmed for D1M* using an in-house computer code based on a quasilocal approximation to the HF exchange potential. We confirm that the most affected quantity in the coordinate space calculation is the spatial density at the origin, but integrated quantities like binding energies or radii show a plateau against the number of iterations, where they are consistent with the values from the HO basis calculation, before diverging for a larger number of iterations. A connection between the last occupied s-orbital in the nucleus and the appearance of instabilities in coordinate space is observed., Comment: 12 pages, 22 figures

The characteristic physical properties of rotating neutron stars under the r-mode oscillation are evaluated using the finite-range simple effective interaction. Emphasis is given on examining the influence of the stiffness of both the symmetric and asymmetric parts of the nuclear equation of state on these properties. The amplitude of the r-mode at saturation is calculated using the data of particular neutron stars from the considerations of "spin equilibrium" and "thermal equilibrium". The upper limit of the r-mode saturation amplitude is found to lie in the range 10^{-8}-10^{-6}, in agreement with the predictions of earlier work., Comment: 28 pages, 7 figures, 5 tables

The D1 family of parametrizations of the Gogny interaction commonly suffers from a rather soft neutron matter equation of state that leads to maximal masses of neutron stars well below the observational value of two solar masses. We propose a reparametrization scheme that preserves the good properties of the Gogny force but allows one to tune the density dependence of the symmetry energy, which, in turn, modifies the predictions for the maximum stellar mass. The scheme works well for D1M, and leads to a new parameter set, dubbed D1M*. In the neutron-star domain, D1M* predicts a maximal mass of two solar masses and global properties of the star in harmony with those obtained with the SLy4 Skyrme interaction. By means of a set of selected calculations in finite nuclei, we check that D1M* performs comparably well to D1M in several aspects of nuclear structure in nuclei., Comment: 6 pages, 5 figures

Relations between the nuclear symmetry energy coefficient and its density derivatives are derived. The relations hold for a class of interactions with quadratic momentum dependence and a power-law density dependence. The structural connection between the different symmetry energy elements as obtained seems to be followed by almost all reasonable nuclear energy density functionals, both relativistic and non-relativistic, suggesting a universality in the correlation structure. This, coupled with known values of some well-accepted constants related to nuclear matter, helps in constraining values of different density derivatives of the nuclear symmetry energy shedding light on the isovector part of the nuclear interaction., Comment: 6 pages, 2 figures

We study the symmetry energy and the core-crust transition in neutron stars using the finite-range Gogny nuclear interaction and examine the deduced crustal thickness and crustal moment of inertia. We start by analyzing the second-, fourth- and sixth-order coefficients of the Taylor expansion of the energy per particle in powers of the isospin asymmetry for Gogny forces. These coefficients provide information about the departure of the symmetry energy from the widely used parabolic law. The neutron star core-crust transition is evaluated by looking at the onset of thermodynamical instability of the liquid core. The calculation is performed with the exact (i.e., without Taylor expansion) Gogny EoS for the core, and also with its Taylor expansion in order to assess the influence of isospin expansions on locating the inner edge of neutron star crusts. It is found that the properties of the core-crust transition derived from the exact EoS differ from the predictions of the Taylor expansion even when the expansion is carried through sixth order in the isospin asymmetry. Gogny forces, using the exact EoS, predict the ranges $0.094 \text{ fm}^{-3} \lesssim \rho_t \lesssim 0.118\text{ fm}^{-3}$ for the transition density and $0.339 \text{ MeV fm}^{-3} \lesssim P_t \lesssim 0.665 \text{ MeV fm}^{-3}$ for the transition pressure. The transition densities show an anticorrelation with the slope parameter $L$ of the symmetry energy. The transition pressures are not found to correlate with $L$. Neutron stars obtained with Gogny forces have maximum masses below $1.74M_\odot$ and relatively small moments of inertia. The crustal mass and moment of inertia are evaluated and comparisons are made with the constraints from observed glitches in pulsars., Comment: 24 pages, 15 figures, discussions and bibliography updated, to appear in Physical Review C

Deformed and spherical even-even nuclei are studied using a finite range simple effective interaction within the Hartree-Fock-Bogoliubov mean field approach. Different parameter sets of the interaction, corresponding to different incompressibility, are constructed by varying the exponent gamma of the density in the traditional density-dependent term. Ten of the twelve parameters of these interactions are determined from properties of asymmetric nuclear matter and spin polarized pure neutron matter. The two remaining parameters are fitted to reproduce the experimental binding energies known in 620 even-even nuclei using several variants of the rotational energy correction. The rms deviations for the binding energy depend on the value of gamma and the way the rotational energy correction is treated but they can be as low as 1.56 MeV, a value competitive with other renowned effective interactions of Skyrme and Gogny type. Charge radii are compared to the experimental values of 313 even-even nuclei and the rms deviation is again comparable and even superior to the one of popular Skyrme and Gogny forces. Emphasis is given to the deformation properties predicted with these interactions by analyzing the Potential Energy Surfaces for several well deformed nuclei and the fission barriers of some nuclei. Comparison of the results with the experimental information, where available, as well as with the results of the Gogny D1S force shows satisfactory agreement., Comment: 20 pages, 6 figures

In this paper we review the semiclassical extended Thomas-Fermi theory for describing the ground-state properties of nuclei. The binding energies calculated in this approach do not contain shell effects and, in this sense, they are analogous to those obtained from the mass formula. We discuss some techniques for incorporating the shell effects which are missing in the semiclassical calculation, such as the so-called expectation value method and the Kohn-Sham scheme. We present numerical applications for effective zero-range Skyrme forces and finite-range Gogny forces., Comment: Contribution to the book "Quantum Theory: in honour of Vladimir A. Fock" pt. 2, eds. Y. Novozhilov, V. Novozhilov, Publishing group of St. Petersburg University, 1998, p. 118

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., Comment: 35 pages including 4 Tables & 12 figures; The article has been accepted for publication in J.Phys. G

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(\rho)$ 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., Comment: 10 pages including 6 figures

The Barcelona-Catania-Paris-Madrid (BCPM) functional recently proposed to describe nuclear structure properties of finite nuclei is generalized as to include a realistic effective mass. The resulting functional is as good as the previous one in describing binding energies, radii, deformation properties, etc and, in addition, the description of Giant Quadrupole Resonance energies is greatly improved., Comment: 8 pages, 5 figures

The information on the symmetry energy and its density dependence is deduced by comparing the available data on the electric dipole polarizability $\alpha_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 $\alpha_D$ are used to validate different correlations involving $\alpha_D$, the symmetry energy at the saturation density $J$, the corresponding slope parameter $L$ and the neutron skin thickness $\Delta 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 $\Delta r_{\!np}$. The resulting intervals are: $J\!=\!30 \text{-}35$ MeV, $L\!=\!20 \text{-} 66$ MeV; and the values for $\Delta 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., Comment: 12 pages, 2 tables and 6 figures

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\mu$ 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., Comment: 23 pages, 17 figures, revised version accepted for publication in Astronomy & Astrophysics

Using HF+BCS method we study light nuclei with nuclear charge in the range $2 \leq Z \leq 8$ and lying near the neutron drip line. The HF method uses effective Skyrme forces and allows for axial deformations. We find that the neutron drip line forms stability peninsulas at $^{18}$He and $^{40}$C. These isotopes are found to be stable against one neutron emission and possess the highest known neutron to proton ratio in stable nuclei.

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., Comment: 31 pages, 13 figures, to appear in Journal of Physics G

Within the context of Hartree-Fock-Bogoliubov theory, we study the behavior of superfluid Fermi systems when they pass from a small to a large container. Such systems can be now realized thanks to recent progress in experimental techniques. It will allow to better understand pairing properties at overflow and in general in rapidly varying external potentials.

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., Comment: 17 pages, 14 figures

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., Comment: 12 pages, including 11 figures, appearing in special issue of EPJ-A on Nuclear Symmetry Energy

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., Comment: 14 pages, 6 figures, accepted version

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., Comment: 17 pages, 10 figures, submitted to EPJA special volume on Nuclear Symmetry Energy

We study the electric dipole polarizability $\alpha_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 $\alpha_D$ and other isovector observables. Insights from the droplet model suggest that the product of $\alpha_D$ and the nuclear symmetry energy at saturation density $J$ is much better correlated with the neutron skin thickness $\Delta r_{np}$ of ${}^{208}$Pb than the polarizability alone. Correlations of $\alpha_D J$ with $\Delta r_{np}$ and with the symmetry energy slope parameter $L$ suggest that $\alpha_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 $\alpha_D$ and the parity-violating asymmetry in ${}^{208}$Pb. We find that the recent experimental determination of $\alpha_D$ in ${}^{208}$Pb in combination with the range for the symmetry energy at saturation density $J=[31\pm (2)_{\rm est.}]$\,MeV suggests $\Delta 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.

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)]., Comment: Proceedings - International Nuclear Physics Conference (INPC), Firenze 2 - 7 June 2013

Exotic and drip-line nuclei as well as nuclei immersed in a low density gas of neutrons in the outer crust of neutron stars are systematically investigated with respect to their neutron pairing properties. This is done using Skyrme density-functional and different pairing forces such as a density-dependent contact interaction and a separable form of a finite-range Gogny interaction. Hartree-Fock-Bogoliubov and BCS theories are compared. It is found that neutron pairing is reduced towards the drip line while overcast by strong shell effects. Furthermore resonances in the continuum can have an important effect counterbalancing the tendency of reduction and leading to a persistence of pairing at the drip line. It is also shown that in these systems the difference between HFB and BCS approaches can be qualitatively large.

The mean field properties and equation of state for asymmetric nuclear matter are studied by using a simple effective interaction which has a single finite range Gaussian term. The study of finite nuclei with this effective interaction is done by means of constructing a quasilocal energy density functional for which the single particle equations take the form of Skryme-Hartree-Fock equations. The predictions of binding energies and charge radii of spherical nuclei are found to be compatible with the results of standard models as well as experimental data.

Given the promising features of the recently proposed Barcelona-Catania-Paris (BCP) functional \cite{Baldo.08}, it is the purpose of this paper to still improve on it. It is, for instance, shown that the number of open parameters can be reduced from 4-5 to 2-3, i.e. by practically a factor of two. One parameter is tightly fixed by a fine-tuning of the bulk, a second by the surface energy. The third is the strength of the spin-orbit potential on which the final result does not depend within the scatter of the values used in Skyrme and Gogny like functionals. An energy rms value of 1.58 MeV is obtained from a fit of these three parameters to the 579 measured masses reported in the Audi and Waspra 2003 compilation. This rms value compares favorably with the one obtained using other successful mean field theories. Charge radii are also well reproduced when compared with experiment. The energies of some excited states, mostly the isoscalar giant monopole resonances, are studied within this model as well., Comment: 23 pages, 12 figures

The momentum and density dependence of mean fields in symmetric and asymmetric nuclear matter are analysed using the simple density dependent finite range effective interaction containing a single Gaussian term alongwith the zero-range terms. Within the formalism developed, it is possible to reproduce the various diverging predictions on the momentum and density dependence of isovector part of the mean field in asymmetric matter. The finite nucleus calculation is formulated for the simple Gaussian interaction in the framework of quasilocal density functional theory. The prediction of energies and charge radii of the interaction for the spherical nuclei compares well with the results of other effective theories., Comment: 10 pages, 3 figures, 2 tables, To appear in the Proceedings of the 11th International Conference on Nucleus-Nucleus Collisions (NN2012), May 27-June 1, 2012, San Antonio, Texas, USA

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., Comment: 12 pages, 11 figures, Phys. Rev.C (in press)

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., Comment: 13 pages; 19 figures

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\"oller 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., Comment: Sequel to A. Bhagwat, X. Vi\~nas, M. Centelles, P. Schuck and R. Wyss, Phys. Rev. C 81, 044321 (2010)

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., Comment: 12 pages, 5 figures, Presented at XXXII Mazurian Lakes Conference on Physics, Piaski, Poland, September 11-18, 2011

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., Comment: Contribution to the 2nd Iberian Nuclear Astrophysics Meeting on Compact Stars proceedings

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$\delta$ is presented which includes four mesons $\sigma$, $\omega$, $\delta$, and $\rho$ 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 $\delta$-meson and its density dependence., Comment: Published Phys. Rev. C (significant digits added to Table I)

Two issues are treated in this work: (i) the generic fact that if a fermionic superfluid in the BCS regime overflows from a narrow container into a much wider one, pairing is much suppressed at the overflow point. Physical examples where this feature may play an important role are discussed (ii) A Thomas-Fermi (TF) approach to inhomogeneous superfluid Fermi-systems is presented and shown that it works well in cases where the Local Density Approximation (LDA) breaks down., Comment: 4 pages, 5 figures

Using HF+BCS method with Skyrme forces we analyze the neutron drip line. It is shown that around magic and new magic numbers the drip line may form stability peninsulas. It is shown the location of these peninsulas does not depend on the choice of Skyrme forces. It is found that the size of the peninsulas is sensitive to the choice of Skyrme forces and the most extended peninsulas appear with the SkI2 set.