Pulsar Research Articles

Accretion tori around rotating neutron stars I. Structure, shape, and size

We present a full general relativistic analytic solution for a radiation-pressure-supported equilibrium fluid torus orbiting a rotating neutron star (NS). We applied previously developed analytical methods that include the effects of both the NS's angular momentum and quadrupole moment in the Hartle-Thorne geometry. The structure, size, and shape of the torus are explored, with a particular focus on the critically thick solution -- the cusp tori. For the astrophysically relevant range of NS parameters, we examined how our findings differ from those obtained for the Schwarzschild space-time. The solutions for rotating stars display signatures of an interplay between relativistic and Newtonian effects where the impact of the NS angular momentum and quadrupole moment are almost counterbalanced at a given radius. Nevertheless, the space-time parameters still strongly influence the size of tori, which can be shown in a coordinate-independent way. Finally, we discuss the importance of the size of the central NS which determines whether or not a surrounding torus exists. We provide a set of tools in a Wolfram Mathematica code, which establishes a basis for further investigation of the impact of the NSs’ super-dense matter equation of state on the spectral and temporal behaviour of accretion tori.

Accretion tori around rotating neutron stars II. Oscillations and precessions

The four characteristic oscillation frequencies of accretion flows (in addition to the Keplerian orbital frequency) are often discussed in the context of the time variability of black hole and neutron star (NS) low-mass X-ray binaries (LMXBs). These four frequencies are the frequencies of the axisymmetric radial and vertical epicyclic oscillations, and the frequencies of non-axisymmetric oscillations corresponding to the periastron (radial) and Lense-Thirring (vertical) precessions. In this context, we investigated the effect of the quadrupole moment of a slowly rotating NS and provide complete formulae for calculating these oscillation and precession frequencies, as well as convenient approximations. Simple formulae corresponding to the geodesic limit of a slender torus (and test-particle motion) and the limit of a marginally overflowing torus (a torus exhibiting a critical cusp) are presented, and more general approximate formulae are included to allow calculations for arbitrarily thick tori. We provide the Wolfram Mathematica code used for our calculations together with the C++ and PYTHON codes for calculating the frequencies. Our formulae can be used for various calculations regarding the astrophysical signatures of the NS super-dense matter equation of state. For instance, we demonstrate that even for a given fixed number of free parameters, a model that accounts for fluid flow precession matches the frequencies of twin-peak quasiperiodic oscillations observed in NS LMXBs better than a model that uses geodesic precession.

Irradiation effect and binary radio pulsars with giant companions

Context. Pulsars are neutron stars that rotate rapidly. Most of the pulsars in binary systems tend to spin faster than those in isolation. According to binary evolution theory, radio pulsars in binary systems can have various types of companion stars. However, binary radio pulsars with giant stars, helium stars, and black holes as companions are lacking so far. Aims. We aim to investigate the possible parameter space of binary radio pulsars with giant companions. Methods. We used the MESA stellar evolution code to consider the effect of irradiation in binary evolution and evolved a series of binary models. Results. We present the potential physical properties of binary radio pulsars with giant star companions in the framework of the classical recycling scenario and the irradiation model. We found that the parameter space of binary radio pulsars with giant companions can be greatly expanded by the effect of irradiation. Moreover, when irradiation is strong enough, submillisecond pulsars might be accompanied by giant companion stars. We also demonstrate a correlation between the timescale of the binary system as a low-mass X-ray binary and the irradiation efficiency. Conclusions. The birthrate problem between millisecond pulsars and low-mass X-ray binaries might not be resolved by the irradiation effect. A more detailed binary population synthesis is necessary. Our findings may offer guidance for observers that wish to locate binary radio pulsars with giant companions or submillisecond pulsars.

Discovery of a Millisecond Pulsar Associated with Terzan 6

