Neutron Star Mass Research Articles

Fast Neutrino Cooling in the Accreting Neutron Star MXB 1659-29

Modeling of crust heating and cooling across multiple accretion outbursts of the low mass X-ray binary MXB 1659-29 indicates that the neutrino luminosity of the neutron star core is consistent with direct Urca (dUrca) reactions occurring in ∼1% of the core volume. We investigate this scenario with neutron star models that include a detailed equation of state parametrized by the slope of the nuclear symmetry energy L, and a range of neutron and proton superfluid gaps. We find that the predicted neutron star mass depends sensitively on L and the assumed gaps. We discuss which combinations of superfluid gaps reproduce the inferred neutrino luminosity. Larger values of L ≳ 80 MeV require superfluidity to suppress dUrca reactions in low mass neutron stars, i.e., the proton or neutron gap is sufficiently strong and extends to high enough density. However, the largest gaps give masses near the maximum mass, making it difficult to accommodate colder neutron stars. The heat capacities of our models span the range from fully paired to fully unpaired nucleons meaning that long-term observations of core cooling could distinguish between models. As a route to solutions with a larger emitting volume, which could provide a more natural explanation for the inferred neutrino luminosity, we discuss the possibility of alternative, less efficient, fast cooling processes in exotic cores. To be consistent with the inferred neutrino luminosity, such processes must be within a factor of ∼1000 of dUrca. We discuss the impact of future constraints on neutron star mass, radius, and the density dependence of the symmetry energy.

Constraints on the merging binary neutron star mass distribution and equation of state based on the incidence of jets in the population

A relativistic jet has been produced in the single well-localised binary neutron star (BNS) merger detected to date in gravitational waves (GWs), and the local rates of BNS mergers and short gamma-ray bursts are of the same order of magnitude. This suggests that jet formation is not a rare outcome for BNS mergers, and we show that this intuition can be turned into a quantitative constraint: at least about one-third of GW-detected BNS mergers and at least about one-fifth of all BNS mergers should produce a successful jet (90% credible level). Whether a jet is launched depends on the properties of the merger remnant and of the surrounding accretion disc, which in turn are a function of the progenitor binary masses and equation of state (EoS). The incidence of jets in the population therefore carries information about the binary component mass distribution and EoS. Under the assumption that a jet can only be produced by a black hole remnant surrounded by a non-negligible accretion disc, we show how the jet incidence can be used to place a joint constraint on the space of BNS component mass distributions and EoS. The result points to a broad mass distribution, with particularly strong support for masses in the 1.3 − 1.6 M⊙ range. The constraints on the EoS are shallow, but we show how they will tighten as the knowledge on the jet incidence improves. We also discuss how to extend the method to include future BNS mergers, with possibly uncertain jet associations.

GW170817 4.5 Yr After Merger: Dynamical Ejecta Afterglow Constraints

GW170817 is the first binary neutron star (NS) merger detected in gravitational waves (GWs) and photons, and so far remains the only GW event of its class with a definitive electromagnetic counterpart. Radio emission from the structured jet associated with GW170817 has faded below the sensitivity achievable via deep radio observations with the most sensitive radio arrays currently in operation. Hence, we now have the opportunity to probe the radio re-brightening that some models predict, which should emerge at late times from the interaction of the dynamically stripped merger ejecta with the interstellar medium. Here we present the latest results from our deep radio observations of the GW170817 field with the Karl G. Jansky Very Large Array (VLA), 4.5 yr after the merger. Our new data at 3 GHz do not show any compelling evidence for emission in excess to the tail of the jet afterglow (<3.3 μJy), confirming our previous results. We thus set new constraints on the dynamical ejecta afterglow models. These constraints favor single-speed ejecta with energies ≲1050 erg (for an ejecta speed of β 0 = 0.5), or steeper energy–speed distributions of the kilonova ejecta. Our results also suggest larger values of the cold, nonrotating maximum NS mass in equal-mass scenarios. However, without a detection of the dynamical ejecta afterglow, obtaining precise constraints on the NS equation of state remains challenging.

Inferring the Neutron Star Maximum Mass and Lower Mass Gap in Neutron Star–Black Hole Systems with Spin

