Search

We present analyses of Suzaku XIS light curves and spectra of the BL Lac object OJ 287 with observations positioned primarily around proposed recurrent optical outbursts. The first two observations were performed in 2007 April 10 - 13 (epoch 1) and 2007 November 7 - 9 (epoch 2) that respectively correspond to a low and a high optical state and which, within the binary supermassive black hole model for OJ 287, precede and follow the impact flare. The last three observations, made consecutively during 2015 May 3 - 9 (epoch 3), were during the post-impact state of the 2013 disc impact and are the longest continuous X-ray observation of OJ 287 taken before the optical outburst in 2015 December. Intraday variability is found in both the soft (0.5 - 2 keV) and hard (2 - 10 keV) bands. The discrete correction function analysis of the light curves in both bands peaks at zero lag during epochs 2 and 3, indicating that the emission in both bands was cospatial and emitted from the same population of leptons. Power spectral densities of all three light curves are red noise dominated, with a rather wide range of power spectrum slopes. These X-ray spectra are overall consistent with power-laws but with significantly different spectral indices. In the 2015 observations the X-ray spectrum softens during the flare, showing an obvious soft X-ray excess that was not evident in the 2007 observations. We discuss the implications of these observations on the jet, the possible accretion disc, and the binary supermassive black hole model proposed for the nearly periodic optical flaring of OJ 287., Comment: 15 pages,9 figures,9 tables

To investigate the short-term variability and determine the size of the optical continuum emitting region of intermediate-mass black holes (IMBHs), we carried out high-cadence, multiband photometric monitoring of a Seyfert 1 galaxy J0249-0815 across two nights, together with a one-night single-band preliminary test. The presence of the broad Ha component in our target was confirmed by recent Paloma spectroscopic observations, 23 years after the Sloan Digital Sky Survey, ruling out the supernovae origin of the broad Ha line. The photometric experiment was primarily conducted utilizing four-channel imagers MuSCAT 3 and 4 mounted on 2 m telescopes within the Las Cumbres Observatory Global Telescope Network. Despite the expectation of variability, we observed no significant variation (<1.4%) on timescales of 6 to 10 hr. This nondetection is likely due to substantial host galaxy light diluting the subtle active galactic nucleus (AGN) variability. Additionally, we cannot rule out that the target was in a relatively quiescent state without intranight variability during our monitoring, owing to the stochastic nature of AGN variations. To enhance the possibility of detecting subtle variability signals and lag in future IMBH reverberation campaigns, it may be beneficial to select targets with a higher AGN-to-host flux ratio, and conduct dual-band preliminary tests and tailored simulations., Comment: 11 pages, 6 figures

We report the study of a huge optical intraday flare on November 12, 2021, at 2 am UT, in the blazar OJ287. In the binary black hole model it is associated with an impact of the secondary black hole on the accretion disk of the primary. Our multifrequency observing campaign was set up to search for such a signature of the impact, based on a prediction made eight years earlier. The first I-band results of the flare have already been reported by \cite{2024ApJ...960...11K}. Here we combine these data with our monitoring in the R-band. There is a big change in the R-I spectral index by $1.0\pm0.1$ between the normal background and the flare, suggesting a new component of radiation. The polarization variation during the rise of the flare suggests the same. The limits on the source size place it most reasonably in the jet of the secondary black hole. We then ask why we have not seen this phenomenon before. We show that OJ287 was never before observed with sufficient sensitivity on the night when the flare should have happened according to the binary model. We also study the probability that this flare is just an oversized example of intraday variability, using the Krakow-dataset of intense monitoring between 2015 and 2023. We find that the occurrence of a flare of this size and rapidity is unlikely. In the Appendix, we give the full orbit-linked historical light curve of OJ287 as well as the dense monitoring sample of Krakow., Comment: to Appear in Astrophysical Journal Letters

The extreme TeV blazar 1ES 0229+200 is a high-frequency-peaked BL Lacertae object. It has not shown intraday variability in extensive optical and X-ray observations. Nor has it shown any significant variability on any measurable timescale in the 1-100 GeV energy range over a 14-year span, but variations in the source flux around its average are present in the energy range above 200 GeV. We searched for intraday optical variability in 1ES 0229+200 as part of an ongoing project to search for variability and quasi-periodic oscillations in the high-cadence, nearly uniformly sampled optical light curves of blazars provided by the Transiting Exoplanet Survey Satellite (TESS). 1ES 0229+200 was monitored by TESS in its Sectors 42, 43 and 44. We analysed the data of all these sectors both with the TESS provided lightkurve software and the eleanor reduction pipeline. We detected a strong, essentially symmetric flare that lasted ~6 hours in Sector 42. We fit the flare's rising and declining phases to exponential functions. We also analysed the light curve of Sector 42 using the Lomb-Scargle periodogram (LSP) and continuous auto-regressive moving average (CARMA) methods. The Sector 42 light curve displayed in the present work provides the first evidence of a strong, rapid, short-lived optical flare on the intraday timescale in 1ES 0229+200. The variability timescale of the flare provides the upper limit for the size of the emission region to be within (3.3\pm0.2 - 8.3\pm0.5)x10^{15}cm. Away from the flare, the slope of the periodogram's power spectrum is fairly typical of many blazars (\alpha<2), but the nominal slopes for the flaring regions are very steep (\alpha~4.3), which may indicate the electron distribution undergoes a sudden change. We discuss possible emission mechanisms that could explain this substantial and rapid flare., Comment: 7 pages, 4 figures, 3 tables, Accepted for publication in A&A

This contribution presents the recent research developments within the Cosmic Ray Extremely Distributed Observatory (CREDO) in the search for resolution of various scientific puzzles, ranging from fundamental physical questions to applications like the determination of earthquake precursors. The state-of-the art theoretical, numerical and computational aspects of these phenomena are addressed, as well as recent experimental developments for detection., Comment: 9 pages, 7 figures. Based on the contribution to the XVIII Mexican Workshop on Particles and Fields in Puebla, Mexico, November 2022. M\'exico

We present an extensive recent multi-band optical photometric observations of the blazar S5 0716+714 carried out over 53 nights with two telescopes in India, two in Bulgaria, one in Serbia, and one in Egypt during 2019 November -- 2022 December. We collected 1401, 689, 14726, and 165 photometric image frames in B, V, R, and I bands, respectively. We montiored the blazar quasi-simultaneously during 3 nights in B, V, R, and I bands; 4 nights in B, V, and R; 2 nights in V, R, and I; 5 nights in B and R; and 2 nights in V and R bands. We also took 37 nights of data only in R band. Single band data are used to study intraday flux variability and two or more bands quasi-simultaneous observations allow us to search for colour variation in the source. We employ the power-enhanced F-test and the nested ANOVA test to search for genuine flux and color variations in the light curves of the blazar on intraday timescales. Out of 12, 11, 53, and 5 nights observations, intraday variations with amplitudes between ~3% and ~20% are detected in 9, 8, 31 and 3 nights in B, V, R, and I bands, respectively, corresponding to duty cycles of 75%, 73%, 58% and 60%. These duty cycles are lower than those typically measured at earlier times. On these timescales color variations with both bluer-when-brighter and redder-when-brighter are seen, though nights with no measurable colour variation are also present. We briefly discuss possible explanations for this observed intraday variability., Comment: 19 pages, 5 figures, 4 tables, Accepted for Publication in MNRAS

