The ANTARES detector is an undersea neutrino telescope in the Mediterranean Sea. The search for point-like neutrino sources is one of the main goals of the ANTARES telescope, requiring a reliable method to evaluate the detector angular resolution and pointing accuracy. This work describes the study of the Sun "shadow" effect with the ANTARES detector. The shadow is the deficit in the atmospheric muon flux in the direction of the Sun caused by the absorption of the primary cosmic rays. This analysis is based on the data collected between 2008 and 2017 by the ANTARES telescope. The observed statistical significance of the Sun shadow detection is $3.9\sigma$, with an estimated angular resolution of $0.45^\circ\pm0.12^\circ$ for downward-going muons. The pointing accuracy is found to be consistent with the expectations and no evidence of systematic pointing shifts is observed.

Context:Spider pulsars are binary systems containing an energetic millisecond pulsar that intensely irradiates a closely orbiting low-mass companion. Modelling their companion's optical light curves is essential to the study of the orbital properties of the binary, including the determination of the pulsar mass, characterising the pulsar wind and the star itself. We aim to generalise the traditional direct heating model of irradiation, whereby energy deposited by the pulsar wind into the stellar envelope is locally re-emitted, by introducing heat redistribution via diffusion and convection within the outer stellar envelope. We approximate the irradiated stellar envelope as a two-dimensional shell. This allows us to propose an effective equation of energy conservation that can be solved at a reduced computational cost. We then implement this model in the \texttt{Icarus} software and use evidence sampling to determine the most likely convection and diffusion laws for the light curve of the redback companion of PSR J2215+5135. Redistribution effects concentrate near the terminator line of pulsar irradiation, and can create apparent hot and cold spots. Among the models tested for PSR J2215+5135, we find that all models with heat redistribution are more likely than symmetric direct heating. The best-fitting redistribution model involves diffusion together with a uniformly rotating envelope. However, we caution that all models still present serious systematic effects, and that prior knowledge from pulsar timing, spectroscopy and distance are key to determine with certainty the most accurate redistribution law. We propose an extension of the direct heating framework that allows for exploring a variety of heat redistribution effects. Future work is necessary to determine the relevant laws from first principles and empirically using complementary observations.

Between May 2016 and September 2018, the intermediate polar (IP) FO Aquarii exhibited two distinct low states and one failed low state. We present optical spectroscopy of FO Aquarii throughout this period, making this the first detailed study of an accretion disc during a low state in any IP. Analysis of these data confirm that the low states are the result of a drop in the mass transfer rate between the secondary star and the magnetic white dwarf primary, and are characterised by a decrease in the system's brightness coupled with a change of the system's accretion structures from an accretion disc-fed geometry to a combination of disc-fed and ballistic stream-fed accretion, and that effects from accretion onto both magnetic poles become detectable. The failed low state only displays a decrease in brightness, with the accretion geometry remaining primarily disc-fed. We also find that the WD appears to be exclusively accretion disc-fed during the high state. There is evidence for an outflow close to the impact region between the ballistic stream and the disc which is detectable in all of the states. Finally, there is marginal evidence for narrow high velocity features in the H$\alpha$ emission line during the low states which may arise due to an outflow from the WD. These features may be evidence of a collimated jet, a long predicted yet elusive feature of cataclysmic variables.

An offline search for a neutrino counterpart to gravitational-wave (GW) events detected during the second observation run (O2) of Advanced-LIGO and Advanced-Virgo performed with ANTARES data is presented. In addition to the search for long tracks induced by $\nu_\mu$ ($\bar{\nu}_{\mu}$) charged current interactions, a search for showering events induced by interactions of neutrinos of any flavour is conducted. The severe spatial and time coincidence provided by the gravitational-wave alert allows regions above the detector horizon to be probed, extending the ANTARES sensitivity over the entire sky. The results of this all-neutrino-flavour and all-sky time dependent analysis are presented. The search for prompt neutrino emission within $\pm$500~s around the time of six GW events yields no neutrino counterparts. Upper limits on the neutrino spectral fluence and constraints on the isotropic radiated energy are set for each GW event analysed.

The discovery that at least some Fast Radio Bursts (FRBs) repeat has ruled out cataclysmic events as the progenitors of these particular bursts. FRB~121102 is the most well-studied repeating FRB but despite extensive monitoring of the source, no underlying pattern in the repetition has previously been identified. Here, we present the results from a radio monitoring campaign of FRB~121102 using the 76-m Lovell telescope. Using the pulses detected in the Lovell data along with pulses from the literature, we report a detection of periodic behaviour of the source over the span of five years of data. We predict that the source is currently `off' and that it should turn `on' for the approximate MJD range $59002-59089$ (2020-06-02 to 2020-08-28). This result, along with the recent detection of periodicity from another repeating FRB, highlights the need for long-term monitoring of repeating FRBs at a high cadence. Using simulations, we show that one needs at least 100 hours of telescope time to follow-up repeating FRBs at a cadence of 0.5--3 days to detect periodicities in the range of 10--150 days. If the period is real, it shows that repeating FRBs can have a large range in their activity periods that might be difficult to reconcile with neutron star precession models.

The KM3NeT research infrastructure is currently under construction at two locations in the Mediterranean Sea. The KM3NeT/ORCA water-Cherenkov neutrino detector off the French coast will instrument several megatons of seawater with photosensors. Its main objective is the determination of the neutrino mass ordering. This work aims at demonstrating the general applicability of deep convolutional neural networks to neutrino telescopes, using simulated datasets for the KM3NeT/ORCA detector as an example. To this end, the networks are employed to achieve reconstruction and classification tasks that constitute an alternative to the analysis pipeline presented for KM3NeT/ORCA in the KM3NeT Letter of Intent. They are used to infer event reconstruction estimates for the energy, the direction, and the interaction point of incident neutrinos. The spatial distribution of Cherenkov light generated by charged particles induced in neutrino interactions is classified as shower- or track-like, and the main background processes associated with the detection of atmospheric neutrinos are recognized. Performance comparisons to machine-learning classification and maximum-likelihood reconstruction algorithms previously developed for KM3NeT/ORCA are provided. It is shown that this application of deep convolutional neural networks to simulated datasets for a large-volume neutrino telescope yields competitive reconstruction results and performance improvements with respect to classical approaches.

Binaries harbouring millisecond pulsars enable a unique path to determine neutron star masses: radio pulsations reveal the motion of the neutron star, while that of the companion can be characterised through studies in the optical range. PSR J1012+5307 is a millisecond pulsar in a 14.5-h orbit with a helium-core white dwarf companion. In this work we present the analysis of an optical spectroscopic campaign, where the companion star absorption features reveal one of the lightest known white dwarfs. We determine a white dwarf radial velocity semi-amplitude of K_2 = 218.9 +- 2.2 km/s, which combined with that of the pulsar derived from the precise radio timing, yields a mass ratio of q=10.44+- 0.11. We also attempt to infer the white dwarf mass from observational constraints using new binary evolution models for extremely low-mass white dwarfs, but find that they cannot reproduce all observed parameters simultaneously. In particular, we cannot reconcile the radius predicted from binary evolution with the measurement from the photometric analysis (R_WD=0.047+-0.003 Rsun). Our limited understanding of extremely low-mass white dwarf evolution, which results from binary interaction, therefore comes as the main factor limiting the precision with which we can measure the mass of the white dwarf in this system. Our conservative white dwarf mass estimate of M_WD = 0.165 +- 0.015 Msun, along with the mass ratio enables us to infer a pulsar mass of M_NS = 1.72 +- 0.16 Msun. This value is clearly above the canonical 1.4 Msun, therefore adding PSR J1012+5307 to the growing list of massive millisecond pulsars.

We present the first measurement of the gravitational quadrupole moment of the companion star of a spider pulsar, namely the black widow PSR J2051$-$0827. To this end we have re-analysed radio timing data using a new model which is able to account for periastron precession caused by tidal and centrifugal deformations of the star as well as by general relativity. The model allows for a time-varying component of the quadrupole moment, thus self-consistently accounting for the ill-understood orbital period variations observed in these systems. Our analysis results in the first detection of orbital precession in a spider system at $\dot{\omega} = -68.6_{-0.5}^{+0.9}$ deg/yr and the most accurate determination of orbital eccentricity for PSR J2051$-$0827 with $e = (4.2 \pm 0.1) \times 10^{-5}$. We show that the variable quadrupole component is about 100 times smaller than the average quadrupole moment $\bar{Q} = -2.2_{-1}^{+0.6} \times 10^{41}$ kg.m$^2$. We discuss how accurate modelling of high precision optical light curves of the companion star will allow its apsidal motion constant to be derived from our results.

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200 Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$ Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma & Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers.

