Open Access
Issue
A&A
Volume 694, February 2025
Article Number A291
Number of page(s) 14
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202450003
Published online 20 February 2025

© The Authors 2025

Licence Creative CommonsOpen 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

Blazars are among the brightest sources in the Universe (Angel & Stockman 1980; Hovatta & Lindfors 2019). The preferential orientation of blazar jets is close to the line of sight (for instance, < 10 deg; Readhead et al. 1978; Blandford & Königl 1979; Hovatta et al. 2009; Liodakis et al. 2018a), resulting in rapid variability, relativistically boosted emission, and superluminal jet motions (Blandford et al. 1977; Dondi & Ghisellini 1995; Fuhrmann et al. 2014). Blazars are traditionally classified as flat-spectrum radio quasars (FSRQs) and BL Lac objects (BL Lacs). The FSRQs are high-power blazars with prominent broad emission lines in their optical spectra, whereas BL Lacs are low-power ones with weak or absent optical emission lines (Urry & Padovani 1995; Blandford et al. 2019).

In more than 15 years of observations of the extragalactic sky with the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope (Atwood et al. 2009), blazars have constituted the largest fraction of the γ-ray sources in the extragalactic sky (e.g. Abdo et al. 2010b; Nolan et al. 2012; Acero et al. 2015; Abdollahi et al. 2020). Many blazars exhibit extreme flaring activity in γ-rays, the causes of which are still a matter of debate (e.g. Tavecchio et al. 1998; Ghisellini et al. 2005; Sikora et al. 2009; MacDonald et al. 2015). The γ-ray emission from relativistic jets is generally attributed to inverse Compton (IC) scattering of low-energy photons by synchrotron emitting electrons in the jet (Boettcher et al. 2012; Madejski & Sikora 2016; Hovatta & Lindfors 2019). The target photons may come from synchrotron radiation by the same population of relativistic electrons or from radiation fields external to the jet, for example, the broad-line region (BLR) and dusty torus.

For the former case, the IC scattering is referred to as synchrotron self-Compton (SSC; e.g. Maraschi et al. 1992), and for the latter one, it is called external Compton (EC; e.g. Dermer et al. 1992). Multi-wavelength flux variations in radio, optical, or X-rays are often observed during γ-ray flares. Significant correlations between γ-ray and radio light curves have been reported by many studies (e.g. Pushkarev et al. 2010; Agudo et al. 2011a; Rani et al. 2013; Fuhrmann et al. 2014; Max-Moerbeck et al. 2014a; Liodakis et al. 2018b). These correlations, along with Very Long Baseline Interferometry (VLBI) observations conducted at angular scales of milliarcseconds (mas), suggest that γ-ray flares most likely occur within parsec-scale jets (e.g. Jorstad et al. 2010; Ramakrishnan et al. 2014; Lisakov et al. 2017; Algaba et al. 2018; Kim et al. 2022a).

CTA 102 (z = 1.037) is a γ-ray bright FSRQ, hosting a supermassive black hole (SMBH) of 8.5 × 108 M (Zamaninasab et al. 2014). This source exhibits multi-wavelength variability from radio to γ-rays (e.g. Casadio et al. 2015; Raiteri et al. 2017; D’Ammando et al. 2019). Previous studies investigated the variability of the source in the 15 GHz radio band in 2012–2018 observations, finding variability on time scales of the order of months (Kim et al. 2022b).

The Fermi-LAT observations of CTA 102 have detected γ-ray activity. It has been suggested that the 2012 γ-ray flare occurred a few parsecs from the black hole and coincided with the ejection of a new superluminal jet component from the 43 GHz core (Casadio et al. 2015). Additionally, the source exhibited similar behaviour in the optical and γ-ray bands during the γ-ray flare, which is consistent with the framework of the SSC process (Larionov et al. 2016).

Many studies have investigated γ-ray flares from CTA 102 during 2016–2017 (Zacharias et al. 2017; Li et al. 2018; Shukla et al. 2018; Kaur & Baliyan 2018; Gasparyan et al. 2018; Prince et al. 2018; Casadio et al. 2019; D’Ammando et al. 2019; Sahakyan 2020; Geng et al. 2022; Sahakyan et al. 2022). A γ-ray flare observed in 2016 January was localised to 5.7–16.7 pc from the central engine (Li et al. 2018). From late 2016 to early 2017, the source was in a prolonged period of activity, reaching its highest γ-ray flux in 2016 December, marking one of the brightest γ-ray flares in blazars (Ciprini 2016). Multi-wavelength studies of the 2016 December γ-ray flare found coincident optical and X-ray flares, suggesting a common origin for the emission in these bands (Kaur & Baliyan 2018; Gasparyan et al. 2018; D’Ammando et al. 2019). The multi-wavelength flux variations during the flare have been interpreted as arising due to inhomogeneities in a curved jet that cause different emission regions to change their orientations and hence their Doppler factors (Raiteri et al. 2017; D’Ammando et al. 2019). One scenario suggests that the flare resulted from the ablation of a gas cloud passing through the jet (Zacharias et al. 2017, 2019). Moreover, investigations of the spectral energy distribution (SED) of CTA 102 during the flare demonstrate that the flare can be explained by leptonic processes (Prince et al. 2018; Gasparyan et al. 2018; Sahakyan 2020; Sahakyan et al. 2022). Alternatively, Casadio et al. (2019) explain the flare during 2016–2017 as the interaction between a new jet component and a recollimation shock at ∼0.1 mas from the core. Magnetic reconnection has also been proposed as a possible physical mechanism due to the extreme γ-ray variability of the source and the proximity (around a few parsecs) of the flaring region to the central black hole (Geng et al. 2022).

Additional strong γ-ray flares were observed from CTA 102 in 2018 (e.g. Geng et al. 2022; Sahakyan et al. 2022). Geng et al. (2022) used γ-ray flaring time scales, the bulk Lorentz factor, and the molecular torus luminosity to calculate an upper limit on the distance between the flaring region and the central black hole of 7.22 pc. The γ-ray emission in 2018 is well described by a combination of the SSC and EC models (Sahakyan et al. 2022). However, while several studies have explored the 2018 γ-ray flares, none searched for correlations between γ-ray and radio activity (e.g. Geng et al. 2022; Sahakyan et al. 2022). We note that no such correlations were found for the major γ-ray flares that occurred in 2012–2017 (e.g. Casadio et al. 2015; D’Ammando et al. 2019). For the period of 2018, Kim et al. (2022b) identify radio activity. Our study is the first to report a correlation between γ-ray and radio activity during this period. Moreover, a correlation study between γ-ray and radio light curves during flaring activity is crucial to constrain the locations of the γ-ray regions, which may give us useful information about the nature of γ-ray flares (e.g. Fuhrmann et al. 2014; Max-Moerbeck et al. 2014b; Algaba et al. 2018; Kim et al. 2022a).

This study aims at understanding the physical mechanisms for generating γ-ray flares and the locations of the emission regions in the jet of CTA 102 during the period 2016–2018. We present results from this work, including multi-wavelength correlations with respect to the γ-ray light curve and kinematics of the radio jet related to γ-ray flares. We provide the first results of examining a correlation between γ-ray and radio flares during 2017–2018. In Sect. 2, we describe the observations and data from radio to γ-rays. Multi-wavelength light curves are analysed statistically in Sect. 3. We present the cross-correlation results in Sect. 4. In Sect. 5, we present the parsec-scale jet kinematics of CTA 102 at 43 GHz. We discuss the physical origins of γ-ray flares based on our results in Sect. 6. A summary of our results is given in Sect. 7.

Throughout this paper, we adopt a flat ΛCDM cosmology with Ωm = 0.315, ΩΛ = 0.685, and H0 = 67.4 km s−1 Mpc−1 (Planck Collaboration VI 2020). The luminosity distance at the redshift of the source (i.e. z = 1.037) is DL = 7112.9 Mpc and the linear scale is 8.311 pc mas−1, calculated with the Cosmology Calculator1 (Wright 2006). A proper motion of 1 mas yr−1 corresponds to 55.2c.

2. Observations and data

2.1. Radio data

We used the multi-wavelength radio data from Kim et al. (2022b) as well as observations from the Submillimeter Array (SMA) taken at 340 GHz (for details and references for the data, see Kim et al. 2022b). The multi-wavelength data cover a period from 2012 November 20 to 2018 September 23 (modified Julian date, MJD 56251–58384) with a frequency range of 15–343 GHz. They are complemented by data from the Owens Valley Radio Observatory (OVRO) 40-m radio telescope at 15 GHz (Richards et al. 2011), the Metsähovi 13.7-m radio telescope at 37 GHz (Teraesranta et al. 1998), the IRAM 30-m telescope at 86 and 230 GHz, the SMA sat 230 GHz (Gurwell et al. 2007), and the Atacama Large millimeter/submillimeter Array (ALMA) at 91, 233, and 343 GHz. The 86/230 GHz IRAM data were obtained under the POLAMI (Polarimetric Monitoring of AGN at Millimetre Wavelengths) IRAM Large Program2 (Agudo et al. 2018a,b; Thum et al. 2018). We note that in combining data at similar frequencies, we used the 89 GHz (IRAM + ALMA), 230 GHz (IRAM + SMA + ALMA), and 343 GHz (SMA + ALMA) GHz data for the analyses discussed in this paper. More detailed descriptions of the multi-wavelength data can be found in Kim et al. (2022b).

The VLBA-BU-BLAZAR monitoring program3 has observed γ-ray bright blazars, including CTA 102, with the Very Long Baseline Array (VLBA) at 43 GHz. We used the calibrated VLBA data4 observed from 2016 January 1 to 2018 August 26 (MJD 57388–58356), with a total of 27 epochs. Details of the data calibration are described in Jorstad et al. (2017). The imaging results using the CLEAN task in the interferometric imaging software DIFMAP (Shepherd et al. 1997) are provided on the VLBA-BU-BLAZAR monitoring program website. Further analysis (e.g. Gaussian model fitting) has been performed on the data as described in Sect. 5.

2.2. Optical data

