Pulsar Age Research Articles

A binary supernova OB-runaway candidate inside Berkeley 97

OB-runaway stars ejected by the binary supernova mechanism can be found near young open star clusters. In this paper, we present an OB-runaway candidate as a pre-SN binary companion to the progenitor of the pulsar PSR J2238$+$5903 inside the young open star cluster Berkeley 97. We tried to find a kinematic outlier based on Gaia DR3 proper motions and parallaxes to be the pre-supernova binary companion to the progenitor of the pulsar. We took the spectra of two bright early B-type stars of the cluster, determined their effective temperature and surface gravity, and updated the parameters of the cluster. Through isochrone fitting of the color-magnitude diagram of the star cluster, we identified the members and determined the stellar parameters of the runaway star. Two bright members of the cluster, HD 240015 and HD 240016, are massive stars with spectral types of B0.5II and B1.5II and effective temperatures of $T_ eff K and $T_ eff K, respectively, as well as surface gravities of $ g cm/s^2 We find that Berkeley 97 is a star cluster with an age of $ age yr )=7.1$, an uncertainty of $<0.1$ dex, and an interstellar extinction of $A_ V mag. The runaway star has an effective temperature of $T_ eff K with a surface gravity of $ g cm/s^2 (B8V type star). By tracing back the proper motion of the runaway star, the explosion center was found for different possible pulsar ages of 10, 20, and 26.6 kyr. The pulsar moving out from the 20 kyr position must have a space velocity of $ km $, which is consistent with the general pulsar velocity distribution. This supports the idea that the pulsar originated from the cluster as a result of a binary supernova. Despite its young age, $ kyr, the supernova remnant is not visible.

Pulsar Glitch Activities: The Spin Parameters Approach

Glitch activity refers to the mean increase in pulsar spin frequency per year due to rotational glitches. It is an important tool for studying super-nuclear matter using neutron star interiors as templates. Glitch events are typically observed in the spin frequency (ν) and frequency derivative ( ν̇ ) of pulsars. The rate of glitch recurrence decreases as the pulsar ages, and the activity parameter is usually measured by linear regression of cumulative glitches over a given period. This method is effective for pulsars with multiple regular glitch events. However, due to the scarcity of glitch events and the difficulty of monitoring all known pulsars, only a few have multiple records of glitch events. This limits the use of the activity parameter in studying neutron star interiors with multiple pulsars. In this study, we examined the relationship between the activity parameters and pulsar spin parameters (spin frequency, frequency derivative, and pulsar characteristic age). We found that a quadratic function provides a better fit for the relationship between activity parameters and spin parameters than the commonly used linear functions. Using this information, we were able to estimate the activity parameters of other pulsars that do not have records of glitches. Our analysis shows that the relationship between the estimated activity parameters and pulsar spin parameters is consistent with that of the observed activity parameters in the ensemble of pulsars.

Confronting the Neutron Star Population with Inverse Cascades

The origin and evolution of magnetic fields of neutron stars from birth have long been a source of debate. Here, motivated by recent simulations of the Hall cascade with magnetic helicity, we invoke a model where the large-scale magnetic field of neutron stars grows as a product of small-scale turbulence through an inverse cascade. We apply this model to a simulated population of neutron stars at birth and show how this model can account for the evolution of such objects across the diagram, explaining both pulsar and magnetar observations. Under the assumption that small-scale turbulence is responsible for large-scale magnetic fields, we place a lower limit on the spherical harmonic degree of the energy-carrying magnetic eddies of ≈40. Our results favor the presence of a highly resistive pasta layer at the base of the neutron star crust. We further discuss the implications of this paradigm on direct observables, such as the nominal age and braking index of pulsars.

MeerKAT caught a Mini Mouse: serendipitous detection of a young radio pulsar escaping its birth site