The gSeaGen code is a GENIE-based application developed to efficiently generate high statistics samples of events, induced by neutrino interactions, detectable in a neutrino telescope. The gSeaGen code is able to generate events induced by all neutrino flavours, considering topological differences between track-type and shower-like events. Neutrino interactions are simulated taking into account the density and the composition of the media surrounding the detector. The main features of gSeaGen are presented together with some examples of its application within the KM3NeT project.

We present an optical, X-ray, and $\gamma$-ray study of 1SXPS J042749.2-670434, an eclipsing X-ray binary which has an associated $\gamma$-ray counterpart, 4FGL J0427.8-6704. This association has led to the source being classified as a transitional millisecond pulsar (tMSP) in an accreting state. We analyse 10.5 years of Fermi LAT data, and detect a $\gamma$-ray eclipse at the same phase as optical and X-ray eclipses at the >5$\sigma$ level, a significant improvement on the 2.8$\sigma$level of the previous detection. The confirmation of this eclipse solidifies the association between the X-ray source and the $\gamma$-ray source, strengthening the tMSP classification. However, analysis of several optical data sets and an X-ray observation do not reveal a change in the source's median brightness over long timescales or a bi-modality on short timescales. Instead, the light curve is dominated by flickering which has a correlation time of 2.6 min alongside a potential quasi-periodic oscillation at $\sim$21 min. The mass of the primary and secondary star are constrained to be $M_1=1.43^{+0.33}_{-0.19}$ M$_{\odot}$ and $M_2=0.3^{+0.17}_{-0.12}$ M$_{\odot}$ through modelling of the optical light curve. While this is still consistent with a white dwarf primary, we favour the transitional millisecond pulsar in a low accretion state classification due to the significance of the $\gamma$-ray eclipse detection.

We present the results of the first combined dark matter search targeting the Galactic Centre using the ANTARES and IceCube neutrino telescopes. For dark matter particles with masses from 50 to 1000 GeV, the sensitivities on the self-annihilation cross section set by ANTARES and IceCube are comparable, making this mass range particularly interesting for a joint analysis. Dark matter self-annihilation through the $\tau^+\tau^-$, $\mu^+\mu^-$, $b\bar{b}$ and $W^+W^-$ channels is considered for both the Navarro-Frenk-White and Burkert halo profiles. In the combination of 2,101.6 days of ANTARES data and 1,007 days of IceCube data, no excess over the expected background is observed. Limits on the thermally-averaged dark matter annihilation cross section $\langle\sigma_A\upsilon\rangle$ are set. These limits present an improvement of up to a factor of two in the studied dark matter mass range with respect to the individual limits published by both collaborations. When considering dark matter particles with a mass of 200 GeV annihilating through the $\tau^+\tau^-$ channel, the value obtained for the limit is $7.44 \times 10^{-24} \text{cm}^{3}\text{s}^{-1}$ for the Navarro-Frenk-White halo profile. For the purpose of this joint analysis, the model parameters and the likelihood are unified, providing a benchmark for forthcoming dark matter searches performed by neutrino telescopes.

We present a comparative study of the low-frequency eclipses of spider (compact, irradiating binary) PSRs B1957+20 and J1816+4510. Combining these data with those of three other eclipsing systems we study the frequency dependence of the eclipse duration. PSRs B1957+20 and J1816+4510 have similar orbital properties, but the companions to the pulsars have masses that differ by an order of magnitude. A dedicated campaign to simultaneously observe the pulsed and imaged continuum flux densities throughout the eclipses reveals many similarities between the excess material within the two binaries, irrespective of the companion star properties. The observations show that the pulsar fluxes are removed from the line of sight throughout the main body of the eclipses. For PSR J1816+4510 we present the first direct evidence of an eclipse mechanism that transitions from one that removes the pulsar flux from the line of sight to one that merely smears out pulsations, and claim that this is a consequence of scattering in a tail of material flowing behind the companion. Inferred mass loss rates from the companion stars are found to be $\dot{M}_{\text{C}} \sim 10^{-12}~M_\odot$~yr$^{-1}$ and $\dot{M}_{\text{C}} \sim 2 \times 10^{-13}~M_\odot$~yr$^{-1}$ for PSR B1957+20 and PSR J1816+4510, respectively; seemingly too low to evaporate the stars within Hubble time. Measurements of eclipse durations over a wide range of radio-frequencies show a significant dependence of eclipse duration on frequency for all pulsars, with wider eclipses at lower-frequencies. These results provide a marked improvement in the observational constraints available for theoretical studies of the eclipse mechanisms.

Fast Radio Bursts (FRBs) are bright, extragalactic radio pulses whose origins are still unknown. Until recently, most FRBs have been detected at frequencies greater than 1 GHz with a few exceptions at 800 MHz. The recent discoveries of FRBs at 400 MHz from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope has opened up possibilities for new insights about the progenitors while many other low frequency surveys in the past have failed to find any FRBs. Here, we present results from a FRB survey recently conducted at the Jodrell Bank Observatory at 332 MHz with the 76-m Lovell telescope for a total of 58 days. We did not detect any FRBs in the survey and report a 90$\%$ upper limit of 5500 FRBs per day per sky for a Euclidean Universe above a fluence threshold of 46 Jy ms. We discuss the possibility of absorption as the main cause of non-detections in low frequency (< 800 MHz) searches and invoke different absorption models to explain the same. We find that Induced Compton Scattering alone cannot account for absorption of radio emission and that our simulations favour a combination of Induced Compton Scattering and Free-Free Absorption to explain the non-detections. For a free-free absorption scenario, our constraints on the electron density are consistent with those expected in the post-shock region of the ionized ejecta in Super-Luminous SuperNovae (SLSNe).

A search for point-like and extended sources of cosmic neutrinos using data collected by the ANTARES and IceCube neutrino telescopes is presented. The data set consists of all the track-like and shower-like events pointing in the direction of the Southern Sky included in the nine-year ANTARES point-source analysis, combined with the through-going track-like events used in the seven-year IceCube point-source search. The advantageous field of view of ANTARES and the large size of IceCube are exploited to improve the sensitivity in the Southern Sky by a factor $\sim$2 compared to both individual analyses. In this work, the Southern Sky is scanned for possible excesses of spatial clustering, and the positions of preselected candidate sources are investigated. In addition, special focus is given to the region around the Galactic Centre, whereby a dedicated search at the location of SgrA* is performed, and to the location of the supernova remnant RXJ 1713.7-3946. No significant evidence for cosmic neutrino sources is found and upper limits on the flux from the various searches are presented.

Spider millisecond pulsars are, along with some eclipsing post-common envelope systems and cataclysmic variables, part of an expanding category of compact binaries with low-mass companions for which puzzling timing anomalies have been observed. The most prominent type of irregularities seen in them are orbital period variations, a phenomenon which has been proposed to originate from changes in the gravitational quadrupole moment of the companion star. A physically sound modelling of the timing of these systems is key to understanding their structure and evolution. In this paper we argue that a complete timing model must account for relativistic corrections as well as rotationally and tidally induced quadrupole distortions. We solve for the resulting orbital dynamics using perturbation theory and derive the corresponding timing model in the low eccentricity limit. We find that the expected strong quadrupole deformation of the companion star results in an effective minimum orbital eccentricity. It is accompanied by a fast periastron precession which, if not taken into account, averages out any measurement of the said eccentricity. We show that, with our model, detection of both eccentricity and precession is likely to be made in many if not all spider pulsar systems. Combined with optical light curves, this will allow us to measure the apsidal motion constant, connecting the quadrupole deformation to the internal structure, and thus opening a new window into probing the nature of their exotic stellar interiors. Moreover, more accurate timing may eventually lead spider pulsars to be used for high-precision timing experiments such as pulsar timing arrays.

