arXiv:2105.04028v1 [astro-ph.HE] 9 May 2021 aao J27wsfis eetdby detected first was 287 OJ of data atr ta.2005 al. et Raiteri 1995 Witzel & Wagner al. et Valtonen eo.Cagsoe ag fmntst esta a (e.g. day a than less to minutes of range a di over three Changes into below. fits blazars in ability di in even 287 source OJ variable of active ( variability exceptionally ago an decades as five known was It nrrdfeunywssuiduigsadr H photomet JHK standard using studied was frequency infrared e.g. (LTV, variations are variations term years term to long months short as several d fined as over of known variations the timescale commonly while are a (STV); months on few those a microvariability); to va or intra-day IDV, or i.e. (INV, variability ity intra-night as defined are J27 LLctp cieglci ulu,lctda a at located nucleus, ( 1967 galactic in active discovered type was 0.306, Lac of BL redshift a 287, OJ Introduction 1. a 1 2021 11, May Astronomy ⋆ -al [email protected] E-mail: hr iesaevraiiyo hsB a beti h near the in object Lac BL this of variability scale time Short pia,Xryadgmarydt fO 8 o h eido 2 of period the for 287 OJ of data gamma-ray and X-ray optical, iie nfiesget,rfre sA ,C ,&Ei hspape this in E & D, C, B, varia A, observed as referred the segments, on five in Based divided period. this during frequencies Results. Methods. rar very are systems binary hole Aims. black the jets, scale parsec Context. itiuin n hshlsu oepoetentr fthis zo of words. single nature Key A the state. explore to high us very helps in this was and source distributions sho the has when source 2020 The mid observed. and been has variability flux though wa time variability fast that with along and states, various 6 a 1 2021 11, May h aueo t aigactivities. flaring its of nature the 1 2 aeeghSDmdln spromdt nwmr bu h p the about more know variability. fast to the performed and emission is modeling SED wavelength 3 4 5 Argentina nttt eAtoíiad aPaa(C aPlata-CONICET- La (CCT Plata La de Astrofísica de Instituto etrfrTertclPyis oihAaeyo Sciences of Academy Polish Physics, Theoretical for Center aa eerhIsiue aahvngr aglr 5600 Bangalore Sadashivanagar, Institute, Research Raman aaIsiueo ula hsc,HN,Klaa etBeng West Kolkata, HBNI, Physics, Nuclear of Institute Saha autdd inisAtoóia efscs Universid Geofísicas, y Astronómicas Ciencias de Facultad opeoAtoóio“lLoct”(ALO,CONICET-UNL (CASLEO), Leoncito” “El Astronómico Complejo ut-aeeghAayi n oeigo J27During 287 OJ of Modeling and Analysis Multi-wavelength h LLcO 8 sa neetn betfrmulti-waveleng for object interesting an is 287 OJ Lac BL The & h lzrO 8 a enpooe sbnr lc oesyste hole black binary as proposed been has 287 OJ blazar The h em A bevtosso oeaeflxlvlo hss this of level flux moderate a show observations LAT Fermi The ( hav we source this in variability temporal the understand To Astrophysics 2006 ; aais cie am as aais niiul:O 28 OJ individuals: galaxies; rays: gamma active; galaxies: a Prince Raj nrwe l 1971 al. et Andrew gra ta.2017 al. et Agarwal .Itaa aiblt nrdoadoptical and radio in variability Intraday ). ; imn1975 Kinman ff aucitn.aanda no. manuscript rn iesae a efudin found be can scales time erent 1 ⋆ dt Agarwal Aditi , ff rn aeoisa described as categories erent atoae al. et Valtaoja ; .Apoe td fthe of study proper A ). ). etr&Prmn2003 Perlman & Rector nrco ta.2011 al. et Andruchow egoA Cellone A. Sergio ikle l 1967 al. et Dickel 2 ( aatr Gupta Nayantara , 1985 2017-2020 .Vari- ). riabil- loetmtdt nesadtentr fvraiiy Fu variability. of nature the understand to estimated also s LLcwt iaysprmsiebakholes. black super-massive binary with Lac BL n ec hssuc svr neetn ostudy. to interesting very is source this hence and e ABSTRACT iiyi pia n -a rqece h niepro 20 period entire the frequencies X-ray and optical in bility 0 India 80, ays lLtio 32 Al.Lotnikow , de- nasrn u aiblt nXry pia,adU uige during UV and optical, X-ray, in variability flux strong a wn ry, l 004 India 700064, al, dNcoa eL lt,PsodlBsu,B90W,L Plata La B1900FWA, Bosque, del Paseo Plata, La de Nacional ad ). .Adtie eprladseta nlssi efre to performed is analysis spectral and temporal detailed A r. NP,L lt,Argentina Plata, La UNLP), ) - ; 4 1 00 hr r eea ihsae notclU n X- and optical-UV in states high several are There 2020. – 017 , eSCeiso oe scniee omdlteseta ene spectral the model to considered is model emission SSC ne 5 n eido 3mnh ( obser months the 23 over of mag period 0.7 ing amplitude of variability showed which h rdcino 2ya yl a confirmed. of was cycle year model 12 of hole prediction the black binary the to according th and hole. hole black primary black the secondary of fla the disc pericen optical accretion between during model, accretion collision another enhanced enhancement or from In passage the results jet. 287 of the OJ result of in the ing factor be Doppler precession thus the pr to would in the due manner of flare regular jet optical a the The in One of changes hole. orientation hole the black black that supermassive mary assumes binary model a of kind of f motion optical periodic ( in earlier 287 the OJ discussed of been outburst have periodic quency the for models ferent ( years. ( 11 others of by interval verified the at were in flux the results has minimum at These and it outbursts years repeated that 11.65 shows of inferred curve val was light the ( it that holes 287, found OJ black of supermassive curve binary light optical term .Andruchow I. , 1988 -N-NJ a un Argentina Juan, San P-UNC-UNSJ, yia rcse epnil o h iutnosbroadb simultaneous the for responsible processes hysical hsuydet t eidcotuss ehv nlsdthe analysed have We outbursts. periodic its to due study th 7 i aefo J27wspeitdt apni 1994 in happen to predicted was 287 OJ from flare big A tde h nr-a,adfatoa aiblt o all for variability fractional and intra-day, the studied e ae nisproi pia ubrt mn lzr with blazars Among outburst. optical periodic its on based m 2 rtkMajumdar Pratik , .Ti a bevdby observed was This ). uc ngmaryfeunytruhu hsperiod, this throughout frequency gamma-ray in ource / 6 asw Poland Warsaw, 46, 4 , 6 3 oeztie l 1989 al. et Lorenzetti Bo˙zena Czerny , ilnä tal. et Sillanpää ilnä ta.1988 al. et Sillanpää e ta.2019 al. et Dey ril ubr ae1o 14 of 1 page number, Article igre l 1992 al. et Kidger te,temulti- the rther, 1 et Valtonen & Lehto ( 7–22 is 2020 – 17 .Fo h long- the From ). 1996a understand ilnä tal. et Sillanpää ry2017 arly © ,involving ), S 2021 ESO n thus and ) and rgy ray the .They ). .Dif- ). , ter- re- ter v- r- i- e . A&A proofs: manuscript no. aanda
(1996) proposed that the reason for the flares is an impact of the model. The first zone gives an LBL SED, and the second zone secondary black hole on the accretion disk of the primary, which gives an HBL SED. In their model, the second zone is located at means that there has to be two such flares during each orbital cy- a distance of parsec-scale from the central engine. cle. This is a unique property of the model, not easily accounted A hadronic model to explain the X-ray and gamma-ray out- for in other proposals. The model predicted the time of the sec- burst of November 2015 outburst of OJ 287 has been given by ond flare in November 1995 within a two week time window Rodríguez-Ramírez et al. (2020). They have used the binary su- (Valtonen 1996), and subsequently Sillanpää et al. (1996b) ob- permassive black hole model, where the initial trigger comes served the flare and confirmed the prediction. Sundelius et al. from the impact of the secondary black hole on the accretion (1996, 1997) calculated the flare arising from tides in this bi- disc of the primary black hole. An idealized spherical outflow is nary model and predicted the next big impact flare in 2005, a generated from this impact. A shock is formed when this spher- year earlier than would be expected from strict periodicity. It ical outflow – containing cosmic rays and thermal ions – in- was reported by Valtonen et al. (2006). The flares come sooner teracts with the AGN wind of the primary black hole. In their than in the strictly periodic models due to precession, as is well model, the cosmic rays are shock accelerated due to the colli- stated in their paper. Finally, the observation of the 2015 flare sion of the outflow with the AGN wind of the primary black confirmed this shift, which by then was 3 years (Valtonen et al. hole. The cosmic ray protons interact with the thermal ions, and 2016). This paper also found the signature of disk impacts, the as a result, secondary leptons, photons are produced in proton- thermal nature of the flare. Valtonen et al. (2019) updated the proton interactions. The optical flare is explained by combining model of Lehto & Valtonen (1996) and determined the disk pa- the jet emission from Kushwaha et al. (2013) and the thermal rameters using time delays calculated in Dey et al. (2018). bremsstrahlung emission in the outflow. The photon field pro- During the phase 2008 – 2010 tidal flares were expected duced as a result of thermal bremsstrahlung acts as target for according to the model by Sundelius et al. (1996, 1997). The inverse Compton emission by the secondary leptons. They have gamma-ray light curve of OJ 287 during 2008 August – 2010 explained the X-ray and gamma-ray data by this inverse Comp- January was studied by Neronov & Vovk (2011). They found the ton emission of the secondary electrons. variability time scale is lower than 3.2 hours. They inferred that Recently, Komossa et al. (2020b) reported the detection of a the observedgamma-rayemission was from the jet of the smaller very bright outburst of OJ 287 covering X-ray, UV, and optical mass black hole. Detection of gamma-rays of energy higher than frequency from April to June of 2020. They concluded that the 10 GeV constrained the lower limit of the Doppler factor to 4. outburst is jet-driven and consistent with the binary supermas- The broadband spectrum of the major gamma-ray flare in sive black hole model. In this model, the impact of the secondary 2009 was studied by Kushwaha et al. (2013). They explained the black hole on the disk of the primary triggers an after-flare. This multi-wavelength spectral energy distribution (SED) by combin- impact enhances the accretion activity of the primary black hole, ing synchrotron, synchrotron self-Compton (SSC), and external which results in enhancedjet emission by the primary black hole. Compton (EC) processes. They suggested that the emission re- In this paper, we have analyzed the multi-wavelength data of gion in the jet is surroundedby a bath of photonsat 250 K. They OJ 287 for the period of 2017 to 2020, which includes the out- also inferred that the location of this emission region is 9 pc burst discussed in the paper by Komossa et al. (2020b). The total away from the central engine. The high activity of OJ 287 dur- period 2017 – 2020 considered in our work has been divided into ing December2015 – April 2016 was studied by Kushwaha et al. five segments after analyzing the variability time scale in optical (2018a), and the authors inferred simultaneous multi-wavelength and X-ray data. We have done the modeling of the SEDs with emission. They explained the optical bump as accretion disc a time-dependent leptonic model, which includes synchrotron, emission associated with the primary black hole. The smaller SSC. The data analysis is discussed in section 2. Our results of bump feature in optical-UV appeared to be consistent with line data analysis and the modelling of SEDs are discussed in sec- emission. They explained the gamma-ray emission with inverse tion 3. The discussions and conclusions of our study are given in Compton scattering of photons from the line emission. section 4. The flux and polarisation variability at optical bands of OJ 287 during the December 2015 to February 2016 outburst was studied by Rakshit et al. (2017). The intra-night optical variabil- 2. Multiwavelength Observations and Data Analysis ity data was analyzed, and the shortest variability time scale was estimated as 142 ± 38 minutes. This constrained the lower limit 2.1. Fermi-LAT on the value of the Doppler factor to 1.17 and the upper limit Fermi-LAT is an excellent space-based telescope to explore the on the value of the magnetic field to 3.8 Gauss. The size of the extragalactic and Galactic objects in the gamma-ray sky. It uses 14 emission region was constrained to less than 2.28 × 10 cm. the pair conversion method to detect gamma-rays in the energy The multi-band optical variability from September 2015 range of 20 MeV – 500 GeV. It has a wide field of view (FoV) to May 2016 was studied by Gupta et al. (2016) using nine of about 2.4 sr (Atwood et al. 2009), which scans 20% of the ground-based optical telescopes. They detected a large op- sky at any time. The total scanning period of the entire sky with tical outburst in December 2015 and a second comparably this telescope is around three hours. OJ 287 was observed in the strong flare in March 2016. The long term optical, ultravi- brightest flaring state in X-ray when monitored by the Swift- olet and X-ray variability in different activity states of OJ telescope (Atel 10043) in 2017, and, soon, flares in other fre- 287 was studied using UVOT and XRT instruments of Swift quency bands were also detected. Fermi-LAT is continuously (Siejkowski & Wierzcholska 2017). They did not find any clear monitoring the source OJ 287 since 2008. We have collected the relation between optical/UV and X-ray emission during quies- data from January 2017 to May 2020, and it is found that the cence state or outbursts. source is in a moderate flux state within this period. We have The strong activity in optical to X-ray frequency during July analyzed the gamma-ray data following the standard data reduc- 2016 to July 2017 was studied by Kushwaha et al. (2018b). The tion and analysis procedure described by Science Tools1. The de- daily gamma-ray fluxes during this time are consistent with no variability. They modeled the SEDs with a two-zone leptonic 1 https://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/
Article number, page 2 of 14 Prince et al.: Temporal and Spectral study of OJ 287 tails of the method of this analysis are discussed in Prince et al. Table 1. Log of photometric observations for the blazar OJ 287. (2018).
2.2. X-ray Observations Dateof Telescope Number of data points On February 3, 2017, an X-ray flare was observed by the Swift- observations telescope, and the results were reported in Atel 10043. It is re- (yyyy mm dd) BVRI ported as the brightest flare ever detected since the monitoring 20190215 C 01 6 1 started by the Swift-telescope. After that, many multiple flares 20190227 A 3 2 30 2 were observed in X-rays until May 2020, and this whole period 20190301 A 2 9 27 1 has been studied in this paper. Swift is a space-based telescope 20190302 A 0 9 22 1 with three instruments onboard, observing all kinds of Galactic 20190303 A 1 17 16 2 and extragalactic sources in soft & hard X-rays, Optical, andUV 20190309 C 11 1 1 simultaneously. The working energy range of Swift-XRT is 0.3- 20190311 C 11 1 1 10.0 keV. The BL Lac OJ 287 was observed by Swift-XRT tele- 20190312 C 01 1 0 scope during the multiple flaring episodes in X-ray frequencies 20190313 C 11 1 1 between the period January 2017 to May 2020. We have ana- 20190326 A 2 1 24 2 lyzed all the observations done during this period, and the pro- 20190405 C 1 1 17 1 2 cessing of the raw data is done by using the task ‘xrtpipeline ’ 20190406 C 1 1 10 1 and cleaned events file are produced for each observation. The 20190407 C 1 2 20 2 CALDB version 20160609 is used while processing the raw 20190409 A 1 16 16 1 data. Our analysis is focused on only the Photon Counting mode 20190410 A 2 10 10 2 observations, and the task ‘xselect’ is used for source and back- 20190410 A 3 4 12 2 ground selection. We have selected a region of 12 arc seconds 20190411 A 2 12 14 2 around the source and away from the source for the source and 20191217 A 01 9 2 background, respectively, in our data analysis. The task ‘xse- 20200103 A 2 2 120 2 lect’ is also used to extract the spectrum and light curve, and the 20200127 B 04 4 1 modeling of the spectrum is done in ‘Xspec’(Arnaud 1996). For modeling the spectra, we have used a single power-law model. 20 −2 The Galactic absorption column density nH = 1.10×10 cm is fixed from Kalberla et al. (2005). The modeling is done for an The preliminary data reduction includes bias correction, energy range of 0.3 – 10.0 keV. flat fielding, and cosmic-ray removal, which was performed with IRAF3 software. We then processed the cleaned CCD im- ages using the Dominion Astronomical Observatory Photometry 2.3. Optical and UV Observations (DAOPHOT II) software (Stetson 1987, 1992) using the aper- Having the Swift Ultraviolet/Optical Telescope (UVOT, ture photometry technique through which we obtained instru- Roming et al. 2005) on board with Swift-XRT has the advantage mental magnitudes for our target and four standard stars located of getting simultaneous observations in Optical and UV bands. in the same field. A more detailed and comprehensive descrip- Swift-UVOT has also observed the OJ 287 in all of the available tion of data reduction methods used is given in Section 2 of six filters U, V, B, W1, M2, and W2, simultaneously with Agarwal et al. (2019). Finally, to extract the instrumental dif- the X-ray observations. The source instrumental magnitudes ferential light curves (LCs), we selected two non-variable stan- are extracted following the uvotsource procedure. We have dards having magnitude and color very similar to that of the considered the region of 5 arcsec around the source and away blazar. The calibrated LCs were obtained using the star 10 of from it as the source and the background region, respec- Fiorucci & Tosti (1996). After constructing the calibrated LCs tively, in our data analysis. The magnitudes are corrected for of our source, we carefully inspected the LCs for any outliers. galactic extinction by using the reddening E(B-V) = 0.0241 A handful of such suspicious data points were detected and cor- from Schlafly & Finkbeiner (2011) and zero points from rected. Breeveld et al. (2011). Moreover, the magnitudes are converted into flux by multiplying by the conversion factor estimated by Poole et al. (2008) and the ratios of extinction to reddening from 2.4. Radio data at 15 GHz Giommi et al. (2006). Owens Valley Radio Observatory(OVRO; Richards et al. 2011) In the period between Feb 2019 - Jan 2020, we performed is one of the observatories that monitors the bright Fermi de- observations of OJ 287 using five different telescopes around tected blazars. It is a 40-meter single-dish antenna working at the globe, which are: 2.15m Jorge Sahade telescope (JS, tele- a frequency of 15 GHz. A large number of Fermi blazars are scope A) and 60 cm Helen Sawyer Hogg telescope (HSH, tele- continuously monitored by OVRO twice a week. Our candidate scope B), CASLEO, Argentina; 1.3m JC Bhattacharya telescope source, OJ 287, is also part of the OVRO monitoring program, (JCBT; telescope C) at the Vainu Bappu Observatory (VBO), In- and we have collected the data from September 2017 to July dia. The technical descriptions of the above telescopes are sum- 2020. marized in Table 1 of Agarwal et al. (2019) and Agarwal et al. (2021). The number of observations made in each band on a particular date during our monitoring campaign is provided in 3 IRAF is distributed by the National Optical Astronomy Observato- Table 1. ries, which are operated by the Association of Universities for Research in Astronomy Inc., under a cooperative agreement with the National 2 https://heasarc.gsfc.nasa.gov/ftools/caldb/help/xrtpipeline.html Science Foundation.
Article number, page 3 of 14 A&A proofs: manuscript no. aanda