Gravitational-wave (GW) detections of merging neutron star–black hole (NSBH) systems probe astrophysical neutron star (NS) and black hole (BH) mass distributions, especially at the transition between NS and BH masses. Of particular interest are the maximum NS mass, minimum BH mass, and potential mass gap between them. While previous GW population analyses assumed all NSs obey the same maximum mass, if rapidly spinning NSs exist, they can extend to larger maximum masses than nonspinning NSs. In fact, several authors have proposed that the ∼2.6 M ⊙ object in the event GW190814—either the most massive NS or least massive BH observed to date—is a rapidly spinning NS. We therefore infer the NSBH mass distribution jointly with the NS spin distribution, modeling the NS maximum mass as a function of spin. Using four LIGO–Virgo NSBH events including GW190814, if we assume that the NS spin distribution is uniformly distributed up to the maximum (breakup) spin, we infer the maximum nonspinning NS mass is (90% credibility), while assuming only nonspinning NSs, the NS maximum mass must be >2.53 M ⊙ (90% credibility). The data support the mass gap’s existence, with a minimum BH mass at . With future observations, under simplified assumptions, 150 NSBH events may constrain the maximum nonspinning NS mass to ±0.02 M ⊙, and we may even measure the relation between the NS spin and maximum mass entirely from GW data. If rapidly rotating NSs exist, their spins and masses must be modeled simultaneously to avoid biasing the NS maximum mass.

Impact of the equation of state on f - and p - mode oscillations of neutron stars

We investigate the impact of the neutron-star matter equation of state on the $f$- and $p_1$-mode oscillations of neutron stars obtained within the Cowling approximation and linearized general relativity. The $f$- and $p_1$-mode oscillation frequencies, and their damping times are calculated using representative sets of Skyrme Hartree-Fock and relativistic mean-field models, all of which reproduce nuclear systematics and support $2M_\odot$ neutron stars. Our study shows strong correlations between the frequencies of $f$- and $p_1$-modes and their damping times with the pressure of $\beta$-equilibrated matter at densities equal to or slightly higher than the nuclear saturation density $\rho_0$. Such correlations are found to be almost independent of the composition of the stars. The frequency of the $p_1$-mode of $1.4M_\odot$ star is strongly correlated with the slope of the symmetry energy $L_0$ and $\beta$-equilibrated pressure at density $\rho_0$. Compared to GR calculations, the error in the Cowling approximation for the $f$-mode is about 30\% for neutron stars of low mass, whereas it decreases with increasing mass. The accuracy of the $p_1$-mode is better than 15\% for neutron stars of maximum mass, and improves for lower masses and higher number of radial nodes.

Exploring the Fate of Stellar Core Collapse with Supernova Relic Neutrinos

Core collapse of massive stars leads to different fates for various physical factors, which gives different spectra of the emitted neutrinos. We focus on the supernova relic neutrinos (SRNs) as a probe to investigate the stellar collapse fate. We present the SRN fluxes and event rate spectra at a detector for three resultant states after stellar core collapse, the typical mass neutron star, the higher mass neutron star, or the failed supernova forming a black hole, based on different nuclear equations of state. Then possible SRN fluxes are formed as mixtures of the three components. We also show the expected sensitivities at the next-generation water-based Cherenkov detectors, SK-Gd and Hyper-Kamiokande, as constraining the mixture fractions. This study provides a practical example of extracting astrophysical constraints through SRN measurement.

Quasi-stationary sequences of hyper-massive neutron stars with exotic equations of state

In this work, we study the effect of differential rotation, finite temperature and strangeness on the quasi stationary sequences of hyper massive neutron stars (HMNS). We generate constant rest mass sequences of differentially rotating and uniformly rotating stars. The nucleonic matter relevant to the star interior is described within the framework of the relativistic mean field model with the DD2 parameter set. We also consider the strange $\Lambda$ hyperons using the BHB$\Lambda\phi$ equation of state (EoS). Additionally, we probe the behaviour of neutron stars (NS) with these compositions at different temperatures. We report that the addition of hyperons to the EoS produces a significant boost to the spin-up phenomenon. Moreover, increasing the temperature can make the spin-up more robust. We also study the impact of strangeness and thermal effects on the T/W instability. Finally, we analyse equilibrium sequences of a NS following a stable transition from differential rotation to uniform rotation. The decrease in frequency relative to angular momentum loss during this transition is significantly smaller for EoS containing hyperons, compared to nucleonic EoS.

Ultrawide Black Hole—Neutron Star Binaries as a Possible Source for Gravitational Waves and Short Gamma-Ray Bursts

