Neutron Star Research Articles

Follow-up Timing of 12 Pulsars Discovered in Commensal Radio Astronomy FAST Survey

We present phase-connected timing ephemerides, polarization pulse profiles, and Faraday rotation measurements of 12 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in the Commensal Radio Astronomy FAST Survey. The observational data for each pulsar span at least 1 yr. Among them, PSR J1840+2843 shows subpulse drifting, and five pulsars are detected to exhibit pulse nulling phenomena. PSR J0640−0139 and PSR J2031−1254 are isolated millisecond pulsars (MSPs) with stable spin-down rates ( Ṗ ) of 4.8981(6) × 10−20 s s−1 and 6.01(2) × 10−21 s s−1, respectively. Additionally, one pulsar (PSR J1602−0611) is in a neutron star−white dwarf (WD) binary system with an 18.23-day orbit and a companion of ≤0.65 M ⊙. PSR J1602−0611 has a spin period, companion mass, and orbital eccentricity that are consistent with the theoretical expectations for MSP−helium WD (He WD) systems. Therefore, we believe that it might be an MSP−He WD binary system. The locations of PSR J1751−0542 and PSR J1840+2843 on the P−Ṗ diagram are beyond the traditional death line. This indicates that FAST has discovered some low- Ė pulsars, contributing new samples for testing pulsar radiation theories. We estimated the distances of these 12 pulsars based on NE2001 and YMW16 electron density models, and our work enhances the data set for investigating the electron density model of the Galaxy.

Electromagnetic signatures of black hole clusters in the center of super-Eddington galaxies

Supermassive black holes (SMBHs) at the centers of active galaxies are fed by accretion disks that radiate from the infrared or optical to the X-ray bands. Several types of objects can orbit SMBHs, including massive stars, neutron stars, clouds from the broad- and narrow-line regions, and X-ray binaries. Isolated black holes with a stellar origin (BHs of ∼10 M⊙) should also be present in large numbers within the central parsec of the galaxies. These BHs are expected to form a cluster around the SMBH as a result of the enhanced star formation rate in the inner galactic region and the BH migration caused by gravitational dynamical friction. However, except for occasional microlensing effects on background stars or gravitational waves from binary BH mergers, the presence of a BH population is hard to verify. In this paper, we explore the possibility of detecting electromagnetic signatures of a central cluster of BHs when the accretion rate onto the central SMBH is greater than the Eddington rate. In these supercritical systems, the accretion disk launches powerful winds that interact with the objects orbiting the SMBH. Isolated BHs can capture matter from this dense wind, leading to the formation of small accretion disks around them. If jets are produced in these single microquasars, they could be sites of particle acceleration to relativistic energies. These particles in turn are expected to cool by various radiative processes. Therefore, the wind of the SMBH might illuminate the BHs through the production of both thermal and nonthermal radiation.

Swiftly Chasing Gravitational Waves across the Sky in Real Time

We introduce a new capability of the Neil Gehrels Swift Observatory, dubbed “continuous commanding,” that achieves 10 s latency response time on orbit to unscheduled target-of-opportunity requests received on the ground. We show that this will allow Swift to respond to premerger (early-warning) gravitational-wave (GW) detections, rapidly slewing the Burst Alert Telescope (BAT) across the sky to place the GW origin in the BAT field of view at or before merger time. This will dramatically increase the GW/gamma-ray burst (GRB) codetection rate and enable prompt arcminute localization of a neutron star merger. We simulate the full Swift response to a GW early-warning alert, including input sky maps produced at different early-warning times, a complete model of the Swift attitude control system, and a full accounting of the latency between the GW detectors and the spacecraft. 60 s of early warning can double the rate of a prompt GRB detection with arcminute localization, and 140 s guarantees observation anywhere on the unocculted sky, even with localization areas ≫1000 deg2. While 140 s is beyond current GW detector sensitivities, 30–70 s is achievable today. We show that the detection yield is now limited by the latency of LIGO/Virgo cyberinfrastructure and motivate a focus on its reduction. Continuous commanding has been integrated as a general capability of Swift, significantly increasing its versatility in response to the growing demands of time-domain astrophysics. We demonstrate this potential on an externally triggered fast radio burst (FRB), slewing 81° across the sky, and collecting X-ray and UV photons from the source position <150 s after the trigger was received from the Canadian Hydrogen Intensity Mapping Experiment, thereby setting the earliest and deepest such constraints on high-energy activity from nonrepeating FRBs. The Swift Team invites the community to consider and propose novel scientific applications of ultra-low-latency UV, X-ray, and gamma-ray observations.