Observations show that globular clusters (GCs) might be among the best places to find millisecond pulsars (MSPs). However, the GC Terzan 6 seems to be an exception without any pulsar discovered, although its high stellar encounter rate suggests that it harbors dozens of them. We report the discovery of the first radio pulsar, PSR J1750–3116A, likely associated with Terzan 6 in a search of C-band (4–8 GHz) data from the Robert C. Byrd Green Bank Telescope with a spin period of 5.33 ms and dispersion measure (DM) ≃ 383 pc cm−3. The mean flux density of this pulsar is approximately 3 μJy. The DM agrees well with predictions from the Galactic free electron density model, assuming a distance of 6.7 kpc for Terzan 6. PSR J1750–3116A is likely an isolated MSP, potentially formed through dynamical interactions, considering the core-collapsed classification and the exceptionally high stellar encounter rate of Terzan 6. This is the highest radio frequency observation that has led to the discovery of a pulsar in a GC to date. While L-band (1–2 GHz) observations of this cluster are unlikely to yield significant returns due to propagation effects, we predict that further pulsar discoveries in Terzan 6 will be made by existing radio telescopes at higher frequencies.

Coherent Radio Emission from “Main-sequence Radio Pulse Emitters”: A New Stellar Diagnostic to Probe 3D Magnetospheric Structures

Main-sequence radio pulse emitters (MRPs) are magnetic early-type stars that produce coherent radio emission observed in the form of periodic radio pulses. The emission mechanism behind this is the electron-cyclotron maser emission (ECME). Among all kinds of magnetospheric emission, ECME is unique due to its high directivity and intrinsically narrow bandwidth. The emission is also highly circularly polarized and the sign of polarization is opposite for the two magnetic hemispheres. This combination of properties makes ECME highly sensitive to the three-dimensional structures in the stellar magnetospheres. This is especially significant for late-B and A-type magnetic stars that do not emit other types of magnetospheric emission such as Hα, the key probe used to trace magnetospheric densities. In this paper, we use an ultra-wideband observation (0.4–2 GHz) of a late B-type MRP HD 133880 to demonstrate how we can extract information on plasma distribution from ECME. We achieve this by examining the differences in pulse arrival times (“lags”) as a function of frequencies and qualitatively comparing those with lags obtained by simulating ECME ray paths in hot stars’ magnetospheres. This reveals that the stellar magnetosphere has a disk-like overdensity inclined to the magnetic equator with a centrally concentrated density that primarily affects the intermediate frequencies (400–800 MHz). This result, which is consistent with the recent density model proposed for hotter centrifugally supported magnetospheres, lends support to the idea of a unifying model for magnetospheric operations in early-type stars, and also provides further motivation to fully characterize the ECME phenomenon in large-scale stellar magnetospheres.

Coherence of Multidimensional Pair Production Discharges in Polar Caps of Pulsars

We report on the first self-consistent multidimensional particle-in-cell numerical simulations of nonhomogeneous pair discharges in polar caps of rotation-powered pulsars. By introducing strong inhomogeneities in the initial plasma distribution in our simulations, we analyze the degree of self-consistently emerging coherence of discharges across magnetic field lines. In 2D, we study discharge evolution for a wide range of physical parameters and boundary conditions corresponding to both the absent and free escape of charged particles from the surface of a neutron star. We also present the results of the first 3D simulations of discharges in a polar cap with a distribution of the global magnetospheric current appropriate for a pulsar with 60° inclination angle. For all parameters, we find the coherence scale of pair discharges across magnetic field lines to be of the order of the gap height. We also demonstrate that the popular “spark” model of pair discharges is incompatible with the universally adopted force-free magnetosphere model: intermittent discharges fill the entire zone of the polar cap that allows pair cascades, leaving no space for discharge-free regions. Our findings disprove the key assumption of the spark model about the existence of isolated distinct discharge columns.

Simulation-based Inference of Radio Millisecond Pulsars in Globular Clusters

