Since the discovery of the first radio pulsar fifty years ago, the population of neutron stars in our Galaxy has grown to over 2,600. A handful of these sources, exclusively seen in X-rays, show properties that are not observed in normal pulsars. Despite their scarcity, they are key to understanding aspects of the neutron star phenomenology and evolution. The forthcoming all-sky survey of eROSITA will unveil the X-ray faint end of the neutron star population at unprecedented sensitivity; therefore, it has the unique potential to constrain evolutionary models and advance our understanding of the sources that are especially silent in the radio and γ-ray regimes. In this contribution I discuss the expected role of eROSITA, and the challenges it will face, at probing the galactic neutron star population.

Whether high-mass stars (M > 7M⊙) emerge from a scaled-up version of the low-mass star formation scenario, i. e. through disk-mediated accretion, is still debated. We present the first results of an observational programme aimed to map the innermost regions of high-mass stellar objects by combining together high-spatial resolution maser and radio continuum observations, and near-infrared imaging.

Phased VLA observations of the Galactic center magnetar J1745-2900 over 8-12 GHz reveal rich single pulse behavior. The average profile is comprised of several distinct components and is fairly stable over day timescales and GHz frequencies. The average profile is dominated by the jitter of relatively narrow pulses. The pulses in each of the four profile components are uncorrelated in phase and amplitude, although the occurrence of pulse components 1 and 2 appear to be correlated. Using a collection of the brightest individual pulses, we verify that the index of the dispersion law is consistent with the expected cold plasma value of 2. The scattering time is weakly constrained, but consistent with previous measurements, while the dispersion measure DM = 1763+3−10 pc cm−3 is lower than previous measurements, which could be a result of time variability in the line-of-sight column density or changing pulse profile shape over time or frequency.

Using the proper-motion as a proxy for a pulsar’s rotation-axis orientation, the emission polarization can be related to the magnetic field direction in the emission region, providing a much more physical foundation for understanding orthogonal mode radiation.

The observed polarization of radio pulsars involves several peculiar effects, such as comparable amount of orthogonal polarization modes (OPMs) which often bear the same handedness of circular polarisation V. In the average profiles of B1913+16 and B1933+16, orthogonal jumps of polarization angle (PA) are observed to occur at the maximum V, instead of V = 0. High levels of V are also observed in core components (eg. in B1237+25), where they are accompanied by strong distortions of PA from the rotating vector model. In weakly polarized emission, PA jumps by 45° are observed in B1919+21 and B0823+26. It is shown that all these peculiarities can be interpreted in a model which assumes coherent addition of waves in natural propagation modes.

The radial velocity of a pulsar induces the Doppler effect on its intrinsic spin properties. In particular, it can generate a contribution to the frequency second derivative. We estimated this effect for each of the International Pulsar Timing Array pulsars. We also assessed the possibility of measuring the frequency second derivative in the observational data.

In 1975 the carousel model was proposed by Ruderman & Sutherland to explain the beautiful phenomenon of drifting subpulses. However the simultaneous appearance of subpulse bands which drift in opposing directions - a feature now found in two pulsars - is difficult to reconcile with this model, both geometrically and physically. Here we propose a geometric resolution of this problem which also may shed light on a range of previously baffling phenomena. The model places significant constraints on the underlying physics of pulsar emission.

Outstanding problems concerning mass-loss from evolved stars include initial wind acceleration and what determines the clumping scale. Reconstructing physical conditions from maser data has been highly uncertain due to the exponential amplification. ALMA and e-MERLIN now provide image cubes for five H2O maser transitions around VY CMa, at spatial resolutions comparable to the size of individual clouds or better, covering excitation states from 204 to 2360 K. We use the model of Gray et al. 2016, to constrain variations of number density and temperature on scales of a few au, an order of magnitude finer than is possible with thermal lines, comparable to individual cloud sizes or locally almost homogeneous regions. We compare results with the models of Decin et al. 2006 and Matsuura et al. 2014 for the circumstellar envelope of VY CMa; in later work this will be extended to other maser sources.