We present the results of our study of cross-correlations between long-term multi-band observations of the radio variability of the blazar 3C 279. More than a decade (2008-2022) of radio data were collected at seven different frequencies ranging from 2 GHz to 230 GHz. The multi-band radio light curves show variations in flux, with the prominent flare features appearing first at higher-frequency and later in lower-frequency bands. This behavior is quantified by cross-correlation analysis, which finds that the emission at lower-frequency bands lags that at higher-frequency bands. Lag versus frequency plots are well fit by straight lines with negative slope, typically ~-30 day/GHz. We discuss these flux variations in conjunction with the evolution of bright moving knots seen in multi-epoch VLBA maps to suggest possible physical changes in the jet that can explain the observational results. Some of the variations are consistent with the predictions of shock models, while others are better explained by a changing Doppler beaming factor as the knot trajectory bends slightly, given a small viewing angle to the jet., Comment: Submitted revised version to MNRAS journal, 11 pages, 6 figures, 4 tables

We study the optical flux and polarization variability of the binary black hole blazar OJ 287 using quasi-simultaneous observations from 2015 to 2023 carried out using telescopes in the USA, Japan, Russia, Crimea, and Bulgaria. This is one of the most extensive quasi-simultaneous optical flux and polarization variability studies of OJ 287. OJ 287 showed large amplitude, ~3.0 mag flux variability, large changes of ~37% in degree of polarization, and a large swing of ~215 degrees in the angle of the electric vector of polarization. During the period of observation, several flares in flux were detected. Those flares are correlated with a rapid increase in the degree of polarization and swings in electric vector of polarization angle. A peculiar behavior of anticorrelation between flux and polarization degree, accompanied by a nearly constant polarization angle, was detected from JD 2,458,156 to JD 2,458,292. We briefly discuss some explanations for the flux and polarization variations observed in OJ 287., Comment: 14 pages, 4 figures, 2 tables, Published in Astrophysical Journal Letters, 2023 November 1

We have analyzed the optical light curves of the blazar OJ 287 obtained with the Transiting Exoplanet Survey Satellite (TESS) over about 80 days from 2021 October 13 to December 31, with an unprecedented sampling of 2 minutes. Although significant variability has been found during the entire period, we have detected two exceptional flares with flux nearly doubling and then nearly tripling over 2 days in the middle of 2021 November. We went through the light curves analysis using the excess variance, generalized Lomb-Scargle periodogram, and Continuous Auto-Regressive Moving Average (CARMA) methods, and estimated the flux halving/doubling timescales. The most probable shortest variability timescale was found to be 0.38 days in the rising phase of the first flare. We briefly discuss some emission models for the variability in radio-loud active galactic nuclei that could be capable of producing such fast flares., Comment: 10 pages, 6 figures. Accepted for publication in Astrophysical Journal

We report the first extensive optical flux and spectral variability study of the TeV blazar TXS 0506+056 on intra-night to long-term timescales using BVRI data collected over 220 nights between January 21, 2017 to April 9, 2022 using 8 optical ground-based telescopes. In our search for intraday variability (IDV), we have employed two statistical analysis techniques, the nested ANOVA test and the power enhanced F-test. We found the source was variable in 8 nights out of 35 in the R-band and in 2 of 14 in the V-band yielding Duty Cycles (DC) of 22.8% and 14.3%, respectively. Clear colour variation in V - R was seen in only 1 out of 14 observing nights, but no IDV was found in the more limited B, I, and B - I data. During our monitoring period the source showed a 1.18 mag variation in the R-band and similar variations are clearly seen at all optical wavelengths. We extracted the optical (BVRI) SEDs of the blazar for 44 nights when observations were carried out in all four of those wavebands. The mean spectral index (\alpha) was determined to be 0.897+-0.171, Comment: 13 pages, 8 figures, 8 tables, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2212.09016

Aims. We aim to investigate the extreme variability properties of the TeV blazar S4 0954+65 using optical photometric and polarisation observations carried out between 2017 and 2023 using three ground-based telescopes. Methods. We examined an extensive dataset comprised of 138 intraday (observing duration shorter than a day) light curves (LCs) of S4 0954+65 for flux, spectral, and polarisation variations on diverse timescales. For the variable LCs, we computed the minimum variability timescales. We investigated flux-flux correlations and colour variations to look for spectral variations on long (several weeks to years) timescales. Additionally, we looked for connections between optical R-band flux and polarisation degree. Results. We found significant variations in 59 out of 138 intraday LCs. We detected a maximum change of 0.58$\pm$0.11 in V-band magnitude within $\sim$2.64 h and a corresponding minimum variability timescale of 18.21$\pm$4.87 mins on 2017 March 25. During the course of our observing campaign, the source brightness changed by $\sim$4 magnitudes in V and R bands; however, we did not find any strong spectral variations. The slope of the relative spectral energy distribution was 1.37$\pm$0.04. The degree of polarisation varied from $\sim$ 3% to 39% during our monitoring. We observed a change of $\sim$120 degrees in polarisation angle (PA) within $\sim$3 h on 2022 April 13. No clear correlation was found between optical flux and the degree of polarisation. Conclusions. The results of our optical flux, colour, and polarisation study provide hints that turbulence in the relativistic jet could be responsible for the intraday optical variations in the blazar S4 0954+65. However, the long-term flux variations may be caused by changes in the Doppler factor., Comment: 9 pages, 10 figures, 4 tables, 4 appendix, Astronomy & Astrophysics journal (in press)

In this contribution we study the possibility of the formation of cosmic ray ensembles (CRE) created by the interaction of ultra-high energy (UHE) photons with the magnetic field of the Sun. The lack of observation of those UHE and the difficulties for their identification given the current methodologies motivates this study. We performed simulations using the PRESHOWER program in order to simulate the expected extensive air showers which might be spatially correlated generated upon entering the Earth's atmosphere. We found characteristic features like very thing and extremely elongates cascades of secondary photons with their corresponding energies spanning the entire cosmic range spectrum. Shower footprints are as large as hundreds of kilometres. An application of this study is the scenario of gamma-ray emission from the vicinity of the Sun as a result of ultra-high energy photon cascading in the solar magnetic field in order to understand recent observations made by the HAWC and Fermi-LAT observatories., Comment: 8 pages, 4 figures, 38th International Cosmic Ray Conference (ICRC2023) 26 July - 3 August, 2023 Nagoya, Japan

We present a variability study of the blazar PKS 0346-27 from December 2018 to January 2022 in its archival $\gamma$-ray observation by Fermi-LAT. We use the Lomb-Scargle periodogram and the weighted wavelet transform methods in order to detect the presence of periodicity/quasi-periodicity and localize this feature in time and frequency space. The significance of the periodicity feature has been estimated using the Monte-Carlo simulation approach. We have also determined the global significance of the periodicity to test the robustness of our claim. To explore the most probable scenario, we modeled the light curve with both a straight jet and a curved jet model. We detect a periodicity feature of $\sim$ 100 days duration for the entire period of observation with a statistical significance of $3\sigma$, which amounts to a 99.7\% confidence level. The global significance of this feature is found to be 96.96\%. Based on the Akaike Information Criteria, the most probable explanation is that the observed emission is enhanced due to the helical motion of a blob within a curved jet. The origin of this QPO is very likely a region of enhanced emission moving helically inside a curved jet. This work presents strong evidence for jet curvature in the source and an independent (albeit a little serendipitous) procedure to estimate the curvature in blazar jets., Comment: 7 pages, 4 figures, Accepted for publication in A&A, in Press

