In this proceeding, we are going to introduce our recent work about the dependence of pulsar death line on the equation of state (Zhou et al. 2017). The results show that the equation of state affects the location of the pulsar death line in the P-Ṗ diagram. We offer another point of view to understand the pulsar death line and multiple observational facts would help us to reveal the nature of pulsars.

Pulsar polarization has been a fruitful area of study since the first discovery of pulsars 50 years ago. Polarization gives information on the geometry of the star, the location of the radio emission in the magnetosphere, the physics behind the radio emission mechanism and a plethora of phenomenology. Here, I will restrict myself to a brief outline of recent work in pulsar polarization using observations taken with the Parkes radio telescope over the past decade.

Sub-pulse drift related profile mode-changes provide a useful probe of pulsar radio emission. Measurements on modal profiles of PSRs B0031–07, J1822–2256 and B2319+60 are presented for different drift modes in this communication. The width of profile increases with the drift rate, while no such trend is seen for the pulsed flux density for all these pulsars. A brief discussion of implications for models is presented.

In-spiraling supermassive black holes should emit gravitational waves, which would produce characteristic distortions in the time of arrival residuals from millisecond pulsars. Multiple national and regional consortia have constructed pulsar timing arrays by precise timing of different sets of millisecond pulsars. An essential aspect of precision timing is the transfer of the times of arrival to a (quasi-)inertial frame, conventionally the solar system barycenter. The barycenter is determined from the knowledge of the planetary masses and orbits, which has been refined over the past 50 years by multiple spacecraft. Within the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), uncertainties on the solar system barycenter are emerging as an important element of the NANOGrav noise budget. We describe what is known about the solar system barycenter, touch upon how uncertainties in it affect gravitational wave studies with pulsar timing arrays, and consider future trends in spacecraft navigation.

Radio continuum surveys are equally sensitive to all pulsars, not affected by dispersion measure smearing, scattering or orbital modulation of spin periods, and therefore allow us to search for extreme pulsars, such as sub-millisecond pulsars, pulsar-black hole systems and pulsars in the Galactic Centre. As we move towards the Square Kilometre Array (SKA) era, searching for pulsars in continuum images will complement conventional pulsar searches, and make it possible to find extreme objects.

Ejection activities in S255IR-SMA1 and AFGL 5142 were investigated by multi-epoch VLBI observations of 22 GHz water masers, tracing bowshocks leading collimated jets. The history of ejections, revealed by the 3D maser motions and supplemented by the literature, suggests that these massive stars formed by episodic accretion, inferred via the accretion-ejection connection. This contribution centers on the role of episodic accretion in overcoming the radiation pressure problem of massive star formation - with maser VLBI and single-dish observations providing essential observational tools.

Radio astrometric campaigns using VLBI have provided distances and proper motions for masers associated with young massive stars (BeSSeL survey). The ongoing BAaDE project plans to obtain astrometric information of SiO maser stars located in the inner Galaxy. These stars are associated with evolved, mass-losing stars. By overlapping optical (Gaia), infrared (2MASS, MSX and WISE) and radio (BAaDE) sources, we expect to obtain important clues on the intrinsic properties and population distribution of late-type stars. Moreover, a comparison of the Galactic parameters obtained with Gaia and VLBI can be done using radio observations on different targets: young massive stars (BeSSeL) and evolved stars (BAaDE).

The final stages of low-mass stellar evolution are characterized by significant mass loss due to stellar pulsations during the AGB phase, which lead to the development of planetary nebulae. Molecular masers of H2O, SiO, and ground state OH transitions are commonly detected in oxygen-rich late-type stars (OH/IR objects). In contrast, excited OH maser transitions are rare. We discuss our study of the carbon-rich pre-planetary nebula CRL618 (a prototypical post-AGB star). Observations conducted in May 2008 with the 305m Arecibo Telescope resulted in the first detection of a 4765MHz OH maser line in a late-type stellar object; the detection was confirmed a few months later also with Arecibo. Subsequent observations in 2015 and 2017 resulted in non-detection of the 4765MHz OH line. Our observations indicate that the 4765MHz OH maser in CRL 618 is highly variable, possibly tracing a short-lived phenomenon during the development of a pre-planetary nebula.

Fifty years of pulsars also mean fifty years of using them as tools to probe other phenomena and physics. One prominent example is the usage of pulsars to test theories of gravity. Probing the quasi-stationary strong-field regime, pulsars allow high precision tests that will maintain their importance even in the era of gravitation wave observations with ground-based detectors. This contribution summarise the methods and status of the field and provides a brief outlook into the future.

A statistical study of the glitch population and the behavior of the glitch activity across the known population of neutron stars is presented. A constant ratio between the glitch activity and the spin-down rate $\dot{\nu }_{\rm {g}}$/|$\dot{\nu }$| = 0.010 ± 0.001 is consistent with the behavior of all rotation-powered pulsars and magnetars. This relation is dominated by large glitches (Δν ≳ 10 μ Hz), which occur at a rate directly proportional to |$\dot{\nu }$|. The only exception are the rotation-powered pulsars with the highest values of |$\dot{\nu }$|, such as the Crab pulsar and PSR B0540–69, which exhibit a much smaller glitch activity, intrinsically different from each other and from the rest of the population. This contribution is based on the work done by Fuentes et al. (2017) “The glitch activity of neutron stars”, accepted for publication in A&A.

Maser theory continues to be driven by advances in observational techniques. Here, I consider the responses to VLBI with space-Earth baselines and cross-correlation spectroscopy (a re-consideration of coherence properties), routine observation in full-Stokes polarization (a re-casting of the polarization transfer equations), and long-term variability monitoring (3-D modelling of irregular domains).