We analyze 7.3 years of ANTARES high-energy neutrino and Fermi LAT {\gamma}-ray data in search of cosmic neutrino + {\gamma}-ray ({\nu}+{\gamma}) transient sources or source populations. Our analysis has the potential to detect either individual {\nu}+{\gamma} transient sources (durations {\delta}t < 1000~s), if they exhibit sufficient {\gamma}-ray or neutrino multiplicity, or a statistical excess of {\nu}+{\gamma} transients of lower multiplicities. Treating ANTARES track and cascade event types separately, we establish detection thresholds by Monte Carlo scrambling of the neutrino data, and determine our analysis sensitivity by signal injection against scrambled datasets. We find our analysis is sensitive to {\nu}+{\gamma} transient populations responsible for $>$5\% of the observed gamma-coincident neutrinos in the track data at 90\% confidence. Applying our analysis to the unscrambled data reveals no individual {\nu}+{\gamma} events of high significance; two ANTARES track + Fermi {\gamma}-ray events are identified that exceed a once per decade false alarm rate threshold ($p=17\%$). No evidence for subthreshold {\nu}+{\gamma} source populations is found among the track ($p=39\%$) or cascade ($p=60\%$) events. While TXS 0506+056, a blazar and variable (non-transient) Fermi {\gamma}-ray source, has recently been identified as the first source of high-energy neutrinos, the challenges in reconciling observations of the Fermi {\gamma}-ray sky, the IceCube high-energy cosmic neutrinos, and ultra-high energy cosmic rays using only blazars suggest a significant contribution by other source populations. Searches for transient sources of high-energy neutrinos remain interesting, with the potential for neutrino clustering or multimessenger coincidence searches to lead to discovery of the first {\nu}+{\gamma} transients.

The KM3NeT Collaboration runs a multi-site neutrino observatory in the Mediterranean Sea. Water Cherenkov particle detectors, deep in the sea and far off the coasts of France and Italy, are already taking data while incremental construction progresses. Data Acquisition Control software is operating off-shore detectors as well as testing and qualification stations for their components. The software, named Control Unit, is highly modular. It can undergo upgrades and reconfiguration with the acquisition running. Interplay with the central database of the Collaboration is obtained in a way that allows for data taking even if Internet links fail. In order to simplify the management of computing resources in the long term, and to cope with possible hardware failures of one or more computers, the KM3NeT Control Unit software features a custom dynamic resource provisioning and failover technology, which is especially important for ensuring continuity in case of rare transient events in multi-messenger astronomy. The software architecture relies on ubiquitous tools and broadly adopted technologies and has been successfully tested on several operating systems.

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 report on $\sim10$ years of observations of PSR J2051$-$0827, at radio frequencies in the range 110--4032 MHz. We investigate the eclipse phenomena of this black widow pulsar using model fits of increased dispersion and scattering of the pulsed radio emission as it traverses the eclipse medium. These model fits reveal variability in dispersion features on timescales as short as the orbital period, and previously unknown trends on timescales of months--years. No clear patterns are found between the low-frequency eclipse widths, orbital period variations and trends in the intra-binary material density. Using polarisation calibrated observations we present the first available limits on the strength of magnetic fields within the eclipse region of this system; the average line of sight field is constrained to be $10^{-4}$ G $\lesssim B_{||} \lesssim 10^2$ G, while for the case of a field directed near-perpendicular to the line of sight we find $B_{\perp} \lesssim 0.3$ G. Depolarisation of the linearly polarised pulses during the eclipse is detected and attributed to rapid rotation measure fluctuations of $\sigma_{\text{RM}} \gtrsim 100$ rad m$^{-2}$ along, or across, the line of sights averaged over during a sub-integration. The results are considered in the context of eclipse mechanisms, and we find scattering and/or cyclotron absorption provide the most promising explanation, while dispersion smearing is conclusively ruled out. Finally, we estimate the mass loss rate from the companion to be $\dot{M}_{\text{C}} \sim 10^{-12} M_\odot$ yr$^{-1}$, suggesting that the companion will not be fully evaporated on any reasonable timescale.

The KM3NeT research infrastructure being built at the bottom of the Mediterranean Sea will host water-Cherenkov telescopes for the detection of cosmic neutrinos. The neutrino telescopes will consist of large volume three-dimensional grids of optical modules to detect the Cherenkov light from charged particles produced by neutrino-induced interactions. Each optical module houses 31 3-inch photomultiplier tubes, instrumentation for calibration of the photomultiplier signal and positioning of the optical module and all associated electronics boards. By design, the total electrical power consumption of an optical module has been capped at seven watts. This paper presents an overview of the front-end and readout electronics system inside the optical module, which has been designed for a 1~ns synchronization between the clocks of all optical modules in the grid during a life time of at least 20 years.

The Photon Transfer Curve (PTC) of a CCD depicts the variance of uniform images as a function of their average. It is now well established that the variance is not proportional to the average, as Poisson statistics would indicate, but rather flattens out at high flux. This "variance deficit", related to the brighter-fatter effect, feeds correlations between nearby pixels, that increase with flux, and decay with distance. We propose an analytical expression for the PTC shape, and for the dependence of correlations with intensity, and relate both to some more basic quantities related to the electrostatics of the sensor, that are commonly used to correct science images for the brighter-fatter effect. We derive electrostatic constraints from a large set of flat field images acquired with a CCD e2v 250, and eventually question the generally-admitted assumption that boundaries of CCD pixels shift by amounts proportional to the source charges. Our results show that the departure of flat field statistics from Poisson law is entirely compatible with charge redistribution during the drift in the sensor.

KM3NeT will be a network of deep-sea neutrino telescopes in the Mediterranean Sea. The KM3NeT/ARCA detector, to be installed at the Capo Passero site (Italy), is optimised for the detection of high-energy neutrinos of cosmic origin. Thanks to its geographical location on the Northern hemisphere, KM3NeT/ARCA can observe upgoing neutrinos from most of the Galactic Plane, including the Galactic Centre. Given its effective area and excellent pointing resolution, KM3NeT/ARCA will measure or significantly constrain the neutrino flux from potential astrophysical neutrino sources. At the same time, it will test flux predictions based on gamma-ray measurements and the assumption that the gamma-ray flux is of hadronic origin. Assuming this scenario, discovery potentials and sensitivities for a selected list of Galactic sources and to generic point sources with an $E^{-2}$ spectrum are presented. These spectra are assumed to be time independent. The results indicate that an observation with $3\sigma$ significance is possible in about six years of operation for the most intense sources, such as Supernovae Remnants RX\,J1713.7-3946 and Vela Jr. If no signal will be found during this time, the fraction of the gamma-ray flux coming from hadronic processes can be constrained to be below 50\% for these two objects.

In the past years, the IceCube Collaboration has reported in several analyses the observation of astrophysical high-energy neutrino events. Despite a compelling evidence for the first identification of a neutrino source, TXS 0506+056, the origin of the majority of these events is still unknown. In this paper, a possible transient origin of the IceCube astrophysical events is searched for using neutrino events detected by the ANTARES telescope. The arrival time and direction of 6894 track-like and 160 shower-like events detected over 2346 days of livetime are examined to search for coincidences with 54 IceCube high-energy track-like neutrino events, by means of a maximum likelihood method. No significant correlation is observed and upper limits on the one-flavour neutrino fluence from the direction of the IceCube candidates are derived. The non-observation of time and space correlation within the time window of 0.1 days with the two most energetic IceCube events constrains the spectral index of a possible point-like transient neutrino source, to be harder than $-2.3$ and $-2.4$ for each event, respectively.

The shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities from the LOFAR MSSS survey and literature flux densities at various frequencies we have obtained integrated radio spectra for 106 galaxies. The spectra are explained through the use of a three-dimensional model of galaxy radio emission, and radiation transfer dependent on the galaxy viewing angle and absorption processes. Spectra of our galaxies are generally flatter at lower compared to higher frequencies but as there is no tendency for the highly inclined galaxies to have more flattened low-frequency spectra, we argue that the observed flattening is not due to thermal absorption, contradicting the suggestion of Israel & Mahoney (1990). According to our modelled radio maps for M51-like galaxies, the free-free absorption effects can be seen only below 30MHz and in the global spectra just below 20MHz, while in the spectra of starburst galaxies, like M82, the flattening due to absorption is instead visible up to higher frequencies of about 150MHz. Locally, within galactic disks, the absorption effects are distinctly visible in M51-like galaxies as spectral flattening around 100-200MHz in the face-on objects, and as turnovers in the edge-on ones, while in M82-like galaxies there are strong turnovers at frequencies above 700MHz, regardless of viewing angle. Our modelling suggests that the weak spectral flattening observed in the nearby galaxies studied here results principally from synchrotron spectral curvature due to cosmic ray energy losses and propagation effects.