The third observing run of the LIGO/Virgo/KARGA collaboration reported a few neutron star–black hole (NSBH) merger events. While NSBH mergers have yet to receive extensive theoretical attention, they may have a promising electromagnetic signature in the form of short gamma-ray bursts. Here we show that NSBH dynamical mergers can naturally form from ultrawide binaries in the field. Flyby gravitational interactions with other neighbors in the galaxy in these ultrawide systems may result in high eccentricity that drives the binary into a merger. We show that this process can result in a merger rate at the order of ∼10 Gpc−3 yr−1 (∼1 Gpc−3 yr−1) for elliptical (spiral) galaxies. This channel predicts a higher merger rate with a higher velocity dispersion of the host galaxy, a delay time distribution that is steeper than uniform but shallower than 1/t, and a higher merger rate for a lower black hole to neutron star mass ratio.

The role of supernova convection for the lower mass gap in the isolated binary formation of gravitational wave sources

ABSTRACT Understanding the astrophysical phenomena involving compact objects requires an insight about the engine behind the core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of the growing populations of the detected black hole–black hole, black hole–neutron star, and neutron star–neutron star mergers. The time-scale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. We adopt new formulas for the relation between the pre-SN star properties and their remnants and check how they impact the population of double compact object (DCO) mergers formed via the isolated binary evolution. The new formulas give one the ability to test a wide spectrum of assumptions on the convection growth time. In particular, the different variants allow for a smooth transition between having a deep mass gap and a remnant mass distribution filled by massive neutron stars and low-mass black holes. We present the distribution of masses, mass ratios, and the local merger rate densities of DCO for the different variants of new formulas and test them together with different approaches to other highly uncertain processes. We find that the mass distribution of DCO mergers is sensitive to the adopted assumption on the SN convection growth time-scale up to $m_1+m_2 \lesssim 35 \, \mathrm{M}_{\odot }$. Between the two extreme tested variants the probability of compact object formation within the mass gap may differ by up to approximately two orders of magnitude.

Kaon–baryon coupling schemes and kaon condensation in hyperon-mixed matter

Abstract The possible coexistence of kaon condensation and hyperons in highly dense matter [the (Y + K) phase] is investigated on the basis of the relativistic mean-field theory combined with the effective chiral Lagrangian. Two coupling schemes for the s-wave kaon–baryon interaction are compared regarding the onset density of kaon condensation in hyperon-mixed matter and the equation of state for the developed (Y + K) phase. One is the contact interaction scheme related to the nonlinear effective chiral Lagrangian. The other is the meson exchange scheme, where the interaction vertices between the kaon field and baryons are described by an exchange of mesons (σ, σ* mesons for scalar coupling, and ω, ρ, ϕ mesons for vector coupling). It is shown that in the meson exchange scheme, the contribution from the nonlinear scalar self-interaction gives rise to a repulsive effect for the kaon effective energy, pushing up the onset density of kaon condensation as compared with the contact interaction scheme. In general, the difference in kaon–baryon dynamics between the contact interaction scheme and the meson exchange scheme relies on the specific forms of the nonlinear self-interacting meson terms. They generate many-baryon forces through the equations of motion for the meson mean fields. However, they should have a definite effect on the ground state properties of nuclear matter only around the saturation density. It is shown that the nonlinear self-interacting term is not relevant to repulsive energy leading to stiffening of the equation of state at high densities, and that it cannot be compensated with a large attractive energy due to the appearance of the (Y + K) phase in the case of the contact interaction scheme. We also discuss what effects are necessary in the contact interaction scheme to make the equation of state with (Y + K) phase stiff enough to be consistent with recent observations of massive neutron stars.

Nearly model-independent constraints on dense matter equation of state in a Bayesian approach

We apply a 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 (${\mathrm{N}}^{3}\mathrm{LO}$) calculation in chiral effective field theory. The tidal deformability and radius of a NS with mass $1--2\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ are found to be strongly correlated with the pressure of $\ensuremath{\beta}$-equilibrated matter at densities higher than the saturation density (${\ensuremath{\rho}}_{0}=0.16\text{ }\text{ }{\mathrm{fm}}^{\ensuremath{-}3}$) in a nearly model-independent manner. These correlations are employed to parametrize the pressure for $\ensuremath{\beta}$-equilibrated matter, around $2{\ensuremath{\rho}}_{0}$, as a function of neutron star mass and the corresponding tidal deformability. The maximum mass of the neutron star is also found to be strongly correlated with the pressure of $\ensuremath{\beta}$-equilibrated matter at densities $\ensuremath{\sim}4.5{\ensuremath{\rho}}_{0}$.

Relativistic mean field model parametrizations in the light of GW170817, GW190814, and PSR J0740+6620

