Galactic Pulsars Research Articles

Constraining positron emission from pulsar populations with AMS-02 data

The cosmic-ray flux of positrons is measured with high precision by the space-borne particle spectrometer AMS-02. The hypothesis that pulsar wind nebulae (PWNe) can significantly contribute to the excess of the positron (e+) cosmic-ray flux has been consolidated after the observation of a γ-ray emission at TeV energies of a few degree size around Geminga and Monogem PWNe.In this work we undertake massive simulations of galactic pulsars populations, adopting different distributions for their position in the Galaxy, intrinsic physical properties, pair emission models, in order to overcome theincompleteness of the ATNF catalog. We fit the e+ AMS-02 data together with a secondary component due to collisions of primary cosmic rays with the interstellar medium.We find that several mock galaxies have a pulsar population able to explain the observed e+ flux, typically by few, bright sources. We determine the physical parameters of the pulsars dominating the e+ flux, and assess the impact of different assumptions on radial distributions, spin-down properties, Galactic propagation scenarios and e+ emission time.

Constraining positron emission from pulsar populations with AMS-02 data

The cosmic-ray flux of positrons is measured with high precision by the space-borne particle spectrometer AMS-02. The hypothesis that pulsar wind nebulae (PWNe) can significantly contribute to the excess of the positron (e+) cosmic-ray flux has been consolidated after the observation of a γ-ray emission at TeV energies of a few degree size around Geminga and Monogem PWNe. In this work we undertake massive simulations of galactic pulsars populations, adopting different distributions for their position in the Galaxy, intrinsic physical properties, pair emission models, in order to overcome the incompleteness of the ATNF catalogue. We fit the e + AMS-02 data together with a secondary component due to collisions of primary cosmic rays with the interstellar medium. We find that several mock galaxies have a pulsar population able to explain the observed e+ flux, typically by few, bright sources. We determine the physical parameters of the pulsars dominating the e + flux, and assess the impact of different assumptions on radial distributions, spin-down properties, Galactic propagation scenarios and e + emission time.

The thousand-pulsar-array programme on MeerKAT VII: polarisation properties of pulsars in the Magellanic Clouds

ABSTRACT The Magellanic Clouds are the only external galaxies known to host radio pulsars. The dispersion and rotation measures of pulsars in the Clouds can aid in understanding their structure, and studies of the pulsars themselves can point to potential differences between them and their Galactic counterparts. We use the high sensitivity of the MeerKAT telescope to observe 17 pulsars in the Small and Large Magellanic Clouds in addition to five foreground (Galactic) pulsars. We provide polarisation profiles for 18 of these pulsars, improved measurements of their dispersion and rotation measures, and derive the mean parallel magnetic field along the lines of sight. The results are broadly in agreement with expectations for the structure and strength of the magnetic field in the Large and Small Magellanic Clouds. The Magellanic Cloud pulsars have profiles which are narrower than expected from the period-width relationship and we show this is due to selection effects in pulsar surveys rather than any intrinsic difference between the population of Galactic and Magellanic objects.

Particle-in-Cell Simulations of Astrophysical Relativistic Jets

Astrophysical relativistic jets in active galactic nuclei, gamma-ray bursts, and pulsars is the main key subject of study in the field of high-energy astrophysics, especially regarding the jet interaction with the interstellar or intergalactic environment. In this work, we review studies of particle-in-cell simulations of relativistic electron–proton (e−−p+) and electron–positron (e±) jets, and we compare simulations that we have conducted with the relativistic 3D TRISTAN-MPI code for unmagnetized and magnetized jets. We focus on how the magnetic fields affect the evolution of relativistic jets of different compositions, how the jets interact with the ambient media, how the kinetic instabilities such as the Weibel instability, the kinetic Kelvin–Helmholtz instability and the mushroom instability develop, and we discuss possible particle acceleration mechanisms at reconnection sites.

Investigating γ -ray halos around three HAWC bright sources in Fermi -LAT data