We have observed a large glitch in the Crab pulsar (PSR B0531+21). The glitch occurred around MJD 58064 (2017 November 8) when the pulsar underwent an increase in the rotation rate of $\Delta \nu = 1.530 \times 10^{-5}$ Hz, corresponding to a fractional increase of $\Delta \nu / \nu = 0.516 \times 10^{-6}$ making this event the largest glitch ever observed in this source. Due to our high-cadence and long-dwell time observations of the Crab pulsar we are able to partially resolve a fraction of the total spin-up of the star. This delayed spin-up occurred over a timescale of $\sim$1.7 days and is similar to the behaviour seen in the 1989 and 1996 large Crab pulsar glitches. The spin-down rate also increased at the glitch epoch by $\Delta \dot{\nu} / \dot{\nu} = 7 \times 10^{-3}$. In addition to being the largest such event observed in the Crab, the glitch occurred after the longest period of glitch inactivity since at least 1984 and we discuss a possible relationship between glitch size and waiting time. No changes to the shape of the pulse profile were observed near the glitch epoch at 610 MHz or 1520 MHz, nor did we identify any changes in the X-ray flux from the pulsar. The long-term recovery from the glitch continues to progress as $\dot{\nu}$ slowly rises towards pre-glitch values. In line with other large Crab glitches, we expect there to be a persistent change to $\dot{\nu}$. We continue to monitor the long-term recovery with frequent, high quality observations.

The population of fast radio bursts (FRBs) will continue to diverge into two groups depending on their method of discovery: those which can be localised, and those which cannot. Events potentially less useful for astronomical and cosmological purposes due to limited localisation will accumulate with the advent of new facilities and continued efforts by, e.g., the SUPERB collaboration, which may require afterglows or multi-wavelength counterparts for sub-arcsecond localisation. It is important to exploit these sources to their maximum scientific potential. We perform analysis of FRB dispersion measures (DMs), considering different theoretical FRB progenitors with view to place more rigorous constraints on FRB redshifts, in particular for large statistical samples, via their DMs. We review FRB DM components, and build redshift-scalable probability distributions corresponding to different progenitor scenarios. We combine these components into a framework for obtaining FRB DM probabilities given their redshifts. Taking into account different possibilities for the evolution of progenitors across cosmic time we invert this model, thus deriving redshift constraints. Effects of varying FRB progenitor models are illustrated. While, as expected, host galaxy DM contributions become decreasingly important with increasing redshift, for AGN-like progenitor scenarios they could remain significant out to redshift 3. Constraints are placed on redshifts of catalogued FRBs with various models and increasingly realistic models may be employed as general understanding of FRBs improves. For localised FRBs, we highlight future prospects for disentangling host and intergalactic medium DM components using their respective redshift scaling. We identify a use for large samples of unlocalised FRBs resulting from upcoming flux-limited surveys, such as with CHIME, in mapping out the Milky Way contribution to the DM.

The majority of fast radio bursts (FRBs) are poorly localised, hindering their potential scientific yield as galactic, intergalactic, and cosmological probes. LOFT-e, a digital backend for the U.K.'s e-MERLIN seven-telescope interferometer will provide commensal search and real-time detection of FRBs, taking full advantage of its field of view (FoV), sensitivity, and observation time. Upon burst detection, LOFT-e will store raw data offline, enabling the sub-arcsecond localisation provided by e-MERLIN and expanding the pool of localised FRBs. The high-time resolution backend will additionally introduce pulsar observing capabilities to e-MERLIN.

For 80 days in 2017, the Kepler Space Telescope continuously observed the transitional millisecond pulsar system PSR J1023+0038 in its accreting state. We present analyses of the 59 second cadence data, focusing on investigations of the orbital light curve of the irradiated companion star, and of flaring activity in the neutron star's accretion disc. The underlying orbital modulation from the companion star retains a similar amplitude and asymmetric heating profile as seen in previous photometric observations of the system in its radio pulsar state, suggesting that the heating mechanism has not been affected by the state change. We also find tentative evidence that this asymmetry may vary with time. The light curve also exhibits "flickering" activity, evident as short time-scale flux correlations throughout the observations, and periods of rapid mode-switching activity on time scales shorter than the observation cadence. Finally, the system spent 20% of the observations in a flaring state, with the length of these flares varying from 2 minutes up to several hours. The flaring behaviour is consistent with a self-organised criticality mechanism, most likely related to the build up and release of mass at the inner edge of the accretion disc.

We present LOFAR high-band data over the frequency range 115-189 MHz for the X-ray binary SS 433, obtained in an observing campaign from 2013 February - 2014 May. Our results include a deep, wide-field map, allowing a detailed view of the surrounding supernova remnant W 50 at low radio frequencies, as well as a light curve for SS 433 determined from shorter monitoring runs. The complex morphology of W 50 is in excellent agreement with previously published higher-frequency maps; we find additional evidence for a spectral turnover in the eastern wing, potentially due to foreground free-free absorption. Furthermore, SS 433 is tentatively variable at 150 MHz, with both a debiased modulation index of 11 per cent and a $\chi^2$ probability of a flat light curve of $8.2 \times 10^{-3}$. By comparing the LOFAR flux densities with contemporaneous observations carried out at 4800 MHz with the RATAN-600 telescope, we suggest that an observed $\sim$0.5-1 Jy rise in the 150-MHz flux density may correspond to sustained flaring activity over a period of approximately six months at 4800 MHz. However, the increase is too large to be explained with a standard synchrotron bubble model. We also detect a wealth of structure along the nearby Galactic plane, including the most complete detection to date of the radio shell of the candidate supernova remnant G 38.7-1.4. This further demonstrates the potential of supernova remnant studies with the current generation of low-frequency radio telescopes.

