We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties of the pulsars, which showed that they have, on average, larger characteristic ages. We present the pulse profiles of the pulsars across the three observing bands, where PSR J1643$+$1338 showed profile evolution that appears not to be well-described by the radius-to-frequency-mapping model. Furthermore, we modelled the spectra of the pulsars across the same observing bands, using a simple power law, and found an average spectral index of $-1.9 \pm 0.5$. Amongst the pulsars studied here, PSR J1657$+$3304 showed large flux density variations of a factor of 10 over 300 days, as well as mode changing and nulling on timescales of a few minutes. We modelled the rotational and orbital properties of PSR J1658$+$3630, which has a spin period of 33 ms in a binary orbit of 3.0 days with a companion of minimum mass of 0.87$M_{\odot}$, likely a Carbon-Oxygen or Oxygen-Neon-Magnesium type white dwarf. PSR J1658$+$3630 has a dispersion measure of 3.0 pc cm$^{-3}$, making it possibly one of the closest binary pulsars known.

The Low-Frequency Array radio telescope discovered the $707$ Hz binary millisecond pulsar (MSP) J0952$-$0607 in a targeted radio pulsation search of an unidentified $\textit{Fermi}$ gamma-ray source. This source shows a weak energy flux of $F_\gamma = 2.6 \times 10^{-12}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$ in the energy range between $100\,\text{MeV}$ and $100\,\text{GeV}$. Here we report the detection of pulsed gamma-ray emission from PSR$\,$J0952$-$0607 in a very sensitive gamma-ray pulsation search. The pulsar's rotational, binary, and astrometric properties are measured over seven years of $\textit{Fermi}$-Large Area Telescope data. For this we take into account the uncertainty on the shape of the gamma-ray pulse profile. We present an updated radio-timing solution now spanning more than two years and show results from optical modeling of the black-widow-type companion based on new multi-band photometric data taken with HiPERCAM on the Gran Telescopio Canarias on La Palma and ULTRACAM on the New Technology Telescope at ESO La Silla. PSR$\,$J0952$-$0607 is now the fastest-spinning pulsar for which the intrinsic spin-down rate has been reliably constrained ($\dot{P}_\text{int} \lesssim 4.6 \times 10^{-21}\,\text{s}\,\text{s}^{-1}$). The inferred surface magnetic field strength of $B_\text{surf} \lesssim 8.2 \times 10^{7}\,\text{G}$ is among the ten lowest of all known pulsars. This discovery is another example of an extremely fast spinning black-widow pulsar hiding within an unidentified $\textit{Fermi} gamma-ray source. In the future such systems might help to pin down the maximum spin frequency and the minimum surface magnetic field strength of MSPs.

In this paper, we present the results of timing observations of PSRs J1949+3106 and J1950+2414, two binary millisecond pulsars discovered in data from the Arecibo ALFA pulsar survey (PALFA). The timing parameters include precise measurements of the proper motions of both pulsars, which show that PSR J1949+3106 has a transversal motion very similar to that of an object in the local standard of rest. The timing also includes measurements of the Shapiro delay and the rate of advance of periastron for both systems. Assuming general relativity, these allow estimates of the masses of the components of the two systems; for PSR J1949+3106, the pulsar mass is $M_p \, = \, 1.34^{+0.17}_{-0.15} \, M_{\odot}$ and the companion mass $M_c \, = \, 0.81^{+0.06}_{-0.05}\, M_{\odot}$; for PSR J1950+2414 $M_p \, = \, 1.496 \, \pm \, 0.023\, M_{\odot}$ and $M_c \, = \, 0.280^{+0.005}_{-0.004}\, M_{\odot}$ (all values 68.3 % confidence limits). We use these masses and proper motions to investigate the evolutionary history of both systems: PSR J1949+3106 is likely the product of a low-kick supernova; PSR J1950+2414 is a member of a new class of eccentric millisecond pulsar binaries with an unknown formation mechanism. We discuss the proposed hypotheses for the formations of these systems in light of our new mass measurements.

Dispersive delays due to the Solar wind introduce excess noise in high-precision pulsar timing experiments, and must be removed in order to achieve the accuracy needed to detect, e.g., low-frequency gravitational waves. In current pulsar timing experiments, this delay is usually removed by approximating the electron density distribution in the Solar wind either as spherically symmetric, or with a two-phase model that describes the contributions from both high- and low-speed phases of the Solar wind. However, no dataset has previously been available to test the performance and limitations of these models over extended timescales and with sufficient sensitivity. Here we present the results of such a test with an optimal dataset of observations of pulsar J0034-0534, taken with the German stations of LOFAR. We conclude that the spherical approximation performs systematically better than the two-phase model at almost all angular distances, with a residual root-mean-square (rms) given by the two-phase model being up to 28% larger than the result obtained with the spherical approximation. Nevertheless, the spherical approximation remains insufficiently accurate in modelling the Solar-wind delay (especially within 20 degrees of angular distance from the Sun), as it leaves timing residuals with rms values that reach the equivalent of 0.3 microseconds at 1400 MHz. This is because a spherical model ignores the large daily variations in electron density observed in the Solar wind. In the short term, broadband observations or simultaneous observations at low frequencies are the most promising way forward to correct for Solar-wind induced delay variations.

Context. High-precision pulsar-timing experiments are affected by temporal variations of the Dispersion Measure (DM), which are related to spatial variations in the interstellar electron content. Correcting for DM variations relies on the cold-plasma dispersion law which states that the dispersive delay varies with the squared inverse of the observing frequency. This may however give incorrect measurements if the probed electron content (and therefore the DM) varies with observing frequency, as is predicted theoretically. Aims. We study small-scale density variations in the ionised interstellar medium. These structures may lead to frequency-dependent DMs in pulsar signals and could inhibit the use of lower-frequency pulsar observations to correct time-variable interstellar dispersion in higher-frequency pulsar-timing data. Methods. We used high-cadence, low-frequency observations with three stations from the German LOng-Wavelength (GLOW) consortium, which are part of the LOw Frequency ARray (LOFAR). Specifically, 3.5 years of weekly observations of PSR J2219+4754 are presented. Results. We present the first detection of frequency-dependent DMs towards any interstellar object and a precise multi-year time-series of the time- and frequency-dependence of the measured DMs. The observed DM variability is significant and may be caused by extreme scattering events. Potential causes for frequency-dependent DMs are quantified and evaluated. Conclusions. We conclude that frequency-dependence of DMs has been reliably detected and is caused by small-scale (up to 10s of AUs) but steep density variations in the interstellar electron content. We find that long-term trends in DM variability equally affect DMs measured at both ends of our frequency band and hence the negative impact on long-term high-precision timing projects is expected to be limited.

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.