Galactic Neutron Stars Research Articles

High Frequency Gravitational Wave bounds from galactic neutron stars

High-Frequency Gravitational Waves (HFGWs) constitute a unique window on the early Universe as well as exotic astrophysical objects. While the current gravitational wave experiments are more dedicated to the low frequency regime, the graviton conversion into photons in a strong magnetic field constitutes a powerful tool to probe HFGWs. In this paper, we show that neutron stars, due to their extreme magnetic field, are a perfect laboratory to study the conversion of HFGWs into photons. Using realistic models for the galactic neutron star population, we calculate for the first time the expected photon flux induced by the conversion of an isotropic stochastic gravitational wave background in the magnetosphere of the ensemble of neutron stars present in the Milky Way. We compare this photon flux to the observed one from several telescopes and derive upper limits on the stochastic gravitational wave background in the frequency range 108 Hz–1025 Hz. We find our limits to be competitive in the frequency range 108 Hz–1012 Hz.

Do Neutron Star Ultraluminous X-Ray Sources Masquerade as Intermediate-mass Black Holes in Radio and X-Ray?

Ultraluminous X-ray sources (ULXs) were once largely believed to be powered by super-Eddington accretion onto stellar-mass black holes, although in some rare cases, ULXs also serve as potential candidates for (sub-Eddington) intermediate-mass black holes. However, a total of eight ULXs have now been confirmed to be powered by neutron stars, thanks to observed pulsations, and may act as contaminants for the radio/X-ray selection of intermediate-mass black holes. Here, we present the first comprehensive radio study of seven known neutron star ULXs using new and archival data from the Karl G. Jansky Very Large Array and the Australia Telescope Compact Array, combined with the literature. Across this sample, there is only one confident radio detection, from the Galactic neutron star ULX Swift J0243.6+6124. The other six objects in our sample are extragalactic, and only one has coincident radio emission, which we conclude is most likely contamination from a background H ii region. We conclude that with current facilities, neutron star ULXs do not produce significant enough radio emission to cause them to be misidentified as radio-/X-ray-selected intermediate-mass black hole candidates. Thus, if background star formation has been properly considered, the current study indicates that a ULX with a compact radio counterpart is not likely to be a neutron star.

Kinematic constraints on the ages and kick velocities of Galactic neutron star binaries

The systems creating binary neutron stars (BNSs) experience systemic kicks when one of the components goes supernova. The combined magnitude of these kicks is still a topic of debate and has implications for the eventual location of the transient resulting from the merger of the binary. For example, the offsets of short-duration gamma-ray bursts resulting from BNS mergers depend on BNS kicks. We investigated Galactic BNSs and traced their motion through the Galaxy. This enabled us to estimate their kinematic ages and construct a BNS kick distribution based on their Galactic trajectories. We used the pulsar periods and their derivatives to estimate the characteristic spin-down ages of the binaries. Moreover, we used a Monte Carlo estimation of their present-day velocity vector in order to trace back their trajectory and estimate their kinematic ages. These trajectories, in turn, were used to determine the eccentricity of their Galactic orbit. Based on simulations of kicked objects in the Galactic potential, we investigated the relationship between this eccentricity and the kick velocity in order to constrain the kicks imparted to the binaries at birth. We find that the Galactic BNSs are likely older than $ Myr, which means their current (scalar) galactocentric speeds are not representative of their initial kicks. However, we find a close relationship between the eccentricity of a Galactic trajectory and the experienced kick. Using this relation, we constrained the kicks of the Galactic BNSs, depending on the kind of isotropy assumed in estimating their velocity vectors. These kick velocities are well described by a log-normal distribution peaking around $ km/s and coincide with the peculiar velocities of the binaries at their last disc crossing. We conclude that BNSs receive kicks following a distribution that peaks at kick velocities lower than found in isolated pulsars. However, we find no tension between this distribution and literature on short-duration gamma-ray burst offsets.

Deceleration of kicked objects due to the Galactic potential