We have observed and analysed the eclipses of the black widow pulsar J1810+1744 at low radio frequencies. Using LOw-Frequency ARray (LOFAR) and Westerbork Synthesis Radio Telescope observations between 2011--2015 we have measured variations in flux density, dispersion measure and scattering around eclipses. High-time-resolution, simultaneous beamformed and interferometric imaging LOFAR observations show concurrent disappearance of pulsations and total flux from the source during the eclipses, with a $3\sigma$ upper limit of 36 mJy ($<10\%$ of the pulsar's averaged out-of-eclipse flux density). The dispersion measure variations are highly asymmetric, suggesting a tail of material swept back due to orbital motion. The egress deviations are variable on timescales shorter than the 3.6 hr orbital period and are indicative of a clumpy medium. Additional pulse broadening detected during egress is typically $<20\%$ of the pulsar's spin period, showing no evidence of scattering the pulses beyond detectability in the beamformed data. The eclipses, lasting $\sim13\%$ of the orbit at 149 MHz, are shown to be frequency-dependent with total duration scaling as $\propto\nu^{-0.41\pm0.03}$. The results are discussed in the context of the physical parameters of the system, and an examination of eclipse mechanisms reveals cyclotron-synchrotron absorption as the most likely primary cause, although non-linear scattering mechanisms cannot be quantitatively ruled out. The inferred mass loss rate is a similar order-of-magnitude to the mean rate required to fully evaporate the companion in a Hubble time.

Black widows and redbacks are binary systems consisting of a millisecond pulsar in a close binary with a companion having matter driven off of its surface by the pulsar wind. X-rays due to an intra-binary shock have been observed from many of these systems, as well as orbital variations in the optical emission from the companion due to heating and tidal distortion. We have been systematically studying these systems in radio, optical and X-rays. Here we will present an overview of X-ray and optical studies of these systems, including new XMM-Newton and NuStar data obtained from several of them, along with new optical photometry.

Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most $\sim 0.02\;M_{\odot}$. We investigate DNS masses, spins and velocities, and in particular correlations between spin period, orbital period and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes which may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short $\gamma$-ray burst.

We report the discovery of PSR J0952$-$0607, a 707-Hz binary millisecond pulsar which is now the fastest-spinning neutron star known in the Galactic field (i.e., outside of a globular cluster). PSR J0952$-$0607 was found using LOFAR at a central observing frequency of 135 MHz, well below the 300 MHz to 3 GHz frequencies typically used in pulsar searches. The discovery is part of an ongoing LOFAR survey targeting unassociated Fermi Large Area Telescope $\gamma$-ray sources. PSR J0952$-$0607 is in a 6.42-hr orbit around a very low-mass companion ($M_\mathrm{c}\gtrsim0.02$ M$_\odot$) and we identify a strongly variable optical source, modulated at the orbital period of the pulsar, as the binary companion. The light curve of the companion varies by 1.6 mag from $r^\prime=22.2$ at maximum to $r^\prime>23.8$, indicating that it is irradiated by the pulsar wind. Swift observations place a 3-$\sigma$ upper limit on the $0.3-10$ keV X-ray luminosity of $L_X < 1.1 \times 10^{31}$ erg s$^{-1}$ (using the 0.97 kpc distance inferred from the dispersion measure). Though no eclipses of the radio pulsar are observed, the properties of the system classify it as a black widow binary. The radio pulsed spectrum of PSR J0952$-$0607, as determined through flux density measurements at 150 and 350 MHz, is extremely steep with $\alpha\sim-3$ (where $S \propto \nu^{\alpha}$). We discuss the growing evidence that the fastest-spinning radio pulsars have exceptionally steep radio spectra, as well as the prospects for finding more sources like PSR J0952$-$0607.

Linares et al. (2016) obtained quasi-simultaneous g', r' and i-band light curves and an absorption line radial velocity curve of the secondary star in the redback system 3FGL J0212.1+5320. The light curves showed two maxima and minima primarily due to the secondary star's ellipsoidal modulation, but with unequal maxima and minima. We fit these light curves and radial velocities with our X-ray binary model including either a dark solar-type star spot or a hot spot due to off-centre heating from an intrabinary shock, to account for the unequal maxima. Both models give a radial velocity semi-amplitude and rotational broadening that agree with the observations. The observed secondary star's effective temperature is best matched with the value obtained using the hot spot model, which gives a neutron star and secondary star mass of $M_{\rm 1}$=1.85$^{+0.32}_{-0.26}$ $M_{\odot}$and $M_{\rm 2}$=0.50$^{+0.22}_{-0.19}$ $M_{\odot}$, respectively.

We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) at 140 MHz. We find a scaling relation $P_{\rm cav} \propto L_{148}^{\beta}$, with a logarithmic slope of $\beta = 0.51 \pm 0.14$, which is in good agreement with previous results based on data at 327 MHz. The large scatter present in this correlation confirms the conclusion reached at higher frequencies that the total radio luminosity at a single frequency is a poor predictor of the total jet power. We show that including measurements at 148 MHz alone is insufficient to reliably compute the bolometric radio luminosity and reduce the scatter in the correlation. For a subset of four well-resolved sources, we examine the detected extended structures at low frequencies and compare with the morphology known from higher frequency images and Chandra X-ray maps. In Perseus we discuss details in the structures of the radio mini-halo, while in the 2A 0335+096 cluster we observe new diffuse emission associated with multiple X-ray cavities and likely originating from past activity. For A2199 and MS 0735.6+7421, we confirm that the observed low-frequency radio lobes are confined to the extents known from higher frequencies. This new low-frequency analysis highlights the fact that existing cavity power to radio luminosity relations are based on a relatively narrow range of AGN outburst ages. We discuss how the correlation could be extended using low frequency data from the LOFAR Two-metre Sky Survey (LoTSS) in combination with future, complementary deeper X-ray observations.

We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispersion of $215.86 \pm 6.34$ km/s. From archival radio observations we see that this galaxy hosts a compact flat-spectrum radio source, and we conclude that it is the active galactic nucleus (AGN) responsible for generating the radio lobes. The radio luminosity distribution of the jets, and the broad-line classification of the host AGN, indicate this GRG is orientated well out of the plane of the sky, making its physical size one of the largest known for any GRG. Analysis of the infrared data suggests that the host is a lenticular type galaxy with a large stellar mass ($\log~\mathrm{M}/\mathrm{M}_\odot = 11.56 \pm 0.12$), and a moderate star formation rate ($1.2 \pm 0.3~\mathrm{M}_\odot/\mathrm{year}$). Spatially smoothing the SDSS images shows the system around UGC 9555 to be significantly disturbed, with a prominent extension to the south-east. Overall, the evidence suggests this host galaxy has undergone one or more recent moderate merger events and is also experiencing tidal interactions with surrounding galaxies, which have caused the star formation and provided the supply of gas to trigger and fuel the Mpc-scale radio lobes.

We present first results from a LOFAR census of non-recycled pulsars. The census includes almost all such pulsars known (194 sources) at declinations ${\rm Dec}> 8^\circ$ and Galactic latitudes $|{\rm Gb}| > 3^\circ$, regardless of their expected flux densities and scattering times. Each pulsar was observed for $\geq 20$ minutes in the contiguous frequency range of 110--188 MHz. Full-Stokes data were recorded. We present the dispersion measures, flux densities, and calibrated total intensity profiles for the 158 pulsars detected in the sample. The median uncertainty in census dispersion measures ($1.5 \times 10^{-3}$ pc cm$^{-3}$) is ten times smaller, on average, than in the ATNF pulsar catalogue. We combined census flux densities with those in the literature and fitted the resulting broadband spectra with single or broken power-law functions. For 48 census pulsars such fits are being published for the first time. Typically, the choice between single and broken power-laws, as well as the location of the spectral break, were highly influenced by the spectral coverage of the available flux density measurements. In particular, the inclusion of measurements below 100 MHz appears essential for investigating the low-frequency turnover in the spectra for most of the census pulsars. For several pulsars, we compared the spectral indices from different works and found the typical spread of values to be within 0.5--1.5, suggesting a prevailing underestimation of spectral index errors in the literature. The census observations yielded some unexpected individual source results, as we describe in the paper. Lastly, we will provide this unique sample of wide-band, low-frequency pulse profiles via the European Pulsar Network Database.

Brown dwarfs and giant planets orbiting close to a host star are subjected to significant irradiation that can modify the properties of their atmospheres. In order to test the atmospheric models that are used to describe these systems, it is necessary to obtain accurate observational estimates of their physical properties (masses, radii, temperatures, albedos). Interacting compact binary systems provide a natural laboratory for studying strongly irradiated sub-stellar objects. As the mass-losing secondary in these systems makes a critical, but poorly understood transition from the stellar to the sub-stellar regime, it is also strongly irradiated by the compact accretor. In fact, the internal and external energy fluxes are both expected to be comparable in these objects, providing access to an unexplored irradiation regime. However, the atmospheric properties of such donors have so far remained largely unknown. Here, we report the direct spectroscopic detection and characterisation of an irradiated sub-stellar donor in an accreting white dwarf binary system. Our near-infrared observations allow us to determine a model-independent mass estimate for the donor of $M_2=0.055\pm0.008M_{\odot}$ and an average spectral type of ${\rm L1}\pm{\rm1}$, supporting both theoretical predictions and model-dependent observational constraints. Our time-resolved data also allow us to estimate the average irradiation-induced temperature difference between the day and night sides on the sub-stellar donor, $\Delta{\rm T}\simeq57$~K, and the maximum difference between the hottest and coolest parts of its surface, of $\Delta{\rm T}_{max}\simeq200$~K. The observations are well described by a simple geometric reprocessing model with a bolometric (Bond) albedo of $A_B<0.54$ at 2-$\sigma$ confidence level, consistent with high reprocessing efficiency, but poor lateral heat redistribution in the donor's atmosphere.