Three parametrizations DOPS1, DOPS2, and DOPS3 (named after the Department of Physics Shimla) of the relativistic mean field model have been proposed with the inclusion of all possible self and mixed interactions between the scalar-isoscalar ($\ensuremath{\sigma}$), vector-isoscalar ($\ensuremath{\omega}$), and vector-isovector ($\ensuremath{\rho}$) mesons up to quartic order. The generated parameter sets are in harmony with the finite and bulk nuclear matter properties. A set of equations of state (EOSs) composed of pure hadronic (nucleonic) matter and nucleonic with quark matter (hybrid EOSs) for superdense hadron-quark matter in $\ensuremath{\beta}$ equilibrium is obtained. The quark matter phase is calculated by using the three-flavor Nambu-Jona-Lasinio (NJL) model. The maximum mass of a nonrotating neutron star with DOPS1 parametrization is found to be around $2.6M\ensuremath{\bigodot}$ for the pure nucleonic matter, which satisfies the recent gravitational wave analysis of GW190814 [Abbott et al., Astrophys. J. Lett. 896, L44 (2020)] with possible maximum mass constraint indicating that the secondary component of GW190814 could be a nonrotating heaviest neutron star composed of pure nucleonic matter. EOSs computed with the DOPS2 and DOPS3 parametrizations satisfy the x-ray observational data [Steiner et al., Astrophys. J. 722, 33 (2010)] and the recent observations of GW170817 maximum mass constraint of a stable nonrotating neutron star in the range $2.01\ifmmode\pm\else\textpm\fi{}0.04--2.16\ifmmode\pm\else\textpm\fi{}0.03M\ensuremath{\bigodot}$ [Rezzolla et al., Astrophys. J. Lett. 852, L25 (2018)] and also in good agreement with constraints on mass and radius measurement for PSR $\mathrm{J}0740+6620$ (NICER) [Riley et al., Astrophys. J. Lett. 918, L27 (2021); Miller et al., Astrophys. J. Lett. 918, L28 (2021)]. The hybrid EOSs obtained with the NJL model also satisfy astrophysical constraints on the maximum mass of a neutron star from PSR J1614-2230 [Demorest et al., Nature (London) 467, 1081 (2010)]. We also present the results for dimensionless tidal deformability, $\mathrm{\ensuremath{\Lambda}}$ which are consistent with the waveform models analysis of GW170817.

MAXI J1957+032: a new accreting millisecond X-ray pulsar in an ultra-compact binary

ABSTRACT The detection of coherent X-ray pulsations at ∼314 Hz (3.2 ms) classifies MAXI J1957+032 as a fast-rotating, accreting neutron star. We present the temporal and spectral analysis performed using NICER observations collected during the latest outburst of the source. Doppler modulation of the X-ray pulsation revealed the ultra-compact nature of the binary system characterized by an orbital period of ∼1 h and a projected semimajor axis of 14 lt-ms. The neutron star binary mass function suggests a minimum donor mass of 1.7 × 10−2 M⊙, assuming a neutron star mass of 1.4 M⊙ and a binary inclination angle lower than 60 deg. This assumption is supported by the lack of eclipses or dips in the X-ray light curve of the source. We characterized the 0.5–10 keV energy spectrum of the source in outburst as the superposition of a relatively cold black-body-like thermal emission compatible with the emission from the neutron star surface and a Comptonization component with photon index consistent with a typical hard state. We did not find evidence for iron K α lines or reflection components.

Unified neutron star EOSs and neutron star structures in RMF models

In the framework of the Thomas-Fermi approximation, we systematically study the EOSs and microscopic structures of neutron star matter in a vast density range with n b ≈ 10−10-2 fm−3, where various covariant density functionals are adopted, i.e., those with nonlinear self couplings (NL3, PK1, TM1, GM1, MTVTC) and density-dependent couplings (DD-LZ1, DDME-X, PKDD, DD-ME2, DD2, TW99). It is found that the EOSs generally coincide with each other at n b ≲ 10−4 fm−3 and 0.1 fm−3 ≲ n b ≲ 0.3 fm−3, while in other density regions they are sensitive to the effective interactions between nucleons. By adopting functionals with a larger slope of symmetry energy L, the curvature parameter K sym and neutron drip density generally increases, while the droplet size, proton number of nucleus, core-crust transition density, and onset density of non-spherical nuclei, decrease. All functionals predict neutron stars with maximum masses exceeding the two-solar-mass limit, while those of DD2, DD-LZ1, DD-ME2, and DDME-X predict optimum neutron star radii according to the observational constraints. Nevertheless, the corresponding skewness coefficients J are much larger than expected, while only the functionals MTVTC and TW99 meet the start-of-art constraints on J. More accurate measurements on the radius of PSR J0740 + 6620 and the maximum mass of neutron stars are thus essential to identify the functional that satisfies all constraints from nuclear physics and astrophysical observations. Approximate linear correlations between neutron stars’ radii at M = 1.4M ⊙ and 2M ⊙, the slope L and curvature parameter K sym of symmetry energy are observed as well, which are mainly attributed to the curvature-slope correlations in the functionals adopted here. The results presented here are applicable for investigations of the structures and evolutions of compact stars in a unified manner.