Numerous extended sources around Galactic pulsars have shown significant $\gamma$-ray emission from GeV to TeV energies, revealing hundreds of TeV energy electrons scattering off of the underlying photon fields through inverse Compton scattering (ICS). HAWC TeV gamma-ray observations of few-degree extended emission around the pulsars Geminga and Monogem, and LAT GeV emission around Geminga, suggest that systems older than 100 kyr have multi-TeV $e^\pm$ propagating beyond the SNR-PWN system into the interstellar medium. Following the discovery of few $\gamma$-ray sources by HAWC at energies E$>100$~TeV, we investigate the presence of an extended $\gamma$-ray emission in Fermi-LAT data around the three brightest sources detected by HAWC up to 100 TeV. We find an extended emission of $\theta_{68} = 1.00^{+0.05}_{-0.07}$ deg around eHWC J1825-134 and $\theta_{68} = 0.71\pm0.10$ deg eHWC J1907+063. The analysis with ICS templates on Fermi-LAT data point to diffusion coefficient values which are significantly lower than the average Galactic one. When studied along with HAWC data, the $\gamma$-ray Fermi-LAT data provide invaluable insight into the very high-energy electron and positron parent populations.

First Discovery of New Pulsars and RRATs with CHIME/FRB

Abstract We report the discovery of seven new Galactic pulsars with the Canadian Hydrogen Intensity Mapping Experiment’s Fast Radio Burst (CHIME/FRB) backend. These sources were first identified via single pulses in CHIME/FRB, then followed up with CHIME/Pulsar. Four sources appear to be rotating radio transients, pulsar-like sources with occasional single-pulse emission with an underlying periodicity. Of those four sources, three have detected periods ranging from 220 ms to 2.726 s. Three sources have more persistent but still intermittent emission and are likely intermittent or nulling pulsars. We have determined phase-coherent timing solutions for the latter two. These seven sources are the first discovery of previously unknown Galactic sources with CHIME/FRB and highlight the potential of fast radio burst detection instruments to search for intermittent Galactic radio sources.

Identification of point sources in gamma rays using U-shaped convolutional neural networks and a data challenge

Context. At GeV energies, the sky is dominated by the interstellar emission from the Galaxy. With limited statistics and spatial resolution, accurately separating point sources is therefore challenging. Aims. Here we present the first application of deep learning based algorithms to automatically detect and classify point sources from gamma-ray data. For concreteness we refer to this approach as AutoSourceID. Methods. To detect point sources, we utilized U-shaped convolutional networks for image segmentation and k-means for source clustering and localization. We also explored the Centroid-Net algorithm, which is designed to find and count objects. Using two algorithms allows for a cross check of the results, while a combination of their results can be used to improve performance. The training data are based on 9.5 years of exposure from The Fermi Large Area Telescope (Fermi-LAT) and we used source properties of active galactic nuclei (AGNs) and pulsars (PSRs) from the fourth Fermi-LAT source catalog in addition to several models of background interstellar emission. The results of the localization algorithm are fed into a classification neural network that is trained to separate the three general source classes (AGNs, PSRs, and FAKE sources). Results. We compared our localization algorithms qualitatively with traditional methods and find them to have similar detection thresholds. We also demonstrate the robustness of our source localization algorithms to modifications in the interstellar emission models, which presents a clear advantage over traditional methods. The classification network is able to discriminate between the three classes with typical accuracy of ∼70%, as long as balanced data sets are used in classification training. We published online our training data sets and analysis scripts and invite the community to join the data challenge aimed to improve the localization and classification of gamma-ray point sources.

Reconciling hints on axion-like-particles from high-energy gamma rays with stellar bounds

It has been recently claimed by two different groups that the spectral modulation observed in gamma raysfrom Galactic pulsars and supernova remnants can be due to conversion of photons into ultra-light axion-like-particles (ALPs) in large-scale Galactic magnetic fields.While we show the required best-fit photon-ALP coupling, gaγ ∼ 2 × 10-10 GeV-1, to be consistent with constraints from observations of photon-ALPs mixing in vacuum, this is in conflict with other bounds, specifically from the CAST solar axion limit, from the helium-burning lifetime in globular clusters, and from the non-observations of gamma rays in coincidencewith SN 1987A. In order to reconcile these different results, we propose that environmental effects in matter would suppress the ALP production in dense astrophysical plasma, allowing to relax previous bounds and make them compatible with photon-ALP conversions in the low-density Galactic medium. If this explanation is correct, the claimed ALP signal would be on the reach of next-generations laboratory experiments such as ALPS II.

Fallback Supernova Assembly of Heavy Binary Neutron Stars and Light Black Hole–Neutron Star Pairs and the Common Stellar Ancestry of GW190425 and GW200115

