In the past decade, the properties of annihilating dark matter models were examined by various kinds of data, including the data of gamma rays, radio waves, X-ray, positrons, electrons, antiprotons and neutrinos. In particular, most of the studies focus on the data of our Galaxy, nearby galaxies (e.g. M31 galaxy) or nearby galaxy clusters (e.g. Fornax cluster). In this article, we examine the archival radio continuum spectral data of a relatively high-redshift galaxy cluster (A697 cluster) to constrain the properties of annihilating dark matter. We find that leptophilic annihilation channels ($e^+e^-$, $\mu^+\mu^-$ and $\tau^+\tau^-$) can give very good fits to the radio continuum spectrum of the A697 cluster.

Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.

Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen to repeat. A leading model for repeating FRBs is that they are extragalactic magnetars, powered by their intense magnetic fields. However, a challenge to this model has been that FRBs must have radio luminosities many orders of magnitude larger than those seen from known Galactic magnetars. Here we report the detection of an extremely intense radio burst from the Galactic magnetar SGR 1935+2154 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project. The fluence of this two-component bright radio burst and the estimated distance to SGR 1935+2154 together imply a 400-800 MHz burst energy of $\sim 3 \times 10^{34}$ erg, which is three orders of magnitude brighter than those of any radio-emitting magnetar detected thus far. Such a burst coming from a nearby galaxy would be indistinguishable from a typical FRB. This event thus bridges a large fraction of the radio energy gap between the population of Galactic magnetars and FRBs, strongly supporting the notion that magnetars are the origin of at least some FRBs.

Using commissioning data from the first year of operation of the Canadian Hydrogen Intensity Mapping Experiment's (CHIME) Pulsar backend system, we conduct a systematic analysis of the Faraday Rotation Measure (RM) of the northern hemisphere pulsars detected by CHIME. We present 55 new RMs as well as obtain improved RM uncertainties for 25 further pulsars. CHIME's low observing frequency and wide bandwidth between 400-800 MHz contribute to the precision of our measurements, whereas the high cadence observation provide extremely high signal-to-noise co-added data. Our results represent a significant increase of the pulsar RM census, particularly regarding the northern hemisphere. These new RMs are for sources that are located in the Galactic plane out to 10 kpc, as well as off the plane to a scale height of ~16 kpc. This improved knowledge of the Faraday sky will contribute to future Galactic large-scale magnetic structure and ionosphere modelling.

We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.

We report on the detection of seven bursts from the periodically active, repeating fast radio burst (FRB) source FRB 180916.J0158+65 in the 300-400-MHz frequency range with the Green Bank Telescope (GBT). Emission in multiple bursts is visible down to the bottom of the GBT band, suggesting that the cutoff frequency (if it exists) for FRB emission is lower than 300 MHz. Observations were conducted during predicted periods of activity of the source, and had simultaneous coverage with the Low Frequency Array (LOFAR) and the FRB backend on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. We find that one of the GBT-detected bursts has potentially associated emission in the CHIME band (400-800 MHz) but we detect no bursts in the LOFAR band (110-190 MHz), placing a limit of $\alpha > -1.0$ on the spectral index of broadband emission from the source. We also find that emission from the source is severely band-limited with burst bandwidths as low as $\sim$40 MHz. In addition, we place the strictest constraint on observable scattering of the source, $<$ 1.7 ms, at 350 MHz, suggesting that the circumburst environment does not have strong scattering properties. Additionally, knowing that the circumburst environment is optically thin to free-free absorption at 300 MHz, we find evidence against the association of a hyper-compact HII region or a young supernova remnant (age $<$ 50 yr) with the source.

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

The first-year results from the Parker Solar Probe (PSP) reveal a gradual decrease of F-coronal dust from distances of $D=0.166-0.336$ AU (or the inner elongations of $\sim 9.22- 18.69~R_{\odot}$) to the Sun, which cannot be explained by the dust sublimation scenario that predicts a dust-free-zone $\lesssim 4-5R_{\odot}$. In this paper, we attempt to explain this puzzle using our newly introduced mechanism of dust destruction so-called Radiative Torque Disruption (RATD) mechanism. We demonstrate that RATD rapidly breaks large grains into nanoparticles so that they can be efficiently destroyed by nonthermal sputtering induced by bombardment of energetic protons from slow solar winds. This joint effect extends the dust-free-zone established by thermal sublimation to $R_{dfz}\sim 8R_{\odot}$, which is called a new dust-free-zone. Beyond this new dust-free-zone, we find that the dust mass decreases gradually from $R\sim 42R_{\odot}$ toward the Sun due to partial removal of nanodust by nonthermal sputtering. This feature can successfully reproduce the gradual decrease of the F-corona between $19-9R_{\odot}$ observed by the PSP. Finally, the RATD mechanism can efficiently produce nanoparticles usually observed in the inner solar system.

We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11-year data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary mass limit for our median-sensitivity pulsar within this period range is 1 M_moon (P / 100 days)^(-2/3).

A handful of fast radio bursts (FRBs) are now known to repeat. However, the question remains --- do they all? We report on an extensive observational campaign with the Australian Square Kilometre Array Pathfinder (ASKAP), Parkes, and Robert C. Byrd Green Bank Telescope, searching for repeat bursts from FRBs detected by the Commensal Real-time ASKAP Fast Transients survey. In 383.2 hr of follow-up observations covering 27 FRBs initially detected as single bursts, only two repeat bursts from a single FRB, FRB 171019, were detected, which have been previously reported by Kumar et al. We use simulations of repeating FRBs that allow for clustering in burst arrival times to calculate new estimates for the repetition rate of FRB 171019, finding only slight evidence for incompatibility with the properties of FRB 121102. Our lack of repeat bursts from the remaining FRBs set limits on the model of all bursts being attributable to repeating FRBs. Assuming a reasonable range of repetition behaviour, at most 60% (90% C.L.) of these FRBs having an intrinsic burst distribution similar to FRB~121102. This result is shown to be robust against different assumptions on the nature of repeating FRB behaviour, and indicates that if indeed all FRBs repeat, the majority must do so very rarely.

This paper presents the first detection of Shapiro delay from the binary millisecond pulsar PSR J1811-2405. We report a 11-${\sigma}$ measurement of the orthometric amplitude, h$_3$=6.8(6)x10$^{-7}$, and a 16-${\sigma}$ measurement of the orthometric ratio, ${\varsigma}$= 0.81(5). Given the relatively high orbital inclination, i= 79(2) deg, of this binary system, we obtain constraints on the companion mass of m$_c$= 0.31$\substack{+0.08 \\ -0.06}M_\odot$. The pulsar mass is currently less well constrained, with a value of 2.0$\substack{+0.8 \\ -0.5}M_\odot$. The companion mass and the orbital period are in agreement with the prediction made by previous numerical calculations of the evolution of compact binary systems. From a study of the polarization, we find that the orbital inclination angle is ~100 deg and that PSR J1811-2405 is an orthogonal rotator. In addition, the ${\mu}$s-level timing precision together with its narrow profile make PSR J1811-2405 a good candidate for inclusion in the pulsar timing arrays being used to detect nHz Gravitational waves.

We report on the discovery and analysis of bursts from nine new repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 195 to 1380 pc cm$^{-3}$. We detect two bursts from three of the new sources, three bursts from four of the new sources, four bursts from one new source, and five bursts from one new source. We determine sky coordinates of all sources with uncertainties of $\sim$10$^\prime$. We detect Faraday rotation measures for two sources, with values $-20(1)$ and $-499.8(7)$ rad m$^{-2}$, that are substantially lower than the RM derived from bursts emitted by FRB 121102. We find that the DM distribution of our events, combined with the nine other repeaters discovered by CHIME/FRB, is indistinguishable from that of thus far non-repeating CHIME/FRB events. However, as previously reported, the burst widths appear statistically significantly larger than the thus far non-repeating CHIME/FRB events, further supporting the notion of inherently different emission mechanisms and/or local environments. These results are consistent with previous work, though are now derived from 18 repeating sources discovered by CHIME/FRB during its first year of operation. We identify candidate galaxies that may contain FRB 190303.J1353+48 (DM = 222.4 pc cm$^{-3}$).

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments.

We present the results of processing an additional 44% of the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey, the most sensitive blind pulsar survey of the southern Galactic plane to date. Our partially-coherent segmented acceleration search pipeline is designed to enable the discovery of pulsars in short, highly-accelerated orbits, while our 72-min integration lengths will allow us to discover pulsars at the lower end of the pulsar luminosity distribution. We report the discovery of 40 pulsars, including three millisecond pulsar-white dwarf binary systems (PSRs J1537-5312, J1547-5709 and J1618-4624), a black-widow binary system (PSR J1745-23) and a candidate black-widow binary system (PSR J1727-2951), a glitching pulsar (PSR J1706-4434), an eclipsing binary pulsar with a 1.5-yr orbital period (PSR J1653-45), and a pair of long spin-period binary pulsars which display either nulling or intermittent behaviour (PSRs J1812-15 and J1831-04). We show that the total population of 100 pulsars discovered in the HTRU-S LowLat survey to date represents both an older and lower-luminosity population, and indicates that we have yet to reach the bottom of the luminosity distribution function. We present evaluations of the performance of our search technique and of the overall yield of the survey, considering the 94% of the survey which we have processed to date. We show that our pulsar yield falls below earlier predictions by approximately 25% (especially in the case of millisecond pulsars), and discuss explanations for this discrepancy as well as future adaptations in RFI mitigation and searching techniques which may address these shortfalls.

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

We report the detection of repeat bursts from the source of FRB 171019, one of the brightest fast radio bursts (FRBs) detected in the Australian Square Kilometre Array Pathfinder (ASKAP) fly's eye survey. Two bursts from the source were detected with the Green Bank Telescope in observations centered at 820 MHz. The repetitions are a factor of $\sim 590$ fainter than the ASKAP-discovered burst. All three bursts from this source show no evidence of scattering and have consistent pulse widths. The pulse spectra show modulation that could be evidence for either steep spectra or patchy emission. The two repetitions were the only ones found in an observing campaign for this FRB totaling 1000 hr, which also included ASKAP and the 64-m Parkes radio telescope, over a range of frequencies (720$-$2000 MHz) at epochs spanning two years. The inferred scaling of repetition rate with fluence of this source agrees with the other repeating source, FRB 121102. The detection of faint pulses from FRB 171019 shows that at least some FRBs selected from bright samples will repeat if follow-up observations are conducted with more sensitive telescopes.

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

The Canadian Hydrogen Observatory and Radio-transient Detector (CHORD) is a next-generation radio telescope, proposed for construction to start immediately. CHORD is a pan-Canadian project, designed to work with and build on the success of the Canadian Hydrogen Intensity Mapping Experiment (CHIME). It is an ultra-wideband, "large-N, small-D" telescope, consisting of a central array of 512x6-m dishes, supported by a pair of distant outrigger stations, each equipped with CHIME-like cylinders and a 64-dish array. CHORD will measure the distribution of matter over a huge swath of the Universe, detect and localize tens of thousands of Fast RadioBursts (FRBs), and undertake cutting-edge measurements of fundamental physics.

We report on the discovery of eight repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 103.5 to 1281 pc cm$^{-3}$. They display varying degrees of activity: six sources were detected twice, another three times, and one ten times. These eight repeating FRBs likely represent the bright and/or high-rate end of a distribution of infrequently repeating sources. For all sources, we determine sky coordinates with uncertainties of $\sim$10$^\prime$. FRB 180916.J0158+65 has a burst-averaged DM = $349.2 \pm 0.3$ pc cm$^{-3}$ and a low DM excess over the modelled Galactic maximum (as low as $\sim$20 pc cm$^{-3}$); this source also has a Faraday rotation measure (RM) of $-114.6 \pm 0.6$ rad m$^{-2}$, much lower than the RM measured for FRB 121102. FRB 181030.J1054+73 has the lowest DM for a repeater, $103.5 \pm 0.3$ pc cm$^{-3}$, with a DM excess of $\sim$ 70 pc cm$^{-3}$. Both sources are interesting targets for multi-wavelength follow-up due to their apparent proximity. The DM distribution of our repeater sample is statistically indistinguishable from that of the first 12 CHIME/FRB sources that have not repeated. We find, with 4$\sigma$ significance, that repeater bursts are generally wider than those of CHIME/FRB bursts that have not repeated, suggesting different emission mechanisms. Our repeater events show complex morphologies that are reminiscent of the first two discovered repeating FRBs. The repetitive behavior of these sources will enable interferometric localizations and subsequent host galaxy identifications.

We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$\alpha$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing - the process of accounting for every rotation of a pulsar over long time periods - can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system (Lorimer & Kramer 2005). One of these, called relativistic Shapiro delay (Shapiro 1964), can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-year data set with recent orbital-phase-specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be $\mathbf{2.14^{+0.10}_{-0.09}}$ solar masses (68.3% credibility interval; 95.4% credibility interval is $\mathbf{2.14^{+0.20}_{-0.18}}$ solar masses). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

Photodissociation by ultraviolet radiation is the key destruction pathway for CS in photon-dominated regions, such as diffuse clouds. However, the large uncertainties of photodissociation cross sections and rates of CS, resulting from a lack of both laboratory experiments and theoretical calculations, limit the accuracy of calculated abundances of S-bearing molecules by modern astrochemical models. Here we show a detailed \textit{ab initio} study of CS photodissociation. Accurate potential energy curves of CS electronic states were obtained by choosing an active space CAS(8,10) in MRCI+Q/aug-cc-pV(5+d)Z calculation with additional diffuse functions, with a focus on the \(B\) and \(C\,^1\Sigma^+\) states. Cross sections for both direct photodissociation and predissociation from the vibronic ground state were calculated by applying the coupled-channel method. We found that the \(C-X\) \((0-0)\) transition has extremely strong absorption due to a large transition dipole moment in the Franck-Condon region and the upper state is resonant with several triplet states via spin-orbit couplings, resulting in predissociation to the main atomic products C \((^3P)\) and S \((^1D)\). Our new calculations show the photodissociation rate under the standard interstellar radiation field is \(2.9\ee{-9}\)\,s\(^{-1}\), with a 57\% contribution from \(C-X\) \((0-0)\) transition. This value is larger than that adopted by the Leiden photodissociation and photoionization database by a factor of 3.0. Our accurate \textit{ab initio} calculations will allow more secure determination of S-bearing molecules in astrochemical models.

We present X-ray and radio monitoring observations of the gamma-ray binary PSR J2032+4127/MT91 213 during its periastron passage in late 2017. Dedicated Chandra, XMM-Newton,NuSTAR X-ray observations and VLA radio observations of this long orbit (50 years), 143 ms pulsar/Be star system clearly revealed flux and spectral variability during the passage. The X-ray spectrum hardened near periastron, with a significant decrease in the power-law photon index from \Gamma ~ 2 to 1.2 and evidence of an increased absorption column density. We identified a possible spectral break at a few keV in the spectrum that suggests synchrotron cooling. A coincident radio and X-ray flare occurred one week after periastron, which is possibly the result of the pulsar wind interacting with the Be stellar disk and generating synchrotron radiation. However, a multi-wavelength comparison indicate that the X-ray and radio spectra cannot be simply connected by a single power-law component. Hence, the emission in these two energy bands must originate from different particle populations.

DEAP-3600 is a single-phase liquid argon (LAr) direct-detection dark matter experiment, operating 2 km underground at SNOLAB (Sudbury, Canada). The detector consists of 3279 kg of LAr contained in a spherical acrylic vessel. This paper reports on the analysis of a 758 tonne\cdot day exposure taken over a period of 231 live-days during the first year of operation. No candidate signal events are observed in the WIMP-search region of interest, which results in the leading limit on the WIMP-nucleon spin-independent cross section on a LAr target of $3.9\times10^{-45}$ cm$^{2}$ ($1.5\times10^{-44}$ cm$^{2}$) for a 100 GeV/c$^{2}$ (1 TeV/c$^{2}$) WIMP mass at 90\% C. L. In addition to a detailed background model, this analysis demonstrates the best pulse-shape discrimination in LAr at threshold, employs a Bayesian photoelectron-counting technique to improve the energy resolution and discrimination efficiency, and utilizes two position reconstruction algorithms based on PMT charge and photon arrival times.

We report the detection of a single burst from the first-discovered repeating Fast Radio Burst source, FRB 121102, with CHIME/FRB, which operates in the frequency band 400-800 MHz. The detected burst occurred on 2018 November 19 and its emission extends down to at least 600 MHz, the lowest frequency detection of this source yet. The burst, detected with a significance of 23.7$\sigma$, has fluence 12$\pm$3 Jy ms and shows complex time and frequency morphology. The 34 ms width of the burst is the largest seen for this object at any frequency. We find evidence of sub-burst structure that drifts downward in frequency at a rate of -3.9$\pm$0.2 MHz ms$^{-1}$. Our best fit tentatively suggests a dispersion measure of 563.6$\pm$0.5 pc cm$^{-3}$, which is ${\approx}$1% higher than previously measured values. We set an upper limit on the scattering time at 500 MHz of 9.6 ms, which is consistent with expectations from the extrapolation from higher frequency data. We have exposure to the position of FRB 121102 for a total of 11.3 hrs within the FWHM of the synthesized beams at 600 MHz from 2018 July 25 to 2019 February 25. We estimate on the basis of this single event an average burst rate for FRB 121102 of 0.1-10 per day in the 400-800 MHz band for a median fluence threshold of 7 Jy ms in the stated time interval.

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

Radio interferometers consisting of identical antennas arranged on a regular lattice permit fast Fourier transform beamforming, which reduces the correlation cost from $\mathcal{O}(n^2)$ in the number of antennas to $\mathcal{O}(n\log n)$. We develop a formalism for describing this process and apply this formalism to derive a number of algorithms with a range of observational applications. These include algorithms for forming arbitrarily pointed tied-array beams from the regularly spaced Fourier-transform formed beams, sculpting the beams to suppress sidelobes while only losing percent-level sensitivity, and optimally estimating the position of a detected source from its observed brightness in the set of beams. We also discuss the effect that correlations in the visibility-space noise, due to cross-talk and sky contributions, have on the optimality of Fourier transform beamforming, showing that it does not strictly preserve the sky information of the $n^2$ correlation, even for an idealized array. Our results have applications to a number of upcoming interferometers, in particular the Canadian Hydrogen Intensity Mapping Experiment--Fast Radio Burst (CHIME/FRB) project.