Context. Various stellar objects experience a velocity kick at some point in their evolution. These include neutron stars and black holes at their birth, or binary systems when one of the two components goes supernova. For most of these objects, the magnitude of the kick and its impact on the object dynamics remains a topic of debate. Aims. We investigate how kicks alter the velocity distribution of objects born in the Milky Way disc, both immediately after the kick and at later times, and whether these kicks are encoded in the observed population of Galactic neutron stars. Methods. We simulated the Galactic trajectories of point masses on circular orbits in the disc after being perturbed by an isotropic kick, with a Maxwellian distribution where σ = 265 km s−1. Then, we simulated the motion of these point masses for 200 Myr. These trajectories were then evaluated, either for the Milky Way population as a whole or for those passing within two kiloparsecs of the Sun, to get the time evolution of the velocities. Results. During the first 20 Myr, the bulk velocity of kicked objects becomes temporarily aligned with the cylindrical radius, implying an anisotropy in the velocity orientations. Beyond this age, the velocity distribution shifts towards lower values and settles to a median of ∼200 km s−1. Around the Sun, the distribution also loses its upper tail, primarily due to unbound objects escaping the Galaxy. We compared this to the velocities of Galactic pulsars and find that pulsars show a similar evolution with characteristic age. Conclusions. The shift in the velocity distribution is due to bound objects spending most of their orbits at larger radii after the kick. They are, therefore, decelerated by the Galactic potential. We find the same deceleration for nearby objects and the total population, and conclude that it is also observed in Galactic pulsars. Because of this effect, the (scalar) speeds of old neutron stars provide little information about their kicks at birth.

Constraining mass transfer models with galactic neutron star–white dwarf binaries as gravitational wave sources

ABSTRACT Neutron star–white dwarf (NSWD) binaries are one of the most abundant sources of gravitational waves (GWs) in the Milky Way. These GW sources are the evolutionary products of primordial binaries that experienced many processes of binary interaction. We employ a binary population synthesis method to investigate the properties of Galactic NSWD binaries detectable by the Laser Interferometer Space Antenna (LISA). In this paper, only the NSWD systems with a COWD or ONeWD component are included. We consider various models related to mass-transfer efficiencies during primordial binary evolution, supernova explosion mechanisms at NS formation, common envelope ejection efficiencies, and critical WD masses that determine the stability of mass transfer between WDs and NSs. Based on our calculations, we estimate that tens to hundreds of LISA NSWD binaries exist in the Milky Way. We find that the detection of LISA NSWD binaries is able to provide profound insights into mass-transfer efficiencies during the evolution of primordial binaries and critical WD masses during mass transfer from a WD to an NS.

Thallium-208: A Beacon of InSitu Neutron Capture Nucleosynthesis.

We demonstrate that the well-known 2.6MeV gamma-ray emission line from thallium-208 could serve as a real-time indicator of astrophysical heavy element production, with both rapid (r) and intermediate (i) neutron capture processes capable of its synthesis. We consider the r process in a Galactic neutron star merger and show Tl-208 to be detectable from ∼12 hours to ∼ten days, and again ∼1-20 years postevent. Detection of Tl-208 represents the only identified prospect for a direct signal of lead production (implying gold synthesis), arguing for the importance of future MeV telescope missions which aim to detect Galactic events but may also be able to reach some nearby galaxies in the Local Group.

Mass Distribution and Maximum Mass of Neutron Stars: Effects of Orbital Inclination Angle

Matter at ultra-high densities finds a physical realization inside neutron stars. One key property is their maximum mass, which has far-reaching implications for astrophysics and the equation of state of ultra dense matter. In this work, we employ Bayesian analysis to scrutinize the mass distribution and maximum mass threshold of galactic neutron stars. We compare two distinct models to assess the impact of assuming a uniform distribution for the most important quantity, the cosine of orbital inclination angles (i), which has been a common practice in previous analyses. This prevailing assumption yields a maximum mass of 2.25 M⊙ (2.15–3.32 M⊙ within 90% confidence), with a strong peak around the maximum value. However, in the second model, which indirectly includes observational constraints of i, the analysis supports a mass limit of 2.56−0.58+0.87M⊙ (2σ uncertainty), a result that points in the same direction as some recent results gathered from gravitational wave observations, although their statistics are still limited. This work stresses the importance of an accurate treatment of orbital inclination angles, and contributes to the ongoing debate about the maximum neutron star mass, further emphasizing the critical role of uncertainties in the individual neutron star mass determinations.

