Binary black hole (BH) central engine description for the unique blazar OJ 287 predicted that the next secondary BH impact-induced bremsstrahlung flare should peak on 2019 July 31. This prediction was based on detailed general relativistic modeling of the secondary BH trajectory around the primary BH and its accretion disk. The expected flare was termed the Eddington flare to commemorate the centennial celebrations of now-famous solar eclipse observations to test general relativity by Sir Arthur Eddington. We analyze the multi-epoch Spitzer observations of the expected flare between 2019 July 31 and 2019 September 6, as well as baseline observations during 2019 February-March. Observed Spitzer flux density variations during the predicted outburst time display a strong similarity with the observed optical pericenter flare from OJ 287 during 2007 September. The predicted flare appears comparable to the 2007 flare after subtracting the expected higher base-level Spitzer flux densities at 3.55 and 4.49 $\mu$m compared to the optical R-band. Comparing the 2019 and 2007 outburst lightcurves and the previously calculated predictions, we find that the Eddington flare arrived within 4 hours of the predicted time. Our Spitzer observations are well consistent with the presence of a nano-Hertz gravitational wave emitting spinning massive binary BH that inspirals along a general relativistic eccentric orbit in OJ 287. These multi-epoch Spitzer observations provide a parametric constraint on the celebrated BH no-hair theorem.

Polarization plays an important role in various time-domain astrophysics to understand the magnetic fields, geometry, and environments of spatially unresolved variable sources. In this paper we present the results of laboratory and on-sky testing of a novel dual-beam, dual-camera optical imaging polarimeter (MOPTOP) exploiting high sensitivity, low-noise CMOS technology and designed to monitor variable and transient sources with low systematic errors and high sensitivity. We present a data reduction algorithm that corrects for sensitivity variations between the cameras on a source-by-source basis. Using our data reduction algorithm, we show that our dual-beam, dual-camera technique delivers the benefits of low and stable instrumental polarization ($<0.05$\% for lab data and $<0.25$\% for on sky data) and high throughput while avoiding the additional sky brightness and image overlap problems associated with dual-beam, single-camera polarimeters.

M87 is one of the closest (z=0.00436) extragalactic sources emitting at very-high-energies (VHE, E > 100 GeV). The aim of this work is to locate the region of the VHE gamma-ray emission and to describe the observed broadband spectral energy distribution (SED) during the low VHE gamma-ray state. The data from M87 collected between 2012 and 2015 as part of a MAGIC monitoring programme are analysed and combined with multi-wavelength data from Fermi-LAT, Chandra, HST, EVN, VLBA and the Liverpool Telescope. The averaged VHE gamma-ray spectrum can be fitted from 100GeV to 10TeV with a simple power law with a photon index of (-2.41 $\pm$ 0.07), while the integral flux above 300GeV is $(1.44 \pm 0.13) \times 10^{-12} cm^{-2} s^{-1}$. During the campaign between 2012 and 2015, M87 is generally found in a low emission state at all observed wavelengths. The VHE gamma-ray flux from the present 2012-2015 M87 campaign is consistent with a constant flux with some hint of variability ($\sim3\sigma$) on a daily timescale in 2013. The low-state gamma-ray emission likely originates from the same region as the flare-state emission. Given the broadband SED, both a leptonic synchrotron self Compton and a hybrid photo-hadronic model reproduce the available data well, even if the latter is preferred. We note, however, that the energy stored in the magnetic field in the leptonic scenario is very low suggesting a matter dominated emission region.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.

In the binary black hole model of OJ 287 the secondary black hole orbits a much more massive primary, and impacts on the primary accretion disk at predictable times. We update the parameters of the disk, the viscosity $\alpha$ and the mass accretion rate $\dot m$. We find $\alpha=0.26 \pm 0.1$ and $\dot m = 0.08 \pm 0.04$ in Eddington units. The former value is consistent with Coroniti (1981) and the latter with Marscher and Jorstad (2011). Predictions are made for the 2019 July 30 superflare in OJ 287. We expect that it will take place simultaneously at the Spitzer infrared channels as well as in the optical and that therefore the timing of the flare in optical can be accurately determined from Spitzer observations. We also discuss in detail the light curve of the 2015 flare and find that the radiating volume has regions where bremsstrahlung dominates as well as regions that radiate primarily in synchrotron radiation. The former region produces the unpolarised first flare while the latter region gives rise to a highly polarized second flare.