The DEAP-3600 experiment is searching for WIMP dark matter with a 3.3 tonne single phase liquid argon (LAr) target, located 2.1 km underground at SNOLAB. The experimental signature of dark matter interactions is keV-scale $^{40}$Ar nuclear recoils (NR) producing 128 nm LAr scintillation photons observed by PMTs. The largest backgrounds in DEAP-3600 are electronic recoils (ER) induced by $\beta$ and $\gamma$-rays originating from internal and external radioactivity in the detector material. A background model of the ER interactions in DEAP-3600 was developed and is described in this work. The model is based on several components which are expected from radioisotopes in the LAr, from ex-situ material assay measurements, and from dedicated independent in-situ analyses. This prior information is used in a Bayesian fit of the ER components to a 247.2 d dataset to model the radioactivity in the surrounding detector materials. While excellent discrimination between ERs and NRs is reached with pulse shape discrimination, utilizing the large difference between fast and slow components of LAr scintillation light, detailed knowledge of the ER background and activity of detector components, sets valuable constraints on other key types of backgrounds in the detector: neutrons and alphas. In addition, the activity of $^{42}$Ar in LAr in DEAP-3600 is determined by measuring the daughter decay of $^{42}$K. This cosmogenically activated trace isotope is a relevant background at higher energies for other rare event searches using atmospheric argon e.g. DarkSide-20k, GERDA or LEGEND. The specific activity of $^{42}$Ar in the atmosphere is found to be $40.4 \pm 5.9$ $\mu$Bq/kg of argon.

Through high-precision radio timing observations, we show that five recycled pulsars in the direction of the Galactic Centre (GC) have anomalous spin period time derivative ($\dot P$) measurements -- PSRs J1748$-$3009, J1753$-$2819, J1757$-$2745, and J1804$-$2858 show negative values of $\dot P$ and PSR J1801$-$3210 is found to have an exceptionally small value of $\dot P$. We attribute these observed $\dot P$ measurements to acceleration of these pulsars along their lines-of-sight (LOSs) due to the Galactic gravitational field. Using models of the Galactic mass distribution and pulsar velocities, we constrain the distances to these pulsars, placing them on the far-side of the Galaxy, providing the first accurate distance measurements to pulsars located in this region and allowing us to consider the electron density along these LOSs. We find the new electron density model YMW16 to be more consistent with these observations than the previous model NE2001. The LOS dynamics further constrain the model-dependent intrinsic $\dot P$ values for these pulsars and they are consistent with measurements for other known pulsars. In the future, the independent distance measurements to these and other pulsars near the GC would allow us to constrain the Galactic gravitational potential more accurately.

Magnetars are strongly magnetized pulsars and they occasionally show violent radiative outbursts. They also often exhibit glitches which are sudden changes in the spin frequency. It was found that some glitches were associated with outbursts but their connection remains unclear. We present a systematic study to identify possible correlations between them. We find that the glitch size of magnetars likely shows a bimodal distribution, different from the distribution of the Vela-like recurrent glitches but consistent with the high end of that of normal pulsars. A glitch is likely a necessary condition for an outburst but not a sufficient condition because only 30\% of glitches were associated with outbursts. In the outburst cases, the glitches tend to induce larger frequency changes compared to those unassociated ones. We argue that a larger glitch is more likely to trigger the outburst mechanism, either by reconfiguration of the magnetosphere or deformation of the crust. A more frequent and deeper monitoring of magnetars is necessary for further investigation of their connection.

Magnetars are neutron stars with extremely high surface magnetic fields. They show diverse X-ray pulse profiles in the quiescent state. We perform a systematic Fourier analysis of their soft X-ray pulse profiles. We find that most magnetars have a single-peaked profile and hence have low amplitudes of the second Fourier harmonics (A2). On the other hand, the pulsed fraction (PF) spreads over a wide range. We compared the results with theoretical profiles assuming various surface hotspot asymmetries, viewing geometries, and beaming functions. We found that a single value of the intensity ratio r between two antipodal hotspots is unable to reproduce the observed distribution of A2 and PF for all magnetars. The inferred r is probably anticorrelated with the thermal luminosity, implying that high-luminosity magnetars tend to have two symmetric hotspots. Our results are consistent with theoretical predictions, for which the existence of an evolving toroidal magnetic field breaks the symmetry of the surface temperature.

X-ray photometry and optical spectra are presented covering the periastron passage of the highly-eccentric, ~50 year binary system PSR J2032+4127 in November 2017. This system consists of a 143 ms pulsar in orbit around a massive OB star, MT 91-213. The data show dramatic changes during the encounter as the pulsar wind collided with the stellar wind. The X-ray flux rose on the approach to periastron, then underwent a major dip in the few days around periastron, and then gradually declined over the next few weeks. The optical spectroscopy revealed a steady decline in the H? line strength on the approach to periastron (from an Equivalent Width of -15A to -7A) implying a truncation of the OB star's circumstellar disk by the approaching neutron star. Smooth Particle Hydrodynamic (SPH) modelling is used here to model the system within the context of the observed behaviour and predict the geometrical configuration of the circumstellar disk with respect to the pulsar's orbit.

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years.

We report on the latest six pulsars discovered through our standard pipeline in the intermediate-latitude region (|b| < 15 deg) of the Parkes High Time Resolution Universe Survey (HTRU). We also present timing solutions for the new discoveries and for 19 further pulsars for which only discovery parameters were previously published. Highlights of the presented sample include the isolated millisecond pulsar J1826-2415, the long-period binary pulsar J1837-0822 in a mildly eccentric 98-day orbit with a > 0.27 M_sun companion, and the nulling pulsar J1638-4233, detected only 10% of the time. Other interesting objects are PSR J1757-1500, exhibiting sporadic mode changes, and PSR J1635-2616 showing one glitch over 6 years. The new discoveries bring the total count of HTRU intermediate-latitude pulsars to 113, 25% of which are recycled pulsars. This is the higest ratio of recycled over ordinary pulsars discoveries of all recent pulsar surveys in this region of the sky. Among HTRU recycled pulsars, four are isolated objects. Comparing the characteristics of Galactic fully-recycled isolated MSPs with those of eclipsing binaries ('spiders'), from which the former are believed to have formed, we highlight a discrepancy in their spatial distribution. This may reflect a difference in the natal kick, hence, possibly, a different formation path. On the other hand, however, isolated fully-recycled MSPs spin periods are, on average, longer than those of spiders, in line with what one would expect, from simple magnetic-dipole spin-down, if the former were indeed evolved from the latter.

The discovery of a repeating Fast Radio Burst (FRB) source, FRB 121102, eliminated models involving cataclysmic events for this source. No other repeating FRB has yet been detected in spite of many recent FRB discoveries and follow-ups, suggesting repeaters may be rare in the FRB population. Here we report the detection of six repeat bursts from FRB 180814.J0422+73, one of the 13 FRBs detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project during its pre-commissioning phase in July and August 2018. These repeat bursts are consistent with originating from a single position on the sky, with the same dispersion measure (DM), ~189 pc cm-3. This DM is approximately twice the expected Milky Way column density, and implies an upper limit on the source redshift of 0.1, at least a factor of ~2 closer than FRB 121102. In some of the repeat bursts, we observe sub-pulse frequency structure, drifting, and spectral variation reminiscent of that seen in FRB 121102, suggesting similar emission mechanisms and/or propagation effects. This second repeater, found among the first few CHIME/FRB discoveries, suggests that there exists -- and that CHIME/FRB and other wide-field, sensitive radio telescopes will find -- a substantial population of repeating FRBs.

Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio flashes likely arriving from far outside the Milky Way galaxy. This phenomenon was discovered at radio frequencies near 1.4 GHz and to date has been observed in one case at as high as 8 GHz, but not below 700 MHz in spite of significant searches at low frequencies. Here we report detections of FRBs at radio frequencies as low as 400 MHz, on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) using the CHIME/FRB instrument. We present 13 FRBs detected during a telescope pre-commissioning phase, when our sensitivity and field-of-view were not yet at design specifications. Emission in multiple events is seen down to 400 MHz, the lowest radio frequency to which we are sensitive. The FRBs show a variety of temporal scattering behaviours, with the majority significantly scattered, and some apparently unscattered to within measurement uncertainty even at our lowest frequencies. Of the 13 reported here, one event has the lowest dispersion measure yet reported, implying it is among the closest yet known, and another has shown multiple repeat bursts, as described in a companion paper. Our low-scattering events suggest that efforts to detect FRBs at radio frequencies below 400 MHz will eventually be successful. The overall scattering properties of our sample suggest that FRBs as a class are preferentially located in environments that scatter radio waves more strongly than the diffuse interstellar medium (ISM) in the Milky Way.

We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres.

Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21 and J2145-0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643-1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145-0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643-1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible.

We have performed a new search for radio pulsars in archival data of the intermediate and high Galactic latitude parts of the Southern High Time Resolution Universe pulsar survey. This is the first time the entire dataset has been searched for binary pulsars, an achievement enabled by GPU-accelerated dedispersion and periodicity search codes nearly 50 times faster than the previously used pipeline. Candidate selection was handled entirely by a Machine Learning algorithm, allowing for the assessment of 17.6 million candidates in a few person-days. We have also introduced an outlier detection algorithm for efficient radio-frequency interference (RFI) mitigation on folded data, a new approach that enabled the discovery of pulsars previously masked by RFI. We discuss implications for future searches, particularly the importance of expanding work on RFI mitigation to improve survey completeness. In total we discovered 23 previously unknown sources, including 6 millisecond pulsars and at least 4 pulsars in binary systems. We also found an elusive but credible redback candidate that we have yet to confirm.

Fast radio bursts (FRBs) are millisecond pulses of radio emission of seemingly extragalactic origin. More than 50 FRBs have now been detected, with only one seen to repeat. Here we present a new FRB discovery, FRB 110214, which was detected in the high latitude portion of the High Time Resolution Universe South survey at the Parkes telescope. FRB 110214 has one of the lowest dispersion measures of any known FRB (DM = 168.9$\pm$0.5 pc cm$^{-3}$), and was detected in two beams of the Parkes multi-beam receiver. A triangulation of the burst origin on the sky identified three possible regions in the beam pattern where it may have originated, all in sidelobes of the primary detection beam. Depending on the true location of the burst the intrinsic fluence is estimated to fall in the range of 50 -- 2000 Jy ms, making FRB 110214 one of the highest-fluence FRBs detected with the Parkes telescope. No repeating pulses were seen in almost 100 hours of follow-up observations with the Parkes telescope down to a limiting fluence of 0.3 Jy ms for a 2-ms pulse. Similar low-DM, ultra-bright FRBs may be detected in telescope sidelobes in the future, making careful modeling of multi-beam instrument beam patterns of utmost importance for upcoming FRB surveys.

Low-frequency gravitational-wave experiments require the highest timing precision from an array of the most stable millisecond pulsars. Several known sources of noise on short timescales in single radio-pulsar observations are well described by a simple model of three components: template-fitting from a finite signal-to-noise ratio, pulse phase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellar scattering variations. Currently template-fitting errors dominate, but as radio telescopes push towards higher signal-to-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics of jitter becomes crucial for properly characterizing arrival-time uncertainties. We characterize the radio-frequency dependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for its frequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations of rms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter is typically ~1% of pulse phase. The average frequency dependence for all pulsars using a power-law model has index -0.42. We investigate the jitter variations for the interpulse of PSR~B1937+21 and find no significant deviations from the main pulse rms jitter. We also test the time-variation of jitter in two pulsars and find that systematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly model jitter as a significant component of the noise within the detector.

We report on updated radio imaging observations of the radio remnant of Supernova 1987A (SN 1987A) at 9 GHz, taken with the Australia Telescope Compact Array (ATCA), covering a 25-year period (1992-2017). We use Fourier modeling of the supernova remnant to model its morphology, using both a torus model and a ring model, and find both models show an increasing flux density, and have shown a continuing expansion of the remnant. As found in previous studies, we find the torus model most accurately fits our data, and has shown a change in the remnant expansion at Day 9,300 $\pm$210 from 2,300 $\pm$200 km/s to 3,610 $\pm$240 km/s. We have also seen an increase in brightness in the western lobe of the remnant, although the eastern lobe is still the dominant source of emission, unlike what has been observed at contemporary optical and X-ray wavelengths. We expect to observe a reversal in this asymmetry by the year $\sim$2020, and note the south-eastern side of the remnant is now beginning to fade, as has also been seen in optical and X-ray data. Our data indicate that high-latitude emission has been present in the remnant from the earliest stages of the shockwave interacting with the equatorial ring around Day 5,000. However, we find the emission has become increasingly dominated by the low-lying regions by Day 9,300, overlapping with the regions of X-ray emission. We conclude that the shockwave is now leaving the equatorial ring, exiting first from the south-east region of the remnant, and is re-accelerating as it begins to interact with the circumstellar medium beyond the dense inner ring.

PSR J1119-6127 is a radio pulsar that behaved magnetar like bursts, and we performed a comprehensive investigation of this pulsar using the archival high-energy observations obtained after its outburst in 2016 July. After the 2016 outburst, specific regions on the neutron star surface were heated up to > 0.3 and > 1 keV from about 0.2 keV. A hard non-thermal spectral component with a photon index < 0.5 related to the magnetospheric emission can be resolved from the NuSTAR spectra above 10 keV. We find that the thermal emitting regions did not cool down and gradually shrank by about (20-35)% four months after the outburst. Hard X-ray pulsations were detected with NuSTAR immediately after the outburst, at 5 sigma confidence level and with a background-subtracted pulsed fraction of (40+-10)%. However, the signal became undetectable after a few days. Using Fermi data, we found that the gamma-ray emission in 0.5-300 GeV was suppressed along with the disappearance of the radio pulsations. This is likely caused by a reconfiguration of the magnetic field. We also discovered that the timing noise evolved dramatically, and the spin-down rate significantly increased after the 2016 glitch. We proposed that post-outburst temporal and spectral behaviors from radio to gamma-ray bands were caused by changes of the magnetosphere structure, pair plasma injection and the shrinking emission sites on the neutron star.

Transient ultraluminous X-ray sources (ULXs) provide an important link bridging transient low-mass X-ray binaries and ULXs. Here we report the first discovery of both a canonical sub-Eddington outburst and an ultraluminous super-Eddington outburst from an unusual transient ULX, NGC 7793 P9 with a variability factor higher than $10^3$. Its X-ray spectrum switches between the typical high/soft state and the steep power-law state during the canonical outburst. The inner radius of the accretion disk and the disk temperature--luminosity correlation suggest that P9 harbors a stellar-mass black hole (BH). At the beginning of the ultraluminous outburst, we observe a cool outer disk with a hard Comptonized spectrum, implying a transition to the ULX regime. When the luminosity increases to $L\gtrsim3\times10^{39}$ erg/s, P9 shows a significantly curved spectrum that can be described by either a slim disk or a strongly curved Comptonized corona. The phenomenological model suggests that the hot disk observed near the peak of the ultraluminous outburst is coincidentally consistent with the extension of the thermal track. Utilizing more physical Comptonized disk models, we suggest that the corona cools down and the apparent disk-like spectrum is dominated by soft Comptonization. The significant variability above 1 keV supports this two-component scenario. The spectral evolution can also be interpreted with the supercritical accretion model. All these indicate that a canonical black hole X-ray binary can show properties of a ULX.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind which produces correlated DM fluctuations in all pulsars. We develop and apply techniques for extracting the imprint of the solar wind from the full collection of DM measurements in the recently released NANOGrav 11-yr data set. We filter out long time scale DM fluctuations attributable to structure in the interstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate the correlated signature of the wind from signals unique to individual lines of sight. When treating the solar wind as spherically symmetric and constant in time, we find the electron number density at 1 A.U. to be $7.9\pm0.2$ cm$^{-3}$. Our data shows little evidence of long-term variation in the density of the wind. We argue that our techniques paired with a high cadence, low radio frequency observing campaign of near-ecliptic pulsars would be capable of mapping out large-scale latitudinal structure in the wind.