ABSTRACT In MeerKAT observations pointed at a Galactic X-ray binary located on the Galactic plane, we serendipitously discovered a radio nebula with cometary-like morphology. The feature, which we named ‘the Mini Mouse’ based on its similarity with the previously discovered ‘Mouse’ nebula, points back towards the previously unidentified candidate supernova remnant G45.24+0.18. We observed the location of the Mini Mouse with MeerKAT in two different observations, and we localized with arcsecond precision the 138-ms radio pulsar PSR J1914+1054g, recently discovered by the FAST telescope, to a position consistent with the head of the nebula. We confirm a dispersion measure of about 418 pc cm−3 corresponding to a distance between 7.8 and 8.8 kpc based on models of the electron distribution. Using our accurate localization and two period measurements spaced 90 d apart, we calculate a period derivative of (2.7 ± 0.3) × 10 −14 s s−1. We derive a characteristic age of approximately 82 kyr and a spin-down luminosity of 4 × 1035 erg s−1. For a pulsar age comparable with the characteristic age, we find that the projected velocity of the neutron star is between 320 and 360 km s−1 if it was born at the location of the supernova remnant. The size of the proposed remnant appears small if compared with the pulsar characteristic age; however, the relatively high density of the environment near the Galactic plane could explain a suppressed expansion rate and thus a smaller remnant.

Pulsars do not produce sharp features in the cosmic-ray electron and positron spectra

Pulsars are considered to be the leading explanation for the excess in cosmic-ray positrons detected by PAMELA and AMS-02. A notable feature of standard pulsar models is the sharp spectral cutoff produced by the increasingly efficient cooling of very-high-energy electrons by synchrotron and inverse-Compton processes. This spectral break has been employed to: (1) constrain the age of pulsars that contribute to the excess, (2) argue that a large number of pulsars must significantly contribute to the positron flux, and (3) argue that spectral cutoffs cannot distinguish between dark matter and pulsar models. We prove that this spectral feature does not exist -- it appears due to approximations that treat inverse-Compton scattering as a continuous, instead of as a discrete and catastrophic, energy-loss process. Astrophysical sources do not produce sharp spectral features via cooling, reopening the possibility that such a feature would provide incontrovertible evidence for dark matter.

Evolution of the angles between magnetic moments and rotation axes in radio pulsars

ABSTRACT Distribution of angles β between magnetic moments and rotation axes for radio pulsars with periods in the interval 0.1 s < P < 2 s has been obtained. This distribution is rather wide with a mean value of 23°. For the bulk of pulsars, the inclination of axes is characterized by moderate values of angles β. About 60 per cent of pulsars considered have angles in the interval from 20° to 45°. The useful correlation between kinematic tkin and characteristic ages τ has been detected. The relationship between tkin and τ can be used to make a conversion of catalogue values of τ to more real pulsar ages tkin. It will be important in further pulsar investigations. It is shown that inclination angles for pulsars with 0.1 s < P < 2 s decrease with their ages. This decrease is the best described by the power law. Values of angles β have been calculated for pulsars, which are at the moment in SNRs. These angles are not changed markedly during 105 yr. Their decrease begins much later.

On the Evolution of Angles between the Magnetic Moment and the Axis of Rotation in Radio Pulsars

Three methods were used to calculate the angles β between the magnetic moment and the axis of rotation of central neutron stars for 307 known radio pulsars. There is no explicit statistical dependence of the calculated angles β on the kinematic age of pulsars. It follows from this that either the indicated angle does not change noticeably with the age of the pulsar, or pulsars are born with an arbitrary inclination angle.

The Eel Pulsar Wind Nebula: A PeVatron-candidate Origin for HAWC J1826−128 and HESS J1826−130

HAWC J1826−128 is one of the brightest Galactic TeV γ-ray sources detected by the High Altitude Water Cherenkov (HAWC) observatory, with photon energies extending up to nearly ∼100 TeV. This HAWC source spatially coincides with the H.E.S.S. TeV source HESS J1826−130 and the “Eel” pulsar wind nebula (PWN), which is associated with the GeV pulsar PSR J1826−1256. In the X-ray band, Chandra and XMM-Newton revealed that the Eel PWN is composed of both a compact nebula (∼15″) and diffuse X-ray emission (∼6′ × 2′) extending away from the pulsar. Our NuSTAR observation detected hard X-ray emission from the compact PWN up to ∼20 keV and evidence of the synchrotron burn-off effect. In addition to the spatial coincidence between HESS J1826−130 and the diffuse X-ray PWN, our multiwavelength spectral energy distribution (SED) analysis using X-ray and γ-ray data establishes a leptonic origin of the TeV emission associated with the Eel PWN. Furthermore, our evolutionary PWN SED model suggests (1) a low PWN B-field of ∼1 μG, (2) a significantly younger pulsar age (t ∼ 5.7 kyr) than the characteristic age (τ = 14.4 kyr), and (3) a maximum electron energy of PeV. The low B-field, as well as the putative supersonic motion of the pulsar, may account for the asymmetric morphology of the diffuse X-ray emission. Our results suggest that the Eel PWN may be a leptonic PeVatron particle accelerator powered by the ∼6 kyr old pulsar PSR J1826−1256 with a spin-down power of 3.6 × 1036 erg s−1.

