Showing 1–50 of 50 results for author: Gulminelli, F

Inference of neutron-star properties with unified crust-core equations of state for parameter estimation

Authors: P. J. Davis, H. Dinh Thi, A. F. Fantina, F. Gulminelli, M. Oertel, L. Suleiman

Abstract: Relating different global neutron-star (NS) properties, such as tidal deformability and radius, or mass and radius, requires an equation of state (EoS). Determining the NS EoS is therefore not only the science goal of a variety of observational projects, but it also enters in the analysis process; for example, to predict a NS radius from a measured tidal deformability via gravitational waves (GW)… ▽ More Relating different global neutron-star (NS) properties, such as tidal deformability and radius, or mass and radius, requires an equation of state (EoS). Determining the NS EoS is therefore not only the science goal of a variety of observational projects, but it also enters in the analysis process; for example, to predict a NS radius from a measured tidal deformability via gravitational waves (GW) during the inspiral of a binary NS merger. To this aim, it is important to estimate the theoretical uncertainties on the EoS, one of which is the possible bias coming from an inconsistent treatment of the low-density region; that is, the use of a so called non-unified NS crust. We propose a numerical tool allowing the user to consistently match a nuclear-physics informed crust to an arbitrary high-density EoS describing the core of the star. We introduce an inversion procedure of the EoS close to saturation density that allows users to extract nuclear-matter parameters and extend the EoS to lower densities in a consistent way. For the treatment of inhomogeneous matter in the crust, a standard approach based on the compressible liquid-drop (CLD) model approach was used in our work. A Bayesian analysis using a parametric agnostic EoS representation in the high-density region is also presented in order to quantify the uncertainties induced by an inconsistent treatment of the crust. We show that the use of a fixed, realistic-but-inconsistent model for the crust causes small but avoidable errors in the estimation of global NS properties and leads to an underestimation of the uncertainties in the inference of NS properties. Our results highlight the importance of employing a consistent EoS in inference schemes. The numerical tool that we developed to reconstruct such a thermodynamically consistent EoS, CUTER, has been tested and validated for use by the astrophysical community. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 16 pages, 11 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 687, A44 (2024)

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

General predictions of neutron star properties using unified relativistic mean-field equations of state

Authors: Luigi Scurto, Helena Pais, Francesca Gulminelli

Abstract: In this work we present general predictions for the static observables of neutron stars (NSs) under the hypothesis of a purely nucleonic composition of the ultra-dense baryonic matter, using Bayesian inference on a very large parameter space conditioned by both astrophysical and nuclear physics constraints. The equation of states are obtained using a unified approach of the NS core and inner crust… ▽ More In this work we present general predictions for the static observables of neutron stars (NSs) under the hypothesis of a purely nucleonic composition of the ultra-dense baryonic matter, using Bayesian inference on a very large parameter space conditioned by both astrophysical and nuclear physics constraints. The equation of states are obtained using a unified approach of the NS core and inner crust within a fully covariant treatment based on a relativistic mean-field Lagrangian density with density dependent couplings. The posterior distributions are well compatible with the ones obtained by semi-agnostic meta-modelling techniques based on non-relativistic functionals, that span a similar portion of the parameter space in terms of nuclear matter parameters, and we confirm that the hypothesis of a purely nucleonic composition is compatible with all the present observations. We additionally show that present observations do not exclude the existence of very massive neutron stars with mass compatible with the lighter partner of the gravitational event GW190814 measured by the LIGO-Virgo collaboration. Some selected representative models, that respect well all the constraints taken into account in this study, and approximately cover the residual uncertainty in our posterior distributions, will be uploaded in the CompOSE database for use by the community. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 28 pages, 21 figures, 9 tables. Submitted to Phys. Rev. D

Nuclear physics inputs for dense-matter modelling in neutron stars. The nuclear equation of state

Authors: A. F. Fantina, F. Gulminelli

Abstract: In this contribution, we briefly present the equation-of-state modelling for application to neutron stars and discuss current constraints coming from nuclear physics theory and experiments. To assess the impact of model uncertainties, we employ a nucleonic meta-modelling approach and perform a Bayesian analysis to generate posterior distributions for the equation of state with filters accounting f… ▽ More In this contribution, we briefly present the equation-of-state modelling for application to neutron stars and discuss current constraints coming from nuclear physics theory and experiments. To assess the impact of model uncertainties, we employ a nucleonic meta-modelling approach and perform a Bayesian analysis to generate posterior distributions for the equation of state with filters accounting for both our present low-density nuclear physics knowledge and high-density neutron-star physics constraints. The global structure of neutron stars thus predicted is discussed in connection with recent astrophysical observations. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 6 pages, 4 figures. Contribution to the proceedings of the 28th International Nuclear Physics Conference (INPC 2022)

Journal ref: J. Phys.: Conf. Ser. 2586, 012112 (2023)

Light clusters in the liquid proto-neutron star inner crust

Authors: H. Dinh Thi, A. F. Fantina, F. Gulminelli

Abstract: Being born hot from core-collapse supernova, the crust of the proto-neutron star is expected to be made of a Coulomb liquid and composed of an ensemble of different nuclear species. In this work, we study the beta-equilibrated proto-neutron-star crust in the liquid phase in a self-consistent multi-component approach, employing a compressible liquid-drop description of the ions including the ion ce… ▽ More Being born hot from core-collapse supernova, the crust of the proto-neutron star is expected to be made of a Coulomb liquid and composed of an ensemble of different nuclear species. In this work, we study the beta-equilibrated proto-neutron-star crust in the liquid phase in a self-consistent multi-component approach, employing a compressible liquid-drop description of the ions including the ion centre-of-mass motion. Particular care is also devoted to the calculation of the rearrangement term, thus ensuring thermodynamic consistency. We compare the results of the multi-component plasma calculations with those obtained within a one-component (single-nucleus) approach, showing that important differences arise between the predictions of the two treatments. In particular, the abundances of helium clusters become important using a complete multi-component plasma approach, and eventually dominate the whole distribution at higher temperature in the crust. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: Submitted to the European Physical Journal A (EPJA) for the Topical Collection "The Nuclear Many-Body Problem"

Journal ref: Eur. Phys. J. A 59, 292 (2023)

The proto-neutron star inner crust in a multi-component plasma approach

Authors: H. Dinh Thi, A. F. Fantina, F. Gulminelli

Abstract: Proto-neutron stars (PNS) are born hot, with temperatures exceeding a few times $10^{10}$ K. In these conditions, the PNS crust is expected to be made of a Coulomb liquid composed of an ensemble of different nuclear species. We perform a study of the beta-equilibrated PNS crust in the liquid phase in a self-consistent multi-component plasma (MCP) approach, thus allowing us to consistently calculat… ▽ More Proto-neutron stars (PNS) are born hot, with temperatures exceeding a few times $10^{10}$ K. In these conditions, the PNS crust is expected to be made of a Coulomb liquid composed of an ensemble of different nuclear species. We perform a study of the beta-equilibrated PNS crust in the liquid phase in a self-consistent multi-component plasma (MCP) approach, thus allowing us to consistently calculate the impurity parameter, often taken as a free parameter in cooling simulations. We developed a self-consistent MCP approach at finite temperature using a compressible liquid-drop description of the ions, with surface parameters adjusted to reproduce experimental masses. The treatment of the ion centre-of-mass motion was included through a translational free-energy term accounting for in-medium effects. The results of self-consistent MCP calculations are systematically compared with those performed in a perturbative and in the one-component plasma treatment. We show that the inclusion of non-linear mixing terms arising from the ion centre-of-mass motion leads to a breakdown of the ensemble equivalence between the one-component and MCP approach. Our findings illustrate that the abundance of light nuclei becomes important, eventually dominating the distribution at higher density and temperature. This is reflected in the impurity parameter, which, in turn, may have a potential impact on NS cooling. For practical applications, we also provide a fitting formula for the impurity parameter in the PNS inner crust. Our results obtained within a self-consistent MCP approach show important differences in the prediction of the PNS composition with respect to those obtained with a one-component or a perturbative MCP approximation, particularly in the deeper region of the crust. This highlights the importance of a full, self-consistent MCP calculation for reliable predictions of the PNS crust composition. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 16 pages, 15 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 677, A174 (2023)

Nuclear physics constraints from binary neutron star mergers in the Einstein Telescope era

Authors: Francesco Iacovelli, Michele Mancarella, Chiranjib Mondal, Anna Puecher, Tim Dietrich, Francesca Gulminelli, Michele Maggiore, Micaela Oertel

