X-ray Properties of 3C 111: Separation of Primary Nuclear Emission and Jet Continuum

: 3C111 is BLRG with signatures of both FSRQ and Sy1 in X-ray spectrum. The signiﬁcant X-ray observational dataset was collected for it by INTEGRAL, XMM-Newton, SWIFT, Suzaku and others. The overall X-ray spectrum of 3C 111 shows signs of a peculiarity with the large value of the high-energy cut-off typical rather for RQ AGN, probably due to the jet contamination. Separating the jet counterpart in the X-ray spectrum of 3C 111 from the primary nuclear counterpart can answer the question is this nucleus truly peculiar or this is a fake “peculiarity” due to a signiﬁcant jet contribution. In view of this question, our aim is to estimate separately the accretion disk/corona and non-thermal jet emission in the 3C 111 X-ray spectra within different observational periods. To separate the disk/corona and jet contributions in total continuum, we use the idea that radio and X-ray spectra of jet emission can be described by a simple power-law model with the same photon index. This additional information allows us to derive rather accurate values of these contributions. In order to test these results, we also consider relations between the nuclear continuum and the line emission.


Introduction
The fact that active galactic nuclei (AGN) can manifest itself as radio loud (RL) or radio quiet (RQ) is clearly known from 1989 when in Reference [1] the criterion of radio loudness was formulated. This criterion is based on radio-to-optical brightness relation R = L R /L V where L R is the AGN luminosity in radio frequencies (usually at 5 GHz) and L V is its luminosity in blue optical light band (near 689 THz) R ≤ 1 means that the object is RQ; R ≥ 10 corresponds to an RL AGN. R = 1 − 10 can be referred to as radio moderate (RM). The last class with strong variations of radio loudness from typically RL to RQ values and unstable jets is becoming more common in the last decade is radio transitional (RT), see Reference [2]. RT AGN possibly have variable spectral shape similarly to X-ray binaries [3][4][5]. An object of our investigation 3C 111 likely belongs to this class rather than to RL one, as well as several other AGNs including the famous 3C 120 [6].
Since that time we know that RQ AGN are 10-20 times more frequent than RL ones. However, the question of the physical difference in AGN "central engine" leading to particular radio properties, is still quite away from the final answer.
Historically, there were several attempts to interpret the RL/RQ dichotomy in different terms. Firstly, the geometrical scheme was proposed for this in Reference [7] similarly to the "unification First X-ray observations of 3C 111 by ASCA [32] detected weak Fe Kα line. Later Suzaku and XMM-Newton observations [33] disclosed emission/absorption variability and weak reflection features. They also detected the Fe K complex. The jet component was probably dominating in the hard X-ray emission during their observable periods, as the high-energy exponential cut-off was not found in the Suzaku/HETG spectrum.
The spectral analysis of the Suzaku/XIS data was performed also in Reference [34] who detected the variability of an ultra-fast outflow on ≈7 days timescale, an ionized Fe emission line in the first observation, and a blueshifted absorption line in the second one, when the flux is ≈30% higher. They stated that the location of this ionized material is constrained at ≈0.006 pc and it belongs to the accretion disc at ≈20-100 R g from the central black hole and a highly ionized and massive ultrafast outflow with velocity ≈0.1 c, respectively.
Also, in Reference [31] it was performed a detailed investigation of a connection between the accretion disk and events in the jet in 3C 111. They pointed out that the X-ray continuum flux and Fe K line intensity are strongly correlated with a time lag less than 90 days.
An attempt to distinguish between the thermal vs. non-thermal high energy emission of 3C 111 has been performed in References [25,35] relying on the 0.4-200 keV composite spectrum from Suzaku, INTEGRAL and Swift/BAT data, as well as Fermi-LAT γ-ray data and radio together with the infrared data. The results of the pexrav spectral model [36]: photon index Γ = 1.68 ± 0.03, high-energy cut-off E C = 227 +143 −67 keV, reflection parameter R = 0.7 ± 0.3; the Gaussian emission line near 6.4 keV was characterized by equivalent width EW = 85 ± 11 eV, that is, the thermal part was dominating in the X-ray spectrum during their observational period. Their radio to γ-ray spectral energy distribution was successfully fit to the single-zone synchrotron self-Compton (SSC) model.
Recently the combined Chandra/HETG and XMM-Newton/RGS observational data on 3C 111 were analyzed in Reference [37]. They has confirmed an excess of neutral absorption with respect to that estimated from 21 cm radio surveys of atomic HI in the Galaxy [38]; they also detected a weak (EW = 16 ± 10 eV) Fe Kα emission line at E = 6.4 keV.
Analysis of the data obtained by VLBA and IRAM within the MOJAVE, UMRAO and F-GAMMA programs revealed the outburst activity and parsec-scale structural variability [39]. Following the results of the mm-VLBI, the object 3C 111 passed in 2007 the flux density outburst stage [40]. On the parsec scales, the outburst is resolved into a number of plasma components forming an intriguing bent structure. The spectacular wide-band radio-IR-X-ray observations of the kpc-scale jet in 3C 111 [41] revealed important differences between the knot morphology in the radio, X-ray and near-IR bands and interpreted them as a sign of two-component synchrotron model. Recently polarimetric VLBA observations of 3C 111 were analyzed in Reference [42] and the propagating polarized features in 3C 111 were interpreted as the result of the interaction of shocks in the underlying flow. Recent analysis of GMVA 86 GHz data from Reference [43] disclosed remarkable radio outburst activity in 3C 111 in 2007.
In this paper we make use of all X-ray data on 3C 111 available publicly in SWIFT, INTEGRAL, XMM-Newton and Suzaku data archives by the date 5 March 2019. To fix the parameters of the synchrotron emission of the jet base in 3C 111 we also use the Planck spectrum within the wavelength range of 10 7 -10 8 Å. The paper structure is as follows: In Section 2 we describe the X-ray observational data used in our work and its reduction. In Section 3 we describe two methods we propose to distinguish between the spectral components induced by jet base and accretion disk/corona. In the Section 4 we present the results of spectral fitting. Finally, in the Sections 5 and 6 we discuss our results and draw out the conclusions.