This paper reports the first results of a direct dark matter search with the DEAP-3600 single-phase liquid argon (LAr) detector. The experiment was performed 2 km underground at SNOLAB (Sudbury, Canada) utilizing a large target mass, with the LAr target contained in a spherical acrylic vessel of 3600 kg capacity. The LAr is viewed by an array of PMTs, which would register scintillation light produced by rare nuclear recoil signals induced by dark matter particle scattering. An analysis of 4.44 live days (fiducial exposure of 9.87 tonne-days) of data taken with the nearly full detector during the initial filling phase demonstrates the detector performance and the best electronic recoil rejection using pulse-shape discrimination in argon, with leakage $<1.2\times 10^{-7}$ (90% C.L.) between 16 and 33 keV$_{ee}$. No candidate signal events are observed, which results in the leading limit on WIMP-nucleon spin-independent cross section on argon, $<1.2\times 10^{-44}$ cm$^2$ for a 100 GeV/c$^2$ WIMP mass (90% C.L.).

Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun ($6\times 10^{-29}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in optical if our line-of-sight is free of ejecta dust, and $10^{36}$ erg s$^{-1}$ ($2\times 10^{-30}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in 2-10 keV X-rays. Our X-ray limits are an order of magnitude less constraining than previous limits because we use a more realistic ejecta absorption model based on three-dimensional neutrino-driven SN explosion models (presented in an accompanying article). The allowed bolometric luminosity of the compact object is 22 Lsun if our line-of-sight is free of ejecta dust, or 138 Lsun if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star to <4-8 MK and do not exclude models, which typically are in the range 3-4 MK. For the simplest accretion model, the accretion rate for an efficiency $\eta$ is limited to $< 10^{-11} \eta^{-1}$ Msun yr$^{-1}$, which excludes most predictions. For pulsar activity modeled by a rotating magnetic dipole in vacuum, the limit on the magnetic field strength ($B$) for a given spin period ($P$) is $B < 10^{14} P^2$ G s$^{-2}$. By combining information about radiation reprocessing and geometry, it is likely that the compact object is a dust-obscured thermally-emitting neutron star, which may appear as a region of higher-temperature ejecta dust emission.

The frequency dependence of radio pulse arrival times provides a probe of structures in the intervening media. Demorest et al. 2013 was the first to show a short-term (~100-200 days) reduction in the electron content along the line of sight to PSR J1713+0747 in data from 2008 (approximately MJD 54750) based on an apparent dip in the dispersion measure of the pulsar. We report on a similar event in 2016 (approximately MJD 57510), with average residual pulse-arrival times of approximately 3.0,-1.3, and -0.7 microseconds at 820, 1400, and 2300 MHz, respectively. Timing analyses indicate possible departures from the standard nu^-2 dispersive-delay dependence. We discuss and rule out a wide variety of potential interpretations. We find the likeliest scenario to be lensing of the radio emission by some structure in the interstellar medium, which causes multiple frequency-dependent pulse arrival-time delays.

We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSRs B1937+21, J1713+0747, and J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability, and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.

Supernova 1987A in the Large Magellanic Cloud has proven a unique laboratory to investigate particle acceleration in young supernova remnants. Here we report the first detection of linear polarization of the supernova's synchrotron emission from imaging observations at frequencies spanning from 20 to 50 GHz, carried out with the Australia Telescope Compact Array between October 2015 and May 2016. The direction of the radio polarization, corrected for Faraday rotation, points to a primarily radial magnetic field across the inner ring, encompassing both the reverse and forward shocks. The magnetic field strength peaks over the high-emissivity eastern sites, where efficient cosmic ray acceleration likely takes place under quasi-parallel shocks at high Mach numbers. The mean fraction of polarized emission in the brightest sites is $2.7\pm0.2\%$ at 22 GHz and $3.5\pm0.7\%$ at 44 GHz. In the inner remnant, non-radial components of the polarized emission appear to be more prevalent. However, the low significance detection in the central regions limits interpretation.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

The young energetic pulsar J1747-2958 ($\tau=26$ kyr, $\dot{E}=2.5\times 10^{36}$ erg s$^{-1}$) powers the Mouse pulsar wind nebula (PWN), famous for its spectacular tail spanning 45$"$ in X-rays and 12$'$ in radio ($d\sim5$ kpc). We present the results of Chandra observations of the PWN and the analysis of archival lower-frequency data. The Chandra HRC image reveals a point-like source at the pulsar position, $\approx$1$'$ behind the bow shock apex of the PWN. The flattened appearance of the compact nebula is consistent with an equatorial outflow deformed by the ram pressure, implying that the angle between the pulsar's spin axis and line of sight is $\sim70^\circ$ (in agreement with the radio and $\gamma$-ray pulse profiles). The spatially-resolved spectroscopy with Chandra ACIS shows that the power-law (PL) spectrum steepens from $\Gamma=1.65\pm0.02$ to $3.0\pm0.1$ over the 45$'$ extent of the X-ray tail. However, the tail's X-ray spectrum integrated over its 45$"$ length fits a single absorbed PL with $\Gamma=2.09\pm0.03$. We also found the Mouse PWN in 150 MHz GMRT data, and a possible counterpart in 24 ${\rm \mu}$m Spitzer data. The multiwavelength data suggest that, at low frequencies, the spectrum of the X-ray-emitting tail region can be described by a broken PL with at least one turnover between radio and X-rays. This is consistent with synchrotron cooling of electrons injected at the termination shock (with an SED slope of 2.2) in an equipartition magnetic field $B\sim200$ ${\rm \mu G}$ and a bulk flow speed $v\sim4000$ km s$^{-1}$.

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.

The Dark matter Experiment using Argon Pulse-shape discrimination (DEAP) has been designed for a direct detection search for particle dark matter using a single-phase liquid argon target. The projected cross section sensitivity for DEAP-3600 to the spin-independent scattering of Weakly Interacting Massive Particles (WIMPs) on nucleons is $10^{-46}~\rm{cm}^{2}$ for a 100 GeV/$c^2$ WIMP mass with a fiducial exposure of 3 tonne-years. This paper describes the physical properties and construction of the DEAP-3600 detector.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a novel transit radio telescope operating across the 400-800-MHz band. CHIME is comprised of four 20-m x 100-m semi-cylindrical paraboloid reflectors, each of which has 256 dual-polarization feeds suspended along its axis, giving it a >200 square degree field-of-view. This, combined with wide bandwidth, high sensitivity, and a powerful correlator makes CHIME an excellent instrument for the detection of Fast Radio Bursts (FRBs). The CHIME Fast Radio Burst Project (CHIME/FRB) will search beam-formed, high time-and frequency-resolution data in real time for FRBs in the CHIME field-of-view. Here we describe the CHIME/FRB backend, including the real-time FRB search and detection software pipeline as well as the planned offline analyses. We estimate a CHIME/FRB detection rate of 2-42 FRBs/sky/day normalizing to the rate estimated at 1.4-GHz by Vander Wiel et al. (2016). Likely science outcomes of CHIME/FRB are also discussed. CHIME/FRB is currently operational in a commissioning phase, with science operations expected to commence in the latter half of 2018.

We present the Neil Gehrels Swift Observatory (Swift), Fermi Large Area Telescope (Fermi-LAT), and Karl G. Jansky Very Large Array (VLA) observations of the gamma-ray binary PSR J2032+4127/MT91 213, of which the periastron passage has just occurred in November 2017. In the Swift X-ray light curve, the flux was steadily increasing before mid-October 2017, however, a sharp X-ray dip on a weekly time-scale is seen during the periastron passage, followed by a post-periastron X-ray flare lasting for ~20 days. We suggest that the X-ray dip is caused by (i) an increase of the magnetization parameter at the shock, and (ii) the suppression due to the Doppler boosting effect. The 20-day post-periastron flare could be a consequence of the Be stellar disk passage by the pulsar. An orbital GeV modulation is also expected in our model, however, no significant variability is seen in the Fermi-LAT light curve. We suspect that the GeV emission resulted from the interaction between the binary's members is hidden behind the bright magnetospheric emission of the pulsar. Pulsar gating technique would be useful to remove the magnetospheric emission and recover the predicted GeV modulation, if an accurate radio timing solution over the periastron passage is provided in the future.

Radio pulsars found in binary systems with short orbital periods are usually fast spinning as a consequence of recycling via mass transfer from their companion stars; this process is also thought to decrease the magnetic field of the neutron star being recycled. Here, we report on timing observations of the recently discovered binary PSR J1755$-$2550 and find that this pulsar is an exception: with a characteristic age of 2.1 Myr, it is relatively young; furthermore, with a spin period of 315 ms and a surface magnetic field strength at its poles of 0.88$\times$10$^{12}$ G the pulsar shows no sign of having been recycled. Based on its timing and orbital characteristics, the pulsar either has a massive white dwarf (WD) or a neutron star (NS) companion. To distinguish between these two cases, we searched radio observations for a potential recycled pulsar companion and analysed archival optical data for a potential WD companion. Neither work returned conclusive detections. We apply population synthesis modelling and find that both solutions are roughly equally probable. Our population synthesis also predicts a minimum mass of 0.90 M$_{\odot}$ for the companion star to PSR J1755$-$2550 and we simulate the systemic runaway velocities for the resulting WDNS systems which may merge and possibly produce Ca-rich supernovae. Whether PSR J1755$-$2550 hosts a WD or a NS companion star, it is certainly a member of a rare subpopulation of binary radio pulsars.

We have investigated the mode-changing properties of PSR B0329+54 using 31 epochs of simultaneous 13 cm/3 cm single-pulse observations obtained with Shanghai Tian Ma 65 m telescope. The pulsar was found in the abnormal emission mode 17 times, accounting for ~13% of the 41.6 hours total observation time. Single pulse analyses indicate that mode changes took place simultaneously at 13 cm/3 cm within a few rotational periods. We detected occasional bright and narrow pulses whose peak flux densities were 10 times higher than that of the integrated profile in both bands. At 3 cm, about 0.66% and 0.27% of single pulses were bright in the normal mode and abnormal mode respectively, but at 13 cm the occurrence rate was only about 0.007%. We divided the pulsar radiation window into three components (C1, C2 and C3) corresponding to the main peaks of the integrated profile. The bright pulses preferentially occurred at pulse phases corresponding to the peaks of C2 and C3. Fluctuation spectra showed that C2 had excess red noise in the normal mode, but broad quasi-periodic features with central frequencies around 0.12 cycles/period in the abnormal mode. At 3 cm, C3 had a stronger quasi-periodic modulation centered around 0.06 cycles/period in the abnormal mode. Although there were some asymmetries in the two-dimensional fluctuation spectra, we found no clear evidence for systematic subpulse drifting. Consistent with previous low-frequency observations, we found a very low nulling probability for B0329+54 with upper limits of 0.13% and 1.68% at 13 cm/3 cm respectively.

Pulsars in relativistic binary systems have emerged as fantastic natural laboratories for testing theories of gravity, the most prominent example being the double pulsar, PSR J0737$-$3039. The HTRU-South Low Latitude pulsar survey represents one of the most sensitive blind pulsar surveys taken of the southern Galactic plane to date, and its primary aim has been the discovery of new relativistic binary pulsars. Here we present our binary pulsar searching strategy and report on the survey's flagship discovery, PSR J1757$-$1854. A 21.5-ms pulsar in a relativistic binary with an orbital period of 4.4 hours and an eccentricity of 0.61, this double neutron star (DNS) system is the most accelerated pulsar binary known, and probes a relativistic parameter space not yet explored by previous pulsar binaries.

We report the discovery of PSR J1757$-$1854, a 21.5-ms pulsar in a highly-eccentric, 4.4-h orbit around a neutron star (NS) companion. PSR J1757$-$1854 exhibits some of the most extreme relativistic parameters of any known pulsar, including the strongest relativistic effects due to gravitational-wave (GW) damping, with a merger time of 76 Myr. Following a 1.6-yr timing campaign, we have measured five post-Keplerian (PK) parameters, yielding the two component masses ($m_\text{p}=1.3384(9)\,\text{M}_\odot$ and $m_\text{c}=1.3946(9)\,\text{M}_\odot$) plus three tests of general relativity (GR), which the theory passes. The larger mass of the NS companion provides important clues regarding the binary formation of PSR J1757$-$1854. With simulations suggesting 3-$\sigma$ measurements of both the contribution of Lense-Thirring precession to the rate of change of the semi-major axis and the relativistic deformation of the orbit within $\sim7-9$ years, PSR J1757$-$1854 stands out as a unique laboratory for new tests of gravitational theories.

Magnetars are neutron stars having extreme magnetic field strengths. Study of their emission properties in quiescent state can help understand effects of a strong magnetic field on neutron stars. SGR 0501+4516 is a magnetar that was discovered in 2008 during an outburst, which has recently returned to quiescence. We report its spectral and timing properties measured with new and archival observations from the Chandra X-Ray Observatory, XMM-Newton, and Suzaku. We found that the quiescent spectrum is best fit by a power-law plus two blackbody model, with temperatures of kT_low~0.26 keV and kT_high~0.62 keV. We interpret these two blackbody components as emission from a hotspot and the entire surface. The hotspot radius shrunk from 1.4 km to 0.49 km since the outburst, and there was a significant correlation between its area and the X-ray luminosity, which agrees well with the prediction by the twisted magnetosphere model. We applied the two temperature spectral model to all magnetars in quiescence and found that it could be a common feature among the population. Moreover, the temperature of the cooler blackbody shows a general trend with the magnetar field strength, which supports the simple scenario of heating by magnetic field decay.

We present the discovery of a binary millisecond pulsar (MSP), PSR J2322$-$2650, found in the Southern section of the High Time Resolution Universe survey. This system contains a 3.5-ms pulsar with a $\sim10^{-3}$ M$_{\odot}$ companion in a 7.75-hour circular orbit. Follow-up observations at the Parkes and Lovell telescopes have led to precise measurements of the astrometric and spin parameters, including the period derivative, timing parallax, and proper motion. PSR J2322$-$2650 has a parallax of $4.4\pm1.2$ mas, and is thus at an inferred distance of $230^{+90}_{-50}$ pc, making this system a candidate for optical studies. We have detected a source of $R\approx26.4$ mag at the radio position in a single $R$-band observation with the Keck Telescope, and this is consistent with the blackbody temperature we would expect from the companion if it fills its Roche lobe. The intrinsic period derivative of PSR J2322$-$2650 is among the lowest known, $4.4(4)\times10^{-22}$ s s$^{-1}$, implying a low surface magnetic field strength, $4.0(4)\times10^7$ G. Its mean radio flux density of 160 $\mu$Jy combined with the distance implies that its radio luminosity is the lowest ever measured, $0.008(5)$ mJy kpc$^2$. The inferred population of these systems in the Galaxy may be very significant, suggesting that this is a common MSP evolutionary path.

The CHIME telescope (the Canadian Hydrogen Intensity Mapping Experiment) recently built in Penticton, Canada, is currently being commissioned. Originally designed as a cosmology experiment, it was soon recognized that CHIME has the potential to simultaneously serve as an incredibly useful radio telescope for pulsar science. CHIME operates across a wide bandwidth of 400-800 MHz and will have a collecting area and sensitivity comparable to that of the 100-m class radio telescopes. CHIME has a huge field of view of ~250 square degrees. It will be capable of observing 10 pulsars simultaneously, 24-hours per day, every day, while still accomplishing its missions to study Baryon Acoustic Oscillations and Fast Radio Bursts. It will carry out daily monitoring of roughly half of all pulsars in the northern hemisphere, including all NANOGrav pulsars employed in the Pulsar Timing Array project. It will cycle through all pulsars in the northern hemisphere with a range of cadence of no more than 10 days.

We describe the Survey for Pulsars and Extragalactic Radio Bursts (SUPERB), an ongoing pulsar and fast transient survey using the Parkes radio telescope. SUPERB involves real-time acceleration searches for pulsars and single-pulse searches for pulsars and fast radio bursts. We report on the observational setup, data analysis, multi-wavelength/messenger connections, survey sensitivities to pulsars and fast radio bursts and the impact of radio frequency interference. We further report on the first 10 pulsars discovered in the project. Among these is PSR~J1306$-$40, a millisecond pulsar in a binary system where it appears to be eclipsed for a large fraction of the orbit. PSR~J1421$-$4407 is another binary millisecond pulsar; its orbital period is $30.7$ days. This orbital period is in a range where only highly eccentric binaries are known, and expected by theory; despite this its orbit has an eccentricity of $10^{-5}$.

Supernova remnants (SNRs) are powerful particle accelerators. As a supernova (SN) blast wave propagates through the circumstellar medium (CSM), electrons and protons scatter across the shock and gain energy by entrapment in the magnetic field. The accelerated particles generate further magnetic field fluctuations and local amplification, leading to cosmic ray production. The wealth of data from Supernova 1987A is providing a template of the SN-CSM interaction, and an important guide to the radio detection and identification of core-collapse SNe based on their spectral properties. Thirty years after the explosion, radio observations of SNR 1987A span from 70 MHz to 700 GHz. We review extensive observing campaigns with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA), and follow-ups with other radio telescopes. Observations across the radio spectrum indicate rapid changes in the remnant morphology, while current ATCA and ALMA observations show that the SNR has entered a new evolutionary phase.

We present an X-ray imaging and spectroscopic study of the molecular cloud interacting mixed-morphology (MM) supernova remnant (SNR) G346.6-0.2 using XMM-Newton. The X-ray spectrum of the remnant is well described by a recombining plasma that most likely arises from adiabatic cooling, and has sub-solar abundances of Mg, Si, and S. Our fits also suggest the presence of either an additional power-law component with a photon index of $\sim$2, or an additional thermal component with a temperature of $\sim$2.0 keV. We investigate the possible origin of this component and suggest that it could arise from either the Galactic ridge X-ray emission, an unidentified pulsar wind nebula or X-ray synchrotron emission from high-energy particles accelerated at the shock. However, deeper, high resolution observations of this object are needed to shed light on the presence and origin of this feature. Based on its morphology, its Galactic latitude, the density of the surrounding environment and its association with a dense molecular cloud, G346.6-0.2 most likely arises from a massive progenitor that underwent core-collapse.

We report a systematic analysis of the spin, orbital, and superorbital modulations of 4U 0114+650, a high-mass X-ray binary consists of one of the slowest spinning neutron stars. Utilizing the dynamic power spectrum, we found that the spin period varied dramatically and anti-correlated with the long-term X-ray flux variation that can be observed using RXTE ASM, Swift BAT, and MAXI. The spin-up rate over the entire dataset is consistent with previously reported values; however, the local spin-up rate is considerably higher. The corresponding local spin-up timescale is comparable to the local spin-up rate of OAO 1657-415, indicating that 4U 0114+650 could also have a transient disk. Moreover, the spin period evolution shows two ~1000-day spin-down/random-walk epochs that appeared together with depressions of the superorbital modulation amplitude. This implies that the superorbital modulation was closely related to the presence of the accretion disk, which is not favored in the spin-down/random-walk epochs because the accretion is dominated by the direct wind accretion.The orbital period is stable during the entire time span; however, the orbital profile changes with time significantly. We found that the depth of the dip near the inferior conjunction of the companion is highly variable, which disfavors the eclipsing scenario. Moreover, the dip was less obvious during the spin-down/random-walk epochs, indicating its correlation with the accretion disk. Further monitoring in both X-ray and optical bands could reveal the establishment of the accretion disk in this system.

Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud (LMC) is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array (ALMA) are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks ("nickel heating"). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained.

In 8 years of operation, the Large Area Telescope (LAT) on the Fermi satellite has impacted our understanding of gamma-ray pulsars dramatically. The LAT now sees over two hundred pulsars: the largest class of GeV sources in the Milky Way. They are diverse -- radio loud versus quiet, young versus millisecond, in evolving binary systems versus isolated, and so on. Relatively few of the GeV pulsars have also been seen in soft gamma rays. After an overview, we present 10 new radio pulsars, six young and four recycled, for which we detect gamma-ray pulsations.

A glitch of a pulsar is known as a sudden increase in the spin frequency and spin-down rate (frequency time derivative), and it can be caused by a sudden rel\ ease of the stress built up in the solid crust of the star or pinned vortices in the superfluid interior. PSR J2021+4026 is the first pulsar that shows a significant change in the gamma-ray flux and pulse profile at the glitch that occurred around 2011 October 16. We report the results of timing and spectral analysis of PSR~J2021+4026 using $\sim$ 8~yr Fermi-LAT data. We find that the pulsar stayed at a high spin-down rate ($\sim 4\%$ higher than the pre-glitch value) and a low gamma-ray state ($\sim 18\%$ lower) for about 3~yr after the glitch. Around 2014 December, the spin-down rate and gamma-ray flux gradually returned to pre-glitch values within a time scale of a few months. The phase-resolved spectra and pulse profiles after the relaxation are also consistent with those before the glitch. The observed long-term evolution of the spin-down rate and the gamma-ray flux indicates that the glitch triggered a mode change in the global magnetosphere. We speculate that the glitch changed the local magnetic field structure around the polar cap and/or the inclination angle of the dipole axis, leading to a change in the electric current circulating in the magnetosphere.

We report the discovery of a synchrotron nebula, G283.1-0.59, associated with PSR J1015-5719. Radio observations using the Molonglo Observatory Synthesis Telescope and the Australia Telescope Compact Array at 36, 16, 6, and 3 cm reveal a complex morphology. The pulsar is embedded in the "head" of the nebula with fan-shaped diffuse emission. This is connected to a circular bubble of 20" radius and a collimated tail extending over 1'. Polarization measurements show a highly ordered magnetic field in the nebula. It wraps around the edge of the head and shows an azimuthal configuration near the pulsar, then switches direction quasi-periodically near the bubble and in the tail. Together with the flat radio spectrum observed, we suggest that this system is most plausibly a pulsar wind nebula (PWN), with the head as a bow shock that has a low Mach number and the bubble as a shell expanding in a dense environment. The bubble could act as a magnetic bottle trapping the relativistic particles. A comparison with other bow-shock PWNe with higher Mach numbers shows similar structure and B-field geometry, implying that pulsar velocity may not be the most critical factor in determining the properties of these systems. We also derive analytic expressions for the projected standoff distance and shape of an inclined bow shock. It is found that the projected distance is always larger than the true distance in three dimensions. On the other hand, the projected shape is not sensitive to the inclination after rescaling with the projected standoff distance.

The Type IIn supernova (SN) 1978K in the nearby galaxy NGC 1313 has remained bright in multiwavelengths for more than $\sim 25$ yr. The archival data of SN 1978K collected with ROSAT, ASCA, XMM-Newton, and Chandra from 1990 to 2006 show no significant variation of the soft X-ray emission but a hint of flux decrease in hard X-rays. In this work, we perform a detailed analysis using more than 15 yr of XMM-Newton observations. Both the 0.5-2 keV and 2-10 keV light-curves decline as $t^{-1}$ from 2000 to 2015. The transition of light curve profiles can be explained in a way that the reverse shock was radiative at an early phase and then became adiabatic at late times. Such a scenario is also supported by the spectral analysis results. We also found a decrease in the absorption column density, which indicates the presence of a cool shell during the radiative phase.

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae.

We report on a timing and spectral analysis of the young, high magnetic field rotation-powered pulsar (RPP) B1509-58 using Chandra continuous-clocking mode observation. The pulsar's X-ray light curve can be fit by the two Gaussian components and the pulsed fraction shows moderate energy dependence over the Chandra band. The pulsed X-ray spectrum is well described by a power law with a photon index 1.16(4), which is harder than the values measured with RXTE PCA and NuSTAR. This result supports the log-parabolic model for the broadband X-ray spectrum. With the unprecedented angular resolution of Chandra, we clearly identified off-pulse X-ray emission from the pulsar, and its spectrum is best fit by a power law plus blackbody model. The latter component has a temperature of ~0.14 keV with a bolometric luminosity comparable to the luminosities of other young and high magnetic field RPPs, and it lies between the temperature of magnetars and typical RPPs. In addition, we found that the nonthermal X-ray emission of PSR B1509-58 is significantly softer in the off-pulse phase than in the pulsed phase, with the photon index varying between 1.0 and 1.8 and anticorrelated with the flux. This is similar to the behavior of three other young pulsars. We interpreted it as different contributions of pair-creation processes at different altitudes from the neutron star surface according to the outer-gap model.

The gamma-ray binary 1FGL J1018.6-5856 is known to be composed of a massive star with a compact object with an orbital period of 16.54 days. We study mechanisms of high-energy emissions (X-ray, GeV, and TeV $\gamma$-rays) from this binary system. The paper is composed of two parts. In the first part, we re-analyze the GeV emission from 1FGL J1018.6-5856 using 7 yrs of the observational data of \textit{Fermi}/LAT. The data analysis of phase-resolved spectrum suggests that the 0.1-10 GeV emissions may contain two different components: (1) a simple power-law with exponential cutoff component in the LOW state; (2) an additional component extended to $\sim $10 GeV in the HIGH state. In the second part, we develop a model to explain the multi-wavelength emissions of this binary. Assuming that 1FGL J1018.6-5856 includes a pulsar, the GeV spectrum in the LOW state can be well explained by the two-layer outer gap model of pulsar magnetosphere. As the pulsar is orbiting around the star, the pulsar is considered to be wrapped by two kinds of termination shock: a bow shock due to the interaction between pulsar wind and stellar wind, and a back shock due to the effect of orbital motion. We proposed that X-rays are mainly produced by the synchrotron radiation at the apex of the bow shock, and TeV $\gamma$-rays are mainly produced by the inverse-Compton scattering at the back shock. The sharp peak in X-ray and TeV light curves are produced by in the shock tail due to the Doppler boosting effect at inferior-conjunction while another broad sinusoidal modulations are possibly produced by the apex of the bow shock and the back shock due to variations of shock conditions.

We have developed FFT beamforming techniques for the CHIME radio telescope, to search for and localize the astrophysical signals from Fast Radio Bursts (FRBs) over a large instantaneous field-of-view (FOV) while maintaining the full angular resolution of CHIME. We implement a hybrid beamforming pipeline in a GPU correlator, synthesizing 256 FFT-formed beams in the North-South direction by four formed beams along East-West via exact phasing, tiling a sky area of ~250 square degrees. A zero-padding approximation is employed to improve chromatic beam alignment across the wide bandwidth of 400 to 800 MHz. We up-channelize the data in order to achieve fine spectral resolution of $\Delta\nu$=24 kHz and time cadence of 0.983 ms, desirable for detecting transient and dispersed signals such as those from FRBs.

PSR J2032+4127 is a radio-loud gamma-ray-emitting pulsar; it is orbiting around a high-mass Be type star with a very long orbital period of 25-50years, and is approaching periastron, which will occur in late 2017/early 2018. This system comprises with a young pulsar and a Be type star, which is similar to the so-called gamma-ray binary PSR~B1259-63/LS2883. It is expected therefore that PSR J2032+4127 shows an enhancement of high-energy emission caused by the interaction between the pulsar wind and Be wind/disk around periastron. Ho et al. recently reported a rapid increase in the X-ray flux from this system. In this paper, we also confirm a rapid increase in the X-ray flux along the orbit, while the GeV flux shows no significant change. We discuss the high-energy emissions from the shock caused by the pulsar wind and stellar wind interaction and examine the properties of the pulsar wind in this binary system. We argue that the rate of increase of the X-ray flux observed by Swift indicates (1) a variation of the momentum ratio of the two-wind interaction region along the orbit, or (2) an evolution of the magnetization parameter of the pulsar wind with the radial distance from the pulsar. We also discuss the pulsar wind/Be disk interaction at the periastron passage, and propose the possibility of formation of an accretion disk around the pulsar. We model high-energy emissions through the inverse-Compton scattering process of the cold-relativistic pulsar wind off soft photons from the accretion disk.

NGC 7793 P13 is an ultraluminous X-ray source harboring an accreting pulsar. We report on the detection of a ~65 d period X-ray modulation with Swift observations in this system. The modulation period found in the X-ray band is P=65.05+/-0.10 d and the profile is asymmetric with a fast rise and a slower decay. On the other hand, the u-band light curve collected by Swift UVOT confirmed an optical modulation with a period of P=64.24+/-0.13 d. We explored the phase evolution of the X-ray and optical periodicities and propose two solutions. A superorbital modulation with a period of ~2,700-4,700 d probably caused by the precession of a warped accretion disk is necessary to interpret the phase drift of the optical data. We further discuss the implication if this ~65d periodicity is caused by the superorbital modulation. Estimated from the relationship between the spin-orbital and orbital-superorbital periods of known disk-fed high-mass X-ray binaries, the orbital period of P13 is roughly estimated as 3-7 d. In this case, an unknown mechanism with a much longer time scale is needed to interpret the phase drift. Further studies on the stability of these two periodicities with a long-term monitoring could help us to probe their physical origins.

By comparing the properties of non-recycled radio-loud $\gamma-$ray pulsars and radio-quiet $\gamma-$ray pulsars, we have searched for the differences between these two populations. We found that the $\gamma-$ray spectral curvature of radio-quiet pulsars can be larger than that of radio-loud pulsars. Based on the full sample of non-recycled $\gamma-$ray pulsars, their distributions of the magnetic field strength at the light cylinder are also found to be different. We notice that this might be resulted from the observational bias. In re-examining the previously reported difference of $\gamma-$ray-to-X-ray flux ratios, we found the significance can be hampered by their statistical uncertainties. In the context of outer gap model, we discuss the expected properties of these two populations and compare with the possible differences identified in our analysis.

We report on six new Chandra observations of the Geminga pulsar wind nebula (PWN). The PWN consists of three distinct elongated structures - two $\approx 0.2 d_{250}$ pc long lateral tails and a segmented axial tail of $\approx 0.05 d_{250}$ pc length, where $d_{250}=d/(250 {\rm pc})$. The photon indices of the power law spectra of the lateral tails, $\Gamma \approx 1$, are significantly harder than those of the pulsar ($\Gamma \approx 1.5$) and the axial tail ($\Gamma \approx 1.6$). There is no significant diffuse X-ray emission between the lateral tails -- the ratio of the X-ray surface brightness between the south tail and this sky area is at least 12. The lateral tails apparently connect directly to the pulsar and show indication of moving footpoints. The axial tail comprises time-variable emission blobs. However, there is no evidence for constant or decelerated outward motion of these blobs. Different physical models are consistent with the observed morphology and spectra of the Geminga PWN. In one scenario, the lateral tails could represent an azimuthally asymmetric shell whose hard emission is caused by the Fermi acceleration mechanism of colliding winds. In another scenario, the lateral tails could be luminous, bent polar outflows, while the blobs in the axial tail could represent a crushed torus. In a resemblance to planetary magnetotails, the blobs of the axial tail might also represent short-lived plasmoids which are formed by magnetic field reconnection in the relativistic plasma of the pulsar wind tail.

The radio and gamma-ray pulsar PSR J2032+4127 was recently found to be in a decades-long orbit with the Be star MT91 213, with the pulsar moving rapidly towards periastron. This binary shares many similar characteristics with the previously unique binary system PSR B1259-63/LS 2883. Here, we describe radio, X-ray, and optical monitoring of PSR J2032+4127/MT91 213. Our extended orbital phase coverage in radio, supplemented with Fermi gamma-ray data, allows us to update and refine the orbital period to 45-50 yr and time of periastron passage to 2017 November. We analyse archival and recent Chandra and Swift observations and show that PSR J2032+4127/MT91 213 is now brighter in X-rays by a factor of ~70 since 2002 and ~20 since 2010. While the pulsar is still far from periastron, this increase in X-rays is possibly due to collisions between pulsar and Be star winds. Optical observations of the Halpha emission line of the Be star suggest that the size of its circumstellar disc may be varying by ~2 over timescales as short as 1-2 months. Multiwavelength monitoring of PSR J2032+4127/MT91 213 will continue through periastron passage, and the system should present an interesting test case and comparison to PSR B1259-63/LS 2883.

We report on Chandra X-ray Observatory (CXO) observations of the pulsar wind nebula (PWN) associated with PSR B0355+54 (eight observations with a 395 ks total exposure, performed over an 8 month period). We investigated the spatial and spectral properties of the emission coincident with the pulsar, compact nebula (CN), and extended tail. We find that the CN morphology can be interpreted in a way that suggests a small angle between the pulsar spin axis and our line-of-sight, as inferred from the radio data. On larger scales, emission from the 7' (2 pc) tail is clearly seen. We also found hints of two faint extensions nearly orthogonal to the direction of the pulsar's proper motion. The spectrum extracted at the pulsar position can be described with an absorbed power-law + blackbody model. The nonthermal component can be attributed to magnetospheric emission, while the thermal component can be attributed to emission from either a hot spot (e.g., a polar cap) or the entire neutron star surface. Surprisingly, the spectrum of the tail shows only a slight hint of cooling with increasing distance from the pulsar. This implies either a low magnetic field with fast flow speed, or particle re-acceleration within the tail. We estimate physical properties of the PWN and compare the morphologies of the CN and the extended tail with those of other bow shock PWNe observed with long CXO exposures.

PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z \approx -0.9\,$ dex) main sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.