The highly variable blazar OJ~287 is commonly discussed as an example of a binary black hole system. The 130 year long optical light curve is well explained by a model where the central body is a massive black hole of 18.35$\times$10$^9$ solar mass that supports a thin accretion disc. The secondary black hole of 0.15$\times$10$^9$ solar mass impacts the disc twice during its 12 year orbit, and causes observable flares. Recently, it has been argued that an accretion disc with a typical AGN accretion rate and above mentioned central body mass should be at least six magnitudes brighter than OJ~287's host galaxy and would therefore be observationally excluded. Based on the observations of OJ~287's radio jet, detailed in Marscher and Jorstad (2011), and up-to-date accretion disc models of Azadi et al. (2022), we show that the V-band magnitude of the accretion disc is unlikely to exceed the host galaxy brightness by more than one magnitude, and could well be fainter than the host. This is because accretion power is necessary to launch the jet as well as to create electromagnetic radiation, distributed across many wavelengths, and not concentrated especially on the optical V-band. Further, we note that the claimed V-band concentration of accretion power leads to serious problems while interpreting observations of other Active Galactic Nuclei. Therefore, we infer that the mass of the primary black hole and its accretion rate do not need to be smaller than what is determined in the standard model for OJ~287.

We present a summary of the results of the OJ 287 observational campaign, which was carried out during the 2021/2022 observational season. This season is special in the binary model because the major axis of the precessing binary happens to lie almost exactly in the plane of the accretion disc of the primary. This leads to pairs of almost identical impacts between the secondary black hole and the accretion disk in 2005 and 2022. In 2005, a special flare called "blue flash" was observed 35 days after the disk impact, which should have also been verifiable in 2022. We did observe a similar flash and were able to obtain more details of its properties. We describe this in the framework of expanding cloud models. In addition, we were able to identify the flare arising exactly at the time of the disc crossing from its photo-polarimetric and gamma-ray properties. This is an important identification, as it directly confirms the orbit model. Moreover, we saw a huge flare that lasted only one day. We may understand this as the lighting up of the jet of the secondary black hole when its Roche lobe is suddenly flooded by the gas from the primary disk. Therefore, this may be the first time we directly observed the secondary black hole in the OJ 287 binary system., Comment: 24 pages, 13 figures, published in Galaxies as part of the Special Issue Distant Glowing Objects: Quest for Quasars

We perform correlation and periodicity search analyses on long-term multi-band light curves of the FSRQ 1510-089 observed by the space-based Fermi--Large Area Telescope in gamma-rays, the SMARTS and Steward Observatory telescopes in optical and near-infrared (NIR) and the 13.7 m radio telescope in Metsahovi Radio Observatory between 2008 and 2018. The z-transform discrete correlation function method is applied to study the correlation and possible time lags among these multi band light curves. Among all pairs of wavelengths, the gamma-ray vs. optical/NIR and optical vs. NIR correlations show zero time lags; however, both the gamma-ray and optical/NIR emissions precede the radio radiation. The Generalized Lomb-Scargle periodogram, Weighted Wavelet Z-transform, and REDFIT techniques are employed to investigate the unresolved-core-emission dominated 37 GHz light curve and yield evidence for a quasi-period around 1540 days, although given the length of the whole data set it cannot be claimed to be significant. We also investigate the optical/NIR color variability and find that this source shows a simple redder-when-brighter behavior over time, even in the low flux state., Comment: Accepted for publication in ApJ; 20 pages, 9 figures, 4 tables

In this work, we report the presence of rapid intra-night optical variations in both -- flux and polarization of the blazar BL Lacertae during its unprecedented 2020--2021 high state of brightness. The object showed significant flux variability and some color changes, but no firmly detectable time delays between the optical bands. The linear polarization was also highly variable in both -- polarization degree and angle (EVPA). The object was observed from several observatories throughout the world, covering in a total of almost 300 hours during 66 nights. Based on our results, we suggest, that the changing Doppler factor of an ensemble of independent emitting regions, travelling along a curved jet that at some point happens to be closely aligned with the line of sight can successfully reproduce our observations during this outburst. This is one of the most extensive variability studies of the optical polarization of a blazar on intra-night timescales., Comment: 23 pages,7 figures, 5 Tables (2 as appendix). Accepted for publication in MNRAS

Here we present the first optical photometric monitoring results of a sample of twelve newly discovered blazars from the ICRF - Gaia CRF astrometric link. The observations were performed from April 2013 until August 2019 using eight telescopes located in Europe. For a robust test for the brightness and colour variability, we use Abbe criterion and F-test. Moreover, linear fittings are performed to investigate the relation in the colour-magnitude variations of the blazars. Variability was confirmed in the case of 10 sources; two sources, 1429+249 and 1556+335 seem to be possibly variable. Three sources (1034+574, 1722+119, and 1741+597) have displayed large amplitude brightness change of more than one magnitude. We found that the seven sources displayed bluer-when-brighter variations, and one source showed redder-when-brighter variations. We briefly explain the various AGN emission models which can explain our results., Comment: 14 pages, 8 figures, 11 tables, 6 supplement figures can be provided on request, MNRAS in press

The bright blazar OJ~287 routinely parades high brightness bremsstrahlung flares, which are explained as being a result of a secondary supermassive black hole (SMBH) impacting the accretion disc of a more massive primary SMBH in a binary system. The accretion disc is not rigid but rather bends in a calculable way due to the tidal influence of the secondary. Below we refer to this phenomenon as a variable disc level. We begin by showing that these flares occur at times predicted by a simple analytical formula, based on general relativity inspired modified Kepler equation, which explains impact flares since 1888. The 2022 impact flare, namely flare number 26, is rather peculiar as it breaks the typical pattern of two impact flares per 12-year cycle. This is the third bremsstrahlung flare of the current cycle that follows the already observed 2015 and 2019 impact flares from OJ~287. It turns out that the arrival epoch of flare number 26 is sensitive to the level of primary SMBH's accretion disc relative to its mean level in our model. We incorporate these tidally induced changes in the level of the accretion disc to infer that the thermal flare should have occurred during July-August 2022, when it was not possible to observe it from the Earth. Thereafter, we explore possible observational evidence for certain pre-flare activity by employing spectral and polarimetric data from our campaigns in 2004/05 and 2021/22. We point out theoretical and observational implications of two observed mini-flares during January-February 2022., Comment: 29 pages, 6 figures, 1 table. arXiv admin note: text overlap with arXiv:2209.08360

We present a long-term and intraday variability study on optical multiwaveband ($U\!BV\!RI$) data from the blazar AO 0235+164 collected by various telescopes for $\sim$44 years (1975--2019). The blazar was found to be significantly variable over the years in all wavebands with a variation of about six magnitudes between its low and active states. The variations in the different wavebands are highly correlated without any time lag. We did not observe any significant trend in color variation with time, but we observed a bluer-when-brighter trend between the $B-I$ color index and the $R$-magnitude. Optical $BV\!R$-band spectral energy distributions always show a convex shape. Significant intraday variability was frequently seen in the quasi-simultaneous observations of AO\,0235+164 made on 22 nights in $R$ and $V$-bands by the CASLEO and CAHA telescopes during 1999--2019. We also estimated the central supermassive black-hole mass of $7.9\times10^7 M_{\odot}$ by analyzing the broad Mg II emission line in AO 0235+164's spectrum. We briefly explore the probable physical scenarios responsible for the observed variability., Comment: 27 pages, 10 figures. Accepted for publication in The Astrophysical Journal Supplement Series (ApJS)