Observations and Data Reduction
In our analysis we used the XMM-Newton/EPIC, Suzaku/XIS, SWIFT (BAT and XRT), and INTEGRAL (ISGRI) datasets ( Figure 1). In this section we describe the reduction performed to the initial data to obtain the spectra and lightcurves. All the data periods are shown on the Figure 1.

Planck
The Planck source spectrum of 3C111 used here is based on The Early Release Compact Source Catalog (ERCSC) of compact sources derived from the data acquired by Planck during its first whole sky coverage, between 13 August 2009 and 6 June 2010. The Planck spectrum in the range 24-240 GHz is publicly available on the HEAVENS webpage 1 .

XMM-Newton
3C 111 was being observed by the XMM-Newton mission three times. Table A1 shows LOG of the XMM-Newton observation data concerning this object, including the exposures, count rates and signal-to-noise rations during all the periods of the observations. First XMM observation of 3C 111, performed in 2001 (OBSID 0065940101) was reduced separately to obtain the complex EPIC spectrum. The resulting files of the next two, being performed closely in time (both in February, 2009), were merged together to give the one common spectrum for two. During the time of all the three observations all EPIC cameras were operated in partial-covering mode. The EPIC data were processed using the standard software package XMM SAS ver. 14.0. The standard SAS chains epproc and emproc were applied for primary data reduction. The single-and double-photon events were taken into account (i.e., the PATTERN ≤ 4 option was applied). To exclude bad pixels and near-CCD-egde events from our consideration, the filter FLAG = 0 was also applied. The source counts were extracted from the 30 sec-radii circular regions around the source centre, and to extract the background counts were chosen the empty regions on the same CCD chip with 50-70 sec-radii. From the cleaned event files for spectra the periods of time of comparably high flaring activity were removed using the tabgtigen task. To construct the ancillary files and response matrices, we used the standard chains arfgen and rmfgen.