PSR J1509-5850 is a middle-aged pulsar with the period P ~ 89 ms, spin-down power Edot = 5.1 x 10^35 erg/s, at a distance of about 3.8 kpc. We report on deep Chandra X-ray Observatory observations of this pulsar and its pulsar wind nebula (PWN). In addition to the previously detected tail extending up to 7' southwest from the pulsar (the southern outflow), the deep images reveal a similarly long, faint diffuse emission stretched toward the north (the northern outflow) and the fine structure of the compact nebula (CN) in the pulsar vicinity. The CN is resolved into two lateral tails and one axial tail pointing southwest (a morphology remarkably similar to that of the Geminga PWN), which supports the assumption that the pulsar moves towards the northeast. The luminosities of the southern and northern outflows are about 1 x 10^33 and 4 x 10^32 erg/s, respectively. The spectra extracted from four regions of the southern outflow do not show any softening with increasing distance from the pulsar. The lack of synchrotron cooling suggests a high flow speed or in-situ acceleration of particles. The spectra extracted from two regions of the northern outflow show a hint of softening with distance from the pulsar, which may indicate slower particle propagation. We speculate that the northern outflow is associated with particle leakage from the bow shock apex into the ISM, while the southern outflow represents the tail of the shocked pulsar wind behind the moving pulsar. We estimate the physical parameters of the observed outflows and compare the J1509-5850 PWN with PWNe of other supersonically moving pulsars.

Pulsars are rapidly spinning and highly magnetized neutron stars, with highly stable rotational period and gradual spin-down over a long timescale due to the loss of radiation. Glitches refer to the events that suddenly increase the rotational speed of a pulsar. The exact causes of glitches and the resulting processes are not fully understood. It is generally believed that couplings between the normal matter and the superfluid components, and the starquakes, are the common causes of glitches. In this study, one famous glitching pulsar, PSR~J2021+4026, is investigated. PSR~J2021+4026 is the first variable gamma-ray pulsar observed by Fermi. From the gamma-ray observations, it is found that the pulsar experienced a significant flux drop, an increase in the spin-down rate, a change in the pulse profile and a shift in the spectral cut-off to a lower energy, simultaneously around 2011 October 16. To explain these effects on the high-energy emissions by the glitch of PSR~J2021+4026, we hypothesized the glitch to be caused by the rearrangement of the surface magnetic field due to the crustal plate tectonic activities on the pulsar which is triggered by a starquake. In this glitch event, the inclination angle of the magnetic dipole axis is slightly shifted. This proposition is then tested by numerical modeling using a three-dimensional two-layer outer gap model. The simulation results indicate that a modification of the inclination angle can affect the pulse profile and the spectral properties, which can explain the observation changes after the glitch.

This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.

Pulsar wind nebulae (PWNe) are suggested to be acceleration sites of cosmic rays in the Galaxy. While the magnetic field plays an important role in the acceleration process, previous observations of magnetic field configurations of PWNe are rare, particularly for evolved systems. We present a radio polarization study of the "Snail" PWN inside the supernova remnant G327.1-1.1 using the Australia Telescope Compact Array. This PWN is believed to have been recently crushed by the supernova (SN) reverse shock. The radio morphology is composed of a main circular body with a finger-like protrusion. We detected a strong linear polarization signal from the emission, which reflects a highly ordered magnetic field in the PWN and is in contrast to the turbulent environment with a tangled magnetic field generally expected from hydrodynamical simulations. This could suggest that the characteristic turbulence scale is larger than the radio beam size. We built a toy model to explore this possibility, and found that a simulated PWN with a turbulence scale of about one-eighth to one-sixth of the nebula radius and a pulsar wind filling factor of 50--75% provides the best match to observations. This implies substantial mixing between the SN ejecta and pulsar wind material in this system.

We report the results from a detailed $\gamma-$ray investigation in the field of two "dark accelerators", HESS J1745-303 and HESS J1741-302, with $6.9$ years of data obtained by the Fermi Large Area Telescope. For HESS J1745-303, we found that its MeV-GeV emission is mainly originated from the "Region A" of the TeV feature. Its $\gamma-$ray spectrum can be modeled with a single power-law with a photon index of $\Gamma\sim2.5$ from few hundreds MeV to TeV. Moreover, an elongated feature, which extends from "Region A" toward northwest for $\sim1.3^{\circ}$, is discovered for the first time. The orientation of this feature is similar to that of a large scale atomic/molecular gas distribution. For HESS J1741-302, our analysis does not yield any MeV-GeV counterpart for this unidentified TeV source. On the other hand, we have detected a new point source, Fermi J1740.1-3013, serendipitously. Its spectrum is apparently curved which resembles that of a $\gamma-$ray pulsar. This makes it possibly associated with PSR B1737-20 or PSR J1739-3023.

The detection of five new fast radio bursts (FRBs) found in the High Time Resolution Universe high latitude survey is presented. The rate implied is 6$^{+4}_{-3}\times~10^3$ (95%) FRBs sky$^{-1}$ day$^{-1}$ above a fluence of between 0.13 and 5.9 Jy ms for FRBs between 0.128 and 262 ms in duration. One of these FRBs has a clear two-component profile, each component is similar to the known population of single component FRBs and are separated by 2.4(4) ms. All the FRB components appear to be unresolved following deconvolution with a scattering tail and accounting for intra-channel smearing. The two-component FRB also has the highest dispersion measure (1629 pc cm$^{-3}$) of any FRB to-date. Many of the proposed models to explain FRBs use a single high energy event involving compact objects (such as neutron star mergers) and therefore cannot easily explain a two-component FRB. Models that are based on extreme versions of flaring, pulsing or orbital events however could produce multiple component profiles. The compatibility of these models and the FRB rate implied by these detections is discussed.

We explore a non-stationary outer gap scenario for gamma-ray emission process in pulsar magnetosphere. Electrons/positrons that migrate along the magnetic field line and enter the outer gap from the outer/inner boundaries activate the pair-creation cascade and high-energy emission process. In our model, the rate of the particle injection at the gap boundaries is key physical quantity to control the gap structure and properties of the gamma-ray spectrum. Our model assumes that the injection rate is time variable and the observed gamma-ray spectrum are superposition of the emissions from different gap structures with different injection rates at the gap boundaries. The calculated spectrum superposed by assuming power law distribution of the particle injection rate can reproduce sub-exponential cut-off feature in the gamma-ray spectrum observed by Fermi-LAT. We fit the phase-averaged spectra for 43 young/middle-age pulsars and 14 millisecond pulsars with the model. Our results imply that (1) a larger particle injection at the gap boundaries is more frequent for the pulsar with a larger spin down power and (2) outer gap with an injection rate much smaller than the Goldreich-Julian value produces observe $>10$GeV emissions. Fermi-LAT gamma-ray pulsars show that (i) the observed gamma-ray spectrum below cut-off energy tends to be softer for the pulsar with a higher spin down rate and (ii) the second peak is more prominent in higher energy bands. Based on the results of the fitting, we describe possible theoretical interpretations for these observational properties. We also briefly discuss Crab-like millisecond pulsars that show phase-aligned radio and gamma-ray pulses.

The DEAP-3600 experiment is located 2 km underground at SNOLAB, in Sudbury, Ontario. It is a single-phase detector that searches for dark matter particle interactions within a 1000-kg fiducial mass target of liquid argon. A first generation prototype detector (DEAP-1) with a 7-kg liquid argon target mass demonstrated a high level of pulse-shape discrimination (PSD) for reducing $\beta$/$\gamma$ backgrounds and helped to develop low radioactivity techniques to mitigate surface-related $\alpha$ backgrounds. Construction of the DEAP-3600 detector is nearly complete and commissioning is starting in 2014. The target sensitivity to spin-independent scattering of Weakly Interacting Massive Particles (WIMPs) on nucleons of 10$^{-46}$ cm$^2$ will allow one order of magnitude improvement in sensitivity over current searches at 100 GeV WIMP mass. This paper presents an overview and status of the DEAP-3600 project and discusses plans for a future multi-tonne experiment, DEAP-50T.