Constraining neutron star properties with a new equation of state insensitive approach

Instead of parameterizing the pressure-density relation of a neutron star (NS), one can parameterize its macroscopic properties such as mass ($M$), radius ($R$), and dimensionless tidal deformability ($\Lambda$) to infer the equation of state (EoS) combining electromagnetic and gravitational wave (GW) observations. We present a new method to parameterize $R(M)$ and $\Lambda(M)$ relations, which approximate the candidate EoSs with accuracy better than 5\% for all masses and span a broad region of $M-R-\Lambda$ plane. Using this method we combine the $M-\Lambda$ measurement from GW170817 and GW190425, and simultaneous $M-R$ measurement of PSR J0030+0451 and PSR J0740+6620 to place joint constraints on NS properties. At 90 \% confidence, we infer $R_{1.4}=12.05_{-0.87}^{+0.98}$ km and $\Lambda_{1.4}=372_{-150}^{+220}$ for a $1.4 M_{\odot}$ NS, and $R_{2.08}=12.65_{-1.46}^{+1.36}$ km for a $2.08 M_{\odot}$ NS. Furthermore, we use the inferred values of the maximum mass of a nonrotating NS $M_{\rm max}=2.52_{-0.29}^{+0.33} M_{\odot}$ to investigate the nature of the secondary objects in three potential neutron star-black hole merger (NSBH) system.

Solving Tolman–Oppenheimer–Volkoff equations in f(T) gravity: A novel approach applied to some realistic equations of state

There are many ways to probe alternative theories of gravity, namely, via experimental tests at solar system scale, cosmological data and models, gravitational waves and compact objects. In this paper, we consider a model of gravity with torsion [Formula: see text] applied to compact objects such as neutron stars (NSs) for a couple of realistic equations of state (EOS). To do so, we follow our previous articles, in which we show how to model compact stars in [Formula: see text] gravity by obtaining its corresponding Tolman–Oppenheimer–Volkof equations. In this modeling of NS in [Formula: see text] gravity presented here, we calculate, among other things, the maximum mass allowed for a given realistic EOS, which would also allow us to evaluate which models are in accordance with the observations. The results already known to General Relativity must be reproduced to some extent and, eventually, we can find models that allow higher maximum masses for NSs than Relativity itself, which could explain, for example, the secondary component of the event GW190814, if this star is a massive NS.

The Hadron-quark Crossover in Neutron Star within Gaussian Process Regression Method

The equations of state of the neutron star at the hadron-quark crossover region are interpolated with the Gaussian process regression (GPR) method, which can reduce the randomness of present interpolation schemes. The relativistic mean-field (RMF) model and Nambu–Jona-Lasinio (NJL) model are employed to describe the hadronic phase and quark phase, respectively. In the RMF model, the coupling term between ω and ρ mesons is considered to control the density-dependent behaviors of symmetry energy, i.e., the slope of symmetry energy L. Furthermore, the vector interaction between quarks is included in the NJL model to obtain the additional repulsive contributions. Their coupling strengths and the crossover windows are discussed in the present framework under the constraints on the neutron star from gravitational-wave detections, massive neutron star measurements, mass–radius simultaneous observation of the NICER Collaboration, and the neutron skin thickness of 208Pb from PREX-II. It is found that the slope of symmetry energy, L, should be around 50−90 MeV and the crossover window is (0.3, 0.6) fm−3 with these observables. Furthermore, the uncertainties of neutron star masses and radii in the hadron-quark crossover regions are also predicted by the GPR method.

Neural network reconstruction of the dense matter equation of state from neutron star observables