Swift
The Swift/XRT spectra were obtained using the online software provided by Dept. of Physics and Astronomy of the University of Leichester (XROA) 2 . We used single-pass centroid with the maximum of 10 attempts and 6 search radius. Event collections of several observations performed within the same months, were merged together (for observational windows 00036367002 and 00036367003, as well as for 00032659002, 00032659005 and 00032659006, and for 00032659008, 00032659009, and 00032659010).
The Swift/BAT event and survey data were processed to obtain the spectra in 14-195 keV energy range for every available period of SWIFT/BAT observations of 3C 111 using the BAT subpackage of the NASA HeaSoft version 6.27.2 software for astronomical data processing and analysis 3 .
All the event data were processed using the batbinevt, batupdatephakw and batphasyserr procedures to create the spectra and generate the error bars (the whole recipe can be found at the UK Swift Science Data Centre web page 4 ). The corresponding responses were generated using the batdrmgen procedure of the same subpackage. To create spectra from detector plane histograms of the BAT survey data the batdph2pha was applied 5 . Table A2 shows the Swift/BAT and XRT observations LOG. The data for the observational windows merged are shown together in the same table cell. Total observational time exposures and count rates, and signal-to-noise rations as well are shown also for these merged periods of the observations.

Suzaku
Suzaku/XIS spectra were reduced through online multi-mission web-interface UDON2 (Universe via Darts ON-line) for quick analysis and processing of X-ray astronomical satellite data archived in DARTS (Data Archives and Transmission System) provided by Institute of Space and Astronautical Science (ISAS) by JAXA (Japan Aerospace Exploration Agency). UDON2 and archived data of Suzaku observations can be found publicly at the corresponding web page 6 . After the images and spectra processing by UDON2, the resulting spectra obtained by all four XIS spectrometers during the same observational period were merged into one using the add_asca_spec routine included in the HEASOFT 6.25 X-ray data analysis and processing off-line software package with corresponding response matrices and ancillary files. Table A3 shows the Suzaku/XIS observations LOG, including the exposures, count rates and signal-to-noise rations during all the periods of the observations.

INTEGRAL
The ISGRI dataset of the INTEGRAL observations analyzed in the present paper includes all the observational data within the period 2003-2018, that is, spacecraft revolutions from 0055 to 1924. The total ISGRI exposure time of the dataset which had been used is 2.9 Ms (including all the observations when the object was at the angle less than 10 • off-axis). We performed INTEGRAL IBIS/ISGRI data analysis with version 10.1 of the Off-line Scientific Analysis software (OSA). We used standard recipes of spectral extraction for IBIS/ISGRI and OSA software. All the spectra were extracted individually for every science window and then summed up into the whole set. The source is detected by ISGRI up to ∼150 keV. The INTEGRAL/ISGRI observations LOG is shown in Table A4, including the exposures, count rates and detection significance during all the periods of the observations.

