Age Determination of NGC 1952 and PSR B0531+21 Using Radio Observations

Abstract

Supernova remnants (SNRs) and pulsars (PSRs) are highly related phenomena: Both can be produced from the same cosmic explosion. NGC 1952 and B0531+21, or the Crab SNR and PSR, are a verified pair. This research uses two strategies to estimate the ages of young SNRs and evaluates their reliability. The first strategy relies on the remnant’s expansion rate and radius to determine a convergence date. We operated the Green Bank Observatory 20-meter radio telescope in high-resolution mode to obtain spectra over a high-pass filter at a central frequency of 1420 MHz. While the data depicted a clear 21 cm hydrogen line, the Doppler calculation resulted in a meager expansion rate, which we found inadequate for subsequent age calculations. Second, due to the pulsar’s gradual loss of rotational energy over time, its period and period derivatives allowed us to extrapolate a formation date. We adjusted the telescope to low resolution at 1395 MHz. Age determinations reveal stages in SNR life cycles and support or disprove SNR-PSR associations. Several factors, including background radio interference, could impact the spectrum scans. However, the pulsar data produced an actual age within 50 years of the supernova’s verified explosion, confirming this method’s relative accuracy.

Keywords:Supernova RemnantsPulsarsAge Determination

Introduction

Supernova remnants (SNRs) have been astronomical phenomena of interest since ancient times due to their brilliance in the visible spectrum. SNR events are evidenced in the artifacts of early civilizations. Records are believed to have begun in 185 C.E. Common type II supernovae occur due to the discontinuation of fusion in massive stars’ iron or nickel cores. However, in type III1 supernovae, the atomic nuclei of oxygen-neon-magnesium cores capture electrons. Both processes drive stellar collapse because of imbalances between gravity and pressure. Protons and electrons meld, forming a small, dense star supported by neutron degeneracy (more massive stars continue to collapse into black holes). Finally, this process produces a shock wave that expands inward and outward–expelling an array of gas and dust in the interstellar medium. Pulsars (PSRs), a form of a neutron star, rotate at extremely high velocities. Although they conserve much of their angular momentum, their periods gradually decline. Pulsars are detectable in gamma, x-ray, ultraviolet, or radio wavelengths due to the emission of synchrotron radiation from their magnetic poles.

A PSR’s characteristic age can be deduced from its period and derivative. The characteristic age assumes significant reductions in the initial period and an absence of magnetic field decay. Its accuracy indicates the relevance of these characteristics. The actual age, meanwhile, uses a pulsar’s braking magnitude and reveals whether multiple sources affect its spin-down.

The Crab Nebula, a remnant of a supernova and a type III plerion located within the constellation Taurus, and its associated pulsar are the subjects of this research. Age determinations assist historians in connecting SNRs to Chinese records of “guest stars.” However, we focus on the Crab because of sufficient evidence from the year it was first observed, proving the accuracy of our results. Additionally, ages confirm pulsar associations: A pulsar may gradually drift from a remnant’s center Condon, 2017. We compare two methods in their effective determination of the Crab SNR’s age: The first utilizes the remnant’s expansion, and the second relies on the spin properties of its pulsar.

An optical image of the NGC 1952 taken by the Hubble Space Telescope in 2005. The central glow highlights the pulsar.

Figure 1:An optical image of the NGC 1952 taken by the Hubble Space Telescope in 2005. The central glow highlights the pulsar.

Methods

Initially, our research relied on the operation of the Pisgah Astronomical Research Institute’s 12-m radio telescope. We sought to measure the diameter of the remnant using multiple continuum scans. However, the telescope could not produce these scans due to the SNR’s uneven distribution of hydrogen gas, nor was its data acquisition system designed to detect fast pulsars.

Radio Observations

The Green Bank 20-m radio telescope with an L-band receiver (1.15-1.73 GHz) was used for all original observations.

Table 1:A summary of completed observing runs for the supernova remnant NGC 1952 and pulsar B0531+21.

RunObjectCoordinates and Apparent MagnitudeTime of Observation (UT)Date
1NGC 1952RA 05:35:54.67 DEC 22:01:49.1 Mag. 8.417:51:052023-06-27
2NGC 1952RA 05:35:54.67 DEC 22:01:49.1 Mag. 8.421:02:372023-07-01
3B0531+21RA 05:34:31.9 DEC 22:00:52.1 Mag. 16.519:00:042023-07-24
4B0531+21RA 05:34:31.9 DEC 22:00:52.1 Mag. 16.520:01:562023-07-25
5B0531+21RA 05:34:31.9 DEC 22:00:52.1 Mag. 16.519:09:292023-07-26

Telescope Settings: SNR Data