Abstract The detection of the unusually heavy binary neutron star merger GW190425 marked a stark contrast to the mass distribution from known Galactic pulsars in double neutron star binaries and gravitational-wave source GW170817. We suggest here a formation channel for heavy binary neutron stars and light black hole–neutron star binaries in which massive helium stars, which had their hydrogen envelope removed during a common-envelope phase, remain compact and avoid mass transfer onto the neutron star companion, possibly avoiding pulsar recycling. We present three-dimensional simulations of the supernova explosion of the massive stripped helium star and follow the mass fallback evolution and the subsequent accretion onto the neutron star companion. We find that fallback leads to significant mass growth in the newly formed neutron star. This can explain the formation of heavy binary neutron star systems such as GW190425, as well as predict the assembly of light black hole–neutron star systems such as GW200115. This formation avenue is consistent with the observed mass–eccentricity correlation of binary neutron stars in the Milky Way. Finally, avoiding mass transfer suggests an unusually long spin-period population of pulsar binaries in our Galaxy.

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.

Particles accelerated in Hertzian magnetic dipole field of a pulsar * *Supported by grants from NSFC (11675182, 11690022)

Supernova remnants are supposed to be the most possible sources of cosmic rays. However, alternative sources of cosmic rays, such as an active galactic nucleus, gamma-ray bursts, and pulsars, have not be excluded. In this study, we investigate the possibility of cosmic rays being generated by pulsars. The pulsar is simply described as a rotational magnetic dipole, the so-called Hertzian magnetic dipole, an exact solution of the d'Alembert equations. In the rotational magnetic dipole field, charged particles experience an accelerated electric field with their radiation reaction. The particles, which are initially static out of the light cylinder radius, can be accelerated up to a high energy.

Timing observations of three Galactic millisecond pulsars

ABSTRACT We report observed and derived timing parameters for three millisecond pulsars (MSPs) from observations collected with the Parkes 64-m telescope, Murriyang. The pulsars were found during reprocessing of archival survey data by Mickaliger et al. One of the new pulsars (PSR J1546–5925) has a spin period P = 7.8 ms and is isolated. The other two (PSR J0921–5202 with P = 9.7 ms and PSR J1146–6610 with P = 3.7 ms) are in binary systems around low-mass (${\gt}0.2\, {\rm M}_{\odot }$) companions. Their respective orbital periods are 38.2 and 62.8 d. While PSR J0921–5202 has a low orbital eccentricity e = 1.3 × 10−5, in keeping with many other Galactic MSPs, PSR J1146–6610 has a significantly larger eccentricity, e = 7.4 × 10−3. This makes it a likely member of a group of eccentric MSP–helium white dwarf binary systems in the Galactic disc whose formation is poorly understood. Two of the pulsars are co-located with previously unidentified point sources discovered with the Fermi satellite’s Large Area Telescope, but no γ-ray pulsations have been detected, likely due to their low spin-down powers. We also show that, particularly in terms of orbital diversity, the current sample of MSPs is far from complete and is subject to a number of selection biases.

Classification of Fermi-LAT sources with deep learning using energy and time spectra

ABSTRACT Despite the growing number of gamma-ray sources detected by the Fermi-Large Area Telescope (LAT), about one-third of the sources in each survey remains of uncertain type. We present a new deep neural network approach for the classification of unidentified or unassociated gamma-ray sources in the last release of the Fermi-LAT catalogue (4FGL-DR2) obtained with 10 yr of data. In contrast to previous work, our method directly uses the measurements of the photon energy spectrum and time series as input for the classification, instead of specific, human-crafted features. Dense neural networks, and for the first time in the context of gamma-ray source classification recurrent neural networks, are studied in depth. We focus on the separation between extragalactic sources, i.e. active galactic nuclei, and Galactic pulsars, and on the further classification of pulsars into young and millisecond pulsars. Our neural network architectures provide powerful classifiers, with a performance that is comparable to previous analyses based on human-crafted features. Our benchmark neural network predicts that of the sources of uncertain type in the 4FGL-DR2 catalogue, 1050 are active galactic nuclei and 78 are Galactic pulsars, with both classes following the expected sky distribution and the clustering in the variability–curvature plane. We investigate the problem of sample selection bias by testing our architectures against a cross-match test data set using an older catalogue, and propose a feature selection algorithm using autoencoders. Our list of high-confidence candidate sources labelled by the neural networks provides a set of targets for further multiwavelength observations addressed to identify their nature. The deep neural network architectures we develop can be easily extended to include specific features, as well as multiwavelength data on the source photon energy and time spectra coming from different instruments.

Gravitational wave signature from phase transition of a combusting neutron star to quark star