Millisecond pulsars (MSPs) are abundant in globular clusters (GCs), which offer favorable environments for their creation. While the advent of recent, powerful facilities led to a rapid increase in MSP discoveries in GCs through pulsation searches, detection biases persist. In this work, we investigate the ability of current and future detections in GCs to constrain the parameters of the MSP population in GCs through a careful study of their luminosity function. Parameters of interest are the number of MSPs hosted by a GC, as well as the mean and the width of their luminosity function, which are typically affected by large uncertainties. While, as we show, likelihood-based studies can lead to ill-behaved posteriors on the size of the MSP population, we introduce a novel, likelihood-free analysis, based on marginal neural ratio estimation, which consistently produces well-behaved posteriors. We focus on the GC Terzan 5 (or Ter 5), which currently counts 48 detected MSPs. We find that 158 −104+294 MSPs should be hosted in this GC, but the uncertainty on this number remains large. We explore the performance of our new method on simulated Terzan 5-like data sets mimicking possible future observational outcomes. We find that significant improvement on the posteriors can be obtained by adding a reliable measurement of the diffuse radio emission of the GC to the analysis or by improving the detection threshold of current radio pulsation surveys by at least a factor of 2.

Multiwavelength identification of millisecond pulsar candidates in the Galactic bulge

Context. The existence of a population of millisecond pulsars in the Galactic bulge is supported, along with other evidence, by the Fermi GeV excess, an anomalous γ-ray emission detected almost 15 years ago in the direction of the Galactic center. However, radio surveys searching for pulsations have not yet revealed bulge millisecond pulsars. Aims. Identifying promising bulge millisecond pulsar candidates is key to motivating pointed radio pulsation searches. Candidates are often selected among steep-spectrum or polarized radio sources, but multiwavelength information can also be exploited: The aim of this work is to pinpoint strong candidates among the yet unidentified X-ray sources. Methods. We investigated the multiwavelength counterparts of sources detected by the Chandra X-ray observatory that have spectral properties expected for millisecond pulsars in the Galactic bulge. We considered that ultraviolet, optical, and strong infrared counterparts indicate that an X-ray source is not a bulge pulsar, while a radio or a faint infrared counterpart makes it a promising candidate. Results. We identify a large population of more than a thousand X-ray sources without optical, ultraviolet, or strong infrared counterparts. Among them, five are seen for the first time in unpublished radio imaging data from the Very Large Array. We provide the list of promising candidates, for most of which follow-up pulsation searches are ongoing.

The Origins of Narrow Spectra of Fast Radio Bursts

Observations find that some fast radio bursts (FRBs) have extremely narrowband spectra, i.e., Δν/ν 0 ≪ 1. We show that, when the angular size of the emission region is larger than the Doppler beaming angle, the observed spectral width (Δν/ν 0) exceeds 0.58 due to the high-latitude effects for a source outside the magnetosphere, even when the spectrum in the source’s comoving frame is monochromatic. The angular size of the source for magnetospheric models of FRBs can be smaller than the Doppler beaming angle, in which case this geometric effect does not influence the observed bandwidth. We discuss various propagation effects to determine if any could transform a broad-spectrum radio pulse into a narrow spectrum signal at the observer’s location. We find that plasma lensing and scintillation can result in a narrow bandwidth in the observed spectrum. However, the likelihood of these phenomena being responsible for the narrow observed spectra with Δν/ν 0 < 0.58 in the fairly large observed sample of FRBs is exceedingly small.

The Northern Cross Fast Radio Burst project

