The TeV-emitting radio galaxy 3C 264. VLBI kinematics and SED modeling

In March 2018, Mukherjee (2018) reported the detection by VERITAS of very-high-energy emission (VHE;>100 GeV) from 3C 264. This is the sixth, and second most distant, radio galaxy ever detected in the TeV regime. In this article we present a radio and X-ray analysis of the jet in 3C 264. We determine the main physical parameters of the parsec-scale flow and explore the implications of the inferred kinematic structure for radiative models of this gamma-ray emitting jet. The radio data set is comprised of VLBI observations at 15 GHz from the MOJAVE program, and cover a time period of ~2 years. Through a segmented wavelet decomposition method we estimate the apparent displacement of individual plasma features; we then perform a pixel-based analysis of the stacked image to determine the jet shape. The X-ray data set includes all available observations from the Chandra, XMM, and Swift satellites, and is used, together with archival data in the other bands, to build the SED. Proper motion is mostly detected along the edges of the flow, which appears strongly limb-brightened. The apparent speeds increase as a function of distance from the core up to a maximum of ~$11.5$ c. This constrains the jet viewing angle to assume relatively small values ($\theta\lesssim10^{\circ}$). In the acceleration region, extending up to a de-projected distance of ~$4.8\times10^4$ Schwarzschild radii (~$11$ pc), the jet is collimating ($r\propto z^{0.40\pm 0.04}$), as predicted for a magnetically-driven plasma flow. By assuming that the core region is indeed magnetically dominated ($U_B/U_e>1$), the SED and the jet power can be well reproduced in the framework of leptonic models, provided that the high-energy component is associated to a second emitting region. The possibility that this region is located at the end of the acceleration zone, either in the jet layer or in the spine, is explored in the modeling.


Introduction
Over the past 10 years, the gamma-ray sky has been thoroughly probed by the Fermi-LAT space telescope, and active galactic nuclei (AGN) can now be firmly defined as the dominant class among γ-ray emitters . In most cases, the γ-ray detections have been associated to jetted-AGN seen at small angles, i.e., to blazars, while only ∼2% of the LATdetected AGN are classified as radio galaxies. The low detection rate of radio galaxies is due to the relativistic nature of plasma flows in AGN: the high-energy emission is preferentially detected when amplified by relativistic Doppler boosting, i.e., when the narrow boosting cone intercepts our line-of-sight. Radio galaxies which are detected at very high energies (V HE; > 100 GeV) with Cherenkov telescopes are, as expected, even rarer. To date, the class of TeV-emitting AGN comprises only 6 misaligned jets, 1 : Centaurus A, M 87, 3C 84, IC 310, PKS 0625-35, and since March 2018 (Mukherjee 2018), 3C 264, subject of the present study. These nearby galaxies all share radio morphologies and energetics typical of Fanaroff-Riley I sources 1 http://tevcat.uchicago.edu/ (FRIs), although IC 310 and PKS 0625-35 may be rather classified as intermediate FRI-BL Lac sources (Kadler et al. 2012;H.E.S.S. Collaboration et al. 2018). The newly detected 3C 264 is at present the second most distant TeV-emitting radio galaxy, being located at a redshift z = 0.0217.
Modeling the high-energy emission from radio galaxies, while highly instructive for testing the unified scheme, as well as our understanding of the γ-ray emission mechanism in AGN, has shown to be challenging. One zone leptonic models, the synchrotron self-Compton (SSC) (Maraschi et al. 1992;Bloom & Marscher 1996) and the external Compton (EC) (Dermer & Schlickeiser 1993;Sikora et al. 1994;Błażejowski et al. 2000;Celotti et al. 2001), are often adequate to reproduce the broadband spectral energy distribution (SED) in blazars (e.g., Böttcher et al. 2013), but appear insufficient for TeV BL Lacs and for radio galaxies (e.g., Georganopoulos et al. 2005;Błażejowski et al. 2005;Brown & Adams 2011;Petropoulou et al. 2014). One of the main difficulties is that the Lorentz factors Γ and Doppler factors δ inferred from radio studies for these two classes are usually low (Γ, δ 5) (Marscher 1999;Piner & Edwards 2004;Piner et al. 2010;Lister et al. 2016;Jorstad et al. 2017), too low Article number, page 1 of 12 arXiv:1905.06634v1 [astro-ph.HE] 16 May 2019 A&A proofs: manuscript no. main to account for the high-energy component of the SEDs. Several alternative models have been proposed that relax the assumption of a one-zone, homogeneous emission region (see, e.g., Rieger & Levinson 2018, for a recent review). For instance, it has been shown that if velocity gradients exist in the flow, either along the jet axis (see, e.g., the decelerating jet model in Georganopoulos & Kazanas 2003) or transverse to it (see, e.g., the spine-sheath model in Ghisellini et al. 2005), the inverse Compton emission is enhanced by the interaction between the different zones. The existence of transverse velocity gradients is supported by observations, e.g., by VLBI studies of the nearby jets in Mrk 501 (Giroletti et al. 2004), 3C 84 (Nagai et al. 2014;Kim et al. 2019;Giovannini et al. 2018), M 87 Kim et al. 2018) and Cygnus A (Boccardi et al. 2016), and the spine-sheath scenario has been successfully applied to reproduce the TeV emission in 3C 84 (Tavecchio & Ghisellini 2014). Deceleration of FRI jets on sub-kiloparsec and kiloparsec scales is likewise supported by abundant observational evidence (e.g., Hardcastle et al. 1997;Feretti et al. 1999;Laing & Bridle 2002), although its role in the TeV production may be challenged by the fast variability that frequently characterize this emission, implying it originates in very compact regions likely located in the innermost jet.
In addition to the multi-zone leptonic models, hadronic and lepto-hadronic models (e.g., Mannheim & Biermann 1992;Dar & Laor 1997;Böttcher 2005;Mastichiadis et al. 2013), which attribute the high-energy SED component to proton-synchrotron or to proton-photon interactions, provide valid alternative solutions for reproducing the TeV emission in blazars and in radio galaxies (Aharonian 2000;Mücke et al. 2003;Aharonian et al. 2005;Protheroe et al. 2003;Petropoulou et al. 2014). The hypothesis that relativistic protons also exist in AGN jets has been recently corroborated by the identification of several BL Lacs as possible sources of the neutrinos detected by IceCube (Padovani & Resconi 2014;Kadler et al. 2016;IceCube Collaboration et al. 2018), since neutrinos are a natural secondary product of protonphoton interactions.
In this article we investigate the radio and X-ray properties of the jet in 3C 264, and we use the inferred parameters to model the broadband SED. This radio source is well studied on subkiloparsec and larger scales, but the innermost regions of its jet, where the most energetic processes are likely to take place, are still largely unexplored. Through multi-epoch VLBI observations, we examine the kinematic properties, the internal structure and the geometry of the plasma flow, which are fundamental parameters in the investigation of the high-energy emission mechanisms. The paper is organized as follows. In Section 2 we summarize previous studies of the source. In Section 3 we present the radio and X-ray data set. In Section 4 we describe the methods adopted in the analysis of the VLBI images and we discuss the results. In Section 5 we model the broadband SED. Throughout the article we assume a ΛCDM cosmology with H 0 = 70.5 h −1 km s −1 Mpc −1 , Ω M = 0.27, Ω Λ = 0.73 (Hinshaw et al. 2009).