We report the flux and spectral variability of PG 1553+113 on intra-night (IDV) to short-term timescales using BVRI data collected over 91 nights from 28 February to 8 November 2019 employing ten optical telescopes: three in Bulgaria, two each in India and Serbia, and one each in Greece, Georgia, and Latvia. We monitored the blazar quasi-simultaneously for 16 nights in the V and R bands and 8 nights in the V, R, I bands and examined the light curves (LCs) for intra-day flux and colour variations using two powerful tests: the power-enhanced F-test and the nested ANOVA test. The source was found to be significantly (> 99%) variable in 4 nights out of 27 in R-band, 1 out of 16 in V-band, and 1 out of 6 nights in I-band. No temporal variations in the colours were observed on IDV timescale. During the course of these observations the total variation in R-band was 0.89 mag observed. We also investigated the spectral energy distribution (SED) using B, V, R, and I band data. We found optical spectral indices in the range of 0.878+-0.029 to 1.106+-0.065 by fitting a power law to these SEDs of PG 1553+113. We found that the source follows a bluer-when-brighter trend on IDV timescales. We discuss possible physical causes of the observed spectral variability., Comment: 13 pages,8 figures, 7 tables, Accepted for publication in MNRAS

We report the detection of several quasi-periodicities around 0.6--2.5 days in the optical emission of the blazar S4 0954+658. The source was observed by the Transiting Exoplanet Survey Satellite (TESS) in six sectors and it showed these features in all but one of them, with a QPO of 1.52 days apparently present in portions of four of them. We used the generalized Lomb-Scargle periodogram method to search for significant signals and we confirmed them using a weighted wavelet transform for time-frequency domain analyses. We discuss several possible explanations for these rapid quasi-periodic variations and suggest that an origin in the innermost part of the accretion disk is most likely. Within this framework, we provide estimates for the mass of the black hole at the core of this blazar., Comment: 20 pages, 6 figures, 4 tables

BL Lacertae had undergone a series of historical high flux activity over a year, from August 2020 in the optical to VHE $\gamma$-rays. In this paper, we report on optical flux and spectral variability of the first historical maxima outburst event during October-November in g, r and i bands with the 1.26m telescope at Xinglong observatory, China. We detected significant intranight variations with amplitude rising up to $\sim 30$%, when the fastest variability timescale is found to be a few tens of minutes, giving an emitting region size of the order $10^{-3}$ pc, which corresponds to $\sim 100$ Schwarzschild radius of the central black hole, likely coming from some jet mini-structures. Unlike on intranight timescale, a clear frequency dependent pattern along with symmetric timescales ($\sim$ 11d) of flux variation are detected on long timescale. The spectral evolution was predominated by flattening of the spectra with increasing brightness i.e., a bluer-when-brighter trend in 96% of the cases. On the night before the outburst peak, the color indices clustered in two distinct branches in color--magnitude diagram within a period of $\sim$ 6 hours that is connected to a hard-soft-hard spectral evolution trend extracted from time-resolved spectra. Such trend has never seen in BL Lac or any other blazars before to the best of our knowledge. The results obtained in this study can be explained in the context of shock induced particle acceleration or magnetic re-connection in the jet where turbulent processes most likely resulted the asymmetric flux variation on nightly timescale., Comment: 15 pages, 3 tables, 8 figures. Accepted for publication in ApJ

The continuous monitoring capability of Fermi-LAT has enabled the exploration of Quasi-Periodic Oscillations (QPOs) in the $\gamma$-ray light curve of blazar that has given a new perspective to probe these sources over a wide range of time scales. We report the presence of transient QPOs in the long-term $\gamma$-ray light curve of blazars PKS 0244-470 and 4C +38.41. We first identified different flux states using the Bayesian Block algorithm and then explored the possible transient QPOs in the segments of each flux phase where the flux level changes over fairly regular intervals. Combining this with the source's intrinsic variance, we identified two flux phases for PKS 0244-470: one activity (AP-1) and one quiescent phase (QP-1). For 4C+38.41, we similarly identified four activity (AP-1, AP-2, AP-3, AP-4) and two quiescent (QP-1, QP-2) phases. AP-1 phase of PKS 0244-470 shows QPO of $\sim$225 days persisting for 8 cycles ($\sim$ 4.1$\sigma$). In 4C+38.41, AP-1 and AP-2 phases show QPO-like behavior of $\sim$110 days and $\sim$ 60 days, respectively, persisting for 5 cycles. In AP-3, we identified three sub-phases, and all show a $\sim$ week scale possible recurrent rise with five complete cycles, while in QP-1, we could identify 2 sub-phases (Q1, and Q2). Q1 phase shows a period of $\sim$ 104 days with six complete cycles. Q2 phase also shows QPO but with only $\sim$3.7 cycles. We discuss the possible origin and argue that the current driven kink instability and curved jet model seem the most likely cause for shorter and longer QPOs., Comment: Accepted for publication in ApJ

We present an extensive study on the X-ray intraday variability of ten TeV-emitting high synchrotron peaked blazars (HBLs): 1ES 0229+200, 1ES 0414+009, PKS 0548-322, 1ES 1101-232, 1H 1219+301, H 1426+428, Mrk 501, 1ES 1959+650, PKS 2005-489, and 1ES 2344+514 made with twenty-five XMM-Newton pointed observations during its operational period. Intraday variability has been estimated in three energy bands: soft (0.3--2 keV), hard (2--10 keV) and total (0.3--10 keV). Although seven out of these ten TeV HBLs exhibited some intraday variability at three-sigma levels no major variations exceeding six percent were detected. We explored the spectral properties of the sample by extracting the hardness ratio from the soft and hard bands; no significant variations in the hardness ratio were observed in any source. We performed power spectral density analyses on the variable light-curves by fitting power-laws, yielding slopes lying in the range from 1.11 to 2.93 for different HBLs. We briefly discuss possible emission mechanisms and carry out rough estimates for magnetic fields, electron Lorentz factors and emission region sizes for seven of these HBLs., Comment: Accepeted for publication in ApJ, 33 pages, 14 figures, 5 tables

Highly variable Markarian 421 is a bright high synchrotron energy peaked blazar showing wide featureless non-thermal spectrum making it a good candidate for our study of intraday flux and spectral variations over time. We analyse its X-ray observations of over 17 years taken with the EPIC-PN instrument to probe into the intraday variability properties. The photon energy band of 0.3 - 10.0 keV, and its sub-bands, soft 0.3-2.0 keV and hard 2.0-10.0 keV. To examine flux variability, fractional variability amplitude and the minimum variability timescale have been calculated. We also probed into the spectral variability by studying hardness ratio for each observation and the correlation between the two energy bands using discrete correlation function and inspecting the normalized light curves. The parameters obtained from these methods have been studied for any correlation or non-random trends. From this work, we speculate on the constraints for possible particle acceleration and emission processes in the jet, for better understanding of the processes involving a turbulent behaviour except of shocks. A positive discrete correlation function between the two sub-bands indicates the role of the same electron population in the emission of photons in the two bands. The correlation between the parameters of flux variability and parameters of spectral variation and lags in sub-energy bands provide the constraints to be considered for any modelling of emission processes., Comment: Accepted in APJS on 31 May 2022; Online supplementary material of paper placed in ancillary material on arxiv