ABSTRACT Fluctuation at the neutron star centre gives rise to a small deconfined quark core very close to the star centre. The density discontinuity at the quark–hadron boundary initiates a shock wave, which propagates outwards of the star. The shock has enough energy to combust nuclear matter to 2-flavour quark matter in the star. The 2-flavour quark matter is not stable and settles to a stable 3-flavour matter in the weakly interacting time-scale. In this paper, we study the conversion of 2-flavour matter to 3-flavour matter. We set-up a differential equation to convert the excess of down quarks to strange quarks involving weak reaction and diffusion of quarks. Calculating the reaction rate and diffusion, we solve the differential equation to find the velocity of the conversion front. As the conversion front moves out, the density profile changes, bringing about a change in the star’s quadrupole moment and thereby emitting gravitational waves (GWs). As the conversion process occurs, the temperature of the star rises, but the neutrino carries away the heat in some tens of ms. The GW amplitude of a colder star is well within the present detector capability, but the frequency is slightly on the higher side. Relatively hotter stars are on the boundary of present detectors and easily detectable with future detectors, and their frequency is also within the present detectability range. In comparison, phase transition from Galactic pulsars is easily detectable with present detectors.

Analyzing the Galactic Pulsar Distribution with Machine Learning

Abstract We explore the possibility of inferring the properties of the Galactic population of neutron stars through machine learning. In particular, in this paper we focus on their dynamical characteristics and show that an artificial neural network is able to estimate with high accuracy the parameters that control the current positions of a mock population of pulsars. For this purpose, we implement a simplified population-synthesis framework (where selection biases are neglected at this stage) and concentrate on the natal kick-velocity distribution and the distribution of birth distances from the Galactic plane. By varying these and evolving the pulsar trajectories in time, we generate a series of simulations that are used to train and validate a suitably structured convolutional neural network. We demonstrate that our network is able to recover the parameters governing the distribution of kick velocity and Galactic height with a mean relative error of about 10−2. We discuss the limitations of our idealized approach and study a toy problem to introduce selection effects in a phenomenological way by incorporating the observed proper motions of 216 isolated pulsars. Our analysis highlights that by increasing the sample of pulsars with accurate proper-motion measurements by a factor of ∼10, one of the future breakthroughs of the Square Kilometre Array, we might succeed in constraining the birth spatial and kick-velocity distribution of the neutron stars in the Milky Way with high precision through machine learning.

The FAST Galactic Plane Pulsar Snapshot survey: I. Project design and pulsar discoveries ⋆ ⋆News and views on this paper

Discovery of pulsars is one of the main goals for large radio telescopes. The Five-hundred-meter Aperture Spherical radio Telescope (FAST), that incorporates an L-band 19-beam receiver with a system temperature of about 20 K, is the most sensitive radio telescope utilized for discovering pulsars. We designed the snapshot observation mode for a FAST key science project, the Galactic Plane Pulsar Snapshot (GPPS) survey, in which every four nearby pointings can observe a cover of a sky patch of 0.1575 square degrees through beam-switching of the L-band 19-beam receiver. The integration time for each pointing is 300 seconds so that the GPPS observations for a cover can be made in 21 minutes. The goal of the GPPS survey is to discover pulsars within the Galactic latitude of ± 10° from the Galactic plane, and the highest priority is given to the inner Galaxy within ± 5°. Up to now, the GPPS survey has discovered 201 pulsars, including currently the faintest pulsars which cannot be detected by other telescopes, pulsars with extremely high dispersion measures (DMs) which challenge the currently widely used models for the Galactic electron density distribution, pulsars coincident with supernova remnants, 40 millisecond pulsars, 16 binary pulsars, some nulling and mode-changing pulsars and rotating radio transients (RRATs). The follow-up observations for confirmation of new pulsars have polarization-signals recorded for polarization profiles of the pulsars. Re-detection of previously known pulsars in the survey data also leads to significant improvements in parameters for 64 pulsars. The GPPS survey discoveries are published and will be updated at http://zmtt.bao.ac.cn/GPPS/.

Radio Counterpart Candidates to Unassociated 4FGL-DR2 Sources