The Galactic neutron star population – II. Systemic velocities and merger locations of binary neutron stars

ABSTRACT The merger locations of binary neutron stars (BNSs) encode their galactic kinematics and provide insights into their connection to short gamma-ray bursts (SGRBs). In this work, we use the sample of Galactic BNSs with measured proper motions to investigate their kinematics and predict their merger locations. Using a synthetic image of the Milky Way and its Galactic potential we analyse the BNS mergers as seen from an extragalactic viewpoint and compare them to the location of SGRBs on and around their host galaxies. We find that the Galactocentric transverse velocities of the BNSs are similar in magnitude and direction to those of their Local Standards of Rest, which implies that the present-day systemic velocities are not isotropically oriented and the peculiar velocities might be as low as those of BNS progenitors. Both systemic and peculiar velocities fit a lognormal distribution, with the peculiar velocities being as low as ∼22–157 km s−1. We also find that the observed BNS sample is not representative of the whole Galactic population, but rather of systems born around the Sun’s location with small peculiar velocities. When comparing the predicted BNS merger locations to SGRBs, we find that they cover the same range of projected offsets, host-normalized offsets, and fractional light. Therefore, the spread in SGRB locations can be reproduced by mergers of BNSs born in the Galactic disc with small peculiar velocities, although the median offset match is likely a coincidence due to the biased BNS sample.

Continuous Gravitational Waves from Galactic Neutron Stars: Demography, Detectability, and Prospects

Abstract We study the prospects for the detection of continuous gravitational signals from normal Galactic neutron stars, i.e., nonrecycled stars. We use a synthetic population generated by evolving stellar remnants in time, according to several models. We consider the most recent constraints set by all-sky searches for continuous gravitational waves and use them for our detectability criteria. We discuss the detection prospects for the current and the next generation of gravitational-wave detectors. We find that neutron stars whose ellipticity is solely caused by magnetic deformations cannot produce any detectable signal, not even by third-generation detectors. The currently detectable sources all have B ≲ 1012 G and deformations that are not solely due to the magnetic field. For these, we find in fact that the larger the magnetic field, the higher the ellipticity required for the signal to be detectable, and this ellipticity is well above the value induced by the magnetic field. Third-generation detectors such as the Einstein Telescope and Cosmic Explorer will be able to detect up to ≈250 more sources than current detectors. We briefly treat the case of recycled neutron stars with a simplified model. We find that continuous gravitational waves from these objects will likely remain elusive to detection by current detectors, but should be detectable with the next generation of detectors.

Deep Einstein@Home All-sky Search for Continuous Gravitational Waves in LIGO O3 Public Data

We present the results of an all-sky search for continuous gravitational waves in the public LIGO O3 data. The search covers signal frequencies 20.0 Hz ≤ f ≤ 800.0 Hz and a spin-down range down to −2.6 × 10−9 Hz s−1, motivated by detectability studies on synthetic populations of Galactic neutron stars. This search is the most sensitive all-sky search to date in this frequency/spin-down region. The initial search was performed using the first half of the public LIGO O3 data (O3a), utilizing graphical processing units provided in equal parts by the volunteers of the Einstein@Home computing project and by the ATLAS cluster. After a hierarchical follow-up in seven stages, 12 candidates remain. Six are discarded at the eighth stage, by using the remaining O3 LIGO data (O3b). The surviving six can be ascribed to continuous-wave fake signals present in the LIGO data for validation purposes. We recover these fake signals with very high accuracy with our last stage search, which coherently combines all O3 data. Based on our results, we set upper limits on the gravitational-wave amplitude h 0 and translate these into upper limits on the neutron star ellipticity and on the r-mode amplitude. The most stringent upper limits are at 203 Hz, with h 0 = 8.1 × 10−26 at the 90% confidence level. Our results exclude isolated neutron stars rotating faster than 5 ms with ellipticities greater than within a distance d from Earth and r-mode amplitudes for neutron stars spinning faster than 150 Hz.