Abstract: The next generation of ground-based gravitational-wave detectors, Einstein Telescope (ET) and Cosmic Explorer (CE), present a unique opportunity to put constraints on dense matter, among many other groundbreaking scientific goals. In a recent study the science case of ET was further strengthened, studying in particular the performances of different detector designs. In this paper we present a more… ▽ More The next generation of ground-based gravitational-wave detectors, Einstein Telescope (ET) and Cosmic Explorer (CE), present a unique opportunity to put constraints on dense matter, among many other groundbreaking scientific goals. In a recent study the science case of ET was further strengthened, studying in particular the performances of different detector designs. In this paper we present a more detailed study of the nuclear physics section of that work. In particular, focusing on two different detector configurations (the single-site triangular-shaped design and a design consisting of two widely separated "L-shaped" interferometers), we study the detection prospects of binary neutron star (BNS) mergers, and how they can reshape our understanding of the underlying equation of state (EoS) of dense matter. We employ several state-of-the-art EoS models and state-of-the-art synthetic BNS merger catalogs, and we make use of the Fisher information formalism (FIM) to quantify statistical errors on the astrophysical parameters describing individual BNS events. To check the reliability of the FIM method, we further perform a full parameter estimation for a few simulated events. Based on the uncertainties on the tidal deformabilities associated to these events, we outline a mechanism to extract the underlying injected EoS using a recently developed meta-modelling approach within a Bayesian framework. Our results suggest that with $\gtrsim 500$ events with signal-to-noise ratio greater than $12$, we will be able to pin down very precisely the underlying EoS governing the neutron star matter. △ Less

Submitted 10 November, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: 18 pages, 10 figures. v2 matches the version accepted for publication in PRD

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Generalised description of Neutron Star matter with nucleonic Relativistic Density Functional

Authors: Prasanta Char, Chiranjib Mondal, Francesca Gulminelli, Micaela Oertel

Abstract: In this work, we propose a meta-modelling technique to nuclear matter on the basis of a relativistic density functional with density-dependent couplings. Identical density dependence for the couplings both in the isoscalar and isovector sectors is employed. We vary the coupling parameters of the model to capture the uncertainties of the empirical nuclear matter parameters at saturation. Then, we c… ▽ More In this work, we propose a meta-modelling technique to nuclear matter on the basis of a relativistic density functional with density-dependent couplings. Identical density dependence for the couplings both in the isoscalar and isovector sectors is employed. We vary the coupling parameters of the model to capture the uncertainties of the empirical nuclear matter parameters at saturation. Then, we construct a large ensemble of unified equations of state in a consistent manner both for clusterized and uniform matter in $β$-equilibrium at zero temperature. Finally, we calculate neutron star properties to check the consistency with astrophysical observations within a Bayesian framework. Out of the different sets of astrophysical data employed, constraint on tidal deformability from the GW170817 event was found to be the most stringent in the posteriors of different neutron star properties explored in the present study. We demonstrate in detail the impact of the isovector incompressibility ($K_{sym}$) on high-density matter that leads to a considerable variation in the composition of neutron star matter. A couple of selected models with extreme values of $K_{sym}$, which satisfy various modern nuclear physics and neutron star astrophysics constraints, are uploaded in the \textsc{CompOSE} \cite{Typel:2013rza} database for use by the community. △ Less

Submitted 30 October, 2023; v1 submitted 23 July, 2023; originally announced July 2023.

Detectability of a phase transition in neutron star matter with third-generation gravitational wave interferometers

Authors: C. Mondal, M. Antonelli, F. Gulminelli, M. Mancini, J. Novak, M. Oertel

Abstract: Possible strong first-order hadron-quark phase transitions in neutron star interiors leave an imprint on gravitational waves, which could be detected with planned third-generation interferometers. Given a signal from the late inspiral of a binary neutron star (BNS) coalescence, %the possibility of assessing the presence of such a phase transition depends on the precision that can be attained in th… ▽ More Possible strong first-order hadron-quark phase transitions in neutron star interiors leave an imprint on gravitational waves, which could be detected with planned third-generation interferometers. Given a signal from the late inspiral of a binary neutron star (BNS) coalescence, %the possibility of assessing the presence of such a phase transition depends on the precision that can be attained in the determination of the tidal deformability parameter, as well as on the model used to describe the hybrid star equation of state. For the latter, we employ here a phenomenological meta-modelling of the equation of state that largely spans the parameter space associated with both the low density phase and the quark high density compatible with current constraints. We show that with a network of third-generation detectors, a single loud BNS event might be sufficient to infer the presence of a phase transition at low baryon densities with an average Bayes factor $B\approx 100$, up to a luminosity distance ($\mathcal{D}_L \lesssim$ 300 Mpc). △ Less

Submitted 11 May, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Journal ref: Mon.Not.Roy.Astron.Soc. 524 (2023) 3, 3464-3473

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

The proto-neutron star inner crust in the liquid phase

Authors: H. Dinh Thi, A. F. Fantina, F. Gulminelli

Abstract: The crust of a neutron star is known to melt at a temperature that increases with increasing matter density, up to about $10^{10}$ K. At such high temperatures and beyond, the crustal ions are put into collective motion and the associated entropy contribution can affect both the thermodynamic properties and the composition of matter. We studied the importance of this effect in different thermodyna… ▽ More The crust of a neutron star is known to melt at a temperature that increases with increasing matter density, up to about $10^{10}$ K. At such high temperatures and beyond, the crustal ions are put into collective motion and the associated entropy contribution can affect both the thermodynamic properties and the composition of matter. We studied the importance of this effect in different thermodynamic conditions relevant to the inner crust of the proto-neutron star, both at beta equilibrium and in the fixed-proton-fraction regime. To this aim, we solved the hydrodynamic equations for an ion moving in an incompressible, irrotational, and non-viscous fluid, with different boundary conditions, thus leading to different prescriptions for the ion effective mass. We then employed a compressible liquid-drop approach in the one-component plasma approximation, including the renormalisation of the ion mass to account for the influence of the surrounding medium. We show that the cluster size is determined by the competition between the ion centre-of-mass motion and the interface properties, namely the Coulomb, surface, and curvature energies. In particular, including the translational free energy in the minimisation procedure can significantly reduce the optimal number of nucleons in the clusters and lead to an early dissolution of clusters in dense beta-equilibrated matter. On the other hand, we find that the impact of translational motion is reduced in scenarios where the proton fraction is assumed constant and is almost negligible on the inner-crust equation of state. Our results show that the translational degrees of freedom affect the equilibrium composition of beta-equilibrated matter and the density and pressure of the crust-core transition in a non-negligible way, highlighting the importance of its inclusion when modelling the finite-temperature inner crust of the (proto-)neutron star. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 14 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 672, A160 (2023)

Science with the Einstein Telescope: a comparison of different designs

Authors: Marica Branchesi, Michele Maggiore, David Alonso, Charles Badger, Biswajit Banerjee, Freija Beirnaert, Enis Belgacem, Swetha Bhagwat, Guillaume Boileau, Ssohrab Borhanian, Daniel David Brown, Man Leong Chan, Giulia Cusin, Stefan L. Danilishin, Jerome Degallaix, Valerio De Luca, Arnab Dhani, Tim Dietrich, Ulyana Dupletsa, Stefano Foffa, Gabriele Franciolini, Andreas Freise, Gianluca Gemme, Boris Goncharov, Archisman Ghosh , et al. (51 additional authors not shown)