The radio galaxy 3C 264
3C 264 is a nearby (D L = 94 Mpc) Fanaroff-Riley I radio galaxy whose jet has been extensively studied in the radio, optical, and X-ray bands (Lara et al. 1997;Baum et al. 1997;Lara et al. 2004;Perlman et al. 2010;Meyer et al. 2015, and references therein). The active source is hosted by the elliptical galaxy NGC 3862, located off-center in the rich cluster Abell 1367, and harbouring a supermassive black hole with mass M BH ∼ 4.7 × 10 8 M . This mass was derived by Balmaverde et al. (2008) from the stellar velocity dispersion using the relation by Tremaine et al. (2002). On kilo-parsec scales, 3C 264 presents a complex head-tailed radio structure embedded in a vast region of low surface brightness emission, which may be indicative of a dense environment. The jet and the weak counter-jet emanate from a compact core, and extend respectively in the north-east and south-west direction. In the inner 2.2 arcseconds, the approaching jet is also shining at optical and X-ray energies via the synchrotron mechanism (Perlman et al. 2010). The analysis presented by Perlman et al. (2010) also showed that the radio and the optical jet are highly polarized on sub-kiloparsec scales, with some regions reaching fractional polarization of the order of 50%. The magnetic field vectors are remarkably well aligned with the jet axis, both at the jet boundary and at the center.
The jet orientation with respect to the observer is uncertain. Hubble Space Telescope imaging reveals a jet punctuated by few bright knots, apparently emanating from a circular optical ring of unclear nature (see the discussion below). The fastest optical knot was observed to propagate with an apparent speed β app ∼7 c, and to¨catch up¨with a slower moving knot (β app ∼1.8 c) located farther downstream, and possible site of particles acceleration (Meyer et al. 2015). The kinematic measurement of the maximum apparent speed constrains the jet viewing angle θ to assume relatively small values (θ 16 • ). Based on the analysis of the longitudinal brightness profiles and of the jet-to-counter jet ratio, previous radio studies had determined a rather different jet orientation of ∼50 • (Baum et al. 1997;Lara et al. 2004), which is however difficult to reconcile with the remarkably circular appearance of the optical ring, seen at a projected radius of ∼300 − 400 pc. In the simplest scenario, the ring marks the boundary of a flattened dust disk seen face-on, as observed in several FRI galaxies with optical jets (Sparks et al. 2000); therefore the jet axis, which is most likely perpendicular to the disk, should not be oriented at a very large angle. More complex interpretations of the ring (a puffed-up disk or a spherical bubble) have been proposed by Baum et al. (1997) and Lara et al. (2004) in view of the fact that its location approximately coincides with the region where the jet starts deflecting towards the north, expands faster and becomes dimmer. These elements would suggest an actual physical interaction between the ring and the jet, and thus an alternative geometry.

