Structure and composition of the 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 -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

Woods-Saxon type of mean-field potentials with effective mass derived from the D1S Gogny force

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

Microscopic-macroscopic approach for ground-state energies based on the Gogny force with the Wigner-Kirkwood averaging scheme

In the previous paper [Bhagwat, Centelles, Viñas, and Schuck, Phys. Rev. C 103, 024320 (2021)] 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 paper using the Gogny D1S force within the ETF scheme. The parametrized ETF potentials and densities of the previous paper 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

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

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

Youthful Familicidal Offenders: Targeted Victims, Planned Attacks

A nonrandom national sample of 16 familicides, which involved 19 offenders (ages 14 to 21 years) who either killed or made a serious attempt to kill their families, was studied. The majority of offenders were Caucasian (78.91 %) males (84.21 %) with interpersonal family conflicts due to parental control, substance use, or physical violence. Prior to the murders, 50 % of the offenders reported to others their intent to kill their families. All of the 42 reported victims were specifically targeted and most of the homicides were planned shooting attacks (75 %) rather than spontaneous eruptions. Immediately following the homicides, 75 % of the offenders stole money from their families, and in 50 % of the cases they either called their friends to report the murders or to plan leisure activities. All offenders were immediate suspects and 81.25 % confessed to the homicides. Implications for furthering our understanding of this group of young offenders are offered

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 calculations using the Argonne v18 potential plus three-body forces computed with the Urbana model. To deal with the inhomogeneous structures of matter in the NS crust, we use a recent nuclear energy density functional that is directly based on the same microscopic 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 new 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 0.067 fm−3 and 0.0825 fm−3, where the transition to the core takes place. The NS core is computed from the new nuclear EoS assuming non-exotic constituents (core of npeμ matter). In each region of the star, we discuss the comparison of the new EoS with previous EoSs for the complete NS structure, widely used in astrophysical calculations. The new microscopically derived EoS fulfills at the same time a NS maximum mass of 2 M with a radius of 10 km, and a 1.5 M NS with a radius of 11.6 km

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 initiocalculations 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 nonrelativistic and relativistic Brueckner calculations of symmetric and asymmetric nuclear matter. The isovector effective mass mp*−mn* derived from relativistic Brueckner theory is used to determine the coupling strength of the δ meson and its density dependence

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

Unified equation of state for neutron stars based on the Gogny interaction

t: The effective Gogny interactions of the D1 family were established by D. Gogny more than forty years ago with the aim to describe simultaneously the mean field and the pairing field corresponding to the nuclear interaction. The most popular Gogny parametrizations, namely D1S, D1N and D1M, describe accurately the ground-state properties of spherical and deformed finite nuclei all across the mass table obtained with Hartree-Fock-Bogoliubov (HFB) calculations. However, these forces produce a rather soft equation of state (EoS) in neutron matter, which leads to predict maximum masses of neutron stars well below the observed value of two solar masses. To remove this limitation, we built new Gogny parametrizations by modifying the density dependence of the symmetry energy predicted by the force in such a way that they can be applied to the neutron star domain and can also reproduce the properties of finite nuclei as good as their predecessors. These new parametrizations allow us to obtain stiffer EoS's based on the Gogny interactions, which predict maximum masses of neutron stars around two solar masses. Moreover, other global properties of the star, such as the moment of inertia and the tidal deformability, are in harmony with those obtained with other well tested EoSs based on the SLy4 Skyrme force or the Barcelona-Catania-Paris-Madrid (BCPM) energy density functional. Properties of the core-crust transition predicted by these Gogny EoSs are also analyzed. Using these new Gogny forces, the EoS in the inner crust is obtained with the Wigner-Seitz approximation in the Variational Wigner-Kirkwood approach along with the Strutinsky integral method, which allows one to estimate in a perturbative way the proton shell and pairing corrections. For the outer crust, the EoS is determined basically by the nuclear masses, which are taken from the experiments, wherever they are available, or by HFB calculations performed with these new forces if the experimental masses are not known