Pulsed non-thermal quiescent emission between 10 keV and around 150 keV has been observed in ~10 magnetars. For inner magnetospheric models of such hard X-ray signals, resonant Compton upscattering of soft thermal photons from the neutron star surface is the most efficient radiative process. We present angle-dependent hard X-ray upscattering model spectra for uncooled monoenergetic relativistic electrons. The spectral cut-off energies are critically dependent on the observer viewing angles and electron Lorentz factor. We find that electrons with energies less than around 15 MeV will emit most of their radiation below 250 keV, consistent with the observed turnovers in magnetar hard X-ray tails. Moreover, electrons of higher energy still emit most of the radiation below around 1 MeV, except for quasi-equatorial emission locales for select pulses phases. Our spectral computations use new state-of-the-art, spin-dependent formalism for the QED Compton scattering cross section in strong magnetic fields.

We present observations of massive star-forming regions selected from the IRAS Point Source Catalog. The observations were made with the Very Large Array and the Large Millimeter Telescope to search for Class I methanol masers. We made interferometric observations of 125 massive star-forming regions in the 44 GHz methanol maser transition; 53 of the 125 fields showed emission. The data allow us to demonstrate associations, at arcsecond precision, of the Class I maser emission with outflows, HII regions and shocks traced by 4.5 μm emission. We made single-dish observations toward 38 of the 53 regions with 44 GHz masers detected to search for the methanol transitions at 84.5, 95.1, 96.7, 107.0, and 108.8 GHz. We find detection rates of 74, 55, 100, 3, and 45%, respectively. We used a wide-band receiver which revealed many other spectral lines that are common in star-forming regions.

Many accretion disks surrounding supermassive black holes in nearby AGN are observed to host 22 GHz water maser activity. We have analyzed single-dish 22 GHz spectra taken with the GBT to identify 32 such “Keplerian disk systems,” which we used to investigate maser excitation and explore the possibility of disk reverberation. Our results do not support a spiral shock model for population inversion in these disks, and we find that any reverberating signal propagating radially outwards from the AGN must constitute <10% of the total observed maser variability. Additionally, we have used ALMA to begin exploring the variety of sub-mm water megamasers that are also predicted, and in the case of the 321 GHz transition found, to be present in these accretion disks. By observing multiple masing transitions within a single system, we can better constrain the physical conditions (e.g., gas temperature and density) in the accretion disk.

We report another 6.7 GHz methanol maser burst in the high mass star region G33.641-0.228. The flare is in its second component at vLSR = 59.6 km s−1 and was observed in August-September 2016 by VIRAC radio telescope RT-32 in Irbene, Latvia. Several bursts of the second spatial component of G33.641-0.228 have been reported previously by Fujisawa et al. The maximum peak flux density of the source was measured to reach 343 Jy that is 13 times increase from its ground level. Significant oscillations were discovered during the decay phase indicating a more complex burst mechanism that cannot be explained by a simple heating of the region.

The RadioAstron space-VLBI mission has successfully detected extragalactic H2O MegaMaser emission regions at very long Earth to space baselines ranging between 1.4 and 26.7 Earth Diameters (ED). The preliminary results for two galaxies, NGC 3079 and NGC 4258, at baselines longer than one ED indicate masering environments and excitation conditions in these galaxies that are distinctly different. Further observations of NGC 4258 at even longer baselines are expected to reveal more of the physics of individual emission regions.

VLBI observation of masers is a powerful mean to understand the early evolutionary phase of massive star formation. A few different scenarios of outflow evolution in the massive protostars have been proposed, and cannot be readily examined because the precise timing of appropriate maser phenomena is difficult. In particular, it has been a matter of debate whether a well-collimated or a less-collimated outflow comes first in the very early phase of the massive protostellar evolution. Long-term, multi-epoch VLBI monitoring is probably the most important method to trace the outflow evolution. Such a monitoring of a massive star-forming region W75N(B) has been very successful. Since the first detection of the expanding water maser shell associated with the star-forming region VLA 2 of W75N(B) in 1999, the observations in 2005 and 2007 displayed that the expanding water maser shell has been evolved to well-collimated from a less collimated morphology. Observations in 2012 also confirmed such a transition. It would be a major breakthrough in our knowledge of the formation and evolution of the first stages of massive protostars. We performed multi-epoch VLBI observations in mid-2014. On the contrary to its expansion for 13 years, the maser shell at VLA 2 observed in 2014 is comparable to the size observed in 2012. The quenching of the maser shell size indicates that the previously expanding outflow has been decelerated plausibly due to the interaction with surrounding interstellar medium.