Two Ways to Distinguish between Nuclear and Jet Base Spectral Components
Jet base (of a percent of pc in size) is composed mainly of ultrarelativistic electron-positron (hereafter electron) plasma ejected at relativistic bulk velocity from the AGN central region. Being accelerated at shock fronts inside the jet, these electrons emit synchrotron radiation visible at a wide range of wavelengths from radio to ulraviolet. These accelerated electrons/positrons in the most cases are distributed over the energies following power-law dependency N(E) ∝ E −γ e [44]. Such distribution of charged particles generates self-absorbed synchrotron emission with spectral distribution: where ν 0 > ν is the break frequency at which an optical depth is τ = 1.
If the source is transparent enough or if we consider the frequency range above the break frequency of our object, ν >> ν 0 , this formula can be approximated with the simple power-law dependency with the slope α = (γ e − 1)/2 = Γ + 1, where Γ is the photon index [45][46][47]: At frequencies well below 1 GHz free-free synchrotron self-absorption flattens the overall spectrum. For 3C 111 in our situation we have 24-240 GHz Planck working range and thus we can use the approximate power-law Formula (2). As it was shown by Reference [47], the jet X-ray self-Comptonized synchrotron (SSC) or Inverse Compton (IC) spectra have the same photon index as the synchrotron emission of the jet at radio wavelengths. Thus one can apply the same power-law photon index value to fit the X-ray SSC/IC and radio synchrotron emission spectra of the jet base. Taking into account that at radio frequencies the jet emission of an RL AGN is many times more intensive than the emission of other AGN components their contribution to the AGN radio spectrum is also negligible in comparison with the jet one. This makes the radio band more promising for estimation of the parameters of the synchrotron emission of the jet base. The following spectral model of 3C 111 jet emission from radio to X-rays was shown in Reference [25]; here we show their one zone SSC model with radio (Planck) and X-rays (Swift/XRT, XMM-Newton/EPIC and Suzaku/XIS) working ranges to illustrate the method we use here to determine the jet spectral counterparts ( Figure 2). The radio-to-X-rays multiwavelength spectral energy distribution of 3C 111 within the one-zone SSC model [25], shown by black solid line. The frequency range of Planck is shown by blue zone and X-rays range corresponding to Swift/XRT, XMM-Newton/EPIC and Suzaku/XIS is shown by red.
In this view we propose to separate the contribution of the jet base and disk/corona based on the equality of synchrotron and SSC/IC slopes. Namely, we use for the jet base component of our X-ray spectral model the frozen values of the photon index Γ f r that is equal to that obtained from the Planck spectrum of 3C 111 in the range of 24-240 GHz. The Planck radio spectrum was fitted with a simple power-law model with the photon index Γ f r = 1.66 ± 0.06. Anyway, we here make the hypothesis about the slope and show that it allows us to reduce the tension concerning the values of the high energy cutoff.
The second way we discuss here, is an use some additional information obtained from in Fe K line emission. Fe K lines emission formed in the accretion disk of an AGN is correlated with the primary nuclear continuum emitted by accretion disk/corona. Here we presume a linear dependence between Fe-Kα line and primary nuclear emission fluxes.
In case of better quality data this information can serve as the basis for an alternative approach to separate the nuclear flux from the jet contributions, which presupposes that some relativistically broadened lines be present in the AGN spectrum and their equivalent width is variable.