Abstract: The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where in each arm there is a `xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogeni… ▽ More The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where in each arm there is a `xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogenic temperature. Here, we examine the scientific perspectives under possible variations of this reference design. We perform a detailed evaluation of the science case for a single triangular geometry observatory, and we compare it with the results obtained for a network of two L-shaped detectors (either parallel or misaligned) located in Europe, considering different choices of arm-length for both the triangle and the 2L geometries. We also study how the science output changes in the absence of the low-frequency instrument, both for the triangle and the 2L configurations. We examine a broad class of simple `metrics' that quantify the science output, related to compact binary coalescences, multi-messenger astronomy and stochastic backgrounds, and we then examine the impact of different detector designs on a more specific set of scientific objectives. △ Less

Submitted 17 June, 2023; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 197 pages, 73 figures. v2: corrections in the part on the sensitivity to stochastic backgrounds. Accepted in JCAP

Report number: ET-0084A-23

Confronting the nucleonic hypothesis with current neutron star observations from GW170817 and PSR J0740+6620

Authors: Hoa Dinh Thi, Chiranjib Mondal, Francesca Gulminelli

Abstract: The nuclear matter equation of state is relatively well constrained at sub-saturation densities thanks to the knowledge from nuclear physics. However, studying its behavior at supra-saturation densities is a challenging task. Fortunately, the extraordinary progress recently made in observations of neutron stars and neutron star mergers has provided us with unique opportunities to unfold the proper… ▽ More The nuclear matter equation of state is relatively well constrained at sub-saturation densities thanks to the knowledge from nuclear physics. However, studying its behavior at supra-saturation densities is a challenging task. Fortunately, the extraordinary progress recently made in observations of neutron stars and neutron star mergers has provided us with unique opportunities to unfold the properties of dense matter. Under the assumption that nucleons are the only constituents of neutron star cores, we perform a Bayesian inference using the so-called meta-modeling technique with a nuclear-physics-informed prior. The latest information from the GW170817 event by the LIGO-Virgo Collaboration (LVC) and from the radius measurement of the heaviest known neutron star PSR J0740+6620 by the Neutron Star Interior Composition Explorer (NICER) telescope and X-ray Multi-Mirror (XMM-Newton) are taken into account as likelihoods in the analysis. The impacts of different constraints on the equation of state as well as on the predictions of neutron star properties are discussed. The obtained posterior reveals that all the current observations are fully compatible with the nucleonic hypothesis. Strong disagreements between our results with future data can be identified as a signal for the existence of exotic degrees of freedom. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: Contribution to the "Journees de Rencontre des Jeunes Chercheurs (JRJC) 2021" proceedings

Anisotropic electron transport in the nuclear pasta phase

Authors: Mateus R. Pelicer, Marco Antonelli, Débora P. Menezes, Francesca Gulminelli

Abstract: The presence of nuclear pasta is expected to modify the transport properties in the mantle of neutron stars. The non-spherical geometry of the pasta nuclear clusters leads to anisotropies in the collision frequencies, impacting the thermal and electrical conductivity. We derive analytical expressions for the anisotropic collision frequencies using the Boltzmann equation in the relaxation time appr… ▽ More The presence of nuclear pasta is expected to modify the transport properties in the mantle of neutron stars. The non-spherical geometry of the pasta nuclear clusters leads to anisotropies in the collision frequencies, impacting the thermal and electrical conductivity. We derive analytical expressions for the anisotropic collision frequencies using the Boltzmann equation in the relaxation time approximation. The average parallel, perpendicular and Hall electrical conductivities are computed in the high-temperature regime above crustal melting, considering incoherent elastic electron-pasta scattering and randomly oriented pasta structures. Numerical values are obtained at different densities and temperatures by using the IUFSU parametrization of the non-linear Walecka model to determine the crustal structure. We find that the anisotropy of the collision frequencies grows with the length of the pasta structures and, independently of the magnetic field, the presence of rod and slab phases decreases the conductivity by more than one order of magnitude. Our numerical results indicate that, even if the pasta structures might survive above the crustal melting point, no strong anisotropies are to be expected in the conduction properties in this temperature regime, even in the presence of a very high magnetic field. △ Less

Submitted 3 May, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 15 pages, 12 figures

Strong magnetic fields and pasta phases revisited

Authors: Luigi Scurto, Helena Pais, Francesca Gulminelli

Abstract: In this work, we compute the structure and composition of the inner crust of a neutron star in the presence of a strong magnetic field, such as it can be found in magnetars. To determine the geometry and characteristics of the crust inhomogeneities, we consider the compressible liquid drop model, where surface and Coulomb terms are included in the variational equations, and we compare our results… ▽ More In this work, we compute the structure and composition of the inner crust of a neutron star in the presence of a strong magnetic field, such as it can be found in magnetars. To determine the geometry and characteristics of the crust inhomogeneities, we consider the compressible liquid drop model, where surface and Coulomb terms are included in the variational equations, and we compare our results with previous calculations based on more approximate treatments. For the equation of state (EoS), we consider two non-linear relativistic mean-field models with different slopes of the symmetry energy, and we show that the extension of the inhomogeneous region inside the star core due to the magnetic field strongly depends on the behavior of the symmetry energy in the crustal EoS. Finally, we argue that the extended spinodal instability observed in previous calculations can be related to the presence of small amplitude density fluctuations in the magnetar outer core, rather than to a thicker solid crust. The compressible liquid drop model formalism, while in overall agreement with the previous calculations, leads to a systematic suppression of the metastable solutions, thus allowing a more precise estimation of the crust-core transition density and pressure, and therefore a better estimation of the crustal radius. △ Less

Submitted 28 April, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 11 pages, 5 figures, 3 tables. Version published on Phys. Rev. C

Journal ref: Phys. Rev. C 107, 045806 (2023)

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1680 additional authors not shown)

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

Nucleonic metamodelling in light of multimessenger, PREX-II and CREX data

Authors: Chiranjib Mondal, Francesca Gulminelli

Abstract: The need of reconciling our understanding of the behavior of hadronic matter across a wide range of densities, especially at the time when data from multimessenger observations and novel experimental facilities are flooding in, has provided new challenges to the nuclear models. Particularly, the density dependence of the isovector channel of the nuclear energy functionals seems hard to pin down if… ▽ More The need of reconciling our understanding of the behavior of hadronic matter across a wide range of densities, especially at the time when data from multimessenger observations and novel experimental facilities are flooding in, has provided new challenges to the nuclear models. Particularly, the density dependence of the isovector channel of the nuclear energy functionals seems hard to pin down if experiments like PREX-II (or PREX) and CREX are required to be taken on the same footing. We put to test this anomaly in a semi-agnostic modelling technique, by performing a full Bayesian analysis of static properties of neutron stars, together with global properties of nuclei as binding energy, charge radii and neutron skin calculated at the semi-classical level. Our results show that the interplay between bulk and surface properties, and the importance of high order empirical parameters that effectively decouple the subsaturation and the supersaturation density regime, might partially explain the tension between the different measurements and observations. If the surface behaviors, however, are decoupled from the bulk properties, we found a rather harmonious situation among experimental and observational data. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 13 pages, 12 figures

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)

Properties of pasta phases in catalyzed neutron stars

Authors: H. Dinh Thi, A. F. Fantina, F. Gulminelli

Abstract: Exotic non-spherical configurations of nuclei, known as ``pasta" phases, are expected to be present at the bottom of the inner crust of a neutron star. We study the properties of these configurations in catalyzed neutron stars within a compressible liquid-drop model approach, with surface parameters optimized to reproduce experimental nuclear masses. Our results show that the properties of the pas… ▽ More Exotic non-spherical configurations of nuclei, known as ``pasta" phases, are expected to be present at the bottom of the inner crust of a neutron star. We study the properties of these configurations in catalyzed neutron stars within a compressible liquid-drop model approach, with surface parameters optimized to reproduce experimental nuclear masses. Our results show that the properties of the pasta phases exhibit strong model dependence. To estimate the model uncertainties, a Bayesian analysis is performed, combining information from nuclear physics experiments and chiral perturbation theoretical calculations with astrophysical observations. The inferred posterior distributions are discussed, with particular focus on the effect of the low-density energy functional on the predictions. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: Proceedings of the International Workshop on Multi-facets of EOS and Clustering 2021; published in Il Nuovo Cimento C

Journal ref: IL NUOVO CIMENTO 45 C (2022) 62

A Bayesian analysis of the properties of hybrid stars with the NJL model

Authors: Antoine Pfaff, Hubert Hansen, Francesca Gulminelli

Abstract: We implement a Bayesian analysis of the properties of non-rotating hybrid stars at equilibrium with quark matter cores, as described by the SU(3) Nambu--Jona-Lasinio (NJL) model. The hadronic phase is described by a unified metamodeling approach, with a prior parameter space covering the present uncertainties on nuclear matter properties with nucleonic degrees of freedom. The parameter space of th… ▽ More We implement a Bayesian analysis of the properties of non-rotating hybrid stars at equilibrium with quark matter cores, as described by the SU(3) Nambu--Jona-Lasinio (NJL) model. The hadronic phase is described by a unified metamodeling approach, with a prior parameter space covering the present uncertainties on nuclear matter properties with nucleonic degrees of freedom. The parameter space of the NJL model includes vector-isoscalar and vector-isovector couplings and, additionally, an effective bag constant for the quark pressure is introduced as a free parameter. The phase transition is assumed to be first order with charge neutral phases, following the Maxwell construction. Our Bayesian framework includes filters on the experimental and theoretical low-density nuclear physics knowledge (atomic masses, ab initio calculations of the EoS) and high density constraints from astrophysical observations (maximum mass of J0348+0432, binary tidal deformability of the GW170817 event). We find that microscopic vector interactions play an important role in quark matter in order to stiffen the equation of state sufficiently to reach high star masses, in agreement with previous studies. Even within a very large prior for both the hadronic phase and the quark phase and the important freedom brought by the effective bag constant, our posterior quark cores tend to be relatively small and only appear in very heavy stars ($M\gtrsim 2 M_\odot$). Coincidentally, the inclusion of the nucleon-quark transition (deconfinement transition) only weakly affects the radii of compact stars, foreshadowing a very low observability of a possible phase transition using x-ray radii measurements. △ Less

Submitted 15 March, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Can we decipher the composition of the core of a neutron star?

Authors: Chiranjib Mondal, Francesca Gulminelli

Abstract: General relativity guarantees a unique one-to-one correspondence between static observables of neutron stars (NSs) accessible by multi-messenger astronomy, such as mass-radius or tidal deformability, and the equation of state (EoS) of beta equilibrated matter. It is well known that these static properties are not enough to discern conventional NSs from hybrid stars. However, if one assumes that ha… ▽ More General relativity guarantees a unique one-to-one correspondence between static observables of neutron stars (NSs) accessible by multi-messenger astronomy, such as mass-radius or tidal deformability, and the equation of state (EoS) of beta equilibrated matter. It is well known that these static properties are not enough to discern conventional NSs from hybrid stars. However, if one assumes that hadrons present in the neutron star core are only neutrons and protons, the lepton fraction is in principle determined unequivocally by the condition of chemical equilibrium. Using a simple analytical method based on a polynomial expansion of the EoS, we show that multiple solutions are possible to the beta-equilibrium equation, leading to a characteristic indetermination on the composition of the interiors of NSs, even in the purely nucleonic hypothesis. We further show that additional empirical information on symmetric matter at high densities are not very efficient to pin down the composition, if uncertainties on measurements are accounted for. We conclude that constraints on the symmetry energy at high densities only, can make meaningful impact to decipher the composition of neutron star core. Our results give a lower limit to the uncertainty on the NS core composition that can be obtained with astrophysical and terrestrial experiments. △ Less

Submitted 4 March, 2022; v1 submitted 8 November, 2021; originally announced November 2021.

Comments: 8 pages, restructured, new figures added

Journal ref: Phys. Rev. D 105, 083016 (2022)

The effect of the energy functional on the pasta-phase properties of catalysed neutron stars

Authors: H. Dinh Thi, A. F. Fantina, F. Gulminelli

Abstract: Nuclear pasta, that is an inhomogeneous distribution of nuclear matter characterised by non-spherical clustered structures, is expected to occur in a narrow spatial region at the bottom of the inner crust of neutron stars, but the width of the pasta layer is strongly model dependent. In the framework of a compressible liquid-drop model, we use Bayesian inference to analyse the constraints on the s… ▽ More Nuclear pasta, that is an inhomogeneous distribution of nuclear matter characterised by non-spherical clustered structures, is expected to occur in a narrow spatial region at the bottom of the inner crust of neutron stars, but the width of the pasta layer is strongly model dependent. In the framework of a compressible liquid-drop model, we use Bayesian inference to analyse the constraints on the sub-saturation energy functional and surface tension imposed by both ab-initio chiral perturbation theory calculations and experimental measurements of nuclear masses. The posterior models are used to obtain general predictions for the crust-pasta and pasta-core transition with controlled uncertainties. A correlation study allows extracting the most influential parameters for the calculation of the pasta phases. The important role of high-order empirical parameters and the surface tension is underlined. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 15 pages, 11 figures. Contribution to the EPJ A topical issue "CompOSE: a repository for Neutron Star Equations of State and Transport Properties"

Journal ref: Eur. Phys. J. A 57, 296 (2021)

The equation of state of neutron stars and the role of nuclear experiments

Authors: F. Gulminelli, A. F. Fantina

Abstract: Neutron stars are unique laboratories to probe matter in extreme conditions, not accessible in terrestrial laboratories. Here, we discuss the modelling of the neutron-star equation of state, particularly in connection with recent constraints coming from both nuclear physics (experiments and ab-initio calculations) and astrophysical observations. Neutron stars are unique laboratories to probe matter in extreme conditions, not accessible in terrestrial laboratories. Here, we discuss the modelling of the neutron-star equation of state, particularly in connection with recent constraints coming from both nuclear physics (experiments and ab-initio calculations) and astrophysical observations. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: 8 pages, 3 figures

Journal ref: Nuclear Physics News, Vol. 31, No 2, p. 9-13 (2021)

Uncertainties in the pasta-phase properties of catalysed neutron stars

Authors: H. Dinh Thi, T. Carreau, A. F. Fantina, F. Gulminelli

Abstract: The interior of a neutron star is expected to exhibit different states of matter. In particular, complex non-spherical configurations known as `pasta' phases may exist at the highest densities in the inner crust, potentially having an impact on different neutron-star phenomena. We study the properties of the pasta phase and the uncertainties in the pasta observables which are due to our incomplete… ▽ More The interior of a neutron star is expected to exhibit different states of matter. In particular, complex non-spherical configurations known as `pasta' phases may exist at the highest densities in the inner crust, potentially having an impact on different neutron-star phenomena. We study the properties of the pasta phase and the uncertainties in the pasta observables which are due to our incomplete knowledge of the nuclear energy functional. To this aim, we employed a compressible liquid-drop model approach with surface parameters optimised either on experimental nuclear masses or theoretical calculations. To assess the model uncertainties, we performed a Bayesian analysis by largely varying the model parameters using uniform priors, and generating posterior distributions with filters accounting for both our present low-density nuclear physics knowledge and high-density neutron-star physics constraints. Our results show that the nuclear physics constraints, such as the neutron-matter equation of state at very low density and the experimental mass measurements, are crucial in determining the crustal and pasta observables. Accounting for all constraints, we demonstrate that the presence of pasta phases is robustly predicted in an important fraction of the inner crust. We estimate the relative crustal thickness associated with pasta phases as $R_{\rm pasta}/R_{\rm crust}=0.128\pm 0.047$ and the relative moment of inertia as $I_{\rm pasta}/I_{\rm crust}=0.480\pm 0.137$. Our findings indicate that the surface and curvature parameters are more influential than the bulk parameters for the description of pasta observables. We also show that using a surface tension that is inconsistent with the bulk functional leads to an underestimation of both the average values and the uncertainties in the pasta properties, thus highlighting the importance of a consistent calculation of the nuclear functional. △ Less