Several theories exist to explain the source of the bright, millisecond duration pulses known as fast radio bursts (FRBs). If the progenitors of FRBs are non-cataclysmic, such as giant pulses from pulsars, pulsar-planet binaries, or magnetar flares, FRB emission may be seen to repeat. We have undertaken a survey of the fields of eight known FRBs from the High Time Resolution Universe survey to search for repeating pulses. Although no repeat pulses were detected the survey yielded the detection of a new FRB, described in Petroff et al. (2015a). From our observations we rule out periodic repeating sources with periods P $\leq$ 8.6 hours and rule out sources with periods 8.6 < P < 21 hours at the 90% confidence level. At P $\geq$ 21 hours our limits fall off as ~1/P. Dedicated and persistent observations of FRB source fields are needed to rule out repetition on longer timescales, a task well-suited to next generation wide-field transient detectors.

Central compact objects (CCOs) are a handful of young neutron stars found at the center of supernova remnants (SNRs). They show high thermal X-ray luminosities but no radio emission. Spin-down rate measurements of the three CCOs with X-ray pulsations indicate surface dipole fields much weaker than those of typical young pulsars. To investigate if CCOs and known radio pulsars are objects at different evolutionary stages, we carried out a census of all weak-field (<1e11 G) isolated radio pulsars in the Galactic plane to search for CCO-like X-ray emission. None of the 12 candidates is detected at X-ray energies, with luminosity limits of 1e32-1e34 erg/s. We consider a scenario in which the weak surface fields of CCOs are due to rapid accretion of supernova materials and show that as the buried field diffuses back to the surface, a CCO descendant is expected to leave the P-Pdot parameter space of our candidates at a young age of a few times 10kyr. Hence, the candidates are likely to be just old ordinary pulsars in this case. We suggest that further searches for orphaned CCO, which are aged CCOs with parent SNRs dissipated, should include pulsars with stronger magnetic fields.

We present initial results from the low-latitude Galactic plane region of the High Time Resolution Universe pulsar survey conducted at the Parkes 64-m radio telescope. We discuss the computational challenges arising from the processing of the terabyte-sized survey data. Two new radio interference mitigation techniques are introduced, as well as a partially-coherent segmented acceleration search algorithm which aims to increase our chances of discovering highly-relativistic short-orbit binary systems, covering a parameter space including potential pulsar-black hole binaries. We show that under a constant acceleration approximation, a ratio of data length over orbital period of ~0.1 results in the highest effectiveness for this search algorithm. From the 50 per cent of data processed thus far, we have re-detected 435 previously known pulsars and discovered a further 60 pulsars, two of which are fast-spinning pulsars with periods less than 30ms. PSR J1101-6424 is a millisecond pulsar whose heavy white dwarf (WD) companion and short spin period of 5.1ms indicate a rare example of full-recycling via Case A Roche lobe overflow. PSR J1757-27 appears to be an isolated recycled pulsar with a relatively long spin period of 17ms. In addition, PSR J1244-6359 is a mildly-recycled binary system with a heavy WD companion, PSR J1755-25 has a significant orbital eccentricity of 0.09, and PSR J1759-24 is likely to be a long-orbit eclipsing binary with orbital period of the order of tens of years. Comparison of our newly-discovered pulsar sample to the known population suggests that they belong to an older population. Furthermore, we demonstrate that our current pulsar detection yield is as expected from population synthesis.

We obtained six observations of PSR J1741-2054 using the $Chandra$ ACIS-S detector totaling $\sim$300 ks. By registering this new epoch of observations to an archival observation taken 3.2 years earlier using X-ray point sources in the field of view, we have measured the pulsar proper motion at $\mu =109 \pm 10 {\rm mas\ yr}^{-1}$ in a direction consistent with the symmetry axis of the observed H$\alpha$ nebula. We investigated the inferred past trajectory of the pulsar but find no compelling association with OB associations in which the progenitor may have originated. We confirm previous measurements of the pulsar spectrum as an absorbed power law with photon index $\Gamma$=2.68$\pm$0.04, plus a blackbody with an emission radius of (4.5$^{+3.2}_{-2.5})d_{0.38}$ km, for a DM-estimated distance of $0.38d_{0.38}$ kpc and a temperature of $61.7\pm3.0$ eV. Emission from the compact nebula is well described by an absorbed power law model with a photon index of $\Gamma$ = 1.67$\pm$0.06, while the diffuse emission seen as a trail extending northeast of the pulsar shows no evidence of synchrotron cooling. We also applied image deconvolution techniques to search for small-scale structures in the immediate vicinity of the pulsar, but found no conclusive evidence for such structures.

Neutron stars lose the bulk of their rotational energy in the form of a pulsar wind: an ultra-relativistic outflow of predominantly electrons and positrons. This pulsar wind significantly impacts the environment and possible binary companion of the neutron star, and studying the resultant pulsar wind nebulae is critical for understanding the formation of neutron stars and millisecond pulsars, the physics of the neutron star magnetosphere, the acceleration of leptons up to PeV energies, and how these particles impact the interstellar medium. With the SKA1 and the SKA2, it could be possible to study literally hundreds of PWNe in detail, critical for understanding the many open questions in the topics listed above.

We report on two years of flux and spin evolution monitoring of 1E 1048.1$-$5937, a 6.5-s X-ray pulsar identified as a magnetar. Using {\it Swift} XRT data, we observed an X-ray outburst consisting of an increase in the persistent 1--10 keV flux by a factor of 6.3$\pm$0.2, beginning on 2011 December 31 (MJD 55926). Following a delay of $\sim100$ days, the magnetar entered a period of large torque variability, with $\dot{\nu}$ reaching a factor of $4.55\pm0.05$ times the nominal value, before decaying in an oscillatory manner over a time scale of months. We show by comparing to previous outbursts from the source that this pattern of behavior may repeat itself with a quasi-period of $\sim1800$ days. We compare this phenomenology to periodic torque variations in radio pulsars, finding some similarities which suggest a magnetospheric origin for the behavior of 1E 1048.1$-$5937.

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$\sigma$) circularly polarized on the leading edge with a 1-$\sigma$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

We present a comprehensive spectral and morphological analysis of the remnant of Supernova (SN) 1987A with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA). The non-thermal and thermal components of the radio emission are investigated in images from 94 to 672 GHz ($\lambda$ 3.2 mm to 450 $\mu$m), with the assistance of a high-resolution 44 GHz synchrotron template from the ATCA, and a dust template from ALMA observations at 672 GHz. An analysis of the emission distribution over the equatorial ring in images from 44 to 345 GHz highlights a gradual decrease of the east-to-west asymmetry ratio with frequency. We attribute this to the shorter synchrotron lifetime at high frequencies. Across the transition from radio to far infrared, both the synchrotron/dust-subtracted images and the spectral energy distribution (SED) suggest additional emission beside the main synchrotron component ($S_{\nu}\propto\nu^{-0.73}$) and the thermal component originating from dust grains at $T\sim22$ K. This excess could be due to free-free flux or emission from grains of colder dust. However, a second flat-spectrum synchrotron component appears to better fit the SED, implying that the emission could be attributed to a pulsar wind nebula (PWN). The residual emission is mainly localised west of the SN site, as the spectral analysis yields $-0.4\lesssim\alpha\lesssim-0.1$ across the western regions, with $\alpha\sim0$ around the central region. If there is a PWN in the remnant interior, these data suggest that the pulsar may be offset westward from the SN position.

We present the discovery of a further five recycled pulsar systems in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods ranging from 2 ms to 66 ms, and four are in binary systems with orbital periods between 10.8 hours and 9.0 days. Three of these binary systems are particularly interesting; PSR J1227-6208 has a pulse period of 34.5 ms and the highest mass function of all pulsars with near-circular orbits. The circular orbit suggests that the companion is not another neutron star, so future timing experiments may reveal one of the heaviest white dwarfs ever found ($>$ 1.3 M$_\odot$). Timing observations of PSR J1431$-$4715 indicate that it is eclipsed by its companion which has a mass indicating it belongs to the redback class of eclipsing millisecond pulsars. PSR J1653-2054 has a companion with a minimum mass of only $0.08$ M$_\odot$, placing it among the class of pulsars with low-mass companions. Unlike the majority of such systems, however, no evidence of eclipses is seen at 1.4 GHz.

The First \fermi-LAT Catalog of Sources Above 10 GeV reported evidence of pulsed emission above 25 GeV from 12 pulsars, including the Vela pulsar, which showed evidence of pulsation at $>37$ GeV energy bands. Using 62 months of \fermi-LAT data, we analyzed the gamma-ray emission from the Vela pulsar and searched for pulsed emission above 50 GeV. Having confirmed the significance of the pulsation in 30-50 GeV with the H-test (p-value $\sim10^{-77}$), we extracted its pulse profile using the Bayesian block algorithm and compared it with the distribution of the 5 observed photons above 50 GeV using the likelihood ratio test. Pulsation was significantly detected for photons above 50 GeV with p-value $=3\times10^{-5}$ ($4.2\sigma$). The detection of pulsation is significant above $4\sigma$ at $>79$ GeV and above $3\sigma$ at $>90$ GeV energy bands, making this the highest energy pulsation significantly detected by the LAT. We explore non-stationary outer gap scenario of the very high-energy emissions from the Vela pulsar.

We report the discovery of PSR J1101-6101, a 62.8 ms pulsar in IGR J11014-6103, a hard X-ray source with a jet and a cometary tail that strongly suggests it is moving away from the center of the supernova remnant (SNR) MSH 11-61A at v>1000 km/s. Two XMM-Newton observations were obtained with the EPIC pn in small window mode, resulting in the measurement of its spin-down luminosity E_dot = 1.36e36 erg/s, characteristic age Tau_c = 116 kyr, and surface magnetic field strength B_s = 7.4e11 G. In comparison to Tau_c, the 10-30 kyr age estimated for MSH 11-61A suggests that the pulsar was born in the SNR with initial period in the range 54 < P_0 < 60 ms. PSR J1101-6101 is the least energetic of the 15 rotation-powered pulsars detected by INTEGRAL, and has a high efficiency of hard X-ray radiation and jet power. We examine the shape of the cometary nebula in a Chandra image, which is roughly consistent with a bow shock at the velocity inferred from the SNR age and the pulsar's E-dot. However, its structure differs in detail from the classic bow shock, and we explore possible reasons for this.

The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There are still a large number of outstanding questions, such the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove & McKee progenitor with an envelope mass of $10 M_{\sun}$ and an energy of $1.5 \times 10^{44} J$. A termination shock in the progenitor's stellar wind at a distance of $0\farcs43-0\farcs51$ provides a good fit to the turn on of radio emission around day 1200. For the H\textsc{ii} region, a minimum distance of $0\farcs63\pm0\farcs01$ and maximum particle number density of $(7.11\pm1.78) \times 10^7$ m$^{-3}$ produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring.

The Fermi source 2FGL~J0523.3$-$2530 has recently been identified as a candidate millisecond pulsar binary with an orbital period of 16.5 hrs. We have carried out detailed studies of the source's emission properties by analyzing data taken with the Fermi Large Area Telescope in the 0.2--300 GeV energy range. Long-term, yearly variability from the source has been found, with a factor of 4 flux variations in 1--300 GeV. From spectral analysis, we find an extra spectral component at 2--3 GeV that causes the source brightening. While no orbital modulations have been found from the Fermi data over the whole period of 2008--2014, orbital modulation in the source's $>$2 GeV emission is detected during the last 1.5 yrs of the Fermi observation. Our results support the millisecond pulsar binary nature of 2FGL~J0523.3$-$2530. Multi-wavelength observations of the source are warranted in order to find any correlated flux variations and thus help determine the origin of the long-term variability, which currently is not understood.

There is a growing new class of young spin-down powered pulsars called GeV-quiet soft gamma-ray pulsar; (1) spectral turnover appears around~10MeV, (2) the X-ray spectra of below 20 keV can be described by power law with photon index around 1.2 and (3) the light curve in X-ray/soft gamma-ray bands shows single broad pulse. Their emission properties are distinct from the normal gamma-ray pulsars, for which the spectral peak in $\nu F_{\nu}$ appears in GeV energy bands and the X-ray/gamma-ray light curves show sharp and double (or more) peaks. In this paper, we discuss that X-ray/soft gamma-ray emissions of the GeV-quiet soft gamma-ray pulsars are caused bythe synchrotron radiation of the electron/positron pairs, which are created by the magnetic pair-creation process near the stellar surface. In our model, the viewing geometry is crucial factor to discriminate between the normal gamma-ray pulsars and soft gamma-ray pulsars. Our model suggests that the difference between the magnetic inclination angle ($\alpha$) and the Earth viewing angle ($\beta$) of the soft gamma-ray pulsars is small, so that the synchrotron emissions from the high magnetic field region around the polar cap region dominates in the observed emissions. Furthermore, the inclination angle of the soft gamma-ray pulsar is relatively small, $\alpha\leq 30$~degree, and our line of sight is out of the gamma-ray beam emitted via the curvature radiation process in the outer gap. We also analysis the six year $Fermi$ data for four soft gamma-ray pulsars to determine the upper limit of the GeV flux.

Central compact objects (CCOs) constitute a population of radio-quiet, slowly-spinning ($\ge$100 ms) young neutron stars with anomalously high thermal X-ray luminosities. Their spin-down properties imply weak dipole magnetic fields ($\sim$$10^{10-11}$ G) and characteristic ages much greater than the ages of their host supernova remnants (SNRs). However, CCOs may possess strong "hidden" internal magnetic fields that may re-emerge on timescales $\gtrsim$10 kyr, with the neutron star possibly activating as a radio pulsar in the process. This suggests that the immediate descendants of CCOs may be masquerading as slowly spinning "old" radio pulsars. We present an X-ray survey of all ordinary radio pulsars within 6 kpc that are positionally coincident with Galactic SNRs in order to test the possible connection between the supposedly old, but possibly very young pulsars, and the SNRs. None of the targets exhibits anomalously high thermal X-ray luminosity, suggesting that they are genuine old ordinary pulsars unrelated to the superposed SNRs. This implies that CCOs are either latent radio pulsars that activate long after their SNRs dissipate or they remain permanently radio-quiet. The true descendants of CCOs remain at large.

We report on our searches for debris disks around seven relatively nearby radio pulsars, which are isolated sources and were carefully selected as the targets on the basis of our deep $K_s$-band imaging survey. The $K_s$ images obtained with the 6.5\,m Baade Magellan Telescope at Las Campanas Observatory are analyzed together with the \textit{Spitzer}/IRAC images at 4.5 and 8.0~$\mu$m and the \textit{WISE} images at 3.4, 4.6, 12 and 22~$\mu$m. No infrared (IR) counterparts of these pulsars are found, with flux upper limits of $\sim \mu$Jy at near-infrared ($\lambda<10\ \mu$m) and $\sim$10--1000\,$\mu$Jy at mid-infrared wavelengths ($\lambda>10\ \mu$m). The results of this search are discussed in terms of the efficiency of converting the pulsar spin-down energy to thermal energy and X-ray heating of debris disks, with comparison made to the two magnetars 4U~0142+61 and 1E~2259+586 which are suggested to harbor a debris disk.

Fast radio bursts (FRBs) are an emerging class of bright, highly dispersed radio pulses. Recent work by Thornton et al. (2013) has revealed a population of FRBs in the High Time Resolution Universe (HTRU) survey at high Galactic latitudes. A variety of progenitors have been proposed including cataclysmic events at cosmological distances, Galactic flare stars, and terrestrial radio frequency interference. Here we report on a search for FRBs at intermediate Galactic latitudes ($-15^{\circ}$ $< b <$ 15$^{\circ}$) in data taken as part of the HTRU survey. No FRBs were discovered in this region. Several effects such as dispersion, scattering, sky temperature and scintillation decrease the sensitivity by more than 3$\sigma$ in $\sim$20\% of survey pointings. Including all of these effects, we exclude the hypothesis that FRBs are uniformly distributed on the sky with 99\% confidence. This low probability implies that additional factors -- not accounted for by standard Galactic models -- must be included to ease the discrepancy between the detection rates at high and low Galactic latitudes. A revised rate estimate or another strong and heretofore unknown selection effect in Galactic latitude would provide closer agreement between the surveys' detection rates. The dearth of detections at low Galactic latitude disfavors a Galactic origin for these bursts.

We report on X-ray and gamma-ray observations of the millisecond pulsar B1937+21 taken with the Chandra X-ray Observatory, XMM-Newton, and the Fermi Large Area Telescope. The pulsar X-ray emission shows a purely non-thermal spectrum with a hard photon index of 0.9+/-0.1, and is nearly 100% pulsed. We found no evidence of varying pulse profile with energy as previously claimed. We also analyzed 5.5 yr of Fermi survey data and obtained much improved constraints on the pulsar's timing and spectral properties in gamma-rays. The pulsed spectrum is adequately fitted by a simple power-law with a photon index of 2.38+/-0.07. Both the gamma-ray and X-ray pulse profiles show similar two-peak structure and generally align with the radio peaks. We found that the aligned profiles and the hard spectrum in X-rays seem to be common properties among millisecond pulsars with high magnetic fields at the light cylinder. We discuss a possible physical scenario that could give rise to these features.

We present a long-term X-ray flux and spectral analysis for 1RXS J170849.0-400910 using Swift/XRT spanning over 8 years from 2005-2013. We also analyze two observations from Chandra and XMM in the period from 2003-2004. In this 10-yr period, 1RXS J170849.0-400910 displayed several rotational glitches. Previous studies have claimed variations in the X-ray emission associated with some of the glitches. From our analysis we find no evidence for significant X-ray flux variations and evidence for only low-level spectral variations. We also present an updated timing solution for 1RXS J170849.0-400910, from RXTE and Swift observations, which includes a previously unreported glitch at MJD 56019. We discuss the frequency and implications of radiatively quiet glitches in magnetars.