Context. Fast radio bursts (FRBs) are energetic, millisecond-duration radio pulses observed at extragalactic distances and whose origins are still a subject of heated debate. A fraction of the FRB population have shown repeating bursts, however it’s still unclear whether these represent a distinct class of sources. Aims. We investigated the bursting behaviour of FRB 20220912A, one of the most active repeating FRBs known thus far. In particular, we focused on its burst energy distribution, linked to the source energetics, and its emission spectrum, with the latter directly related to the underlying emission mechanism. Methods. We monitored FRB 20220912A at 408 MHz with the Northern Cross radio telescope and at 1.4 GHz using the 32-m Medicina Grueff radio telescope. Additionally, we conducted 1.2 GHz observations taken with the upgraded Giant Meter Wave Radio Telescope (uGMRT) searching for a persistent radio source coincident with FRB 20220912A, and included high energy observations in the 0.3–10 keV, 0.4–100 MeV and 0.03–30 GeV energy range. Results. We report 16 new bursts from FRB 20220912A at 408 MHz during the period between October 16th 2022 and December 31st 2023. Their cumulative spectral energy distribution follows a power law with slope αE = −1.3 ± 0.2 and we measured a repetition rate of 0.19 ± 0.03 hr−1 for bursts having a fluence of ℱ ≥ 17 Jy ms. Furthermore, we report no detections at 1.4 GHz for ℱ ≥ 20 Jy ms. These non-detections imply an upper limit of β &lt; −2.3, with β being the 408 MHz – 1.4 GHz spectral index of FRB 20220912A. This is inconsistent with positive β values found for the only two known cases in which an FRB has been detected in separate spectral bands. We find that FRB 20220912A shows a decline of four orders of magnitude in its bursting activity at 1.4 GHz over a timescale of one year, while remaining active at 408 MHz. The cumulative spectral energy distribution (SED) shows a flattening for spectral energy Eν ≥ 1031 erg Hz−1, a feature seen thus far in only two hyperactive repeaters. In particular, we highlight a strong similarity between FRB 20220912A and FRB 20201124A, with respect to both the energy and repetition rate ranges. We also find a radio continuum source with 240 ± 36 μJy flux density at 1.2 GHz, centered on the FRB 20220912A coordinates. Finally, we place an upper limit on the γ to radio burst efficiency η to be η &lt; 1.5 × 109 at 99.7% confidence level, in the 0.4–30 MeV energy range. Conclusions. The strong similarity between the cumulative energy distributions of FRB 20220912A and FRB 20201124A indicate that bursts from these sources are generated via similar emission mechanisms. Our upper limit on β suggests that the spectrum of FRB 20220912A is intrinsically narrow-band. The radio continuum source detected at 1.2 GHz is likely due to a star formation environment surrounding the FRB, given the absence of a source compact on millisecond scales brighter than 48 μJy beam−1. Finally, the upper limit on the ratio between the γ and radio burst fluence disfavours a giant flare origin for the radio bursts unlike observed for the Galactic magnetar SGR 1806-20.

Follow-up Observations of Apparently One-off Sources from the Parkes Telescope

A small fraction of fast radio bursts (FRBs) have been observed with multiple bursts, whereas most Galactic sources emitting radio pulses are known to repeat. Here we present the results of follow-up observations of two FRBs and four rotating radio transients (RRATs). Among these, only one RRAT has been observed with repeating pulses, with an estimated period of around 1.297047 s. For comparison, we reanalysed the Parkes archival follow-up observations in CSIRO’s data archive for all apparently one-off sources discovered by the Parkes telescopes, including 13 RRATs and 29 FRBs. In total, 3 RRATs are suggested to be repeaters, but no repeating signals were detected from the other sources. Reporting details of the nondetection observations for the apparently one-off sources would help investigate their origins, and catastrophic scenarios are worth proposing for both extragalactic and Galactic sources.

Optimal radar ranging pulse to resolve two reflectors

Previous work established fundamental bounds on subwavelength resolution for the radar range resolution problem, called superradar []. In this work, we identify the optimal waveforms for distinguishing the range resolution between two reflectors of identical strength, leveraging results in quantum metrology. We discuss both the unnormalized optimal waveform as well as the best square-integrable pulse and their variants. Using orthogonal function theory, we give an explicit algorithm to optimize the wave pulse in finite time to have the best performance. We also explore range resolution estimation with unnormalized waveforms with multiparameter methods to also independently estimate loss and time of arrival. These results are consistent with the earlier single parameter approach of range resolution only and give deeper insight into the ranging estimation problem. Experimental results are presented using radio pulse reflections inside coaxial cables, showing robust range resolution smaller than a tenth of the inverse bandlimit, with uncertainties close to the derived Cramér-Rao bound. Published by the American Physical Society 2024

Long-term evolutionary links between the isolated neutron star populations