We produced two spectrum scans of the Crab SNR to estimate its expansion rate. This approach relied on measuring the Doppler shift of the hydrogen 21 cm emission line (with a precise frequency of 1420.406 MHz) and assumed that the shell had undergone no acceleration while dispersing from a point. The pulsar’s energy input, combined with the shell’s evolving interaction with the interstellar medium, would cause numerous complications Brandt & Williamson, 1979. We pulled frequency values from the minimum coordinates of the spectral features displayed by our scans and averaged the corresponding pairs. We then substituted these values into Equation (1) to find a velocity difference for the slightly redshifted and blueshifted radio waves (due to the nearer approaching side of the shell and the further receding side). Finally, with measurements for the expansion velocity, we aimed to estimate a convergence date from the radio size of the remnant.

Δv=Δvv×C\Delta v = \frac{\Delta v}{v}\times C

Telescope Settings: PSR Data

We conducted three scans of PSR B0531+21, as described in Table 2, to determine the approximate period. Green Bank’s systems produced the data in processed Prepfold Plots Figure 2. Due to the spin-down of the pulsar, the period, PP, gradually decreases over many years. The period derivative, P˙\dot{P}, describes this change. We found P˙P\dot{P}P using the Prepfolds, as well as a recorded period from the Jodrell Bank Observatory in 1982 Lyne & Tong, n.d..

Table 2:Control parameters for all telescope observations.

SettingNGC 1952B0 531+21
Receiver ModeHighresLowres
Central Frequency (MHz)14201395
FilterHIHI
PatternTrackTrack
Duration (sec)60180
Integration Time (sec)0.50.0001
Number of Channels10241024
Min. Sun Separation (degrees)120
Min. Elevation (degrees)2020

Characteristic and True Ages

We calculated the characteristic age, denoted τc\tau_c, using Equation (2), which required only the period and period derivative.

τc=P2P˙\tau_c = \frac{P}{2\dot{P}}

With Equation (3), we then estimated the true age, τ. Unlike the characteristic age, the true age of a pulsar incorporates nn, the magnetic braking index or the loss of rotational energy due to ionized material traveling along the pulsar’s magnetic field lines, and P0P_0, the initial period upon the formation of the pulsar.

τ=P(n1)P˙×[1(P0P)n1]\tau = \frac{P}{(n-1)\dot{P}} \times \left[1-\left(\frac{P_0}{P}\right)^{n-1}\right]
P=kP2nP = kP^{2-n}

A differentiated and rearranged form of Equation (4), which describes the relationship between the period, the constant, kk, and the braking index, allowed us to solve for nn: n=2PP¨PP˙P2n=2-\frac{P \ddot{P}{P}} {\dot{P}{P^2}}. The derivations of P¨P\ddot{P}{P} and P0P_0 require complex models actively being explored. The process is further complicated by numerous glitches or random and disruptive spin-up events experienced by the Crab PSR. However, recent studies have given estimates for P0P_0 Zhang et al., 2022, as well as Malov, 2017 between glitches: values which we use here.

Results

SNR Rate of Expression

The Crab remnant is composed primarily of hydrogen gas. Neutral hydrogen produces the 21 cm line, which is evident as a sharp emission feature around 1420 MHz Figure 2. The 21 cm undergoes slight dispersion in energy and quickly passes through the Earth’s atmosphere, making it observable with little interference. The averages of the spectral features’ minima, 1420.16 MHz and 1420.67 MHz, resulted in a total expansion velocity of 108 km/s.

However, the data may have been skewed by various factors. For example, it is possible that the sharp peak was produced by background sources and not the composition of the remnant itself. The Crab SNR has a relatively small angular size (4-8 arcmin) compared with the telescope’s wide field of view. We did not produce further scans due to these significant uncertainties, nor did we carry out subsequent age calculations. Yang & Chevalier (2015) estimate that the velocities of most shell features range approximately between 1260-1700 km/s. “Daisy” telescope scans process the sky surrounding an object, allowing the observer to discern radio-emitting sources within a set radius. We may implement this method in future trials to recognize sources of error.

A spectrum scan from observing run 1 displaying the linearly polarized components of the signal (XX2 and YY2), plotted as power (K) versus frequency (MHz).

Figure 2:A spectrum scan from observing run 1 displaying the linearly polarized components of the signal (XX2XX2 and YY2YY2), plotted as power (KK) versus frequency (MHz).

PSR Period, Braking, and Age

The pulse profile (top left of Figure 3) represents the period produced from folded data and reveals recurring signals. The three measurements resulted in an average period of 0.0338149930 seconds. Thus, with the February 15, 1982 record of 0.0332676584 seconds, we found the period derivative to equal 4.18551013ss14. 1855 \cdot 10^{-13} ss^{-1}.

A pulsar dataset (compiled by Green Bank into a Prepfold Plot) depicting a regular, well-defined signal from observing run 3.