We used data taken in the V and R bands by the optical spectropolarimetric monitoring program of the Steward Observatory at the University of Arizona. This program utilises either the 2.3-m Bok Telescope on Kitt Peak or the 1.54-m Kuiper Telescope on Mt. Bigelow (see Smith et al. 2009, for a more detailed description). The data span a period from 2012 October 11 to 2018 July 7 (MJD 56211–58306). The reduced photometry data are publicly available5. The magnitudes of the optical data have been corrected for Galactic extinction with values in the NASA/IPAC Extragalactic Database (NED)6. The effective wavelengths for the V and R bands are 0.55 and 0.64 μm, respectively (Schlafly & Finkbeiner 2011). We converted the corrected magnitudes into flux densities (in units of mJy) using zero-magnitude flux densities (Mead et al. 1990).

2.3. Hard X-ray data

The Neil Gehrels Swift Observatory (Gehrels et al. 2004) detected the source during 157 months from 2012 September 29 to 2017 December 30 (MJD 56200–58118). The observations were performed with the Burst Alert Telescope (BAT; Barthelmy et al. 2005) and will be included in the upcoming Swift-BAT 157-Month Hard X-ray Survey (Lien et al., in prep.)7. We obtained the Crab-weighted monthly binned light curve at 14–195 keV. We flagged the data when the count rate was below zero.

2.4. γ-ray data

We have analysed the γ-ray data from the Fermi-LAT during the period from 2012 September 29 to 2018 October 8 (MJD 56200–58400), following the standard unbinned likelihood procedure. We performed the analysis with the SCIENCETOOLS8 software package version v11r5p3. We selected the SOURCE class events (evclass = 128, evtype = 3) at 100 MeV–300 GeV. To ensure that we used good quality data, we used the filter parameters DATA_Qual > 0 && LAT_CONFIG = = 1. We applied the maximum zenith angle of 90 deg to reduce contamination from Earth limb γ-rays. We defined a Region of Interest (ROI) of 15 deg radius centred at the location of CTA 102. We used P8R3_SOURCE_V1 as the instrument response function. In order to model the background, we also used the Galactic diffuse emission model gll_iem_v07 and the isotropic background model iso_P8R3_SOURCE_V3_v1.

The spectral analysis was performed using the Science tool gtlike. The source model included all sources from the fourth Fermi Large Area Telescope catalogue (4FGL; Abdollahi et al. 2020) that appeared within 30 deg of the centre of the ROI. For 4FGL sources, we parameterised their spectra using a power law (PL), a log-parabola (LP), or a PL with a super-exponential cut-off (see e.g. Prince et al. 2018, for detailed equations). Model parameters for the sources within 10 deg of the centre of the ROI and for those whose detection significances are ≥15 were set free. For CTA 102, the PL model was used, in the form of dN(E)/dE = N0(E/E0)−Γ, where N0 is the prefactor, Γ is the photon index, and E0 is the pivot energy. We left the prefactor and the photon index free to vary and used the pivot energy E0 = 414 MeV from 4FGL (Abdollahi et al. 2020). The normalisation factors of variable sources (those with variability indexes ≥ 18.48 in the 4FGL catalogue) that were within 10 deg of CTA 102 were also left free to vary. For the other sources, the normalisation factor and the spectral shape parameters were fixed. We evaluated the significance of the γ-ray signals from the sources using a maximum-likelihood analysis and setting a test statistic (TS) value of 10 as the detection limit.

We performed the first maximum-likelihood analysis over the whole period. The second maximum-likelihood analysis was performed after removing the sources with TS < 10. Applying this result to the source model, we generated a weekly binned light curve for the entire period in this work. One of the solar arrays on the Fermi-LAT experienced an anomaly on 2018 March 169. We did not make use of the data from mid-March to early April 2018. The fit results for integrating over the whole period yield a TS value of 340 798 in the 100 MeV–300 GeV energy range. The integrated average flux is (129.4 ± 0.5)×10−8 ph cm−2 s−1, and the photon index is 2.239 ± 0.003.

3. Multi-wavelength light curves

Fig. 1 presents multi-wavelength light curves from CTA 102 in the γ-ray, X-ray, optical, and radio bands. The light curves cover a time range from 2012 September 30 to 2018 October 9 (MJD 56200–58400). We define the γ-ray brightness state (i.e. a quiescent state and an active state) of the source based on the weighted mean flux of ⟨Fω⟩ = 0.7 × 10−6 ph cm−2 s−1 and the weighted standard deviation of σω = 0.9 × 10−6 ph cm−2 s−1 (Williamson et al. 2014). The source is quiescent with flux densities at the level of Fγ < ⟨Fω⟩, while active with those at the level of Fγ > ⟨Fω⟩ + σω. We regard a threshold as ⟨Fω⟩ + 3σω = 3.4 × 10−6 ph cm−2 s−1 to define a flaring state. Until late 2015, the γ-ray brightness of the source was quiescent with flux densities at the level of Fγ ≤ 1.6 × 10−6 ph cm−2 s−1. After 2016, the source became active in γ-rays.

thumbnail Fig. 1.

Multi-wavelength light curves from CTA 102. The entire period is from 2012 September 30 to 2018 October 9 (MJD 56200–58400). From top to bottom, the panels show light curves in γ-rays, X-rays, optical, and radio with different colours and symbols indicating different frequencies. The red triangles in the γ-ray light curve indicate upper limits on the flux density. The optical and radio light curves are displayed with arbitrary offsets for clarity. The grey horizontal lines for the γ-ray light curve denote the flux density threshold (see the text for details). The black dotted lines indicate the peaks of the γ-ray flares (1–6). The grey shaded areas mark three periods (P1–3) for the cross-correlation analysis.

For the period from 2015 October 5 to 2018 October 9 (MJD 57300–58400), we identified six γ-ray flares, labelled (1–6) in Fig. 1. Flare 3 is the brightest γ-ray flare, reaching a peak of (1.415 ± 0.018)×10−5 ph cm−2 s−1 on 2016 December 28 (MJD 57750). D’Ammando et al. (2019) studied this source for the period of 2013 January 1–2017 February 9 (MJD 56293–57793), which overlaps with the period studied here. Three of the flares identified in this work (Flares 1, 2, and 3) were also identified by that group (denoted Flares IV, V, and VI in their work).

We employed a Bayesian block analysis (Scargle et al. 2013) to identify the local maxima of multi-wavelength light curves. This method divides the light curves into blocks so that the flux measurements within each block are statistically consistent with a constant value. In this analysis, we set a false positive rate of 0.01 (the probability of falsely reporting the detection of a change point) for all energy bands. We identify a local maximum if the flux density of the block exceeds those of both previous and subsequent blocks. Fig. 2 presents the Bayesian block representation of the light curves in the γ-ray, X-ray, optical, and radio bands, zooming in on the period that includes the γ-ray flares (see above). Our analysis revealed the following number of local maxima: 24 in the γ-ray light curve, 2 in the X-ray light curve, 31 in the optical V-band light curve, 21 in the optical R-band light curve, and 4–11 in the radio light curves (4 at 343 GHz, 5 at 15/37 GHz, 6 at 89 GHz, and 11 at 230 GHz). We note a larger number of local maxima in the optical and γ-ray light curves compared to the radio band, suggesting more rapid variability at higher energies. The small number of local maxima in X-rays is mainly related to the relatively low sensitivity of Swift-BAT. Kim et al. (2022b) have reported that the rise time of radio flares in the jet of the source is < 250 days.

thumbnail Fig. 2.

Bayesian block representation of the multi-wavelength light curves during the flaring period. The time range of the period is from 2015 October 5 to 2018 October 9 (MJD 57300–58400). The units of the flux densities and vertical lines are the same as in Fig. 1. The horizontal arrows indicate the time ranges of flaring periods (P1–3).

To explore the cross-correlations between γ-ray and lower-energy light curves in Sect. 4, we split the period into three parts, P1, P2, and P3, based on the Bayesian block analysis results. The period P2 is defined to include the strongest γ-ray flare (Flare 3) and its neighboring flares (Flares 2, 4, and 5), as well as the most prominent local maxima in the X-ray and optical light curves. This period extends from MJD 57550 to MJD 57940. The beginning and ending epochs of the period are aligned with the local minima identified in the Bayesian block analysis of the well-sampled 15 GHz light curve. The periods preceding and following P2 within the overall flaring period were determined as P1 (MJD 57300–57550) and P3 (MJD 57940–58400), respectively. Table 1 provides the duration of the periods and the peak information of the γ-ray flares.

Table 1.

γ-ray flares in 2016–2018 (MJD 57300–58400)

We searched for variations in the γ-ray photon index over the flaring period divided into seven-day time bins. Fig. 3 presents the γ-ray light curve in the upper panel and the γ-ray photon index in the lower panel. The black dotted lines (1–6) indicate γ-ray flares that are the same as in Fig. 1. The γ-ray photon index decreases when the source becomes brighter, indicative of spectral hardening (also known as the harder-when-brighter trend). We note that the γ-ray spectrum also becomes harder when small flux enhancements (e.g. between Flares 1 and 2 or between Flares 5 and 6) are observed below the γ-ray flaring threshold.

thumbnail Fig. 3.

γ-ray light curve (top) and time series for the γ-ray photon index (bottom). Black dotted lines (1–6) indicate the peaks of the γ-ray flares.

Fig. 4 displays the γ-ray flux density versus the γ-ray photon index over the flaring period. As in Fig. 3, a clear harder-when-brighter trend is seen in Fig. 4. The γ-ray photon index varies between 1.98 and 2.62 with an average of 2.31 ± 0.01. Consistent results were obtained by D’Ammando et al. (2019) who studied variations in the γ-ray photon index occurring over a time period of 2013 January 1 to 2017 February 9 (MJD 56293–57793) divided into 30-day time bins. Assuming the LP model, they found a spectral slope that varied between 1.88 and 2.97 with an average of 2.30 ± 0.09. We performed a Pearson correlation analysis to evaluate the correlation between the γ-ray flux density and the γ-ray photon index. In quantifying the correlation between photon index and flux density, we calculated the Pearson correlation coefficient, obtaining a value of rp = −0.68 for a confidence level greater than 99.99%. The negative value for rp suggests an anti-correlation. Therefore, we conclude that the harder-when-brighter trend is seen in the γ-ray band (see also Algaba et al. 2018; Rani et al. 2018; Prince 2020).