The Equation of State (EoS) of strongly interacting cold and hot ultra-dense QCD matter remains a major challenge in the field of nuclear astrophysics. With the advancements in measurements of neutron star masses, radii, and tidal deformabilities, from electromagnetic and gravitational wave observations, neutron stars play an important role in constraining the ultra-dense QCD matter EoS. In this work, we present a novel method that exploits deep learning techniques to reconstruct the neutron star EoS from mass-radius (M-R) observations. We employ neural networks (NNs) to represent the EoS in a model-independent way, within the range of ∼1-7 times the nuclear saturation density. The unsupervised Automatic Differentiation (AD) framework is implemented to optimize the EoS, so as to yield through TOV equations, an M-R curve that best fits the observations.We demonstrate that this method works by rebuilding the EoS on mock data, i.e., mass-radius pairs derived from a randomly generated polytropic EoS. The reconstructed EoS fits the mock data with reasonable accuracy, using just 11 mock M-R pairs observations, close to the current number of actual observations.

Upgraded GMRT survey for pulsars in globular clusters

Context. Globular clusters (GCs) contain a unique pulsar population, with many exotic systems that can form only in their dense stellar environments. Such systems are potentially very interesting for new tests of gravity theories and neutron-star mass measurements. Aims. The leap in sensitivity of the upgraded Giant Metrewave Radio Telescope (uGMRT) in India, especially at low radio frequencies (&lt; 1 GHz), motivated a new search for radio pulsars in a group of eight southern GCs. We aim to image these clusters in order to have independent measurements of the radio fluxes of known pulsars and the identification of bright radio sources that could be pulsars missed by pulsation search pipelines due to their inherent limitations. Methods. The observations were conducted at 650 MHz (Band 4 receivers) on Terzan 5, NGC 6441, NGC 6440, and NGC 6544, and at 400 MHz (Band 3 receivers) on NGC 6652, NGC 6539, NGC 1851, and M 30. Segmented acceleration and jerk searches were performed on the data. Simultaneously, we obtained interferometric data on these clusters, which were later converted into radio images. Results. We discovered PSR J1835−3259B, a 1.83-ms pulsar in NGC 6652; this is in a near-circular wide orbit of 28.7-h with an unidentified low-mass (∼0.2 M⊙) companion, likely a helium white dwarf. We derived a ten-year timing solution for this system. We also present measurements of scattering, flux densities, and spectral indices for some of the previously known pulsars in these GCs. Conclusions. A significant fraction of the pulsars in these clusters have steep spectral indices. Additionally, we detected eight radio point sources not associated with any known pulsar positions in the radio images. There are four newly identified sources, three in NGC 6652 and one in NGC 6539, as well as one previously identified source in NGC 1851, NGC 6440, NGC 6544, and Terzan 5. Surprisingly, our images show that our newly discovered pulsar, PSR J1835−3259B, is the brightest pulsar in all GCs we have imaged; like other pulsars with broad profiles (Terzan 5 C and O), its flux density in the radio images is much larger than in its pulsations. This indicates that their pulsed emission is only a fraction of their total emission. The detection of radio sources outside the core radii but well within the tidal radii of these clusters show that future GC surveys should complement the search analysis by using the imaging capability of interferometers, and preferentially synthesise large number of search beams in order to obtain a larger field of view.

Thermodynamics of Hot Neutron Stars and Universal Relations

Over the last few years, the detection of gravitational waves from binary neutron star systems has rekindled our hopes for a deeper understanding of the unknown nature of ultradense matter. In particular, gravitational wave constraints on the tidal deformability of a neutron star can be translated into constraints on several neutron star properties using a set of universal relations. Apart from binary neutron star mergers, supernova explosions are also important candidates for the detection of multimessenger signals. Such observations may allow us to impose significant constraints on the binding energy of neutron stars. The purpose of the present study is twofold. Firstly, we investigate the agreement of finite temperature equations of state with established universal relations. Secondly, we examine the possible existence of a universal relation between the binding energy and the dimensionless tidal deformability, which are the bulk properties connected to the most promising sources for multimessenger signals. We find that hot equations of state are not always compatible with accepted universal relations. Therefore, the use of such expressions for probing general relativity or imposing constraints on the structure of neutron stars would be inconclusive (when thermal effects are present). Additionally, we show that the binding energy and the dimensionless tidal deformability exhibit a universal trend at least for moderate neutron star masses. The latter allows us to set bounds on the binding energy of a 1.4 M⊙ neutron star using data from the GW170817 event. Finally, we provide a relation between the compactness, the binding energy and the dimensionless tidal deformability of a neutron star that is accurate for cold and hot isentropic equations of state.