Radio data
The radio data set comprises 7 VLBA (Very Long Baseline Array) observations at 15 GHz (2 cm) performed between September 2016 and October 2018 within the MOJAVE survey (Lister et al. 2018). The calibrated data were downloaded from the MO-JAVE public archive, and the CLEAN images were produced in DIFMAP (Differential Mapping, Shepherd et al. 1994) through the CLEAN algorithm after minor editing. In Table 1 we report the log of observations and the characteristics of the seven maps. A stacked image, created after convolving each map with a common beam of 0.85 × 0.30 mas, 25 • , is shown in Fig. 1. The beam position angle was set to coincide with the mean jet direction in order to more clearly show the internal structure of the flow, which is characterized by a prominent limb-brightening. The jet extends approximately for ∼8 mas, corresponding, at the source redshift and for the adopted cosmology, to ∼3.5 projected parsecs (1 mas∼0.44 parsec). The faint jet emission appears smooth, with a single, main brightness enhancement (knot) at a radial distance of ∼6 mas. The radio core emission is moderately variable 3.2. X-ray data 3C 264 was observed by the Chandra, XMM-Newton and Swift X-ray satellites (see Table 2). The XMM-Newton observation of 2001 was analyzed by Donato et al. (2004) while the Chandra observation of 2004 was presented by Evans et al. (2006). In both analyses the source was characterized by a steep power law photon index Γ∼2.4 and a luminosity of ∼ 10 42 erg s −1 between 2-10 keV. Thanks to an accurate analysis of the Chandra image from the 2004 observations, the X-ray counterpart of the radiooptical jet was later revealed by Perlman et al. (2010).
Here we present a comprehensive X-ray study of the nuclear region of 3C 264, collecting all the available data from the XMM-Newton, Chandra, and Swift/XRT archives. In three occasions, 3C 264 was observed more than 8' off-axis, being the X-ray centroid of the cluster Abell 1367 the primary target. On the contrary, Swift pointed directly at the source and performed 17 snapshots (less than 4 kiloseconds each) in 2018.
All Chandra observations were performed using CCD/ACIS-S. The data were reprocessed by using CIAO version 4.10 with calibration database CALDB version 4.7.8, and applying the standard procedures. The radius r of the extraction region depends on the distance of the source from the aim-points. We chose a circle of radius r = 5" when the source was on-axis and r = 27" for an off-axis position. This ensures the collection of at least 90% of counts in the selected regions 2 . The background spectra were extracted near the source in the same CCD. Data were then grouped to a minimum of 15 counts per bin over the energy range 0.3 − 7 keV.
We checked if 3C 264 was affected by severe pile-up in 2000 (when it was the main target of the pointing). As a 1/8 sub-array configuration was adopted during the observation (time frame 0.41 s), we estimated 3 a negligible (5%) pile-up. However, as a further check, we also extracted a spectrum from an annulus to exclude the inner region where a higher signal degradation is expected. The two spectra (circle and annulus) were fitted with the same model. The spectral parameters were completely consistent, although less statistically constrained in the annulus case.
In the two XMM-Newton observations, 3C 264 was off-axis and therefore no pile-up effect is expected. These data were analysed using SAS version 14.0 and the latest calibration files. Periods of high particle background were screened by computing light curves above 10 keV. In 2009, only the MOS1 spectrum could be analyzed. The pn camera was strongly affected by high background, while in MOS2 the source fell exactly in a CCD gap. Different radii were selected for the extraction circle of the source, depending on its position in the different cameras. In at least three cases, 3C 264 was near the edge of a CCD and it was necessary to assume small regions. Spectra from the first observation (2001) were obtained by extracting events in circular regions of r = 40" (pn), r = 19" (MOS1), and r = 22" (MOS2). For the only spectra available in 2009 a circle of radius r = 22" was adopted. Data were then grouped to a minimum of 25 counts per bin over the energy range 0.3 − 7 keV.
Swift/XRT data were reduced using the data product generator of the UK Swift Science Data Center 4 . Source spectra for each observation were extracted from a circular region with a radius which depends on the count rate of the source. No pile-up effect reduced the quality of the data. The background was taken from an annulus with an inner radius of 142" and outer radius of 260" centered on the source. Spectra were grouped in order to have at least 1 count per bin (see Evans et al. 2009, for details). All Swift observations were also combined together in order to A&A proofs: manuscript no. main α -Net exposure time. β -Flux in units of 10 −13 erg cm −2 s −1 in the 2-10 keV range. ζ -Luminosity in units of 10 41 erg s −1 in the 2-10 keV range. κ -pn, MOS1 and MOS2 are fitted together using the same model but different power law normalisation. Exposure time. flux, and luminosity refer to MOS1. -Column density in units of 10 20 cm −2 . †c indicates C-statistic. § -When a thermal emission is included in the model, the χ 2 value refers to the global fit. produce an average spectrum with higher signal to noise ratio. This could be grouped to a minimum of 15 counts in the 0.3-7 keV range allowing the application of the χ 2 -statistic. Spectra were then fitted with XSPEC12.10 (Arnaud et al. 1999). The XMM spectra of the 2001 observation were fitted together using the same model. Only the flux was allowed to vary in order to take into account possible mis-calibration among the instruments. For Chandra and XMM-Newton, a χ 2 -statistic could be applied, while for the single Swift observations the Cstat-statistic was generally preferred, considering the short exposure and the limited number of counts. However, in at least 4 observations data could be grouped to a minimum of 15 counts per bin and the χ 2 -statistic adopted.
We considered a power law model to fit the nuclear continuum and a thermal emission (APEC model) to reproduce the gas emission associated to the cluster Abell 1367. The Galactic column density was fixed to N H =1.82×10 20 cm −2 (Kalberla et al. 2005). The presence of an intrinsic N H,z , responsible for the attenuation of the nuclear emission, was also tested and required by all the Chandra and XMM-Newton observations. The N H,z intercepted by our line of sight is of a few 10 20 cm −2 , more probably related to the galaxy rather than to the nuclear region. Hot emission due to the gaseous cluster environment was detected by Chandra and Swift/XRT (considering the average spectrum). The measured gas temperatures (kT∼0.3) and luminosities (L 0.5−2∼keV ∼ 10 41 erg · s −1 ) are consistent within the uncertainties. The large point spread function of the MOS off-axis and the short XRT integration time of each single observation prevent the detection of any low surface brightness extended emission. The results are reported in Tables 2 and 3. The parameter uncertainties correspond to 90% confidence for a single parameter, i.e., ∆χ 2 = +2.7 (Avni 1976).
The light curve, shown in Figure 2 (upper panel), collects the 2 − 10 keV fluxes measured by all the satellites from 2000 to 2018. The Swift flux refers to an average value obtained from the collection of all the 2018 observations. It is evident that the source luminosity changed by about a factor 3 from 2000 to 2018. In the lower panel the photon indices measured in the same period are shown. The investigation of a possible correlation between flux and photon index could, in principle, provide further insights on the physical mechanisms driving the variability. For instance, a "hardening when brightening" spectral behaviour is often observed in jets (e.g., Ahnen et al. 2017), and could be expected as a consequence of the injection of fresh energetic particles in the flow. However, due to the bias introduced by integrating the spectrum to calculate the flux (see Massaro & Trevese 1996), the present data are insufficient to attest the statistical significance of a correlation between Γ and flux. This may be better explored in the future based on a denser monitoring of the source.  In order to statistically test the high state of the source in that period, all the observations performed in a single month were collected together. Each spectra could be grouped in order to have at least 20 counts per bin and the chi-statistic applied. A power law absorbed by Galactic column density was a good fit in all the three cases (see Table 2). The average January, March, April fluxes are shown in Figure 3 as red points, with error-bars corresponding to 1-σ (68% confidence level for two parameters of interest i.e. ∆χ 2 =2.3). A standard χ 2 test, applied to the average fluxes (red points in Figure  3), confirms the variability of 3C 264, being the probability to be constant less than 10 −3 . Interestingly, the X-ray enhancement was observed two days after the end of the VERITAS campaign that detected for the first time 3C 264 in the TeV band. The two bands could then be related to each other and the TeV photons produced by up-scattering of keV photons in the jet.