The extreme conditions found in and around pulsars make them fantastic natural laboratories, providing insights to a rich variety of fundamental physics and astronomy. To discover more pulsars we have begun the High Time Resolution Universe (HTRU) survey: a blind survey of the northern sky with the 100-m Effelsberg radio telescope in Germany and a twin survey of the southern sky with the 64-m Parkes radio telescope in Australia. The HTRU is an international collaboration with expertise shared among the MPIfR in Germany, ATNF/CASS and Swinburne University of Technology in Australia, University of Manchester in the UK and INAF in Italy. The HTRU survey uses multi-beam receivers and backends constructed with recent advancements in technology, providing unprecedentedly high time and frequency resolution, allowing us to probe deeper into the Galaxy than ever before. While a general overview of HTRU has been given by Keith at this conference, here we focus on three further aspects of HTRU discoveries and highlights. These include the `Diamond-planet pulsar' binary PSR J1719-1438 and a second similar system recently discovered. In addition, we provide specifications of the HTRU-North survey and an update of its status. In the last section we give an overview of the search for highly-accelerated binaries in the Galactic plane region. We discuss the computational challenges arising from the processing of the petabyte-sized HTRU survey data. We present an innovative segmented search technique which aims to increase our chances of discovering highly accelerated relativistic binary systems, potentially including pulsar-black-hole binaries.

We report on the discovery of four millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) pulsar survey being conducted at the Parkes 64-m radio telescope. All four MSPs are in binary systems and are likely to have white dwarf companions. In addition, we present updated timing solutions for 12 previously published HTRU MSPs, revealing new observational parameters such as five proper motion measurements and significant temporal dispersion measure variations in PSR J1017-7156. We discuss the case of PSR J1801-3210, which shows no significant period derivative after four years of timing data. Our best-fit solution shows a period derivative of the order of $10^{-23}$, an extremely small number compared to that of a typical MSP. However, it is likely that the pulsar lies beyond the Galactic Centre, and an unremarkable intrinsic period derivative is reduced to close to zero by the Galactic potential acceleration. Furthermore, we highlight the potential to employ PSR J1801-3210 in the strong equivalence principle test due to its wide and circular orbit. In a broader comparison with the known MSP population, we suggest a correlation between higher mass functions and the presence of eclipses in `very low-mass binary pulsars', implying that eclipses are observed in systems with high orbital inclinations. We also suggest that the distribution of the total mass of binary systems is inversely-related to the Galactic height distribution. Finally, we report on the first detection of PSRs J1543-5149 and J1811-2404 as gamma-ray pulsars.