We present the most extensive and well-sampled long-term multi-band near-infrared (NIR) temporal and spectral variability study of OJ 287, considered to be the best candidate binary supermassive black hole blazar. These observations were made between December 2007 and November 2021. The source underwent ~ 2 -- 2.5 magnitude variations in the J, H, and Ks NIR bands. Over these long-term timescales there were no systematic trends in either flux or spectral evolution with time or with the source's flux states. However, on shorter timescales, there are significant variations in flux and spectra indicative of strong changes during different activity states. The NIR spectral energy distributions show diverse facets at each flux state, from the lowest to the highest. The spectra are, in general, consistent with a power-law spectral profile (within 10%) and many of them indicate minor changes (observationally insignificant) in the shift of the peak. The NIR spectra generally steepens during bright phases. We briefly discuss these behaviors in the context of blazar emission scenarios/mechanisms, OJ 287's well-known traditional behavior, and implications for models of the source central engine invoked for its long-term optical semi-periodic variations., Comment: 17 pages, 6 figures, 5 tables. Accepted for publication in ApJS

We present a long term optical $R$ band light curve analysis of the gravitationally lensed blazar AO 0235+164 in the time span 1982 - 2019. Several methods of analysis lead to the result that there is a periodicity of ~8.13 years present in these data. In addition, each of these five major flares are apparently double-peaked, with the secondary peak following the primary one by ~2 years. Along with the well known system, OJ 287, our finding constitutes one of the most secure cases of long term quasi-periodic optical behaviour in a blazar ever found. A binary supermassive black hole system appears to provide a good explanation for these results., Comment: 7 pages, 3 figures, 1 table, Accepted for publication in MNRAS

Old, digitized astronomical images taken before the human spacefaring age offer a unique view of the sky devoid of known artificial satellites. In this paper, we have carried out the first optical searches ever for non-terrestrial artifacts near the Earth following the method proposed in Villarroel et al. (2022). We use images contained in the First Palomar Sky Survey to search for simultaneous (during a plate exposure time) transients that in addition to being point-like, are aligned. We provide a shortlist of the most promising candidates of aligned transients, that must be examined with the help of a microscope to separate celestial sources from plate defects with coincidentally star-like brightness profiles. We further explore one possible, but not unique, interpretation in terms of fast reflections off high-albedo objects in geosynchronous orbits around Earth. If a future study rules out each multiple transient candidate, the estimated surface density becomes an upper limit of $<10^{-9}$ objects km$^{-2}$ non-terrestrial artifacts in geosynchronous orbits around Earth. Finally, we conclude that observations and analysis of multiple, simultaneously appearing and vanishing light sources on the sky merit serious further attention, regardless of their origin., Comment: Submitted to Astronomical Journal

Cosmic Ray Ensembles (CRE) are yet not observed groups of cosmic rays with a common primary interaction vertex or the same parent particle. One of the processes capable of initiating identifiable CRE is an interaction of an ultra-high energy (UHE) photon with the solar magnetic field which results in an electron pair production and the subsequent synchrotron radiation. The resultant electromagnetic cascade forms a very characteristic line-like front of a very small width ($\sim$ meters), stretching from tens of thousands to even many millions of kilometers. In this contribution we present the results of applying a toy model to simulate detections of such CRE at the ground level with an array of ideal detectors of different dimensions. The adopted approach allows us to assess the CRE detection feasibility for a specific configuration of a detector array. The process of initiation and propagation of an electromagnetic cascade originated from an UHE photon passing near the Sun, as well as the resultant particle distribution on ground, were simulated using the CORSIKA program with the PRESHOWER option, both modified accordingly. The studied scenario results in photons forming a cascade that extends even over tens of millions of kilometers when it arrives at the top of the Earth's atmosphere, and the photon energies span practically the whole cosmic ray energy spectrum. The topology of the signal consists of very extended CRE shapes, and the characteristic, very much elongated disk-shape of the particle distribution on ground illustrates the potential for identification of CRE of this type., Comment: 8 pages, 3 figures. Proceedings of the International Cosmic Rays Conference 2021, 12-23 July, Berlin, Germany

We present periodicity search analyses on the $\gamma$-ray lightcurve of the TeV blazar PKS 1510-089 observed by the \textit{Fermi} Large Area Telescope. We report the detection of two transient quasi-periodic oscillations: a 3.6-day QPO during the outburst in 2009 that lasted five cycles (MJD 54906--54923); and a periodicity of 92 days spanning over 650 days from 2018 to 2020 (MJD 58200--58850), which lasted for seven cycles. We employed the Lomb-Scargle periodogram, Weighted Wavelet Z-transform, REDFIT, and the Monte Carlo lightcurve simulation techniques to find any periodicity and the corresponding significance. The 3.6-day QPO was detected at a moderate significance of $\sim$3.5$\sigma$, while the detection significance of the 92-day QPO was $\sim$7.0$\sigma$. We explore a few physical models for such transient QPOs including a binary black hole system, precession of the jet, a non-axisymmetric instability rotating around the central black hole near the innermost stable circular orbit, the presence of quasi-equidistant magnetic islands inside the jet, and a geometric model involving a plasma blob moving helically inside a curved jet., Comment: 10 pages, 4 figures, submitted to MNRAS

We present X-ray intraday variability and power spectral density (PSD) analyses of the longest 23 pointed XMM-Newton observations of the blazar 3C 273 that were taken during 2000 -- 2021. These good time intervals contained between 5 and 24.6 hours of data. Variability has been estimated in three energy bands: 0.2 -- 2 keV (soft), 2 -- 10 keV (hard), and 0.2 -- 10 keV (total). Nine of the 23 observations exhibited some variability, though no major variations exceeding 5 per cent were detected. Typical timescales for variability were $\sim 1$ ks. For those variable light curves we find that a power-law model provides good fits to each PSD, with most of the slopes between $-1.7$ and $-2.8$. Although no variations of hardness ratio could be measured in any individual observation, an anti-correlation in flux and hardness ratio is found in long term data that indicates a harder when brighter trend. Our flux and spectral analyses indicate that both particle acceleration and synchrotron cooling processes make an important contribution to the emission from this blazar., Comment: 12 pages, 4 tables, 10 figures, Complete Fig. 1 - 6 will be published as supplementary material and can be provided on request; Accepted for Publication in MNRAS