An Overview of Compact Star Populations and Some of Its Open Problems

The study of compact object populations has come a long way since the determination of the mass of the Hulse–Taylor pulsar, and we now count on more than 150 known Galactic neutron stars and black hole masses, as well as another 180 objects from binary mergers detected from gravitational-waves by the Ligo–Virgo–KAGRA Collaboration. With a growing understanding of the variety of systems that host these objects, their formation, evolution and frequency, we are now in a position to evaluate the statistical nature of these populations, their properties, parameter correlations and long-standing problems, such as the maximum mass of neutron stars and the black hole lower mass gap, to a reasonable level of statistical significance. Here, we give an overview of the evolution and current state of the field and point to some of its standing issues. We focus on Galactic black holes, and offer an updated catalog of 35 black hole masses and orbital parameters, as well as a standardized procedure for dealing with uncertainties.

XMM-Newton and SRG/eROSITA observations of the isolated neutron star candidate 4XMM J022141.5−735632

We report the results of follow-up investigations of a possible new thermally emitting isolated neutron star (INS), 4XMM J022141.5−735632, using observations from XMM-Newton and Spectrum Roentgen Gamma (SRG) eROSITA. The analysis is complemented by Legacy Survey imaging in the optical and near-infrared wavelengths. The X-ray source, the first to be targeted by XMM-Newton in an effort to identify new INS candidates from the fourth generation of the XMM-Newton serendipitous source catalogue Data Release 9 (4XMM-DR9), shows a remarkably soft energy distribution and a lack of catalogued counterparts; the very high X-ray-to-optical flux ratio virtually excludes any other identification than an INS. Within current observational limits, no significant flux variation or change of spectral state is registered over nearly ten years. Future dedicated observations, particularly to search for pulsations, are crucial to shed further light on the nature of the X-ray source and relations to other Galactic neutron stars.

A Robust Test of the Existence of Primordial Black Holes in Galactic Dark Matter Halos

If very low mass primordial black holes (PBH) within the asteroid/moon-mass range indeed reside in galactic dark matter halos, they must necessarily collide with galactic neutron stars (NSs). These collisions must, again necessarily, form light black holes (LBHs) with masses of typical NSs, M LBH ≈ 1–2 M ⊙. LBHs may be behind events already detected by ground-based gravitational-wave detectors (GW170817, GW190425, and others such as a mixed stellar black hole–NS-mass event GW191219_163120), and most recently by microlensing (OGLE-BLG-2011-0462). Although the status of these observations as containing LBHs is not confirmed, there is no question that gravitational-wave detectors and microlensing are in principle and in practice capable of detecting LBHs. We have calculated the creation rate of LBHs resulting from these light primordial black hole (PBH) collisions with NSs. On this basis, we claim that if improved gravitational-wave detectors and microlensing statistics of the LBH events would indicate that the number of LBHs is significantly lower that what follows from the calculated creation rate, then this would be an unambiguous proof that there is no significant light PBH contribution to the galactic dark matter halos. Otherwise, if observed and calculated numbers of LBHs roughly agree, then the hypothesis of primordial black hole existence gets strong observational support, and in addition their collisions with NSs may be considered a natural creation channel for the LBHs, solving the problem of their origin, as it is known that they cannot be a product of standard stellar evolution.

Burst timescales and luminosities as links between young pulsars and fast radio bursts