Fitting the Spectra and Separation of the Jet Base Contribution
The model of 3C 111 3.0-300 keV X-ray spectrum consists of such main components: We perform our fitting in the two following steps: • fitting the high-energy continuum excluding the 5-7 keV area where the Fe-K emission lines can be visible, with the model, including only the jet and nuclear continuum emission; • fitting the 5-7 keV area with the Gaussian emission line and, in the case of presence of such, absorption line also (on this step we also use the jet+nuclear continuum model with parameters frozen to the values obtained at the first step).
Keeping in mind the spin paradigm for the radio-loud AGN 3C 111, we assume that the X-ray power-law tail above 80 keV of this object is mainly due to the jet base F jet . In this connection, we carried out the power-law fit of the total spectrum for ISGRI and BAT spectra above 80 keV to get the photon index Γ f it ; the resulting values are shown in Table 1 (in the bottom line of each cell of the second column 7 ). These values are mainly consistent within error limits with other calculations (except for observation period 4, where the match is marginal).
The proper and Galactic absorption N H0 was modelled using tbabs model. The Galactic absorption was included in the model as the lowest possible value of 7.4 × 10 21 cm −2 , following Reference [37].
Additional intercalibration constants were used in the simultaneous modelling of XMM-Newton/EPIC, Suzaku/XIS or SWIFT/XRT with INTEGRAL/ISGRI or Swift/BAT spectra and its value always was in the range 0.8-1.2. The errors, lower and upper limits indicated below correspond to a 90% confidence level (∆χ 2 = 2.71).
Thus, in our analysis, we fitted all the components listed above with the model: pexrav All the spectral fits were performed using XSpec v.12.13 of HEASOFT package, version 6.27. Swift/BAT and INTEGRAL spectra were treated within the observational time intervals overlapping with shorter intervals of XMM-Newton/EPIC, Suzaku/XIS or SWIFT/XRT observations; they were processed together to trace the variability details.
Observational data are divided into 12 observational periods described in Table A5. The spectra and models are shown in the Figures A1-A12. The initial datasets involve the combined soft 7 For the first period we have not obtained a reliable value.
(XMM-Newton/EPIC, SUZAKU/XIS, or SWIFT/XRT) X-ray spectra and hard (INTEGRAL/ISGRI or SWIFT/BAT) ones collected during the same month. The intercalibration constants between these pairs of spectra are also shown. Table 1. Hard X-ray continuum parameters with marginal errors (90% probability) for 3C 111. First lines in each cell include nuclear continuum parameters estimated with the frozen value of the jet base emission index Γ f r :index Γ nucl , exponential cutoff energy E c and reflection coefficient R of the nuclear emission. The second lines present analogous fit with frozen Γ f r and frozen reflection coefficient R f r = 0.7 as it was obtained by Reference de Jong et al. [25]. The lower lines of the first column show the jet emission indexes Γ f it (free) from the independent power-law fit of the spectral tail for E > 80 keV. First column: observational period; second to fifths columns: photon index, high-energy cut-off and reflection of primary nuclear emission; fifth column: absorbing column density excess in comparison with the Galactic absorption; sixth column: chi-squared/degrees-of-freedom.  Table 1 shows fitting parameters of the continuum spectra of 3C 111: the photon index Γ nucl , exponential high-energy cut-off E C of the primary nuclear continuum emission, its reflection parameter R, the absorption excess N H and discrepancies χ 2 /d.o. f . of the spectral fits. Here and in Tables below,  the observation periods are the same as in Table A5. The fitting has been performed with frozen value Γ f r = 1.66 of the jet base power-law SSC/IC emission. For comparison, the first column of Table 1 shows (in the bottom line) the indexes Γ f it derived independently from the power-law fit of the high-energy tail. The results of different fittings are in satisfactory agreement within errors.
Additionally, taking into account the results of References [25,35], we performed also the fitting with the frozen value of R = R f r = 0.7 as it was shown there. As a result we can see that there are periods of higher and lower values of the reflection coefficient within the error levels (i.e., periods 1, 5-8). We can explain these variations of the reflection fraction in terms of the changes of the inner radius of the accretion disc (or its radiatively efficient zone). Table 2 shows fitting parameters of emission spectral lines near 6 keV. In the second and third columns of this Table we show the line energy and width, and in the last one we show the equivalent widths calculated using the XSPEC12 standard procedure eqwidth. The eqwidth command calculates the equivalent width as a relation between the integrated photon flux produced by an additive model component (including in our case the line, jet and nuclear continuum), at the location of the peak of the photon spectrum (i.e., near 6.4 keV), and the flux (in units of photons per keV) at that energy of the continuum (i.e., jet + nuclear continuum in our case). The equivalent width is then defined in units of keV. The continuum (jet + nuclear) model is integrated from E − ∆E to E + ∆E. The initial value of ∆E = 0.3 keV was used. Table A6 shows the absorption line parameters, which were detected during 3, 5 and 12 periods. We also estimated the jet and primary nuclear component fluxes and the Fe-Kα emission line as well basing on this model; these fluxes are shown in Table 3.