Capability of Searching for Kilonova Associated with a Short Gamma-Ray Burst by SVOM

In spite of the importance of studying the cosmic generation of heavy elements through the r-process, the detection of a kilonova resulting from the merger of a neutron star binary is still a challenging task. In this paper, we show that the Visible Telescope (VT) onboard the ongoing SVOM space mission is powerful for identifying kilonova candidates associated with short gamma-ray bursts up to a distance of 600 Mpc. A significant color variation, turning blue and then turning red, is revealed by calculating the light curves in both red and blue channels of VT by a linear combination of an afterglow and an associated kilonova. The maximum color variation is as high as ∼0.5–1 mag, which is far larger than the small photometry error of ∼0.2 mag of VT for a point source with a brightness of 23 mag. Up to a distance of 600 Mpc, ∼1–2 kilonova candidates per year are predicted to be identified by VT.

Interpreting Mass and Radius Measurements of Neutron Stars with Dark Matter Halos

The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

Energy-dependent and Energy-integrated Two-moment General-relativistic Neutrino Transport Simulations of a Hypermassive Neutron Star

We compare two-moment-based energy-dependent and three variants of energy-integrated neutrino transport general-relativistic magnetohydrodynamics simulations of a hypermassive neutron star. To study the impacts due to the choice of the neutrino transport schemes, we perform simulations with the same setups and input neutrino microphysics. We show that the main differences between energy-dependent and energy-integrated neutrino transport are found in the disk and ejecta properties, as well as in the neutrino signals. The properties of the disk surrounding the neutron star and the ejecta in energy-dependent transport are very different from the ones obtained using energy-integrated schemes. Specifically, in the energy-dependent case, the disk is more neutron-rich at early times and becomes geometrically thicker at later times. In addition, the ejecta is more massive and, on average, more neutron-rich in the energy-dependent simulations. Moreover, the average neutrino energies and luminosities are about 30% higher. Energy-dependent neutrino transport is necessary if one wants to better model the neutrino signals and matter outflows from neutron star merger remnants via numerical simulations.

The Meissner effect in neutron stars

Abstract We present the first model aimed at understanding how the Meissner effect in a young neutron star affects its macroscopic magnetic field. In this model, field expulsion occurs on a dynamical timescale, and is realised through two processes that occur at the onset of superconductivity: fluid motions causing the dragging of field lines, followed by magnetic reconnection. Focussing on magnetic fields weaker than the superconducting critical field, we show that complete Meissner expulsion is but one of four possible generic scenarios for the magnetic-field geometry, and can never expel magnetic flux from the centre of the star. Reconnection causes the release of up to ∼5 × 1046 erg of energy at the onset of superconductivity, and is only possible for certain favourable early-phase dynamics and for pre-condensation fields 1012 G ≲ B ≲ 5 × 1014 G. Fields weaker or stronger than this are predicted to thread the whole star.

X-ray and optical observations of the millisecond pulsar binary PSR J1431–4715