Submitted 27 October, 2021; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: 12 pages, 8 figures, accepted for publication in Astronomy and Astrophysics. Typos corrected

Journal ref: A&A 654, A114 (2021)

The nuclear matter density functional under the nucleonic hypothesis

Authors: Hoa Dinh Thi, Chiranjib Mondal, Francesca Gulminelli

Abstract: A Bayesian analysis of the possible behaviors of the dense matter equation of state informed through recent LIGO-Virgo as well as NICER measurements reveals that all the present observations are compatible with a fully nucleonic hypothesis for the composition of dense matter, even in the core of the most massive pulsar PSR J0740+6620. Under the hypothesis of a nucleonic composition, we extract the… ▽ More A Bayesian analysis of the possible behaviors of the dense matter equation of state informed through recent LIGO-Virgo as well as NICER measurements reveals that all the present observations are compatible with a fully nucleonic hypothesis for the composition of dense matter, even in the core of the most massive pulsar PSR J0740+6620. Under the hypothesis of a nucleonic composition, we extract the most general behavior of the energy per particle of symmetric matter and density dependence of the symmetry energy, compatible with the astrophysical observations as well as our present knowledge of low energy nuclear physics from effective field theory predictions and experimental nuclear mass data. These results can be used as a null hypothesis to be confronted with future constraints on dense matter to search for possible exotic degrees of freedom. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: Submitted to Universe, invited contribution to the special Issue: Properties and Dynamics of Neutron Stars and Proto-Neutron Stars, V.Dexheimer and R.Negreiros eds., 2021

Journal ref: Universe 2021, 7(10), 373

Binding energy shifts from heavy-ion experiments in a nuclear statistical equilibrium model

Authors: S. Mallik, H. Pais, F. Gulminelli

Abstract: Chemical constants extracted from $^{124}$Xe+ $^{124}$Sn collisions at 32 AMeV are compared to the predictions of an extended Nuclear Statistical Equilibrium model including mean-field interactions and in-medium binding energy shifts for the light ($Z\leq 2$) clusters. The ion species and density dependence of the in-medium modification is directly extracted from the experimental data. We show tha… ▽ More Chemical constants extracted from $^{124}$Xe+ $^{124}$Sn collisions at 32 AMeV are compared to the predictions of an extended Nuclear Statistical Equilibrium model including mean-field interactions and in-medium binding energy shifts for the light ($Z\leq 2$) clusters. The ion species and density dependence of the in-medium modification is directly extracted from the experimental data. We show that the shift increases with the mass of the cluster and the density of the medium, and we provide a simple linear fit for future use in astrophysical simulations in the framework of the CompOSE data base. The resulting mass fractions are computed in representative thermodynamic conditions relevant for supernova and neutron star mergers. A comparison to the results of a similar analysis of the same data performed in the framework of a relativistic mean-field model shows a good agreement at low density, but significant discrepancies close to the Mott dissolution of clusters in the dense medium. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 7 pages, 4 figures

Journal ref: Eur. Phys. J. A (2021) 57:262

Fluctuations in the pasta phase

Authors: Mateus R. Pelicer, Débora P. Menezes, Celso C. Barros Jr, Francesca Gulminelli

Abstract: Baryonic matter close to the saturation density is very likely to present complex inhomogeneous structures collectively known under the name of pasta phase. At finite temperature, the different geometric structures are expected to coexist, with potential consequences on the neutron star crust conductivity and neutrino transport in supernova matter. In the framework of a statistical multi-component… ▽ More Baryonic matter close to the saturation density is very likely to present complex inhomogeneous structures collectively known under the name of pasta phase. At finite temperature, the different geometric structures are expected to coexist, with potential consequences on the neutron star crust conductivity and neutrino transport in supernova matter. In the framework of a statistical multi-component approach, we calculate the composition of matter in the pasta phase considering density, proton fraction, and geometry fluctuations. Using a realistic energy functional from relativistic mean field theory and a temperature and isospin dependent surface tension fitted from Thomas-Fermi calculations, we show that different geometries can coexist in a large fraction of the pasta phase, down to temperatures of the order of the crystallization temperature of the neutron star crust. Quantitative estimates of the charge fluctuations are given. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 6 pages, 3 figures, 1 table