The study of astrophysical maser formation provides a useful probe of the chemical composition and physical conditions of the sources they are observed in. This exploration requires continuously solving the SE equations for the populations of the energy levels in search of conditions that will produce an inversion. After evaluation of available implementations applying the Escape Probability approximation, the masers solver was developed to provide an efficient and robust matrix inversion calculation. This open source package is hosted at https://bitbucket.org/ruby_van_rooyen/masers.

The Neutron Star Interior Composition Explorer (NICER) presents an exciting new capability for exploring the modulation properties of X-ray emitting neutron stars, including large area, low background, extremely precise absolute event time stamps, superb low-energy response and flexible scheduling. The Pulsation Searches and Multiwavelength Coordination working group has designed a 2.5 Ms observing program to search for emission and characterize the modulation properties of about 30 known or suspected neutron star sources across a number of source categories. A key early goal will be to search for pulsations from millisecond pulsars that might exhibit thermal pulsations from the surface suitable for pulse profile modeling to constrain the neutron star equation of state. In addition, we will search for pulsations from transitional millisecond pulsars, isolated neutron stars, low-mass X-ray binaries (LMXBs), accretion-powered millisecond pulsars, central compact objects and other sources. We present our science plan and initial results from the first months of the NICER mission, including the discovery of pulsations from the millisecond pulsar J1231–1411.

This talk commented on our progress in understanding high Ė pulsar output, with the photon power dominated by GeV radiation and the total power dominated by the e±/B wind. We are increasingly appreciating the anisotropy in these outflows, with high energy pulsar beaming probed by the distribution of γ-ray pulse profiles and wind anisotropy mapped by synchrotron images of PWNe. Possible hemispheric asymmetry and the prospects for additional probes of pulsar spindown, particularly from compact binaries in the black widow class, are briefly mentioned.

Fifty years of pulsar data has led to the discovery of emission and rotation variability on timescales of months and years; we have developed techniques to identify this long timescale variability. Individual observations may be too noisy to identify subtle changes in a pulse profile; we use Gaussian process regression to model noisy observations and produce a continuous map of pulse profile variability. Generally, multiple observing epochs are required to obtain the pulsar spin frequency derivative. Gaussian process regression is, therefore, also used to monitor this rate of spindown. We have applied variability detection techniques to both millisecond and long period pulsar datasets. I will discuss the techniques used and present the most interesting results from the pulsars analysed.

The aim of our project is to search for ways to best extract information on pulsar profiles and the interstellar medium (ISM), using the wide frequency bands that are typical of radio telescopes today. Pulsar profiles typically show a strong dependence on frequency. This depends both on the intrinsic radio emission mechanism, and the interaction of the radio waves with the ISM that lies between the pulsars and our detectors on Earth, due mostly to the effects of dispersion and scattering. In this work, we make use of radio pulsar beam models from the existing literature, to generate simulated pulse profiles, observed across various bands (centre frequencies and bandwidths), for each beam model. For all the chosen geometric parameters of the pulsar beam, observed in any frequency band, the simulated profiles manifest a relative shift in phase in their observed components, as a result of the intrinsic profile evolution. This relative shift in phase could be interpreted as an additional component to the ISM induced dispersion measure (DM). This additional DM component due to profile evolution is frequency dependent. We discuss the systematics introduced to pulsar data due to this effect.

The Galway Astronomical Stokes Polarimeter (GASP) is a high time resolution, full-Stokes imaging polarimeter which utilises a Fresnel rhomb prism as a beam-splitter and quarter wave retarder. The ability to perform optical photometry and polarimetry at high frame rates enables more detailed studies of a multitude of compact objects including; optical pulsars, magnetars and Active Galactic Nuclei (AGN). We present a brief discussion of the instruments’ current capabilities and hardware design. GASP is currently configured with two Electron Multiplying Charge-Coupled Devices (EMCCDs), which offer sub-millisecond frame readout speeds over a defined region of interest on the sensor. We will report results from an observational campaign at the William Herschel Telescope (WHT) in December, 2015. During this campaign GASP was used to study the Crab Pulsar, V404 Cygni and polarisation standards. As a subset of our analysis we have identified a significant contribution from the Interstellar Medium (ISM) and the interaction of a varying polarised source with the ISM, the subsequent conversion of linearly polarised light to circular, and its dependence on the angle of the source emission electric field orientation. Further to the presentation of results we will discuss future observational work, which is planned for September 2017 and subsequent improvements to increase the temporal resolution of the detectors.