We present the first X-ray observation of the energetic millisecond pulsar binary PSR J1431−4715, performed with XMM-Newton and complemented with fast optical multi-band photometry acquired with the ULTRACAM instrument at ESO-NTT. It is found as a faint X-ray source without a significant orbital modulation. This contrasts with the majority of systems that instead display substantial X-ray orbital variability. The X-ray spectrum is dominated by non-thermal emission and, due to the lack of orbital modulation, does not favour an origin in an intrabinary shock between the pulsar and companion star wind. While thermal emission from the neutron star polar cap cannot be excluded in the soft X-rays, the dominance of synchrotron emission favours an origin in the pulsar magnetosphere that we describe at both X-ray and gamma-ray energies with a synchro-curvature model. The optical multi-colour light curve folded at the 10.8 h orbital period is double-humped and dominated by ellipsoidal effects, but also affected by irradiation. The ULTRACAM light curves are fit with several models encompassing direct heating and a cold spot, or heat redistribution after irradiation either through convection or convection plus diffusion. Despite the inability to constrain the best irradiation models, the fits provide consistent system parameters, giving an orbital inclination of 59 ± 6° and a distance of 3.1 ± 0.3 kpc. The companion is found to be an F-type star, underfilling its Roche lobe (fRL = 73 ± 4%) with a mass of 0.20 ± 0.04 M⊙, confirming the redback status, but hotter than the majority of redbacks. The stellar dayside and nightside temperatures of 7500 K and 7400 K, respectively, indicate a weak irradiation effect on the companion, likely due to its high intrinsic luminosity. Although the pulsar mass cannot be precisely derived, a heavy (1.8−2.2 M⊙) neutron star is favoured.

A census of the Sun’s ancestors and their contributions to the Solar System chemical composition

In this work, we compute the rates and numbers of different types of stars and phenomena (supernovae, novae, white dwarfs, merging neutron stars, black holes) that contributed to the chemical composition of the Solar System. During the Big Bang, only light elements formed, while all the heavy ones, from carbon to uranium and beyond, have since been created inside stars. Stars die and release the newly formed elements into the interstellar gas. This process is called ‘chemical evolution’. In particular, we analyse the death rates of stars of all masses, whether they die quiescently or explosively. These rates and total star numbers are computed in the context of a revised version of the two-infall model for the chemical evolution of the Milky Way, which reproduces the observed abundance patterns of several chemical species, the global solar metallicity, and the current gas, stellar, and total surface mass densities relatively well. We also compute the total number of stars ever born and still alive as well as the number of stars born up to the formation of the Solar System with mass and metallicity like those of the Sun. This latter number accounts for all the possible existing Solar systems that can host life in the solar vicinity. We conclude that, among all the stars (from 0.8 to 100 M⊙) that were born and died from the Big Bang up until the Solar System formation epoch and that contributed to its chemical composition, 93.00% were stars that died as single white dwarfs (without interacting significantly with a companion star) and originated in the mass range of 0.8–8 M⊙, while 5.24% were neutron stars and 0.73% were black holes, both originating from core-collapse supernovae (M &gt; 8 M⊙); 0.64% were Type Ia supernovae and 0.40% were nova systems, both originating from the same mass range as the white dwarfs. The number of stars similar to the Sun born from the Big Bang up until the formation of the Solar System, with metallicity in the range 12+log(Fe/H)= 7.50 ± 0.04 dex, is ~31•107, and in particular our Sun is the ~2.61• 107-th star of this kind.

Two-dimensional cylindrical magnetosonic shock waves in a relativistic degenerated plasma

Abstract In this paper, the characteristics of two-dimensional magnetosonic (MS) shock waves have been studied in a nonplanar relativistic degenerate collisional magnetoplasma whose constituents are non-degenerate warm ions fluid and relativistic degenerated electrons. Employing fluid model equations for such plasma along with Maxwell equations, a set of magnetohydrodynamic (MHD) model equations is obtained. Based on the newly obtained MHD equations, a Burgers-Kadomstev-Petviashvili equation (which describes shock wave structures) is derived in cylindrical geometry using the reductive perturbation technique. The considered plasma system was investigated under the impacts of spin-magnetization, relativistic degeneracy, cylindrical geometry, and dissipation. Numerical results revealed that the relativistic degeneracy, dissipation, and electron spin-magnetization as well as nonplanar geometry significantly altered the MS shock wave properties. Interestingly, it’s found that the shock nature changes and turns into new structures attributed to the transverse perturbation and cylindrical geometry. The implications of our investigation may be applicable to relativistic nonplanar quantum plasmas in astrophysical environments, particularly neutron stars and white dwarfs.&amp;#xD;