Journal ref: Phys. Rev. C 104, 022801 (2021)

Parametrization of the surface energy in the ETF approximation

Authors: Ubiratãn José Furtado, Francesca Gulminelli

Abstract: We perform extended Thomas-Fermi calculations in the Wigner-Seitz cell below and above neutron drip with realistic functionals. The resulting energy density is decomposed as a sum of bulk terms and a surface term, and a compressible liquid drop analytical formula is used to fit the surface tension. The effect of curvature terms and neutron skin is studied in detail. A very good reproduction of the… ▽ More We perform extended Thomas-Fermi calculations in the Wigner-Seitz cell below and above neutron drip with realistic functionals. The resulting energy density is decomposed as a sum of bulk terms and a surface term, and a compressible liquid drop analytical formula is used to fit the surface tension. The effect of curvature terms and neutron skin is studied in detail. A very good reproduction of the microscopic data is obtained using an expression that depends only on the mass and charge of the cluster, showing that the in-medium modifications of the nuclear energy in the presence of an external neutron gas can be effectively accounted for in the isospin dependence of the surface tension. In this first application, aimed at establishing the fitting protocol, we concentrate on the Sly4 energy functional, but the study can be easily generalized to different functionals, and the resulting parametrizations can be used for direct applications in pasta calculations in neutron star crusts and supernova matter. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 44 pages, 17 figures, 5 tables

Light clusters in warm stellar matter: calibrating the cluster couplings

Authors: Tiago Custódio, Alexandre Falcão, Helena Pais, Constança Providência, Francesca Gulminelli, Gerd Röpke

Abstract: The abundances of light clusters within a formalism that considers in-medium effects are calculated using several relativistic mean-field models, with both density-dependent and density-independent couplings. Clusters are introduced as new quasiparticles, with a modified coupling to the scalar meson field. A comparison with experimental data from heavy ion collisions allows settling the model depe… ▽ More The abundances of light clusters within a formalism that considers in-medium effects are calculated using several relativistic mean-field models, with both density-dependent and density-independent couplings. Clusters are introduced as new quasiparticles, with a modified coupling to the scalar meson field. A comparison with experimental data from heavy ion collisions allows settling the model dependence of the results and the determination of the couplings of the light clusters to the meson fields. We find that extra experimental constraints at higher density are needed to convincingly pin down the density associated to the melting of clusters in the dense nuclear medium. The role of neutron rich clusters, such as $^6$He, in asymmetric matter is discussed. △ Less

Submitted 25 November, 2020; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: 10 pages, 8 figures, and 1 table; accepted for publication in the EPJA Topical Issue "Light Clusters in Nuclei and Nuclear Matter: Nuclear Structure and Decay, Heavy Ion Collisions, and Astrophysics"

Journal ref: Eur. Phys. J. A 56, 295 (2020)

Improved method for the experimental determination of in-medium effects from heavy-ion collisions

Authors: Helena Pais, Rémi Bougault, Francesca Gulminelli, Constança Providência, Eric Bonnet, Bernard Borderie, Abdelouahad Chbihi, John D. Frankland, Emmanuelle Galichet, Diégo Gruyer, Maxime Henri, Nicolas Le Neindre, Olivier Lopez, Loredana Manduci, Marian Pârlog, Giuseppe Verde

Abstract: The equation of state with light clusters for nuclear and stellar matter is determined using chemical equilibrium constants evaluated from the analysis of the recently published (Xe$+$Sn) heavy ion data, corresponding to three reactions with different isotopic contents of the emission source. The measured multiplicities are used to extract the thermodynamic properties, and an in-medium correction… ▽ More The equation of state with light clusters for nuclear and stellar matter is determined using chemical equilibrium constants evaluated from the analysis of the recently published (Xe$+$Sn) heavy ion data, corresponding to three reactions with different isotopic contents of the emission source. The measured multiplicities are used to extract the thermodynamic properties, and an in-medium correction to the ideal gas internal partition function of the clusters is included in the analysis. This in-medium correction and its respective uncertainty are calculated via a Bayesian analysis, with the unique hypothesis that the different nuclear species in a given sample must correspond to a unique common value for the density of the expanding source. Different parameter sets for the correction are tested, and the effect of the radius of the clusters on the thermodynamics and on the chemical equilibrium constants is also addressed. It is shown that the equilibrium constants obtained are almost independent of the isospin content of the analysed systems. Finally, a comparison with a relativistic mean field model proves that data are consistent with a universal in-medium correction of the scalar $σ$-meson coupling for nucleons bound in clusters. The obtained value, $g_s/g_s^0 = 0.92 \pm 0.02$, is larger than that obtained in a previous study not including in-medium effects in the data analysis. This result implies a smaller effect on the binding energy of the clusters and, as a consequence, larger melting densities, and an increased cluster contribution in supernova matter. △ Less

Submitted 8 September, 2020; v1 submitted 12 June, 2020; originally announced June 2020.

Comments: 29 pages, 16 figures, accepted for publication in J. Phys. G

Journal ref: J. Phys. G: Nucl. Part. Phys. 47, 105204 (2020)

Inner crust of a neutron star at crystallization in a multi-component approach

Authors: T. Carreau, A. F. Fantina, F. Gulminelli

Abstract: The possible presence of amorphous and heterogeneous phases in the inner crust of a neutron star is expected to reduce the electrical conductivity of the crust, with potentially important consequences on the magneto-thermal evolution of the star. In cooling simulations, the disorder is quantified by an impurity parameter which is often taken as a free parameter. We aim to give a quantitative predi… ▽ More The possible presence of amorphous and heterogeneous phases in the inner crust of a neutron star is expected to reduce the electrical conductivity of the crust, with potentially important consequences on the magneto-thermal evolution of the star. In cooling simulations, the disorder is quantified by an impurity parameter which is often taken as a free parameter. We aim to give a quantitative prediction of the impurity parameter as a function of the density in the crust,performing microscopic calculations including up-to-date microphysics of the crust. A multi-component approach is developed at finite temperature using a compressible liquid drop description of the ions with an improved energy functional based on recent microscopic nuclear models and optimized on extended Thomas-Fermi calculations. Thermodynamic consistency is ensured by adding a rearrangement term and deviations from the linear mixing rule are included in the liquid phase. The impurity parameter is consistently calculated at the crystallization temperature as determined in the one-component plasma approximation for the different functionals. Our calculations show that at the crystallization temperature the composition of the inner crust is dominated by nuclei with charge number around $Z \approx 40$, while the range of the $Z$ distribution varies from about 20 near the neutron drip to about 40 closer to the crust-core transition. This reflects on the behavior of the impurity parameter that monotonically increases with density up to around 40 in the deeper regions of the inner crust. Our study shows that the contribution of impurities is non-negligible, thus potentially having an impact on the transport properties in the neutron-star crust. The obtained values of the impurity parameter represent a lower limit; larger values are expected in the presence of non-spherical geometries and/or fast cooling dynamics. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Comments: 10 pages, 6 figures

Crystallization of the outer crust of a non-accreting neutron star

Authors: A. F. Fantina, S. De Ridder, N. Chamel, F. Gulminelli

Abstract: The interior of a neutron star (NS) is usually assumed to be made of cold catalyzed matter. However, the outer layers are unlikely to remain in full equilibrium during the formation of the star and its cooling, especially after crystallization. We study the cooling and equilibrium composition of the outer layers of a NS down to crystallization. Here the impurity parameter, usually a free parameter… ▽ More The interior of a neutron star (NS) is usually assumed to be made of cold catalyzed matter. However, the outer layers are unlikely to remain in full equilibrium during the formation of the star and its cooling, especially after crystallization. We study the cooling and equilibrium composition of the outer layers of a NS down to crystallization. Here the impurity parameter, usually a free parameter in cooling simulations, is calculated self-consistently using a microscopic model for which a unified equation of state has recently been determined. We follow the evolution of the nuclear distributions of the multi-component Coulomb liquid plasma (MCP) fully self-consistently, adapting a general formalism originally developed for the description of supernova cores. We calculate the impurity parameter at the crystallization as determined in the one-component plasma (OCP) approximation. Our analysis shows that the sharp changes in composition obtained in the OCP approximation are smoothed out when a full nuclear distribution is allowed. The Coulomb coupling parameter at melting is found to be reasonably close to the canonical value of 175, except for specific pressures for which supercooling occurs in the OCP approximation. Our MCP treatment leads to non-monotonic variations of the impurity parameter with pressure. Its values can change by several orders of magnitude reaching about 50, suggesting that the crust may be composed of an alternation of pure (highly conductive) and impure (highly resistive) layers. The results presented here complement the recent unified equation of state obtained with the same model. Our self-consistent approach to hot MCP shows that the presence of impurities in the outer crust of a NS is non-negligible and may have a sizeable impact on transport properties. In turn, this may have important implications for the cooling of NS and their magneto-rotational evolution. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: 15 pages, 8 figures, accepted for publication in Astronomy and Astrophysics

