Evolution of the Radio Remnant of Supernova 1987A: Morphological Changes from Day 7000

We present radio imaging observations of supernova remnant 1987A at 9 GHz, taken with the Australia Telescope Compact Array over 21 years from 1992 to 2013. By employing a Fourier modeling technique to fit the visibility data, we show that the remnant structure has evolved significantly since day 7000 (mid-2006): the emission latitude has gradually decreased, such that the overall geometry has become more similar to a ring structure. Around the same time, we find a decreasing trend in the east-west asymmetry of the surface emissivity. These results could reflect the increasing interaction of the forward shock with material around the circumstellar ring, and the relative weakening of the interaction with the lower-density material at higher latitudes. The morphological evolution caused an apparent break in the remnant expansion measured with a torus model, from a velocity of 4600+150-200 km/s between day 4000 and 7000 to 2400+100-200 km/s after day 7000. However, we emphasize that there is no conclusive evidence for a physical slowing of the shock at any given latitude in the expanding remnant, and that a change of radio morphology alone appears to dominate the evolution. This is supported by our ring-only fits which show a constant expansion of 3890+/-50 km/s without deceleration between days 4000 and 9000. We suggest that once the emission latitude no longer decreases, the expansion velocity obtained from the torus model should return to the same value as that measured with the ring model.Comment: 11 pages, 5 figures, accepted for publication in ApJ, Figure 1 has been scaled dow

High-resolution radio observations of SNR 1987A at high frequencies

We present new imaging observations of the remnant of Supernova (SN) 1987A at 44 GHz, performed in 2011 with the Australia Telescope Compact Array (ATCA). The 0\farcs35\times0\farcs23 resolution of the diffraction-limited image is the highest achieved to date in high-dynamic range. We also present a new ATCA image at 18 GHz derived from 2011 observations, which is super-resolved to 0\farcs25. The flux density is 40$\pm$2 mJy at 44 GHz and 81$\pm$6 mJy at 18 GHz. At both frequencies, the remnant exhibits a ring-like emission with two prominent lobes, and an east-west brightness asymmetry that peaks on the eastern lobe. A central feature of fainter emission appears at 44 GHz. A comparison with previous ATCA observations at 18 and 36 GHz highlights higher expansion velocities of the remnant eastern side. The 18-44 GHz spectral index is $\alpha=-0.80$ ($S_{\nu}\propto\nu^{\alpha}$). The spectral index map suggests slightly steeper values at the brightest sites on the eastern lobe, whereas flatter values are associated with the inner regions. The remnant morphology at 44 GHz generally matches the structure seen with contemporaneous X-ray and H$\alpha$ observations. Unlike the H$\alpha$ emission, both the radio and X-ray emission peaks on the eastern lobe. The regions of flatter spectral index align and partially overlap with the optically-visible ejecta. Simple free-free absorption models suggest that emission from a pulsar wind nebula or a compact source inside the remnant may now be detectable at high frequencies, or at low frequencies if there are holes in the ionised component of the ejecta.Comment: References updated. High resolution version may be found at http://ict.icrar.org/store/staff/gio/Papers/Zanardo_2013.pd

Constraining gamma-ray pulsar gap models with a simulated pulsar population

With the large sample of young gamma-ray pulsars discovered by the Fermi Large Area Telescope (LAT), population synthesis has become a powerful tool for comparing their collective properties with model predictions. We synthesised a pulsar population based on a radio emission model and four gamma-ray gap models (Polar Cap, Slot Gap, Outer Gap, and One Pole Caustic) normalizing to the number of detected radio pulsars in select group of surveys. The luminosity and the wide beams from the outer gaps can easily account for the number of Fermi detections in 2 years of observations. The wide slot-gap beams requires an increase by a factor of ~10 of the predicted luminosity to produce a reasonable number of gamma-ray pulsars. Such large increases in the luminosity may be accommodated by implementing offset polar caps. The narrow polar-cap beams contribute at most only a handful of LAT pulsars. Standard distributions in birth location and pulsar spin-down power (Edot) fail to reproduce the LAT findings: all models under-predict the number of LAT pulsars with high Edot, and they cannot explain the high probability of detecting both the radio and gamma-ray beams at high Edot. The beaming factor remains close to 1 over 4 decades in Edot evolution for the slot gap whereas it significantly decreases with increasing age for the outer gaps. The evolution of the slot-gap luminosity with Edot is compatible with the large dispersion of gamma-ray luminosity seen in the LAT data. The stronger evolution predicted for the outer gap, which is linked to the polar cap heating by the return current, is apparently not supported by the LAT data. The LAT sample of gamma-ray pulsars therefore provides a fresh perspective on the early evolution of the luminosity and beam width of the gamma-ray emission from young pulsars, calling for thin and more luminous gaps.Comment: 23 pages, 21 figures, accepted for publication in A&

The asymmetric radio remnant of SN 1987A