ABSTRACT We have investigated the evolutionary connections of the isolated neutron star (NS) populations including radio pulsars (RPs), anomalous X-ray pulsars (AXPs), soft gamma repeaters (SGRs), dim isolated NSs (XDINs), ‘high-magnetic field’ RPs (‘HBRPs’), central compact objects (CCOs), rotating radio transients (RRATs), and long-period pulsars (LPPs) in the fallback disc model. The model can reproduce these NS families as a natural outcome of different initial conditions (initial period, disc mass, and dipole moment, μ) with a continuous μ distribution in the $\sim 10^{27} - 5 \times 10^{30}$ G cm$^3$ range. Results of our simulations can be summarized as follows: (1) A fraction of ‘HBRPs’ with relatively high μ evolve into the persistent AXP/SGR properties, and subsequently become LPPs. (2) Persistent AXP/SGRs do not have evolutionary links with CCOs, XDINs, and RRATs. (3) For a wide range of μ, most RRATs evolve passing through RP or ‘HBRP’ properties during their early evolutionary phases. (4) A fraction of RRATs which have the highest estimated birth rate seem to be the progenitors of XDINs. (5) LPPs, whose existence was predicted by the fallback disc model, are the sources evolving in the late stage of evolution before the discs become inactive. These results provide concrete support to the ideas proposing evolutionary connections between the NS families to account for the ‘birth rate problem’, the discrepancy between the cumulative birth rate estimated for these systems and the core-collapse supernova rate.

Diffraction efficiency in acousto-optical interaction with short radio pulses

Development of acousto-optic modulators with high diffraction efficiency, which is achieved by using multi-frequency acoustic signals, is of great interest. In this regard, it seems interesting to study the possibility of implementing highly efficient diffraction for wideband signals in the form of radio pulses with short duration. The goal of presented article is to develop a mathematical model based on a quantum description of the interaction between a photon and a phonon, simulate various operating modes of an acousto-optical modulator and compare theoretical results with experimental ones. A comparison between experimental and theoretical results was made, demonstrating the possibility of obtaining high diffraction efficiency (more than 90%) for wideband signals. The proposed method makes it possible to determine the conditions under which maximum diffraction efficiency is observed, and it also seems possible to determine the parameters of the pulse sequence.

Sudden Polarization Angle Jumps of the Repeating Fast Radio Burst FRB 20201124A

We report the first detection of polarization angle orthogonal jumps, a phenomenon previously only observed from radio pulsars, from a fast radio burst (FRB) source FRB 20201124A. We find three cases of orthogonal jumps in over 2000 bursts, all resembling those observed in pulsar single pulses. We propose that the jumps are due to the superposition of two orthogonal emission modes that could only be produced in a highly magnetized plasma, and they are caused by the line of sight sweeping across a rotating magnetosphere. The shortest jump timescale is of the order of 1 millisecond, which hints that the emission modes come from regions smaller than the light cylinder of most pulsars or magnetars. This discovery provides convincing evidence that FRB emission originates from the complex magnetosphere of a magnetar, suggesting an FRB emission mechanism that is analogous to radio pulsars despite a huge luminosity difference between two types of objects.

Neutron stars in f(R,T) theory: slow rotation approximation

In this paper, we study the slowly rotating neutron stars in f(R, T) gravity based on Hartle-Thorne formalism. We first consider the simplest matter-geometry coupled modified gravity, namely f(R, T) = R + 2χ T. We compute the mass, radius, moment of inertia, change in radius, and binding energy due to rotation, eccentricity, quadrupole moment, and the tidal love number. The quantities, which are of the second order in angular velocity, like change in radius and binding energy due to rotation, eccentricity, and quadrupole moment, deviate more from their corresponding general relativistic counterparts in lighter neutron stars than heavier ones. Whereas the moment of inertia, which is of the first order in angular velocity, in f(R, T) = R + 2χ T modified gravity, barely diverges from the general relativistic one. The Equation of state-independent I-Love-Q relation retains in this f(R, T) modified gravity, and it coincides with the general relativistic ones within less than one percent even for the maximum allowed coupling parameters. We also study the slowly rotating neutron star in f(R, T) = R + αR 2 + 2χT up to first order their angular velocity. We calculate the mass, radius, and moment of inertia of neutron stars in this modified gravity. The results show that the impact of the matter-geometric coupling parameter is greater on lighter neutron stars in both of these modified gravity models.