γ-ray data
In the GeV band, 3C 264 is associated with 3FGL J1145.1+1935, a source reported in the Third Fermi Large Area Telescope Source Catalog (Acero et al. 2015). It was detected at the 5.7σ significance level with a flux in the 1-100 GeV band of (2.6 ± 0.7) × 10 −10 ph · cm −2 · s −1 and a power law photon index Γ 3FGL =1.98±0.2. The GeV spectrum of 3C 264 is quite hard compared to the other radio galaxies with a γ-ray counterpart. Indeed, the source is also listed in the third Fermi-LAT catalog of hard-spectrum sources (Ajello et al. 2017), containing all objects detected in seven years of survey above 10 GeV. In the hard regime, its power law index is Γ >10GeV = 1.65 ± 0.33 and the flux between 10 GeV and 1 TeV F 3FHL = (2.5 ± 1.0) × 10 −11 ph · cm −2 · s −1 . No significant flux variation on time scales of one month is reported for this source in the 3FGL catalog.

Kinematics
The kinematic analysis of the radio jet was performed using the the publicly available Wavelet Image Segmentation and Evaluation (WISE) code (Mertens & Lobanov 2015. By exploiting a segmented wavelet decomposition method (SWD), the WISE code enables a multi-scale detection of statistically significant structural patterns (SSP) in the examined astronomical image. Given a multi-epoch set of images, the significant patterns are then tracked in time by running a multi-scale cross-correlation algorithm (MCC) between pairs of consecutive epochs, which yields a measurement of the pattern displacement.  The chosen approach is ideal for the kinematic study of the jet in 3C 264, which is transversely resolved and lacks bright and well defined knots. Traditional methods, like the MODELFIT subroutine in DIFMAP (Shepherd et al. 1994), rely on a priori assumptions concerning the shape of the structures contributing to the brightness distribution. For instance, in the case of jets, a set of two-dimensional Gaussian components is typically fitted to the data. The wavelet decomposition, on the other hand, enables arbitrarily shaped structures to be identified.
The maps considered in the SWD analysis were slightly super-resolved in the direction perpendicular to the jet axis, being restored with a common beam of 0.85 × 0.30 mas, 25 • (the same used to create the stacked image in Fig.1). These parameters were chosen to facilitate the investigation of the jet kinematic structure in the transverse direction. To cross-check the results obtained, the analysis was also performed after restoring all images with the average natural beam, which is approximately 0.85 × 0.40 mas, 5 • .

Image decomposition
The segmented wavelet decomposition was carried out for each epoch on a set of five spatial scales (0.12 mas, 0.24 mas, 0.48 mas, 0.72 mas, 1.44 mas), plus three intermediate scales (0.18 mas, 0.36 mas, 0.96 mas). The intermediate wavelet decomposition (IWD), described in , is aimed at recovering structural information in optically thin, stratified flows, where multiple patterns can overlap without inducing full obscuration. The complex limb-brightened jet structure shown in Fig. 1 suggests that this is most likely the case in 3C 264.
In Figure 4 we present an example of the multi-scale wavelet decomposition of the jet in 3C 264, using data from July 2017. While at the largest spatial scale of 1.44 mas the algorithm detects only 2 components, the smaller scale decomposition obviously yields an increasingly rich set of information, with the detection of 18 components at the smallest spatial scale of 0.12 mas. In the specific observation presented in Fig. 4 the jet shows a quite complex structure, being composed by four filaments which appear to possibly intersect and cross each other. Such pattern is however not clearly seen in other maps. As shown in Fig. 5, a smaller number of filaments and/or features is detected in other epochs. This is most likely due to variations of the local signal to noise ratio with time, with only some portions of the jet being illuminated in a single epoch.

Cross-identification and proper motions
Once the significant patterns are identified, the multi-scale crosscorrelation algorithm matches them by examining pairs of adjacent epochs. The matching of features detected on the larger scales is used as a reference for the robust cross-identification of features on the smaller scales. In the MCC analysis, we assumed a correlation threshold of 0.6 and we enabled the displacements to vary between 0 and 12 mas/year in the longitudinal direction and between -6 and 6 mas/year transversely. About 59% of the total features detected at all scales was successfully crossidentified, with the normalized cross-correlation coefficient being systematically larger than 0.75. When considering only the results obtained at the smallest scale of 0.12 mas, the success rate of the matching algorithm decreases to 39%.
In Figure 6 we present the displacement vectors of the patterns that were successfully cross-identified in at least three epochs, overlaid on the stacked image. While we show the re- sults obtained at the finest scale, we remark that the displacements determined at larger scales are very consistent. The features are observed to move on complex trajectories, mainly following the jet direction but also showing a significant transverse component in some cases. In the inner jet, up to a distance of ∼ 3 mas, most of the features that could be matched belong to the eastern limb. In the outer jet, instead, the proper motion is detected on both limbs. In addition, thanks to the fact that the jet transverse structure is better resolved at larger distances, in the outer jet we can have a closer look at kinematic properties near the jet axis. The displacement of features detected in the central body of the jet is characterized by a larger transverse component. In particular, the jet limbs appear to clearly intersect and cross each other at a projected distance from the core of ∼ 4.5 mas, in the proximity of the brightest jet knot.
From the proper motions, we have then computed the apparent velocity β app of each feature in units of the speed of light c. At the source redshift and for the adopted cosmology, a proper motion of 1 mas/year corresponds to an apparent speed of 1.52 c. The uncertainties on the velocities are calculated as described in  by accounting for the positional error of each SSP, which in turn depends on the signal-to-noise ratio of the SSP detection. In Figure 7, we report the obtained values of β app as a function of radial distance from the jet core z. In this plot, the color assigned to each feature matches the color of the displacement vectors shown in Fig. 6. The velocities are distributed in a broad interval between ∼0.5 c and ∼11.5 c. The presence of a gap is apparent between a lower envelope of speeds reaching a maximum value of ∼5 c, and an upper envelope peaking at ∼11 c. Both the upper and the lower envelope show a trend of increasing velocity with distance up to ∼3.5 − 4.5 mas (or ∼1.5 − 2 projected parsecs) from the core, followed by a possible flattening at larger separation. Similar acceleration profiles have been recently obtained in the kinematic studies of M 87  and Cygnus A (Boccardi et al. 2016), in agreement with the theoretical prediction that magnetic acceleration of the bulk flow proceeds slowly and extends on parsec scales (e.g., Vlahakis & Königl 2004).  In Figure 8, we also examine how the apparent velocities are distributed in the transverse direction. To this aim, we plot the apparent speeds as a function of position angle of the features. The jet axis is oriented at an average position angle of ∼0.44 radiants. The SSP are detected in a range between ∼0.2 and ∼0.65 radians and, as already discussed above, they tend to cluster at the edges of this interval, i.e. on the jet limbs. The high velocities, above ∼5 c, are found either close to the jet axis or along the eastern limb. Therefore from our analysis we do not infer a clear speed stratification, e.g., of the spine-sheath kind. It is possible, however, that such a velocity gradient does exist in the flow, but the low and the high speed components appear mixed due to the overlap of optically thin filaments seen in projection. This scenario was proposed previously for the interpretation of the kinematic properties of the M 87 jet , where a spine-sheath structure is also not directly and unambiguously inferred from the kinematic study, even though the plasma flow appears prominently stratified.
As mentioned in Sect. 4.1, the SWD analysis was performed a second time by considering maps convolved with the average natural beam, which is approximately 0.85 × 0.40 mas, 5 • , with the aim of cross-checking the robustness of the kinematic results. While the algorithm cannot successfully detect and/or match all of the features in Fig. 7 due to the lower resolution in the transverse direction, both the displacement map and the range of apparent speeds inferred though this second approach are very consistent.

Jet orientation, Lorentz and Doppler factors
The detection of moderate-to-high superluminal speeds in the jet of 3C 264 is in agreement with previous optical kinematic studies of the source (Meyer et al. 2015(Meyer et al. , 2017, and constrains the jet viewing angle θ to assume relatively small values. Since we measure a maximum apparent speed of ∼11.5 c, the equation that relates the apparent to the intrinsic speed, limits the viewing angle to a maximum value θ max ∼ 10 • . On the other hand, given the absence of prominent radio flux density variability (see Table 1), and considering the quite extended large scale morphology of 3C 264, the viewing angle is unlikely to assume very small values. Therefore, in the following calculations we will assume the upper limit θ = 10 • . For this orientation, an apparent speed β app = 10 c, which is the average speed in the upper envelope in Fig. 7, implies an intrinsic speed β = 0.9978 and a Lorentz factor Γ = 1/ (1 − β 2 ) ∼ 15. For the lower speed envelope with β app = 5 c, instead, we obtain β = 0.9808 and Γ ∼ 5. The assumption of the above parameters finally enables us to calculate the Doppler factor δ of the two envelopes, which is a function of both the intrinsic speed and the orientation: For the high and the low speed values we obtain respectively δ ∼ 4 and δ ∼ 6. We note that if the slow envelope is actually reflecting the speed of the limbs and the fast one is associated to regions close to the jet axis, the difference between these Doppler factors naturally explains the observed limb-brightening. Moreover, the jet orientation is such that a large difference in Γ (Γ ∼ 15 versus Γ ∼ 5) translates into a moderate difference in δ, which enables the fast filaments to be also visible. In other words, the central filaments are not de-boosted (δ < 1), but are less boosted than the limbs.

Collimation profile
As a final step in the study of the innermost radio jet in 3C 264, we have examined its collimation properties. The analysis of the jet shape complements the kinematic study presented in Sect. De-projected distance z [R s ] d ∝ z 0.40 ± 0.04 Fig. 9. Collimation profile of the jet in 3C 264. The flow has a close-to-parabolic shape on parsec scales. A recollimation is observed to start at a distance of ∼11 pc (dashed vertical line). This location may coincide with the end of the acceleration zone, as seen in Fig. 7. A fit performed by excluding the recollimation region, i.e. in the interval 0 < z < 11 pc, yields d ∝ z 0.66±0.02 .
4.1, as it allows us to further test our finding that bulk acceleration is taking place on the examined VLBI scales. In fact, according to theoretical models, acceleration through MHD processes cannot take place at large distances from the black hole if the jet is freely expanding, i.e., it has a conical shape; therefore a jet which is still accelerating on parsec scales should undergo active collimation in the same regions (see, e.g., Beskin & Nokhrina 2006;Komissarov et al. 2007;Lyubarsky 2009).
With the aim of determining the jet width at each pixel and investigating its evolution with distance from the core, we have created a new stacked image after restoring all maps with a circular beam of 0.6 × 0.6 mas, corresponding to the average equivalent beam b eq = √ b min × b max . The choice of a slightly larger restoring beam in the transverse direction makes the analysis less sensitive to the the jet transverse stratification. The stacked image was sliced pixel by pixel (1px=0.06 mas) in the direction perpendicular to the the jet axis using the AIPS task SLICE. Through the task SLFIT, the transverse intensity distribution in each slice was then modeled by fitting a Gaussian profile, in order to extract its width (FWHM) and its peak flux density. Figure 9 shows the dependence of the de-convolved jet width d on the de-projected distance from the core z. The distance is now expressed in linear scale in units of parsecs and of Schwarzschild radii R S , and was de-projected assuming a viewing angle θ = 10 • , as determined in Sect. 4.2. In the plot we only report measurements obtained for transverse profiles with a signal-to-noise ratio higher than 10. The error bars on the width are conservatively assumed equal to one fifth of the convolved FWHM. The jet appears to expand quite smoothly up to a distance of ∼11 pc (∼4 projected mas), where a recollimation occurs. By fitting a power law of the form d = a · z b to the whole profile, we obtain d ∝ z 0.40±0.04 . Therefore, as expected for an accelerating jet, the flow is collimating and has a close-to-parabolic shape. This conclusion is still valid if, in the fit, we consider only the inner regions of the jet, excluding the recollimation at ∼11 pc. In this case, we obtain d ∝ z 0.66±0.02 . We note that the distance at which the recollimation occurs coincides with the distance, in Fig. 7, where the speeds profile starts flattening, and where the jet limbs intersect (Fig. 6). Ultimately, our results support a scenario in which the acceleration and collimation of the jet in 3C 264 take place in the inner ∼4 projected milliarcseconds, corresponding to a de-projected distance of ∼11 parsecs, or ∼4.8 × 10 4 R S (assuming M BH ∼ 4.7 × 10 8 M ).

SED modeling
In light of the properties inferred from the analysis of the radio and X-ray data, we modeled the nuclear SED of 3C 264 (Fig.  10). The radio to γ-ray data were collected from publicly available databases (NED 5 , ASDC 6 ). In addition to our estimated Xray fluxes, the average γ-ray flux reported in the 3FGL (lightgrey bow tie in Fig. 10) and the 3FHL data points are shown. While a detailed analysis of the source behavior in the GeV-TeV regime is beyond the scope of this work, for guidance in the modeling, we included in the SED an estimate of its VHE flux. We used the integrated (>300 GeV) flux, (1.3 ± 0.2)×10 −12 ph cm −2 s −1 (Mukherjee 2018), and assumed a power-law spectral shape with a hard and a steep photon index (Γ V HE =2.3 and 2.8, see e.g. the review by Rieger & Levinson 2018).
The SED displays the characteristic double-hump shape, with the low-energy component extending up to the X-rays and showing no evidence of features related to the accretion process (such as a bump in the UV regime and broad emission lines in the optical band, see also Buttiglione et al. 2009Buttiglione et al. , 2010. Such properties are in agreement with the expectations for a low-power AGN jet which, according to our estimate of the viewing angle, lies at the boundary between FR I and BL Lac sources. The fact that 3C 264 is not always significantly detected on monthly timescales by Fermi-LAT and the recent detection at TeV energies suggest that the source is a relatively faint γ-ray emitter with possible, low-amplitude variability (Rieger & Levinson 2018). As discussed in Sect. 3.2, flux variability is clearly observed at Table 4. Model parameters (see Sec. 5 in the text) and jet powers of the soft low state (SLS) and hard high state (HHS) of 3C 264: L rad is the total (synchrotron and SSC) radiative power, L B the Poynting flux, L e and L p the kinetic powers in emitting electrons and in protons (assuming one cold proton per electron), respectively, and L kin =L e +L p . X-rays, a property which is also evident at the respective range in the SED.

SLS
In the modeling presented in the following we explore a scenario in which the source switches between two states: a soft/low state (SLS) characterized by the lowest flux and softest X-ray spectrum, and a hard/high state (HHS) when the source is the brightest and has a hard X-ray spectrum. Based on the results obtained in our kinematic study, the modeling relies on two working hypotheses: 1. the SLS corresponds to the core's average SED and 2. the core is located at the base of the acceleration region, where the plasma still has a mildly relativistic bulk speed (Γ b ∼2).
In this framework, we will discuss different scenarios for the origin of the HHS. We assume that the emission is produced in a spherical blob of radius R, moving with a bulk Lorentz factor Γ b , and seen at a viewing angle θ. The electrons follow a broken power-law energy distribution: p 1 , p 2 are the indices below and above the spectral break, γ min , γ max and γ break the minimum, maximum Lorentz factors and the Lorentz factor at the spectral break, respectively. By interacting with the magnetic field (B), the electrons radiate via the synchrotron mechanism, which accounts for the low-energy component of the SED. The synchrotron photons are then up-scattered by the same electron population (Synchrotron-Self-Compton, SSC), giving rise to the highenergy hump.
We first tested a scenario in which the SLS and the HHS emission are both produced in the core (Model 1 - Figure 10, upper panel), but the HHS originates in a more compact and more magnetized region. The model parameters for the two states are shown in Table 4. The compact volume assumed for HHS implies that the emission can vary on timescales as short as few days (t var ∼R(1 + z)/(cδ)≈2.7 days), in agreement with those suggested by the Swift monitoring (Sect. 3.2, Fig. 3). In the HHS the blob dissipates almost 60% of its total power (L jet = L kin + L B ∼ 8 × 10 41 erg s −1 ) in the form of radiation. The model parameters of the SLS are rather similar, suggesting that the latter state could directly derive from the former. The SLS could be the result of strong radiative losses plus adiabatic losses due to the mild volume expansion of the same emitting region (the radius increases by a factor of 2 from the HHS to the SLS). However, the energetic budget of the two states disfavor this hypothesis, given that the estimated total power in the SLS (∼ 8.7 × 10 41 erg s −1 ) exceeds the one left in the HHS after dissipation. Alternatively, as proposed by some jet models (see e.g. Marscher 2014), the two states could originate in different cells within the radio core region. In the latter case, the estimated powers carried by these emitting regions would be only a fraction of the total jet power, since the jet still has to propagate and form the large scale radio structure. Indeed, a higher jet power (L jet = 1.8 × 10 44 erg s −1 ) was estimated by Balmaverde et al. (2008), based on a semiempirical relation considering the radio core luminosity (Heinz et al. 2007).
This simple one zone model provides a reasonably good description of the broadband SED in 3C 264, but it requires the emitting region to be very weakly magnetized. Indeed, as shown in Table 4, the ratio between the energy density of the magnetic field and the leptons, U B and U e respectively, is much smaller than unity both in the SLS (U B /U e ∼0.13) and the HHS (U B /U e ∼0.02). Similar ratios have been obtained for the misaligned jet in IC 310 (Ahnen et al. 2017), which is also detected at TeV energies, as well as in a subclass of TeV detected BL Lacs (Tavecchio et al. 2010) characterized by higher Doppler factors (δ 15). In the case of 3C 264, this result is in contrast with our initial hypothesis that acceleration of the bulk flow is taking place on the examined VLBI scales. In fact, according to theoretical models and simulations of magnetic jet launching (see e.g., Meier 2012, and references therein), in the jet acceleration zone most of the energy should be stored in the magnetic field, and gradually converted into kinetic energy of the bulk flow. In a highly magnetized jet, acceleration of particles radiating up to very high energies may be triggered by magnetic dissipation 15 Quillen et al. (2003) for the radio and microwave data, Leipski et al. (2009) for infrared data, and Chiaberge et al. (1999Chiaberge et al. ( , 2002 for the optical-UV. The X-ray points are from this paper, while the γ-ray data are from the 3FGL and 3FHL catalogs Ajello et al. 2017). The TeV spectrum was derived assuming a power law photon index ranging between 2.3 and 2.8 (Rieger & Levinson 2018) and a flux above 300 GeV as reported by Mukherjee (2018). Upper panel: the SLS (solid red line for the synchrotron curve and blue line for the SSC component) and HHS (dashed red and blue lines) are both produced by an emitting region in the core (Model 1), still within the acceleration region. A large particle dominance is required to reach the observed γ-ray fluxes. Lower panel: the emission of the SLS and HHS are produced in two different regions, in the core and in the layer, respectively (Model 2). This time, the core region is assumed to be Poynting flux-dominated (U B /U e =37). The SSC emission of the core component (solid blue line) is significantly below the observed γ-ray emission, and remains such when equipartition between particles and magnetic field is assumed (dot-dashed line). The HHS originates in a region in the layer at the end of the acceleration region, moving with a bulk motion Γ bulk ∼5. through reconnection, a process which requires the field and the emitting particles to be at least in equipartition (Giannios et al. 2009;Sironi et al. 2015). While a production of high-energy particles in the vicinity of the black hole magnetosphere may still be possible, as suggested by recent magnetospheric models drawn by pulsar-like acceleration mechanisms (see Aleksić et al. 2014;Hirotani & Pu 2016;Ahnen et al. 2017), in the following we test the acceleration scenario by imposing that U B /U e ≥ 1 in the core region (Model 2, Figure 10, bottom panel).
With this constrain, the synchrotron emission from the core can still reproduce the radio to X-ray observed SED. However, for B 0.1 G the SSC curve significantly under predicts the emission at higher energies, even in the most favorable case of equipartition between U e and U B . Therefore, we consider a case in which the low-and high-energy components originate in two separate emitting regions. The low-energy synchrotron emission is still produced in an accelerating, Poynting flux-dominated region in the core and accounts for the SLS in X-rays. The variable high-energy emission is instead associated to a blob at the end of the acceleration region (z∼11 pc), either in the layer (Γ b =5, Model 2) or in the spine (Γ b =8-10, Model 3). In Model 3, we assume a blob of the spine that is temporarily pointing towards us (θ ≤5 • ). This possibility finds support in the results of kinematics study, showing that some jet features move along curved trajectories, and could account for the fact that the source is only sporadically detected in the GeV to TeV band, as well as in Xrays in the hard state. Variations of the Doppler factor due to a change of θ have been proved able to consistently reproduce the long-term flux variability of other blazars (e.g. 4C 38.41, CTA 102 Raiteri et al. 2012Raiteri et al. , 2017, and references therein). In Figure  10, we show that this scenario can explain even the VERITAS detection in the TeV band and, with some variation of the layer (or spine) parameters the average GeV flux measured by Fermi.
The total power in the layer (or spine) component is in rough agreement with that of the core (a few 10 43 erg s −1 ), but the Poynting flux, which is energetically dominant in the latter, has been converted to kinetic power in the former. Here, for comparison with previous works, we have included in the jet energy budget a proton component (L p ) and assumed one cold proton for radiating electron. However, L p is not energetically dominant in any of the models, nor strictly necessary, hence a "light" electron-positron jet could be equally possible.
The total power estimated in Model 2 and Model 3 is consistent with the one independently derived by Balmaverde et al. (2008).

Summary
Motivated by the recent detection of 3C 264 in the TeV regime (Mukherjee 2018), we have conducted a multi-band study of its misaligned jet, combining VLBI radio images with high frequency data collected by different satellites (Chandra, XMM, Swift, Fermi-LAT) and adopting the inferred properties to model the broadband SED.
-By employing a multi-scale wavelet decomposition method (WISE code, Mertens & Lobanov 2015, we have investigated the bulk motion within the radio jet, which appears strongly limb-brightened. Superluminal apparent speeds up to ∼11.5 c are determined in our analysis. Such high speeds, in agreement with findings from optical kinematics on kilo-parsec scales (Meyer et al. 2015(Meyer et al. , 2017, limit the viewing angle to relatively small values (θ 10 • ). A trend of increasing speed with distance from the core is suggested on a de-projected scale of ∼11 pc (or ∼4.8 × 10 4 R S ). The plasma features are observed to move mainly along the jet limbs, but significant transverse motion is also detected when the jet filaments appear to intersect and cross each other. In agreement with theoretical models, which predict active collimation to take place along the acceleration region, the jet has a parabolic shape on the examined VLBI scales (r ∝ z 0.40±0.04 ).
-We have analyzed all the available X-ray data, including those obtained in 2018 through a dedicated Swift monitoring. The X-ray lightcurves indicate the presence of both long and short time variability. The luminosity increased almost by a factor of ∼3 from 2000 to 2018, while the spectrum became harder. The Swift data reveal day-scale variability by a factor of ∼2 between January and April 2018, and a possible flaring event two days after the end of the VERITAS campaign, in March 2018. -Based on the radio and X-ray results, we have modeled the broad-band emission assuming that the source switches between a low/soft and a high/hard state. A simple onezone synchrotron self-Compton model provides a reasonably good description of the SED in both states, but requires the core region to be strongly particle-dominated (U B /U e ∼0.02 − 0.13). Since we expect the plasma to be still magnetically dominated in the acceleration and collimation region, we have then tested a second scenario by imposing that U B /U e ≥ 1 in the core. The SED, as well as the expected total jet power, can be well reproduced with this assumption, provided that the high-energy component is associated to a second emission region. This could be located, for instance, at the end of the acceleration region, either in the spine or in the sheath.
A location of the high-energy emitting region in the jet acceleration zone has been considered in previous studies (e.g., Ghisellini 2012) but, for the assumed jet geometry and kinematics, the GeV-TeV emission was found to be significantly suppressed by γγ-pair production. The results of our kinematic study, however, indicate that the jet acceleration region in 3C 264 is relatively extended, which allows us to release the constraints on the size and location of the emitting region. There is, on the other hand, increasing evidence, based on VLBI observations, for the existence of such extended acceleration regions. For instance, about half of the AGN jets in the MOJAVE sample presents features moving on accelerating trajectories (Lister et al. 2016), on linear scales reaching up to hundreds of parsecs. In the case of 3C 264, its proximity and the relatively large black hole mass and viewing angle enable us to probe the innermost jet with high resolution in units of Schwarzschild radii, on scales where the bulk of the MHD processes is still taking place.