We present the results of the first detection of submillimeter water maser emission toward water-fountain nebulae. Using APEX we found emission at 321.226 GHz toward two sources: IRAS 18043−2116, and IRAS 18286−0959. The submillimeter H2O masers exhibit expansion velocities larger than those of the OH masers, suggesting that these masers, similarly to the 22 GHz masers, originate in fast bipolar outflows. The 321 GHz masers in IRAS 18043−2116 and IRAS 18286−0959, which figure among the sources with the fastest H2O masers, span a velocity range similar to that of the 22 GHz masers, indicating that they probably coexist. The intensity of the submillimeter masers is comparable to the 22 GHz masers, implying that the kinetic temperature of the region where the masers originate is Tk>1000 K. We propose a simple model invoking the passage of two shocks through the same gas that creates the conditions for explaining the strong high-velocity 321 GHz masers coexisting with the 22 GHz masers in the same region.

We present the results of the linear polarisation observations of methanol masers at 44 and 95 GHz towards 39 massive star forming regions (Kang et al. 2016). These two lines are observed simultaneously with the 21-m Korean VLBI Network (KVN) telescope in single dish mode. About 60% of the observed showed fractional polarisation of a few percents at least at one of the two transition lines. We note that the linear polarisation of the 44 GHz methanol maser is first detected in this study including single dish and interferometer observations. We find the polarisation properties of these two lines are similar as expected, since they trace similar regions. As a follow-up study, we have carried out the VLBI polarisation observations toward some 44 GHz maser targets using the KVN telescope. We present preliminary VLBI polarisation results of G10.34-0.14, which show consistent polarisation properties in multiple epoch observations.

Radio pulsars have been responsible for many astonishing astrophysical and fundamental physics breakthroughs since their discovery 50 years ago. In this review I will discuss many of the highlights, most of which were only possible because of the provision of large-scale observing facilities. The next 50 years of pulsar astronomy can be very bright, but only if our governments properly plan and fund the infrastructure necessary to enable future discoveries. Being a small sub-field of astronomy places an onus on the pulsar community to have an open-source/open access approach to data, software, and major observing facilities to enable new groups to emerge to keep the field vibrant.

PSR J1913+1102 is a double neutron star system (DNS) discovered in the Pulsar Arecibo L-band Feed Array survey. We have now very precisely measured the rate of advance of periastron for the system and the Einstein delay. From general relativity, this results in precise mass measurements: 1.65 ± 0.05 and 1.24 ± 0.05 M⊙ for the pulsar and neutron-star companion, respectively. This makes PSR J1913+1102 both the most massive double neutron star system known, and the most asymmetric in mass among compact DNS binaries. This asymmetry will allow for stringent limits on the effects of dipolar gravitational-wave radiation, predicted by alternative theories of gravity, as well as insight into heavy-element production from the eventual merger of this system and others like it. Further observations will also tighten constraints on formation and evolution models; this is crucial for understanding the DNS population, for which there are relatively few mass measurements.

In recent years, surprise discoveries of pulsed emission from the Crab and Vela pulsars above 100 GeV have drawn renewed attention to this largely unexplored region of the energy range. In this paper, we discuss example light curves due to curvature emission, with good resolution in the different energy bands. Continued light curve modelling may help to discriminate between different emission mechanisms, as well as constrain the location where emission is produced within the pulsar magnetosphere, including regions beyond the light cylinder.

In 2014 we conducted a survey for 6.7 GHz methanol masers with the Arecibo Telescope toward far infrared sources selected from the Hi-GAL catalog of massive cores. We found a number of sources with weak 6.7 GHz methanol masers, possibly indicating regions in early stages of star formation. Here we describe the results of follow-up observations that were conducted with the Very Large Array in New Mexico to characterize this new population of “weak” 6.7 GHz methanol masers.

We have conducted astrometric observations toward a 22 GHz water maser source associated with the Sgr B2 complex in the Galactic center region with VERA (VLBI exploration of Radio Astrometry). We measured a trigonometric parallax and absolute proper motion of the Sgr B2 complex with respect to an extra-galactic source by observing the water maser source at 10 epochs from 2014 to 2017. The measured distance was 7.52+3.01−1.67 kpc for the Sgr B2M region.

We also succeeded to measure internal motions of maser spots in Sgr B2M, and N region. The number of spots which we could measure the internal motions is about 400. The distribution of the maser spots shows that the maser spots are associated with envelope of HII region seen in radio continuum image obtained with VLA and ALMA. We discuss relative motions between Sgr B2M, and N by using the internal motion.

The brightness of maser features are fascinating and give valuable insight for circumstellar physics of oxygen-rich, intermediate-mass stars, in particular the final evolution of circumstellar envelopes (CSEs). The variety of accompanying masers such as SiO, H2O, and OH in the CSEs may provide unique probes into different stages of rapid CSE evolution. However, with only sparse monitoring of these masers one can sometimes find it difficult to accurately interpret their spatio-kinematics, origins and excitation mechanisms. Examples can be seen in the variety of proposed models for water masers associated with “water fountains” and for silicon-monoxide masers. In order to better understand these issues, one needs to consider continuous monitoring of the individual maser gas clumps over a few stellar cycles or episodic ejection events. Here I present our previous long-term monitoring observations, especially for the water fountain source W43A. Our current efforts involve programs of intensive monitoring observations of circumstellar maser sources over decadal time periods. These programs with the East Asia VLBI Network observe H2O and SiO maser lines simultaneously mapped at high cadence (2–8 weeks) with VLBI observations.