Kinematics, structure and abundances of supernova remnant 0540-69.3

Aims. Our goal is to investigate the structure, elemental abundances, physical conditions, and the immediate surroundings of supernova remnant 0540-69.3 in the Large Magellanic Cloud. Methods. Imaging in [O III] and spectroscopic studies through various slits were carried out using European Souther Observatory’s Very Large and New Technology Telescopes. Densities, temperatures, and abundances were estimated applying nebular analysis for various parts of the remnant. Results. Several new spectral lines are identified, both from ejecta embedded in the pulsar-wind nebula, and in interstellar clouds shocked by the supernova blast wave. For the filaments in the pulsar-wind nebula, all lines are redshifted by 440 ± 80 km s−1 with respect to the rest frame of the host galaxy, and a 3D representation of the [O III] emission displays a symmetry axis of ring-like structures which could indicate that the pulsar shares the same general redshift as the central supernova ejecta. We note that [O II], [S II], [Ar III], and Hβ share a common more compact structure than [O III], and possibly [Ne III]. The average [O III] temperature for the filaments in the pulsar-wind nebula is 23 500 ± 1800 K, and the electron density derived from [S II] is typically ∼ 103 cm−3. By mass, the relative elemental abundances of the shocked ejecta in the pulsar-wind nebula are O : Ne : S : Ar ≈ 1 : 0.07 : 0.10 : 0.02, consistent with explosion models of 13 − 20 M⊙ progenitors, and similar to that of SN 1987A, as is also the explosive mixing of hydrogen and helium into the center. From Hβ and He Iλ5876, the mass ratio of He/H in the center is estimated to be in excess of ∼0.8. The rapid cooling of the shocked ejecta could potentially cause variations in the relative abundances if the ejecta are not fully microscopically mixed, and this is highlighted for S/O for the period 1989–2006. Also, [O III] is seen in presumably freely coasting photoionized ejecta outside the pulsar-wind nebula at inferred velocities out to well above 2000 km s−1, and in projection, [O III] is seen out to ∼10″ from the pulsar. This was used to estimate that the pulsar age is ≈1200 years. The freely coasting [O III]-emitting ejecta have a strictly nonspherical distribution, and their mass is estimated to be ∼0.12 M⊙. A possible outer boundary of oxygen-rich ejecta is seen in [O II] λλ3726,3729 at ∼2000 − 2100 km s−1. Four filaments of a shocked interstellar medium are identified, and there is a wide range in the degree of ionization of iron, from Fe+ to Fe13+. One filament belongs to a region also observed in X-rays, and another one has a redshift of 85 ± 30 km s−1 relative to the host. From this we estimate that the electron density of the [O III]-emitting gas is ∼ 103 cm−3, and that the line of the most highly ionized ion, [Fe XIV] λ5303, comes from an evaporation zone in connection with the radiatively cooled gas emitting, for example, [O III], and not from immediately behind the blast wave. We do not find evidence for nitrogen-enriched ejecta in the southwestern part of the remnant, as was previously suggested. Emission in this region is instead from a severely reddened H II-region.

Evidences of low-diffusion bubbles around Galactic pulsars

