We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year project which aims to observe (a) more than 1000 pulsars to obtain high-fidelity pulse profiles, (b) some 500 of these pulsars over multiple epochs, (c) long sequences of single-pulse trains from several hundred pulsars. The scientific outcomes from the programme will include determination of pulsar geometries, the location of the radio emission within the pulsar magnetosphere, the connection between the magnetosphere and the crust and core of the star, tighter constraints on the nature of the radio emission itself as well as interstellar medium studies. First results presented here include updated dispersion measures, 26 pulsars with Faraday rotation measures derived for the first time and a description of interesting emission phenomena observed thus far.

We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.

Here, we report on the detection and verification of Fast Radio Burst FRB 180301, which occurred on UTC 2018 March 1 during the Breakthrough Listen observations with the Parkes telescope. Full-polarization voltage data of the detection were captured--a first for non-repeating FRBs--allowing for coherent de-dispersion and additional verification tests. The coherently de-dispersed dynamic spectrum of FRB 180301 shows complex, polarized frequency structure over a small fractional bandwidth. As FRB 180301 was detected close to the geosynchronous satellite band during a time of known 1-2 GHz satellite transmissions, we consider whether the burst was due to radio interference emitted or reflected from an orbiting object. Based on the preponderance of our verification tests, we find that FRB 180301 is likely of astrophysical origin, but caution that anthropogenic sources cannot conclusively be ruled out.

$\textit{Herschel}$ observations of far infrared N$^+$ emission lines have demonstrated that dense plasma, with $n_e\sim30\,{\rm cm^{-3}}$, is ubiquitous in the inner Galactic plane. By combining the information from $\textit{Herschel}$ with other tracers of ionised gas, we build a picture of this dense plasma. We adopt a collisional ionisation model, so the analysis is not tied to a specific energisation mechanism. We find that the dense plasma is concentrated in a disk that is $130\;{\rm pc}$ thick, and makes a significant contribution to radio pulsar dispersion measures in the inner Galactic plane. The strength of the far infrared N$^+$ emission requires high temperatures in the plasma, with $T \simeq 19{,}000\,{\rm K}$ indicated both by the ratio of N$^+$ to C$^+$, and by the ratio of N$^+$ to microwave bremsstrahlung in the inner Galactic plane. This parallels the situation at high Galactic latitudes, where strong optical emission is observed from N$^+$ (and S$^+$), relative to both H$\alpha$ and microwave bremsstrahlung, and suggests a common origin. If so, the same gas provides a natural explanation for the extreme radio-wave scattering phenomena that are sometimes observed in pulsars and quasars. We therefore propose a new picture of the warm ionised medium as seen in emission, in which the plasma is dense, hot, and localised in numerous structures of size $\sim 10^2\,{\rm AU}$ that are clustered around stars.

Fast Radio Bursts (FRBs) are a class of short-duration transients at radio wavelengths with inferred astrophysical origin. The prototypical FRB is a broadband signal that occurs over the extent of the receiver frequency range, is narrow in time, and is highly dispersed, following a $\nu^{-2}$ relation. However, some FRBs appear band-limited, and show apparent scintillation, complex frequency-dependent structure, or multi-component pulse shapes. While there is sufficient evidence that FRBs are indeed astrophysical, their one-off nature necessitates extra scrutiny when reporting a detection as bona fide and not a false positive. Currently, there is no formal validation framework for FRBs, rather a set of community practices. In this article, we discuss potential sources of false positives, and suggest a framework in which FRB-like events can be evaluated as real or otherwise. We present examples of false-positive events in data from the Arecibo, LOFAR, and Nanshan telescopes, which while FRB-like, are found to be due to instrumental variations, noise, and radio-frequency interference. Differentiating these false-positive detections from astrophysical events requires knowledge and tests beyond thresholded single-pulse detection. We discuss post-detection analyses, verification tests, and datasets which should be provided when reporting an FRB detection.

We have observed a complex and continuous change in the integrated pulse profile of PSR B2217+47, manifested as additional components trailing the main peak. These transient components are detected over 6 years at $150$ MHz using the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous Lovell observations at $1.5$ GHz. We argue that propagation effects in the ionized interstellar medium (IISM) are the most likely cause. The putative structures in the IISM causing the profile variation are roughly half-way between the pulsar and the Earth and have transverse radii $R \sim 30$ AU. We consider different models for the structures. Under the assumption of spherical symmetry, their implied average electron density is $\overline{n}_e \sim 100$ cm$^{-3}$. Since PSR B2217+47 is more than an order of magnitude brighter than the average pulsar population visible to LOFAR, similar profile variations would not have been identified in most pulsars, suggesting that subtle profile variations in low-frequency profiles might be more common than we have observed to date. Systematic studies of these variations at low frequencies can provide a new tool to investigate the proprieties of the IISM and the limits to the precision of pulsar timing.

