Engineering topics
Somasundaram, Rahul
Publications and source records attributed to Somasundaram, Rahul.
Impact of dark matter on tidal signatures in neutron star mergers with the Einstein Telescope
If dark matter (DM) accumulates inside neutron stars (NS), it changes their internal structure and causes a shift of the tidal deformability from the value predicted by the dense-matter equation of state (EOS). In principle, this shift could be observable in the gravitational-wave (GW) signal of binary neutron star (BNS) mergers. We investigate the effect of fermionic, noninteracting DM when observing a large number of GW events from DM-admixed BNSs with the precision of the proposed Einstein telescope (ET). Specifically, we study the impact on the recovery of the baryonic EOS and whether DM properties can be constrained. For this purpose, we create event catalogs of BNS mock events with DM fraction up to 1%, from which we reconstruct the posterior uncertainties with the Fisher matrix approach. Using this data, we perform joint Bayesian inference on the baryonic EOS, DM particle mass, and DM particle fraction in each event. Here, our results reveal that when falsely ignoring DM effects, the EOS posterior is biased toward softer EOSs, though the offset is rather small. Further, we find that within our assumptions of our DM model and population, ET will likely not be able to test the presence of DM in BNSs, even when combining many events and adding Cosmic Explorer (CE) to the next-generation detector network. Likewise, the potential constraints on the DM particle mass will remain weak because of degeneracies with the fraction and EOS.
Classification of compact objects and model comparison using EOS knowledge
Nuclear theory and experiments, alongside astrophysical observations, constrain the equation of state (EOS) of supranuclear-dense matter. Conversely, knowledge of the EOS allows an improved interpretation of nuclear or astrophysical data. In this article, we use several established constraints on the EOS and the new NICER measurement of PSR J0437-4715 to comment on the nature of the primary companion in GW230529 and the companion of PSR J0514-4002E. We find that, with a probability of ≳84% and ≳68%, respectively, both objects are black holes. These likelihoods increase to above 95% when one uses GW170817’s remnant as an upper limit on the TOV mass. We also demonstrate that the current knowledge of the EOS substantially disfavors high masses and radii for PSR J0030+0451, inferred recently when combining NICER with XMM-Newton background data and using particular hot-spot models. Lastly, we also use our obtained EOS knowledge to comment on measurements of the nuclear symmetry energy, finding that the large value predicted by the PREX-II measurement displays some mild tension with other constraints on the EOS.
Measuring neutron star radius with second and third generation gravitational wave detector networks
The next generation of ground-based interferometric gravitational wave detectors will observe mergers of black holes and neutron stars throughout cosmic time. A large number of the binary neutron star merger events will be observed with extreme high fidelity, and will provide stringent constraints on the equation of state of nuclear matter. In this paper, we investigate the systematic improvement in the measurability of the equation of state with increase in detector sensitivity by combining constraints obtained on the radius of a 1.4 M ⊙ neutron star from a simulated source population. Since the measurability of the equation of state depends on its stiffness, we consider a range of realistic equations of state that span the current observational constraints. We show that a single 40 km Cosmic Explorer detector can pin down the neutron star radius for a soft, medium and stiff equation of state with a precision of 10 m within a decade, whereas the current generation of ground-based detectors like the Advanced LIGO-Virgo network would take $\mathcal{O}(10_5)$ years to do so for a soft equation of state.
Equation of state at neutron-star densities and beyond from perturbative QCD
We explore the consequences of imposing robust thermodynamic constraints arising from perturbative quantum chromodynamics (QCD) when inferring the dense-matter equation-of-state (EOS). We find that the termination density, up to which the EOS modeling is performed in an inference setup, strongly affects the constraining power of the QCD input. This sensitivity in the constraining power arises from EOSs that have a specific form, with drastic softening immediately above the termination density followed by a strong stiffening. Here we also perform explicit modeling of the EOS down from perturbative-QCD densities to construct a new QCD likelihood function that incorporates additional perturbative-QCD calculations of the sound speed and is insensitive to the termination density, which we make publicly available.
Probing quarkyonic matter in neutron stars with the Bayesian nuclear-physics multimessenger astrophysics framework
The interiors of neutron stars contain matter at the highest densities realized in our Universe. Interestingly, theoretical studies of dense matter, in combination with the existence of two-solar-mass neutron stars, indicate that the speed of sound $c_s$ has to increase to values well above the conformal limit ($c_s^2$ = 1/3) before decreasing again at higher densities. Further, the decrease could be explained by either a strong first-order phase transition or a crossover transition from hadronic to quark matter. The latter scenario leads to a pronounced peak in the speed of sound, reaching values above the conformal limit, naturally explaining the inferred behavior. In this work, we use the nuclear-physics multimessenger astrophysics (NMMA) framework to compare predictions of the quarkyonic matter model with astrophysical observations of neutron stars, with the goal of constraining model parameters. Assuming quarkyonic matter to be realized within neutron stars, we find that there can be a significant amount of quarks inside the cores of neutron stars with masses in the two-solar-mass range, amounting to up to ≈0.13$M$ ⊙ , contributing ≈ 5.9% of the total mass. Furthermore, for the quarkyonic matter model investigated here, the radius of a 1.4$M$ ⊙ neutron star would be $13.44_{–1.54}^{+1.69}(13. 54_{–1.04}^{+1.02})$ km, at 95% credibility, without (with) the inclusion of AT2017gfo.
An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers
The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M ⊙ neutron star to be $R$ $=$ $11.98^{+0.35}_{-0.04}$ km.
Relativistic Hartree–Fock chiral Lagrangians with confinement, nucleon finite size and short-range effects
A relativistic Hartree–Fock Lagrangian including a chiral potential and nucleon polarisation is investigated in hopes of providing a better description of dense nuclear matter. Here, we fully consider the contribution of the exchange Fock term to the energy and the self-energies, and in addition we investigate the nucleon’s compositeness and finite size effects (confinement and form factors) and short range correlations modeled by a Jastrow ansatz. These effects are added step by step, such that their impact on the dense matter properties can be analysed in details. The parameters of the model are adjusted to reproduce fundamental properties related to the QCD theory at low energy, such as the chiral symmetry breaking, nucleon’s quark substructure and Lattice-QCD predictions, as well as two empirical properties at saturation: the binding energy and the density. All other empirical parameters, e.g., symmetry energy and its slope, incompressibility modulus, effective mass, as well as spin–isospin Landau–Midgal parameter are predictions of the models and can be used to evaluate the gain of the different approximation schemes in describing nuclear properties. Bayesian statistics is employed in order to propagate parameter uncertainties into predictions for the nuclear matter properties. We show that the splitting of the effective Landau mass is largely influenced by the value of the ρ T coupling, and we show that the fit to the symmetry energy, which induces an increase of the coupling constant g ρ by about 20–25% compared to the case where it is fixed by the quark model, provides a very good EoS compatible with the present nuclear physics knowledge.
Impact of O4 future detections on the determination of the dense matter equations of state
Here, in view of the next LIGO-Virgo-KAGRA Observing period O4 (to start in Spring 2023), we address the question of the ability of the interferometers network to discriminate among different neutron stars equation of states better than what was possible with the observation of the binary neutron stars merger GW170817. We show that the observation of an event similar to GW170817 during O4 would allow us to resolve the dimensionless effective tidal deformability $\tilde{Λ}$ within an uncertainty 7 times better than the one obtained in O2. Thanks to the expected increase in sensitivities, we show that any GW170817-like single-event within a distance of 100 Mpc would imply significantly improved constraints of the neutron stars equations of state. We also illustrate the important impact of the noise in the analysis of the signal, showing how it can impact the effective tidal deformability probability density function for large signal-to-noise ratio.
Perturbative QCD and the neutron star equation of state
Here, we construct a physics-agnostic approach to the neutron star (NS) equation of state (EoS) based on a sound speed model, which connects both low-density information from nuclear theory and high-density constraints from perturbative quantum chromodynamics (pQCD). Using this approach, we study the impact of pQCD calculations on NS EoS that have been constrained by astrophysical observations. We find that pQCD affects the EoS mainly beyond the densities realized in NS. Furthermore, we observe an interesting interplay between pQCD and astrophysical constraints, with pQCD preferring softer EoS for the heaviest NS while recent NICER observations suggest the EoS to be stiffer. We explore the sensitivity of our findings to pQCD uncertainties and study the constraining power of pQCD if future observations of heavy NS were to suggest radii larger than 13 km.
Confronting a set of Skyrme and $\chi _{EFT}$ predictions for the crust of neutron stars
With the improved accuracy of neutron star observational data, it is necessary to derive new equation of state where the crust and the core are consistently calculated within a unified approach. For this purpose we describe non-uniform matter in the crust of neutron stars employing a compressible liquid-drop model, where the bulk and the neutron fluid terms are given from the same model as the one describing uniform matter present in the core. We then generate a set of fifteen unified equations of state for cold catalyzed neutron stars built on realistic modelings of the nuclear interaction, which belongs to two main groups: the first one derives from the phenomenological Skyrme interaction and the second one from χEFT Hamiltonians. The confrontation of these model predictions allows us to investigate the model dependence for the crust properties, and in particular the effect of neutron matter at low density. Here, the new set of unified equations of state is available at the CompOSE repository.