Nuclear Emission vs. Fe-Kα Flux
In the previous subsection we have calculated the flux levels of the jet and primary nuclear counterparts in the spectra, basing on the spectral fitting model with the frozen value of the photon index of the jet counterpart. Here we propose another method to estimate the jet and nuclear fluxes using the results of the Fe-Kα emission line fitting. The advantage of this second method is that its result depends on the emission line parameters only and not connected directly with the continuum parameters (except the flux in 5-7 keV energy range), such as photon indices of jet or primary nuclear emission. This advantage makes this method independent from the radio spectral fitting results (i.e., jet photon index).
For most measurement periods it is rather difficult to isolate Fe Kα against the noise background and this is the reason of large errors of F line 5−7 in Table 3 (the index 5-7 shows the energy range in keV); see also Figure 3 representing the equivalent widths of Fe Kα line near 6.4 keV. Nevertheless, one can use these data in order to test some prospects for future use of F line 5−7 for separation of F jet and F nucl provided that more accurate data be available.
In Figure 3 we arranged equivalent widths of Fe Kα line vs. F jet 5−7 in order to estimate the dependence on the jet base continuum emission levels. If we remove the unreliable data with large error bars from Figure 3, we clearly see that the larger equivalent widths correspond to lowest F jet 5−7 . The highest value of equivalent width derived is about 150 eV (3-rd period, the left point in Figure 3 corresponding to zero jet base emission), which is compatible with the typical values of EW of RQ AGN [48]. This can be considered in favour of the similarity of the "central engine" of RQ and RL AGN. The reason leading to lower observed EW values in RL AGN than in RQ is the presence of an additional continuum component due to the jet base emission (cf. References [25,49]). Now we consider a relation between the X-ray continuum in 5-7 keV range and flux in the Fe Kα line. Reference [31] found out that the X-ray continuum flux and maximum flux in the Fe Kα line are correlated with the time lag consistent with zero. For a rough estimate, we suppose that the total (energy-integrated) observed flux F line 5−7 in the Fe Kα line is approximately proportional to the total nucleus (disk/corona) flux in the 5-7 keV range Then for the total continuum flux in the 5-7 keV range we have First of all, we can use the results of Section 4.1 in order to derive F line 5−7 , F jet 5−7 and κ, and then to check whether κ is the same for different observation periods. If there are periods of zero-compatible variations of the jet emission, then the alternative estimation of κ is possible. Namely, if one can assume for adjacent observation periods ∆F jet ≈ 0, so we could derive κ from variations of the fluxes to compare with the previous results of Section 1.
In case of X-ray 3C 111 X-ray spectrum, these 5-7 keV fluxes for the three model components (power-law jet emission, primary nuclear (disk/corona) emission and Fe Kα 6.4 keV emission line) are shown in Table 3. We have 3-rd period of zero-compatible jet emission with fairly good data on F nucl 5−7 and F line 5−7 , which can be used for determination of k yielding the value k = 33 ± 3. We used this value of k to estimate primary nuclear fluxes F nucl 5−7 from values of F line 5−7 using Formula (3) (Figure 4, lower panel). Despite the large errors one can see that the measured values F nucl 5−7 fall into the interval determined the error bars of estimated F nucl 5−7 , and the correlation coefficient between them is 0.85. This gives reason to hope that Formulas (3) and (4) can be used for independent determination of the nuclear flux in future measurements with more accurate data on the Fe K emission lines.  Table 3 (blue balls, "fitted model").