Abstract For the duration of the Fermi Gamma-Ray Space Telescope's mission, approximately one-third of the point sources detected have been noted as “unassociated,” meaning that they seem to have no known counterpart at any other wavelength/frequency. This mysterious part of the gamma-ray sky is perhaps one of the largest unknowns in current astronomical pursuits and as such has been probed extensively by various techniques at various frequencies. Radio frequencies have been perhaps one of the most fruitful, producing a large fraction of the identified and associated active galactic nuclei and pulsars noted in each update of the point-source catalogs. Here we present a catalog of 7432 radio counterpart candidates for unassociated gamma-ray fields in the Second Data Release of the Fourth Fermi Point Source Catalog (4FGL-DR2). A description of the catalog and source types is provided followed by a discussion that demonstrates how the results of this work will aid new associations and identifications. As part of this work, we also present a first catalog derived from “quicklook” images of the Very Large Array Sky Survey.

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.

Research progress of Dark Matter Particle Explorer

Quite a lot of astronomical and cosmological observations reveal that dark matter widely exists in our Universe. It occupies nearly 1/4 of the total energy budget of the Universe, and is about 5 times more abundant than the normal matter. The nature of dark matter is one of the most important fundamental questions of modern physics and astronomy. To probe dark matter particles is also highly competitive among many countries around the world. China has launched the Dark Matter Particle Explorer (DAMPE; also known as “Wukong”) to indirectly detect dark matter via high-precision observations of high-energy electrons and gamma-rays. As a high-energy particle detector, the DAMPE data can also be useful for probing the origin and propagation of cosmic rays, which is another century-long problem, and the energetic activities of extreme objects such as black holes and neutron stars. DAMPE has the highest energy resolution for electron and gamma-ray measurements, and the highest electron-proton discrimination capability above TeV energies. Using the data recorded by DAMPE, we obtain by far the most precise measurements of the spectra of cosmic ray electrons and protons, and reveal new spectral features which were not discovered or clarified before. For the total electron plus positron spectrum, the DAMPE result covers an energy range of 25 GeV–4. 6 TeV, which reaches the highest energy among all the space direct detections. A spectral hardening around 50 GeV, and a softening around 0.9 TeV have been revealed by DAMPE. Some hint of fine structures of the spectrum is also shown which requires more data for confirmation or exclusion. The spectral features of the DAMPE electron plus positron spectrum have very interesting and important implications on the origin and distribution of high-energy electron sources. Particularly, the 0.9 TeV break may reveal the fact that sources of TeV electrons distribute discretely in the Galaxy. The energy range coverage of the DAMPE proton spectrum is from 40 GeV to 100 TeV. For the first time, DAMPE measures the proton spectrum up to 100 TeV with a high statistics. Besides the hundreds of GeV spectral hardening, the DAMPE measurement further reveals a clear spectral softening around 14 TeV. The softening feature of the proton spectrum may be an imprint of nearby cosmic ray acceleration sources. Alternatively, the spectral structures of protons may be due to the superposition of multiple source populations. High-energy gamma-rays have also been robustly identified from the DAMPE data. The analysis shows the detection of luminous diffuse emission from the Galactic plane and nearly 250 pointlike sources in the gamma-ray sky, most of which are active galactic nuclei and pulsars. The DAMPE results offer important data in constraining/understanding the properties of dark matter and cosmic rays. We review the design, operation and data analysis of DAMPE in this work, with a focus on the analyses of spectra of cosmic ray electrons and protons.

An excess radio signal in the Abell 4038 cluster

ABSTRACT In the past decade, various instruments, such as the Large Area Telescope (LAT) on the Fermi Gamma Ray Space Telescope, the Alpha Magnetic Spectrometer (AMS) and the Dark Matter Particle Explorer(DAMPE), have been used to detect the signals of annihilating dark matter in our Galaxy. Although some excesses of gamma rays, antiprotons and electrons/positrons have been reported and are claimed to be dark matter signals, the uncertainties of the contributions of Galactic pulsars are still too large to confirm the claims. In this paper, we report on a possible radio signal of annihilating dark matter manifested in the archival radio continuum spectral data of the Abell 4038 cluster. By assuming a thermal annihilation cross-section and comparing the dark matter annihilation model with the null hypothesis (cosmic ray emission without dark matter annihilation), we obtain very large test statistic (TS) values, TS > 45, for four popular annihilation channels, which correspond to more than 6σ statistical preference. This reveals a possible potential signal of annihilating dark matter. In particular, our results are also consistent with the recent claims of dark matter mass, m ≈ 30–50 GeV, annihilating via the $\rm b\bar{b}$ quark channel with the thermal annihilation cross-section. However, at this time, we cannot exclude the possibility that a better background cosmic ray model could explain the spectral data without recourse to dark matter annihilations.