thumbnail Fig. 4.

γ-ray photon index (Γ) versus γ-ray flux density (Fγ). The Pearson correlation coefficient rp is presented with the confidence level.

4. Correlation analysis

4.1. Method

We carried out a correlation analysis that is applicable to unevenly sampled data sets. For this analysis, we calculated the discrete correlation function (DCF) in order to test for correlations between the various multi-wavelength light curves and to determine their time lags. Details on the DCF are described in Edelson & Krolik (1988). A positive (negative) time lag indicates that the higher-frequency time series leads (lags) the lower-frequency time series. For our calculations, we used the PYTHON implementation of the public DCF code developed by Robertson et al. (2015).

In order to compute the DCF, there are a few effects that we consider. Peterson et al. (2004) found that detrending does not improve the accuracy of estimates of time lags in unevenly sampled time-series data; rather, detrending results in large errors. Thus, we did not detrend our light curves, following Max-Moerbeck et al. (2014b). For each frequency pair, we equated the bin width of the time lag (Δtbin) with the mean cadence (tint) of the less frequently sampled light curve. The optimal bin width was determined in a range from tint to 2tint. For computing the time lags, we selected a range that covers at least two-thirds of the light curves (e.g. Liodakis et al. 2018b).

To determine the significance of the DCF, we simulated artificial light curves with statistical properties, that is, the power spectral density (PSD) and the probability density function (PDF), similar to the real light curves for each frequency pair. For these simulations, we used the public code developed by Connolly (2015) that implements the model described in Emmanoulopoulos et al. (2013). We calculated the PSDs using the normalised Scargle periodogram (Scargle 1982). In order to fit the periodogram (Pf), we used a power-law spectrum, given by P f f β $ P_{\mathrm{f}} \propto {f}^{-\beta} $, where β is the slope of the power-law spectrum (Uttley et al. 2002; Vaughan et al. 2003). We performed Monte Carlo simulations of light curves to find the best-fit slopes of the observed power spectra (see Park & Trippe 2014, 2017, for more details on the method). We also calculated the PDFs of the multi-wavelength light curves, finding that the optical and γ-ray light curves follow gamma distributions, while the radio and X-ray light curves follow log-normal distributions (see Algaba et al. 2018, for more details on the method). We then calculated the DCFs for a large sample of simulated light curves, obtaining 10,000 DCFs, which we used to estimate 1σ, 2σ, and 3σ (68%, 95%, and 99.7%, respectively) confidence intervals.

To measure the uncertainty in a DCF peak (DCFpeak) and its corresponding time lag (τDCF), we used the flux randomisation and random subset selection (FRRSS) method (Peterson et al. 1998). Using this method, we simulated 10 000 artificial light curves randomly sampled from each light curve. Then, we performed the DCF analysis and obtained 10 000 DCFpeak and τDCF values for each frequency pair. We fit the DCFpeak and τDCF distributions with Gaussian functions with means and standard deviations corresponding to the best-fit values for the parameters and their uncertainties, respectively.

4.2. Results

We carried out the cross-correlation analysis over the entire period from 2012 September 30 to 2018 October 9 (MJD 56200–58400) in order to characterise the overall correlation between the γ-ray light curve and those of the other bands. The analysis reveals significant correlations between the optical/X-ray and γ-ray light curves above 99.7% confidence levels, implying a time coincidence and hence co-spatiality of the optical/X-ray and γ-ray emission regions, given in Fig. 5. The 3σ correlations come from the cross-match of Flare 3 with local maxima in the optical and X-ray bands during P2. Cross-correlation analyses during the flaring period (MJD 57300–58400) confirmed a time coincidence between the optical/X-ray and γ-ray light curves. D’Ammando et al. (2019) also reported a correlation with a zero time lag between optical/X-ray and γ-ray light curves, consistent with our results.

thumbnail Fig. 5.

Discrete correlation functions for γ-ray versus optical (R and V bands) and X-ray over the entire period, from top left to bottom right. The red solid and green dashed lines indicate 2σ (95%) and 3σ (99.7%) confidence intervals, respectively. Positive time lags indicate that the flux density in γ-rays leads the flux density in lower-energy bands.

The γ-ray light curve exhibited more complicated correlations with the radio light curves than with the optical band. Fig. 6 presents the DCF results for γ-ray versus radio (15 and 89 GHz) over the entire period. The analysis reveals double 2σ correlations (for 15 and 37 GHz) with positive time lags of a few hundred days (top in Fig. 6), as well as 2σ correlations (for 89, 230, and 343 GHz) with a broad range of positive time lags (bottom in Fig. 6). The former comes from the matching of the 15/37 GHz local maxima in P3 with Flare 1 (or Flare 3) based on the peak-to-peak time intervals (see Fig. 2). The latter occurs because local maxima in the radio (i.e. 89–343 GHz) coincide with the γ-ray flares (e.g. Flare 5 or Flare 6), while others lag behind the γ-ray flares (e.g. Flare 1 or Flare 2).

thumbnail Fig. 6.

Discrete correlation functions for γ-ray versus radio (15 and 89 GHz) over the entire period. Horizontal lines are the same as in Fig. 5.

These correlations are consistent with those observed during the flaring period. As noted in the Bayesian block analysis (Sect. 3), the different variability behaviours between the γ-ray and radio emissions could affect the correlation results. To further explore these correlations, we investigated cross-correlations between the γ-ray and lower-energy light curves during three distinct periods (P1–3).

We found no significant correlations for P1. However, significant correlations (above 99.7% confidence level) were observed between the γ-ray and radio light curves for P2 and P3. Fig. 7 presents the observed negative time lags of <  − 50 days in the 3σ correlation between the γ-ray and 15/37 GHz radio light curves during P2. This correlation is based on matching the 15/37 GHz local maxima in P2 with Flare 3, suggesting that the γ-ray emission lagged behind the radio emission.

thumbnail Fig. 7.

Discrete correlation functions for γ-ray versus radio (15 and 37 GHz) and optical during P2. Horizontal lines are the same as in Fig. 5.

Flare 6, the most recent γ-ray flare in P3, is significantly correlated (> 99.7% confidence levels) with the local maxima in the radio band, yielding time lags of less than a month (see Fig. 8). The positive time lags between the γ-ray and radio light curves during P3 imply that the γ-ray emission preceded the radio emission, including the strongest local maxima at 15–230 GHz frequencies. Table 2 summarises the cross-correlation results for P2 and P3.

thumbnail Fig. 8.

Discrete correlation functions for γ-ray versus radio (15, 37, 89, and 230 GHz) and optical (R and V bands) in P3 from top left to bottom right. Horizontal lines are the same as in Fig. 5.

Table 2.

Cross-correlation results for γ-ray versus multi-wavelength (radio and optical) pairs during P2 and P3 (MJD 57550–58400).

In the γ-radio correlation analysis, we observed both positive and negative time lags, indicating the existence of multiple γ-ray emission regions in the jet. Positive time lags are a common finding in studies of correlations between radio and γ-ray activity in blazar jets (e.g. Max-Moerbeck et al. 2014b; Ramakrishnan et al. 2015; Algaba et al. 2018; Hodgson et al. 2018; Liodakis et al. 2018b). Kim et al. (2022b) examined the synchrotron spectrum of CTA 102 during P3, finding that the synchrotron radiation is self-absorbed with a turnover frequency of around 110 GHz. Opacity effects due to synchrotron self-absorption can result in frequency-dependent time lags with shorter time lags at higher frequencies (Fuhrmann et al. 2014).

We note negative time lags in the correlation between radio and γ-ray light curves during P2. Previous studies have reported that this phenomenon occurs due to the interaction of a traveling shock with a stationary shock in a relativistic jet (e.g. Agudo et al. 2011a; Ramakrishnan et al. 2014). We conducted a jet kinematics analysis to explore the possibility of the shock-shock interaction in Sect. 5. In summary, our analysis suggests the presence of radio counterparts for two γ-ray flares (Flare 3 and Flare 6).

We found 3σ correlations between the γ-ray and optical light curves in P2 and P3. We observed zero time lags during P2 (Fig. 7), which is consistent with the γ-optical correlations over the entire period. However, we found negative time lags of ∼10 days between the γ-ray and optical light curves during P3 (Fig. 8). This implies that the optical emission preceded the γ-ray emission. A detailed discussion of the interpretation of negative time lags is presented in Sect. 6.2.

5. Jet kinematics

For the analysis of jet kinematics, we followed the formalism of Jorstad et al. (2017) and Weaver et al. (2022). To parameterise the compact bright features in the jet of CTA 102, we performed a model-fitting analysis of VLBA 43 GHz data taken from 2016 January 1 to 2018 August 26 (MJD 57388–58356), using the MODELFIT task in DIFMAP. We fitted circular two-dimensional Gaussians to CLEAN images, finding multiple jet components for each epoch. The fits also yield jet parameters for the components, including the flux density, the relative distance from the core (R), the size, and the relative position angle with respect to the core (PA). We calculated uncertainties of jet parameters using formulas given by Jorstad et al. (2017). An example is provided in Fig. 9, which displays the fit of Gaussian components superimposed on the CLEAN image taken on 2016 January 31. The core (labelled C0) is the brightest region presumed to be stationary at the upper end of the jet (i.e. the origin of each map). We labelled the jet components close to the core as C1–3 and the others as J1–4. The model fitting parameters of the core and jet components and their uncertainties are available on Zenodo10.

thumbnail Fig. 9.

VLBA 43 GHz image of CTA 102 taken on 2016 January 31. The contours increase by a factor of two from the lowest level of 3.92 mJy beam−1. The map peak is 1.60 Jy beam−1. The grey ellipse in the bottom right corner shows the restoring beam of 0.35 × 0.15 mas at −10 deg. The circles with crosses indicate the positions and sizes of the core (C0) and the jet components (C1 and J1–4).