Figure 3:A pulsar dataset (compiled by Green Bank into a Prepfold Plot) depicting a regular, well-defined signal from observing run 3.

Characteristic Age

We determined the characteristic age to be 4.039510104.0395 \cdot 10^{10} seconds or 1280.9 years. For millisecond pulsars, the characteristic age can be many times greater than the actual age, which indicates that P0P_0 is similar to PP. Although the approximation above is reasonable, the actual age demonstrates higher accuracy. Using Equation (4) and a period second derivative of 2.71024-2.7 \cdot 10^{-24}, nn was found to be 2.52. The braking index is slightly less than the expected value of 3 if the spin-down were due solely to magnetic dipole radiation; therefore, other factors contribute to the pulsar’s loss of rotational velocity Kou & Tong, 2015.

True Age

Utilizing Equation (3) and a value of 0.0183 seconds for P0P_0, the actual age of B0531+21 is approximately 1020 years. The Crab Supernova was first recorded in July of 1054 C.E., resulting in an actual age of 969.67 years. Lastly, an experimental error of 5% was computed, indicating a relatively accurate result. In comparison, the characteristic age produced an error of 32.09%32.09\%.

Chinese Sung Dynasty astronomers described the remnant in the Sung Shih annals, and the Japanese documented it in the texts Mei Getsuki and Ichidai, stating, “It shone like a comet and was as bright as Jupiter” Brandt & Williamson, 1979. It is also evidenced in a pictographic representation of the ruins of the Anasazi’s Peñasco Blanco.

Discussion and Conclusion

Age determinations derived from spectrum scans of the Crab SNR’s radio expansion velocity are impacted by the remnant’s competing traits and the likelihood that most of the nebular ejecta does not move directly toward or away from Earth. Considering the possibility of radio frequency interference (RFI), the expansion velocity could be measured more accurately in the visible spectrum. Published high-resolution photographs taken across several decades could also account for the velocities of individual knots. However, Nugent (1998) notes that their optical velocity closely agreed with other published research, indicating a convergence date 76 years after the 1054 C.E. supernova. Meanwhile, Bietenholz & Nugent (2015) determined that the Crab’s inner “synchrotron-emitting bubble,” detectable on the radio, expands faster than the outward optical filaments. They reached this conclusion and accounted for the remnant’s less prominent features by applying unique scaling techniques to high-resolution radio and optical photographs. The results of the spectrum scans presented here support the need for advanced observing techniques that capture the remnant’s asymmetrical geometries and dynamic physical properties.

Meanwhile, despite undergoing glitches, pulsars’ compact structures and the benefits of folded data in distinguishing recurring signals from RFI contribute to very accurate time of arrival measurements. We encounter problems considering that many young SNRs do not contain detectable pulsars – perhaps because these pulsars do not emit signals that can be discerned among the remnant’s intense radiation. However, at least 20 pairs of plerions and pulsars are known with certainty to exist. Our research conveys that pulsars provide valuable insight into the age and evolution of SNRs. In contrast to expansive nebulae and supernova ejecta, the properties of pulsars make them a very reliable tool for astronomers.

Acknowledgments

The research presented in this paper was made possible through the generous allocation of resources and equipment by the Green Bank Observatory. I thank Tim DeLisle and Melanie Crowson at the Pisgah Astronomical Research Institute for their invaluable guidance in operating radio telescopes. Additionally, I appreciate James Happer, a physics instructor at The North Carolina School of Science and Mathematics, for his assistance during the initial stages of brainstorming research topics and for providing valuable suggestions in analyzing spectrum data. I also thank Jacob Brown for his insightful advice regarding potential sources of error.

References
  1. Condon, J. J. (2017). Essential radio astronomy. Contemporary Physics, 58(3), 278–279. 10.1080/00107514.2017.1312538
  2. Brandt, J. C., & Williamson, J. C. (1979). The 1054 supernova and native American rock art (Vol. 10). Journal for the History of Astronomy, Archaeoastronomy Supplement.
  3. Lyne, A., & Tong, H. (n.d.). CRABTIME - Crab Pulsar Timing. https://heasarc.gsfc.nasa.gov/W3Browse/all/crabtime.html
  4. Zhang, C.-M., Cui, X.-H., Li, D., Wang, D.-H., Wang, S.-Q., Wang, N., Zhang, J.-W., Peng, B., Zhu, W.-W., Yang, Y.-Y., & Pan, Y.-Y. (2022). Evolution of Spin Period and Magnetic Field of the Crab Pulsar: Decay of the Braking Index by the Particle Wind Flow Torque. 10.48550/ARXIV.2212.04674
  5. Malov, I. (2017). On the second derivatives of the spin periods and braking indices in radio pulsars. Monthly Notices of the Royal Astronomical Society, 468(3), 2713–2718. 10.1093/mnras/stx619