Crystallization of the inner crust of a neutron star and the influence of shell effects

Authors: T. Carreau, F. Gulminelli, N. Chamel, A. F. Fantina, J. M. Pearson

Abstract: Context. In the cooling process of a non-accreting neutron star, the composition and properties of the crust are thought to be fixed at the finite temperature where nuclear reactions fall out of equilibrium. A lower estimation for this temperature is given by the crystallization temperature, which can be as high as $\approx 7\times 10^9$ K in the inner crust, potentially leading to sizeable differ… ▽ More Context. In the cooling process of a non-accreting neutron star, the composition and properties of the crust are thought to be fixed at the finite temperature where nuclear reactions fall out of equilibrium. A lower estimation for this temperature is given by the crystallization temperature, which can be as high as $\approx 7\times 10^9$ K in the inner crust, potentially leading to sizeable differences with respect to the simplifying cold-catalyzed matter hypothesis. Aims. We extend the recent work by Fantina et al. (2019) on the outer crust, to the study of the crystallization of the inner crust and the associated composition in the one-component plasma approximation. Methods. The finite temperature variational equations for non-uniform matter in both the liquid and the solid phases are solved using a compressible liquid-drop approach with parameters optimized on four different microscopic models which cover the present uncertainties in nuclear modeling. Results. We consider separately the effect of the different nuclear ingredients with their associated uncertainties, namely the nuclear equation of state, the surface properties in the presence of a uniform gas of dripped neutrons, and the proton shell effects arising from the ion single-particle structure. Our results suggest that the highest source of model dependence comes from the smooth part of the nuclear functional. Conclusions. We show that shell effects play an important role at the lowest densities close to the outer crust, but the most important physical ingredient to be settled for a quantitative prediction of the inner crust properties is the surface tension at extreme isospin values. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Comments: 9 pages, 9 figures

Journal ref: A&A 635, A84 (2020)

Low-density in-medium effects on light clusters from heavy-ion data

Authors: Helena Pais, Rémi Bougault, Francesca Gulminelli, Constança Providência, Eric Bonnet, Bernard Borderie, Abdelouahad Chbihi, John D. Frankland, Emmanuelle Galichet, Diégo Gruyer, Maxime Henri, Nicolas Le Neindre, Olivier Lopez, Loredana Manduci, Marian Parlôg, Giuseppe Verde

Abstract: The modification of the ground state properties of light atomic nuclei in the nuclear and stellar medium is addressed, using chemical equilibrium constants evaluated from a new analysis of the intermediate energy heavy-ion (Xe$+$Sn) collision data measured by the INDRA collaboration. Three different reactions are considered, mainly differing by the isotopic content of the emission source. The ther… ▽ More The modification of the ground state properties of light atomic nuclei in the nuclear and stellar medium is addressed, using chemical equilibrium constants evaluated from a new analysis of the intermediate energy heavy-ion (Xe$+$Sn) collision data measured by the INDRA collaboration. Three different reactions are considered, mainly differing by the isotopic content of the emission source. The thermodynamic conditions of the data samples are extracted from the measured multiplicities allowing for a parametrization of the in-medium modification, determined with the single hypothesis that the different nuclear species in a given sample correspond to a unique common value for the density of the expanding source. We show that this correction, which was not considered in previous analyses of chemical constants from heavy ion collisions, is necessary, since the observables of the analyzed systems show strong deviations from the expected results for an ideal gas of free clusters. This data set is further compared to a relativistic mean-field model, and seen to be reasonably compatible with a universal correction of the attractive $σ$-meson coupling. △ Less

Submitted 12 June, 2020; v1 submitted 25 November, 2019; originally announced November 2019.

Comments: 5 pages, 2 figures, 1 table, accepted for publication in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 125, 012701 (2020)

Impact of electron capture rates on nuclei far from stability on core-collapse supernovae

Authors: Aurélien Pascal, Simon Giraud, Anthea Fantina, Francesca Gulminelli, Jerome Novak, Micaela Oertel, Adriana Raduta

Abstract: The impact of electron-capture (EC) cross sections on neutron-rich nuclei on the dynamics of core-collapse during infall and early post-bounce is studied performing spherically symmetric simulations in general relativity using a multigroup scheme for neutrino transport and full nuclear distributions in extended nuclear statistical equilibrium models. We thereby vary the prescription for EC rates o… ▽ More The impact of electron-capture (EC) cross sections on neutron-rich nuclei on the dynamics of core-collapse during infall and early post-bounce is studied performing spherically symmetric simulations in general relativity using a multigroup scheme for neutrino transport and full nuclear distributions in extended nuclear statistical equilibrium models. We thereby vary the prescription for EC rates on individual nuclei, the nuclear interaction for the EoS, the mass model for the nuclear statistical equilibrium distribution and the progenitor model. In agreement with previous works, we show that the individual EC rates are the most important source of uncertainty in the simulations, while the other inputs only marginally influence the results. A recently proposed analytic formula to extrapolate microscopic results on stable nuclei for EC rates to the neutron rich region, with a functional form motivated by nuclear-structure data and parameters fitted from large scale shell model calculations, is shown to lead to a sizable (16%) reduction of the electron fraction at bounce compared to more primitive prescriptions for the rates, leading to smaller inner core masses and slower shock propagation. We show that the EC process involves $\approx$ 130 different nuclear species around 86 Kr mainly in the N = 50 shell closure region, and establish a list of the most important nuclei to be studied in order to constrain the global rates. △ Less

Submitted 2 September, 2019; v1 submitted 12 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. C 101, 015803 (2020)

Full distribution of clusters with universal couplings and in-medium effects

Authors: Helena Pais, Francesca Gulminelli, Constança Providência, Gerd Röpke

Abstract: Light and heavy clusters are calculated for asymmetric warm nuclear matter in a relativistic mean-field approach. In-medium effects, introduced via a universal cluster-meson coupling, and a binding energy shift contribution, calculated in a Thomas-Fermi approximation, were taken into account. This work considers, besides the standard lightest bound clusters $^4$He, $^3$He, $^3$H, and $^2$H, also s… ▽ More Light and heavy clusters are calculated for asymmetric warm nuclear matter in a relativistic mean-field approach. In-medium effects, introduced via a universal cluster-meson coupling, and a binding energy shift contribution, calculated in a Thomas-Fermi approximation, were taken into account. This work considers, besides the standard lightest bound clusters $^4$He, $^3$He, $^3$H, and $^2$H, also stable and unstable clusters with higher number of nucleons, in the range $5\leq A \leq 12$, as it is natural that heavier clusters also form in core-collapse supernova matter, before the pasta phases set in. We show that these extra degrees of freedom contribute with non-negligible mass fractions to the composition of nuclear matter, and may prevail over deuterons and $α$ particles at high density in strongly asymmetric matter, and not too high temperatures. The presence of the light clusters reduces the contribution of heavy clusters to a much smaller density range, and to a smaller mass fraction. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: 17 pages, 14 figures, accepted for publication in PRC

Journal ref: Phys. Rev. C 99, 055806 (2019)

Quantifying the uncertainties on spinodal instability in stellar matter through meta-modeling

Authors: Sofija Antic, Debarati Chatterjee, Thomas Carreau, Francesca Gulminelli

Abstract: The influence of the uncertainties of the equation of state empirical parameters on the neutron stars crust-core phase transition is explored within a meta-modeling approach, in which the energy per particle is expanded as a Taylor series in density and asymmetry around the saturation point. The phase transition point is estimated from the intersection of the spinodal instability region for dynami… ▽ More The influence of the uncertainties of the equation of state empirical parameters on the neutron stars crust-core phase transition is explored within a meta-modeling approach, in which the energy per particle is expanded as a Taylor series in density and asymmetry around the saturation point. The phase transition point is estimated from the intersection of the spinodal instability region for dynamical fluctuations with the chemical equilibrium curve. Special attention is paid to the inclusion of high-order parameters of the Taylor series and their influence on the transition point. An uncorrelated prior distribution is considered for the empirical parameters, with bulk properties constrained through effective field theory predictions, while the surface parameters are controlled from a fit of nuclear masses using the extended Thomas Fermi approximation. The results show that the isovector compressibility $K_{sym}$ and skewness $Q_{sym}$ have the most significant correlations with the transition point, along with the previously observed influence of the $L_{sym}$ parameter. The estimated density and pressure of the crust-core transition are $n_t = (0.071 \pm 0.011) fm^{-3}$ and $P_t = (0.294 \pm 0.102) MeV fm^{-3}$. △ Less

Submitted 13 January, 2019; originally announced January 2019.

Comments: 11 pages, 7 figures

Estimating magnetar radii with an empirical meta-model

Authors: Debarati Chatterjee, Francesca Gulminelli, Debora P. Menezes