We examined the jet components close to the core, at R < 0.5 mas. The top panel of Fig. 10 displays the separations of jet components from the core as a function of time. Over the course of the data period, a jet component C1 showed a quasi-stationary at ∼0.1 mas from the core, which has been consistently reported in multiple studies (e.g. Jorstad et al. 2005; Fromm et al. 2013; Casadio et al. 2015; Jorstad et al. 2017; Casadio et al. 2019; Weaver et al. 2022). Additionally, we found that two of the jet components, C2 and C3, moved farther downstream from the core. The bottom panel of Fig. 10 presents the light curves of the core and jet components of the source. To determine the motions and physical parameters for each jet component, we fitted the measurements of R for components C2 and C3 with a linear function and obtained the best-fit solution of the proper motion μ in mas/yr. The apparent speed, βapp, was obtained by using the relation that a proper motion of 1 mas/yr corresponds to 55.2c (see Sect. 1). The ejection time of a jet component, t0, is defined as the time when the centroid of the jet component coincides with the centroid of the core. We determined t0 by extrapolating to R = 0 from the fitted linear motion. The uncertainty in t0 was computed by propagating errors from the fit results. We added a minimum uncertainty of 0.005/μ yr to the uncertainty in t0, where 0.005 mas is the half-size of a pixel in the VLBA images (Weaver et al. 2022). Table 3 presents the motions of the moving jet components we calculated.

thumbnail Fig. 10.

Kinematics and light curves of the 43 GHz jet components. Top: separations of the jet components from the core (R < 0.5 mas). The grey vertical lines indicate the peaks of the γ-ray flares. The blue dotted line shows the position of C0. The cyan and red solid lines represent linear fits to the positions of C2 and C3, respectively. The cyan and red horizontal lines correspond to the uncertainties in the times of ejection of C2 and C3, respectively. Bottom: Light curves of the 43 GHz core and jet components. The cyan and red vertical lines indicate the ejection times of C2 and C3, respectively.

Table 3.

Motions of the moving jet components (MJD 57388–58356).

We derived the physical parameters of the moving jet components. To estimate the variability time scales of τvar, we followed Jorstad et al. (2017) in using a model with an exponential decay in the form of lnS(t) ∝ −t. Radio flares from blazar jets, including CTA 102, have been shown to follow exponential profiles (e.g. Terasranta & Valtaoja 1994; Rani et al. 2013; Liodakis et al. 2018a; Kim et al. 2022b). Individual exponential flares are related to the emergence of new jet components in VLBI observations, and the decaying profile of the flux density is similar to that of the total flux density (Valtaoja et al. 1999; Savolainen et al. 2002). We note that Kim et al. (2022b) identified seven exponential flares in the 15 GHz light curve, and two of them are closely connected with the emergence of VLBI jet components whose decaying flux density profile resembles that of the 15 GHz flux density (see below). Thus, the decay in flux of a jet component can be approximated as ln(S(t)/So) = k(t − tmax), where S(t) is the flux density as a function of time, tmax(≥t) is the time corresponding to the maximum flux density of Smax, S0 is the flux density of the fit to the light curve decay at tmax, and k is the slope of the fit. The variability time scales were computed to be τvar = |1/k|. To calculate the uncertainty in τvar, we used EMCEE, a PYTHON package that employs the Markov Chain Monte Carlo technique (Foreman-Mackey et al. 2013). The variability Doppler factor is given by

δ var = KsD τ var ( 1 + z ) , $$ \begin{aligned} {\delta _{\rm var} = \frac{K s D}{\tau _{\rm var} (1+z)},} \end{aligned} $$(1)

where K = 15.8 yr mas−1 Gpc−1, D is the luminosity distance in Gpc, τvar is in the unit of years, and s is the angular size of the jet component in mas, which is 1.6 times the full width at half maximum of the Gaussian modelled jet component assuming that the actual shape of the jet component is similar to a uniform face-on disk. Using the apparent speed and the variability Doppler factor, the Lorentz factor (Γvar) and the viewing angle (θvar) with respect to the line of sight are given by

Γ var = β app 2 + δ var 2 + 1 2 δ var , $$ \begin{aligned} {\Gamma _{\rm var}}&= \frac{\beta _{\rm app}^2 + \delta _{\rm var}^2 + 1}{2 \delta _{\rm var}}, \end{aligned} $$(2)

θ var = atan ( 2 β app β app 2 + δ var 2 1 ) . $$ \begin{aligned} {\theta _{\rm var}}&= \mathrm{atan}(\frac{2 \beta _{\rm app}}{\beta _{\rm app}^2 + \delta _{\rm var}^2 - 1}). \end{aligned} $$(3)

Table 4 presents the physical parameters of the jet components.

Table 4.

Physical parameters of the moving jet components (MJD 57388–58356)

Fig. 11 displays the VLBA 43 GHz images of the source in 2016–2018, including detections of C2 and C3. C2 was first detected on 2016 November 28 (MJD 57720) at a distance of 0.077 ± 0.007 mas from the core. It then travelled down the jet with an apparent speed of βapp = 9.8 ± 0.6. The variability Doppler factor of C2 is δvar = 12 ± 2. C2 was ejected in a 50 day window around 2016 June 5 ( MJD 57545 23 + 24 $ \mathrm{MJD}~57545^{+24}_{-23} $). Our estimate of t0 is roughly consistent with that of Casadio et al. (2019) within the uncertainty.

thumbnail Fig. 11.

VLBA 43 GHz images of CTA 102 in 2016–2018. The epoch of the observations is presented above each contour map. The contours increase by a factor of two from the lowest level of 5.90 mJy beam−1. All maps are restored with a common beam that is the same as in Fig. 9 (illustrated at the bottom right corner). The circles with crosses indicate the positions and sizes of the core (C0) and the jet components (C2, C3, J1, J2, J3, and J4).

The jet component C3 was first detected on 2018 March 10 (MJD 58187) at a distance of 0.042 ± 0.009 mas from the core. Emerging from the core, C3 moved downstream the jet with an apparent speed of β app = 17 . 8 2.7 + 2.0 $ \beta_{\mathrm{app}} = 17.8^{+2.0}_{-2.7} $. The variability Doppler factor of C3 is δvar = 38 ± 6. We note that this high estimate of the Doppler factor is consistent with those in previous studies (e.g. Casadio et al. 2015, 2019). C3 was ejected in a 33 day window around 2018 February 1 ( MJD 58151 16 + 17 $ \mathrm{MJD} 58151^{+17}_{-16} $).

6. Discussion

6.1. γ-ray flares in shock-shock interaction regions

We discuss the origin of the γ-ray flares based on the jet kinematics. The motion of C2, presented in Sect. 5, shows that the passage of C2 through C1 (at ∼0.1 mas from the core) coincides with the peak of Flare 3 on 2016 December 27 (MJD 57750), as well as the peaks of optical and X-ray local maxima during P2. The peak flux density of C2 matches the peak of Flare 3, indicating that the high-energy activity is likely a result of the shock-shock interaction between C1 and C2 (see also, Casadio et al. 2019). The γ-radio correlation analysis yielding negative time lags of <  − 50 days supports the scenario. During P2, the core flux density reached a peak of 3.7 ± 0.2 Jy on MJD 57636 (bottom in Fig. 10), accompanied by the size variation with a factor of ∼5 increase (from ∼0.01 to ∼0.07 mas). The increase in both the flux density and size of the core can be attributed to the ejection of C2 emerging from the core (see Fig. 10). In addition, the variation in the core flux density matches well with the radio local maxima in P2 (Fig. 2). The observed negative time lags (the radio emission leading the γ-ray emission) confirm that Flare 3 occurred downstream of the radio core.

Furthermore, we examined the interaction between C1 and C2 by modelling the trajectory of C2 using its R and PA. Fig. 12 presents the trajectory of C2 and the mean separation from the core for C1 (top) and the γ-ray light curve (bottom). We used the linear-fit result of C2 as presented in Sect. 5 and fitted a second-order polynomial function to the PA of C2. For the size of C2, we took the minimum value. The analysis revealed that the traveling shock (C2), propagating along a helical trajectory (top in Fig. 12), passed through C1 beginning at least MJD 57607 and exited from the region of influence on the shock-shock interaction at MJD 57852. The period of the shock-shock interaction is estimated to be around 250 days, coinciding with three γ-ray flares (Flares 2–4). This implies that the shock-shock interaction likely triggered not only Flare 3 but also other major γ-ray flares during P2. Throughout the interaction between C1 and C2, the separation from the core for C2 varied from ∼0.02 to ∼0.16 mas. Assuming a viewing angle of θ ∼ 4.5 deg taken from C2, the angular separation corresponds to the de-projected distance (from the core) of ∼2–17 pc. Therefore, we suggest the possibility that the γ-ray flares during P2 occurred parsecs downstream of the core. Our results are consistent with results estimated in previous studies of blazars (e.g. León-Tavares et al. 2011; Schinzel et al. 2012; Rani et al. 2015; Kramarenko et al. 2022).

thumbnail Fig. 12.

Motion and position of C2 relative to the core, with γ-ray light curve during the flaring period. Top: vector motion fit and sky position of C2. The position of C2 is relative to the core, denoted by the blue dotted lines. The black solid line represents the vector fit to the separation from the core for C2. The position at each epoch is denoted by filled cyan circles while the positions inferred from the motion of C2 are shown with unfilled large cyan circles. The unfilled purple circle represents the mean separation of C1 from the core. The filled red circles indicate the size of C2 assuming its lower limit, during the period of the shock–shock interaction. The red part in the solid line indicates the interaction region between C1 and C2. Bottom: γ-ray light curve during the flaring period (MJD 57300–58200). The shaded red region represents the period of the interaction between C1 and C2.