Fast radio bursts (FRBs) are extragalactic radio flashes of unknown physical origin. Their high luminosities and short durations require extreme energy densities, such as those found in the vicinity of neutron stars and black holes. Studying the burst intensities and polarimetric properties on a wide range of timescales, from milliseconds down to nanoseconds, is key to understanding the emission mechanism. However, high-time-resolution studies of FRBs are limited by their unpredictable activity levels, available instrumentation and temporal broadening in the intervening ionized medium. Here we show that the repeating FRB 20200120E can produce isolated shots of emission as short as about 60 nanoseconds in duration, with brightness temperatures as high as 3 × 1041 K (excluding relativistic effects), comparable with ‘nano-shots’ from the Crab pulsar. Comparing both the range of timescales and luminosities, we find that FRB 20200120E observationally bridges the gap between known Galactic young pulsars and magnetars and the much more distant extragalactic FRBs. This suggests a common magnetically powered emission mechanism spanning many orders of magnitude in timescale and luminosity. In this Article, we probe a relatively unexplored region of the short-duration transient phase space; we highlight that there probably exists a population of ultrafast radio transients at nanosecond to microsecond timescales, which current FRB searches are insensitive to. The range of timescales and luminosities measured from the nearby fast radio burst FRB 20200120E observationally connects these extreme extragalactic transients with studies of Galactic neutron stars.

Characterizing the Fast Radio Burst Host Galaxy Population and its Connection to Transients in the Local and Extragalactic Universe

Abstract We present the localization and host galaxies of one repeating and two apparently nonrepeating fast radio bursts (FRBs). FRB 20180301A was detected and localized with the Karl G. Jansky Very Large Array to a star-forming galaxy at z = 0.3304. FRB20191228A and FRB20200906A were detected and localized by the Australian Square Kilometre Array Pathfinder to host galaxies at z = 0.2430 and z = 0.3688, respectively. We combine these with 13 other well-localized FRBs in the literature, and analyze the host galaxy properties. We find no significant differences in the host properties of repeating and apparently nonrepeating FRBs. FRB hosts are moderately star forming, with masses slightly offset from the star-forming main sequence. Star formation and low-ionization nuclear emission-line region emission are major sources of ionization in FRB host galaxies, with the former dominant in repeating FRB hosts. FRB hosts do not track stellar mass and star formation as seen in field galaxies (more than 95% confidence). FRBs are rare in massive red galaxies, suggesting that progenitor formation channels are not solely dominated by delayed channels which lag star formation by gigayears. The global properties of FRB hosts are indistinguishable from core-collapse supernovae and short gamma-ray bursts hosts, and the spatial offset (from galaxy centers) of FRBs is mostly inconsistent with that of the Galactic neutron star population (95% confidence). The spatial offsets of FRBs (normalized to the galaxy effective radius) also differ from those of globular clusters in late- and early-type galaxies with 95% confidence.

Population Properties of Neutron Stars in the Coalescing Compact Binaries

Abstract We perform a hierarchical Bayesian inference to investigate the population properties of the coalescing compact binaries involving at least one neutron star (NS). With the current gravitational-wave (GW) observation data, we can rule out none of the double Gaussian, single Gaussian, and uniform NS mass distribution models, though a specific double Gaussian model inferred from the Galactic NSs is found to be slightly more preferred. The mass distribution of black holes (BHs) in the neutron star–black hole (NSBH) population is found to be similar to that in the Galactic X-ray binaries. Additionally, the ratio of the merger rate densities between NSBHs and BNSs is estimated to be ∼3:7. The spin properties of the binaries, though constrained relatively poorly, play a nontrivial role in reconstructing the mass distribution of NSs and BHs. We find that a perfectly aligned spin distribution can be ruled out, while a purely isotropic distribution of spin orientation is still allowed. To evaluate the feasibility of reliably determining the population properties of NSs in the coalescing compact binaries with upcoming GW observations, we perform simulations with a mock population. We find that with 100 detections (including BNSs and NSBHs) the mass distribution of NSs can be well determined, and the fraction of BNSs can also be accurately estimated.

Modeling the Galactic Neutron Star Population for Use in Continuous Gravitational-wave Searches

Abstract Searches for continuous gravitational waves (GW) from unknown Galactic neutron stars provide limits on the shapes of neutron stars. A rotating neutron star will produce GW if asymmetric deformations exist in its structure that are characterized by the star’s ellipticity. In this study, we use a simple model of the spatial and spin distribution of Galactic neutron stars to estimate the total number of neutron stars probed, using GW, to a given upper limit on the ellipticity. This may help optimize future searches with improved sensitivity. The improved sensitivity of third-generation GW detectors may increase the number of neutron stars probed, to a given ellipticity, by factors of 100 to 1000.