Recently, a few-degrees extended $\gamma$-ray halo in the direction of Geminga pulsar has been detected by HAWC, Milagro and Fermi-LAT. These observations can be interpreted with positrons ($e^+$) and electrons ($e^-$) accelerated by Geminga pulsar wind nebula (PWN), released in a Galactic environment with a low diffusion coefficient ($D_0$), and inverse Compton scattering (ICS) with the interstellar radiation fields. We inspect here how the morphology of the ICS $\gamma$-ray flux depends on the energy, the pulsar age and distance, and the strength and extension of the low-diffusion bubble. In particular we show that $\gamma$-ray experiments with a peak of sensitivity at TeV energies are the most promising ones to detect ICS halos. We perform a study of the sensitivity of HAWC, HESS and the future CTA experiment finding that, with efficiencies of the order of a few %, the first two experiments should have already detected a few tens of ICS halos while the latter will increase the number of detections by a factor of 4. We then consider a sample of sources associated to PWNe and detected in the HESS Galactic plane survey and in the second HAWC catalog. We use the information available in these catalogs for the $\gamma$-ray spatial morphology and flux of these sources to inspect the value of $D_0$ around them and the $e^{\pm}$ injection spectrum. All sources are detected as extended with a $\gamma$-ray emission extended about $15-80$ pc. Assuming that most of the $e^{\pm}$ accelerated by these sources have been released in the interstellar medium, the diffusion coefficient is $2-30 \cdot 10^{26}$ cm$^2$/s at 1 TeV, i.e. two orders of magnitude smaller than the value considered to be the average in the Galaxy. These observations imply that Galactic PWNe have low-diffusion bubbles with a size of at least 80 pc.

Discovery of a pulsar wind nebula around B0950 + 08 with the ELWA

ABSTRACT With the Expanded Long Wavelength Array and pulsar binning techniques, we searched for off-pulse emission from PSR B0950 + 08 at 76 MHz. Previous studies suggest that off-pulse emission can be due to pulsar wind nebulae (PWNe) in younger pulsars. Other studies, such as that done by Basu et al., propose that in older pulsars this emission extends to some radius that is on the order of the light cylinder radius, and is magnetospheric in origin. Through imaging analysis, we conclude that this older pulsar with a spin-down age of 17 Myr has a surrounding PWN, which is unexpected since as a pulsar ages its PWN spectrum is thought to shift from being synchrotron to inverse Compton scattering dominated. At 76 MHz, the average flux density of the off-pulse emission is 0.59 ± 0.16 Jy. The off-pulse emission from B0950+08 is ∼ 110 ± 17 arcsec (0.14 ± 0.02 pc) in size, extending well beyond the light cylinder diameter and ruling out a magnetospheric origin. Using data from our observation and the surveys VLSSr, TGSS, NVSS, FIRST, and VLASS, we have found that the spectral index for B0950 + 08 is about −1.36 ± 0.20, while the PWN’s spectral index is steeper than −1.85 ± 0.45.

The X-Ray Outburst of the Galactic Center Magnetar over Six Years of Chandra Observations

Abstract The magnetar SGR J1745−2900, discovered at a distance of parsecs from the Milky Way central black hole, Sagittarius A⋆, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745−2900 has led to interesting debates about the number, age, and nature of pulsars expected in the Galactic center region. In this work, we present extensive X-ray monitoring of the outburst of SGR J1745−2900 using the Chandra X-ray Observatory, the only instrument with the spatial resolution to distinguish the magnetar from the supermassive black hole (2.″4 angular distance). It was monitored from its outburst onset in 2013 April until 2019 August, collecting more than 50 Chandra observations for a total of more than 2.3 Ms of data. Soon after the outburst onset, the magnetar emission settled onto a purely thermal emission state that cooled from a temperature of about 0.9–0.6 keV over 6 yr. The pulsar timing properties showed at least two changes in the period derivative, increasing by a factor of about 4 during the outburst decay. We find that the long-term properties of this outburst challenge current models for the magnetar outbursts.

Optical Identification of the Millisecond Pulsar J0621+2514

AbstractUsing the SDSS and Pan-STARRS1 survey data, we found a likely companion of the recently discovered binary γ-ray radio-loud millisecond pulsar J0621+2514. Its visual brightness is about 22 mag. The broadband magnitudes and colours suggest that this is a white dwarf. Comparing the data with various white dwarfs evolutionary tracks, we found that it likely belongs to a class of He-core white dwarfs with a temperature of about 10 000 K and a mass of ≲ 0.5 M⊙. For a thin hydrogen envelope of the white dwarf, its cooling age is ≲ 0.5 Gyr which is smaller than the pulsar characteristic age of 1.8 Gyr. This may indicate that the pulsar age is overestimated. Otherwise, this may be explained by the presence of a thick hydrogen envelope or a low metallicity of the white dwarf progenitor.