Fig. 13 displays the trajectory of C3 (top) and the γ-ray and 37 GHz light curves with the Bayesian blocks (bottom) during P3. The Bayesian block analysis reveals that the γ-ray light curve during P3 is characterised by two local maxima at MJD 58121 (Flare 6) and MJD 58142. The second local maximum coincides with the ejection time of C3, estimated to be MJD 58151 16 + 17 $ \mathrm{MJD}~58151^{+17}_{-16} $ within the uncertainty. The lack of 43 GHz VLBA data near the ejection time of C3 makes it challenging to study the relationship between the variation in the core flux density and the ejection of C3 from the core. Kim et al. (2022b) reported a strong correlation between the core flux variability and the single-dish flux variability of the jet in CTA 102. During P3, a local maximum at MJD 58158 in the 37 GHz light curve is consistent with the ejection time of C3. The time coincidence of the local maxima with the shock ejection suggests that the interaction between C3 and the core (possibly a standing shock) may have been associated with the γ-ray activity during P3. Similar to C2, C3 propagates down the jet along a helical trajectory due to the large variation in PA.

thumbnail Fig. 13.

Motion and position of C3 relative to the core, with γ-ray and radio light curves during the flaring period. Top: same as Fig. 12 for the jet component C3. The grey arrow indicates the ejection position of C3. Bottom: γ-ray (black) and 37 GHz (cyan) light curves during P3 (MJD 58000–58400). The black and cyan solid lines represent the Bayesian blocks of the γ-ray and 37 GHz light curves, respectively. The black vertical dashed lines indicate peaks of two γ-ray local maxima (including Flare 6) in the Bayesian block analysis. The cyan vertical line denotes the peak of the 37 GHz flux density. The shaded grey region shows the ejection of C3.

In order to constrain the location of the γ-ray flaring region, we use the multi-wavelength time lags between the radio and γ-ray light curves. The distance between the γ-ray emission region and the radio core can be calculated as Δ r radio γ = β app c t radio , γ sour / sin θ $ \Delta r_{\mathrm{radio-\gamma}} = \beta_{\mathrm{app}} c t^{\mathrm{sour}}_{\mathrm{radio, \gamma}} / {\sin{\theta}} $, where t radio , γ sour $ t^{\mathrm{sour}}_{\mathrm{radio, \gamma}} $ is the radio to γ-ray time delay in the source frame (e.g. Pushkarev et al. 2010; Ramakrishnan et al. 2014; Kim et al. 2022a; Cheong et al. 2024). We computed this time delay following the relation of t radio , γ sour = t radio , γ obs ( 1 + z ) 1 $ t^{\mathrm{sour}}_{\mathrm{radio, \gamma}} = t^{\mathrm{obs}}_{\mathrm{radio, \gamma}} (1+z)^{-1} $, where t radio , γ obs $ t^{\mathrm{obs}}_{\mathrm{radio, \gamma}} $ is the radio to γ-ray time delay in the observer frame. Regarding the association of the ejection epoch for C3 with the peak of the γ-ray flare, we used the physical jet parameters (i.e. μ, βapp, and θvar) of C3. We found that the γ-ray flaring regions are located at Δrγ, 15 GHz = 8.1 ± 2.8 pc and Δrγ, 37 GHz = 4.6 ± 2.2 pc upstream of the radio core.

The analysis of the frequency-dependent shift of the VLBI core location by Fromm et al. (2015) indicates that the de-projected separation between the SMBH and the 86 GHz core is 7.1 ± 3.3 pc. Using the results of the core shift analysis from Fromm et al. (2015), the core shift values relative to 86 GHz are Δr86, 15 GHz = 28.3 ± 12.5 pc, Δr86, 37 GHz = 7.2 ± 6.8 pc, and Δr86, 43 GHz = 5.2 ± 6.2 pc. The location of the radio core (from the SMBH) can be estimated to be 35.4 ± 12.9 pc, 14.3 ± 7.6 pc, and 12.3 ± 7.0 pc at 15, 37, and 43 GHz, respectively. These distances correspond to ∼104 − 5 gravitational radii, which are comparable to those estimated for other blazars (e.g. Marscher et al. 2008; Abdo et al. 2010c).

We discuss the seed photons for the physical mechanism that contributes to the γ-ray flares of CTA 102. Given the location of the 43 GHz core relative to the central region and the traveling shock (C2) from the core, the shock-shock interaction region associated with the γ-ray flares in 2016–2017 is located at tens of parsecs from the central engine, which is well beyond the BLR (≪1 pc, Pian et al. 2005) and the dust torus (∼2 pc, Nenkova et al. 2008). This suggests that the sources of seed photons to generate γ-rays may come from the jet itself.

The 2018 γ-ray flaring region from the SMBH is estimated to be 9.7 ± 7.9 pc. This location is far from the BLR, but it could be comparable to the size of the torus, suggesting that the dust torus may have provided the seed photons for the EC mechanism during Flare 6. The core position could vary due to flare events and newly ejected jet components propagating along the jet (e.g. Plavin et al. 2019). The typical variability of the core position has been reported to be ∼2 pc (Plavin et al. 2019), which is not likely to have significant effects on our results.

6.2. γ-ray/optical correlations

We investigated the relationship between the γ-ray and optical fluxes to understand the emission mechanisms of the γ-ray flares. Most cases of correlations between γ-ray and optical light curves were found to be significant with confidence levels > 99.7% (Sect. 4.2). In a simple leptonic scenario (with constant magnetic field strength and viewing angle for a flaring period), it is assumed that synchrotron radiation from the jet is mostly emitted at optical wavelengths, whereas the γ-ray flaring emission is from IC scattering of optical/IR photons by relativistic electrons in the jet (see, e.g. Larionov et al. 2016; Liodakis et al. 2019; Sahakyan et al. 2022). Under this assumption, the relation between the flux density of the IC scattered emission (FIC) and that of the synchrotron radiation (Fsync) is FIC ∝ Fsyncm, where m is a slope of the flux–flux relation. If the dominant γ-ray emission mechanism is EC, then the only component of the synchrotron emission that also contributes to the IC emission is the population of relativistic electrons, in which case, m = 1. If it is SSC, both the relativistic electrons and their associated synchrotron photons contribute to the IC emission, resulting in an additional power (i.e. m = 2) in the flux–flux relation.

D’Ammando et al. (2019) investigated the relationship between the gamma-ray and optical fluxes on a logarithmic space during MJD 56293–57793 (including P1–2), yielding F IC F sync 0.82 $ F_{\mathrm{IC}} \propto F_{\mathrm{sync}}^{0.82} $ with a Pearson correlation coefficient of rp = 0.95. We conducted the analysis during the entire period (MJD 56200–58400). For this, we converted the unit of the γ-ray flux density (ph cm−2 s−1) into mJy to be matched with that of the optical flux density. The optical flux density was shifted close in time with respect to the γ-ray flux density, regarding τDCF = 3 ± 5 days. As seen in Fig. 14, we found a slope of m = 0.87 ± 0.04 with rp = 0.85, which is consistent with the result of D’Ammando et al. (2019). The slope close to the unity corresponds to the EC mechanism under the simple scenario.

thumbnail Fig. 14.

Quasi-simultaneous γ-ray versus optical (R-band) fluxes in log–log space. The data was taken from the six-year observation period in γ-rays (MJD 56200–58400). The best-fit solution of the flux–flux relation shows F IC F sync 0.87 ± 0.04 $ F_{\mathrm{IC}} \propto F_{\mathrm{sync}}^{0.87 \pm 0.04} $ (plotted in black and red).

Considering the γ-ray emission site near the dust torus during Flare 6 (Sect. 6.1), we assume that a γ-ray photon is produced by the upscattering of an infrared (IR) photon by a relativistic electron, expressed as ϵs ≳ γe2ϵi (Boettcher et al. 2012), where ϵs is the scattered γ-ray photon energy, γe is the electron Lorentz factor, and ϵi is the incident IR photon energy. Geng et al. (2022) estimated the maximum γ-ray photon energy during Flare 6 to be ϵs = 30.9 GeV. The IR photon energy is predicted to be ∼102 − 3 eV, corresponding to a wavelength of 1.7–78.6 μm (Malmrose et al. 2011). This suggests γe ≲ 104, which is roughly consistent with the predictions from the SED analysis of blazars (Kang et al. 2014; Pandey et al. 2024).

However, we note that a unity slope is possible within the SSC framework if variations in the magnetic field or the Doppler factor (resulting from changes in the jet viewing angle) are considered (Larionov et al. 2016; Liodakis et al. 2019; D’Ammando et al. 2019; Acciari et al. 2020). The multi-wavelength light curve analysis yielded a viewing angle of ∼8 deg in the 15 GHz band and ∼2 deg in the optical or γ-ray bands (D’Ammando et al. 2019), suggesting that the unity slope for the γ-ray flares in P1–2 (2016–2017) may be due to changes in the Doppler factor resulting from variations in the jet viewing angle. The γ-ray variability during Flare 6 may be associated with these variations. We found a variation in the viewing angle from ∼4.5 deg to ∼1.2 deg (Table 4), as well as the large variation in the PA of C3 from 89.1 ± 5.0 deg to 129.2 ± 4.1 deg. Assuming that the observed flux density (Fν) is amplified due to the Doppler boosting effect, Fν is related to the flux density in the jet frame (Fν) by Fν = δ3 − αFν, where Fν ∝ να, δ is the Doppler factor defined as δ = 1/Γ(1 − β cos θ), α is the optically thin spectral index (α < 0), β is the intrinsic velocity of the jet, and Γ is the bulk Lorentz factor defined as Γ = 1 / 1 β 2 $ \Gamma = 1/\sqrt{1 - \beta^2} $ (Blandford & Königl 1979). The ratio between the maximum and minimum flux densities in the optical during Flare 6 is ∼12. We thus obtain a ratio between the maximum and minimum Doppler factors of δmax/δmin < 2.3. The maximum Doppler factor of C3 is δmax ∼ 44 at θ ∼ 0 deg due to Γvar ≈ 23 (see Table 4). This yields a lower limit of > 19 for the minimum Doppler factor at θ < 3 deg.