In an optical monitoring program to characterize the variability properties of blazar, we observed 10 sources from the Roma-BZCAT catalogue for 26 nights in V and R bands during October 2014 to June 2015 with two telescopes located in India. The sample includes mainly newly discovered BL Lacs where the redshift of some sources are not known yet. We present the results of flux and color variations of the sample on intraday and short time scales obtained by using the power-enhanced F-test and the nested-ANOVA tests, along with their spectral behavior. We find significant intraday variability in the single FSRQ in our sample, having an amplitude of variation ~12%. Although a few of BL Lacs showed probable variation in some nights, none of them passes the variability tests at 99.9% significance level. We find that 78% of the sample showed significant negative colour--magnitude correlations i.e., a redder-when-brighter spectral evolution. Those which do not show strong or clear chromatism, predominantly exhibit a redder-when-brighter trends. Unlike on hourly timescales, the high synchrotron peaked (HSP) blazars in the sample (BZGJ0656+4237, BZGJ0152+0147 and BZBJ1728+5013) show strong flux variation on days to months timescales, where again we detect a decreasing trend of the spectral slope with brightness. We observe a global steepening of the optical spectrum with increasing flux on intranight timescale for the entire blazar sample. Non-variability in the BL Lacs in our sample could be resulted by distinct contribution from the disk as well as from other components in the studied energy range., Comment: 19 pages, 4 figures, 6 tables; accepted for publication in The Astrophysical Journal Supplement Series

Nine point sources appeared within half an hour on a region within $\sim$ 10 arcmin of a red-sensitive photographic plate taken in April 1950 as part of the historic Palomar Sky Survey. All nine sources are absent on both previous and later photographic images, and absent in modern surveys with CCD detectors which go several magnitudes deeper. We present deep CCD images with the 10.4-meter Gran Telescopio Canarias (GTC), reaching brightness $r \sim 26$ mag, that reveal possible optical counterparts, although these counterparts could equally well be just chance projections. The incidence of transients in the investigated photographic plate is far higher than expected from known detection rates of optical counterparts to e.g.\ flaring dwarf stars, Fast Radio Bursts (FRBs), Gamma Ray Bursts (GRBs) or microlensing events. One possible explanation is that the plates have been subjected to an unknown type of contamination producing mainly point sources with of varying intensities along with some mechanism of concentration within a radius of $\sim$ 10 arcmin on the plate. If contamination as an explanation can be fully excluded, another possibility is fast (t $<0.5$ s) solar reflections from objects near geosynchronous orbits. An alternative route to confirm the latter scenario is by looking for images from the First Palomar Sky Survey where multiple transients follow a line., Comment: 21 pages, 3 figures, 1 table. Published in Scientific Reports. Supplementary information can be found on the publishers webpage (open access)

We present variability analyses of twenty pointed XMM-Newton observations of the high energy peaked TeV blazar PG 1553+113 taken during 2010 to 2018. We found intraday variability in the total X-ray energy range (0.3 -- 10 keV) in 16 out of 19 light curves or a duty cycle of ~84%. A discrete correlation function analysis of the intraday light curves in the soft and hard X-ray bands peaks on zero lag, showing that the emission in hard and soft bands are co-spatial and emitted from the same population of leptons. Red-noise dominates the power spectral density (PSD) of all the LCs although the PSDs have a range of spectral slopes from -2.36 to -0.14. On longer timescales, the optical and UV variability patterns look almost identical and well correlated, as are the soft and hard X-ray bands, but the optical/UV variations are not correlated to those in the X-ray band, indicating that the optical/UV and X-ray emissions are emitted by two different populations of leptons. We briefly discuss physical mechanisms which may be capable of explaining the observed flux and spectral variability of PG 1553+113 on these diverse timescales., Comment: 11 pages, 4 tables, 7 figures, The redshift of the blazar has been updated, leading to minor changes in certain quantities but no changes to conclusions. Complete Fig. 1-3 will be published as supplementary material and can be provided on request

We have examined 13 pointed observations of the TeV emitting high synchrotron peak blazar PKS 2155-304, taken by the Suzaku satellite throughout its operational period. We found that the blazar showed large-amplitude intraday variabilities in the soft (0.8 - 1.5 keV) and the hard (1.5 - 8.0 keV) bands in the light curves. Spectral variability on intraday timescales is estimated using the hardness ratio. The blazar usually becomes harder when brighter and vice versa, following the typical behavior of high synchrotron peak blazars. The power spectral density (PSD) analyses of 11 out of 13 light curves in the total energy (0.8 - 8.0 keV) are found to be red-noise dominated, with power-law spectral indices that span a large range, from -2.81 to -0.88. Discrete correlation function analyses of all the 13 light curves between the soft and the hard bands show that they are well correlated and peak at, or very close to, zero lag. This indicates that the emissions in soft and hard bands are probably cospatial and emitted from the same population of leptons. Considering fluxes versus variability timescales, we found no correlation on intraday timescales, implying that X-ray emission from PKS 2155-304 is not dominated by simple changes in the Doppler factor. We briefly discuss the most likely emission mechanisms responsible for the observed flux and spectral variabilities and place constraints on magnetic field strength and Lorentz factors of the electrons emitting the X-rays in the most likely scenario., Comment: 18 pages, 7 figures, 4 tables; Accepted for publication in ApJ

We present a detailed timing and spectral study of an extremely variable narrow-line Seyfert~1 galaxy NGC 4748 using observations in the year 2017 and 2014 performed with AstroSat and XMM-Newton, respectively. Both observations show extremely variable soft and hard X-ray emission that are correlated with each other. In the 2014 data set, the source retains its general behaviour of "softer when brighter" while the 2017 observation exhibits a "harder when brighter" nature. Such changing behaviour is rare in AGNs and is usually observed in the black hole binary systems. The "harder when brighter" is confirmed with the anti-correlation between the photon index and the 0.3-10 keV power-law flux. This suggests a possible change in the accretion mode from standard to the advection-dominated flow. Additionally, both the observations show soft X-ray excess below 2 keV over the power-law continuum. This excess was fitted with a single or multiple blackbody component(s). The origin of soft excess during the 2017 observation is likely due to the cool Comptonization as the photon index changes with time. On the other hand, the broad iron line and delayed UV emission during the 2014 observation strongly suggest that X-ray illumination onto the accretion disk and reflection and reprocessing play a significant role in this AGN., Comment: 13 pages, 9 figures, 2 tables; accepted for publication in the Journal of Astrophysics and Astronomy

We present time series analyses of three-decade long radio observations of the BL Lacertae object AO 0235+164 made at the University of Michigan Radio Astronomical Observatory operating at three central frequencies of 4.8 GHz, 8.0 GHz and 14.5 GHz. We detected a quasi-periodic oscillation of $\sim$965 days in all three frequency bands in the light curve of the effectively simultaneous observations, along with strong signals at $\sim$1950 d, $\sim$1350 d, and $\sim$660 d. The periodicity is analyzed with three methods: Data Compensated Discrete Fourier Transform, Generalized Lomb-Scargle Periodogram and Weighted Wavelet Z-transform. These methods are chosen as they have different analysis approaches toward robust measurement of claimed periodicities. The QPO at $965\pm 50$ days is found to be significant (at least $3.5\sigma$) and is persistent throughout the observation for all three radio frequencies, and the others, which may be harmonics, are comparably significant in at least the 8.0 GHz and 14.5 GHz bands. We briefly discuss plausible explanations for the origin of such long and persistent periodicity., Comment: 12 pages, 4 figures. Accepted for publication in MNRAS