Abstract: The presence of strong magnetic fields in neutron stars, such as in magnetars, may significantly affect their crust-core transition properties and the crust size. This knowledge is crucial in the correct interpretation of astrophysical phenomena involving magnetars, such as glitches in observed rotation frequencies, cooling, bursts and possibly tidal polarizabilities. A recently developed meta-mod… ▽ More The presence of strong magnetic fields in neutron stars, such as in magnetars, may significantly affect their crust-core transition properties and the crust size. This knowledge is crucial in the correct interpretation of astrophysical phenomena involving magnetars, such as glitches in observed rotation frequencies, cooling, bursts and possibly tidal polarizabilities. A recently developed meta-modelling technique allows exploring the model dependence of density functional theory equation of state calculations. In this work, we extend this meta-model to investigate the effect of strong magnetic fields on spinodal instabilities of neutron star matter and the associated crust-core properties. Both Tolman-Oppenheimer-Volkov and a full self-consistent numerical calculations are performed for the neutron star structure, the results being quantitatively different for strong magnetic fields. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 24 pages, 10 figures, 9 tables, regular journal article

General predictions for the neutron star crustal moment of inertia

Authors: Thomas Carreau, Francesca Gulminelli, Jérôme Margueron

Abstract: The neutron star crustal EoS and transition point properties are computed within a unified meta-modeling approach. A Bayesian approach is employed including two types of filters: bulk nuclear properties are controlled from low density effective field theory (EFT) predictions as well as the present knowledge from nuclear experiments, while the surface energy is adjusted on experimental nuclear mass… ▽ More The neutron star crustal EoS and transition point properties are computed within a unified meta-modeling approach. A Bayesian approach is employed including two types of filters: bulk nuclear properties are controlled from low density effective field theory (EFT) predictions as well as the present knowledge from nuclear experiments, while the surface energy is adjusted on experimental nuclear masses. Considering these constraints, a quantitative prediction of crustal properties can be reached with controlled confidence intervals and increased precision with respect to previous calculations: {$\approx 11\%$} dispersion on the crustal width and {$\approx 27\%$} dispersion on the fractional moment of inertia. The crust moment of inertia is also evaluated as a function of the neutron star mass, and predictions for mass and radii are given for different pulsars. The possible crustal origin of Vela pulsar glitches is discussed within the present estimations of crustal entrainment, disfavoring a large entrainment phenomenon if the Vela mass is above $1.4M_\odot$. Further refinement of the present predictions requires a better estimation of the high order isovector empirical parameters, e.g. $K_{sym}$ and $Q_{sym}$, and a better control of the surface properties of extremely neutron rich nuclei. △ Less

Submitted 15 February, 2019; v1 submitted 1 October, 2018; originally announced October 2018.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. C 100, 055803 (2019)

On the effect of high order empirical parameters on the nuclear equation of state

Authors: J. Margueron, F. Gulminelli

Abstract: A quantitative knowledge of the nuclear Equation of State (EoS) requires an accurate estimation of the uncertainties on the EoS parameters and their mutual correlations. Such correlations are empirically observed in a large set of EoS models by different authors, but they are not always fully understood. We show that some of these correlations can be interpreted from basic physical constraints imp… ▽ More A quantitative knowledge of the nuclear Equation of State (EoS) requires an accurate estimation of the uncertainties on the EoS parameters and their mutual correlations. Such correlations are empirically observed in a large set of EoS models by different authors, but they are not always fully understood. We show that some of these correlations can be interpreted from basic physical constraints imposed on a simple Taylor expansion of the binding energy around saturation density. In particular, we investigate the correlations among the following empirical parameters: the symmetry energy $E_{sym}$, the slope and curvature of the symmetry energy $L_{sym}$ and $K_{sym}$, and the curvature and skewness of the binding energy in symmetric matter $K_{sat}$ and $Q_{sat}$. The uncertainties on these correlations is estimated through analytical modelling as well as a meta-modelling analysis of the EoS subject to physical constraints. We show that a huge dispersion of the correlations among low order empirical parameters is induced by the unknown higher order empirical parameters, such as $K_{sym}$ (second order) and $Q_{sat}$ and $Q_{sym}$ (third order). We also propose an explanation of the reason why $Q_{sat}$ is weakly constrained by present experimental data. We conclude that selective observables on high order parameters, such as $K_{sym}$ and $Q_{sat}$, should be better determined before the present uncertainties on the EoS can be further reduced. △ Less

Submitted 4 July, 2018; originally announced July 2018.

Report number: INT-PUB-18-028

Journal ref: Phys. Rev. C 99, 025806 (2019)

French SKA White Book - The French Community towards the Square Kilometre Array

Authors: F. Acero, J. -T. Acquaviva, R. Adam, N. Aghanim, M. Allen, M. Alves, R. Ammanouil, R. Ansari, A. Araudo, E. Armengaud, B. Ascaso, E. Athanassoula, D. Aubert, S. Babak, A. Bacmann, A. Banday, K. Barriere, F. Bellossi, J. -P. Bernard, M. G. Bernardini, M. Béthermin, E. Blanc, L. Blanchet, J. Bobin, S. Boissier , et al. (153 additional authors not shown)

Abstract: The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive p… ▽ More The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration. △ Less

Submitted 28 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferrière, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Deschênes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

Examination of strangeness instabilities and effects of strange meson couplings in dense strange hadronic matter and compact stars

Authors: James R. Torres, Francesca Gulminelli, Debora P. Menezes

Abstract: Background : The emergence of hyperon degrees of freedom in neutron star matter has been associated to first order phase transitions in some phenomenological models, but conclusions on the possible physical existence of an instability in the strangeness sector are strongly model dependent. Purpose : The purpose of the present study is to assess whether strangeness instabilities are related to spec… ▽ More Background : The emergence of hyperon degrees of freedom in neutron star matter has been associated to first order phase transitions in some phenomenological models, but conclusions on the possible physical existence of an instability in the strangeness sector are strongly model dependent. Purpose : The purpose of the present study is to assess whether strangeness instabilities are related to specific values of the largely unconstrained hyperon interactions, and to study the effect of the strange meson couplings on phenomenological properties of neutron stars and supernova matter, once these latter are fixed to fulfill the constraints imposed by hypernuclear data. Method : We consider a phenomenological RMF model sufficiently simple to allow a complete exploration of the parameter space. Results : We show that no instability at supersaturation density exists for the RMF model, as long as the parameter space is constrained by basic physical requirements. This is at variance with a non-relativistic functional, with a functional behavior fitted through ab-initio calculations. Once the study is extended to include the full octet, we show that the parameter space allows reasonable radii for canonical neutron stars as well as massive stars above two-solar mass, together with an important strangeness content of the order of 30\%, slightly decreasing with increasing entropy, even in the absence of a strangeness driven phase transition. Conclusions : We conclude that the hyperon content of neutron stars and supernova matter cannot be established with present constraints, and is essentially governed by the unconstrained coupling to the strange isoscalar meson. △ Less

Submitted 8 February, 2017; v1 submitted 17 August, 2016; originally announced August 2016.

Comments: 13 pages, 8 figures

Journal ref: Phys. Rev. C 95, 025201 (2017)

Neutron star radii and crusts: uncertainties and unified equations of state

Authors: M. Fortin, C. Providencia, A. R. Raduta, F. Gulminelli, J. L Zdunik, P. Haensel, M. Bejger

Abstract: The uncertainties in neutron star (NS) radii and crust properties due to our limited knowledge of the equation of state (EOS) are quantitatively analysed. We first demonstrate the importance of a unified microscopic description for the different baryonic densities of the star. If the pressure functional is obtained matching a crust and a core EOS based on models with different properties at nuclea… ▽ More The uncertainties in neutron star (NS) radii and crust properties due to our limited knowledge of the equation of state (EOS) are quantitatively analysed. We first demonstrate the importance of a unified microscopic description for the different baryonic densities of the star. If the pressure functional is obtained matching a crust and a core EOS based on models with different properties at nuclear matter saturation, the uncertainties can be as large as $\sim 30\%$ for the crust thickness and $4\%$ for the radius. Necessary conditions for causal and thermodynamically consistent matchings between the core and the crust are formulated and their consequences examined. A large set of unified EOS for purely nucleonic matter is obtained based on 24 Skyrme interactions and 9 relativistic mean-field nuclear parametrizations. In addition, for relativistic models 17 EOS including a transition to hyperonic matter at high density are presented. All these EOS have in common the property of describing a $2\;M_\odot$ star and of being causal within stable NS. A span of $\sim 3$ km and $\sim 4$ km is obtained for the radius of, respectively, $1.0\;M_\odot$ and $2.0\;M_\odot$ star. Applying a set of nine further constraints from experiment and ab-initio calculations the uncertainty is reduced to $\sim 1$ km and $2$ km, respectively. These residual uncertainties reflect lack of constraints at large densities and insufficient information on the density dependence of the EOS near the nuclear matter saturation point. The most important parameter to be constrained is shown to be the symmetry energy slope $L$ which exhibits a linear correlation with the stellar radius, particularly for masses $\sim 1.0\;M_\odot$. Potential constraints on $L$, the NS radius and the EOS from observations of thermal states of NS are also discussed. [Abriged] △ Less