Timing of Millisecond Pulsars in NGC 6752. III. On the Presence of Nonluminous Matter in the Cluster’s Core

Millisecond pulsars are subject to accelerations in globular clusters (GCs) that manifest themselves in both the first and second spin period time derivatives, and can be used to explore the mass distribution of the potentials they inhabit. Here we report on over 20 yr of pulsar timing observations of five millisecond radio pulsars in the core of the core-collapse GC NGC 6752 with the Parkes (Murriyang) and MeerKAT radio telescopes, which have allowed us to measure the proper motions, positions, and first and second time derivatives of the pulsars. The pulsar timing parameters indicate that all the pulsars in the core experience accelerations and jerks that can be explained only if an amount of nonluminous mass of at least 2.56 × 103 M ⊙ is present in the core of NGC 6752. On the other hand, our studies highly disfavor the presence of an intermediate-mass black hole at the center of the cluster, with a mass equal to or greater than ∼3000M ⊙.

Behavior of the Position and Ellipticity Angles at Polarization Mode Transitions in Pulsar Radio Emission

Polarization observations of radio pulsars show that abrupt transitions in the polarization vector’s position angle can be accompanied by large excursions in the vector’s ellipticity angle, suggesting the vector passes near the right or left circular pole of the Poincaré sphere. The behavior of the angles can be explained by a transition in dominance of the orthogonal polarization modes or a vector rotation caused by a change in the phase difference between the modes. Four polarization models are examined to quantify and understand the behavior of the angles at a mode transition: coherent polarization modes, partially coherent modes, incoherent modes with nonorthogonal polarization vectors, and incoherent orthogonal modes with an elliptically polarized emission component. In all four models, the trajectory of the mode transition on the Poincaré sphere follows the geodesic that connects the orientations of the mode polarization vectors. The results from the models can be similar, indicating that the interpretation of an observed transition within the context of a particular model is not necessarily unique. The polarization fraction of the emission and the average ellipticity angle depend upon the statistical character of the mode intensity fluctuations. The polarization fraction increases as the fluctuations increase. The excursion in ellipticity angle can be large when the mode intensities are quasi-stable and is suppressed when the intensity fluctuations are large.

A bright burst from FRB 20200120E in a globular cluster of the nearby galaxy M81

Fast radio bursts (FRBs) are immensely energetic millisecond-duration radio pulses. Observations indicate that nearby FRBs can be produced by old stellar populations, as suggested by the localization of the repeating source FRB 20200120E in a globular cluster of M81. Nevertheless, the burst energies of FRB 20200120E are significantly smaller than those of other cosmological FRBs. Here, we report the detection of a bright burst from FRB 20200120E in 1.1 – 1.7 GHz, with a fluence of approximately 30 Jy ms, which is more than 42 times larger than the previously detected bursts near 1.4 GHz frequency. It reaches one-third of the energy of the weakest burst from FRB 20121102A and is detectable at a distance exceeding 200 Mpc. Our finding bridges the gap between nearby and cosmological FRBs and indicates that FRBs hosted in globular clusters can be bright enough to be observable at cosmological distances.

A New Puzzling Periodic Signal in GeV Energies of the γ-Ray Binary LS I+61°303

LS I+61°303 is a high-mass X-ray binary system comprising a massive Be star and a rapidly rotating neutron star. Its spectral energy distribution across multiwavelengths categorizes it as a γ-ray binary system. In our analysis of LS I+61°303 using Fermi Large Area Telescope observations, we not only confirmed the three previously discussed periodicities of orbital, superorbital, and orbital–superorbital beat periods observed in multiwavelength observations, but also identified an additional periodic signal. This newly discovered signal exhibits a period of ∼26.3 days at a ∼7σ confidence level. Moreover, the power spectrum peak of the new signal gradually decreases as the energy increases across the energy ranges of 0.1–0.3, 0.3–1.0, and 1.0–500.0 GeV. Interestingly, a potential signal with a similar period was found in data obtained from the Owens Valley Radio Observatory 40 m telescope. We suggest that the newly discovered periodic signal may originate from a coupling between the orbital period and the retrograde stellar precession period.