Issue |
A&A
Volume 688, August 2024
|
|
---|---|---|
Article Number | L16 | |
Number of page(s) | 5 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/202449409 | |
Published online | 12 August 2024 |
Letter to the Editor
Clues on the nature of the quasi-periodic optical outbursts of the blazar OJ 287
UM-DAE Centre for Excellence in Basic Sciences (CEBS), Vidyanagari, Mumbai 400098, India
e-mail: gopaltani@gmail.com
Received:
30
January
2024
Accepted:
4
July
2024
The fascination with the blazar OJ 287 stems not only from its status as a prominent candidate for a close supermassive black hole (SMBH) binary, but also because of the thermal bremsstrahlung origin proposed for its giant optical outbursts. These outbursts arrive as pairs, quasi-periodically every ∼12 years, based on the unique 130-year-long, well-sampled optical light curve available for this blazar. For its three well-known, large quasi-periodic optical outbursts (QPOOs), observed in 1983, 2007 and 2015, optical photo-polarimetric monitoring has been reported in the literature. For these initially radio-undetected QPOOs, widely acclaimed as ‘bremsstrahlung flares’, we have scrutinised the available measurements of optical polarisation and spectral index during the rising phase. Several inconsistencies of these data with the optical bremsstrahlung interpretation are noted, which point towards a synchrotron-dominated alternative interpretation for all these prominent QPOOs, just as for the optical emission observed between the outbursts. Possible reasons for the radio non-detection of the QPOOs during the initial stage are outlined.
Key words: galaxies: active / BL Lacertae objects: individual: OJ 287 / galaxies: jets / galaxies: nuclei / quasars: general / quasars: supermassive black holes
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Even among the tiny sub-population of active galactic nuclei (AGNs) known as blazars, exhibiting strong variability, OJ 287 (J0854+2006, PG0851+202, z = 0.306) stands out for its highly intriguing optical variability patterns. Its exceptionally long recorded optical light curve, going back to 1888, is dotted with outbursts which occur with a quasi-periodicity of about 12 years (Sillanpaa et al. 1988; Hudec et al. 2013). These outbursts are double-peaked with a peak separation of order of one year (Sillanpaa et al. 1996a). Between such quasi-periodic optical outbursts (QPOOs), the observed optical emission, as in the case of all known blazars, is predominantly of synchrotron origin, as inferred from its generally high polarisation levels (e.g. D’arcangelo et al. 2009; Villforth et al. 2009; Rakshit et al. 2017; Singh 2022; Gupta et al. 2023), strong intranight optical variability (INOV, e.g. Sillanpaa et al. 1991; Sagar et al. 2004; Villforth et al. 2009; Goyal et al. 2013; Rakshit et al. 2017; Gupta et al. 2017; Wehrle et al. 2019) and SED modelling (e.g. Abdo et al. 2009; Kushwaha et al. 2018). The latter demonstrated it to be a low-frequency-peaked blazar (Padovani & Giommi 1995). It had been proposed by Sillanpaa et al. (1988) that the QPOOs originate from a close binary of supermassive black holes (SMBHs) undergoing general relativistic precession, such that the pericentric passage of the secondary SMBH triggers a tidally induced mass flow from the disk to the primary SMBH (see, also Valtaoja et al. 2000). Another major difference from typical blazars is that the dominant starting flare of the QPOOs is not detected at radio wavelengths (e.g. Sillanpaa et al. 1996b; Valtaoja et al. 2000). Moreover, it was claimed to be accompanied by a decreasing percentage optical polarisation, P, (see reviews by Takalo 1994; Valtonen et al. 2016; Kushwaha 2022; Valtonen et al. 2023). The concomitant drop in P, together with the radio non-detection, led to the ‘bremsstrahlung flare’ paradigm for the QPOOs (e.g. Sillanpaa 1991; Lehto & Valtonen 1996; Valtaoja et al. 2000) and to the widely explored ‘disk-piercing’ (DP) model of OJ 287, introduced by Lehto & Valtonen (1996), specifically to explain the arrival of the QPOOs in pair, with a year-like separation. In this DP model, the orbiting secondary SMBH, upon impacting the disk of an ultra-massive primary SMBH (∼1010 M⊙), punches a hole in the disk, twice per its highly eccentric, fast-precessing orbit. An optical outburst (‘impact flare’) of bremsstrahlung radiation is supposed to ensue at each piercing of the disk by the secondary SMBH, as the shock-heated disk plasma is pulled out of the impact site and the plasma bubbles blown supersonically on either side of the disk expand and cool (see, also Ivanov et al. 1998; Pihajoki et al. 2013a; Valtonen et al. 2019). Within this framework, the expected gravitational response of the accretion disk to the orbiting secondary SMBH has been calculated by Sundelius et al. (1997) and Pihajoki et al. (2013a) and it constitutes an essential ingredient to predicting the arrival times of the QPOOs (e.g. Valtonen et al. 2023, and references therein). We recall that the DP model posits a geometrically thin (Shakura-Sunyaev) accretion disk; however, this assumption was recently questioned by Antonucci (2023)(see also, Villforth et al. 2010).
Variants of the binary SMBH scenario for explaining the QPOOs of this blazar, some of which associate the QPOOs directly with the primary SMBH, include those proposed by Sillanpaa et al. (1988), Katz (1997), Villata et al. (1998), Valtaoja et al. (2000), Liu & Wu (2002), Qian (2015), Dey et al. (2018), Britzen et al. (2018), and Britzen et al. (2023). In particular, the latter two papers have explicitly argued for a substantial deterministic component in the observed variability of flux density and of the SED state of OJ 287 (see, Kushwaha et al. 2018) by proposing a model of precession and nutation jet motions, partly constrained by the available four-decade-long sequence of very long baseline interferometry (VLBI) images (see also Cohen 2017; Butuzova & Pushkarev 2020). We note that these models belong to a category where the binary is not required to have an ultra-massive black hole (∼1010 M⊙). To verify the predictions of the different models, a dedicated long-term program of dense monitoring of OJ 287 at over a dozen different frequencies from radio to X-rays has been carried out since late 2015: Multi-wavelength Observations and Modelling of OJ 287 (MOMO, Komossa et al. 2021a). The project makes use of a large number of ground-based observatories worldwide and space observatories. This has made OJ 287 one of the best-monitored blazars, particularly in the X-ray/UV/optical regimes (Komossa et al. 2023a,b).
The goal of the present study is basically limited to investigating the nature of the radiation associated with a set of well-observed, prominent QPOOs of OJ 287, which are much-acclaimed cases of ‘bremsstrahlung flares’ (e.g. Valtonen et al. 2021). In Sect. 2, we carried out consistency checks on three such QPOOs, observed in 1983, 2007 and 2015, with the hypothesis of their bremsstrahlung origin, which predicts: (i) a flattening of the optical spectral index (α) due to dilution by the thermal emission, whose spectrum should be nearly flat; (ii) a dip in the optical fractional polarization (P), again due to dilution by the thermal emission which is expected to be unpolarised; and (iii) general stability (i.e. the absence of large fluctuations and/or swings) of the polarisation angle (θ), since θ cannot be altered by addition of the (unpolarised) bremsstrahlung emission (e.g. see Valtaoja et al. 2000). Thus, these are the tell-tale signs of a bremsstrahlung flare, which can be looked for in the existing optical data. Among these, we deem the θ-test to be the most powerful test of the thermal dilution hypothesis. This is because (unlike the first two tests) the θ-test is practically immune with respect to the uncertainty about the level of the underlying synchrotron emission during a putative bremsstrahlung flare, a factor that weakens the other two tests. For instance, this uncertainty virtually precludes the error estimation for the derived spectral index (see e.g. Valtonen et al. 2023). Furthermore, the error gets accentuated due to the small factor of just ∼2 in frequency, which the optical window affords.
The aforementioned three famous QPOOs chosen here for the purposes of the above tests are an unbiased selection, guided solely by the availability of published optical photo-polarimetric measurements covering several consecutive (> 5) nights around the QPOO peak. We recall that a bremsstrahlung flare in the DP model is attributed to collision of the orbiting secondary SMBH with the disk around the (much more massive) primary SMBH. In this process, a bubble of shock-heated plasma is released out of the disk on each side and envelopes the impact site as it expands and cools adiabatically (see above)1. As the expanding plasma bubble becomes transparent, it is posited to emit a burst of ‘bremsstrahlung’ radiation, often capable of overwhelming the underlying synchrotron radiation. After some further cooling, the thermal plasma bubble would evolve into an H II region (e.g. Valtonen et al. 2021). We also note that the QPOOs are found to be sharply peaked initially and compared to the rising phase, which lasts just a couple of days, the decay is slower, probably due to an observed surging contribution from (polarised) optical synchrotron emission, possibly originating near the impact site and becoming increasingly prominent as the opacity of the thermal plasma bubble declines due to its continued expansion (e.g. Valtonen et al. 2021). For this reason, it is much preferable to look for signatures of thermal dilution by a putative ‘bremsstrahlung’ component during the rising part of the QPOO, where this component is supposed to dominate according to the DP model. This is the approach followed below.
2. Analysis of the photo-polarimetric monitoring observations of the three prominent QPOOs
The three QPOOs of OJ 287, with published optical photo-polarimetric measurements on several consecutive nights around the flux maximum, are associated with: January 1983 (Smith et al. 1987), September 2007 (Valtonen et al. 2008), and December 2015 (Valtonen et al. 2016). For the first two QPOOs, the monitoring data of polarisation position angle (θ) have also been reported, in addition to the fractional polarisation (P) and the flux monitoring data. We examine these three datasets below, proceeding in reverse chronological order.
2.1. The outburst of 5 December 2015
This large QPOO began around mid-November 2015, coinciding with the centenary of Einstein’s general relativity theory; hence, its denomination as a ‘GR’ flare (e.g. Valtonen et al. 2016). It was the brightest QPOO observed since the January 1983 QPOO (Sect. 2.3). The outburst consisted of two main flares, altogether spanning ∼40 days. The first (larger) flare peaked on 5 December 2015 (JD 2457361). The second (strongly polarised) flare, peaking ∼18 days later, is clearly a synchrotron flare, thought to originate from the same impact site (e.g. Valtonen et al. 2019). The published data allow for the following two tests to be carried out:
– The α-test: The available multi-colour photometry (Gupta et al. 2017) has scanty pre-flare coverage. Nonetheless, it is seen from their Fig. 2 that the (V–R), (B–V), and (R–I) colour indices at the flare maximum are very close to their long-term values and any flattening of α at the flare maximum is < 0.05. According to these authors, αVR for the low brightness state was ∼ −2.0, while for the outburst (QPOO) state αVR was ∼ −2.13. Thus, there is no evidence of α flattening during the outburst.
– The P-test: Fig. 3 of Valtonen et al. (2016) displays measurements of P (R-band) during this QPOO. For 5 consecutive days up to the ‘bremsstrahlung’ peak (JD 2457361), the plot shows 8 measurements of P, out of which four were made using 2-meter class telescopes and the remaining four with a 0.6-metre telescope (these latter four values are given below in parentheses). These measurements of P, read from their Fig. 3 are, starting from the day of the peak: [6.3%, (8.8)%], −1 day [(5.6%), 7.4%, 7.0%], −2 day [(10.8%)], −3 day [(2.8%)], and −4 day [6.7%]. The corresponding averaged values of P are 7.6% (on the day of the flare peak), 6.7% (-1 day), 10.8% (-2 day), 2.8% (-3 day), and 6.7% (−4 day). We note that during those 5 days, the flux of the flare rose ∼3 times. If we now confine attention to the four measurements made using 2-m class telescopes, the values of P are: 6.3% (at the flare maximum), 7.2% (−1 day), and 6.7% (−4 day). Thus, either way, no systematic decline in P, in sync with the nearly three-fold steady rise in flux, is evident. Unfortunately, the θ-test for this QPOO is precluded due to non-availability of θ monitoring data2.
2.2. The outburst of 13 September 2007
– The P-test: Polarimetric measurements of this QPOO are displayed in Fig. 2 of Valtonen et al. (2008). These cover three consecutive nights preceding the continuum peak (13 September 2007) and three consecutive nights following the peak. The observed two-fold decline in P (starting from P ∼ 8%), concurrently with the doubling of the flux seems consistent with the thermal dilution hypothesis. However, intriguingly, the maximum dip in P, taking it down to ∼2%, occurred two nights before the continuum peak3 (see, Valtonen et al. 2008; Valtonen & Sillanpää 2011);
– The θ-test: Fig. 2 of Valtonen et al. (2008), shows that the only large inter-night spike in θ (by ∼35°) within the 40 days long plot, occurred precisely on the night of the QPOO peak. Clearly, both these polarimetric results are difficult to reconcile with the bremsstrahlung interpretation of this flare.
2.3. The giant outburst of 8 January 1983
For this QPOO, the available data (Smith et al. 1987; Sillanpaa 1991) permit all three tests of the bremsstrahlung hypothesis:
– The α-test: Based on the V-band light curve of the QPOO, displayed in their Fig. 1, Valtonen et al. (2008) set 1983.00 as the starting point of this QPOO. The flare is seen to peak on 8 January, before fading (albeit, dotted with a few smaller synchrotron flares) back to the pre-flare level in ∼45 days. The rising part of the flare amounts to at least a three-fold increase in flux density (Fig. 1), which Valtonen et al. (2008, 2021) identified as a bremsstrahlung event. Such a massive infusion of thermal flux (whose optical spectrum is expected to be nearly flat) should have resulted in a conspicuous spectral flattening in tandem with the flux rise during this giant flare. However, the observed spectral flattening during this week-long rising phase was very marginal (Δα < 0.05, see Fig. 2(b) of Sillanpaa (1991). This is clearly difficult to reconcile with the thermal dilution, namely, the bremsstrahlung hypothesis.
![]() |
Fig. 1. UBVRI photo-polarimetric measurements of the January 1983 QPOO of OJ287. The data are taken from Smith et al. (1987). In the V-band light curve, the first and the last four data points are taken from Sillanpaa (1991). |
– The P-test: The observed drop in P near the peak of this giant QPOO has long been cited as powerful evidence of a bremsstrahlung-dominated rise of the QPOOs of OJ 287 (e.g. Valtonen et al. 2021, however, see more details below.) From Fig. 1, we can see that the light curves in all five optical bands exhibit a general rise in flux density over the five days preceding the peak on 8 January, whereafter a similar decline is seen over the next three days, becoming more gradual over the following week. At the first sight, the observed dipping of P to a minimum exactly on the night of the continuum peak (8 January) does appear to be in accord with the hypothesis of thermal dilution of the existing synchrotron flux, by the flare. However, we argue below that in quantitative terms, this inference is too simplistic.
From Fig. 1, the sharpest drop in P occurred between January 6–8 and we now examine if this precipitous decline in P was indeed due to thermal dilution. From Table 3 of Smith et al. (1987), we see that the R-band flux densities on these two dates were: 30.2±1.4 mJy and 39.2±1.8 mJy (the peak). The corresponding observed values of P were 15.4±0.1% and 0.4±0.1%. This plunge in P, which must be real since it is seen in all the 5 bands (Fig. 1), took place as ∼9 mJy were concurrently added to the flare (R band). Even if all this added flux were bremsstrahlung emission, and making a neutral assumption that the underlying synchrotron component experienced no large change over the two-day interval, the observed increase in flux from ∼30 to ∼39 mJy, would have only decreased P from the observed 15.4 ± 0.1% on 6 January to ∼12% on 8 January, which is 30 times higher than P = 0.4 ± 0.1% actually measured on 8 January. Clearly, most of the observed fall in P did not occur due to thermal dilution and probably has a different origin, possibly a rapid variation in the flux/polarisation of the underlying synchrotron component (e.g. see Holmes et al. 1984). This would not be so surprising for blazars are known to exhibit large (and often correlated) rapid changes in optical continuum and polarisation, thanks to both arising from the same relativistically beamed synchrotron jet. Thus, once again, this QPOO seems to be a synchrotron-dominated event and any role played by thermal dilution is secondary.
– The θ-test: As mentioned in Sect. 1, this is probably the cleanest test of the thermal dilution hypothesis. Thus, although thermal dilution is expected to leave θ unaffected, we notice from Fig. 1 here, along with the V-band data given in Table 3 of Smith et al. (1987), that a large drop in θ occurred on the night of the flare maximum (8 January): from 71.4±0.5 deg on 7 January to 22.3±5.6 deg on 8 January ( 8.7σ change) and then to 133.2±1.0 deg on the very next night (20σ change). Likewise, between 7-9 January, θ changed by 61.8±1.1 deg (55σ). All this leaves little doubt that the giant QPOO of January 1983 was accompanied by huge swings in θ near the peak of the flare. We also note that these swings were seen in all the five bands. This pattern is extremely hard to reconcile with the notion of a bremsstrahlung-dominated flare.
It is also interesting that the large swing in θ within 2 days appears superposed on a fairly steadily rising θ-profile during that week when the QPOO flare stayed above the background (Fig. 1). This indicates presence of two emission components within the flare. Such a two-component structure is also discernible in the flux profile of this QPOO (Fig. 1), wherein a shoulder-like formation can be seen on 6 January, after which the flare steepens again to reach its maximum on 8 January. A similar hint is present on the declining side of the flare (Fig. 1). Thus, conceivably, a short-duration component in the flux profile, lasting for two to three days, which we associate with the simultaneously observed large to-and-fro swing in θ, is superposed on a medium-duration (∼ 10 days) component that may be associated with the above-mentioned steady rise in θ (Fig. 1). Another interesting point to note is that the large abrupt swing in θ over 6-9 January, particularly the sharp dip, is distinctly chromatic, becoming less pronounced on the redder side (Fig. 1, see also Holmes et al. 1984).
3. Summary and conclusions
In this work, we scrutinise the published optical photo-polarimetric monitoring data for three large and prominent QPOOs of OJ 287, selected in an unbiased manner. These three QPOOs have been widely acclaimed in the literature as bremsstrahlung impact flares. We subjected them to a variety of observational tests, based on the premise that the influence of any bremsstrahlung flare should be reflected in the behaviour of three measured optical parameters, namely α, P, and θ. The data availability (Smith et al. 1987; Sillanpaa 1991) enabled us to conduct all three tests for the giant QPOO of January 1983. For the remaining two large QPOOs (September 2007 and December 2015), the available data (Valtonen et al. 2008, 2016; Gupta et al. 2017) only allowed for two of the three tests to be performed (Sect. 2). In not a single one of these seven instances of testing did we find a consonance with the thermal dilution hypothesis, which rests on the bremsstrahlung-dominated interpretation of the QPOOs during the rising phase. Perhaps the most glaring discrepancy pertains to the θ-test, since for both QPOOs (January 1983 and September 2007) for which the published data have permitted this test, a statistically robust huge swing in θ was observed very close to the continuum peak. Likewise, the expected drop in P in sync with the continuum peak is found to be either undetected, or inexplicably large for the thermal dilution scenario, or temporally offset from the continuum peak. Thus, altogether, the present analysis bolsters the case for an alternative mechanism, most probably synchrotron, even for the radio-undetected flares of the QPOOs of OJ 287 (see below).
For the giant QPOO of January 1983, the present analysis has hinted at two emission components within the purported bremsstrahlung flare. These (likely synchrotron) components, lasting ∼10 days and ∼2–3 days, exhibit very different patterns of θ variability (Fig. 1; Sect. 2.3). Holmes et al. (1984) sought to explain the observed large swing in θ near the flux maximum, in terms of two polarised synchrotron sources of comparable luminosities, moving in a relativistic jet, such that their spectral indices differ by ∼0.5 and their polarisation angles are mutually orthogonal, for instance, a shock forming in a stable jet flow. However, that analysis has somehow remained largely ignored in the extensive literature published on this blazar, (see, e.g. Valtonen et al. 2023; Kushwaha 2022).
According to a conceivable, albeit fairly speculative alternative scenario sketched below, the impact of the secondary SMBH on the accretion disk of the primary SMBH triggers jet activity in the secondary SMBH, possibly via Roche-lobe flooding (e.g. Pihajoki et al. 2013a,b; Valtonen et al. 2023). However, the idea of triggering a short-lived powerful jet in the secondary SMBH during its collision with the putative thin disk would require improved theoretical and observational underpinnings, given the extremely violent circumstances expected to prevail during such collisions (see e.g. Komossa et al. 2022). Nonetheless, with this caveat in mind, it is tempting in the context of the January 1983 QPOO, to identify the secondary jet-triggering hypothesis with the above-mentioned first emission component found to last a few days (see Sect. 2.3). The second emission component (of longer duration) could then be an aggregate contribution coming from synchrotron pockets formed near the SMBH-impacted part of the disk which can be strongly magnetised with B > 104 G (Valtonen et al. 2021). Theoretical studies show that the conditions for synchroron emission can arise in the strong shocks of the SMBH-disk collision (see e.g. Medvedev & Loeb 1999; Giacalone & Jokipii 2007). Their synchrotron flux would remain subdued initially because of synchrotron self-absorption and/or free-free absorption, the latter declining as the surrounding impact-driven plasma bubble causing the absorption expands out of the disk (Ivanov et al. 1998; Valtonen et al. 2021). Starting from the actual impact event, the computation of the actual time of onset of synchrotron flare (proposed to be contributed by the magnetised pockets near the impact site, as well as the nascent secondary jet) could be beset with a good deal of uncertainty of astrophysical origin. This problem is compounded for the standard DP model since for predicting the times of appearance of the impact flares (QPOOs), it invokes a bremsstrahlung-dominated post-impact astrophysics, a premise that has now been cast into serious doubt by the present analysis of three famous, well-observed QPOOs that have hitherto been proclaimed to be outstanding cases of ’bremsstrahlung flares’. We note that these cautionary remarks become specially relevant since a sub-day precision has sometimes been estimated or claimed for a predicted time interval between the outbursts observed a decade apart (see e.g. Valtonen et al. 2023; Laine et al. 2020). In general, such an accuracy seems overly optimistic, considering the complex post-impact astrophysics involved, as mentioned above. Furthermore, it may be noted that even when the opacity of the expanding plasma cloud has declined enough to allow escape of the optical photons arising from the embedded synchrotron pockets and the young secondary jet, the opacity may still remain too high for their radio-frequency photons to escape. This could well explain the oft-reported radio non-detection of the initial main flare of the QPOOs (Sect. 1). Another possibility is that the very strong magnetic field near the accretion disk may often shift the bulk of the synchrotron output into the optical regime (see Valtonen et al. 2019).
The work presented here covers the period until the end of 2015 when the QPOO known as the ‘GR flare’ was observed, marking 100 years since the completion of Einstein’s theory of general relativity (see Valtonen et al. 2016; Dey et al. 2018; Valtonen et al. 2019). Then, in late 2015, a new observational phase of this blazar commenced under project MOMO (Sect. 1), resulting in an exceptionally dense long-term multi-band monitoring of this blazar. The largest outburst observed under this project so far had a duration of at least 0.5 yr, with brightness in different bands peaking between October 2016 and March 2017 (Komossa et al. 2020, 2023a). Various observational indicators favour a synchrotron origin for this high-amplitude outburst (Komossa et al. 2023a, and references therein). Coming < 1 yr after the QPOO, called the ‘GR flare’ or the ‘centennial flare’ mentioned above, the 2016/17 outburst was initially perceived as jet-flaring activity unrelated to the SMBH binary or a primary jet-driven ‘after flare’ of the GR impact flare (Valtonen et al. 2017; Komossa et al. 2021b), as envisioned in the binary black hole model of OJ 287, wherein the impact of the secondary black hole on the disk could trigger time-delayed accretion and new jet activity in the primary black hole (Sundelius et al. 1997; Valtonen et al. 2009). However, it was recently suggested that in several properties, the pronounced outburst of 2016/17 matches the previous well known double-peaked impact flares of OJ 287 observed during 1970–2000; thus, it may probably be the latest addition to the long sequence of recorded optical outbursts of this blazar (Komossa et al. 2023a), except that its timing is clearly inconsistent with the DP model prediction. If this surmise is correct, the optical-UV-X-ray outburst of 2016–17, which is also accompanied throughout by radio emission and by γ − ray emission in its rising phase, would be the most recent double-peaked outburst of OJ 287, with even its first peak being synchrotron-dominated. In that case, we note that the first peak is identifiable with a UV flare observed around February 2016 and the second peak with a much larger UV flare seen in October 2016, both having radio counterparts, as per (Fig. 9 Komossa et al. 2023a)4 The next impact outburst of OJ 287, which will be verifiable in the near future, is predicted to begin over 2026–28 (Komossa et al. 2023a). Interestingly, this timing agrees with the prediction coming from a totally different approach, namely the ‘precessing jet’ model of this blazar reported in Britzen et al. (2018).
Reverting to the main theme of the present study, we have shown here that the published optical photo-polarimetric measurements of three very prominent and well-established QPOOs of the blazar OJ 287, selected in an unbiased manner, are inconsistent with their much-acclaimed bremsstrahlung interpretation. Instead, the data distinctly favour the alternative of synchrotron mechanism for all these QPOOs. Possible reasons for the non-detection of the QPOOs at radio wavelengths are outlined. The multi-band search for future QPOO outbursts of OJ 287, expected to occur in 2026–28 (see above), should enable strong constraints to be placed on the parameters of the SMBH binary and on the specific physical mechanism behind the origin of the QPOOs.
Tanaka (2013) has suggested an alternative model for the QPOOs, wherein the secondary SMBH creates a cavity in the disk once every orbit and the disk plasma leaks into it. In a radically different approach which does not involve a SMBH binary, Villforth et al. (2010) have proposed a ‘magnetic breathing’ scenario to explain the double-peaked optical outbursts of OJ 287.
Actually, θ light curve for this QPOO has recently been published by that team (see, Fig. 1 of Gupta et al. (2023). However, individual data points in their long-term θ light curve cannot be visually separated in time, because of the heavy compression of the plot along the time axis. Also, the date of the ‘bremsstrahlung’ flare maximum is given by them as 8 December 2015, which is three days after the date of 5 December 2015 given in Valtonen et al. (2016, Sect. 2.1). The difference is more than the rise time-scale of the flare (∼2 days). Origin of this discrepancy is unclear.
Within the synchrotron jet scenario, a possible explanation for such a temporal offset between flux and P has been provided in the ‘swinging jet’ model (Gopal-Krishna & Wiita 1992).
It may be recalled that for the impact flare of OJ 287, termed ‘Eddington flare’, observed in July 2019 in the near-infrared/optical bands, a thermal spectrum has been inferred, based on relatively sparce data (Laine et al. 2020).
Acknowledgments
I thank Drs. Dusmanta Patra, Krishan Chand and Vibhore Negi for discussions and the Indian National Science Academy for a Senior Scientist position.
References
- Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, ApJ, 700, 597 [CrossRef] [Google Scholar]
- Antonucci, R. R. J. 2023, Galaxies, 11, 102 [NASA ADS] [CrossRef] [Google Scholar]
- Britzen, S., Fendt, C., Witzel, G., et al. 2018, MNRAS, 478, 3199 [NASA ADS] [Google Scholar]
- Britzen, S., Zajaček, M., Gopal-Krishna, et al. 2023, ApJ, 951, 106 [CrossRef] [Google Scholar]
- Butuzova, M. S., & Pushkarev, A. B. 2020, Universe, 6, 191 [NASA ADS] [CrossRef] [Google Scholar]
- Cohen, M. H. 2017, Galaxies, 5, 12 [NASA ADS] [CrossRef] [Google Scholar]
- D’arcangelo, F. D., Marscher, A. P., Jorstad, S. G., et al. 2009, ApJ, 697, 985 [CrossRef] [Google Scholar]
- Dey, L., Valtonen, M. J., Gopakumar, A., et al. 2018, ApJ, 866, 11 [NASA ADS] [CrossRef] [Google Scholar]
- Giacalone, J., & Jokipii, J. R. 2007, AIP Conf. Ser., 932, 243 [NASA ADS] [CrossRef] [Google Scholar]
- Gopal-Krishna, & Wiita, P. J. 1992, A&A, 259, 109 [NASA ADS] [Google Scholar]
- Goyal, A., Gopal-Krishna, Wiita, P. J., Stalin, C. S., & Sagar, R. 2013, MNRAS, 435, 1300 [CrossRef] [Google Scholar]
- Gupta, A. C., Agarwal, A., Mishra, A., et al. 2017, MNRAS, 465, 4423 [NASA ADS] [CrossRef] [Google Scholar]
- Gupta, A. C., Kushwaha, P., Valtonen, M. J., et al. 2023, ApJ, 957, L11 [CrossRef] [Google Scholar]
- Holmes, P. A., Brand, P. W. J. L., Impey, C. D., et al. 1984, MNRAS, 211, 497 [NASA ADS] [CrossRef] [Google Scholar]
- Hudec, R., Bašta, M., Pihajoki, P., & Valtonen, M. 2013, A&A, 559, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ivanov, P. B., Igumenshchev, I. V., & Novikov, I. D. 1998, ApJ, 507, 131 [NASA ADS] [CrossRef] [Google Scholar]
- Katz, J. I. 1997, ApJ, 478, 527 [NASA ADS] [CrossRef] [Google Scholar]
- Komossa, S., Grupe, D., Parker, M. L., et al. 2020, MNRAS, 498, L35 [NASA ADS] [CrossRef] [Google Scholar]
- Komossa, S., Grupe, D., Kraus, A., et al. 2021a, Universe, 7, 261 [NASA ADS] [CrossRef] [Google Scholar]
- Komossa, S., Grupe, D., Parker, M. L., et al. 2021b, MNRAS, 504, 5575 [NASA ADS] [CrossRef] [Google Scholar]
- Komossa, S., Grupe, D., Kraus, A., et al. 2022, MNRAS, 513, 3165 [NASA ADS] [CrossRef] [Google Scholar]
- Komossa, S., Kraus, A., Grupe, D., et al. 2023a, ApJ, 944, 177 [NASA ADS] [CrossRef] [Google Scholar]
- Komossa, S., Kraus, A., Grupe, D., et al. 2023b, Astron. Nachr., 344 [Google Scholar]
- Kushwaha, P. 2022, J. Astrophys. Astron., 43, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Kushwaha, P., Gupta, A. C., Wiita, P. J., et al. 2018, MNRAS, 479, 1672 [NASA ADS] [Google Scholar]
- Laine, S., Dey, L., Valtonen, M., et al. 2020, ApJ, 894, L1 [NASA ADS] [CrossRef] [Google Scholar]
- Lehto, H. J., & Valtonen, M. J. 1996, ApJ, 460, 207 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, F. K., & Wu, X. B. 2002, A&A, 388, L48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Medvedev, M. V., & Loeb, A. 1999, ApJ, 526, 697 [NASA ADS] [CrossRef] [Google Scholar]
- Padovani, P., & Giommi, P. 1995, ApJ, 444, 567 [NASA ADS] [CrossRef] [Google Scholar]
- Pihajoki, P., Valtonen, M., & Ciprini, S. 2013a, MNRAS, 434, 3122 [CrossRef] [Google Scholar]
- Pihajoki, P., Valtonen, M., Zola, S., et al. 2013b, ApJ, 764, 5 [NASA ADS] [CrossRef] [Google Scholar]
- Qian, S.-J. 2015, RAA, 15, 687 [NASA ADS] [Google Scholar]
- Rakshit, S., Stalin, C. S., Muneer, S., Neha, S., & Paliya, V. S. 2017, ApJ, 835, 275 [NASA ADS] [CrossRef] [Google Scholar]
- Sagar, R., Stalin, C. S., Gopal-Krishna, & Wiita, P. J. 2004, MNRAS, 348, 176 [NASA ADS] [CrossRef] [Google Scholar]
- Sillanpaa, A. 1991, Ap&SS, 176, 297 [NASA ADS] [CrossRef] [Google Scholar]
- Sillanpaa, A., Haarala, S., Valtonen, M. J., Sundelius, B., & Byrd, G. G. 1988, ApJ, 325, 628 [NASA ADS] [CrossRef] [Google Scholar]
- Sillanpaa, A., Takalo, L. O., Kikuchi, S., Kidger, M., & de Diego, J. A. 1991, AJ, 101, 2017 [NASA ADS] [CrossRef] [Google Scholar]
- Sillanpaa, A., Takalo, L. O., Pursimo, T., et al. 1996a, A&A, 315, L13 [Google Scholar]
- Sillanpaa, A., Takalo, L. O., Valtaoja, E., et al. 1996b, Mem. Soc. Astron. It., 67, 539 [NASA ADS] [Google Scholar]
- Singh, K. P. 2022, JApA, 43, 85 [Google Scholar]
- Smith, P. S., Balonek, T. J., Elston, R., & Heckert, P. A. 1987, ApJS, 64, 459 [NASA ADS] [CrossRef] [Google Scholar]
- Sundelius, B., Wahde, M., Lehto, H. J., & Valtonen, M. J. 1997, ApJ, 484, 180 [NASA ADS] [CrossRef] [Google Scholar]
- Takalo, L. O. 1994, Vistas Astron., 38, 77 [NASA ADS] [CrossRef] [Google Scholar]
- Tanaka, T. L. 2013, MNRAS, 434, 2275 [NASA ADS] [CrossRef] [Google Scholar]
- Valtaoja, E., Teräsranta, H., Tornikoski, M., et al. 2000, ApJ, 531, 744 [NASA ADS] [CrossRef] [Google Scholar]
- Valtonen, M., & Sillanpää, A. 2011, Acta Polytechn., 51, 76 [NASA ADS] [CrossRef] [Google Scholar]
- Valtonen, M. J., Lehto, H. J., Nilsson, K., et al. 2008, Nature, 452, 851 [Google Scholar]
- Valtonen, M. J., Nilsson, K., Villforth, C., et al. 2009, ApJ, 698, 781 [NASA ADS] [CrossRef] [Google Scholar]
- Valtonen, M. J., Zola, S., Ciprini, S., et al. 2016, ApJ, 819, L37 [NASA ADS] [CrossRef] [Google Scholar]
- Valtonen, M., Zola, S., Jermak, H., et al. 2017, Galaxies, 5, 83 [CrossRef] [Google Scholar]
- Valtonen, M. J., Zola, S., Pihajoki, P., et al. 2019, ApJ, 882, 88 [CrossRef] [Google Scholar]
- Valtonen, M. J., Dey, L., Gopakumar, A., et al. 2021, Galaxies, 10, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Valtonen, M. J., Dey, L., Gopakumar, A., et al. 2023, Galaxies, 11, 82 [NASA ADS] [CrossRef] [Google Scholar]
- Villata, M., Raiteri, C. M., Sillanpaa, A., & Takalo, L. O. 1998, MNRAS, 293, L13 [NASA ADS] [CrossRef] [Google Scholar]
- Villforth, C., Nilsson, K., Østensen, R., et al. 2009, MNRAS, 397, 1893 [CrossRef] [Google Scholar]
- Villforth, C., Nilsson, K., Heidt, J., et al. 2010, MNRAS, 402, 2087 [Google Scholar]
- Wehrle, A. E., Carini, M., & Wiita, P. J. 2019, ApJ, 877, 151 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
![]() |
Fig. 1. UBVRI photo-polarimetric measurements of the January 1983 QPOO of OJ287. The data are taken from Smith et al. (1987). In the V-band light curve, the first and the last four data points are taken from Sillanpaa (1991). |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.