The Galactic neutron star population – I. An extragalactic view of the Milky Way and the implications for fast radio bursts

ABSTRACT A key tool astronomers have to investigate the nature of extragalactic transients is their position on their host galaxies. Galactocentric offsets, enclosed fluxes, and the fraction of light statistic are widely used at different wavelengths to help infer the nature of transient progenitors. Motivated by the proposed link between magnetars and fast radio bursts (FRBs), we create a face-on image of the Milky Way using best estimates of its size, structure, and colour. We place Galactic magnetars, pulsars, low-mass, and high-mass X-ray binaries on this image, using the available distance information. Galactocentric offsets, enclosed fluxes, and fraction of light distributions for these systems are compared to extragalactic transient samples. We find that FRBs follow the distributions for Galactic neutron stars closest, with 24 (75 per cent) of the Anderson–Darling tests we perform having a p-value greater than 0.05. This suggests that FRBs are located on their hosts in a manner consistent with Galactic neutron stars on the Milky Way’s light, although we cannot determine which specific neutron star population is the best match. The Galactic distributions are consistent with other extragalactic transients much less often across the range of comparisons made, with type Ia SNe in second place, at only 33 per cent of tests exceeding 0.05. Overall, our results provide further support for FRB models invoking isolated young neutron stars, or binaries containing a neutron star.

Detectability of continuous gravitational waves from isolated neutron stars in the Milky Way

Aims. We estimate the number of pulsars, detectable as continuous gravitational wave sources with the current and future gravitational-wave detectors, assuming a simple phenomenological model of evolving non-axisymmetry of the rotating neutron star. Methods. We employed a numerical model of the Galactic neutron star population, with the properties established by comparison with radio observations of isolated Galactic pulsars. We generated an arbitrarily large synthetic population of neutron stars and evolved their period, magnetic field, and position in space. We used a gravitational wave emission model based on exponentially decaying ellipticity (i.e. non-axisymmetry of the star) with no assumption of the origin of a given ellipticity. We calculated the expected signal in a given detector for a one-year observation, and assumed a detection criterion of the signal-to-noise ratio of 11.4, comparable to a targeted continous wave search. We analysed the detectable population separately in each detector: Advanced LIGO, Advanced Virgo, and the planned Einstein Telescope. In the calculation of the expected signal we neglect the frequency change of the signals due to the source’s spindown and the Earth’s motion with respect to the solar barycentre. Results. With conservative values for the neutron star evolution (a supernova rate of once per 100 years, initial ellipticity ϵ0 ≃ 10−5 with no decay of the ellipticity η = thub ≃ 104 Myr), the expected number of detected neutron stars is 0.15 (based on a simulation of 10 M stars) for the Advanced LIGO detector. A broader study of the parameter space (ϵ0, η) is presented. With the planned sensitivity for the Einstein Telescope, and assuming the same ellipiticity model, the expected detection number is 26.4 pulsars during a one-year observing run.

The Peculiar X-Ray Transient Swift J0840.7−3516: An Unusual Low-mass X-Ray Binary or a Tidal Disruption Event?

Abstract We report on the X-ray properties of the new transient Swift J0840.7−3516, discovered with Swift/BAT in 2020 February, using extensive data from Swift, MAXI, NICER, and NuSTAR. The source flux increased for ∼103 s after the discovery, decayed rapidly over ∼5 orders of magnitude in five days, and then remained almost constant over nine months. Large-amplitude short-term variations on timescales of 1–104 s were observed throughout the decay. In the initial flux rise, the source showed a hard power-law-shaped spectrum with a photon index of ∼1.0 extending up to ∼30 keV, above which an exponential cutoff was present. The photon index increased in the following rapid decay and became ∼2 at the end of the decay. A spectral absorption feature at 3–4 keV was detected in the decay. It is not straightforward to explain all the observed properties by any known class of X-ray sources. We discuss the possible nature of the source, including a Galactic low-mass X-ray binary with multiple extreme properties and a tidal disruption event by a supermassive black hole or a Galactic neutron star.