We measure the effects of interstellar scattering on average pulse profiles from 13 radio pulsars with simple pulse shapes. We use data from the LOFAR High Band Antennas, at frequencies between 110 and 190~MHz. We apply a forward fitting technique, and simultaneously determine the intrinsic pulse shape, assuming single Gaussian component profiles. We find that the constant $\tau$, associated with scattering by a single thin screen, has a power-law dependence on frequency $\tau \propto \nu^{-\alpha}$, with indices ranging from $\alpha = 1.50$ to $4.0$, despite simplest theoretical models predicting $\alpha = 4.0$ or $4.4$. Modelling the screen as an isotropic or extremely anisotropic scatterer, we find anisotropic scattering fits lead to larger power-law indices, often in better agreement with theoretically expected values. We compare the scattering models based on the inferred, frequency dependent parameters of the intrinsic pulse, and the resulting correction to the dispersion measure (DM). We highlight the cases in which fits of extreme anisotropic scattering are appealing, while stressing that the data do not strictly favour either model for any of the 13 pulsars. The pulsars show anomalous scattering properties that are consistent with finite scattering screens and/or anisotropy, but these data alone do not provide the means for an unambiguous characterization of the screens. We revisit the empirical $\tau$ versus DM relation and consider how our results support a frequency dependence of $\alpha$. Very long baseline interferometry, and observations of the scattering and scintillation properties of these sources at higher frequencies, will provide further evidence.

Observations of pulsars across the radio spectrum are revealing a dependence of the characteristic scattering time ($\tau$) on frequency, which is more complex than the simple power law with a theoretically predicted power law index. In this paper we investigate these effects using simulated pulsar data at frequencies below 300 MHz. We investigate different scattering mechanisms, namely isotropic and anisotropic scattering, by thin screens along the line of sight, and the particular frequency dependent impact on pulsar profiles and scattering time scales of each. We also consider how the screen shape, location and offset along the line of sight lead to specific observable effects. We evaluate how well forward fitting techniques perform in determining $\tau$. We investigate the systematic errors in $\tau$ associated with the use of an incorrect fitting method and with the determination of an off-pulse baseline. Our simulations provide examples of average pulse profiles at various frequencies. Using these we compute spectra of $\tau$ and mean flux for different scattering setups. We identify setups that lead to deviations from the simple theoretical picture. This work provides a framework for interpretation of upcoming low frequency data, both in terms of modelling the interstellar medium and understanding intrinsic emission properties of pulsars.

This paper gives a complete characterization of resonant orbits in a Kerr spacetime. A resonant orbit is defined as a geodesic for which the longitudinal and radial orbital frequencies are commensurate. Our analysis is based on expressing the resonance condition in its most symmetric form using Carlson's integrals. We provide a number of concise formulae for the dependence of resonances on the system parameters. Resonant effects may be observable during the in-spiral of a compact object into a super-massive black hole. When the slowly evolving orbital frequencies pass through a series of low-order resonances, rapid changes in the orbital parameters could produce measurable phase shifts in the emitted gravitational radiation (GW). Resonant orbits may also capture dust leading to electromagnetic emission. The KAM theorem indicates that, low order resonant orbits demarcate the regions where the onset of chaos could occur around a perturbed black-hole. We find that the 1/2 and 2/3 resonances occur at ~4 and 5.4 Schwarzschild radii (Rs) from the event horizon. For compact object in-spirals around super-massive black holes, this region lies within the sensitivity band of space-based GW detectors. For Sgr A*, length scales of ~41 and 55 microarcseconds and timescales of 50 and 79 min respectively should be associated with resonant effects, if Sgr A* is non-spinning. Spin decreases these values by up to ~32% and ~28%. These length-scales are potentially resolvable with VLBI measurements. We find that all low-order resonances are localized to the strong field region r < 50 Rs. This fact guarantees the validity of using approximations based on averaging to model the frequency evolution of a test object in region 50 Rs <r <1000 Rs. The systematic determination of the multipole moments of the central object by observing the orbit of a pulsar, free of chaotic effects, is thus possible.

We compute the length and timescales associated with resonant orbits in the Kerr Metric for all orbital and spin parameters. Resonance induced effects are potentially observable when the Event Horizon telescope resolves the inner structure of Sgr A*, space-based gravitational wave detectors record phase-shifts in the waveform during the resonant passage of a neutron star as it spirals in to the black hole and in the frequencies of quasi periodic oscillations for accreting black holes. The onset of geodesic chaos for non-Kerr spacetimes should occur at the resonance locations quantified here.

The formation of massive stellar clusters in turbulent molecular clouds is investigated. We include artificial star formation and energy feedback of newly born stars. The obtained systems are not likely to survive. Case studies to determine conditions necessary for forming bound clusters will be done in the future.

The effect of gas ejection on the structure and binding energy of newly formed stellar clusters is investigated. The star formation efficiency (SFE), necessary for forming a gravitationally bound stellar cluster, is determined. Two sets of numerical N-body simulations are presented: As a first simplified approach we treat the residual gas as an external potential. The gas expulsion is approximated by reducing the gas mass to zero on a given timescale, which is treated as a free parameter. In a second set of simulations we use smoothed particle hydrodynamics (SPH) to follow the dynamics of the outflowing residual gas self-consistently. We investigate cases where gas outflow is induced by an outwards propagating shock front and where the whole gas cloud is heated homogeneously, leading to ejection. If the stars are in virial equilibrium with the gaseous environment initially, bound clusters only form in regions where the local SFE is larger than 50% or where the gas expulsion timescale is long compared to the dynamical timescale. A small initial velocity dispersion of the stars leads to a compaction of the cluster during the expulsion phase and reduces the SFE needed to form bound clusters to less than 10%.