We report the re-emergence of a new broadband emission through a detailed and systematic study of the multi-wavelength spectral and temporal behavior of OJ 287 after its first-ever reported VHE activity in 2017 to date, which includes the second-highest X-ray flux of the source. The source shows high optical to X-ray flux variations, accompanied mainly by strong spectral changes. The optical to X-ray flux variations are correlated and simultaneous except for two durations when they are anti-correlated. The flux variations, however, are anti-correlated with the X-ray spectral state while correlated with optical-UV (ultraviolet). Weekly binned {\it Fermi}-LAT data around the duration of the highest X-ray activity show a few detections with a log-parabola model but none with a power-law; yet the extracted LAT spectral energy distribution (SED) of the high activity duration for both the models is similar and show a hardening above 1 GeV. Further, near-infrared (NIR) data indicate strong spectral change, resembling a thermal component. Overall, the combined optical to gamma-ray broadband spectrum establishes the observed variations to a new high-energy-peaked (HBL) broadband emission component, similar to the one seen during the highest reported X-ray flux state of the source in 2017. The observed activities indicate some peculiar features that seem to be characteristic of this emission component while its appearance, a few years around the claimed \(\sim 12\)-year optical outbursts strongly indicate a connection between the two., Comment: 14 pages, 6 figures, 1 table, ApJ (in press)

The Cosmic Ray Extremely Distributed Observatory (CREDO) is a newly formed, global collaboration dedicated to observing and studying cosmic rays (CR) and cosmic ray ensembles (CRE): groups of a minimum of two CR with a common primary interaction vertex or the same parent particle. The CREDO program embraces testing known CR and CRE scenarios, and preparing to observe unexpected physics, it is also suitable for multi-messenger and multi-mission applications. Perfectly matched to CREDO capabilities, CRE could be formed both within classical models (e.g. as products of photon-photon interactions), and exotic scenarios (e.g. as results of decay of Super Heavy Dark Matter particles). Their fronts might be significantly extended in space and time, and they might include cosmic rays of energies spanning the whole cosmic ray energy spectrum, with a footprint composed of at least two extensive air showers with correlated arrival directions and arrival times. Since CRE are mostly expected to be spread over large areas and, because of the expected wide energy range of the contributing particles, CRE detection might only be feasible when using available cosmic ray infrastructure collectively, i.e. as a globally extended network of detectors. Thus, with this review article, the CREDO Collaboration invites the astroparticle physics community to actively join or to contribute to the research dedicated to CRE, and in particular to share any cosmic ray data useful for the specific CRE detection strategies., Comment: 59 pages, 23 figures

Quasi-periodic fluctuations in the light curves of blazars can provide insight into the underlying emission process. This type of flux modulation hints at periodic physical processes that result in emission. CTA 102, a flat spectrum radio quasar at $z=1.032$, has displayed significant activity since 2016. The multi-waveband light curve of CTA 102 shows signs of quasi-periodic oscillations during the 2016-2017 flare. Our goal is to quantify the presence of any periodicity in the emitted flux during the mentioned period and to explore the causes that can give rise to it. Techniques such as the Lomb-Scargle periodogram and weighted wavelet z-transform were employed to observe the power emitted at different frequencies. To quantify the significance of the dominant period, Monte-Carlo techniques were employed considering an underlying smooth bending power-law model for the power spectrum. In addition, the light curve was modeled using an autoregressive process to analytically obtain the significance of the dominant period. Lastly, the light curve was modeled using a generalized autoregressive integrated moving average process to check whether introducing a periodic component results in a statistically preferable model. Significant, simultaneous quasi-periodic oscillations (QPOs) were observed in the $\gamma$-ray and optical fluxes of blazar CTA 102 during its 2016-2017 flare. The periodic flux modulation had a dominant period of $\sim$ 7.6 days and lasted for $\sim$ 8 cycles (MJD 57710-57770). All of the methods used point toward significant ($>4\sigma$) quasi-periodic modulation in both $\gamma$-ray and optical fluxes. Several possible models were explored while probing the origin of the periodicity, and by extension, the 2016--2017 optical flare. The best explanation for the detected QPO appears to be a region of enhanced emission, moving helically inside the jet., Comment: 9 pages, 4 figures. Accepted for publication in A&A

We report the detection of a probable $\gamma$-ray quasi-periodic oscillation (QPO) of around 314 days in the monthly binned 0.1 -- 300 GeV $\gamma$-ray {\it Fermi}-LAT light curve of the well known BL Lac blazar OJ 287. To identify and quantify the QPO nature of the $\gamma$-ray light curve of OJ 287, we used the Lomb-Scargle periodogram (LSP), REDFIT, and weighted wavelet z-transform (WWZ) analyses. We briefly discuss possible emission models for radio-loud active galactic nuclei (AGN) that can explain a $\gamma$-ray QPO of such a period in a blazar. Reports of changes in the position of quasi-stationary radio knots over a yearly timescale as well as a strong correlation between gamma-ray and mm-radio emission in previous studies indicate that the signal is probably associated with these knots., Comment: 6 pages, 4 figures, in press (MNRAS)

The Vanishing & Appearing Sources during a Century of Observations (VASCO) project investigates astronomical surveys spanning a time interval of 70 years, searching for unusual and exotic transients. We present herein the VASCO Citizen Science Project, which can identify unusual candidates driven by three different approaches: hypothesis, exploratory, and machine learning, which is particularly useful for SETI searches. To address the big data challenge, VASCO combines three methods: the Virtual Observatory, user-aided machine learning, and visual inspection through citizen science. Here we demonstrate the citizen science project and its improved candidate selection process, and we give a progress report. We also present the VASCO citizen science network led by amateur astronomy associations mainly located in Algeria, Cameroon, and Nigeria. At the moment of writing, the citizen science project has carefully examined 15,593 candidate image pairs in the data (ca. 10% of the candidates), and has so far identified 798 objects classified as "vanished". The most interesting candidates will be followed up with optical and infrared imaging, together with the observations by the most potent radio telescopes., Comment: 16 pages, 7 figures, contribution to IAC 2020 conference

We reviewed X-ray flux and spectral variability properties studied to date by various X-ray satellites for Mrk 421 and PKS 2155-304, which are TeV emitting blazars. Mrk 421 and PKS 2155-304 are the most X-ray luminous blazars in the northern and southern hemispheres, respectively. Blazars show flux and spectral variabilities in the complete electromagnetic spectrum on diverse timescales ranging from a few minutes to hours, days, weeks, months and even several years. The flux and spectral variability on different timescales can be used to constrain the size of the emitting region, estimate the super massive black hole mass, find the dominant emission mechanism in the close vicinity of the super massive black hole, search for quasi-periodic oscillations in time series data and~several other physical parameters of blazars. Flux and spectral variability is also a dominant tool to explain jet as well as disk emission from blazars at different epochs of observations., Comment: 15 pages, Published in the special issue "X-Ray Flux and Spectral Variability of Blazars" of Galaxies Journal

We present the results of a temporal and spectral study of the BL Lacertae object OJ 287 in optical, UV, and X-ray bands with observations performed by Swift satellite during September 2019--March 2020. In this period, the source showed moderate variability characterized by variability amplitude of $\sim$22--31% in all the wavelengths on a short timescale, except the hard X-ray band which was variable by only $\sim$8%. We observed that the X-ray flux of the source was significantly dominated by the soft photons below $2$ keV. Soft lags of $\sim$45 days were detected between the optical/UV and soft X-ray emissions, while there is no correlation between the hard X-rays and the lower energy bands indicating the presence of two emission components or electron populations. Although two components contribute to the X-ray emission, most of the 0.3--10 keV spectra were well fitted with an absorbed power-law model which outlines the dominance of synchrotron over inverse Compton (IC) mechanism. The X-ray spectra follow a weak "softer when brighter" trend., Comment: 11 pages, 4 figures, published in the special issue "X-Ray Flux and Spectral Variability of Blazars" of Galaxies journal