We discuss the observed negative time lags between the γ-ray and optical light curves during P3 (Sect. 4.2), which have been reported in previous studies as well (e.g. Abdo et al. 2010a; Agudo et al. 2011b; Hayashida et al. 2012; Gupta et al. 2017; Liodakis et al. 2019). Fig. 15 displays the γ-ray and optical (R-band) light curves with normalised flux densities and the Bayesian block representation during P3. In this period, we identified seven optical and five γ-ray local maxima (including Flare 6). After applying a time shift equal to the observed time lag, the optical local maxima at MJD 58097 and MJD 58142 appear to have no γ-ray counterparts. The negative time lags are likely due to the cross-match of two optical local maxima at MJD 58104 and MJD 58132 with two γ-ray local maxima at MJD 58121 (Flare 6) and MJD 58142. Assuming a co-spatial origin of optical and γ emission, the optical local maxima at MJD 58132 and MJD 58142 seem to overlap with the two γ-ray local maxima. The presence of optical local maxima without γ-ray counterparts may be attributed to magnetic field changes leading to synchrotron emission modulation with steady IC emission (e.g. Chatterjee et al. 2013; Cohen et al. 2014). High-cadence observations and modelling of the broadband SED are needed for a clear understanding of the γ-optical correlation.

thumbnail Fig. 15.

Light curves in the γ-ray and optical (R) bands during P3. The flux densities have been normalised to the maximum value in each band. The solid lines correspond to the Bayesian block representation of the light curves. The vertical dashed lines indicate peaks of the light curves.

The observed negative time lags could be explained in the framework of the one-zone leptonic EC scenario. In this scenario, we assume that the same electron population produces the optical and γ-ray emission by the synchrotron and EC processes in the fast-cooling regime, that is, the electron cooling time scales are much shorter than the source light crossing time scales (e.g. Derishev 2007; Hayashida et al. 2012; Janiak et al. 2012). The synchrotron and EC losses depend on the energy densities of the magnetic field and the external radiation field (e.g. the BLR, dusty torus, etc.), respectively. The energy densities could have different spatial distributions along the jet (refer to Fig. 6 in Janiak et al. 2012). If the magnetic field energy density decreases more rapidly than the external radiation field energy density, the EC (γ-rays) peak could be delayed relative to the synchrotron (optical) peak, resulting in a negative time lag.

7. Summary

The blazar CTA 102 exhibited intense γ-ray flaring activity from 2016 to 2018. We investigated the physical origins of γ-ray flares in CTA 102 by performing a cross-correlation analysis between its γ-ray light curve and lower-energy light curves from multi-wavelength data. Using VLBA observations, we explored the pc-scale jet activity during the flaring periods. We summarise our conclusions as follows:

  1. Based on our criteria, we identified six flares in the γ-ray light curve. Assuming a power-law spectrum in the γ-ray band, we found changes in the spectral index, exhibiting a harder-when-brighter trend, during the flaring period from 2016 to 2018.

  2. We conducted a cross-correlation analysis between the γ-ray light curve and those in lower-energy bands for both the entire observation period and the flaring period. The analysis revealed that the time lags between the γ-ray light curve and the optical and X-ray light curves were nearly zero, implying a time coincidence and thus a co-spatial origin.

  3. The correlations between the γ-ray and radio light curves were more complex over long periods, possibly due to more rapid variability in the γ-ray band compared to the radio band. We divided the flaring period into three parts (P1–3) and explored the correlations between the γ-ray and radio light curves. There are no significant correlations between the γ-ray and radio light curves in P1. The brightest γ-ray flare (Flare 3) showed significant correlations with the radio (15 and 37 GHz) local maxima as well as a negative time lag, indicating that the γ-ray emission lagged behind the radio emission. The most recent γ-ray flare (Flare 6) in P3 exhibited significant correlations with radio flares, with positive time lags between the γ-ray and radio light curves. This implies that the γ-ray emission preceded the radio emission.

  4. We used VLBA maps and found that two jet components (C2 and C3) were newly ejected from the core during P1 and P3, respectively. Right after the ejection of C2, there was an increase in both the flux density and size of the core. The variation in core flux density seems to be linked to the radio local maxima correlated with Flare 3. The interaction between the moving jet component (C2) and a quasi-stationary jet component C1 (at ∼0.1 mas) is associated with Flare 3. Negative time lags in the γ-radio correlation support the shock-shock interaction scenario. In addition to Flare 3, multiple γ-ray flares in P2 could be expected to occur as the traveling shock passes through the stationary one, at ∼2–17 pc downstream of the radio core.

  5. The ejection time of the jet component C3 coincided with the γ-ray activity in P3, suggesting a connection between the two phenomena. The γ-ray flaring region is located parsecs upstream of the radio core as inferred from multi-wavelength time lags between the radio and γ-ray light curves and physical jet parameters of C3. The γ-ray flaring regions are well beyond the BLR, which suggests that the sources of seed photons for γ-rays may come from either the jet itself or the dust torus.

  6. We found a linear relationship between the γ-ray and the optical fluxes, suggesting that the EC mechanism is at play. Alternatively, this relationship may be consistent with the SSC mechanism as well, provided that the jet viewing angle varies as the jet propagates.

Data availability

Appendix A, containing the table used in this publication, is available on Zenodo at: https://doi.org/10.5281/zenodo.14642338.


Acknowledgments