The bright blazar OJ 287 is the best-known candidate for hosting a nanohertz gravitational wave (GW) emitting supermassive binary black hole (SMBBH) in the present observable universe. The binary black hole (BBH) central engine model, proposed by Lehto and Valtonen in 1996, was influenced by the two distinct periodicities inferred from the optical light curve of OJ 287. The current improved model employs an accurate general relativistic description to track the trajectory of the secondary black hole (BH) which is crucial to predict the inherent impact flares of OJ 287. The successful observations of three predicted impact flares open up the possibility of using this BBH system to test general relativity in a hitherto unexplored strong field regime. Additionally, we briefly describe an on-going effort to interpret observations of OJ 287 in a Bayesian framework.

The origin of the luminosity of superluminous supernovae (SLSNe) is an unresolved mystery, and a number of very different physical scenarios (including energy injection from magnetars, collision with a dense circumstellar medium and pair instability-induced explosions) have been invoked. The application of polarimetry to normal SNe has been shown to probe the three-dimensional structure of exploding stars, providing clues to the nature of the explosion mechanism. We report imaging linear polarimetry observations of the Type I SLSN 2017egm, in the galaxy NGC 3191, conducted with the Liverpool Telescope and the RINGO3 instrument. Observations were acquired at four epochs, spanning 4 - 19 days after light-curve maximum, however, polarization was not detected at a level of $>3\sigma$. At +7 and +15 days, and in the average over all epochs, we find a possible polarization signal, detected at a significance of $\approx 2\sigma$ in the "blue" channel. This signal is seen, primarily, in the Stokes $q$ parameter, with a corresponding polarization angle consistent with the orientation of the spiral arm in proximity to the position of SN 2017egm. We interpret this as indicating that any polarization, if present, originates from dust in the host galaxy rather than being intrinsic to the SN itself. Despite its apparent peculiarities, compared to other Type I SLSNe, the polarization characteristics of SN 2017egm are consistent with the previously reported low polarization of other SLSNe of this variety.

Blazar OJ287 exhibits large thermal flares at least twice every 12 years. The times of these flares have been predicted successfully using the model of a quasi-Keplerian eccentric black hole binary where the secondary impacts the accretion disk of the primary, creating the thermal flares. New measurements of the historical light curve have been combined with the observations of the 2015 November/December flare to identify the impact record since year 1886, and to constrain the orbit of the binary. The orbit solution shows that the binary period, now 12.062 year, is decreasing at the rate of 36 days per century. This corresponds to an energy loss to gravitational waves that is 6.5 +- 4 % less than the rate predicted by the standard quadrupolar gravitational wave (GW) emission. We show that the difference is due to higher order gravitational radiation reaction terms that include the dominant order tail contributions.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process.

We present optical polarimetry in the period May-August 2017 of the enigmatic "dipping" star KIC 8462852. During that period three ~1% photometric dips were reported by other observers. We measured the average absolute polarization of the source, and find no excess or unusual polarization compared to a nearby comparison star. We place tight upper limits on any change in the degree of polarization of the source between epochs in- and out-of-dip of <0.1% (8500-Ang.) and <0.2% (7050-Ang. and 5300-Ang.). How our limits are interpreted depends on the specific model being considered. If the whole stellar disk were covered by material with an optical depth of ~0.01 then the fractional polarisation introduced by this material must be less than 10-20%. While our non-detection does not constrain the comet scenario, it predicts that even modest amounts of dust that have properties similar to Solar System comets may be detectable. We note that the sensitivity of our method scales with the depth of the dip. Should a future ~20% photometric dip be observed (as was previously detected by Kepler) our method would constrain any induced polarization associated with any occulting material to 0.5-1.0%.

We present a catalog of early-time (~10^2-10^4s) photometry and polarimetry of all Gamma-Ray Burst (GRB) optical afterglows observed with RINGO2 imaging polarimeter on the Liverpool Telescope. For the 19 optical afterglows observed, the following 9 were bright enough to perform photometry and attempt polarimetry: GRB 100805A, GRB 101112A, GRB 110205A, GRB 110726A, GRB 120119A, GRB 120308A, GRB 120311A, GRB 120326A and GRB 120327A. We present multi-wavelength light curves for these 9 GRBs, together with estimates of their optical polarization degrees and/or limits. We carry out a thorough investigation of detection probabilities, instrumental properties and systematics. Using two independent methods, we confirm previous reports of significant polarization in GRB 110205A and 120308A, and report new detection of $P=6^{+3}_{-2}\%$ in GRB101112A. We discuss the results for the sample in the context of the reverse and forward shock afterglow scenario, and show that GRBs with detectable optical polarization at early time have clearly identifiable signatures of reverse-shock emission in their optical light curves. This supports the idea that GRB ejecta contain large-scale magnetic fields and highlights the importance of rapid-response polarimetry.

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.