Pulsar-wind nebulae meeting the circumstellar media of their progenitors

A significative fraction of high-mass stars sail away through the interstellar medium of the galaxies. Once they evolved and died via a core-collapse supernova, a magnetised, rotating neutron star (a pulsar) is usually left over. The immediate surroundings of the pulsar is the pulsar wind, which forms a nebula whose morphology is shaped by the supernova ejecta and channelled into the circumstellar medium of the progenitor star in the pre-supernova time. Irregular pulsar-wind nebulae display a large variety of radio appearances, screened by their interacting supernova blast wave, or harbour asymmetric up–down emission. Here, we present a series of 2.5-dimensional (2 dimensions for the scalar quantities plus a toroidal component for the vectors) black non-relativistic magneto-hydrodynamical simulations exploring the evolution of the pulsar-wind nebulae generated by a red supergiant and a Wolf-Rayet massive supernova progenitor, moving with Mach number $M=1$ and $M=2$ into the warm phase of the Galactic plane. black In such a simplified approach, the progenitor's direction of motion, the local ambient medium magnetic field, and the progenitor and pulsar axis of rotation, are all aligned; this restricted our study to peculiar pulsar-wind nebula of high-equatorial-energy flux. We find that the reverberation of the termination shock of the pulsar-wind nebulae, when sufficiently embedded into its dead stellar surroundings and interacting with the supernova ejecta, is asymmetric and differs greatly as a function of the past circumstellar evolution of its progenitor, which reflects into their projected radio synchrotron emission. This mechanism is particularly at work in the context of remnants involving slowly moving or very massive stars. We find that the mixing of material in plerionic core-collapse supernova remnants is strongly affected by the asymmetric reverberation in their pulsar-wind nebulae.

Lower-mass-gap Black Holes in Dense Star Clusters

The existence of compact stellar remnants in the mass range 2–5 M ⊙ has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5 M ⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte Carlo N-body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc−3 yr−1.

First-order phase transitions in the cores of neutron stars

First-order phase transitions in the cores of neutron stars

Probing the Propeller Regime With Wide Neutron Star Binaries

ABSTRACTWe model the duration of the propeller stage in wide binary systems with neutron stars and calculate the time of accretion onset for various propeller models. We apply our modeling to the symbiotic X‐ray binary SWIFT J0850.8‐4219. Unless a propeller with a very slow spin‐down is operating, it is very improbable to find a system similar to SWIFT J0850.8‐4219 at the propeller stage. Then we model the evolution of a neutron star in a binary with a solar‐like companion. We calculate for which orbital separations and magnetic fields a neutron star can start to accrete while the companion is still on the Main sequence. We demonstrate that for the magnetic field G neutron stars at the orbital separation AU do not reach the propeller stage. In the case of a slow propeller spin‐down, neutron stars never start to accrete. For the more rapid propeller spin‐down, a neutron star can start to accrete or spend a long time at the propeller stage depending on the parameters.

Editor's Note: The X‐Ray Mysteries of Neutron Stars and White Dwarfs

Editor's Note: The X‐Ray Mysteries of Neutron Stars and White Dwarfs

A Comprehensive Study of Thermonuclear X-Ray Bursts from 4U 1820–30 with NICER: Accretion Disk Interactions and a Candidate Burst Oscillation