Supernova (SN) explosions are crucial engines driving the evolution of galaxies by shock heating gas, increasing the metallicity, creating dust, and accelerating energetic particles. In 2012 we used the Atacama Large Millimeter/Submillimeter Array to observe SN 1987A, one of the best-observed supernovae since the invention of the telescope. We present spatially resolved images at 450um, 870um, 1.4mm, and 2.8mm, an important transition wavelength range. Longer wavelength emission is dominated by synchrotron radiation from shock-accelerated particles, shorter wavelengths by emission from the largest mass of dust measured in a supernova remnant (>0.2Msun). For the first time we show unambiguously that this dust has formed in the inner ejecta (the cold remnants of the exploded star's core). The dust emission is concentrated to the center of the remnant, so the dust has not yet been affected by the shocks. If a significant fraction survives, and if SN 1987A is typical, supernovae are important cosmological dust producers.

We report the discovery of four gamma-ray pulsars, detected in computing-intensive blind searches of data from the Fermi Large Area Telescope (LAT). The pulsars were found using a novel search approach, combining volunteer distributed computing via Einstein@Home and methods originally developed in gravitational-wave astronomy. The pulsars PSRs J0554+3107, J1422-6138, J1522-5735, and J1932+1916 are young and energetic, with characteristic ages between 35 and 56 kyr and spin-down powers in the range $6\times10^{34}$ - $10^{36}$ erg s$^{-1}$. They are located in the Galactic plane and have rotation rates of less than 10 Hz, among which the 2.1 Hz spin frequency of PSR J0554+3107 is the slowest of any known gamma-ray pulsar. For two of the new pulsars, we find supernova remnants coincident on the sky and discuss the plausibility of such associations. Deep radio follow-up observations found no pulsations, suggesting that all four pulsars are radio-quiet as viewed from Earth. These discoveries, the first gamma-ray pulsars found by volunteer computing, motivate continued blind pulsar searches of the many other unidentified LAT gamma-ray sources.

As the remnant of Supernova (SN) 1987A has been getting brighter over time, new observations at high frequencies have allowed imaging of the radio emission at unprecedented detail. We present a new radio image at 44 GHz of the supernova remnant (SNR), derived from observations performed with the Australia Telescope Compact Array (ATCA) in 2011. The diffraction-limited image has a resolution of $349\times225$ mas, which is the highest achieved to date in high-dynamic range images of the SNR. We also present a new image at 18 GHz, also derived from ATCA observations performed in 2011, which is super-resolved to $0./!/!^{\prime\prime}25$. The new 44 and 18 GHz images yield the first high-resolution spectral index map of the remnant. The comparison of the 44 GHz image with contemporaneous X-ray and H$\alpha$ observations allows further investigations of the nature of the remnant asymmetry and sheds more light into the progenitor hypotheses and SN explosion. In light of simple free-free absorption models, we discuss the likelihood of detecting at 44 GHz the possible emission originating from a pulsar wind nebula (PWN) or a compact source in the centre of the remnant.

Supernovae and their remnants are believed to be prodigious sources of Galactic cosmic rays and interstellar dust. Understanding the mechanisms behind their surprisingly high production rate is helped by the study of nearby young supernova remnants. There has been none better in modern times than SN1987A, for which radio observations have been made for over a quarter of a century. We review extensive observations made with the Australia Telescope Compact Array (ATCA) at centimetre wavelengths. Emission at frequencies from 1 to 100 GHz is dominated by synchrotron radiation from an outer shock front which has been growing exponentially in strength from day 3000, and is currently sweeping around the circumstellar ring at about 4000 km s$^{-1}$. Three dimensional models of the propagation of the shock into the circumstellar medium are able to reproduce the main observational features of the remnant, and their evolution. We find that up to 4% of the electrons encountered by the shock are accelerated to relativistic energies. High-frequency ALMA observations will break new ground in the understanding of dust and molecule production.

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.

We present a polarimetric analysis of 49 long-period pulsars discovered as part of the High Time Resolution Universe (HTRU) southern survey. The sources exhibit the typical characteristics of "old" pulsars, with low fractional linear and circular polarisation and narrow, multicomponent profiles. Although the position angle swings are generally complex, for two of the analysed pulsars (J1622-3751 and J1710-2616) we obtained an indication of the geometry via the rotating vector model. We were able to determine a value of the rotation measure (RM) for 34 of the sources which, when combined with their dispersion measures (DM), yields an integrated magnetic field strength along the line of sight. With the data presented here, the total number of values of RM associated to pulsars discovered during the HTRU southern survey sums to 51. The RMs are not consistent with the hypothesis of a counter-clockwise direction of the Galactic magnetic field within an annulus included between 4 and 6 kpc from the Galactic centre. A partial agreement with a counter-clockwise sense of the Galactic magnetic field within the spiral arms is however found in the area of the Carina-Sagittarius arm.

This paper presents the discovery and timing parameters for five millisecond pulsars (MSPs), four in binary systems with probable white dwarf companions and one isolated, found in ongoing processing of the High Time Resolution Universe Pulsar Survey (HTRU). We also present high quality polarimetric data on four of them. These further discoveries confirm the high potential of our survey in finding pulsars with very short spin periods. At least two of these five MSPs are excellent candidates to be included in the Pulsar Timing Array projects. Thanks to the wealth of MSP discoveries in the HTRU survey, we revisit the question of whether the luminosity distributions of isolated and binary MSPs are different. Using the Cordes and Lazio distance model and our new and catalogue flux density measurements, we find that 41 of the 42 most luminous MSPs in the Galactic disk are in binaries and a statistical analysis suggests that the luminosity functions differ with 99.9% significance. We conclude that the formation process that leads to solitary MSPs affects their luminosities, despite their period and period derivatives being similar to those of pulsars in binary systems.

Searches for transient astrophysical sources often reveal unexpected classes of objects that are useful physical laboratories. In a recent survey for pulsars and fast transients we have uncovered four millisecond-duration radio transients all more than 40{\deg} from the Galactic plane. The bursts' properties indicate that they are of celestial rather than terrestrial origin. Host galaxy and intergalactic medium models suggest that they have cosmological redshifts of 0.5 to 1, and distances of up to 3 gigaparsecs. No temporally coincident x- or gamma-ray signature was identified in association with the bursts. Characterization of the source population and identification of host galaxies offers an opportunity to determine the baryonic content of the Universe.

We have used millisecond pulsars (MSPs) from the southern High Time Resolution Universe (HTRU) intermediate latitude survey area to simulate the distribution and total population of MSPs in the Galaxy. Our model makes use of the scale factor method, which estimates the ratio of the total number of MSPs in the Galaxy to the known sample. Using our best fit value for the z-height, z=500 pc, we find an underlying population of MSPs of 8.3(\pm 4.2)*10^4 sources down to a limiting luminosity of L_min=0.1 mJy kpc^2 and a luminosity distribution with a steep slope of d\log N/d\log L = -1.45(\pm 0.14). However, at the low end of the luminosity distribution, the uncertainties introduced by small number statistics are large. By omitting very low luminosity pulsars, we find a Galactic population above L_min=0.2 mJy kpc^2 of only 3.0(\pm 0.7)*10^4 MSPs. We have also simulated pulsars with periods shorter than any known MSP, and estimate the maximum number of sub-MSPs in the Galaxy to be 7.8(\pm 5.0)*10^4 pulsars at L=0.1 mJy kpc^2. In addition, we estimate that the high and low latitude parts of the southern HTRU survey will detect 68 and 42 MSPs respectively, including 78 new discoveries. Pulsar luminosity, and hence flux density, is an important input parameter in the model. Some of the published flux densities for the pulsars in our sample do not agree with the observed flux densities from our data set, and we have instead calculated average luminosities from archival data from the Parkes Telescope. We found many luminosities to be very different than their catalogue values, leading to very different population estimates. Large variations in flux density highlight the importance of including scintillation effects in MSP population studies.

Magnetars are neutron stars showing dramatic X-ray and soft $\gamma$-ray outbursting behaviour that is thought to be powered by intense internal magnetic fields. Like conventional young neutron stars in the form of radio pulsars, magnetars exhibit "glitches" during which angular momentum is believed to be transferred between the solid outer crust and the superfluid component of the inner crust. Hitherto, the several hundred observed glitches in radio pulsars and magnetars have involved a sudden spin-up of the star, due presumably to the interior superfluid rotating faster than the crust. Here we report on X-ray timing observations of the magnetar 1E 2259+586 which we show exhibited a clear "anti-glitch" -- a sudden spin down. We show that this event, like some previous magnetar spin-up glitches, was accompanied by multiple X-ray radiative changes and a significant spin-down rate change. This event, if of origin internal to the star, is unpredicted in models of neutron star spin-down and is suggestive of differential rotation in the neutron star, further supporting the need for a rethinking of glitch theory for all neutron stars.

Forecasting large Solar Energetic Particle (SEP) events associated with shocks driven by fast coronal mass ejections (CME) pose a major difficulty in the field of Space Weather. Besides issues associated with CME initiation, the SEP intensities are difficult to predict, spanning 3 orders of magnitude at any given CME speed. Many lines of indirect evidence point to the pre-existence of suprathermal seed particles for injection into the acceleration process as a key ingredient limiting the SEP intensity of a given event. This paper outlines the observational and theoretical basis for the inference that a suprathermal particle population is present prior to large SEP events, explores various scenarios for generating seed particles and their observational signatures, and explains how such suprathermals could be detected through measuring the wings of the H I Ly-alpha line.

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.

Pulsars are remarkable objects that emit across the entire electromagnetic spectrum, providing a powerful probe of the interstellar medium. In this study, we investigate the relation between dispersion measure (DM) and X-ray absorption column density NH using 68 radio pulsars detected at X-ray energies with the Chandra X-ray Observatory or XMM-Newton. We find a best-fit empirical linear relation of NH (10^20 cm^-2)= 0.30+0.13-0.09 DM (pc cm^-3), which corresponds to an average ionization of 10+4-3%, confirming the ratio of one free electron per ten neutral hydrogen atoms commonly assumed in the literature. We also compare different NH estimates and note that some NH values obtained from X-ray observations are higher than the total Galactic HI column density along the same line of sight, while the optical extinction generally gives the best NH predictions.

In a recent numerical study [Ng et al., Astrophys. J. {\bf 747}, 109, 2012], with a three-dimensional model of coronal heating using reduced magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate versus Lundquist number based on a series of runs in which random photospheric motions are imposed for hundreds to thousands of \al time in order to obtain converged statistical values. The heating rate found in these simulations saturate to a level that is independent of the Lundquist number. This scaling result was also supported by an analysis with the assumption of the Sweet-Parker scaling of the current sheets, as well as how the width, length and number of current sheets scale with Lundquist number. In order to test these assumptions, we have implemented an automated routine to analyze thousands of current sheets in these simulations and return statistical scalings for these quantities. It is found that the Sweet-Parker scaling is justified. However, some discrepancies are also found and require further study.

High-magnetic-field radio pulsars are important transition objects for understanding the connection between magnetars and conventional radio pulsars. We present a detailed study of the young radio pulsar J1119-6127, which has a characteristic age of 1900yr and a spin-down-inferred magnetic field of 4.1e13G, and its associated supernova remnant G292.2-0.5, using deep XMM-Newton and Chandra X-ray Observatory exposures of over 120ks from each telescope. The pulsar emission shows strong modulation below 2.5keV, with a single-peaked profile and a large pulsed fraction of 0.48+/-0.12. Employing a magnetic, partially ionized hydrogen atmosphere model, we find that the observed pulse profile can be produced by a single hot spot of temperature 0.13keV covering about one third of the stellar surface, and we place an upper limit of 0.08keV for an antipodal hot spot with the same area. The nonuniform surface temperature distribution could be the result of anisotropic heat conduction under a strong magnetic field, and a single-peaked profile seems common among high-B radio pulsars. For the associated remnant G292.2-0.5, its large diameter could be attributed to fast expansion in a low-density wind cavity, likely formed by a Wolf-Rayet progenitor, similar to two other high-B radio pulsars.

On 2011 July 14, a transient X-ray source, Swift J1822.3-1606, was detected by Swift BAT via its burst activities. It was subsequently identified as a new magnetar upon the detection of a pulse period of 8.4 s. Using follow-up RXTE, Swift, and Chandra observations, we have determined a spin-down rate of $\dot{P}\sim3\times10^{-13}$, implying a dipole magnetic field of $\sim5\times10^{13}$\,G, second lowest among known magnetars, although our timing solution is contaminated by timing noise. The post-outburst flux evolution is well modelled by surface cooling resulting from heat injection in the outer crust, although we cannot rule out other models. We measure an absorption column density similar to that of the open cluster M17 at 10\arcmin\ away, arguing for a comparable distance of $\sim$1.6\,kpc. If confirmed, this could be the nearest known magnetar.

Swift J1822.3-1606 was discovered on 2011 July 14 by the Swift Burst Alert Telescope following the detection of several bursts. The source was found to have a period of 8.4377 s and was identified as a magnetar. Here we present a phase-connected timing analysis and the evolution of the flux and spectral properties using RXTE, Swift, and Chandra observations. We measure a spin frequency of 0.1185154343(8) s$^{-1}$ and a frequency derivative of $-4.3\pm0.3\times10^{-15}$ at MJD 55761.0, in a timing analysis that include significant non-zero second and third frequency derivatives that we attribute to timing noise. This corresponds to an estimated spin-down inferred dipole magnetic field of $B\sim5\times10^{13}$ G, consistent with previous estimates though still possibly affected by unmodelled noise. We find that the post-outburst 1--10 keV flux evolution can be characterized by a double-exponential decay with decay timescales of $15.5\pm0.5$ and $177\pm14$ days. We also fit the light curve with a crustal cooling model which suggests that the cooling results from heat injection into the outer crust. We find that the hardness-flux correlation observed in magnetar outbursts also characterizes the outburst of Swift J1822.3-1606. We compare the properties of Swift J1822.3-1606 with those of other magnetars and their outbursts.

We present 75 pulsars discovered in the mid-latitude portion of the High Time Resolution Universe survey, 54 of which have full timing solutions. All the pulsars have spin periods greater than 100 ms, and none of those with timing solutions are in binaries. Two display particularly interesting behaviour; PSR J1054-5944 is found to be an intermittent pulsar, and PSR J1809-0119 has glitched twice since its discovery. In the second half of the paper we discuss the development and application of an artificial neural network in the data-processing pipeline for the survey. We discuss the tests that were used to generate scores and find that our neural network was able to reject over 99% of the candidates produced in the data processing, and able to blindly detect 85% of pulsars. We suggest that improvements to the accuracy should be possible if further care is taken when training an artificial neural network; for example ensuring that a representative sample of the pulsar population is used during the training process, or the use of different artificial neural networks for the detection of different types of pulsars.

We report on a new Chandra exposure of PSR J1809-2332, the recently discovered pulsar powering the bright EGRET source 3EG J1809-2328. By registration of field X-ray sources in an archival exposure, we measure a significant proper motion for the pulsar point source over an ~11 year baseline. The shift of 0.30+/-0.06" (at PA= 153.3+/-18.4) supports an association with proposed SNR parent G7.5-1.7. Spectral analysis of diffuse emission in the region also supports the interpretation as a hard wind nebula trail pointing back toward the SNR.

In a recent computational campaign [Ng et al., Astrophys. J. 747, 109, 2012] to investigate a three-dimensional model of coronal heating using reduced magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate versus Lundquist number based on a series of runs in which random photospheric motions are imposed for hundreds to thousands of Alfv\'en time in order to obtain converged statistical values. Using this collection of numerical data, we have performed additional statistical analysis related to the formation of current sheets and heating events, or nanoflares [Parker, Astrophys. J. 330, 474, 1988]. While there have been many observations of the energy distribution of solar flares, there have not been many results based on large-scale three-dimensional direct simulations due to obvious numerical difficulties. We will present energy distributions and other statistics based on our simulations, calculated using a method employed in [Dmitruk & G\'omez, Astrophys. J., 484, L83, 1997]. We will also make comparisons of our results with observations.

The proximity of core-collapse Supernova 1987A (SN1987A) in the Large Magellanic Cloud (LMC) and its rapid evolution make it a unique case study of the development of a young supernova remnant. We aim at resolving the remnant of SN1987A for the first time in the 3-mm band (at 94 GHz). We observed the source at 3-mm wavelength with a 750-m configuration of the Australia Telescope Compact Array (ATCA). We compare the image with a recent 3-cm image and with archival X-ray images. We present a diffraction-limited image with a resolution of 0.7", revealing the ring structure seen at lower frequencies and at other wavebands. The emission peaks in the eastern part of the ring. The 3-mm image bears resemblance to early X-ray images (from 1999-2000). We place an upper limit of 1 mJy (2 \sigma) on any discrete source of emission in the centre (inside of the ring). The integrated flux density at 3 mm has doubled over the six years since the previous observations at 3 mm. At 3 mm - i.e. within the operational domain of the Atacama Large Millimeter/submillimeter Array (ALMA) - SN1987A appears to be dominated by synchrotron radiation from the inner rim of the equatorial ring, characterised by moderately-weak shocks. There is no clear sign of emission of a different nature, but the current limits do not rule out such component altogether.

Parker's model is one of the most discussed mechanisms for coronal heating and has generated much debate. We have recently obtained new scaling results in a two-dimensional (2D) version of this problem suggesting that the heating rate becomes independent of resistivity in a statistical steady state [Ng and Bhattacharjee, Astrophys. J., 675, 899 (2008)]. Our numerical work has now been extended to 3D by means of large-scale numerical simulations. Random photospheric footpoint motion is applied for a time much longer than the correlation time of the motion to obtain converged average coronal heating rates. Simulations are done for different values of the Lundquist number to determine scaling. In the high-Lundquist number limit, the coronal heating rate obtained so far is consistent with a trend that is independent of the Lundquist number, as predicted by previous analysis as well as 2D simulations. In the same limit the average magnetic energy built up by the random footpoint motion tends to have a much weaker dependence on the Lundquist number than that in the 2D simulations, due to the formation of strong current layers and subsequent disruption when the equilibrium becomes unstable. We will present scaling analysis showing that when the dissipation time is comparable or larger than the correlation time of the random footpoint motion, the heating rate tends to become independent of Lundquist number, and that the magnetic energy production is also reduced significantly.

The Frying Pan (G315.9-0.0) is a radio supernova remnant with a peculiar linear feature (G315.78-0.23) extending 10' radially outward from the rim of the shell. We present radio imaging and polarization observations obtained from the Molonglo Observatory Synthesis Telescope and the Australia Telescope Compact Array, confirming G315.78-0.23 as a bow-shock pulsar wind nebula (PWN) powered by the young pulsar J1437-5959. This is one of the longest pulsar tails observed in radio and it has a physical extent over 20 pc. We found a bow-shock stand-off distance of 0.002 pc, smallest among similar systems, suggesting a large pulsar velocity over 1000 km/s and a high Mach number ~200. The magnetic field geometry inferred from radio polarimetry shows a good alignment with the tail orientation, which could be a result of high flow speed. There are also hints that the postshock wind has a low magnetization and is dominated by electrons and positrons in energy. This study shows that PWNe can offer a powerful probe of their local environment, particularly for the case of a bow shock where the parent supernova shell is also detected.

We present an analysis of the extended emission around the anomalous X-ray pulsar 1E 1547.0-5408 using four XMM-Newton observations taken with the source in varying states of outburst as well as in quiescence. We find that the extended emission flux is highly variable and strongly correlated with the flux of the magnetar. Based on this result, as well as on spectral and energetic considerations, we conclude that the extended emission is dominated by a dust-scattering halo and not a pulsar wind nebula (PWN), as has been previously argued. We obtain an upper limit on the 2-10 keV flux of a possible PWN of 4.7e-14 erg/s/cm^2, three times less than the previously claimed value, implying an efficiency for conversion of spin-down energy into nebular luminosity of <9e-4 (assuming a distance of 4 kpc). We do, however, find strong evidence for X-ray emission from the supernova remnant shell surrounding the pulsar, as previously reported.

We have recently carried out a computational campaign to investigate a model of coronal heating in three-dimensions using reduced magnetohydrodynamics (RMHD). Our code is built on a conventional scheme using the pseudo-spectral method, and is parallelized using MPI. The current investigation requires very long time integrations using high Lundquist numbers, where the formation of very fine current layers challenge the resolutions achievable even on massively parallel machines. We present here results of a port to Nvidia CUDA (Compute Unified Device Architecture) for hardware acceleration using graphics processing units (GPUs). In addition to a brief discussion of our general strategy, we will report code performance on several machines which span a variety of hardware configurations and capabilities. These include a desktop workstation with commodity hardware, a dedicated research workstation equipped with four Nvidia C2050 GPUs, as well as several large-scale GPU accelerated distributed memory machines: Lincoln/NCSA, Dirac/NERSC, and Keeneland/NICS.

On 2011 July 14, a new magnetar candidate, Swift J1822.3-1606, was identified via a rate trigger on the Swift/Burst Alert Telescope. Here we present an initial analysis of the X-ray properties of the source, using data from the Rossi X-ray Timing Explorer, Swift, and the Chandra X-ray Observatory, spanning 2011 July 16--September 22. We measure a precise spin period of P=8.43771963(5) s and a spin-down rate of 2.97(28)E-13, at MJD 55761.0, corresponding to an inferred surface dipole magnetic field strength of B=5.1E13 G, the second lowest thus far measured for a magnetar, though similar to 1E~2259+586 as well as to several high-magnetic field radio pulsars. We show that the pulsed X-ray flux decay in the 2--10 keV band is best fit by an exponential with a time constant of 16.4+/-0.3 days. After increasing from ~35% during the first week after the onset of the outburst, the pulsed fraction in the 2--10 keV band remained constant at ~45. We argue that these properties confirm this source to be a new member of the class of objects known as magnetars.

Electrostatic structures have been observed in many regions of space plasmas, including the solar wind, the magnetosphere, the auroral acceleration region, and in association with shocks, turbulence, and magnetic reconnection. Due to potentially large amplitude of electric fields within these structures, their effects on particle heating, scattering, or acceleration can be important. One possible theoretical description of some of these structures is the concept of Bernstein-Greene-Kruskal (BGK) modes, which are exact nonlinear solutions of the Vlasov-Poisson system of equations in collisionless kinetic theory. BGK modes have been studied extensively for many decades, predominately in one dimension (1D), although there have been observations showing that some of these structures have clear 3D features. While there have been approximate solutions of higher dimensional BGK modes, an exact 3D BGK mode solution in a finite magnetic field has not been found yet. Recently we have constructed exact solutions of 2D BGK modes in a magnetized plasma with finite magnetic field strength in order to gain insights of the ultimate 3D theory [Ng, Bhattacharjee, and Skiff, Phys. Plasmas 13, 055903 (2006)]. Based on the analytic form of these solutions, as well as Particle-in-Cell (PIC) simulations, we will present numerical studies of their stability for different levels of background magnetic field strength.

An analytical and numerical treatment is given of a constrained version of the tectonics model developed by Priest, Heyvaerts, & Title [2002]. We begin with an initial uniform magnetic field ${\bf B} = B_0 \hat{\bf z}$ that is line-tied at the surfaces $z = 0$ and $z = L$. This initial configuration is twisted by photospheric footpoint motion that is assumed to depend on only one coordinate ($x$) transverse to the initial magnetic field. The geometric constraints imposed by our assumption precludes the occurrence of reconnection and secondary instabilities, but enables us to follow for long times the dissipation of energy due to the effects of resistivity and viscosity. In this limit, we demonstrate that when the coherence time of random photospheric footpoint motion is much smaller by several orders of magnitude compared with the resistive diffusion time, the heating due to Ohmic and viscous dissipation becomes independent of the resistivity of the plasma. Furthermore, we obtain scaling relations that suggest that even if reconnection and/or secondary instabilities were to limit the build-up of magnetic energy in such a model, the overall heating rate will still be independent of the resistivity.

The solar wind and the interstellar medium are permeated by large-scale magnetic fields that render magnetohydrodynamic (MHD) turbulence anisotropic. In the weak-turbulence limit in which three-wave interactions dominate, analytical and high-resolution numerical results based on random scattering of shear-Alfv\'en waves propagating parallel to a large-scale magnetic field, as well as direct simulations demonstrate rigorously an anisotropic energy spectrum that scales as $k^{-2}_\perp$, instead of the famous Iroshnikov-Kraichnan (IK) spectrum of $k^{-3/2}$ for the isotropic case. Even in the absence of a background magnetic field, anisotropy is found to develop with respect to the local magnetic field, although the energy spectrum is globally isotropic and is found to be consistent with a $k^{-3/2}$ scaling. It is also found in direct numerical simulations that the energy cascade rate is much closer to IK scaling than a Kolmogorov scaling. Recent observations in the solar wind on cascade rates (as functions of the proton temperature and solar wind speed at 1 AU) seem to support this result [Vasquez et al. 2007].

High-magnetic-field pulsars represent an important class of objects for studying the relationship between magnetars and radio pulsars. Here we report on four Chandra observations of the high-magnetic-field pulsar J1718-3718 ($B=7.4\times10^{13}$\,G) taken in 2009 as well as on a re-analysis of 2002 Chandra observations of the region. We also report an improved radio position for this pulsar based on ATCA observations. We detect X-ray pulsations at the pulsar's period in the 2009 data, with pulsed fraction 52%$\pm$13% in the 0.8--2.0 keV band. We find that the X-ray pulse is aligned with the radio pulse. The data from 2002 and 2009 show consistent spectra and fluxes; a merged overall spectrum is well fit by a blackbody of temperature $186^{+19}_{-18}$\,eV, slightly higher than predicted by standard cooling models, however, the best-fit neutron star atmosphere model is consistent with standard cooling. We find the bolometric luminosity $L^{\infty}_{\rm bb}=4^{+5}_{-2} \times10^{32}$\,erg\,s$^{-1}\sim0.3\dot{E}$, for a distance of 4.5 kpc. We compile measurements of the temperatures of all X-ray detected high-$B$ pulsars as well as those of low-$B$ radio pulsars and find evidence for the former being on average hotter than the latter.

Recently, secondary island formation due to the tearing instability of the Sweet-Parker current sheet was identified as a possible mechanism that can lead to fast reconnection (less sensitive dependence on Lundquist number $S$) both in numerical simulations using Particle-in-Cell (PIC) method [Daughton et al. 2009], as well as using resistive magnetohydrodynamics (MHD) [Lapenta 2008; Bhattacharjee et al. 2009]. This instability is thought to appear when $S$ is greater than a certain threshold. These recent results prompt us to perform more resistive MHD simulations of a basic reconnection configuration based on the island coalescence instability, using much higher resolutions and larger $S$. Our simulations are based on a fairly standard pseudo spectral code, which has been tested for accuracy, convergence, and compared well with codes using other methods [Ng et al. 2008]. In our simulations, formation of plasmoids were not found, except when insufficient resolution was used, or when a small amount of noise was added externally. The reconnection rate is found to follow the Sweet-Parker scaling when no noise is added, but increases to a level independent of $S$ with noise, when plasmoids form. Latest results with $S$ up to $2 \times 10^5$ will be presented.

Whether the phenomenology governing MHD turbulence is Kolmogorov or Iroshnikov-Kraichnan (IK) remains an open question, theoretically as well as observationally. The ion heating profile observed in the solar wind provides a quantitative, if indirect, observational constraint on the relevant phenomenology. Recently, a solar wind heating model based on Kolmogorov spectral scaling has produced reasonably good agreement with observations, provided the effect of turbulence generation due to pickup ions is included in the model. Without including the pickup ion contributions, the Kolmogorov scaling predicts a proton temperature profile that decays too rapidly beyond a radial distance of 15 AU. In the present study, we alter the heating model by applying an energy cascade rate based on IK scaling, and show that the model yields higher proton temperatures, within the range of observations, with or without the inclusion of the effect due to pickup ions. Furthermore, the turbulence correlation length based on IK scaling seems to follow the trend of observations better.

In recent years, a phenomenological solar wind heating model based on a turbulent energy cascade prescribed by the Kolmogorov theory has produced reasonably good agreement with observations on proton temperatures out to distances around 70 AU, provided the effect of turbulence generation due to pickup ions is included in the model. In a recent study [Ng et al., J. Geophys. Res., 115, A02101 (2010)], we have incorporated in the heating model the energy cascade rate based on Iroshnikov-Kraichnan (IK) scaling. We showed that the IK cascade rate can also produce good agreement with observations, with or without the inclusion of pickup ions. This effect was confirmed both by integrating the model using average boundary conditions at 1 AU, and by applying a method [Smith et al., Astrophys. J., 638, 508 (2006)] that uses directly observed values as boundary conditions. The effects due to pickup ions is found to be less important for the IK spectrum, which is shallower than the Kolmogorov spectrum. In this paper, we will present calculations of the pickup ions effect in more details, and discuss the physical reason why a shallower spectrum generates less waves and turbulence.

We present a new Chandra observation of supernova remnant (SNR) G350.1-0.3. The high resolution X-ray data reveal previously unresolved filamentary structures and allow us to perform detailed spectroscopy in the diffuse regions of this SNR. Spectral analysis demonstrates that the region of brightest emission is dominated by hot, metal-rich ejecta while the ambient material along the perimeter of the ejecta region and throughout the remnant's western half is mostly low-temperature, shocked interstellar/circumstellar medium (ISM/CSM) with solar-type composition. The data reveal that the emission extends far to the west of the ejecta region and imply a lower limit of 6.6 pc on the diameter of the source (at a distance of 4.5 kpc). We show that G350.1-0.3 is likely in the free expansion (ejecta-dominated) stage and calculate an age of 600-1200 years. The derived relationship between the shock velocity and the electron/proton temperature ratio is found to be entirely consistent with that of other SNRs. We perform spectral fits on the X-ray source XMMU J172054.5-372652, a candidate central compact object (CCO), and find that its spectral properties fall within the typical range of other CCOs. We also present archival 24 um data of G350.1-0.3 taken with the Spitzer Space Telescope during the MIPSGAL galactic survey and find that the infrared and X-ray morphologies are well-correlated. These results help to explain this remnant's peculiar asymmetries and shed new light on its dynamics and evolution.

The ~800 yr-old pulsar J1846-0258 is a unique transition object between rotation-powered pulsars and magnetars: though behaving like a rotation-powered pulsar most of the time, in 2006 it exhibited a distinctly magnetar-like outburst accompanied by a large glitch. Here we present X-ray timing observations taken with the Rossi X-ray Timing Explorer over a 2.2-yr period after the X-ray outburst and glitch had recovered. We observe that the braking index of the pulsar, previously measured to be n=2.65+/-0.01, is now n=2.16+/-0.13, a decrease of 18+/-5%. We also note a persistent increase in the timing noise relative to the pre-outburst level. Despite the timing changes, a 2009 Chandra X-ray Observatory observation shows that the X-ray flux and spectrum of the pulsar and its wind nebula are consistent with the quiescent levels observed in 2000.

We present a detailed analysis of the first very long baseline interferometry (VLBI) detection of the radio remnant of supernova 1987A. The VLBI data taken in 2007 and 2008 at 1.4 and 1.7 GHz, respectively, provide images sensitive to angular scales from 0.1" to 0.7", the highest resolution to date at radio frequencies. The results reveal two extended lobes with an overall morphology consistent with observations at lower resolutions. We find evidence of small-scale features in the radio shell, which possibly consist of compact clumps near the inner surface of the shell. These features have angular extent smaller than 0.2" and contribute less than 13% of the total remnant flux density. No central source is detected in the VLBI images. We place a 3-sigma flux density limit of 0.3 mJy on any pulsar or pulsar wind nebula at 1.7 GHz.

We present results from observations of the magnetar 1E 1547.0-5408 (SGR J1550-5418) taken with the Chandra X-ray Observatory and the Rossi X-ray Timing Explorer (RXTE) following the source's outbursts in 2008 October and 2009 January. During the time span of the Chandra observations, which covers days 4 through 23 and days 2 through 16 after the 2008 and 2009 events, respectively, the source spectral shape remained stable, while the pulsar's spin-down rate in the same span in 2008 increased by a factor of 2.2 as measured by RXTE. The lack of spectral variation suggests decoupling between magnetar spin-down and radiative changes, hence between the spin-down-inferred magnetic field strength and that inferred spectrally. We also found a strong anti-correlation between the phase-averaged flux and the pulsed fraction in the 2008 and 2009 Chandra data, but not in the pre-2008 measurements. We discuss these results in the context of the magnetar model.

Anomalous X-ray pulsars and soft gamma repeaters have recently emerged as a unified class of neutron stars, identified by dramatic X-ray and gamma-ray outbursts and via luminous X-ray pulsations, both thought to be powered by the decay of an enormous internal magnetic field. This "magnetar" hypothesis has raised the question of these objects' physical relationship with conventional rotation-powered pulsars (RPPs). The highest magnetic-field RPPs might therefore be expected to be transition objects between the two populations. The recently reported magnetar-like outburst of PSR J1846-0258, previously thought to be purely rotation-powered, clearly supports this suggestion. Here we review the observational properties of the highest magnetic-field RPPs known, and show some common characteristics that are notable among RPPs, which are plausibly related to their high fields. Using these objects, we consider the evidence for proposed "magneto-thermal evolution" in neutron stars, and argue that while some exists, it is not yet conclusive.

We analysed the XMM-Newton archival observations of 16 neutron star (NS) low-mass X-ray binaries (LMXBs) to study the Fe K emission in these objects. The sample includes all the observations of NS LMXBs performed in EPIC pn Timing mode with XMM-Newton publicly available until September 30, 2009. We performed a detailed data analysis considering pile-up and background effects. The properties of the iron lines differed from previous published analyses due to either incorrect pile-up corrections or different continuum parameterization. 80% of the observations for which a spectrum can be extracted showed significant Fe line emission. We found an average line centroid of 6.67 $\pm$ 0.02 keV and a finite width, $\sigma$, of 0.33 $\pm$ 0.02 keV. The equivalent width of the lines varied between 17 and 189 eV, with an average weighted value of 42 $\pm$ eV. For sources where several observations were available the Fe line parameters changed between observations whenever the continuum changed significantly. The line parameters did not show any correlation with luminosity. Most important, we could fit the Fe line with a simple Gaussian component for all the sources. The lines did not show the asymmetric profiles that were interpreted as an indication of relativistic effects in previous analyses of these LMXBs.

We report radio polarization observations of G319.9-0.7 (MSC 319.9-0.7) at 3 and 6 cm obtained with the Australia Telescope Compact Array. The source shows a highly elongated morphology with the energetic pulsar J1509-5850 located at the tip. We found a flat radio spectrum of index \alpha=-0.26 +/- 0.04 and a high degree of linear polarization. These results confirm G319.9-0.7 as a bow-shock pulsar wind nebula. The polarization maps suggest a helical magnetic field trailing the pulsar, with the symmetry axis parallel to the system's inferred direction of motion. This is the first time such a field geometry has been seen in a bow-shock nebula, and it may be the result of an alignment between the pulsar spin axis and its space velocity. Compared to other bow-shock examples, G319.9-0.7 exhibits very different properties in the field structure and surface brightness distribution, illustrating the large diversity of the population.

We present spectropolarimetric radio images of the supernova remnant (SNR) G296.5+10.0 at frequencies near 1.4 GHz, observed with the Australia Telescope Compact Array. By applying rotation measure (RM) synthesis to the data, a pixel-by-pixel map of Faraday rotation has been produced for the entire remnant. We find G296.5+10.0 to have a highly ordered RM structure, with mainly positive RMs (mean RM of +28 rad/m**2) on the eastern side and negative RMs (mean RM of -14 rad/m**2) on the western side, indicating a magnetic field which is directed away from us on one side and toward us on the other. We consider several possible mechanisms for creating the observed RM pattern. Neither Faraday rotation in foreground interstellar gas nor in a homogeneous ambient medium swept up by the SNR shell can easily explain the magnitude and sign of the observed RM pattern. Instead, we propose that the observed RMs are the imprint of an azimuthal magnetic field in the stellar wind of the progenitor star. Specifically, we calculate that a swept-up magnetized wind from a red supergiant can produce RMs of the observed magnitude, while the azimuthal pattern of the magnetic field at large distances from the star naturally produces the anti-symmetric RM pattern observed. Expansion into such a wind can possibly also account for the striking bilateral symmetry of the SNR's radio and X-ray morphologies.

We present extensive observations of the radio emission from the remnant of SN 1987A made with the Australia Telescope Compact Array (ATCA), since the first detection of the remnant in 1990. The radio emission has evolved in time providing unique information on the interaction of the supernova shock with the circumstellar medium. We particularly focus on the monitoring observations at 1.4, 2.4, 4.8 and 8.6 GHz, which have been made at intervals of 4-6 weeks. The flux density data show that the remnant brightness is now increasing exponentially, while the radio spectrum is flattening. The current spectral index value of -0.68 represents an 18+/-3% increase over the last 8 years. The exponential trend in the flux is also found in the ATCA imaging observations at 9 GHz, which have been made since 1992, approximately twice a year, as well as in the 843 MHz data set from the Molonglo Observatory Synthesis Telescope from 1987 to March 2007. Comparisons with data at different wavelengths (X-ray, H\alpha) are made. The rich data set that has been assembled in the last 22 years forms a basis for a better understanding of the evolution of the supernova remnant.

We report observations of the remnant of Supernova 1987A with the High Resolution Camera (HRC) onboard the Chandra X-ray Observatory. A direct image from the HRC resolves the annular structure of the X-ray remnant, confirming the morphology previously inferred by deconvolution of lower resolution data from the Advanced CCD Imaging Spectrometer. Detailed spatial modeling shows that the a thin ring plus a thin shell gives statistically the best description of the overall remnant structure, and suggests an outer radius 0.96" +/- 0.05" +/- 0.03" for the X-ray-emitting region, with the two uncertainties corresponding to the statistical and systematic errors, respectively. This is very similar to the radius determined by a similar modeling technique for the radio shell at a comparable epoch, in contrast to previous claims that the remnant is 10-15% smaller at X-rays than in the radio band. The HRC observations put a flux limit of 0.010 cts/s (99% confidence level, 0.08-10 keV range) on any compact source at the remnant center. Assuming the same foreground neutral hydrogen column density as towards the remnant, this allows us to rule out an unobscured neutron star with surface temperature T^\infty > 2.5MK observed at infinity, a bright pulsar wind nebula or a magnetar.

The aftermath of supernova (SN) 1987A continues to provide spectacular insights into the interaction between a SN blastwave and its circumstellar en- vironment. We here present 36 GHz observations from the Australia Telescope Compact Array of the radio remnant of SN 1987A. These new images, taken in 2008 Apr and 2008 Oct, substantially extend the frequency range of an ongo- ing monitoring and imaging program conducted between 1.4 and 20 GHz. Our 36.2 GHz images have a diffraction-limited angular resolution of 0.3-0.4 arcseconds, which covers the gap between high resolution, low dynamic range VLBI images of the remnant and low resolution, high dynamic range images at frequencies between 1 and 20 GHz. The radio morphology of the remnant at 36 GHz is an elliptical ring with enhanced emission on the eastern and western sides, similar to that seen previously at lower frequencies. Model fits to the data in the Fourier domain show that the emitting region is consistent with a thick inclined torus of mean radius 0.85 arcsec, and a 2008 Oct flux density of 27 +/- 6 mJy at 36.2 GHz. The spectral index for the remnant at this epoch, determined between 1.4 GHz and 36.2 GHz, is -0.83. There is tentative evidence for an unresolved central source with flatter spectral index.

The faint radio supernova remnant SNR G315.9-0.0 is notable for a long and thin trail that extends outward perpendicular from the edge of its approximately circular shell. In a search with the Parkes telescope we have found a young and energetic pulsar that is located at the tip of this collimated linear structure. PSR J1437-5959 has period P = 61 ms, characteristic age tau_c = 114 kyr, and spin-down luminosity dE/dt = 1.4e36 erg/s. It is very faint, with a flux density at 1.4 GHz of about 75 uJy. From its dispersion measure of 549 pc/cc, we infer d ~ 8 kpc. At this distance and for an age comparable to tau_c, the implied pulsar velocity in the plane of the sky is V_t = 300 km/s for a birth at the center of the SNR, although it is possible that the SNR/pulsar system is younger than tau_c and that V_t > 300 km/s. The highly collimated linear feature is evidently the pulsar wind trail left from the supersonic passage of PSR J1437-5959 through the interstellar medium surrounding SNR G315.9-0.0.

Radio continuum emission from the supernova remnant G296.5+10.0 was observed using the Australia Telescope Compact Array. Using a 104 MHz bandwidth split into 13 x 8 MHz spectral channels, it was possible to produce a pixel-by-pixel image of Rotation Measure (RM) across the entire remnant. A lack of correlation between RM and X-ray surface brightness reveals that the RMs originate from outside the remnant. Using this information, we will characterise the smooth component of the magnetic field within the supernova remnant and attempt to probe the magneto-ionic structure and turbulent scale sizes in the ISM and galactic halo along the line-of-sight.

We present the results of detailed spatial and spectral analysis of the pulsar wind nebula (PWN) in supernova remnant Kes 75 (G29.7-0.3) using a deep exposure with Chandra X-ray observatory. The PWN shows a complex morphology with clear axisymmetric structure. We identified a one-sided jet and two bright clumps aligned with the overall nebular elongation, and an arc-like feature perpendicular to the jet direction. Further spatial modeling with a torus and jet model indicates a position angle $207\arcdeg\pm8 \arcdeg$ for the PWN symmetry axis. We interpret the arc as an equatorial torus or wisp and the clumps could be shock interaction between the jets and the surrounding medium. The lack of any observable counter jet implies a flow velocity larger than 0.4c. Comparing to an archival observation 6 years earlier, some small-scale features in the PWN demonstrate strong variability: the flux of the inner jet doubles and the peak of the northern clump broadens and shifts 2" outward. In addition, the pulsar flux increases by 6 times, showing substantial spectral softening from $\Gamma$=1.1 to 1.9 and an emerging thermal component which was not observed in the first epoch. The changes in the pulsar spectrum are likely related to the magnetar-like bursts of the pulsar that occurred 7 days before the Chandra observation, as recently reported from RXTE observations.

We present detailed Fourier modeling of the radio remnant of Supernova 1987A, using high-resolution 9 GHz and 18 GHz data taken with the Australia Telescope Compact Array over the period 1992 to 2008. We develop a parameterized three-dimensional torus model for the expanding radio shell, in which the emission is confined to an inclined equatorial belt; our model also incorporates both a correction for light travel-time effects and an overall east-west gradient in the radio emissivity. By deriving an analytic expression for the two-dimensional Fourier transform of the projected three-dimensional brightness distribution, we can fit our spatial model directly to the interferometric visibility data. This provides robust estimates to the radio morphology at each epoch. The best-fit results suggest a constant remnant expansion at 4000 +/- 400 km/s over the 16-year period covered by the observations. The model fits also indicate substantial mid-latitude emission, extending to 40 degree on either side of the equatorial plane. This likely corresponds to the extra-planar structure seen in H$\alpha$ and Ly$\alpha$ emission from the supernova reverse shock, and broadly supports hydrodynamic models in which the complex circumstellar environment was produced by a progression of interacting winds from the progenitor. Our model quantifies the clear asymmetry seen in the radio images: we find that the eastern half of the radio remnant is consistently ~40 brighter than the western half at all epochs, which may result from an asymmetry in the ejecta distribution between these two hemispheres.

We present an XMM-Newton observation of the long-overlooked radio source G350.1-0.3. The X-ray spectrum of G350.1-0.3 can be fit by a shocked plasma with two components: a high-temperature (1.5 keV) region with a low ionization time scale and enhanced abundances, plus a cooler (0.36 keV) component in ionization equilibrium and with solar abundances. The X-ray spectrum and the presence of non-thermal, polarized, radio emission together demonstrate that G350.1-0.3 is a young, luminous supernova remnant (SNR), for which archival HI and 12-CO data indicate a distance of 4.5 kpc. The diameter of the source then implies an age of only ~900 years. The SNR's distorted appearance, small size and the presence of 12-CO emission along the SNR's eastern edge all indicate that the source is interacting with a complicated distribution of dense ambient material. An unresolved X-ray source, XMMU J172054.5-372652, is detected a few arcminutes west of the brightest SNR emission. The thermal X-ray spectrum and lack of any multi-wavelength counterpart suggest that this source is a neutron star associated with G350.1-0.3, most likely a "central compact object", as seen coincident with other young SNRs such as Cassiopeia A.

We describe recent Chandra ACIS observations of the Vela-like pulsar PSR J2021+3651 and its pulsar wind nebula (PWN). This `Dragonfly Nebula' displays an axisymmetric morphology, with bright inner jets, a double-ridged inner nebula, and a ~30" polar jet. The PWN is embedded in faint diffuse emission: a bow shock-like structure with standoff ~1' brackets the pulsar to the east and emission trails off westward for 3-4'. Thermal (kT=0.16 +/-0.02 keV) and power law emission are detected from the pulsar. The nebular X-rays show spectral steepening from Gamma=1.5 in the equatorial torus to Gamma=1.9 in the outer nebula, suggesting synchrotron burn-off. A fit to the `Dragonfly' structure suggests a large (86 +/-1 degree) inclination with a double equatorial torus. Vela is currently the only other PWN showing such double structure. The >12 kpc distance implied by the pulsar dispersion measure is not supported by the X-ray data; spectral, scale and efficiency arguments suggest a more modest 3-4 kpc.

A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.

We have applied the torus fitting procedure described in Ng & Romani (2004) to PWNe observations in the Chandra data archive. This study provides quantitative measurement of the PWN geometry and we characterize the uncertainties in the fits, with statistical errors coming from the fit uncertainties and systematic errors estimated by varying the assumed fitting model. The symmetry axis $\Psi$ of the PWN are generally well determined, and highly model-independent. We often derive a robust value for the spin inclination $\zeta$. We briefly discuss the utility of these results in comparison with new radio and high energy pulse measurements

Evidence from pulsar wind nebula symmetry axes and radio polarization observations suggests that pulsar motions correlate with the spin directions. We assemble this evidence for young isolated pulsars and show how it can be used to quantitatively constrain birth kick scenarios. We illustrate by computing several plausible, but idealized, models where the momentum thrust is proportional to the neutrino cooling luminosity of the proto-neutron star. Our kick simulations include the effects of pulsar acceleration and spin-up and our maximum likelihood comparison with the data constrains the model parameters. The fit to the pulsar spin and velocity measurements suggests that: i) the anisotropic momentum required amounts to ~10% of the neutrino flux, ii) while a pre-kick spin of the star is required, the preferred magnitude is small 10-20rad/s, so that for the best-fit models iii) the bulk of the spin is kick-induced with $\bar \Omega$ ~120rad/s and iv) the models suggest that the anisotropy emerges on a timescale $\tau$ ~1-3s.