Around the time of its perihelion passage the observability of 67P/Churyumov-Gerasimenko from Earth was limited to very short windows each morning from any given site, due to the low solar elongation of the comet. The peak in the comet's activity was therefore difficult to observe with conventionally scheduled telescopes, but was possible where service/queue scheduled mode was possible, and with robotic telescopes. We describe the robotic observations that allowed us to measure the total activity of the comet around perihelion, via photometry (dust) and spectroscopy (gas), and compare these results with the measurements at this time by Rosetta's instruments. The peak of activity occurred approximately two weeks after perihelion. The total brightness (dust) largely followed the predictions from Snodgrass et al. 2013, with no significant change in total activity levels from previous apparitions. The CN gas production rate matched previous orbits near perihelion, but appeared to be relatively low later in the year.

We present ~2000 polarimetric and ~3000 photometric observations of 15 gamma-ray bright blazars over a period of 936 days (11/10/2008 - 26/10/2012) using data from the Tuorla blazar monitoring program (KVA DIPOL) and Liverpool Telescope (LT) RINGO2 polarimeters (supplemented with data from SkyCamZ (LT) and Fermi-LAT gamma-ray data). In 11 out of 15 sources we identify a total of 19 electric vector position angle (EVPA) rotations and 95 flaring episodes. We group the sources into subclasses based on their broadband spectral characteristics and compare their observed optical and gamma-ray properties. We find that (1) the optical magnitude and gamma-ray flux are positively correlated, (2) EVPA rotations can occur in any blazar subclass, 4 sources show rotations that go in one direction and immediately rotate back, (3) we see no difference in the gamma-ray flaring rates in the sample; flares can occur during and outside of rotations with no preference for this behaviour, (4) the average degree of polarisation (DoP), optical magnitude and gamma-ray flux are lower during an EVPA rotation compared with during non-rotation and the distribution of the DoP during EVPA rotations is not drawn from the same parent sample as the distribution outside rotations, (5) the number of observed flaring events and optical polarisation rotations are correlated, however we find no strong evidence for a temporal association between individual flares and rotations and (6) the maximum observed DoP increases from ~10% to ~30% to ~40% for subclasses with synchrotron peaks at high, intermediate and low frequencies respectively.

We describe the design, construction and commissioning of LOTUS; a simple, low-cost long-slit spectrograph for the Liverpool Telescope. The design is optimized for near-UV and visible wavelengths and uses all transmitting optics. It exploits the instrument focal plane field curvature to partially correct axial chromatic aberration. A stepped slit provides narrow (2.5x95 arcsec) and wide (5x25 arcsec) options that are optimized for spectral resolution and flux calibration respectively. On sky testing shows a wavelength range of 3200-6300 Angstroms with a peak system throughput (including detector quantum efficiency) of 15 per cent and wavelength dependant spectral resolution of R=225-430. By repeated observations of the symbiotic emission line star AG Peg we demonstrate the wavelength stability of the system is less than 2 Angstroms rms and is limited by the positioning of the object in the slit. The spectrograph is now in routine operation monitoring the activity of comet 67P/Churyumov-Gerasimenko during its current post-perihelion apparition.

OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.

IO:I is a new instrument that has recently been commissioned for the Liverpool Telescope, extending current imaging capabilities beyond the optical and into the near infrared. Cost has been minimised by use of a previously decommissioned instrument's cryostat as the base for a prototype and retrofitting it with Teledyne's 1.7$\mu m$ cutoff Hawaii-2RG HgCdTe detector, SIDECAR ASIC controller and JADE2 interface card. In this paper, the mechanical, electronic and cryogenic aspects of the cryostat retrofitting process will be reviewed together with a description of the software/hardware setup. This is followed by a discussion of the results derived from characterisation tests, including measurements of read noise, conversion gain, full well depth and linearity. The paper closes with a brief overview of the autonomous data reduction process and the presentation of results from photometric testing conducted on on-sky, pipeline processed data.