Statistical Study of Glitch Behaviours of Glitching Pulsars

Detailed long-term timing observations have revealed that the expected smooth spin- down of many pulsars is prone to a variety of discrete disruptions often referred to as glitches. Although the nature and behaviour of small glitches are still poorly understood compared to large glitches, it is widely believed that both originate from some complex dynamical changes within the neutron star interior and their study could provide valuable information about the internal structure and dynamics of the neutron stars. In this paper, the distribution of glitch sizes, glitch patterns and possible relationships between glitch parameters and pulsar rotational parameters were statistically investigated using 482 glitches reported in 168 pulsars. The distribution of glitch sizes showed predominance of large glitches for J0537-6910, J0835-4510, J1341-6220 and J18001-2304; small glitches for J0534+2200, J0631+1036 and J1740-3015 and continuous glitch size distribution for J0534+2200, J1341-6220, J1740-3015 and J1801-2304. PSRs J0537-6910 and J0835-4510 showed specific regular pattern with J1740-3015 showing a quasi-regular pattern. The mean glitch size of these pulsars relates considerably with rotational frequency (ν) and spin down rate ( ) in simple power laws. Similarly, variation of glitch activity with the characteristic age (τ) traces a curve that peaks at τ = 104 yr and decays with age for older pulsars with τ ≥ 104yr. The angular momentum transfer resulting to glitches appears to be maximum at youthful age (≈ 104-106yr) of pulsars when certain rotational properties as well as temperature of the star best supports vortex pinning and unpinning of the superfluid of the star interior.

How young the accretion-powered pulsars could be?

A question about the age of accretion-powered X-ray pulsars has recently been reopened by a discovery of the X-ray pulsar SXP 1062 in the SMC. This High Mass X-ray Binary (HMXB) contains a neutron star rotating with the period of 1062 s and is associated with a supernova remnant of the age ∼ 104 yr. An attempt to explain the origin of this young long-period X-ray pulsar within the traditional scenario of three basic states (ejector, propeller and accretor) encounters difficulties. Even if this pulsar were born as a magnetar the spin-down time during the propeller stage would exceed 104 yr. Here we explore a more circuitous way of the pulsar spin evolution in HMXBs, in which the propeller stage in the evolutionary track is avoided. We find this way to be possible if the stellar wind of the massive companion to the neutron star is magnetized. The geometry of plasma flow captured by the neutron star in this case differs from spherically symmetrical and the magnetospheric radius of the neutron star is smaller than that evaluated in the convention accretion scenarios. We show that the age of an accretion-powered pulsar in this case can be as small as ∼ 104 years without the need of invoking initial magnetic field in excess of 1013 G.

The space velocities of radio pulsars

Known models proposed to explain the high space velocities of pulsars based on asymmetry of the transport coefficients of different sorts of neutrinos or electromagnetic radiation can be efficient only in the presence of high magnetic fields (to 1016 G) or short rotation periods for the neutron stars (of the order of 1 ms). This current study shows that the observed velocities are not correlated with either the pulsar periods or their surface magnetic fields. The initial rotation periods are estimated in a model for the magnetedipolar deceleration of their spin, assuming that the pulsar ages are equal to their kinematic ages. The initial period distribution is bimodal, with peaks at 5 ms and 0.5 s, and similar to the current distribution of periods. It is shown that asymmetry of the pulsar electromagnetic radiation is insufficient to give rise to additional acceleration of pulsars during their evolution after the supernova explosion that gave birth to them. The observations testify to deceleration of the motion, most likely due to the influence of the interstellar medium and interactions with nearby objects. The time scale for the exponential decrease in the magnetic field τD and in the angle between the rotation axis and magnetic moment τs are estimated, yielding τβ = 1.4 million years. The derived dependence of the transverse velocity of a pulsar on the angle between the line of sight and the rotation axis of the neutron star corresponds to the expected dependence for acceleration mechanisms associated with asymmetry of the radiation emitted by the two poles of the star.

The evolution of the magnetic inclination angle as an explanation of the long term red timing-noise of pulsars