We present seven years of radio observations of SN 1987A made with the Australia Telescope Compact Array. At 1.4, 2.4, 4.8 and 8.6 GHz, the flux density of the radio remnant has increased monotonically since emission was redetected 1200 days after the explosion. On day 3200, the remnant was expanding at 2800 +/- 400 km/s, which we interpret as indicating significant deceleration of the fastest moving ejecta. Since day 1787 the spectral index has remained constant at alpha = -0.95 +/- 0.04. These observations are all consistent with the shock having encountered a denser, shocked, component of the progenitor's stellar wind. At the current rate of expansion, the shock is expected to encounter the inner optical ring in the year 2006 +/- 3. Using super-resolution, we have also obtained 9 GHz images of the remnant (resolution approx 0".5) at four epochs. The emission is distributed around the rim of a near-circular shell, but has become increasingly asymmetric with time. There are two "hotspots" to the east and west, aligned along the major axis of the optical ring. This morphology is most likely indicative of an axisymmetric circumstellar medium into which the shock is expanding, consistent with present understanding of the progenitor star and its environment.Comment: 45 pages, LaTeX, including 15 PostScript figures. To appear in "The Astrophysical Journal", volume 479 (20 Apr 1997

Application of the Gaussian mixture model in pulsar astronomy -- pulsar classification and candidates ranking for {\it Fermi} 2FGL catalog

Machine learning, algorithms to extract empirical knowledge from data, can be used to classify data, which is one of the most common tasks in observational astronomy. In this paper, we focus on Bayesian data classification algorithms using the Gaussian mixture model and show two applications in pulsar astronomy. After reviewing the Gaussian mixture model and the related Expectation-Maximization algorithm, we present a data classification method using the Neyman-Pearson test. To demonstrate the method, we apply the algorithm to two classification problems. Firstly, it is applied to the well known period-period derivative diagram, where we find that the pulsar distribution can be modeled with six Gaussian clusters, with two clusters for millisecond pulsars (recycled pulsars) and the rest for normal pulsars. From this distribution, we derive an empirical definition for millisecond pulsars as $\frac{\dot{P}}{10^{-17}} \leq3.23(\frac{P}{100 \textrm{ms}})^{-2.34}$. The two millisecond pulsar clusters may have different evolutionary origins, since the companion stars to these pulsars in the two clusters show different chemical composition. Four clusters are found for normal pulsars. Possible implications for these clusters are also discussed. Our second example is to calculate the likelihood of unidentified \textit{Fermi} point sources being pulsars and rank them accordingly. In the ranked point source list, the top 5% sources contain 50% known pulsars, the top 50% contain 99% known pulsars, and no known active galaxy (the other major population) appears in the top 6%. Such a ranked list can be used to help the future follow-up observations for finding pulsars in unidentified \textit{Fermi} point sources.Comment: 9 pages, 4 figures, accepted by MNRA

Study of the April 20, 2007 CME-Comet Interaction Event with an MHD Model

This study examines the tail disconnection event on April 20, 2007 on comet 2P/Encke, caused by a coronal mass ejection (CME) at a heliocentric distance of 0.34 AU. During their interaction, both the CME and the comet are visible with high temporal and spatial resolution by the STEREO-A spacecraft. Previously, only current sheets or shocks have been accepted as possible reasons for comet tail disconnections, so it is puzzling that the CME caused this event. The MHD simulation presented in this work reproduces the interaction process and demonstrates how the CME triggered a tail disconnection in the April 20 event. It is found that the CME disturbs the comet with a combination of a $180^\circ$ sudden rotation of the interplanetary magnetic field (IMF), followed by a $90^\circ$ gradual rotation. Such an interpretation applies our understanding of solar wind-comet interactions to determine the \textit{in situ} IMF orientation of the CME encountering Encke.Comment: 13 pages, 3 figures, accepted by the ApJ Letter

The Contribution of Fermi Gamma-Ray Pulsars to the local Flux of Cosmic-Ray Electrons and Positrons

We analyze the contribution of gamma-ray pulsars from the first Fermi-Large Area Telescope (LAT) catalogue to the local flux of cosmic-ray electrons and positrons (e+e-). We present new distance estimates for all Fermi gamma-ray pulsars, based on the measured gamma-ray flux and pulse shape. We then estimate the contribution of gamma-ray pulsars to the local e+e- flux, in the context of a simple model for the pulsar e+e- emission. We find that 10 of the Fermi pulsars potentially contribute significantly to the measured e+e- flux in the energy range between 100 GeV and 1 TeV. Of the 10 pulsars, 2 are old EGRET gamma-ray pulsars, 2 pulsars were discovered with radio ephemerides, and 6 were discovered with the Fermi pulsar blind-search campaign. We argue that known radio pulsars fall in regions of parameter space where the e+e- contribution is predicted to be typically much smaller than from those regions where Fermi-LAT pulsars exist. However, comparing the Fermi gamma-ray flux sensitivity to the regions of pulsar parameter space where a significant e+e- contribution is predicted, we find that a few known radio pulsars that have not yet been detected by Fermi can also significantly contribute to the local e+e- flux if (i) they are closer than 2 kpc, and if (ii) they have a characteristic age on the order of one mega-year.Comment: 21 pages, 6 figures, accepted for publication in JCA