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.

The magnetograms of sixteen Active Regions, observed during June - November 2015, were obtained from the Joint Science Operations Center and used to determine the fractal dimensions of various magnetic field groups. The fields were divided into strong (positive and negative) and weak (positive and negative) groupings. The area-perimeter method was used to determine the fractal dimensions for the umbral and penumbral regions of the magnetograms. The fractal dimensions were found to be 1.79 \pm 0.49 and 1.96 \pm 0.29 for strong and weak magnetic fields respectively. When compared to the umbral and penumbral fractal dimensions determined by (Rajkumar, Haque and Hrudey, 2017) using white light images, an inverse relationship was found, with the umbra showing lower fractal dimensions than the penumbra for the magnetic field groups compared to the intensity images. This directly implies that there is greater complexity in the penumbral magnetic field groupings than the umbral region. This has implications for models that constrain how the magnetic field groups are contained in the cooler umbral regions.

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.

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 images of sunspots in sixteen active regions taken at the UCCI Observatory in Grand Cayman during June - November 2015, were used to determine their fractal dimensions using the perimeter-area method for the umbral and the penumbral region. Scale free fractal dimensions of 2.09 $\pm$ 0.42 and 1.72 $\pm$ 0.4 were found respectively. This value was higher than the value determined by Chumak and Chumak (1996) who used a similar method but for the penumbral region only for their sample set. There is a positive correlation r = 0.58 between umbral and the penumbral fractal dimensions for the specific sunspots. Furthermore, a time series analysis was done similarly on eight images of AR 12403, from 21st August 2015 - 28th August 2015 taken from the Debrecen Photoheliographic Data (DPD). The correlation, r = 0.623 between the umbral and penumbral fractal dimensions in the time series indicating that the complexity in morphology indicated by the fractal dimension between the umbra and penumbra followed each other in time as well.

We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies.

A non-minimally coupled scalar field can have, in principle, a negative effective Planck mass squared which depends on the scalar field. Surprisingly, an isotropic and homogeneous cosmological universe with a non-minimally coupled scalar field is perfectly smooth as the rolling scalar field causes the effective Planck mass to change sign and pass through zero. However, we show that any small deviations from homogeneity diverge as the effective Planck mass vanishes, with catastrophic consequences for the cosmology. The physical origin of the divergence is due to the presence of non-zero scalar anisotropic stress from the non-minimally coupled scalar field. Thus, while the homogeneous and isotropic cosmology appears surprisingly sensible when the effective Planck mass vanishes, inhomogeneities tell a different story.

(Abridged) OJ287 is a BL Lac object that has shown double-peaked bursts at regular intervals of ~12 yr during the last ~40 yr. We analyse optical photopolarimetric monitoring data from 2005-2009, during which the latest double-peaked outburst occurred. The aim of this study is twofold: firstly, we aim to analyse variability patterns and statistical properties of the optical polarization light-curve. We find a strong preferred position angle in optical polarization. The preferred position angle can be explained by separating the jet emission into two components: an optical polarization core and chaotic jet emission. The optical polarization core is stable on time scales of years and can be explained as emission from an underlying quiescent jet component. The chaotic jet emission sometimes exhibits a circular movement in the Stokes plane. We interpret these events as a shock front moving forwards and backwards in the jet, swiping through a helical magnetic field. Secondly, we use our data to assess different binary black hole models proposed to explain the regularly appearing double-peaked bursts in OJ287. We compose a list of requirements a model has to fulfil. The list includes not only characteristics of the light-curve but also other properties of OJ287, such as the black hole mass and restrictions on accretion flow properties. We rate all existing models using this list and conclude that none of the models is able to explain all observations. We discuss possible new explanations and propose a new approach to understanding OJ287. We suggest that both the double-peaked bursts and the evolution of the optical polarization position angle could be explained as a sign of resonant accretion of magnetic field lines, a 'magnetic breathing' of the disc.