We report on VLBA astrometry and CXO imaging of PSR J0538+2817 in the supernova remnant S147. We measure a parallax distance of $1.47^{+0.42}_{-0.27}$kpc along with a high-precision proper motion, giving a transverse velocity $V_\perp=400^{+114}_{-73}$km/s. A small extended wind nebula is detected around the pulsar; the symmetry axis of this structure suggests that the spin axis lies $12\arcdeg\pm4\arcdeg$ from the velocity vector (2-D), but the emission is too faint for robust model independent statements. The neutron star is hot, consistent with the young ~40kyr kinematic age. The pulsar progenitor is likely a runaway from a nearby cluster, with NGC 1960 (M36) a leading candidate.

It has been suggested that the Crab pulsar's proper motion is well aligned with the symmetry axis of the pulsar wind nebula. We have re-visited this question, examining over 6 years of F547M WFPC2 chip 3 images to obtain a best-fit value of $\mu_\ast = 14.9 \pm 0.8$mas/yr at PA $278\arcdeg \pm 3\arcdeg$. At $26\arcdeg\pm 3\arcdeg$ to the nebula axis, this substantially relaxes constraints on the birth kick of this pulsar. Such misalignment allows the momentum to be imparted over $\sim$1s timescales.

We report on Chandra ACIS imaging of the pulsar wind nebula (PWN) of the young Vela-like PSR B1706-44, which shows the now common pattern of an equatorial wind and polar jets. The structure is particularly rich, showing a relativistically boosted termination shock, jets with strong confinement, a surrounding radio/X-ray PWN and evidence for a quasi-static `bubble nebula'. The structures trace the pulsar spin geometry and illuminate its possible relation to SNR G343.1-2.3. We also obtain improved estimates of the pulsar flux and nebular spectrum, constraining the system age and energetics.

We report on Chandra ACIS and XMM-Newton MOS/PN imaging observations of two pulsar wind nebulae (K3/PSR J1420-6048 and G313.3+0.1=`the rabbit') associated with the Galactic unidentified gamma-ray source GeV J1417-6100. With the excellent ACIS imaging, the very energetic pulsar PSR J1420-6048 is separated from its surrounding nebula. This nebula has surprisingly little compact structure, although a faint arc is seen near the pulsar. Similarly, two point sources are resolved in the rabbit nebula. The large XMM-Newton collecting area provides useful spectral constraints on the rabbit and the associated point sources. Based on spectra and X-ray morphology, we identify one point source as a plausible pulsar counterpart. Large backgrounds and low source counts limited pulse search sensitivities, but we report pulse upper limits and a candidate 108ms period for the rabbit pulsar based on the XMM-Newton data and an ACIS CC observation. Comparison of the X-ray images with high resolution ATCA radio maps shows that the non-thermal X-ray emission corresponds well with the radio structure.

A remarkable number of pulsar wind nebulae (PWN) are coincident with EGRET gamma-ray sources. X-ray and radio imaging studies of unidentified EGRET sources have resulted in the discovery of at least 6 new pulsar wind nebulae (PWN). Stationary PWN (SPWN) appear to be associated with steady EGRET sources with hard spectra, typical for gamma-ray pulsars. Their toroidal morphologies can help determine the geometry of the pulsar which is useful for constraining models of pulsed gamma-ray emission. Rapidly moving PWN (RPWN) with more cometary morphologies seem to be associated with variable EGRET sources in regions where the ambient medium is dense compared to what is typical for the ISM.

We present results from X-ray and radio observations of the recently discovered young Vela-like pulsar PSR J2021+3651, which is coincident with the EGRET gamma-ray source GeV 2020+3658. A 19.0-ks Chandra ACIS-S observation has revealed a ~20'' x 10'' pulsar wind nebula that is reminiscent of the equatorial tori seen around some young pulsars, along with thermal emission from an embedded point source (kT = 0.15 +/- 0.02 keV). We name the nebula G75.2+0.1. Its spectrum is well fit by an absorbed power-law model with photon index 1.7 +/- 0.3, hydrogen column density nH = 7.8 +/- 1.7 x 10^21 cm^-2, and an unabsorbed 0.3-10.0 keV flux of 1.9 +/- 0.3 x 10^-12 erg cm^-2 s^-1. We have spatially fit G75.2+0.1 with a model that assumes a toroidal morphology, and from this we infer that the torus is highly inclined 83 deg +/- 1 deg to the line of sight. A 20.8-ks Chandra observation in continuous-clocking mode reveals a possible pulse detection, with a pulsed fraction of ~37% and an H-test probability of occuring by chance of 1.2 x 10^-4. Timing observations with the Arecibo radio telescope spanning two years show that PSR J2021+3651 glitched sometime between MJDs 52616 and 52645 with parameters delta(v)/v = (2.587 +/- 0.002) x 10^-6 and delta(dot(v))/v = (6.2 +/- 0.3) x 10^-3, similar to those of the largest glitches observed in the Vela pulsar. PSR J2021+3651 is heavily scattered (T_sc = 17.7 ms +/- 0.9 ms at 1 GHz) and exhibits a significant amount of timing noise.

CXO imaging has shown that equatorial tori, often with polar jets, are very common in young pulsar wind nebulae (PWNe). These structures are interesting both for what they reveal about the relativistic wind itself and for the (nearly) model-independent measurement of the neutron star spin orientation they provide. The later is a particularly valuable probe of pulsar emission models and of neutron star physics.We describe here a procedure for fitting simple 3-D torus models to the X-ray data which provides robust estimates of the geometric parameters. An application to 6 PWN tori gives orientations, PWN shock scales and post-shock wind speeds along with statistical errors. We illustrate the use of these data by commenting on the implications for kick physics and for high energy beaming models.

We find evidence for a faint wind nebula surrounding PSR J0538+2817 in CXO-ACIS imaging. This object is particularly interesting, as the pulsar spindown age is largest for any such X-ray PWN. If interpreted as an equatorial torus, the PWN supports the claimed association with the S147 supernova remnant and implies good alignment between the pulsar spin and space velocity. Comparison of the SNR age, X-ray cooling age and characteristic age suggests a birth spin period of >130ms. In turn, if we accept as causal the alignment of the linear and angular momenta, this places strong constraints on the origin of the `kick' at the neutron star birth.

As the chiral symmetry is widely recognized as an important driver of the strong interaction dynamics, current strange stars models based on MIT bag models do not obey such symmetry. We investigate properties of bare strange stars using the Cloudy Bag Model, in which a pion cloud coupled to the quark-confining bag is introduced such that chiral symmetry is conserved. We find that in this model the decay of pions is a very efficient cooling way. In fact it can carry out most the thermal energy in a few milliseconds and directly convert them into 100MeV photons via pion decay. This may be a very efficient $\gamma$-ray burst mechanism. Furthermore, the cooling behavior may provide a possible way to distinguish a compact object between a neutron star, MIT strange star and Cloudy Bag strange star in observations.

We evaluate the conductivity tensor for molecular gas at densities ranging from 10^4 to 10^15 cm^-3 for a variety of grain models. The Hall contribution to the conductivity has generally been neglected in treatments of the dynamics of molecular gas. We find that it is not important if only 0.1 micron grains are considered, but for a Mathis-Rumpl-Nordsieck grain-size distribution (with or without PAHs) it becomes important for densities between 10^7 and 10^11 cm^-3. If PAHs are included, this range is reduced to 10^9 -- 10^10 cm^-3. The consequences for the magnetic field evolution and dynamics of dense molecular gas are profound. To illustrate this, we consider the propagation of Alfven waves under these conditions. A linear analysis yields a dispersion relation valid for frequencies below the neutral collision frequencies of the charged species. The dispersion relation shows that there is a pair of circularly polarised modes with distinct propagation speeds and damping rates. We note that the gravitational collapse of dense cloud cores may be substantially modified by the Hall term.