The eclipses of certain types of binary millisecond pulsars (i.e. `black widows' and `redbacks') are often studied using high-time-resolution, `beamformed' radio observations. However, they may also be detected in images generated from interferometric data. As part of a larger imaging project to characterize the variable and transient sky at radio frequencies <200 MHz, we have blindly detected the redback system PSR J2215+5135 as a variable source of interest with the Low-Frequency Array (LOFAR). Using observations with cadences of 2 weeks - 6 months, we find preliminary evidence that the eclipse duration is frequency dependent ($\propto \nu^{-0.4}$), such that the pulsar is eclipsed for longer at lower frequencies, in broad agreement with beamformed studies of other similar sources. Furthermore, the detection of the eclipses in imaging data suggests an eclipsing medium that absorbs the pulsed emission, rather than scattering it. Our study is also a demonstration of the prospects of finding pulsars in wide-field imaging surveys with the current generation of low-frequency radio telescopes.

We report on the results of a search for radio transients between 115 and 190\,MHz with the LOw-Frequency ARray (LOFAR). Four fields have been monitored with cadences between 15 minutes and several months. A total of 151 images were obtained, giving a total survey area of 2275 deg$^2$. We analysed our data using standard LOFAR tools and searched for radio transients using the LOFAR Transients Pipeline (TraP). No credible radio transient candidate has been detected; however, we are able to set upper limits on the surface density of radio transient sources at low radio frequencies. We also show that low-frequency radio surveys are more sensitive to steep-spectrum coherent transient sources than GHz radio surveys. We used two new statistical methods to determine the upper limits on the transient surface density. One is free of assumptions on the flux distribution of the sources, while the other assumes a power-law distribution in flux and sets more stringent constraints on the transient surface density. Both of these methods provide better constraints than the approach used in previous works. The best value for the upper limit we can set for the transient surface density, using the method assuming a power-law flux distribution, is 1.3$\cdot$10$^{-3}$ deg$^{-2}$ for transients brighter than 0.3 Jy with a time-scale of 15 min, at a frequency of 150 MHz. We also calculated for the first time upper limits for the transient surface density for transients of different time-scales. We find that the results can differ by orders of magnitude from previously reported, simplified estimates.

The relativistic double neutron star binary PSR J0737-3039 shows clear evidence of orbital phase-dependent wind-companion interaction, both in radio and X-rays. In this paper we present the results of timing analysis of PSR J0737-3039 performed during 2006 and 2011 XMM-Newton Large Programs that collected ~20,000 X-ray counts from the system. We detected pulsations from PSR J0737-3039A (PSR A) through the most accurate timing measurement obtained by XMM-Newton so far, the spin period error being of 2x10^-13 s. PSR A's pulse profile in X-rays is very stable despite significant relativistic spin precession that occurred within the time span of observations. This yields a constraint on the misalignment between the spin axis and the orbital momentum axis Delta_A ~6.6^{+1.3}_{-5.4} deg, consistent with estimates based on radio data. We confirmed pulsed emission from PSR J0737-3039B (PSR B) in X-rays even after its disappearance in radio. The unusual phenomenology of PSR B's X-ray emission includes orbital pulsed flux and profile variations as well as a loss of pulsar phase coherence on time scales of years. We hypothesize that this is due to the interaction of PSR A's wind with PSR B's magnetosphere and orbital-dependent penetration of the wind plasma onto PSR B closed field lines. Finally, the analysis of the full XMM-Newton dataset provided evidences of orbital flux variability (~7%) for the first time, involving a bow-shock scenario between PSR A's wind and PSR B's magnetosphere.

LS I +61 303 is a gamma-ray binary that exhibits an outburst at GHz frequencies each orbital cycle of $\approx$ 26.5 d and a superorbital modulation with a period of $\approx$ 4.6 yr. We have performed a detailed study of the low-frequency radio emission of LS I +61 303 by analysing all the archival GMRT data at 150, 235 and 610 MHz, and conducting regular LOFAR observations within the Radio Sky Monitor (RSM) at 150 MHz. We have detected the source for the first time at 150 MHz, which is also the first detection of a gamma-ray binary at such a low frequency. We have obtained the light-curves of the source at 150, 235 and 610 MHz, all of them showing orbital modulation. The light-curves at 235 and 610 MHz also show the existence of superorbital variability. A comparison with contemporaneous 15-GHz data shows remarkable differences with these light-curves. At 15 GHz we see clear outbursts, whereas at low frequencies we see variability with wide maxima. The light-curve at 235 MHz seems to be anticorrelated with the one at 610 MHz, implying a shift of $\sim$ 0.5 orbital phases in the maxima. We model the shifts between the maxima at different frequencies as due to changes in the physical parameters of the emitting region assuming either free-free absorption or synchrotron self-absorption, obtaining expansion velocities for this region close to the stellar wind velocity with both mechanisms.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOFAR in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths, and flux densities and compare these with higher observing frequencies. The flux-calibrated, multifrequency LOFAR pulse profiles are publicly available via the EPN Database of Pulsar Profiles. We also present average values of dispersion measures (DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered profiles, and the rest are only weakly detected. A qualitative comparison of LOFAR profiles with those at higher radio frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As previously demonstrated this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterising and correcting these systematic effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect gravitational waves.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation.

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

Fast Radio Bursts (FRBs), are millisecond radio signals that exhibit dispersion larger than what the Galactic electron density can account for. We have conducted a 1446 hour survey for Fast Radio Bursts (FRBs) at 145~MHz, covering a total of 4193 sq. deg on the sky. We used the UK station of the LOFAR radio telescope -- the Rawlings Array -- , accompanied for a majority of the time by the LOFAR station at Nan\c{c}ay, observing the same fields at the same frequency. Our real-time search backend, ARTEMIS, utilizes graphics processing units to search for pulses with dispersion measures up to 320 cm$^{-3}$ pc. Previous derived FRB rates from surveys around 1.4~GHz, and favoured FRB interpretations, motivated this survey, despite all previous detections occurring at higher dispersion measures. We detected no new FRBs above a signal-to-noise threshold of 10, leading to the most stringent upper limit yet on the FRB event rate at these frequencies: 29 sky$^{-1}$ day$^{-1}$ for 5~ms-duration pulses above 62~Jy. The non-detection could be due to scatter-broadening, limitations on the volume and time searched, or the shape of FRB flux density spectra. Assuming the latter and that FRBs are standard candles, the non-detection is compatible with the published FRB sky rate, if their spectra follow a power law with frequency ($\propto \nu^{\alpha}$), with $\alpha\gtrsim+0.1$, demonstrating a marked difference from pulsar spectra. Our results suggest that surveys at higher frequencies, including the low frequency component of the Square Kilometre Array, will have better chances to detect, estimate rates and understand the origin and properties of FRBs.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

Current and future astronomical survey facilities provide a remarkably rich opportunity for transient astronomy, combining unprecedented fields of view with high sensitivity and the ability to access previously unexplored wavelength regimes. This is particularly true of LOFAR, a recently-commissioned, low-frequency radio interferometer, based in the Netherlands and with stations across Europe. The identification of and response to transients is one of LOFAR's key science goals. However, the large data volumes which LOFAR produces, combined with the scientific requirement for rapid response, make automation essential. To support this, we have developed the LOFAR Transients Pipeline, or TraP. The TraP ingests multi-frequency image data from LOFAR or other instruments and searches it for transients and variables, providing automatic alerts of significant detections and populating a lightcurve database for further analysis by astronomers. Here, we discuss the scientific goals of the TraP and how it has been designed to meet them. We describe its implementation, including both the algorithms adopted to maximize performance as well as the development methodology used to ensure it is robust and reliable, particularly in the presence of artefacts typical of radio astronomy imaging. Finally, we report on a series of tests of the pipeline carried out using simulated LOFAR observations with a known population of transients.

Impulsive radio bursts that are detectable across cosmological distances constitute extremely powerful probes of the ionized Inter-Galactic Medium (IGM), intergalactic magnetic fields, and the properties of space-time itself. Their dispersion measures (DMs) will enable us to detect the "missing" baryons in the low-redshift Universe and make the first measurements of the mean galaxy halo profile, a key parameter in models of galaxy formation and feedback. Impulsive bursts can be used as cosmic rulers at redshifts exceeding 2, and constrain the dark energy equation-of-state parameter, $w(z)$ at redshifts beyond those readily accessible by Type Ia SNe. Both of these goals are realisable with a sample of $\sim 10^4$ fast radio bursts (FRBs) whose positions are localized to within one arcsecond, sufficient to obtain host galaxy redshifts via optical follow-up. It is also hypothesised that gravitational wave events may emit coherent emission at frequencies probed by SKA1-LOW, and the localization of such events at cosmological distances would enable their use as cosmological standard sirens. To perform this science, such bursts must be localized to their specific host galaxies so that their redshifts may be obtained and compared against their dispersion measures, rotation measures, and scattering properties. The SKA can achieve this with a design that has a wide field-of-view, a substantial fraction of its collecting area in a compact configuration (80\% within a 3\,km radius), and a capacity to attach high-time-resolution instrumentation to its signal path.

We present the highest-quality polarisation profiles to date of 16 non-recycled pulsars and four millisecond pulsars, observed below 200 MHz with the LOFAR high-band antennas. Based on the observed profiles, we perform an initial investigation of expected observational effects resulting from the propagation of polarised emission in the pulsar magnetosphere and the interstellar medium. The predictions of magnetospheric birefringence in pulsars have been tested using spectra of the pulse width and fractional polarisation from multifrequency data. The derived spectra offer only partial support for the expected effects of birefringence on the polarisation properties, with only about half of our sample being consistent with the model's predictions. It is noted that for some pulsars these measurements are contaminated by the effects of interstellar scattering. For a number of pulsars in our sample, we have observed significant variations in the amount of Faraday rotation as a function of pulse phase, which is possibly an artefact of scattering. These variations are typically two orders of magnitude smaller than that observed at 1400 MHz by Noutsos et al. (2009), for a different sample of southern pulsars. In this paper we present a possible explanation for the difference in magnitude of this effect between the two frequencies, based on scattering. Finally, we have estimated the magnetospheric emission heights of low-frequency radiation from four pulsars, based on the phase lags between the flux-density and the PA profiles, and the theoretical framework of Blaskiewicz, Cordes & Wasserman (1991). These estimates yielded heights of a few hundred km; at least for PSR B1133+16, this is consistent with emission heights derived based on radius-to-frequency mapping, but is up to a few times larger than the recent upper limit based on pulsar timing.

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.

The SKA will discover tens of thousands of pulsars and provide unprecedented data quality on these, as well as the currently known population, due to its unrivalled sensitivity. Here, we outline the state of the art of our understanding of magnetospheric radio emission from pulsars and how we will use the SKA to solve the open problems in pulsar magnetospheric physics.

Since their discovery in the late 1960's the population of known neutron stars (NSs) has grown to ~2500. The last five decades of observations have yielded many surprises and demonstrated that the observational properties of NSs are remarkably diverse. The surveys that will be performed with SKA (the Square Kilometre Array) will produce a further tenfold increase in the number of Galactic NSs known. Moreover, the SKA's broad spectral coverage, sub-arraying and multi-beaming capabilities will allow us to characterise these sources with unprecedented efficiency, in turn enabling a giant leap in the understanding of their properties. Here we review the NS population and outline our strategies for studying each of the growing number of diverse classes that are populating the "NS zoo". Some of the main scientific questions that will be addressed by the much larger statistical samples and vastly improved timing efficiency provided by SKA include: (i) the spin period and spin-down rate distributions (and thus magnetic fields) at birth, and the associated information about the SNe wherein they are formed; (ii) the radio pulsar-magnetar connection; (iii) the link between normal radio pulsars, intermittent pulsars and rotating radio transients; (iv) the slowest possible spin period for a radio pulsar (revealing the conditions at the pulsar death-line); (v) proper motions of pulsars (revealing SN kick physics); (vi) the mass distribution of NSs (vii) the fastest possible spin period for a recycled pulsar (constraining magnetosphere-accretion disc interactions, gravitational wave radiation and the equation-of-state); (viii) the origin of high eccentricity millisecond pulsars (MSPs); (ix) the formation channels for recently identified triple systems; and finally (x) how isolated MSPs are formed. We expect that the SKA will break new ground unveiling exotic systems that will challenge... [abridged]

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

We report the optical identification of the companion to the {\it Fermi} black widow millisecond pulsar PSR J1544+4937. We find a highly variable source on Keck LRIS images at the nominal pulsar position, with 2 magnitude variations over orbital period in the B, g, R, and I bands. The nearly achromatic light curves are difficult to explain with a simply irradiated hemisphere model, and suggest that the optical emission is dominated by a nearly isothermal hot patch on the surface of the companion facing the pulsar. We roughly constrain the distance to PSR J1544+4937 to be between 2 and 5 kpc. A more reliable distance measurement is needed in order to constrain the composition of the companion.

We present a study of the quiescent optical counterpart of the Accreting Millisecond X-ray Pulsar XTE J1814-338, carrying out multiband (BVR) orbital phase-resolved photometry using the ESO VLT/FORS2. The optical light curves are consistent with a sinusoidal variability modulated with the orbital period, showing evidence for a strongly irradiated companion star, in agreement with previous findings. The observed colours cannot be accounted for by the companion star alone, suggesting the presence of an accretion disc during quiescence. The system is fainter in all analysed bands compared to previous observations. The R band light curve displays a possible phase offset with respect to the B and V band. Through a combined fit of the multi-band light curves we derive constraints on the companion star and disc fluxes, on the system distance and on the companion star mass. The irradiation luminosity required to account for the observed day-side temperature of the companion star is consistent with the spin-down luminosity of a millisecond radio pulsar. The flux decrease and spectral evolution of the quiescent optical emission observed comparing our data with previous observations, collected over 5 years, cannot be well explained with the contribution of an irradiated companion star and an accretion disc alone. The progressive flux decrease as the system gets bluer could be due to a continuum component evolving towards a lower, bluer spectrum. While most of the continuum component is likely due to the disc, we do not expect it to become bluer in quiescence. Hence we hypothesize that an additional component, such as synchrotron emission from a jet was contributing significantly in the earlier data obtained during quiescence and then progressively fading or moving its break frequency toward longer wavelengths.

The double pulsar system PSR J0737-3039A/B is a double neutron star binary, with a 2.4-hour orbital period, which has allowed measurement of relativistic orbital perturbations to high precision. The low mass of the second-formed neutron star, as well as the low system eccentricity and proper motion, point to a different evolutionary scenario compared to most other known double neutron star systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to be a near-orthogonal rotator, with an average separation between its spin and magnetic axes of 90 +/- 11 +/- 5 deg. Furthermore, we find a mean 95% upper limit on the misalignment between its spin and orbital angular momentum axes of 3.2 deg, assuming that the observed emission comes from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

In the last few years, over 43 millisecond radio pulsars have been discovered by targeted searches of unidentified gamma-ray sources found by the Fermi Gamma-Ray Space Telescope. A large fraction of these millisecond pulsars are in compact binaries with low-mass companions. These systems often show eclipses of the pulsar signal and are commonly known as black widows and redbacks because the pulsar is gradually destroying its companion. In this paper, we report on the optical discovery of four strongly irradiated millisecond pulsar companions. All four sources show modulations of their color and luminosity at the known orbital periods from radio timing. Light curve modelling of our exploratory data shows that the equilibrium temperature reached on the companion's dayside with respect to their nightside is consistent with about 10-30% of the available spin-down energy from the pulsar being reprocessed to increase the companion's dayside temperature. This value compares well with the range observed in other irradiated pulsar binaries and offers insights about the energetics of the pulsar wind and the production of gamma-ray emission. In addition, this provides a simple way of estimating the brightness of irradiated pulsar companions given the pulsar spin-down luminosity. Our analysis also suggests that two of the four new irradiated pulsar companions are only partially filling their Roche lobe. Some of these sources are relatively bright and represent good targets for spectroscopic follow-up. These measurements could enable, among other things, mass determination of the neutron stars in these systems.

We report on time-resolved optical imaging of the X-ray binary SAX J1808.4-3658 during its quiescent state and 2008 outburst. The binary, containing an accretion-powered millisecond pulsar, has a large sinusoidal-like modulation in its quiescent optical emission. We employ a Markov chain Monte Carlo technique to fit our multi-band light curve data in quiescence with an irradiated star model, and derive a tight constraint of $50^{+6}_{-5}$ deg on the inclination angle $i$ of the binary system. The pulsar and its companion are constrained to have masses of $0.97^{+0.31}_{-0.22} M_{\sun}$ and $0.04^{+0.02}_{-0.01} M_{\sun}$ (both 1$\sigma$ ranges), respectively. The dependence of these results on the measurements of the companion's projected radial velocity is discussed. We also find that the accretion disk had nearly constant optical fluxes over a $\sim$500 day period in the quiescent state our data covered, but started brightening 1.5 months before the 2008 outburst. Variations in modulation during the outburst were detected in our four observations made 7-12 days after the start of the outburst, and a sinusoidal-like modulation with 0.2 mag amplitude changed to have a smaller amplitude of 0.1 mag. The modulation variations are discussed. We estimate the albedo of the companion during its quiescence and the outburst, which was approximately 0 and 0.8 (for isotropic emission), respectively. This large difference probably provides additional evidence that the neutron star in the binary turns on as a radio pulsar in quiescence.

The double pulsar PSR J0737-3039A/B displays short, 30 s eclipses that arise around conjunction when the radio waves emitted by pulsar A are absorbed as they propagate through the magnetosphere of its companion pulsar B. These eclipses offer a unique opportunity to probe directly the magnetospheric structure and the plasma properties of pulsar B. We have performed a comprehensive analysis of the eclipse phenomenology using multi-frequency radio observations obtained with the Green Bank Telescope. We have characterized the periodic flux modulations previously discovered at 820 MHz by McLaughlin et al., and investigated the radio frequency dependence of the duration and depth of the eclipses. Based on their weak radio frequency evolution, we conclude that the plasma in pulsar B's magnetosphere requires a large multiplicity factor (~ 10^5). We also found that, as expected, flux modulations are present at all radio frequencies in which eclipses can be detected. Their complex behavior is consistent with the confinement of the absorbing plasma in the dipolar magnetic field of pulsar B as suggested by Lyutikov & Thompson and such a geometric connection explains that the observed periodicity is harmonically related to pulsar B's spin frequency. We observe that the eclipses require a sharp transition region beyond which the plasma density drops off abruptly. Such a region defines a plasmasphere which would be well inside the magnetospheric boundary of an undisturbed pulsar. It is also two times smaller than the expected standoff radius calculated using the balance of the wind pressure from pulsar A and the nominally estimated magnetic pressure of pulsar B.

We present an analysis and interpretation of the Kepler binary system KOI 1224. This is the fourth binary found with Kepler that consists of a thermally bloated, hot white dwarf in a close orbit with a more or less normal star of spectral class A or F. As we show, KOI 1224 contains a white dwarf with Teff = 14400 +/- 1100 K, mass = 0.20 +/- 0.02 Msun, and radius = 0.103 +/- 0.004 Rsun, and an F-star companion of mass = 1.59 +/- 0.07 Msun that is somewhat beyond its terminal-age main sequence. The orbital period is quite short at 2.69802 days. The ingredients that are used in the analysis are the Kepler binary light curve, including the detection of the Doppler boosting effect; the NUV and FUV fluxes from the Galex images of this object; an estimate of the spectral type of the F-star companion; and evolutionary models of the companion designed to match its effective temperature and mean density. The light curve is modelled with a new code named Icarus which we describe in detail. Its features include the full treatment of orbital phase-resolved spectroscopy, Doppler boosting, irradiation effects and transits/eclipses, which are particularly suited to irradiated eclipsing binaries. We interpret the KOI 1224 system in terms of its likely evolutionary history. We infer that this type of system, containing a bloated hot white dwarf, is the direct descendant of an Algol-type binary. In spite of this basic understanding of the origin of KOI 1224, we discuss a number of problems associated with producing this type of system with this short of an short orbital period.

The most massive neutron stars constrain the behavior of ultra-dense matter, with larger masses possible only for increasingly stiff equations of state. Here, we present evidence that the black widow pulsar, PSR B1957+20, has a high mass. We took spectra of its strongly irradiated companion and found an observed radial-velocity amplitude of K_obs=324+/-3 km/s. Correcting this for the fact that, due to the irradiation, the center of light lies inward relative to the center of mass, we infer a true radial-velocity amplitude of K_2=353+/-4 km/s and a mass ratio q=M_PSR/M_2=69.2+/-0.8. Combined with the inclination i=65+/-2 deg inferred from models of the lightcurve, our best-fit pulsar mass is M_PSR=2.40+/-0.12 M_sun. We discuss possible systematic uncertainties, in particular in the lightcurve modeling. Taking an upper limit of i<85 deg based on the absence of radio eclipses at high frequency, combined with a conservative lower-limit to the motion of the center of mass, K_2>343 km/s (q>67.3), we infer a lower limit to the pulsar mass of M_PSR>1.66 M_sun.

We present the evolution of the radio emission from the 2.8-s pulsar of the double pulsar system PSR J0737-3039A/B. We provide an update on the Burgay et al. (2005) analysis by describing the changes in the pulse profile and flux density over five years of observations, culminating in the B pulsar's radio disappearance in 2008 March. Over this time, the flux density decreases by 0.177 mJy/yr at the brightest orbital phases and the pulse profile evolves from a single to a double peak, with a separation rate of 2.6 deg/yr. The pulse profile changes are most likely caused by relativistic spin precession, but can not be easily explained with a circular hollow-cone beam as in the model of Clifton & Weisberg (2008). Relativistic spin precession, coupled with an elliptical beam, can model the pulse profile evolution well. This particular beam shape predicts geometrical parameters for the two bright orbital phases which are consistent and similar to those derived by Breton et al. (2008). However, the observed decrease in flux over time and B's eventual disappearance cannot be easily explained by the model and may be due to the changing influence of A on B.

Among the initial results from Kepler were two striking lightcurves, for KOI 74 and KOI 81, in which the relative depths of the primary and secondary eclipses showed that the more compact, less luminous object was hotter than its stellar host. That result became particularly intriguing because a substellar mass had been derived for the secondary in KOI 74, which would make the high temperature challenging to explain; in KOI 81, the mass range for the companion was also reported to be consistent with a substellar object. We re-analyze the Kepler data and demonstrate that both companions are likely to be white dwarfs. We also find that the photometric data for KOI 74 show a modulation in brightness as the more luminous star orbits, due to Doppler boosting. The magnitude of the effect is sufficiently large that we can use it to infer a radial velocity amplitude accurate to 1 km/s. As far as we are aware, this is the first time a radial-velocity curve has been measured photometrically. Combining our velocity amplitude with the inclination and primary mass derived from the eclipses and primary spectral type, we infer a secondary mass of 0.22+/-0.03 Msun. We use our estimates to consider the likely evolutionary paths and mass-transfer episodes of these binary systems.

This thesis focuses on the study of binary radio pulsars, their evolution and some specific use of their properties to investigate fundamental physics such as general relativity and other gravitational theories. The work that we present here is organized in three main parts. First, we report on the study of PSR J1744-3922, a binary pulsar presenting a peculiar `flickering' flux behavior as well as spin and orbital properties that do not correspond to the expectations of standard evolution scenarios. Second, we conducted an in-depth analysis of the eclipses in the relativistic double pulsar system PSR J0737-3039A/B. From our modeling of the eclipses, we precisely determined the geometry of pulsar B in space and used this information to study the temporal behavior of the eclipses, which revealed that pulsar B precesses around the angular momentum of the system in a way that is consistent with the prediction of general relativity. Third, we searched for the signature of latitudinal aberration in the pulse profile of pulsar A in the double pulsar system. The non-detection of this effect allows us to put an upper limit on its amplitude, which constrains the geometry of pulsar A with respect to our line of sight as well as its emission geometry. (Abridged)

The double pulsar PSR J0737-3039A/B consists of two neutron stars in a highly relativistic orbit that displays a roughly 30-second eclipse when pulsar A passes behind pulsar B. Describing this eclipse of pulsar A as due to absorption occurring in the magnetosphere of pulsar B, we successfully use a simple geometric model to characterize the observed changing eclipse morphology and to measure the relativistic precession of pulsar B's spin axis around the total orbital angular momentum. This provides a test of general relativity and alternative theories of gravity in the strong-field regime. Our measured relativistic spin precession rate of 4.77 (+0.66,-0.65) degrees per year (68% confidence level) is consistent with that predicted by general relativity within an uncertainty of 13%.

The double pulsar system PSR J0737-3039A/B is a highly relativistic double neutron star (DNS) binary, with a 2.4-hour orbital period. The low mass of the second-formed NS, as well the low system eccentricity and proper motion, point to a different evolutionary scenario compared to other known DNS systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to have a low misalignment between its spin and orbital angular momentum axes, with a 68.3% upper limit of 6.1 degrees, assuming emission from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

PSR J1744-3922 is a binary pulsar exhibiting highly variable pulsed radio emission. We report on a statistical multi-frequency study of the pulsed radio flux variability which suggests that this phenomenon is extrinsic to the pulsar and possibly tied to the companion, although not strongly correlated with orbital phase. The pulsar has an unusual combination of characteristics compared to typical recycled pulsars: a long spin period (172 ms); a relatively high magnetic field strength (1.7x10^10 G); a very circular, compact orbit of 4.6 hours; and a low-mass companion (0.08 Msun). These spin and orbital properties are likely inconsistent with standard evolutionary models. We find similarities between the properties of the PSR J1744-3922 system and those of several other known binary pulsar systems, motivating the identification of a new class of binary pulsars. We suggest that this new class could result from either: a standard accretion scenario of a magnetar or a high-magnetic field pulsar; common envelope evolution with a low-mass star and a neutron star, similar to what is expected for ultra-compact X-ray binaries; or, accretion induced collapse of a white dwarf. We also report the detection of a possible K'=19.30(15) infrared counterpart at the position of the pulsar, which is relatively bright if the companion is a helium white dwarf at the nominal distance, and discuss its implications for the pulsar's companion and evolutionary history.