Overview of Our Results
We had shown that there is a possibility to separate the primary nuclear and jet base contributions in the hard X-ray spectra of RL AGN 3C 111 using the observational data obtained by INTEGRAL, SWIFT and Suzaku cosmic missions during the 12 periods of observations in common even taking into account that the X-ray and radio data by Planck mission were not synchronized in time.
In order to fulfill this separation we assumed the same photon index for the X-ray component of the jet base emission and for the radio emission spectra of this object (note that the radio emission is pure jet base synchrotron and the X-ray jet base counterpart is considered as the IC or SSC emission with the same photon index as the synchrotron one).
This method shows quite accurate results and can be recommended for further usage in the other objects, but we should concede that the Planck radio spectrum is averaged upon the time and the synchrotron photon index can vary with time as well. Therefore, it would be better to have the time rows of simultaneous observations in radio-and X-rays to establish more accurate dependence between these spectra and their parameters.
We also used possible correlation between the emission levels in primary nuclear continuum and Fe Kα emission lines formed in the near vicinity of the central black hole of an AGN [25,49]. The errors of the emission line parameters are quite large in most observation periods and thus give very rough estimates of the nuclear flux levels. In the same time, both approaches gave compatible results within the error bars. This gives reason to hope for the successful application of both approaches in the future for other objects, keeping in mind the necessity to have more precise spectral data for determining emission line parameters and synchronized time rows of radio-and X-ray observations. The results of Table 1 for hard X-ray spectrum of 3C 111 show that, for the periods having modest error limits, the spectral parameters of the primary nuclear continuum (i.e., photon index and high-energy exponential cut-off) mostly do not contradict the "spin paradigm" predictions. The exception is the 6-th period, when we have the high-energy exponential cut-off of primary nuclear continuum more than ∼130 keV for the fitting with free reflection R; on the other hand the result on this period are not reliable and also we have a low E c derived with fixed R f r = 0.7 that seems to be more physical. Therefore, 3C 111 deserves further observation and analysis with better data set in the high energy region.
In the same time we have shown that the maximum value of the Fe Kα line equivalent width obtained for the phase of the low level jet base emission is in agreement with the values typical for RQ AGN suggesting similarities of the "central engines" in RQ and RL AGN.

3C 111 in γ-Rays and Further Perspectives
The emission of RL AGN above 100 keV in 3C 111 was observed at higher energies by Fermi-LAT mission [50][51][52]. The simple power-law approximation of the spectrum has been performed by Reference [50] in Fermi-LAT 10 8 -10 11 eV energy range yielding the photon index Γ = 2.7 ± 0.2. This differs significantly from the value obtained from the Planck spectrum of 3C 111; the explanation of this difference is due to different emission mechanism in γ-rays. Since then, the Fermi-LAT observational exposure of 3C 111 has increased significantly and it makes sense to analyze the extended dataset in a more detailed way in future. The higher energy ranges brings the key information about AGN jets composition as the two best known scenarios of jet formation (hadronic and leptonic) leads to the significantly different spectral shape above 10 10 -10 11 eV [53]. Leptonic scenario leads to the exponential fall of the spectrum, whereas the hadronic component produces more complicated spectral shape due to secondary pion decay. Tracing out this difference can be performed more effectively by combining the Fermi-LAT and Cherenkov Telescope Array (CTA) observational data [54].

1.
We separated the "jet base" and "central engine" spectral components of AGN 3C 111 in the range from X-ray to soft γ-ray spectrum. For this purpose, we used the fixed value of the photon index from Planck radio data. We performed some tests showing that our main results are reliable.
In particular, we tested the results under supposition of the linear relationship between Fe Kα emission line and the 5-7 keV fluxes from the "central engine". 2.
The value of the high-energy exponential cut-off E c in the nuclear spectrum of the "central engine", obtained after the separation of the "jet base" spectral components from the total 3C 111 spectra, are mostly below 100 keV or compatible with this within the error limits. This means that spectrum of 3C 111 "central engine" most likely does not contradict the predictions of the spin-or gap-paradigm [15,16] Funding: This work has been supported in part by scientific program "Astronomy and space physics" (project 19BF023-01) of Taras Shevchenko National University of Kyiv.             During this period jet base emission was dominating over nuclear one (the jet base flux in 5-7 keV range was near as twice as higher than the nuclear one); emission in Fe-Kα line was weak; also the tiny signs of the weak absorption line near 5.9 were present.