We present the results obtained from timing and spectral studies of 15 thermonuclear X-ray bursts from 4U 1820–30 observed with the Neutron Star Interior Composition Explorer (NICER) during its 5 yr of observations between 2017 and 2022. All bursts showed clear signs of photospheric radius expansion (PRE), where the neutron star (NS) photosphere expanded more than 50 km above the surface. One of the bursts produced a superexpansion with a blackbody emission radius of 902 km for the first time with NICER. We searched for burst oscillations in all 15 bursts and found evidence of a coherent oscillation at 716 Hz in a burst, with a 2.9σ detection level based on Monte Carlo simulations. If confirmed with future observations, 4U 1820–30 would become the fastest-spinning NS known in X-ray binary systems. The fractional rms amplitude of the candidate burst oscillation was found to be 5.8% in the energy range of 3–10 keV. Following the variable persistent model from burst time-resolved spectroscopy, an anticorrelation is seen between the maximum scaling factor value and the (preburst) persistent flux. We detected a low value of ionization at the peak of each burst based on reflection modeling of burst spectra. A partially interacting inner accretion disk or a weakly ionized outer disk may cause the observed ionization dip during the PRE phase.

Average speed of sound in neutron stars

Average speed of sound in neutron stars

A View of the Long-term Spectral Behavior of Ultracompact X-Ray Binary 4U 0614+091

In this study, we examine 51 archival NICER observations and six archival NuSTAR observations of the neutron star ultracompact X-ray binary (UCXB) 4U 0614+091, which span over 5 yr. The source displays persistent reflection features, so we use a reflection model designed for UCXBs, with overabundant carbon and oxygen (xillverCO) to study how various components of the system vary over time. The flux of this source is known to vary quasiperiodically on a timescale of a few days, so we study how the various model components change as the overall flux varies. The flux of most components scales linearly with the overall flux, while the power law, representing coronal emission, is anticorrelated, as expected. This is consistent with previous studies of the source. We also find that during observations of the high-soft state, the disk emissivity profile as a function of radius becomes steeper. We interpret this as the corona receding to be closer to the compact object during these states, at which point the assumed power law illumination of xillverCO may be inadequate to describe the illumination of the disk.

Binary neutron star merger offsets from their host galaxies. GW 170817 as a case study

The locations of binary neutron star (BNS) mergers within their host galaxies encode the systemic kicks that these systems received in the supernova aftermath. We investigate how the galactic potential and the systemic kicks shape the offset distribution of BNS mergers with a case study of GW 170817 and its host NGC 4993. We derived dynamical constraints on the host potential from integral field spectroscopy with Jeans anisotropic modelling. We evolved the trajectories of synthetic BNSs from the BPASS code in the galactic potential, using two different kick prescriptions to investigate how the observed offsets might differentiate between these two possibilities. The simulation was repeated after swapping the host potential with that of a dwarf galaxy, to test the effect of the potential on the offsets. The location of GW 170817 is entirely consistent with our predictions, regardless of large or small kicks, because the strong potential of NGC 4993 is only diagnostic of very large kicks. In galaxies of similar or greater mass, large offsets can constrain large kicks, while small offsets do not provide much information. In an old dwarf galaxy, on the other hand, small offsets can constrain small kicks, while large offsets would prevent host association.

Accreting Neutron Stars in 3D General-relativistic Magnetohydrodynamic Simulations: Jets, Magnetic Polarity, and the Interchange Slingshot

Abstract Accreting neutron stars differ from black holes by the presence of the star’s own magnetic field, whose interaction with the accretion flow is a central component in understanding these systems’ disk structure, outflows, jets, and spin evolution. It also introduces an additional degree of freedom, as the stellar dipole can have any orientation relative to the inner disk’s magnetic field. We present a suite of 3D general-relativistic magnetohydrodynamic simulations in which we investigate the two extreme polarities, with the dipole field being either parallel or antiparallel to the initial disk field, in both the accreting and propeller states. When the magnetosphere truncates the disk near or beyond the corotation radius, most of the system’s properties, including the relativistic jet power, are independent of the star–disk relative polarity. However, when the disk extends well inside the corotation radius, in the parallel orientation the jet power is suppressed and the inner disk is less dense and more strongly magnetized. We suggest a physical mechanism that may account for this behavior—the interchange slingshot—and discuss its astrophysical implications. When the star is in the rapidly accreting regime, which in most cases will be associated with strong spin-up, we expect large observational differences between the two magnetic orientations. This may be reflected in increased variability as the accretion flow drags in successive magnetic structures of varying polarity.