We study the possibility that the long term red timing-noise in pulsars originates from the evolution of the magnetic inclination angle $\chi$. The braking torque under consideration is a combination of the dipole radiation and the current loss. We find that the evolution of $\chi$ can give rise to extra cubic and fourth-order polynomial terms in the timing residuals. These two terms are determined by the efficiency of the dipole radiation, the relative electric-current density in the pulsar tube and $\chi$. The following observation facts can be explained with this model: a) young pulsars have positive $\ddot{\nu}$; b) old pulsars can have both positive and negative $\ddot{\nu}$; c) the absolute values of $\ddot{\nu}$ are proportional to $-\dot{\nu}$; d) the absolute values of the braking indices are proportional to the characteristic ages of pulsars. If the evolution of $\chi$ is purely due to rotation kinematics, then it can not explain the pulsars with braking index less than 3, and thus the intrinsic change of the magnetic field is needed in this case. Comparing the model with observations, we conclude that the drift direction of $\chi$ might oscillate many times during the lifetime of a pulsar. The evolution of $\chi$ is not sufficient to explain the rotation behavior of the Crab pulsar, because the observed $\chi$ and $\dot{\chi}$ are inconsistent with the values indicated from the timing residuals using this model.

Higgs portals to pulsar collapse

Pulsars apparently missing from the galactic center could have been destroyed by asymmetric fermionic dark matter ($m_X = 1-100$ GeV) coupled to a light scalar ($m_{\phi}= 5-20$ MeV), which mixes with the Higgs boson. We point out that this pulsar-collapsing dark sector can resolve the core-cusp problem and will either be excluded or discovered by upcoming direct detection experiments. Another implication is a maximum pulsar age curve that increases with distance from the galactic center, with a normalization that depends on the couplings and masses of dark sector particles. In addition, we use old pulsars outside the galactic center to place bounds on asymmetric Higgs portal models.

TESTING DISSIPATIVE MAGNETOSPHERE MODEL LIGHT CURVES AND SPECTRA WITHFERMIPULSARS

We explore the emission properties of a dissipative pulsar magnetosphere model introduced by Kalapotharakos et al. (2014), comparing its high energy light curves and spectra, due to curvature radiation, with data collected by the Fermi LAT. The magnetosphere structure is assumed to be near the force-free solution. The accelerating electric field, inside the light-cylinder, is assumed to be negligible, while outside the light-cylinder it rescales with a finite conductivity ({\sigma}). In our approach we calculate the corresponding high energy emission by integrating the trajectories of test particles that originate from the stellar surface, taking into account both the accelerating electric field components and the radiation reaction forces. First we explore the parameter space assuming different value sets for the stellar magnetic field, stellar period, and conductivity. We show that the general properties of the model are in a good agreement with observed emission characteristics of young {\gamma}-ray pulsars, including features of the phase resolved spectra. Second we find model parameters that fit each pulsar belonging to a group of eight bright pulsars that have a published phase-resolved spectrum. The {\sigma} values that best describe each of the pulsars in this group show an increase with the spin-down rate $(\dot{E})$ and a decrease with the pulsar age, expected if pair cascades are providing the magnetospheric conductivity. Finally, we explore the limits of our analysis and suggest future directions for improving such models.

Detecting dark matter with imploding pulsars in the galactic center.

The paucity of old millisecond pulsars observed at the galactic center of the Milky Way could be the result of dark matter accumulating in and destroying neutron stars. In regions of high dark matter density, dark matter clumped in a pulsar can exceed the Schwarzschild limit and collapse into a natal black hole which destroys the pulsar. We examine what dark matter models are consistent with this hypothesis and find regions of parameter space where dark matter accumulation can significantly degrade the neutron star population within the galactic center while remaining consistent with observations of old millisecond pulsars in globular clusters and near the solar position. We identify what dark matter couplings and masses might cause a young pulsar at the galactic center to unexpectedly extinguish. Finally, we find that pulsar collapse age scales inversely with the dark matter density and linearly with the dark matter velocity dispersion. This implies that maximum pulsar age is spatially dependent on position within the dark matter halo of the Milky Way. In turn, this pulsar age spatial dependence will be dark matter model dependent.