-
Calibrating global behaviour of equation of state by combining nuclear and astrophysics inputs in a machine learning approach
Authors:
Sk Md Adil Imam,
Prafulla Saxena,
Tuhin Malik,
N. K. Patra,
B. K. Agrawal
Abstract:
We implemented symbolic regression techniques to identify suitable analytical functions that map various properties of neutron stars (NSs), obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) equations, to a few key parameters of the equation of state (EoS). These symbolic regression models (SRMs) are then employed to perform Bayesian inference with a comprehensive dataset from nuclear physic…
▽ More
We implemented symbolic regression techniques to identify suitable analytical functions that map various properties of neutron stars (NSs), obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) equations, to a few key parameters of the equation of state (EoS). These symbolic regression models (SRMs) are then employed to perform Bayesian inference with a comprehensive dataset from nuclear physics experiments and astrophysical observations. The posterior distributions of EoS parameters obtained from Bayesian inference using SRMs closely match those obtained directly from the solutions of TOV equations. Our SRM-based approach is approximately 100 times faster, enabling efficient Bayesian analyses across different combinations of data to explore their sensitivity to various EoS parameters within a reasonably short time.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
Unraveling the global behavior of equation of state by explicit finite nuclei constraints
Authors:
Anagh Venneti,
Sakshi Gautam,
Sarmistha Banik,
B. K. Agrawal
Abstract:
We obtain posterior distribution of equations of state (EOSs) across a broad range of density by imposing explicitly the constraints from precisely measured fundamental properties of finite nuclei, in combination with the experimental data from heavy-ion collisions and the astrophysical observations of radius, tidal deformability and minimum-maximum mass of neutron stars. The acquired EOSs exhibit…
▽ More
We obtain posterior distribution of equations of state (EOSs) across a broad range of density by imposing explicitly the constraints from precisely measured fundamental properties of finite nuclei, in combination with the experimental data from heavy-ion collisions and the astrophysical observations of radius, tidal deformability and minimum-maximum mass of neutron stars. The acquired EOSs exhibit a distinct global behavior compared to those usually obtained by imposing the finite nuclei constraints implicitly through empirical values of selected key parameters describing symmetric nuclear matter and symmetry energy in the vicinity of the saturation density. The explicit treatment of finite nuclei constraints yields softer EOSs at low densities which eventually become stiffer to meet the maximum mass criteria. The Kullback-Leibler divergence has been used to perform a quantitative comparison of the distributions of neutron star properties resulting from the EOSs obtained from implicit and explicit finite nuclei constraints.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
Implications of comprehensive nuclear and astrophysics data on the equations of state of neutron star matter
Authors:
Sk Md Adil Imam,
Tuhin Malik,
Constança Providência,
B. K. Agrawal
Abstract:
The equations of state (EoSs) governing neutron star (NS) matter obtained for both non-relativistic and relativistic mean-field models are systematically confronted with a diverse set of terrestrial data and astrophysical observations within the Bayesian framework. The terrestrial data, spans from bulk properties of finite nuclei to the heavy-ion collisions, constrain the symmetric nuclear matter…
▽ More
The equations of state (EoSs) governing neutron star (NS) matter obtained for both non-relativistic and relativistic mean-field models are systematically confronted with a diverse set of terrestrial data and astrophysical observations within the Bayesian framework. The terrestrial data, spans from bulk properties of finite nuclei to the heavy-ion collisions, constrain the symmetric nuclear matter EoS and the symmetry energy up to twice the saturation density ($ρ_0$= 0.16 fm$^{-3}$). The astrophysical observations encompass the NS radius, the tidal deformability, and the lower bound on maximum mass. Three distinct posterior distributions of EoSs are generated by gradually updating the priors with different constraints: (i) only the maximum NS mass, (ii) incorporating additional terrestrial data, (iii) combining both the terrestrial data and astrophysical observations. These EoS distributions are then compared using the Kullback-Liebler divergence which highlights the significant constraints imposed on the EoSs by the currently available lower bound of NS maximum mass and terrestrial data. The remaining astrophysical observations marginally refine the EoS within the density range $\sim$ 2-3$ρ_0$. It is observed that the relativistic mean field model yields stiffer EoS around the saturation density, but predict smaller values of the speed of sound and proton fraction in the interior of massive stars.
△ Less
Submitted 16 May, 2024; v1 submitted 11 January, 2024;
originally announced January 2024.
-
Establishing connection between neutron star properties and nuclear matter parameters through a comprehensive multivariate analysis
Authors:
N. K. Patra,
Prafulla Saxena,
B. K. Agrawal,
T. K. Jha
Abstract:
We have attempted to mitigate the challenge of connecting the neutron star (NS) properties with the nuclear matter parameters that describe equations of state (EoSs). The efforts to correlate various neutron star properties with individual nuclear matter parameters have been inconclusive. A Principal Component Analysis is employed as a tool to uncover the connection between multiple nuclear matter…
▽ More
We have attempted to mitigate the challenge of connecting the neutron star (NS) properties with the nuclear matter parameters that describe equations of state (EoSs). The efforts to correlate various neutron star properties with individual nuclear matter parameters have been inconclusive. A Principal Component Analysis is employed as a tool to uncover the connection between multiple nuclear matter parameters and the tidal deformability as well as the radius of neutron stars within the mass range of $1.2-1.8M_\odot$. The essential EOSs for neutron star matter at low densities have been derived using both uncorrelated uniform distributions and minimally constrained joint posterior distributions of nuclear matter parameters. For higher densities ($ρ> 0.32$fm$^{-3}$), the EOSs have been established through a suitable parameterization of the speed of sound, which consistently maintains causality and gradually approaches the conformal limit. Our analysis reveals that in order to account for over 90\% of the variability in NS properties, it is crucial to consider two or more principal components, emphasizing the significance of employing multivariate analysis. To explain the variability in tidal deformability needs a greater number of principal components compared to those for the radius at a given NS mass. The contributions from iso-vector nuclear matter parameters to the tidal deformability and radius of NS decrease by $\sim$ 25\% with the increase in mass of NS from 1.2$M_\odot$ to 1.8$M_\odot$.
△ Less
Submitted 10 December, 2023; v1 submitted 26 August, 2023;
originally announced August 2023.
-
Direct mapping of tidal deformability to the iso-scalar and iso-vector nuclear matter parameters
Authors:
Sk Md Adil Imam,
Arunava Mukherjee,
B. K. Agrawal,
Gourab Banerjee
Abstract:
Background: The equations of state (EoSs) which determine the properties of neutron stars (NSs) are often characterized by the iso-scalar and iso-vector nuclear matter parameters (NMPs). Recent attempts to relate the radius and tidal deformability of a NS to the individual NMPs have been inconclusive. These properties display strong correlations with the pressure of NS matter which depends on seve…
▽ More
Background: The equations of state (EoSs) which determine the properties of neutron stars (NSs) are often characterized by the iso-scalar and iso-vector nuclear matter parameters (NMPs). Recent attempts to relate the radius and tidal deformability of a NS to the individual NMPs have been inconclusive. These properties display strong correlations with the pressure of NS matter which depends on several NMPs. The knowledge of minimal NMPs that determine the NS properties will be necessary to address any connection between NS properties (e.g., tidal deformability) and that of finite nuclei.
Purpose: To identify the important NMPs required to describe the tidal deformability of neutron star for astrophysically relevant range of their gravitational masses (1.2 -- 1.8 M$_\odot$) as encountered in the binary neutron star merger events.
Method: We construct a large set of EoSs using four iso-scalar and five iso-vector NMPs. These EOSs are employed to perform a systematic analysis to isolate the NMPs that predominantly determine the tidal deformability, over a wide range of NS mass. The tidal deformability is then directly mapped to these NMPs.
Results: The tidal deformability of the NS with mass 1.2-1.8 M$_\odot$ can be determined within 10$\%$ directly in terms of four nuclear matter parameters, namely, the incompressibility $K_0$ and skewness $Q_0$ of symmetric nuclear matter, and the slope $L_0$ and curvature parameter $K_{\rm sym,0}$ of symmetry energy.
Conclusion: A function that quickly estimates the value of tidal deformability in terms of minimal nuclear matter parameters is developed. Our method can also be extended to other NS observables.
△ Less
Submitted 24 February, 2024; v1 submitted 18 May, 2023;
originally announced May 2023.
-
CREX- and PREX-II-motivated relativistic interactions and their implications for the bulk properties of nuclear matter and neutron stars
Authors:
Mukul Kumar,
Sunil Kumar,
Virender Thakur,
Raj Kumar,
B. K. Agrawal,
Shashi K. Dhiman
Abstract:
We investigate the implications of parity-violating electron scattering experiment on neutron skin thickness of $^{48}$Ca (CREX) and $^{208}$Pb (PREX-II) data on the bulk properties of finite nuclei, nuclear matter, and neutron stars. The neutron skin thickness from the CREX and PREX-II data is employed to constrain the parameters of relativistic mean field models which includes different non-line…
▽ More
We investigate the implications of parity-violating electron scattering experiment on neutron skin thickness of $^{48}$Ca (CREX) and $^{208}$Pb (PREX-II) data on the bulk properties of finite nuclei, nuclear matter, and neutron stars. The neutron skin thickness from the CREX and PREX-II data is employed to constrain the parameters of relativistic mean field models which includes different non-linear, self and cross-couplings among isoscalar-scalar $σ$, isoscalar-vector $ω$, isovector-scalar $δ$ and isovector-vector $ρ$ meson fields up to the quartic order. Three parametrizations of RMF model are proposed by fitting CREX, PREX-II and both CREX as well as PREX-II data to assess their implications. A covariance analysis is performed to assess the theoretical uncertainties of model parameters and nuclear matter observables along with correlations among them. The RMF model parametrization obtained with the CREX data acquires much smaller value of symmetry energy (J= 28.97$\pm$ 0.99 MeV), its slope parameter (L= 30.61$\pm 6.74$ MeV) in comparison to those obtained with PREX-II data. The neutron star properties are studied by employing the equations of state (EoSs) composed of nucleons and leptons in $β$ equilibrium.
△ Less
Submitted 10 May, 2023;
originally announced May 2023.
-
Neutron Skin Thickness Dependence of Astrophysical $S$-factor
Authors:
T. Ghosh,
Sangeeta,
G. Saxena,
B. K. Agrawal,
Ushasi Datta
Abstract:
Background: The density dependence of nuclear symmetry energy is crucial in determining several properties of finite nuclei to the neutron stars with mass $\sim$ 1.4 $M_\odot$. The values of neutron skin thickness, isovector giant dipole resonances energies and various nuclear reaction cross-sections in asymmetric nuclei have been utilized to determine the slope of symmetry energy ($L_0$) at the s…
▽ More
Background: The density dependence of nuclear symmetry energy is crucial in determining several properties of finite nuclei to the neutron stars with mass $\sim$ 1.4 $M_\odot$. The values of neutron skin thickness, isovector giant dipole resonances energies and various nuclear reaction cross-sections in asymmetric nuclei have been utilized to determine the slope of symmetry energy ($L_0$) at the saturation density. Recent PREX-II and CREX measurements of neutron skin thickness in $^{208}$Pb and $^{48}$Ca nuclei yield very different values of $L_0$ which overlap marginally within 90$\%$ confidence interval.
Purpose: Our objective is to demonstrate the role of symmetry energy on the sub-barrier fusion cross-section and the astrophysical $S$-factor for asymmetric nuclei.
Method: The nucleus nucleus potentials are generated using the double folding model (DFM) for three different nucleon-nucleon interactions. These DFM potentials are used for the calculation of the sub-barrier fusion cross-section and the astrophysical $S$-factor. The nucleon densities required for DFM potentials are generated from different families of non-relativistic and relativistic mean-field models which correspond to a wide range of neutron skin thickness or $L_0$.
Results: We have calculated the sub-barrier fusion cross-section for several asymmetric nuclei involving O, Ca, Ni, and Sn isotopes. The results are presented for the barrier parameters, cross-section, and astrophysical $S$-factor for $^{54}$Ca+$^{54}$Ca and $^{124}$Sn+$^{124}$Sn as a function of neutron skin thickness.
Conclusions: The cross-section for the neutron-rich nuclei show a strong dependence on the behavior of symmetry energy or the neutron skin thickness. The increase in skin thickness lowers the height of the barrier as well as its width which enhances the values of the $S$-factor by more than an order of magnitude.
△ Less
Submitted 21 March, 2023;
originally announced March 2023.
-
Systematic analysis of the impacts of symmetry energy parameters on neutron star properties
Authors:
N. K. Patra,
Anagh Venneti,
Sk Md Adil Imam,
Arunava Mukherjee,
B. K. Agrawal
Abstract:
The impacts of various symmetry energy parameters on the properties of neutron stars (NSs) have been recently investigated, and the outcomes are at variance, as summarized in Table III of Phys. Rev. D 106, 063005 (2022). We have systematically analyzed the correlations of slope and curvature parameters of symmetry energy at the saturation density ($ρ_0=0.16 \text{fm}^{-3}$) with the tidal deformab…
▽ More
The impacts of various symmetry energy parameters on the properties of neutron stars (NSs) have been recently investigated, and the outcomes are at variance, as summarized in Table III of Phys. Rev. D 106, 063005 (2022). We have systematically analyzed the correlations of slope and curvature parameters of symmetry energy at the saturation density ($ρ_0=0.16 \text{fm}^{-3}$) with the tidal deformability and stellar radius of non-spinning neutron stars in the mass range of $1.2 - 1.6 M_\odot$ using a large set of minimally constrained equations of state (EoSs). The EoSs at low densities correspond to the nucleonic matter and are constrained by empirical ranges of a few low-order nuclear matter parameters from the finite nuclei data and the pure neutron matter EoS from chiral effective field theory. The EoSs at high densities ($ρ> 1.5 - 2ρ_0$) are obtained by a parametric form for the speed of sound that satisfies the causality condition. Several factors affecting the correlations between the NS properties and the individual symmetry energy parameters usually encountered in the literature are considered. These correlations are quite sensitive to the choice of the distributions of symmetry energy parameters and their interdependence. But, variations of NS properties with the pressure of $β-$ equilibrated matter at twice the saturation density remain quite robust which maybe due to the fact that the pressure depends on the combination of multiple nuclear matter parameters that describe the symmetric nuclear matter as well as the density dependence of the symmetry energy. Our results are practically insensitive to the behavior of EoS at high densities.
△ Less
Submitted 18 May, 2023; v1 submitted 8 February, 2023;
originally announced February 2023.
-
Inferring the nuclear symmetry energy at supra saturation density from neutrino cooling
Authors:
Tuhin Malik,
B. K. Agrawal,
Constança Providência
Abstract:
An ambitious goal of the astrophysical community is not only to constrain the equation of state (EOS) of neutron star (NS) matter by confronting it with astrophysics observations, but ultimately also to infer the NS composition. Nevertheless, the composition of the NS core is likely to remain uncertain unless we have an accurate determination of the nuclear symmetry energy at supra saturation dens…
▽ More
An ambitious goal of the astrophysical community is not only to constrain the equation of state (EOS) of neutron star (NS) matter by confronting it with astrophysics observations, but ultimately also to infer the NS composition. Nevertheless, the composition of the NS core is likely to remain uncertain unless we have an accurate determination of the nuclear symmetry energy at supra saturation density ($ρ>ρ_0$). We investigate how the nucleonic direct Urca (dUrca) processes can be used as an effective probe to constraint the high density nuclear symmetry energy. A large number of minimally constrained EOSs has been constructed by applying a Bayesian approach to study the correlations of the symmetry energy at different densities with a few selected properties of a NS. The nuclear symmetry energy above the baryon density 0.5 fm$^{-3}$ ($\sim 3 ρ_0$) is found to be strongly correlated with NS mass at which the onset of nucleonic dUrca neutrino cooling takes place in the core. This allows us to constrain the high density behavior of nuclear symmetry energy within narrow bounds. {The pure neutron matter pressure constraint from chiral effective field theory rules out the onset of nucleonic dUrca in stars with a mass $\lesssim$ 1.4 $M_\odot$.} The onset of dUrca inside 1.6 M$_\odot$ to 1.8 M$_\odot$ NS implies a slope of the symmetry energy $L$ at $\sim 2.5~ρ_0$, respectively, between 54 and 48 MeV.
△ Less
Submitted 3 July, 2022; v1 submitted 30 June, 2022;
originally announced June 2022.
-
Astrophysical reaction rates with realistic nuclear level densities
Authors:
Sangeeta,
T. Ghosh,
B. Maheshwari,
G. Saxena,
B. K. Agrawal
Abstract:
Realistic nuclear level densities (NLDs) obtained within the spectral distribution method (SDM) are employed to study nuclear processes of astrophysical interest. The merit of SDM lies in the fact that the NLDs corresponding to many body shell model Hamiltonian consisting of residual interaction can be obtained for the full configurational space without recourse to the exact diagnolization of huge…
▽ More
Realistic nuclear level densities (NLDs) obtained within the spectral distribution method (SDM) are employed to study nuclear processes of astrophysical interest. The merit of SDM lies in the fact that the NLDs corresponding to many body shell model Hamiltonian consisting of residual interaction can be obtained for the full configurational space without recourse to the exact diagnolization of huge matrices. We calculate NLDs and s-wave neutron resonance spacings which agree reasonably well with the available experimental data. By employing these NLDs, we compute reaction cross-sections and astrophysical reaction rates for radiative neutron capture in few Fe-group nuclei, and compare them with experimental data as well as with those obtained with NLDs from phenomenological and microscopic mean-field models. The results obtained for the NLDs from SDM are able to explain the experimental data quite well. These results are of particular importance since the configuration mixing through the residual interaction naturally accounts for the collective excitations. In the mean-field models, the collective effects are included through the vibrational and rotational enhancement factors and their NLDs are further normalized at low energies with neutron resonance data.
△ Less
Submitted 27 March, 2022;
originally announced March 2022.
-
Nearly model-independent constraints on dense matter equation of state in a Bayesian approach
Authors:
N. K. Patra,
Sk Md Adil Imam,
B. K. Agrawal,
Arunava Mukherjee,
Tuhin Malik
Abstract:
We apply Bayesian approach to construct a large number of minimally constrained equations of state (EOSs) and study their correlations with a few selected properties of a neutron star (NS). Our set of minimal constraints includes a few basic properties of saturated nuclear matter and low-density pure neutron matter EOS which is obtained from a precise next-to-next-to-next-to-leading-order (N…
▽ More
We apply Bayesian approach to construct a large number of minimally constrained equations of state (EOSs) and study their correlations with a few selected properties of a neutron star (NS). Our set of minimal constraints includes a few basic properties of saturated nuclear matter and low-density pure neutron matter EOS which is obtained from a precise next-to-next-to-next-to-leading-order (N$^{3}$LO) calculation in chiral effective field theory. The tidal deformability and radius of NS with mass $1-2 M_\odot$ are found to be strongly correlated with the pressure of $β$-equilibrated matter at densities higher than the saturation density ($ρ_0 = 0.16$ fm$^{-3}$) in a nearly model-independent manner. These correlations are employed to parametrize the pressure for $β$-equilibrated matter, around 2$ρ_0$, as a function of neutron star mass and the corresponding tidal deformability. The maximum mass of neutron star is also found to be strongly correlated with the pressure of $β$-equilibrated matter at densities $\sim 4.5ρ_0$.
△ Less
Submitted 9 August, 2022; v1 submitted 16 March, 2022;
originally announced March 2022.
-
Relativistic description of dense matter equation of state and compatibility with neutron star observables: a Bayesian approach
Authors:
Tuhin Malik,
Márcio Ferreira,
B. K. Agrawal,
Constança Providência
Abstract:
The general behavior of the nuclear equation of state (EOS), relevant for the description of neutron stars (NS), is studied within a Bayesian approach applied to a set of models based on a density dependent relativistic mean field description of nuclear matter. The EOS is subjected to a minimal number of constraints based on nuclear saturation properties and the low density pure neutron matter EOS…
▽ More
The general behavior of the nuclear equation of state (EOS), relevant for the description of neutron stars (NS), is studied within a Bayesian approach applied to a set of models based on a density dependent relativistic mean field description of nuclear matter. The EOS is subjected to a minimal number of constraints based on nuclear saturation properties and the low density pure neutron matter EOS obtained from a precise next-to-next-to-next-to-leading order (N$^{3}$LO) calculation in chiral effective field theory ($χ$EFT). The posterior distributions of the model parameters obtained under these minimal constraints are employed to construct the distributions of various nuclear matter properties and NS properties such as radii, tidal deformabilites, central energy densities and speeds of sound etc. We found that 90% confidence interval (CI) for allowed NS mass - radius relationship and tidal deformabilites are compatible with GW170817 and recent NICER observations, without invoking the exotic degrees of freedom. A central speed-of-sound of the order of $\sqrt{2/3}$ $c$ is obtained. The maximum neutron star mass allowed by the model is 2.5 $M_\odot$.
△ Less
Submitted 1 March, 2022; v1 submitted 29 January, 2022;
originally announced January 2022.
-
Bayesian reconstruction of nuclear matter parameters from the equation of state of neutron star matter
Authors:
Sk Md Adil Imam,
N. K. Patra,
C. Mondal,
Tuhin Malik,
B. K. Agrawal
Abstract:
The nuclear matter parameters (NMPs), those underlie in the construction of the equation of state (EoS) of neutron star matter, are not directly accessible. The Bayesian approach is applied to reconstruct the posterior distributions of NMPs from the EoS of neutron star matter. The constraints on lower-order parameters as imposed by the finite nuclei observables are incorporated through appropriate…
▽ More
The nuclear matter parameters (NMPs), those underlie in the construction of the equation of state (EoS) of neutron star matter, are not directly accessible. The Bayesian approach is applied to reconstruct the posterior distributions of NMPs from the EoS of neutron star matter. The constraints on lower-order parameters as imposed by the finite nuclei observables are incorporated through appropriately chosen prior distributions. The calculations are performed with two sets of pseudo data on the EoS whose true models are known. The median values of second or higher order NMPs show sizeable deviations from their true values and associated uncertainties are also larger. The sources of these uncertainties are intrinsic in nature, identified as (i) the correlations among various NMPs and (ii) the variations in the EoS of symmetric nuclear matter, symmetry energy, and the neutron-proton asymmetry in such a way that the neutron star matter EoS remain almost unaffected.
△ Less
Submitted 22 June, 2022; v1 submitted 29 October, 2021;
originally announced October 2021.
-
Unveiling the correlations of tidal deformability with the nuclear symmetry energy parameters
Authors:
Tuhin Malik,
B. K. Agrawal,
Constança Providência,
J. N. De
Abstract:
The chi-squared based covariance approach allows one to estimate the correlations among desired observables related to nuclear matter directly from a set of fit data without taking recourse to the distributions of the nuclear matter parameters (NMPs). Such an approach is applied to study the correlations of tidal deformability of neutron star with the slope and the curvature parameters of nuclear…
▽ More
The chi-squared based covariance approach allows one to estimate the correlations among desired observables related to nuclear matter directly from a set of fit data without taking recourse to the distributions of the nuclear matter parameters (NMPs). Such an approach is applied to study the correlations of tidal deformability of neutron star with the slope and the curvature parameters of nuclear symmetry energy governed by an extensive set of fit data on the finite nuclei together with the maximum mass of the neutron star. The knowledge of the distributions of NMPs consistent with the fit data is implicitly inbuilt in the Hessian matrix which is central to this covariance approach. Comparing our results with those obtained with the explicit use of the distributions of NMPs, we show that the appropriate correlations among NMPs as induced by the fit data are instrumental in strengthening the correlations of the tidal deformability with the symmetry energy parameters, without it, the said correlations tend to disappear. The interplay between isoscalar and isovector NMPs is also emphasized.
△ Less
Submitted 8 August, 2020;
originally announced August 2020.
-
Constraining nuclear matter parameters from correlation systematics:a mean-field perspective
Authors:
B. K. Agrawal,
Tuhin Malik,
J. N. De,
S. K. Samaddar
Abstract:
The nuclear matter parameters define the nuclear equation of state (EoS), they appear as coefficients of expansion around the saturation density of symmetric and asymmetric nuclear matter. We review their correlations with several properties of finite nuclei and of neutron stars within mean-field frameworks. The lower order nuclear matter parameters such as the binding energy per nucleon, incompre…
▽ More
The nuclear matter parameters define the nuclear equation of state (EoS), they appear as coefficients of expansion around the saturation density of symmetric and asymmetric nuclear matter. We review their correlations with several properties of finite nuclei and of neutron stars within mean-field frameworks. The lower order nuclear matter parameters such as the binding energy per nucleon, incompressibility and the symmetry energy coefficients are found to be constrained in narrow limits through their strong ties with selective properties of finite nuclei. From the correlations of nuclear matter parameters with neutron star observables, we further review how precision knowledge of the radii and tidal deformability of neutron stars in the mass range $1 - 2 M_\odot$ may help cast them in narrower bounds. The higher order parameters such as the density slope and the curvature of the symmetry energy or the skewness of the symmetric nuclear matter EoS are, however, plagued with larger uncertainty. From inter-correlation of these higher order nuclear matter parameters with lower order ones, we explore how they can be brought to more harmonious bounds.
△ Less
Submitted 10 June, 2020;
originally announced June 2020.
-
Empirical constraints on the high-density equation of state from multi-messenger observables
Authors:
Márcio Ferreira,
M. Fortin,
Tuhin Malik,
B. K. Agrawal,
Constança Providência
Abstract:
We search for possible correlations between neutron star observables and thermodynamic quantities that characterize high density nuclear matter. We generate a set of model-independent equations of state describing stellar matter from a Taylor expansion around saturation density. Each equation of state which is a functional of the nuclear matter parameters is thermodynamically consistent, causal an…
▽ More
We search for possible correlations between neutron star observables and thermodynamic quantities that characterize high density nuclear matter. We generate a set of model-independent equations of state describing stellar matter from a Taylor expansion around saturation density. Each equation of state which is a functional of the nuclear matter parameters is thermodynamically consistent, causal and compatible with astrophysical observations. We find that the neutron star tidal deformability and radius are strongly correlated with the pressure, the energy density and the sound velocity at different densities. Similar correlations are also exhibited by a large set of mean-field models based on non-relativistic and relativistic nuclear energy density functionals. These model independent correlations can be employed to constrain the equation of state at different densities above saturation from measurements of NS properties with multi-messenger observations. In particular, precise constraints on the radius of PSR J0030+0451 thanks to NICER observations would allow to better infer the properties of matter around two times the nuclear saturation density.
△ Less
Submitted 23 December, 2019;
originally announced December 2019.
-
GW170817: constraining the nuclear matter equation of state from the neutron star tidal deformability
Authors:
Tuhin Malik,
N. Alam,
M. Fortin,
C. Providência,
B. K. Agrawal,
T. K. Jha,
Bharat Kumar,
S. K. Patra
Abstract:
Constraints set on key parameters of the nuclear matter equation of state (EoS) by the values of the tidal deformability, inferred from GW170817, are examined by using a diverse set of relativistic and non-relativistic mean field models. These models are consistent with bulk properties of finite nuclei as well as with the observed lower bound on the maximum mass of neutron star…
▽ More
Constraints set on key parameters of the nuclear matter equation of state (EoS) by the values of the tidal deformability, inferred from GW170817, are examined by using a diverse set of relativistic and non-relativistic mean field models. These models are consistent with bulk properties of finite nuclei as well as with the observed lower bound on the maximum mass of neutron star $\sim 2 ~ {\rm M}_\odot$. The tidal deformability shows a strong correlation with specific linear combinations of the isoscalar and isovector nuclear matter parameters associated with the EoS. Such correlations suggest that a precise value of the tidal deformability can put tight bounds on several EoS parameters, in particular, on the slope of the incompressibility and the curvature of the symmetry energy. The tidal deformability obtained from the GW170817 and its UV/optical/infrared counterpart sets the radius of a canonical $1.4~ {\rm M}_{\odot}$ neutron star to be $11.82\leqslant R_{1.4}\leqslant13.72$ km.
△ Less
Submitted 29 September, 2018; v1 submitted 30 May, 2018;
originally announced May 2018.
-
New relativistic effective interaction for finite nuclei, infinite nuclear matter and neutron stars
Authors:
Bharat Kumar,
S. K. Patra,
B. K. Agrawal
Abstract:
We carry out the study for finite nuclei, infinite nuclear matter and neutron star properties with the newly developed relativistic force named as the Institute Of Physics Bhubaneswar-I(IOPB-I).
Using this force, we calculate the binding energies, charge radii and neutron skin thickness for some selected nuclei. From the ground state properties of superheavy element i.e. Z=120, it is noticed tha…
▽ More
We carry out the study for finite nuclei, infinite nuclear matter and neutron star properties with the newly developed relativistic force named as the Institute Of Physics Bhubaneswar-I(IOPB-I).
Using this force, we calculate the binding energies, charge radii and neutron skin thickness for some selected nuclei. From the ground state properties of superheavy element i.e. Z=120, it is noticed that considerable shell gaps appear at neutron numbers N=172, 184 and 198, showing the magicity of these numbers. The low density behavior of the equation of state for pure neutron matter is compatible with other microscopic models. Along with the nuclear symmetry energy, its slope and curvature parameters at the saturation density are consistent with those extracted from various experimental data. We calculate the neutron star properties with the equation of state composed of nucleons and leptons in $\it beta-equilibrium$ which are in good agreement with the X-ray observations by Steiner and Nättilä. We find that the maximum mass of the neutron star to be 2.15$M_{\odot}$ and stellar radius 11.936 km . Moreover, the radius and tidal deformability of a {\it canonical} neutron star mass 1.4$M_\odot$ come out to be 13.242 km and 3.910$\times$10$^{36}$ g cm$^2$ s$^2$ respectively within this parameter set.
△ Less
Submitted 7 February, 2018; v1 submitted 13 November, 2017;
originally announced November 2017.
-
Warm unstable asymmetric nuclear matter: critical properties and the density dependence of the symmetry energy
Authors:
N. Alam,
H. Pais,
C. Providência,
B. K. Agrawal
Abstract:
The spinodal instabilities in hot asymmetric nuclear matter and some important critical parameters derived thereof are studied using six different families of relativistic mean-field (RMF) models. The slopes of the symmetry energy coefficient vary over a wide range within each family. The critical densities and proton fractions are more sensitive to the symmetry energy slope parameter at temperatu…
▽ More
The spinodal instabilities in hot asymmetric nuclear matter and some important critical parameters derived thereof are studied using six different families of relativistic mean-field (RMF) models. The slopes of the symmetry energy coefficient vary over a wide range within each family. The critical densities and proton fractions are more sensitive to the symmetry energy slope parameter at temperatures much below its critical value ($T_c\sim$14-16 MeV). The spread in the critical proton fraction at a given symmetry energy slope parameter is noticeably larger near $T_c$, indicating that the warm equation of state of asymmetric nuclear matter at sub-saturation densities is not sufficiently constrained. The distillation effects are sensitive to the density dependence of the symmetry energy at low temperatures which tend to wash out with increasing temperature.
△ Less
Submitted 5 May, 2017;
originally announced May 2017.
-
Strong correlations of neutron star radii with the slopes of nuclear matter incompressibility and symmetry energy at saturation
Authors:
N. Alam,
B. K. Agrawal,
M. Fortin,
H. Pais,
C. Providência,
Ad. R. Raduta,
A. Sulaksono
Abstract:
We examine the correlations of neutron star radii with the nuclear matter incompressibility, symmetry energy, and their slopes, which are the key parameters of the equation of state (EoS) of asymmetric nuclear matter. The neutron star radii and the EoS parameters are evaluated using a representative set of 24 Skyrme-type effective forces and 18 relativistic mean field models, and two microscopic c…
▽ More
We examine the correlations of neutron star radii with the nuclear matter incompressibility, symmetry energy, and their slopes, which are the key parameters of the equation of state (EoS) of asymmetric nuclear matter. The neutron star radii and the EoS parameters are evaluated using a representative set of 24 Skyrme-type effective forces and 18 relativistic mean field models, and two microscopic calculations, all describing 2$M_\odot$ neutron stars. Unified EoSs for the inner-crust-core region have been built for all the phenomenological models, both relativistic and non-relativistic. Our investigation shows the existence of a strong correlation of the neutron star radii with the linear combination of the slopes of the nuclear matter incompressibility and the symmetry energy coefficients at the saturation density. Such correlations are found to be almost independent of the neutron star mass in the range $0.6\text{-}1.8M_{\odot}$. This correlation can be linked to the empirical relation existing between the star radius and the pressure at a nucleonic density between one and two times saturation density, and the dependence of the pressure on the nuclear matter incompressibility, its slope and the symmetry energy slope. The slopes of the nuclear matter incompressibility and the symmetry energy coefficients as estimated from the finite nuclei data yield the radius of a $1.4M_{\odot}$ neutron star in the range $11.09\text{-}12.86$ km.
△ Less
Submitted 20 October, 2016;
originally announced October 2016.
-
Diversities in the properties of neutron stars at a fixed neutron-skin thickness in $^{208}$Pb nucleus
Authors:
N. Alam,
A. Sulaksono,
B. K. Agrawal
Abstract:
We study the diversities in the properties of the neutron stars arising due to the different choices for the cross-coupling between various mesons which governs the density dependence of the nuclear symmetry energy in the extended relativistic mean-field(RMF) model. For this purpose, we obtain two different families of the extended RMF model corresponding to different non-linear cross-coupling ter…
▽ More
We study the diversities in the properties of the neutron stars arising due to the different choices for the cross-coupling between various mesons which governs the density dependence of the nuclear symmetry energy in the extended relativistic mean-field(RMF) model. For this purpose, we obtain two different families of the extended RMF model corresponding to different non-linear cross-coupling term in the isovector part of the effective Lagrangian density. The lowest order contributions for the $δ$ mesons are also included. The different models within the same family are so obtained that they yield wide variation in the value of neutron-skin thickness in the $^{208}$Pb nucleus. These models are employed to compute the neutron star properties such as, core-crust transition density, radius and red shift at canonical mass ($1.4M_{\odot}$), tidal polarizability parameter, and threshold mass required for the enhanced cooling through direct Urca process. Most of the neutron star properties considered are significantly different(10\%-40\%) for the different families of models at a smaller neutron-skin thickness ($\sim 0.15$ fm) in the $^{208}$Pb nucleus.
△ Less
Submitted 3 July, 2015;
originally announced July 2015.
-
Constraining the density dependence of symmetry energy from nuclear masses
Authors:
B. K. Agrawal,
J. N. De,
S. K. Samaddar,
G. Colo,
A. Sulaksono
Abstract:
Empirically determined values of the nuclear volume and surface symmetry energy coefficients from nuclear masses are expressed in terms of density distributions of nucleons in heavy nuclei in the local density approximation. This is then used to extract the value of the symmetry energy slope parameter $L$. The density distributions in both spherical and well deformed nuclei calculated within micro…
▽ More
Empirically determined values of the nuclear volume and surface symmetry energy coefficients from nuclear masses are expressed in terms of density distributions of nucleons in heavy nuclei in the local density approximation. This is then used to extract the value of the symmetry energy slope parameter $L$. The density distributions in both spherical and well deformed nuclei calculated within microscopic framework with different energy density functionals give $L = 59.0 \pm 13.0$ MeV. Application of the method also helps in a precision determination of the neutron skin thickness of nuclei that are difficult to measure accurately.
△ Less
Submitted 23 May, 2013;
originally announced May 2013.
-
Determining the density content of symmetry energy and neutron skin: an empirical approach
Authors:
B. K. Agrawal,
J. N. De,
S. K. Samaddar
Abstract:
The density dependence of nuclear symmetry energy remains poorly constrained. Starting from precise empirical values of the nuclear volume and surface symmetry energy coefficients and the nuclear saturation density, we show how in the ambit of microscopic calculations with different energy density functionals, the value of the symmetry energy slope parameter $L$ alongwith that for neutron skin can…
▽ More
The density dependence of nuclear symmetry energy remains poorly constrained. Starting from precise empirical values of the nuclear volume and surface symmetry energy coefficients and the nuclear saturation density, we show how in the ambit of microscopic calculations with different energy density functionals, the value of the symmetry energy slope parameter $L$ alongwith that for neutron skin can be put in tighter bounds. The value of $L$ is found to be $L$= 64$\pm $5 MeV. For $^{208}$Pb, the neutron skin thickness comes out to be 0.188 $\pm $0.014 fm. Knowing $L$, the method can be applied to predict neutron skins of other nuclei.
△ Less
Submitted 3 December, 2012;
originally announced December 2012.
-
Existence of hyperons in the pulsar PSRJ1614-2230
Authors:
A. Sulaksono,
B. K. Agrawal
Abstract:
The possibility of existence of hyperons in the recently measured $2M_\odot$ pulsar PSRJ1614-2230 is explored using a diverse set of nuclear equations of state calculated within the relativistic mean-field models. Our results indicate that the nuclear equations of state compatible with heavy-ion data allow the hyperons to exist in the PSRJ1614-2230 only for significantly larger values for the meso…
▽ More
The possibility of existence of hyperons in the recently measured $2M_\odot$ pulsar PSRJ1614-2230 is explored using a diverse set of nuclear equations of state calculated within the relativistic mean-field models. Our results indicate that the nuclear equations of state compatible with heavy-ion data allow the hyperons to exist in the PSRJ1614-2230 only for significantly larger values for the meson-hyperon coupling strengths. The maximum mass configurations for these cases contain sizable hyperon fractions ($\sim 60%$) and yet masquared their counterpart composed of only nucleonic matter.
△ Less
Submitted 27 September, 2012;
originally announced September 2012.
-
Electric dipole polarizability and the neutron skin
Authors:
J. Piekarewicz,
B. K. Agrawal,
G. Colo,
W. Nazarewicz,
N. Paar,
P. -G. Reinhard,
X. Roca-Maza,
D. Vretenar
Abstract:
The recent high-resolution measurement of the electric dipole (E1) polarizability (alphad) in 208Pb [Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness of this nucleus. The neutron-skin thickness (rskin) of 208Pb is a quantity of critical importance for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model dependence…
▽ More
The recent high-resolution measurement of the electric dipole (E1) polarizability (alphad) in 208Pb [Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness of this nucleus. The neutron-skin thickness (rskin) of 208Pb is a quantity of critical importance for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model dependence of the correlation between alphad and rskin, we carry out systematic calculations for 208Pb, 132Sn, and 48Ca based on the nuclear density functional theory (DFT) using both non-relativistic and relativistic energy density functionals (EDFs). Our analysis indicates that whereas individual models exhibit a linear dependence between alphad and rskin, this correlation is not universal when one combines predictions from a host of different models. By averaging over these model predictions, we provide estimates with associated systematic errors for rskin and alphad for the nuclei under consideration. We conclude that precise measurements of rskin in both 48Ca and 208Pb---combined with the recent measurement of alphad---should significantly constrain the isovector sector of the nuclear energy density functional.
△ Less
Submitted 18 January, 2012;
originally announced January 2012.
-
Equations of state and stability of color-superconducting quark matter cores in hybrid stars
Authors:
B. K. Agrawal
Abstract:
The stable configurations of non-rotating and rotating hybrid stars composed of colour superconducting quark matter core are constructed using several equations of state (EOSs). We use a set of diverse EOSs for the nuclear matter which represents the low density phase. The EOSs at higher densities correspond to the quark matter in the colour superconducting phase and are computed within the NJL-…
▽ More
The stable configurations of non-rotating and rotating hybrid stars composed of colour superconducting quark matter core are constructed using several equations of state (EOSs). We use a set of diverse EOSs for the nuclear matter which represents the low density phase. The EOSs at higher densities correspond to the quark matter in the colour superconducting phase and are computed within the NJL-like model for different values of the scalar diquark and vector current couplings strengths. The phase transition to the quark matter is computed by a Maxwell construction. We find that the stability of the hybrid stars are mainly governed by the behaviour of the EOSs for the colour superconducting quark matter. However the compositions of hybrid star are sensitive to the EOS of the nuclear matter. The value of the critical rotation frequency for the hybrid star depends strongly on the EOS of the nuclear matter as well as that for the colour superconducting quark matter. Our results indicate that the EOS for the colour superconducting quark matter can be obtained, by adjusting the parameters of the NJL model, to yield the stable configurations of the hybrid star having the maximum mass $\sim 1.5M_\odot$ in the non-rotating limit and the critical rotation frequency $\sim 1$ kHz.
△ Less
Submitted 10 January, 2010;
originally announced January 2010.
-
Stable configurations of hybrid stars with colour-flavour-locked core
Authors:
B. K. Agrawal,
Shashi K. Dhiman
Abstract:
We construct static and mass-shedding limit sequences of hybrid stars, composed of colour flavour locked (CFL) quark matter core, for a set of equations of state (EOSs). The EOS for the hadronic matter is obtained using appropriately calibrated extended field theoretical based relativistic mean-field model. The MIT bag model is employed to compute the EOSs of the CFL quark matter for different v…
▽ More
We construct static and mass-shedding limit sequences of hybrid stars, composed of colour flavour locked (CFL) quark matter core, for a set of equations of state (EOSs). The EOS for the hadronic matter is obtained using appropriately calibrated extended field theoretical based relativistic mean-field model. The MIT bag model is employed to compute the EOSs of the CFL quark matter for different values of the CFL gap parameter in the range of $50 - 150\text{MeV}$ with the deconfinement phase transition density ranging from $4ρ_0 - 6ρ_0$ ($ρ_0 = 0.16\text{fm}^{-3}$). We find, depending on the values of the CFL gap parameter and the deconfinement phase transition density, the sequences of stable configurations of hybrid stars either form third families of the compact stars or bifurcate from the hadronic sequence. The hybrid stars have masses $1.0 - 2.1 M_\odot$ with radii $9 - 13.5$ km. The maximum values of mass shedding limit frequency for such hybrid stars are $1 -2$ kHz. For the smaller values of the CFL gap parameter and the deconfinement phase transition density, mass-radius relationships are in harmony with those deduced by applying improved hydrogen atmosphere model to fit the high quality spectra from compact star X7 in the globular cluster 47 Tucanae. We observed for some cases that the third family of compact stars exist in the static sequence, but, disappear from the mass-shedding limit sequence. Our investigation suggests that the third family of compact stars in the mass-shedding limit sequence is more likely to appear, provided they have maximum mass in the static limit higher than their second family counterpart composed of pure hadronic matter.
△ Less
Submitted 20 April, 2009;
originally announced April 2009.
-
Correlations in the properties of static and rapidly rotating compact stars
Authors:
B. K. Agrawal,
Raj Kumar,
S. K. Dhiman
Abstract:
Correlations in the properties of the static compact stars (CSs) and the ones rotating with the highest observed frequency of 1122Hz are studied using a large set of equations of state (EOSs). These EOSs span various approaches and their chemical composition vary from the nucleons to hyperons and quarks in $β$-equilibrium. It is found that the properties of static CS, like, the maximum gravitati…
▽ More
Correlations in the properties of the static compact stars (CSs) and the ones rotating with the highest observed frequency of 1122Hz are studied using a large set of equations of state (EOSs). These EOSs span various approaches and their chemical composition vary from the nucleons to hyperons and quarks in $β$-equilibrium. It is found that the properties of static CS, like, the maximum gravitational mass $M_{\rm max}^{\rm stat}$ and radius $R_{1.4}^{\rm stat}$ corresponding to t he canonical mass and supramassive or non-supramassive nature of the CS rotating at 1122 Hz are strongly correlated. In particular, only those EOSs yield the CS rotating at 1122Hz to be non-supramassive for which $\left (\frac{M_{\rm max}^{\rm stat}}{M_\odot}\right )^{1/2} \left (\frac{10{\rm km}}{R_{1.4}^{\rm stat}})^{3/2} $ is greater than unity. Suitable parametric form which can be used to split the $M_{\rm max}^{\rm stat}$ $-$ $R_{1.4}^{\rm stat}$ plane into the regions of different supramassive nature of the CS rotating at 1122Hz is presented. Currently measured maximum gravitational mass 1.76$M_\odot$ of PSR J0437-4715 suggests that the CS rotating at 1122Hz can be non-supramassive provided $R_{1.4}^{\rm stat} \leqslant 12.4$ km.
△ Less
Submitted 22 February, 2008;
originally announced February 2008.
-
Non-rotating and rotating neutron stars in the extended field theoretical model
Authors:
Shashi K. Dhiman,
Raj Kumar,
B. K. Agrawal
Abstract:
We study the properties of non-rotating and rotating neutron stars for a new set of equations of state (EOSs) with different high density behaviour obtained using the extended field theoretical model. The high density behaviour for these EOSs are varied by varying the $ω-$meson self-coupling and hyperon-meson couplings in such a way that the quality of fit to the bulk nuclear observables, nuclea…
▽ More
We study the properties of non-rotating and rotating neutron stars for a new set of equations of state (EOSs) with different high density behaviour obtained using the extended field theoretical model. The high density behaviour for these EOSs are varied by varying the $ω-$meson self-coupling and hyperon-meson couplings in such a way that the quality of fit to the bulk nuclear observables, nuclear matter incompressibility coefficient and hyperon-nucleon potential depths remain practically unaffected. We find that the largest value for maximum mass for the non-rotating neutron star is $2.1M_\odot$. The radius for the neutron star with canonical mass is $12.8 - 14.1$ km provided only those EOSs are considered for which maximum mass is larger than $1.6M_\odot$ as it is the lower bound on the maximum mass measured so far. Our results for the very recently discovered fastest rotating neutron star indicate that this star is supra massive with mass $1.7 - 2.7M_\odot$ and circumferential equatorial radius $12 - 19$ km.
△ Less
Submitted 26 September, 2007;
originally announced September 2007.