We present the results of optical photometric observations of three extreme TeV blazars, 1ES 0229$+$200, 1ES 0414$+$009, and 1ES 2344$+$514, taken with two telescopes (1.3 m Devasthal Fast Optical Telescope, and 1.04 m Sampuranand Telescope) in India and two (1.4 m Milankovi\'{c} telescope and 60 cm Nedeljkovi\'{c} telescope) in Serbia during 2013--2019. We investigated their flux and spectral variability on diverse timescales. We examined a total of 36 intraday $R-$band light curves of these blazars for flux variations using the power-enhanced {\it F}-test and the nested ANOVA test. No significant intraday variation was detected on 35 nights, and during the one positive detection the amplitude of variability was only 2.26 per cent. On yearly timescales, all three blazars showed clear flux variations in all optical wavebands. The weighted mean optical spectral index ($\alpha_{BR}$), calculated using $B - R$ color indices, for 1ES 0229$+$200 was 2.09 $\pm$ 0.01. We also estimated the weighted mean optical spectral indices of 0.67 $\pm$ 0.01 and 1.37 $\pm$ 0.01 for 1ES 0414$+$009, and 1ES 2344$+$514, respectively, by fitting a single power-law ($F_{\nu} \propto \nu^{-\alpha}$) in their optical ({\it VRI}) spectral energy distributions. A bluer-when-brighter trend was only detected in the blazar 1ES 0414$+$009. We briefly discuss different possible physical mechanisms responsible for the observed flux and spectral changes in these blazars on diverse timescales., Comment: 17 pages, 12 figures, 8 tables, accepted for publication in MNRAS journal

We report the detection ($>4\sigma$) of a Quasi-Periodic Oscillation (QPO) in the $\gamma$-ray light curve of 3C 454.3 along with a simultaneous marginal QPO detection ($>2.4\sigma$) in the optical light curves. Periodic flux modulations were detected in both of these wavebands with a dominant period of $\sim 47$ days. The $\gamma$-ray QPO lasted for over 450 days (from MJD 56800 to 57250) resulting in over nine observed cycles which is among the highest number of periods ever detected in a blazar light curve. The optical light curve was not well sampled for almost half of the $\gamma$-ray QPO span due to the daytime transit of the source, which could explain the lower significance of the optical QPO. Autoregressive Integrated Moving Average (ARIMA) modelling of the light curve revealed a significant, exponentially decaying, trend in the light curve during the QPO, along with the $47$ days periodicity. We explore several physical models to explain the origin of this transient quasi-periodic modulation and the overall trend in the observed flux with a month-like period. These scenarios include a binary black hole system, a hotspot orbiting close to the innermost stable circular orbit of the supermassive black hole, and precessing jets. We conclude that the most likely scenario involves a region of enhanced emission moving helically inside a curved jet. The helical motion gives rise to the QPO and the curvature ($\sim 0.05^{\circ}$ pc$^{-1}$) of the jet is responsible for the observed trend in the light curve., Comment: 13 pages, 10 figures, 1 table

We report the results of our optical (VRI) photometric observations of the TeV blazar 1ES 0806$+$524 on 153 nights during 2011-2019 using seven optical telescopes in Europe and Asia. We investigated the variability of the blazar on intraday as well as on long-term timescales. We examined eighteen intraday light curves for flux and color variations using the most reliable power-enhanced F-test and the nested ANOVA test. Only on one night was a small, but significant, variation found, in both $V$ band and $R$ band light curves. The $V-R$ color index was constant on every one of those nights. Flux density changes of around 80 % were seen over the course of these eight years in multiple bands. We found a weighted mean optical spectral index of 0.639$\pm$0.002 during our monitoring period by fitting a power law ($F_{\nu} \propto \nu^{-\alpha}$) in 23 optical ($VRI$) spectral energy distributions of 1ES 0806$+$524. We discuss different possible mechanisms responsible for blazar variability on diverse timescales., Comment: 14 pages, 6 figures, 7 tables, accepted for publication in ApJ

In this paper we report the current status of a new research program. The primary goal of the "Vanishing & Appearing Sources during a Century of Observations" (VASCO) project is to search for vanishing and appearing sources using existing survey data to find examples of exceptional astrophysical transients. The implications of finding such objects extend from traditional astrophysics fields to the more exotic searches for evidence of technologically advanced civilizations. In this first paper we present new, deeper observations of the tentative candidate discovered by Villarroel et al. (2016). We then perform the first searches for vanishing objects throughout the sky by comparing 600 million objects from the US Naval Observatory Catalogue (USNO) B1.0 down to a limiting magnitude of $\sim 20 - 21$ with the recent Pan-STARRS Data Release-1 (DR1) with a limiting magnitude of $\sim$ 23.4. We find about 150,000 preliminary candidates that do not have any Pan-STARRS counterpart within a 30 arcsec radius. We show that these objects are redder and have larger proper motions than typical USNO objects. We visually examine the images for a subset of about 24,000 candidates, superseding the 2016 study with a sample ten times larger. We find about $\sim$ 100 point sources visible in only one epoch in the red band of the USNO which may be of interest in searches for strong M dwarf flares, high-redshift supernovae or other catagories of unidentified red transients., Comment: SETI meets time domain astronomy. Accepted into the Astronomical Journal

Characterisation of the long-term variations in the broad line region in a luminous blazar, where Comptonisation of broad-line emission within a relativistic jet is the standard scenario for production of gamma-ray emission that dominates the spectral energy distribution. We analysed ten years of optical spectroscopic data from the Steward Observatory for the blazar 3C 454.3, as well as gamma-ray data from the Fermi Large Area Telescope (LAT). The optical spectra are dominated by a highly variable non-thermal synchrotron continuum with a prominent Mg II broad emission line. The line flux was obtained by spectral decomposition including significant contribution from the Fe II pseudo-continuum. Three methods were used to characterise variations in the line flux: (1) stacking of the continuum-subtracted spectra, (2) subtracting the running mean light curves calculated for different timescales, and (3) evaluating potential time delays via the discrete correlation function (DCF). Despite very large variations in the gamma-ray and optical continua, the line flux changes only moderately (< 0.1 dex). The data suggest that the line flux responds to a dramatic change in the blazar activity from a very high state in 2010 to a deep low state in 2012. Two interpretations are possible: either the line flux is anti-correlated with the continuum or the increase in the line luminosity is delayed by ~600 days. If this time delay results from the reverberation of poorly constrained accretion disc emission in both the broad-line region (BLR) and the synchrotron emitting blazar zone within a relativistic jet, we would obtain natural estimates for the BLR radius [R_{BLR,MgII} >~ 0.28 pc] and for the supermassive black hole mass [M_SMBH ~ 8.5x10^8 M_sun]. We did not identify additional examples of short-term 'flares' of the line flux, in addition to the previously reported case observed in 2010., Comment: 8 pages, 7 figures, published in A&A, replaced to match the title