We thank the anonymous referee for useful comments that helped improve this manuscript. This research has made use of data from the OVRO 40-m monitoring program (Richards et al. 2011), supported by private funding from the California Institute of Technology and the Max Planck Institute for Radio Astronomy, and by NASA grants NNX08AW31G, NNX11A043G, and NNX14AQ89G and NSF grants AST-0808050 and AST- 1109911. This publication makes use of data obtained at the Metsähovi Radio Observatory, operated by Aalto University in Finland. The IAA-CSIC group acknowledges financial support from the grant CEX2021-001131-S funded by MCIN/AEI/10.13039/501100011033 to the “Instituto de Astrofísica de Andalucía-CSIC”. Acquisition and reduction of the POLAMI data was supported in part by MICIN through grants PID2019-107847RB-C44 abd PID2022-139117NB-C44. The POLAMI observations were carried out at the IRAM 30 m Telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). CC acknowledges support by the European Research Council (ERC) under the HORIZON ERC Grants 2021 programme under grant agreement No. 101040021. The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. We recognize that Maunakea is a culturally important site for the indigenous Hawaiian people; we are privileged to study the cosmos from its summit. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00001.CAL. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This study makes use of VLBA data from the VLBA-BU Blazar Monitoring Program (BEAM-ME and VLBA-BU-BLAZAR; http://www.bu.edu/blazars/BEAM-ME.html), funded by NASA through the Fermi Guest Investigator Program. The VLBA is an instrument of the National Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated by Associated Universities, Inc. Data from the Steward Observatory spectropolarimetric monitoring project were used. This program is supported by Fermi Guest Investigator grants NNX08AW56G, NNX09AU10G, NNX12AO93G, and NNX15AU81G. The Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat à l’Energie Atomique and the Centre National de la Recherche Scientifique / Institut National de Physique Nucléaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d’Études Spatiales in France. This work performed in part under DOE Contract DE-AC02-76SF00515. This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MIST) (2020R1A2C2009003). This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT; RS-2024-00449206). This research has been supported by the POSCO Science Fellowship of POSCO TJ Park Foundation.

References

  1. Abdo, A. A., Ackermann, M., Agudo, I., et al. 2010a, ApJ, 721, 1425 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010b, ApJS, 188, 405 [NASA ADS] [CrossRef] [Google Scholar]
  3. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010c, Nature, 463, 919 [NASA ADS] [CrossRef] [Google Scholar]
  4. Abdollahi, S., Acero, F., Ackermann, M., et al. 2020, ApJS, 247, 33 [Google Scholar]
  5. Acciari, V. A., Ansoldi, S., Antonelli, L. A., et al. 2020, ApJS, 248, 29 [Google Scholar]
  6. Acero, F., Ackermann, M., Ajello, M., et al. 2015, ApJS, 218, 23 [Google Scholar]
  7. Agudo, I., Jorstad, S. G., Marscher, A. P., et al. 2011a, ApJ, 726, L13 [NASA ADS] [CrossRef] [Google Scholar]
  8. Agudo, I., Marscher, A. P., Jorstad, S. G., et al. 2011b, ApJ, 735, L10 [NASA ADS] [CrossRef] [Google Scholar]
  9. Agudo, I., Thum, C., Molina, S. N., et al. 2018a, MNRAS, 474, 1427 [NASA ADS] [CrossRef] [Google Scholar]
  10. Agudo, I., Thum, C., Ramakrishnan, V., et al. 2018b, MNRAS, 473, 1850 [Google Scholar]
  11. Algaba, J.-C., Lee, S.-S., Kim, D.-W., et al. 2018, ApJ, 852, 30 [Google Scholar]
  12. Angel, J. R. P., & Stockman, H. S. 1980, ARA&A, 18, 321 [NASA ADS] [CrossRef] [Google Scholar]
  13. Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 [CrossRef] [Google Scholar]
  14. Barthelmy, S. D., Barbier, L. M., Cummings, J. R., et al. 2005, Space Sci. Rev., 120, 143 [Google Scholar]
  15. Blandford, R. D., & Königl, A. 1979, ApJ, 232, 34 [Google Scholar]
  16. Blandford, R. D., McKee, C. F., & Rees, M. J. 1977, Nature, 267, 211 [NASA ADS] [CrossRef] [Google Scholar]
  17. Blandford, R., Meier, D., & Readhead, A. 2019, ARA&A, 57, 467 [NASA ADS] [CrossRef] [Google Scholar]
  18. Boettcher, M., Harris, D. E., & Krawczynski, H. 2012, Relativistic Jets from Active Galactic Nuclei (Wiley-VCH) [CrossRef] [Google Scholar]
  19. Casadio, C., Gómez, J. L., Jorstad, S. G., et al. 2015, ApJ, 813, 51 [Google Scholar]
  20. Casadio, C., Marscher, A. P., Jorstad, S. G., et al. 2019, A&A, 622, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Chatterjee, R., Fossati, G., Urry, C. M., et al. 2013, ApJ, 763, L11 [NASA ADS] [CrossRef] [Google Scholar]
  22. Cheong, W. Y., Lee, S.-S., Kim, S.-H., et al. 2024, MNRAS, 527, 882 [Google Scholar]
  23. Ciprini, S. 2016, The Astronomer’s Telegram, 9869, 1 [Google Scholar]
  24. Cohen, D. P., Romani, R. W., Filippenko, A. V., et al. 2014, ApJ, 797, 137 [NASA ADS] [CrossRef] [Google Scholar]
  25. Connolly, S. D. 2015, ArXiv e-prints [arXiv:1503.06676] [Google Scholar]
  26. D’Ammando, F., Raiteri, C. M., Villata, M., et al. 2019, MNRAS, 490, 5300 [Google Scholar]
  27. Derishev, E. V. 2007, Ap&SS, 309, 157 [NASA ADS] [CrossRef] [Google Scholar]
  28. Dermer, C. D., Schlickeiser, R., & Mastichiadis, A. 1992, A&A, 256, L27 [NASA ADS] [Google Scholar]
  29. Dondi, L., & Ghisellini, G. 1995, MNRAS, 273, 583 [NASA ADS] [Google Scholar]
  30. Edelson, R. A., & Krolik, J. H. 1988, ApJ, 333, 646 [Google Scholar]
  31. Emmanoulopoulos, D., McHardy, I. M., & Papadakis, I. E. 2013, MNRAS, 433, 907 [NASA ADS] [CrossRef] [Google Scholar]
  32. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  33. Fromm, C. M., Ros, E., Perucho, M., et al. 2013, A&A, 551, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Fromm, C. M., Perucho, M., Ros, E., Savolainen, T., & Zensus, J. A. 2015, A&A, 576, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Fuhrmann, L., Larsson, S., Chiang, J., et al. 2014, MNRAS, 441, 1899 [NASA ADS] [CrossRef] [Google Scholar]
  36. Gasparyan, S., Sahakyan, N., Baghmanyan, V., & Zargaryan, D. 2018, ApJ, 863, 114 [NASA ADS] [CrossRef] [Google Scholar]
  37. Gehrels, N., Chincarini, G., Giommi, P., et al. 2004, ApJ, 611, 1005 [Google Scholar]
  38. Geng, X., Ding, N., Cao, G., et al. 2022, ApJS, 260, 48 [NASA ADS] [CrossRef] [Google Scholar]
  39. Ghisellini, G., Tavecchio, F., & Chiaberge, M. 2005, A&A, 432, 401 [CrossRef] [EDP Sciences] [Google Scholar]
  40. Gupta, A. C., Mangalam, A., Wiita, P. J., et al. 2017, MNRAS, 472, 788 [NASA ADS] [CrossRef] [Google Scholar]
  41. Gurwell, M. A., Peck, A. B., Hostler, S. R., Darrah, M. R., & Katz, C. A. 2007, in From Z-Machines to ALMA: (Sub)Millimeter Spectroscopy of Galaxies, eds. A. J. Baker, J. Glenn, A. I. Harris, J. G. Mangum, & M. S. Yun, ASP Conf. Ser., 375, 234 [NASA ADS] [Google Scholar]
  42. Hayashida, M., Madejski, G. M., Nalewajko, K., et al. 2012, ApJ, 754, 114 [NASA ADS] [CrossRef] [Google Scholar]
  43. Hodgson, J. A., Rani, B., Lee, S.-S., et al. 2018, MNRAS, 475, 368 [NASA ADS] [CrossRef] [Google Scholar]
  44. Hovatta, T., & Lindfors, E. 2019, New Astron. Rev., 87, 101541P [Google Scholar]
  45. Hovatta, T., Valtaoja, E., Tornikoski, M., & Lähteenmäki, A. 2009, A&A, 494, 527 [CrossRef] [EDP Sciences] [Google Scholar]
  46. Janiak, M., Sikora, M., Nalewajko, K., Moderski, R., & Madejski, G. M. 2012, ApJ, 760, 129 [NASA ADS] [CrossRef] [Google Scholar]
  47. Jorstad, S. G., Marscher, A. P., Lister, M. L., et al. 2005, AJ, 130, 1418 [Google Scholar]
  48. Jorstad, S. G., Marscher, A. P., Larionov, V. M., et al. 2010, ApJ, 715, 362 [Google Scholar]
  49. Jorstad, S. G., Marscher, A. P., Morozova, D. A., et al. 2017, ApJ, 846, 98 [Google Scholar]
  50. Kang, S.-J., Chen, L., & Wu, Q. 2014, ApJS, 215, 5 [NASA ADS] [CrossRef] [Google Scholar]
  51. Kaur, N., & Baliyan, K. S. 2018, A&A, 617, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Kim, D.-W., Kravchenko, E. V., Kutkin, A. M., et al. 2022a, ApJ, 925, 64 [NASA ADS] [CrossRef] [Google Scholar]
  53. Kim, S.-H., Lee, S.-S., Lee, J. W., et al. 2022b, MNRAS, 510, 815 [Google Scholar]
  54. Kramarenko, I. G., Pushkarev, A. B., Kovalev, Y. Y., et al. 2022, MNRAS, 510, 469 [Google Scholar]
  55. Larionov, V. M., Villata, M., Raiteri, C. M., et al. 2016, MNRAS, 461, 3047 [NASA ADS] [CrossRef] [Google Scholar]
  56. León-Tavares, J., Valtaoja, E., Tornikoski, M., Lähteenmäki, A., & Nieppola, E. 2011, A&A, 532, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Li, X., Mohan, P., An, T., et al. 2018, ApJ, 854, 17 [NASA ADS] [CrossRef] [Google Scholar]
  58. Liodakis, I., Hovatta, T., Huppenkothen, D., et al. 2018a, ApJ, 866, 137 [Google Scholar]
  59. Liodakis, I., Romani, R. W., Filippenko, A. V., et al. 2018b, MNRAS, 480, 5517 [NASA ADS] [CrossRef] [Google Scholar]
  60. Liodakis, I., Romani, R. W., Filippenko, A. V., Kocevski, D., & Zheng, W. 2019, ApJ, 880, 32 [NASA ADS] [CrossRef] [Google Scholar]
  61. Lisakov, M. M., Kovalev, Y. Y., Savolainen, T., Hovatta, T., & Kutkin, A. M. 2017, MNRAS, 468, 4478 [Google Scholar]
  62. MacDonald, N. R., Marscher, A. P., Jorstad, S. G., & Joshi, M. 2015, ApJ, 804, 111 [Google Scholar]
  63. Madejski, G. G., & Sikora, M. 2016, ARA&A, 54, 725 [NASA ADS] [CrossRef] [Google Scholar]
  64. Malmrose, M. P., Marscher, A. P., Jorstad, S. G., Nikutta, R., & Elitzur, M. 2011, ApJ, 732, 116 [NASA ADS] [CrossRef] [Google Scholar]
  65. Maraschi, L., Ghisellini, G., & Celotti, A. 1992, ApJ, 397, L5 [CrossRef] [Google Scholar]
  66. Marscher, A. P., Jorstad, S. G., D’Arcangelo, F. D., et al. 2008, Nature, 452, 966 [Google Scholar]
  67. Max-Moerbeck, W., Hovatta, T., Richards, J. L., et al. 2014a, MNRAS, 445, 428 [Google Scholar]
  68. Max-Moerbeck, W., Richards, J. L., Hovatta, T., et al. 2014b, MNRAS, 445, 437 [NASA ADS] [CrossRef] [Google Scholar]
  69. Mead, A. R. G., Ballard, K. R., Brand, P. W. J. L., et al. 1990, A&AS, 83, 183 [Google Scholar]
  70. Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezić, Ž., & Elitzur, M. 2008, ApJ, 685, 160 [Google Scholar]
  71. Nolan, P. L., Abdo, A. A., Ackermann, M., et al. 2012, ApJS, 199, 31 [Google Scholar]
  72. Pandey, A., Kushwaha, P., Wiita, P. J., et al. 2024, A&A, 681, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Park, J.-H., & Trippe, S. 2014, ApJ, 785, 76 [NASA ADS] [CrossRef] [Google Scholar]
  74. Park, J., & Trippe, S. 2017, ApJ, 834, 157 [CrossRef] [Google Scholar]
  75. Peterson, B. M., Wanders, I., Horne, K., et al. 1998, PASP, 110, 660 [Google Scholar]
  76. Peterson, B. M., Ferrarese, L., Gilbert, K. M., et al. 2004, ApJ, 613, 682 [Google Scholar]
  77. Pian, E., Falomo, R., & Treves, A. 2005, MNRAS, 361, 919 [NASA ADS] [CrossRef] [Google Scholar]
  78. Planck Collaboration VI 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Plavin, A. V., Kovalev, Y. Y., Pushkarev, A. B., & Lobanov, A. P. 2019, MNRAS, 485, 1822 [Google Scholar]
  80. Prince, R. 2020, ApJ, 890, 164 [NASA ADS] [CrossRef] [Google Scholar]
  81. Prince, R., Raman, G., Hahn, J., Gupta, N., & Majumdar, P. 2018, ApJ, 866, 16 [NASA ADS] [CrossRef] [Google Scholar]
  82. Pushkarev, A. B., Kovalev, Y. Y., & Lister, M. L. 2010, ApJ, 722, L7 [NASA ADS] [CrossRef] [Google Scholar]
  83. Raiteri, C. M., Villata, M., Acosta-Pulido, J. A., et al. 2017, Nature, 552, 374 [NASA ADS] [CrossRef] [Google Scholar]
  84. Ramakrishnan, V., León-Tavares, J., Rastorgueva-Foi, E. A., et al. 2014, MNRAS, 445, 1636 [NASA ADS] [CrossRef] [Google Scholar]
  85. Ramakrishnan, V., Hovatta, T., Nieppola, E., et al. 2015, MNRAS, 452, 1280 [CrossRef] [Google Scholar]
  86. Rani, B., Krichbaum, T. P., Fuhrmann, L., et al. 2013, A&A, 552, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Rani, B., Krichbaum, T. P., Marscher, A. P., et al. 2015, A&A, 578, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Rani, B., Jorstad, S. G., Marscher, A. P., et al. 2018, ApJ, 858, 80 [Google Scholar]
  89. Readhead, A. C. S., Cohen, M. H., Pearson, T. J., & Wilkinson, P. N. 1978, Nature, 276, 768 [Google Scholar]
  90. Richards, J. L., Max-Moerbeck, W., Pavlidou, V., et al. 2011, ApJS, 194, 29 [Google Scholar]
  91. Robertson, D. R. S., Gallo, L. C., Zoghbi, A., & Fabian, A. C. 2015, MNRAS, 453, 3455 [NASA ADS] [Google Scholar]
  92. Sahakyan, N. 2020, A&A, 635, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Sahakyan, N., Israyelyan, D., Harutyunyan, G., et al. 2022, MNRAS, 517, 2757 [NASA ADS] [CrossRef] [Google Scholar]
  94. Savolainen, T., Wiik, K., Valtaoja, E., Jorstad, S. G., & Marscher, A. P. 2002, A&A, 394, 851 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. Scargle, J. D. 1982, ApJ, 263, 835 [Google Scholar]
  96. Scargle, J. D., Norris, J. P., Jackson, B., & Chiang, J. 2013, ApJ, 764, 167 [Google Scholar]
  97. Schinzel, F. K., Lobanov, A. P., Taylor, G. B., et al. 2012, A&A, 537, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  98. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
  99. Shepherd, M. C. 1997, in Astronomical Data Analysis Software and Systems VI, eds. G. Hunt, & H. Payne, ASP Conf. Ser., 125, 77 [Google Scholar]
  100. Shukla, A., Mannheim, K., Patel, S. R., et al. 2018, ApJ, 854, L26 [NASA ADS] [CrossRef] [Google Scholar]
  101. Sikora, M., Stawarz, Ł., Moderski, R., Nalewajko, K., & Madejski, G. M. 2009, ApJ, 704, 38 [CrossRef] [Google Scholar]
  102. Smith, P. S., Montiel, E., Rightley, S., et al. 2009, ArXiv e-prints [arXiv:0912.3621] [Google Scholar]
  103. Tavecchio, F., Maraschi, L., & Ghisellini, G. 1998, ApJ, 509, 608 [Google Scholar]
  104. Terasranta, H., & Valtaoja, E. 1994, A&A, 283, 51 [NASA ADS] [Google Scholar]
  105. Teraesranta, H., Tornikoski, M., Mujunen, A., et al. 1998, A&AS, 132, 305 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  106. Thum, C., Agudo, I., Molina, S. N., et al. 2018, MNRAS, 473, 2506 [NASA ADS] [CrossRef] [Google Scholar]
  107. Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 [NASA ADS] [CrossRef] [Google Scholar]
  108. Uttley, P., McHardy, I. M., & Papadakis, I. E. 2002, MNRAS, 332, 231 [Google Scholar]
  109. Valtaoja, E., Lähteenmäki, A., Teräsranta, H., & Lainela, M. 1999, ApJS, 120, 95 [NASA ADS] [CrossRef] [Google Scholar]
  110. Vaughan, S., Edelson, R., Warwick, R. S., & Uttley, P. 2003, MNRAS, 345, 1271 [Google Scholar]
  111. Weaver, Z. R., Jorstad, S. G., Marscher, A. P., et al. 2022, ApJS, 260, 12 [NASA ADS] [CrossRef] [Google Scholar]
  112. Williamson, K. E., Jorstad, S. G., Marscher, A. P., et al. 2014, ApJ, 789, 135 [NASA ADS] [CrossRef] [Google Scholar]
  113. Wright, E. L. 2006, PASP, 118, 1711 [NASA ADS] [CrossRef] [Google Scholar]
  114. Zacharias, M., Böttcher, M., Jankowsky, F., et al. 2017, ApJ, 851, 72 [NASA ADS] [CrossRef] [Google Scholar]
  115. Zacharias, M., Böttcher, M., Jankowsky, F., et al. 2019, ApJ, 871, 19 [Google Scholar]
  116. Zamaninasab, M., Clausen-Brown, E., Savolainen, T., & Tchekhovskoy, A. 2014, Nature, 510, 126 [Google Scholar]

All Tables

Table 1.

γ-ray flares in 2016–2018 (MJD 57300–58400)

Table 2.

Cross-correlation results for γ-ray versus multi-wavelength (radio and optical) pairs during P2 and P3 (MJD 57550–58400).

Table 3.

Motions of the moving jet components (MJD 57388–58356).

Table 4.

Physical parameters of the moving jet components (MJD 57388–58356)

All Figures

thumbnail Fig. 1.

Multi-wavelength light curves from CTA 102. The entire period is from 2012 September 30 to 2018 October 9 (MJD 56200–58400). From top to bottom, the panels show light curves in γ-rays, X-rays, optical, and radio with different colours and symbols indicating different frequencies. The red triangles in the γ-ray light curve indicate upper limits on the flux density. The optical and radio light curves are displayed with arbitrary offsets for clarity. The grey horizontal lines for the γ-ray light curve denote the flux density threshold (see the text for details). The black dotted lines indicate the peaks of the γ-ray flares (1–6). The grey shaded areas mark three periods (P1–3) for the cross-correlation analysis.

In the text
thumbnail Fig. 2.

Bayesian block representation of the multi-wavelength light curves during the flaring period. The time range of the period is from 2015 October 5 to 2018 October 9 (MJD 57300–58400). The units of the flux densities and vertical lines are the same as in Fig. 1. The horizontal arrows indicate the time ranges of flaring periods (P1–3).

In the text
thumbnail Fig. 3.

γ-ray light curve (top) and time series for the γ-ray photon index (bottom). Black dotted lines (1–6) indicate the peaks of the γ-ray flares.

In the text
thumbnail Fig. 4.

γ-ray photon index (Γ) versus γ-ray flux density (Fγ). The Pearson correlation coefficient rp is presented with the confidence level.

In the text
thumbnail Fig. 5.

Discrete correlation functions for γ-ray versus optical (R and V bands) and X-ray over the entire period, from top left to bottom right. The red solid and green dashed lines indicate 2σ (95%) and 3σ (99.7%) confidence intervals, respectively. Positive time lags indicate that the flux density in γ-rays leads the flux density in lower-energy bands.

In the text
thumbnail Fig. 6.

Discrete correlation functions for γ-ray versus radio (15 and 89 GHz) over the entire period. Horizontal lines are the same as in Fig. 5.

In the text
thumbnail Fig. 7.

Discrete correlation functions for γ-ray versus radio (15 and 37 GHz) and optical during P2. Horizontal lines are the same as in Fig. 5.

In the text
thumbnail Fig. 8.

Discrete correlation functions for γ-ray versus radio (15, 37, 89, and 230 GHz) and optical (R and V bands) in P3 from top left to bottom right. Horizontal lines are the same as in Fig. 5.

In the text
thumbnail Fig. 9.

VLBA 43 GHz image of CTA 102 taken on 2016 January 31. The contours increase by a factor of two from the lowest level of 3.92 mJy beam−1. The map peak is 1.60 Jy beam−1. The grey ellipse in the bottom right corner shows the restoring beam of 0.35 × 0.15 mas at −10 deg. The circles with crosses indicate the positions and sizes of the core (C0) and the jet components (C1 and J1–4).

In the text
thumbnail Fig. 10.

Kinematics and light curves of the 43 GHz jet components. Top: separations of the jet components from the core (R < 0.5 mas). The grey vertical lines indicate the peaks of the γ-ray flares. The blue dotted line shows the position of C0. The cyan and red solid lines represent linear fits to the positions of C2 and C3, respectively. The cyan and red horizontal lines correspond to the uncertainties in the times of ejection of C2 and C3, respectively. Bottom: Light curves of the 43 GHz core and jet components. The cyan and red vertical lines indicate the ejection times of C2 and C3, respectively.

In the text
thumbnail Fig. 11.

VLBA 43 GHz images of CTA 102 in 2016–2018. The epoch of the observations is presented above each contour map. The contours increase by a factor of two from the lowest level of 5.90 mJy beam−1. All maps are restored with a common beam that is the same as in Fig. 9 (illustrated at the bottom right corner). The circles with crosses indicate the positions and sizes of the core (C0) and the jet components (C2, C3, J1, J2, J3, and J4).

In the text
thumbnail Fig. 12.

Motion and position of C2 relative to the core, with γ-ray light curve during the flaring period. Top: vector motion fit and sky position of C2. The position of C2 is relative to the core, denoted by the blue dotted lines. The black solid line represents the vector fit to the separation from the core for C2. The position at each epoch is denoted by filled cyan circles while the positions inferred from the motion of C2 are shown with unfilled large cyan circles. The unfilled purple circle represents the mean separation of C1 from the core. The filled red circles indicate the size of C2 assuming its lower limit, during the period of the shock–shock interaction. The red part in the solid line indicates the interaction region between C1 and C2. Bottom: γ-ray light curve during the flaring period (MJD 57300–58200). The shaded red region represents the period of the interaction between C1 and C2.

In the text
thumbnail Fig. 13.

Motion and position of C3 relative to the core, with γ-ray and radio light curves during the flaring period. Top: same as Fig. 12 for the jet component C3. The grey arrow indicates the ejection position of C3. Bottom: γ-ray (black) and 37 GHz (cyan) light curves during P3 (MJD 58000–58400). The black and cyan solid lines represent the Bayesian blocks of the γ-ray and 37 GHz light curves, respectively. The black vertical dashed lines indicate peaks of two γ-ray local maxima (including Flare 6) in the Bayesian block analysis. The cyan vertical line denotes the peak of the 37 GHz flux density. The shaded grey region shows the ejection of C3.

In the text
thumbnail Fig. 14.

Quasi-simultaneous γ-ray versus optical (R-band) fluxes in log–log space. The data was taken from the six-year observation period in γ-rays (MJD 56200–58400). The best-fit solution of the flux–flux relation shows F IC F sync 0.87 ± 0.04 $ F_{\mathrm{IC}} \propto F_{\mathrm{sync}}^{0.87 \pm 0.04} $ (plotted in black and red).

In the text
thumbnail Fig. 15.

Light curves in the γ-ray and optical (R) bands during P3. The flux densities have been normalised to the maximum value in each band. The solid lines correspond to the Bayesian block representation of the light curves. The vertical dashed lines indicate peaks of the light curves.

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.