Submitted 7 April, 2016; originally announced April 2016.

Comments: Submitted to Phys. Rev. C. Supplemental material not included

Clusterized nuclear matter in the (proto-)neutron star crust and the symmetry energy

Authors: Ad. R. Raduta, F. Aymard, F. Gulminelli

Abstract: Though generally agreed that the symmetry energy plays a dramatic role in determining the structure of neutron stars and the evolution of core-collapsing supernovae, little is known in what concerns its value away from normal nuclear matter density and, even more important, the correct definition of this quantity in the case of unhomogeneous matter. Indeed, nuclear matter traditionally addressed b… ▽ More Though generally agreed that the symmetry energy plays a dramatic role in determining the structure of neutron stars and the evolution of core-collapsing supernovae, little is known in what concerns its value away from normal nuclear matter density and, even more important, the correct definition of this quantity in the case of unhomogeneous matter. Indeed, nuclear matter traditionally addressed by mean-field models is uniform while clusters are known to exist in the dilute baryonic matter which constitutes the main component of compact objects outer shells. In the present work we investigate the meaning of symmetry energy in the case of clusterized systems and the sensitivity of the proto-neutron star composition and equation of state to the effective interaction. To this aim an improved Nuclear Statistical Equilibrium (NSE) model is developed, where the same effective interaction is consistently used to determine the clusters and unbound particles energy functionals in the self-consistent mean-field approximation. In the same framework, in-medium modifications to the cluster energies due to the presence of the nuclear gas are evaluated. We show that the excluded volume effect does not exhaust the in-medium effects and an extra isospin and density dependent energy shift has to be considered to consistently determine the composition of subsaturation stellar matter. The symmetry energy of diluted matter is seen to depend on the isovector properties of the effective interaction, but its behavior with density and its quantitative value are strongly modified by clusterization. △ Less

Submitted 16 July, 2013; originally announced July 2013.

Comments: A contribution to the upcoming EPJA Special Volume on Nuclear Symmetry Energy

Journal ref: Eur. Phys. J. A (2014) 50:24

Ensemble in-equivalence in supernova matter within a simple model

Authors: F. Gulminelli, Ad. R. Raduta

Abstract: A simple, exactly solvable statistical model is presented for the description of baryonic matter in the thermodynamic conditions associated to the evolution of core-collapsing supernova. It is shown that the model presents a first order phase transition in the grandcanonical ensemble which is not observed in the canonical ensemble. Similar to other model systems studied in condensed matter physics… ▽ More A simple, exactly solvable statistical model is presented for the description of baryonic matter in the thermodynamic conditions associated to the evolution of core-collapsing supernova. It is shown that the model presents a first order phase transition in the grandcanonical ensemble which is not observed in the canonical ensemble. Similar to other model systems studied in condensed matter physics, this ensemble in-equivalence is accompanied by negative susceptibility and discontinuities in the intensive observables conjugated to the order parameter. This peculiar behavior originates from the fact that baryonic matter is subject to attractive short range strong forces as well as repulsive long range electromagnetic interactions, partially screened by a background of electrons. As such, it is expected in any theoretical treatment of nuclear matter in the stellar environment. Consequences for the phenomenology of supernova dynamics are drawn. △ Less

Submitted 10 October, 2011; originally announced October 2011.

Journal ref: Phys. Rev. C 85, 025803 (2012)

Boundary conditions for star matter and other periodic fermionic systems

Authors: Francesca Gulminelli, Takuya Furuta, Olivier Juillet, Christian Leclercq

Abstract: Bulk fermionic matter, as it can be notably found in supernova matter and neutrons stars, is subject to correlations of infinite range due to the antisymmetrisation of the N-body wave function, which cannot be explicitly accounted for in a practical simulation. This problem is usually addressed in condensed matter physics by means of the so-called Twist Averaged Boundary Condition method. A differ… ▽ More Bulk fermionic matter, as it can be notably found in supernova matter and neutrons stars, is subject to correlations of infinite range due to the antisymmetrisation of the N-body wave function, which cannot be explicitly accounted for in a practical simulation. This problem is usually addressed in condensed matter physics by means of the so-called Twist Averaged Boundary Condition method. A different ansatz based on the localized Wannier representation has been proposed in the context of antisymmetrized molecular dynamics. In this paper we work out the formal relation between the two approaches. We show that, while the two coincide when working with exact eigenstates of the N-body Hamiltonian, differences appear in the case of variational approaches, which are currently used for the description of stellar matter. Some model applications with Fermionic Molecular Dynamics are shown. △ Less

Submitted 3 October, 2011; originally announced October 2011.

Statistical description of complex nuclear phases in supernovae and proto-neutron stars

Authors: Ad. R. Raduta, F. Gulminelli

Abstract: We develop a phenomenological statistical model for dilute star matter at finite temperature, in which free nucleons are treated within a mean-field approximation and nuclei are considered to form a loosely interacting cluster gas. Its domain of applicability, that is baryonic densities ranging from about $ρ>10^8$ g $\cdot$ cm$^{-3}$ to normal nuclear density, temperatures between 1 and 20 MeV and… ▽ More We develop a phenomenological statistical model for dilute star matter at finite temperature, in which free nucleons are treated within a mean-field approximation and nuclei are considered to form a loosely interacting cluster gas. Its domain of applicability, that is baryonic densities ranging from about $ρ>10^8$ g $\cdot$ cm$^{-3}$ to normal nuclear density, temperatures between 1 and 20 MeV and proton fractions between 0.5 and 0, make it suitable for the description of baryonic matter produced in supernovae explosions and proto-neutron stars. The first finding is that, contrary to the common belief, the crust-core transition is not first order, and for all subsaturation densities matter can be viewed as a continuous fluid mixture between free nucleons and massive nuclei. As a consequence, the equations of state and the associated observables do not present any discontinuity over the whole thermodynamic range. We further investigate the nuclear matter composition over a wide range of densities and temperatures. At high density and temperature our model accounts for a much larger mass fraction bound in medium nuclei with respect to traditional approaches as Lattimer-Swesty, with sizeable consequences on the thermodynamic quantities. The equations of state agree well with the presently used EOS only at low temperatures and in the homogeneous matter phase, while important differences are present in the crust-core transition region. The correlation among the composition of baryonic matter and neutrino opacity is finally discussed, and we show that the two problems can be effectively decoupled. △ Less

Submitted 8 December, 2010; v1 submitted 12 September, 2010; originally announced September 2010.

Comments: 40 pages, 25 figures

Journal ref: Phys.Rev.C82:065801,2010

Isospin-dependent clusterization of Neutron-Star Matter

Authors: Camille Ducoin, Philippe Chomaz, Francesca Gulminelli

Abstract: Because of the presence of a liquid-gas phase transition in nuclear matter, compact-star matter can present a region of instability against the formation of clusters. We investigate this phase separation in a matter composed of neutrons, protons and electrons, within a Skyrme-Lyon mean-field approach. Matter instability and phase properties are characterized through the study of the free-energy… ▽ More Because of the presence of a liquid-gas phase transition in nuclear matter, compact-star matter can present a region of instability against the formation of clusters. We investigate this phase separation in a matter composed of neutrons, protons and electrons, within a Skyrme-Lyon mean-field approach. Matter instability and phase properties are characterized through the study of the free-energy curvature. The effect of beta-equilibrium is also analyzed in detail, and we show that the opacity to neutrinos has an influence on the presence of clusterized matter in finite-temperature proto-neutron stars. △ Less

Submitted 19 April, 2007; v1 submitted 11 December, 2006; originally announced December 2006.

Comments: To appear in Nuclear Physics A

Journal ref: Nucl.Phys.A789:403-425,2007

Anomalous thermodynamics and phase transitions of neutron-star matter

Authors: P. Chomaz, F. Gulminelli, C. Ducoin, P. Napolitani, K. H. O. Hasnaoui

Abstract: In this letter we show that the presence of the long-range Coulomb force in dense stellar matter implies that the total charge cannot be associated with a chemical potential, even if it is a conserved quantity. As a further consequence, the analytical properties of the partition sum are modified, changing the order of the phase transitions and affecting the possible occurrence of critical behavi… ▽ More In this letter we show that the presence of the long-range Coulomb force in dense stellar matter implies that the total charge cannot be associated with a chemical potential, even if it is a conserved quantity. As a further consequence, the analytical properties of the partition sum are modified, changing the order of the phase transitions and affecting the possible occurrence of critical behaviours. The peculiar thermodynamic properties of the system can be understood introducing a model hamiltonian in which each charge is independently neutralized by a uniform background of opposite charge. △ Less

Submitted 20 June, 2007; v1 submitted 27 July, 2005; originally announced July 2005.

Comments: 4 pages - GANIL, LPC

Report number: GANIL P 05 07

Journal ref: Phys.Rev.C75:065805,2007