Jet collimation in NGC 315 and other nearby AGN

Aims. The collimation of relativistic jets in galaxies is a poorly understood process. Detailed radio studies of the jet collimation region have been performed so far in few individual objects, providing important constraints for jet formation models. However, the extent of the collimation zone as well as the nature of the external medium possibly confining the jet are still debated. Methods. In this article we present a multi-frequency and multi-scale analysis of the radio galaxy NGC 315, including the use of mm-VLBI data up to 86 GHz, aimed at revealing the evolution of the jet collimation profile. We then consider results from the literature to compare the jet expansion profile in a sample of 27 low-redshift sources, mainly comprising radio galaxies and BL Lacs, classified based on the accretion properties as low-excitation (LEG) and high-excitation (HEG) galaxies. Results.The jet collimation in NGC 315 is completed on sub-parsec scales. A transition from a parabolic to conical jet shape is detected at $z_{t}=0.58\pm0.28$ parsecs or ${\sim}5\times 10^3$ Schwarzschild radii ($R_{S}$) from the central engine, a distance which is much smaller than the Bondi radius, $r_{B}{\sim}92$ $\rm pc$, estimated based on X-ray data. The jet in this and in few other LEGs in our sample may be initially confined by a thick disk extending out to ${\sim}10^3$-$10^4$ $R_{S}$. A comparison between the mass-scaled jet expansion profiles of all sources indicates that jets in HEGs are surrounded by thicker disk-launched sheaths and collimate on larger scales with respect to jets in LEGs. These results suggest that disk winds play an important role in the jet collimation mechanism, particularly in high-luminosity sources. The impact of winds for the origin of the FRI/FRII dichotomy in radio galaxies is also discussed.


Introduction
Extragalactic jets are collimated outflows of relativistic plasma emanating from the center of active galaxies (see Blandford et al. 2019, for a recent review). On parsec scales, they are observed to propagate with opening angles as small as a fraction of a degree (e.g., Pushkarev et al. 2009), a striking feature which can be preserved up to kilo-parsec distances from the central engine. This high degree of collimation is thought to be achieved gradually, through physical mechanisms which are, at present, still unclear.
Magneto-hydrodynamic simulations of jet launching (McKinney & Gammie 2004;Hawley & Krolik 2006;McKinney 2006;Tchekhovskoy et al. 2011) predict the formation of a light, relativistic outflow powered by the rotational energy of the black hole, as described in the work of Blandford & Znajek (1977), as well as of a heavier and mildly relativistic wind powered by the accretion disk, as originally proposed by Blandford & Payne (1982). At their base, these outflows are expected to be slow and broad, with opening angles of tens of degrees. Along the collimation region, the bulk flow is accelerated by magnetic pressure gradients, thus the magnetic energy which dominates at the jet base is converted into kinetic energy. This description applies to magnetically-dominated cold jets, in which the acceleration and collimation mechanisms are shown to be necessarily co-spatial (Komissarov 2012, and references therein). In hot jets (internal energy-dominated) thermal acceleration can also occur in a conical flow which expands adiabatically in an external medium with a steep density gradient (e.g., Perucho & Martí 2007). Thermal acceleration could, at least initially, play a role (Vlahakis & Königl 2004), particularly for electron-positron jets and for winds powered by hot accretion flows. The magnetic field, helically wrapped around the jet axis, is thought to contribute not only to the acceleration but also to the collimation, through the "hoop stress" exerted by its toroidal component and/or through the pressure of the poloidal component (see Spruit et al. 1997, for a discussion). However, it has been shown that magnetic selfcollimation may not be sufficiently effective in relativistic flows to account for the observed collimation degree (e.g., Tomimatsu 1994;Beskin et al. 1998), and that some confinement from the ambient medium is required (Komissarov et al. 2007;Lyubarsky 2009). The nuclear environment in active galaxies is complex, and the nature of this confining medium is far from understood. At short distances from the black hole it may be the accretion disk itself or dense gas clouds in its surroundings, while at larger Article number, page 1 of 20 arXiv:2012.14831v1 [astro-ph.HE] 29 Dec 2020 A&A proofs: manuscript no. main distances there could be contributions from the shocked cocoon surrounding the jet or from the interstellar medium. In addition, disk winds are likely to play a crucial role, as they may effectively confine the inner relativistic jet (Bogovalov & Tsinganos 2005;Globus & Levinson 2016).
Several observational studies have been performed with the aim of determining the extension of the acceleration and collimation region, thus providing constraints for theoretical models. Such studies were mostly focused on the analysis of lowluminosity radio galaxies, such as NGC 6251 (Tseng et al. 2016), 3C 84 , 3C 270 (Nakahara et al. 2018), 3C 264 (Boccardi et al. 2019), NGC 1052 (Nakahara et al. 2020) and M 87 (Asada & Nakamura 2012;Mertens et al. 2016;Hada et al. 2016;Kim et al. 2018;Nokhrina et al. 2019), but some constraints for high-luminosity radio galaxies, like Cygnus A (Boccardi et al. 2016b,a;Nakahara et al. 2019), and blazars (Giroletti et al. 2008;Akiyama et al. 2018;Hada et al. 2018;Algaba et al. 2019;Traianou et al. 2020) were also provided. Moreover, part of the MOJAVE sample was investigated by Pushkarev et al. (2017); Kovalev et al. (2020). Asada & Nakamura (2012) have first shown that the inner jet in M 87 expands following a characteristic parabolic profile, and that this shape is preserved up to ∼10 5 Schwarzschild radii (R S ) from the black hole. At larger distances, the collimation stops and the flow assumes a conical shape. The transition between these two regimes occurs in the proximity of the stationary feature HST-1 and of the Bondi radius r B . Assuming that the jet is confined by a spherical, hot accretion flow of the Bondi type (Bondi 1952), it was proposed that the transition, as well as the recollimation at HST-1 (e.g., Levinson & Globus 2017), may be induced by a change in the ambient pressure profile beyond r B . A transition from parabolic to conical expansion at similar distances (∼10 4 −10 6 R S ) was later observed in other sources, most recently in ten nearby objects in the MO-JAVE 1 sample (Kovalev et al. 2020), including both low-power (Fanaroff-Riley I, FRI) and high-power (Fanaroff-Riley II, FRII) radio galaxies, as well as BL Lacs.
If the ambient medium does shape jets, however, one could expect to observe differences in the collimation profile of sources characterized, for instance, by different accretion modes. Most jets in FRII galaxies and flat spectrum radio quasars (FSRQs) are thought to be powered by radiatively efficient, geometrically thin accretion systems fed by cold gas (Baum et al. 1995;Sambruna et al. 1999;Grandi & Palumbo 2004;Grandi et al. 2006;Ballantyne & Fabian 2005;Hardcastle et al. 2007) and surrounded by a torus (Ogle et al. 2006). On the contrary, FRIs and BLLacs are fed by hot, radiatively inefficient accretion flows (Baum et al. 1995;Chiaberge et al. 2000;Balmaverde et al. 2006;Allen et al. 2006;Hardcastle et al. 2007;Yuan & Narayan 2014), and generally lack a torus (Chiaberge et al. 1999;Whysong & Antonucci 2004). The transition from radiatively efficient to inefficient accretion mode occurs around Eddington ratios of ∼0.01-0.1 (Narayan & Yi 1995), and is reflected by a change in the optical spectra, which, below this threshold, tend to lack strong high-excitation lines from the narrow-line region (Laing et al. 1994;Jackson & Rawlings 1997;Best & Heckman 2012;Heckman & Best 2014). The different energy output of low-excitation (LEG) and high-excitation (HEG) galaxies likely result in a different feedback on the environment (jet mode vs. radiative mode, see e.g., Heckman & Best 2014). While cold and hot disks are both expected to launch winds and a mildly relativistic jet sheath (e.g., Hawley & Krolik 2006;McKinney 2006;Ohsuga et al. 2009;Liska et al. 2019), the properties of these disk outflows 1 http://www.physics.purdue.edu/astro/MOJAVE/ and their impact on jet collimation in the two cases are not well constrained.
Ultimately, understanding the role of the environment in the collimation process requires a detailed investigation of the connection between the properties of the jet base and those of the AGN in the vicinity of the black hole in different luminosity classes. Imaging of the innermost jet base has now become possible in several nearby sources thanks to very-long-baseline interferometry (VLBI) at millimeter wavelengths (mm-VLBI, e.g., Boccardi et al. 2017, and references therein); however, studies of the nuclear environment, for instance of the hot X-ray emitting gas within the black hole sphere of influence, are resolutionlimited. The best chance is offered by nearby radio galaxies powered by very massive black holes, like M 87, where both the black hole sphere of influence and the event horizon have a large apparent size. Based on this criterion, our team has identified a first sample of suitable sources for performing high-resolution studies of jet formation in different AGN classes (Boccardi et al. in prep.). In this article we present a detailed analysis of jet collimation in one of the best targets identified so far, the lowluminosity source NGC 315. Being close (z = 0.0165), quite bright (0.15-0.4 Jy in the mm-band), and hosting a billion solar masses black hole (Sect. 3.1), this giant radio galaxy could be imaged through mm-VLBI with a resolution as high as ∼160 R S . In the second part of the paper, the results obtained for NGC 315 are compared with those obtained for other nearby AGN. The article is structured as follows: in Sect. 2 we describe the NGC 315 data set and its analysis; in Sect. 3 we present the jet collimation profile and investigate possible sources of confinement; in Sect. 4 we present a comparison between the properties of the collimation region in NGC 315 and in other nearby jets, based on results obtained for the MOJAVE sample and for other objects in the literature; we summarize our conclusions in Sect. 5. Throughout the article we assume a ΛCDM cosmology with H 0 = 71 km s −1 Mpc −1 , Ω M = 0.27, Ω Λ = 0.73 (Komatsu et al. 2009).

NGC 315: data set and analysis
The VLBI data set considered for the analysis of the parsec and sub-parsec scale jet structure in NGC 315 comprises eighteen observations at eight frequencies, spanning the range between 1.4 GHz and 86 GHz (Table 1). Several arrays were used in these observations: the Very Long Baseline Array (VLBA), the European VLBI Network (EVN), the High Sensitivity Array (HSA), and the Global Mm-VLBI Array (GMVA). Of these eighteen data sets, fourteen were calibrated in AIPS following the standard procedures, while the remaining four, specifically the data from February 2006 (at 8.4 GHz, 12.1 GHz, and 15 GHz) and from March 2005 (at 1.6 GHz), are calibrated data obtained from the MOJAVE archive and from the EVN archive respectively. To examine the large scale expansion profile of the source, two Very Large Array (VLA) calibrated data set at 1.4 GHz and 5 GHz provided by the NRAO VLA Archive 2 were also considered. All the data were imaged using DIFMAP (version 2.5e). The main information on the VLBI observations and clean maps are reported in Table 1. The basic properties of the VLA data sets are presented in Table 2. VLBI images of the source at 22 GHz, 43 GHz, and 86 GHz are published for the first time in this article, while previous analyses at 1.4 GHz, 5 GHz, and 8 GHz considering some of these data were presented by Cotton et al. (1999); Giovannini et al. (2001). All the images are shown in Figures 1-8 of the Appendix, while a multi-scale view of the source, including the Mega-parsec scale structure probed at 1.4 GHz by the VLA, is depicted in Figure 1. On parsec scales, the source presents a bright core and a straight one-sided jet. A faint counter-jet is often, but not always, detected, in agreement with previous results (Cotton et al. 1999;Giovannini et al. 2001

Alignment of maps at different frequencies
In order to combine the data obtained at each frequency and correctly reconstruct the jet expansion profile, it is necessary to refer all the measured distances to a common origin, ideally the central supermassive black hole. In each map, the origin coincides approximately with the position of the emission peak, which is however frequency-dependent due to synchrotron opacity effects at the jet base. Several methods have been developed for determining the opacity shift and aligning VLBI images (see e.g., O'Sullivan & Gabuzda 2009, and references therein). In this study we performed a 2D cross-correlation analysis taking into account optically thin regions of the jet at given pairs of frequencies (Table 3). We selected pairs of images from close-intime observations, in order to minimize the uncertainties arising from flux variability and proper motion of the plasma. Same-day observations were available in all but two cases. The closest in time available observations at 1.4 GHz and 5 GHz are separated by 11 days (November 1994), a time short enough for estimating the core-shift, since structural changes in radio galaxies are observed over relatively long timescales (∼months). The coreshift analysis at 43 GHz and 86 GHz was not performed, since the closest-in-time observations are separated by two months, and the core-shift is expected to be small.
Images in each pair were restored with a common circular beam, corresponding to the average equivalent beam b eq = b min × b maj of the two images, where b min and b maj are the minor and major axis of the natural beam, respectively (see Column 7 in Table 1). The pixel size was set to one tenth of the beam FWHM, which we estimate to be comparable to the resolution limits in the brightest regions of our images. Before performing the cross-correlation, each image was slightly shifted so that the pixel with peak flux density was exactly centered at the origin. Therefore the error on the shift determination in x and y corresponds to the in quadrature sum of the error on the core alignment, equal to one pixel, and the error on the images alignment, also equal to one pixel, since the 2D cross-correlation algorithm cannot determine shifts smaller than this size. The results obtained from the cross-correlation, including the correlation coefficients, are reported in Tables 3 and 4 for each frequency. In Figure 2 we show the dependence of the derived core positions z core relative to the mm-VLBI core (43 GHz and 86 GHz), as a function of frequency. By fitting a power law of the form z core ∝ ν −1/a , we obtain a = 0.84 ± 0.06. According to the Blandford & Koenigl jet model, this result is in agreement with the expectations for a synchrotron self-absorbed conical jet in equipartition, as observed in other extra-galactic jets (e.g., Hada et al. 2011;Sokolovsky et al. 2011). Note that the result of the fit is strongly influenced by the 1.4 GHz data point and, as it will be shown in the following, the core at this frequency is indeed located in the conical jet region.

Modeling
Aiming at measuring the transverse width of the jet as a function of distance from its base, we have modeled the data through the MODELFIT sub-routine in DIFMAP by fitting circular Gaussian components to the visibilities. The derived MODELFIT parameters at each frequency and epoch are reported in Tables 1-10 of the Appendix. After the modeling, two corrections were applied to the position of each feature, reported in Columns 3 and 4 of these Tables. The positions were first shifted so that the brightest core component would be centered exactly at zero; then a second shift was applied to correct for synchrotron opacity, as described in Sect. 2.1 and Tables 3-4. The shifted positions are shown in Columns 5 and 6 of Tables 1-10 in the Appendix. In addition to the integrated flux density, reported in Column 2, the component size d, assumed equal to the full width at half maximum (FWHM) of the Gaussian, is reported in Column 7 for each feature. Since the errors on the size determined by MODELFIT are extremely small, as the dominant systematic errors associated to the calibration and imaging procedure are not taken into account, in the following we will assume a more conservative error equal to one fourth of the FWHM.

Black hole mass
The black hole mass M BH in NGC 315 has been estimated in several works. Based on the correlation between the mass and the stellar velocity dispersion (M-σ relation), Woo & Urry (2002) found M BH ∼8 × 10 8 M , Satyapal et al. (2005) Beifiori et al. (2009). In the following, we will assume an intermediate mass value of M BH ∼1.3 × 10 9 M (Satyapal et al. 2005). Then, for the adopted cosmology, 1 mas corresponds to 0.331 pc and 2662 R S .

Expansion profile
In Figure 3, left panel, we examine the expansion profile of the approaching jet based on the analysis of the VLBI data set, which covers distances from the black hole ranging from the subparsec to the sub-kiloparsec. For the assumed black hole mass M BH ∼1.3 × 10 9 M , these distances translate to a range between ∼10 2 and ∼5 × 10 5 R S . For the de-projection we assume a viewing angle of 38 • , derived by Canvin et al. (2005) based on the large scale radio properties of the source. While a larger value of 49.8 • was obtained in a more recent study by Laing & Bridle (2014), an angle of 38 • is in better agreement with VLBI constraints provided by Giovannini et al. (2001), who derived an interval between 30 • − 40 • .
Not all the jet MODELFIT features (Tables 1-10 of the Appendix) were considered in our analysis. Those which were not fully resolved at a given frequency were filtered out. Specifically, we have excluded components with FWHM smaller than half of the beam minor axis. While the resolution limits may be smaller than this in some cases, this approach is justified by the large number of images and frequencies available, and ensures a reliable determination of the jet width on each scale. The reliability of these measurements is demonstrated by the fact that we obtain similar width values at different frequencies, thus at different resolutions. As explained in Sect. 2.1, the radial distances of each jet feature are relative to the positions of the 43 GHz and 86 GHz cores (origin of the axes in Fig. 3) while, ideally, the jet expansion profile should be described with respect to the black hole position. However, based on our core-shift study, the expected separation between the mm-VLBI core and the black hole is negligible with respect to the range of scales considered.
The presence of possibly two breaks in the jet expansion profile is suggested by the data in Figure 3, left panel. A flattening is observed in the inner jet, in the region described by the 86 GHz and 43 GHz data and, partially, by the 22 GHz data, while a possible steepening is hinted by the two outermost data points at 1 GHz. As we are firstly interested in the properties of the innermost jet regions, we have tested the existence of a transition in the jet shape by excluding the two outermost points and by fitting a broken power law to describe the jet width d as a function of de-projected distance z. This power law has the form: where k u and k d are respectively the upstream and downstream power-law coefficients, d t is the width of the jet at the transition, z t is the distance from the core at which the transition occurs, and h is a parameter regulating the sharpness of the break. The fit was performed based on a nonlinear least-squares (NLLS) Marquardt-Levenberg algorithm, implemented in GNUPLOT. The sharpness parameter was fixed, and several values were tested; the reduced χ 2 (χ 2 /dof) was found to be only weakly dependent on h, improving slightly for larger values. By fixing h to 100 and letting all the other four parameters free, we achieved a reduced χ 2 (χ 2 /dof) of 1.35. A single power-law (d ∝ z k ) fitted for comparison in the same region does not well describe the innermost data points, and the reduced χ 2 is larger (1.44). The broken power-law fit yields a power-law index k u = 0.45 ± 0.15 in the inner jet, which thus has a close to parabolic shape, while a close to conical shape is determined outwards, with k d = 0.84 ± 0.04. The transition occurs at a distance z t = 0.58 ± 0.28 de-projected parsecs, corresponding to ∼5 × 10 3 R S . The two outermost VLBI data points at 1.4 GHz, which were so far excluded by our analysis, hint at a possible further steepening of the expansion profile. In order to examine how the profile evolves on larger scales, in the right panel of Fig. 3 we add the VLA data, which enable us to probe the jet up to a distance of almost one Megaparsec (10 9 − 10 10 R S ) from the central engine. The unresolved VLA MODELFIT components were filtered out following the same criteria described for the VLBI data. In the case of the 1 GHz data, which describe the entire large scale structure of the approaching jet, we have also excluded data points beyond the ninthy-degrees bend of the jet direction towards south ( Figure 1). The analysis of the multi-scale data set confirms our result that the jet collimation is completed on small scales, and also shows that, except for some local oscillations of the jet width (e.g., between 10 4 − 10 5 pc) the global jet shape remains conical for several orders of magnitude in distance after the initial collimation. A single power-law fit performed considering the entire jet after the transition (i.e., excluding the innermost 0.58 pc) yields a power-law index k = 0.99 ± 0.01, only slightly steeper than the one derived based on the VLBI data set. The power-law index does not change significantly if the innermost parabolic region is included in the fit (k = 0.98 ± 0.01), A transition from a parabolical to conical jet shape is observed on sub-parsec scales. The fitted red line is a broken power-law. The dotted lines delimit the uncertainty interval of the transition distance. The latter is much smaller than the estimated Bondi radius. Right panel: VLBI and VLA data. After the transition, the global jet shape remains conical beyond the Bondi radius and up to Megaparsec distance from the black hole.
but the reduced χ 2 becomes larger (2.30 vs 2.07). In Figure 4, the residuals of a single power-law fit performed considering the entire data set as well as the VLBI data only are displayed. In both cases, a single power-law cannot well describe the expansion profile on sub-parsec scales, where the residuals are of positive sign. A broken power-law fit considering the entire VLA plus VLBI data set was also attempted. Even in this case, the results are consistent with a small scale shape transition. The fit becomes more sensitive to the assumed sharpness parameter (which we varied between 1 and 50, as the fit did not converge for larger values) even though the reduced χ 2 does not change significantly (χ 2 /dof∼1.76 − 1.79). The assumption of a soft break (h=1) yields a transition distance z t = 1.29 ± 1.23 pc with k u = 0.48 ± 0.22 and k d = 1.01 ± 0.01, while the assumption of a sharp break (h=50) places the transition at z t = 6.32±2.67 pc with k u = 0.74 ± 0.05 and k d = 1.01 ± 0.01, thus with the inner jet shape deviating significantly from a parabola. By letting all the parameters, including the sharpness h, free to vary, we obtain h = 1.50 ± 1.89 and a transition distance z t = 2.75 ± 2.88 pc with k u = 0.62 ± 0.19 and k d = 1.01 ± 0.02, in agreement with the VLBI results within the uncertainty interval. Overall, our findings match those obtained by Pushkarev et al. (2017) and Kovalev et al. (2020), who performed a detailed analysis of the 15 GHz VLBI stacked image, probing scales comprised between ∼1 − 10 de-projected parsecs. The power-law indices obtained by these authors by fitting a single power-law are k = 0.86 ± 0.01 and k = 1.07±0.05, respectively. Kovalev et al. (2020) also tested the presence of a transition along the 15 GHz jet, with negative results, and indeed the transition we detect is suggested to occur on smaller scales, unresolved at 15 GHz. In conclusion, the present data are compatible with the occurrence of a sub-parsec scale transition from parabolic to conical expansion, after which the global jet shape remains conical.

Jet confinement from the environment
As discussed in Sect. 1, the transition from parabolic to conical shape is observed in M87 and in several other sources on relatively large scales, at distances of the same order as the sphere of influence of the black hole. Like M 87, NGC 315 is a cluster member and, as shown by Worrall et al. (2007), is powered by a radiatively inefficient nucleus. Under the assumption that the hot accretion flow can be well described by the Bondi theory, we can test the hypothesis that the jet shape transition occurs in the proximity of the Bondi radius. The Bondi radius r B is the radius inside which the gravitational influence of the central black hole dominates over the thermal motion of the gas in the accretion flow. Its expression is (see e.g., Frank et al. 2002): where G is the gravitational constant and c s is the sound speed.
Since the sound speed depends on the gas temperature T , assuming an adiabatic index γ = 5/3 the Bondi radius can be also expressed, in convenient units, as (see Russell et al. 2015): The expected pressure profile of the interstellar medium inside the Bondi sphere is p ∝ z −2 and, according to analytical and numerical models (Tchekhovskoy et al. 2008;Komissarov et al. 2009;Beskin et al. 2017), a jet propagating in a medium with such a pressure gradient develops the characteristic parabolic shape. Based on Eq. 3, estimating the Bondi radius requires the knowledge of the black hole mass and of the temperature of the accretion flow in the nuclear region. A detailed analysis of the Xray emission in NGC 315 was carried out by Worrall et al. (2003Worrall et al. ( , 2007 using sensitive Chandra data. The jet, which presents several X-ray-bright knots of synchrotron emission, is embedded in a hot gaseous atmosphere, as inferred from the presence of an Xray thermal emission component in the spectrum. Such a feature is also detected within a circle of 1 arcsecond radius centered around the nucleus, in addition to a dominant, mildly absorbed power-law component possibly associated to the jet. The nuclear hot atmosphere has a temperature kT = 0.44 +0.08 −0.04 keV. Assuming a black hole mass of 1.3 × 10 9 M , we then estimate r B = 92 pc (dark blue vertical line in Fig. 3). Even if the black hole mass was smaller by a factor of several, this radius would still be much larger (two orders of magnitude) than the distance at which we observe the jet shape transition from parabolic to conical, z t = 0.58 ± 0.28 pc. This result would also not be affected by varying the gas temperature within the given small uncertainty. By assuming slightly different parameters, a similar value for the Bondi radius in NGC 315 was recently obtained by Inayoshi et al. (2020), who estimated an uncertainty of 50 percent 3 . Then, if the transition from parabolic to conical shape in NGC 315 is induced by a change in the external pressure gradient, this must occur not in the proximity of the black hole sphere of influence, but in the vicinity of the black hole, on sub-parsec scales. Based on the analysis of the VLA data in Fig. 3, right panel, we also note that no discontinuity is observed in the expansion profile after the jet crosses the Bondi radius: the jet shape is close to conical both inside and outside the Bondi sphere.
It is important to point out that the Bondi theory is likely over-simplified, not taking into account important aspects like viscosity, the presence of magnetic fields or a non-zero angular momentum of the accretion flow. In fact, simulations assuming more realistic physical conditions (e.g., Gaspari et al. 2013) have suggested that, even in the case of radiatively inefficient AGN, the accretion is likely to be cold and chaotic, rather than hot and spherically symmetric. The presence of cold gas in NGC 315 on scales of tens of parsecs was revealed by VLBI HI absorption studies (Morganti et al. 2009). However, there is no evidence for the existence of a dusty torus feeding the nucleus (X-ray data reveal only mild intrinsic absorption), nor of a cold thin disk. The nuclear emission in NGC 315 may be adequately described by more complex models of hot accretion flows, such as advection-dominated accretion flows (ADAF). Wu et al. (2007) showed that the X-ray spectrum in NGC 315 can be well fitted by an ADAF. Interestingly, such thick disks are suggested to have an outer radius of the order of ∼10 3 R S (e.g., Mahadevan 1997), comparable to the distance at which the jet break is observed (∼5 × 10 3 R S ). Thus it is possible that the jet in NGC 315 is initially confined by the thick disk itself, and starts to freely expand beyond the outer disk radius. Alternatively, analytical and semianalytical models (Lyubarsky 2009;Beskin et al. 2017;Kovalev et al. 2020) have shown that the transition to a conical jet shape can occur also in the presence of a single external pressure profile, and that the location of the transition may be strongly dependent on the initial jet magnetization. The external medium may not necessarily be the gas in the accretion flow. Globus & Levinson (2016) have recently extended the work of Lyubarsky (2009) to describe the case in which the Poynting-flux dominated jet is confined by a hydrodynamic wind layer. In this semianalytic model, the wind layer was also shown to be characterized by a pressure profile p ∝ z −2 , thus conferring the characteristic parabolic shape to the central jet filaments. The collimation was found to be effective for sufficiently high ratios (> 0.1) of wind power to jet power, and the extension of the collimation region was shown to be dependent on the thickness and power distribution of the wind layer. The possibility that this and other jets are collimated by accretion disk winds is further discussed in the following.

The parabolic jet in different sources
Substantial observational evidence has accumulated so far concerning the existence of parabolically expanding jet bases in AGN. The study of NGC 315 presented in this paper shows another example of this, but also indicates that the jet collimation can be completed already on sub-parsec scales. In this section we investigate how the properties of the jet collimation region vary in different sources, combining our findings with results from the literature. In particular, one open question which we would like to address concerns how the presence of disk-winds affects jet collimation, and how are the properties of these winds related to the nature of the accretion flow. Millimiter-VLBI observations of M 87 indicate that the jet base is anchored in the vicinity of the ergosphere, having a transverse size of the order of few R S (e.g., Doeleman et al. 2012;Kim et al. 2018). On the other hand, imaging of the two-sided jet in Cygnus A has revealed that its base is quite broad, with a transverse size of ∼230 R S (Boccardi et al. 2016a), thus part of the jet may be launched from the outer regions of the accretion disk. M 87 and Cygnus A clearly differ both in the jet power, of ∼10 43 erg s −1 in the former and ∼10 46 erg s −1 in the latter, and in the type of the accretion disk they host, the former being powered by an ADAF (e.g., Di Matteo et al. 2003) and the latter by a heavily absorbed cold disk (e.g., Young et al. 2002;Reynolds et al. 2015). The two sources also differ in the jet large scale morphology, FR I and FR II respectively. In the following we compare the mass-scaled properties of the jet collimation zone in NGC 315 and in several other sources of both high-and low-power. We take into account results from studies of individual nearby objects, namely Cygnus A (Boccardi et al. 2016b), 3C 264 (Boccardi et al. 2019, and NGC 4261 (Nakahara et al. 2018), as well as results obtained for the MOJAVE sub-sample analyzed by Pushkarev et al. (2017); Kovalev et al. (2020), limiting our analysis to sources with redshift smaller than 0.15 and with known Table 5. Sample of sources considered for the analysis in Section 4, sorted by redshift. Column 1: Source name (B1950). Column 2: Other name. Column 3: Redshift. Column 4: Source classification. RG: Radio galaxy, BL: BL Lac, NLSy1: Narrow-line Seyfert 1. Column 5: Log of the black hole mass M BH , expressed in units of solar masses. Column 6: Ratio between the X-ray luminosity in the 2-10 KeV band (reference in the Table note) and the Eddington luminosity L Edd = 1.3 · 10 38 M BH /M . Column 7: Classification as HEG: High-excitation galaxy or LEG: Low-excitation galaxy. A dividing limit L X /L Edd = 1.1e−3 was assumed. Column 8: Conversion factor, parsec per milliarcsecond. Column 9: Conversion factor, Schwarzschild radii per milliarcsecond. Column 10: Viewing angle. Column 11: Power-law index defining the inner jet shape. Column 12: Transition distance in units of projected milli-arcseconds. References as in Column 11.

B1950
Other  (2002) black hole mass. We exclude BL Lacertae, whose peculiar properties will be discussed in a dedicated paper (C. Casadio et al. in preparation). While observational constraints exist also for highredshift sources based on the analysis of the MOJAVE sample, in this case we prefer to focus on nearby objects, for which several estimates of the black hole mass are given in literature and a higher spatial resolution can be achieved in VLBI observations. A summary of the properties of the twenty-seven sources in our sample, thirteen radio galaxies, thirteen BL Lac objects, and one Narrow-line Seyfert 1, is presented in Table 5. In Column 7 we report the classification as low-excitation (LEG) or highexcitation (HEG) galaxy. As mentioned in Sect. 1, this division is indicative of the AGN accretion regime and is based on the optical spectroscopic properties of the nucleus. Since this optical classification is not available in the literature for all the sources in our sample, in order to adopt a uniform criterion we consider instead the X-ray luminosity L X as a proxy of the accretion activity. Following the work of Macconi et al. (2020), who have inves-tigated the X-ray properties of a large sample of high-excitation (HERG) and low-excitation (LERG) radio galaxies, we assume a ratio L X /L Edd = 1.1e−3 4 as the limit above (below) which the source is classified as HEG (LEG). This method is applied to all the radio galaxies and to the NLSy1; the BL Lacs are instead assumed to be LEG since they form the beamed parent population of LERG. In the following we use the abbreviations HEG and LEG when we generically refer to high-and low-excitation galaxies respectively, and HERG and LERG when we limit our analysis to radio galaxies.
The jet viewing angles, reported in Column 10, are assumed based on results from the literature. In the case of several BL Lacs for which more stringent constraints are unavailable, the adopted angles were used in the cited works to reproduce the broadband spectral energy distribution. The jet shape  is defined by the power-law index k, reported in Column 11. For jets showing a parabolic shape in the innermost regions, k is the power-law index in the parabolic region. About half of the sources in the sample, mostly radio galaxies, show parabolic expansion (0.3 ≤ k ≤ 0.7), while the shape of the other objects, mainly BL Lacs, is either conical or hyperbolic (k ≥ 0.7). For most of these parabolic jets, a possible location for the transition distance has been suggested in the literature through the observation of a jet break in the expansion profile and/or based on the jet kinematic properties. The transition distance expressed in projected milli-arcseconds is reported in Column 12. Note that some sources in the MOJAVE sub-sample (e.g., 3C 84, NGC 1052, NGC 6251) have also been studied by other authors (e.g., Tseng et al. 2016;Giovannini et al. 2018;Nakahara et al. 2020) but we refer to the MOJAVE results for uniformity. The dedicated studies are discussed in the text.

Literature data
For the visualization of the 15 GHz data from the MOJAVE survey we follow a different but complementary approach with respect to the one used by Pushkarev et al. (2017); Kovalev et al. (2020). These authors have analyzed the jet expansion by examining transverse intensity profiles in stacked images, and have determined the jet shape. For the purposes of our comparison, we consider instead the results from Gaussian fitting in the visibility domain, presented by Lister et al. (2019). The availability of many observing epochs for each source implies that several measurements of the jet width at similar distance from the core are provided. Since one of our main goals is to compare the full transverse width of different jets, and the jet cross-section may be not fully visible in a single epoch, we adopt the following empirical method, which was tested against the stacking 5 : for every five measurements of the jet width we select the maximum; then, to reduce the scatter, we compute a 5-point average of the maxima. In the single case of 1ES 1959+650, where the number of points is insufficient, we select the maximum every 3 measurements and we compute a 3-point average. When the jet is two-sided, we consider width measurements in both sides. To ensure that we are taking into account resolved jet regions, we exclude data points in the inner 0.5 mas. For sources showing a parabolic shape on the MOJAVE scales, we also exclude data points beyond the transition distance, whenever this was determined (Column 12 in Table 5). In the case of M 87, and possibly of 3C 84 and Mrk 501, the MOJAVE data describe only part of the parabolic region. The distances of the jet features reported by Lister et al. (2019) are computed with respect to the VLBI core position, while we ideally want to compare the width of different jets at the same distance from the black hole. Even though this offset is usually negligible on the considered scales, we apply the correction when the offset is known (based on the analysis from Kovalev et al. 2020, Column 7 in Table 3). No shift was applied in the case of Mrk 501 and I Zw 187, for which a negative offset was determined. The radio galaxy 3C 264 is also part of the MOJAVE sample, but was not analyzed by Lister et al. (2019). For this source we include data from Boccardi et al. (2019) , Fig. 9, where the jet expansion has been investigated in the image domain by considering the stacked 15 GHz image. We take into account data for the innermost jet region until the recollimation at ∼11 pc and we apply the same method described above, except that, since in the stacked image we already recover most of the jet cross-section, we do not average maxima but the data points directly.
For Cygnus A, we report 43 GHz data from Boccardi et al. (2016b), while for NGC 4261 the data points are not publicly available and we refer to Fig. 8 of Nakahara et al. (2018).

Mass-scaled expansion profiles
The mass-scaled jet expansion profiles are shown in Fig. 5, left and right panels. In the left panel we compare the profiles of fourteen sources showing a parabolic expansion. The data points are coded with a different color and symbol for each source. A power law of the form r ∝ z k is also fitted to the data of M 87, Cygnus A, and 3C 120 (magenta, green, and violet line, respectively in Fig. 5, left panel), with the coefficient k fixed to the value reported in Table 5 for each source. Indeed, our aim is not to determine the jet shape, an information that already exists in the literature, but to examine where these profiles lie with respect to each other. The comparison shows that the sources are not all aligned along the same profile. NGC315 and M 87 present the "thinnest" jet, and sources like NGC 6251, 3C 264, 3C 371, Mrk 501, and Mrk 180 lie on a similar or slightly upshifted profile. On the contrary, the rest of the radio galaxies, as well as the narrow-line Seyfert 1 1H 0323+342, present "thicker" jets, with the maximum width being observed in 3C 120 and in 1H 0323+342. None of the radio galaxies classified as HERG lies on the M 87 profile, but two of the LERG (3C 84, NGC 1052) are aligned with the profiles of Cygnus A or 3C 111. While, in Fig. 5, we do not report data for the LERG NGC 4261, Nakahara et al. (2018) showed that this object presents a jet width which is intermediate between the one of Cygnus A and M 87.
The relation between the source classification and the collimation profile is further explored in the right panel, where all the sources in Table 5, not only those with parabolic shape, are compared. As evident from the figure, the BL Lacs are on average probed on larger scales than radio galaxies, which is the most likely reason for the observation of mainly conical shapes (i.e., the jet collimation occurred on unresolved scales). By expanding faster, the BL Lac jets reach, at the largest distances (∼10 7 R S ), widths which are comparable with those of the HERG (green points) and of the NLSy1 1H 0323+342 (dark blue points). However, the BL Lacs conical profiles smoothly connect to the parabolic profiles of M 87 and other LERG (magenta points) observed on smaller scales. If we assume that the jet expansion in the collimation zone is described by a single parabolic profile from the launching up to the transition region, as observed for M 87 (e.g., Nokhrina et al. 2019), than these results suggest that the HERG and the NLSy1 in our sample have jets launched at larger disk radii than jets in BL Lacs and in most of the LERG.
Concerning the two "outliers", NGC 1052 and 3C 84, we note that these are both peculiar LERG. NGC 1052 presents features which are unusual for a low-luminosity AGN, since the nucleus is obscured by a high column density torus (Kadler et al. 2004), and a broad iron line of unclear origin has been detected based on X-ray observations (Brenneman et al. 2009). Reb et al. (2018) proposed that that NGC 1052 hosts a hybrid accretion disk in a transition regime, i.e. formed by a central ADAF and an outer thin disk (e.g., Esin et al. 1997). A disk rather than a spherical accretion flow has been proposed also for 3C 84 (Plambeck et al. 2014) based on a study of Faraday rotation in the nucleus, and space-VLBI observations with the RadioAstron telescope have shown that the jet base is indeed wide, and the jet may be anchored to the outer accretion disk . As discussed further in Sect. 4.5, both these sources show signs of young and/or restarted activity.

Caveats
Before discussing further the possible implications of the results in Fig. 5, we comment on some caveats to be taken into account, and that are also relevant for the comparison in the next Sect. 4.3. Firstly, there exists the possibility that some of the profiles in Fig. 5 are misplaced. The main source of uncertainty is in this case the black hole mass. The assumed values may be incorrect, and a different mass will cause a given profile to shift along the x and y axis by the same amount. Since the Schwarzschild radius is directly proportional to the mass, the shift will also be directly proportional to the change in mass, and the profiles will shift along parallel lines. Another source of uncertainty is related to the jet viewing angle. While the mass-scaled transverse widths (y axis) only depend on the mass, the de-projected distances on the x axis depend, in addition, on the assumed θ, being inversely proportional to sin θ. Our choice to limit the sample to nearby and well studied objects is aimed at mitigating the impact of the uncertainty on the mass and jet orientation, but incorrect assumptions are possible. In the case of 1H 0323+342, in particular, León Tavares et al. (2014) have proposed a mass larger by one order of magnitude; the radio properties of some of the BL Lacs in the sample are also compatible with a more misaligned jet orientation (Giroletti et al. 2004b(Giroletti et al. , 2006 Concerning the conclusions we can draw on the jet origin, based on a back-extrapolation of the jet profiles, it should be noted that the assumption that the jet expands following a single parabolic profile in the collimation region may not be valid. While this is observed in the best studied case of M 87, we cannot exclude that jets in other sources present more complex profiles at their base. In the case of NGC 1052, one of the two LERG placed on the HERG region, high resolution VLBI observations with the GMVA (Baczko et al. 2016) do not clearly resolve the jet base, thus it is plausible that this jet was initially narrow and has experienced a fast expansion at larger distances. Moreover, Nakahara et al. (2020) did not confirm a parabolic expansion for this source, suggesting instead a transition from cylindrical to conical shape at a distance of ∼10 4 R S . 3C 84, i.e., the other LERG in the HERG region, shows evidence for a more complex profile as well. In this case the jet base was well resolved by space-VLBI observations , and the inner-jet shape was found to also approach a cylinder.
Finally, another element to consider is the impact of relativistic effects. If a wide and mildly relativistic, disk-launched component was present in some of the BL Lacs, its detection could be prevented by the strong relativistic Doppler boosting, which amplifies the emission from the narrow and fast spine. However the good alignment between the BL Lacs profiles and those of their misaligned parent population, i.e. the low-power radio galaxies, suggests that BL Lacs intrinsically miss this component. Similarly, the jet orientation could influence the observed difference between HERGs and LERG, as a mildly relativistic jet sheath could appear more or less prominent in radio galaxies seen at different angles. This also seems to be not a concern, given that different widths are observed in radio galaxies seen at similar angles. For instance, the jets in M 87 and 3C 120 have very similar orientation (θ ∼ 16 • − 19 • ), but the mass-scaled jet widths differ by one order of magnitude. In fact, unlike M 87, 3C 120 hosts a powerful classic disk with high Eddington ratio Vasudevan & Fabian 2009;Cowperthwaite & Reynolds 2012), and a disk origin for this jet was suggested based on a combined X-ray and radio monitoring by Marscher et al. (2002). To test further the impact of relativistic effects, it would be interesting to examine where FSRQs, not included in our sample due to the redshift cut, would be placed in Fig. 5. The only high-excitation source with a blazar-like jet orientation is 1H 0323+342, whose jet profile appears to be well aligned with that of 3C 120. This would indicate that a jet orientation towards the line of sight does not prevent the observer from detecting an extended jet sheath, if this is present.
Ultimately, a more solid investigation of all these aspects will be possible in the future by considering larger samples at the highest possible resolution, in order to reduce the impact of the uncertainties related to the mass and viewing angle and to the back-extrapolation of the profiles, and to test how the appearance of the jet structure depends on relativistic effects.

Transition distance
Limiting our analysis to sources which show a transition in the jet shape and/or in the kinematic properties of the jet (Col. 11 in Table 5), we examine how the transition distances are distributed based on the AGN classification as HEG or LEG. With the exception of Mrk 180 and 3C 371, this excludes most of the BL Lacs, which are either already freely expanding or show a purely parabolic profile on the examined scales (this is the case for Mrk 501 only). For the latter reason, 3C 84 is also excluded.
Results are presented in the histograms in Fig. 6. In the left panel, the transition distance is expressed in de-projected parsecs. When considering the low-excitation sources, it is evident that NGC 315 is not the only object showing a subparsec scale transition, which is observed also in other two lowluminosity jets, NGC 1052 and NGC 4261. At least in the case of NGC 4261, this distance is also much smaller than the Bondi radius, estimated by Balmaverde et al. (2008) to be ∼32 pc, but it may match the actual extent of the hot gas phase as inferred from a spectral fitting performed by Kaburaki et al. (2010). None of the jets in HEGs presents a sub-parsec scale transition, and the transition distances seem to be shifted to larger radii than in LEGs. However, a two-sample Kolmogorov-Smirnov (KS) test does not allow us to conclude that the two classes are drawn from a different population, as the p-value 6 is p ks = 0.32. In the right panel of Fig. 6, the transition distance is expressed in deprojected Schwarzschild radii. With the exception of Cygnus A, all the jets in HEGs show a transition on scales > 10 6 R S , while all the LEGs are below this limit. In this case, the separation between the two classes becomes more evident, and the KStest indicates that they are drawn from a different population (p ks = 0.02).
Additional results from the literature for sources not included in our sample are consistent with this picture. In the LERG Centaurus A (Müller et al. 2011), the jet was observed to be freely expanding already on scales of ∼0.1 pc, or ∼10 4 R S for M BH = 5 × 10 7 M (Neumayer 2010), in agreement with our 6 Here the p-value is the largest probability that, given two random samples drawn from identical populations (null hypothesis), the difference between the two cumulative distributions is at least as extreme as the observed one. The null hypothesis is assumed to be rejected for p < 0.05. See, however, the discussion on p-values by Wasserstein & Lazar (2016). findings for NGC 315 and other LERG. On the other hand, studies of the jet shape in some high redshift FSRQs (Algaba et al. 2019;Traianou et al. 2020) locate the shape transition at distances > 10 6 R S , as found for the HEGs in our sample.
The different extent of the collimation zone in low-and high-luminosity sources is well matched by the findings of Potter & Cotter (2015), who have investigated the extent of the acceleration region in a sample of blazars. FSRQs were shown to reach their terminal Lorentz factor on scales larger than 10 5 R S , while shorter distances were inferred for BL Lacs. This also supports the idea that, both in high-excitation and in low-excitation sources, the jet acceleration and collimation processes are co-spatial, as expected based on theoretical predictions for magnetically-driven cold outflows (Komissarov et al. 2007;Tchekhovskoy et al. 2008;Lyubarsky 2009).

Implications for the disk-jet connection
Keeping in mind the caveats discussed in Sect. 4.2, a relation between the properties of the jet collimation region and the properties of the accretion disk is suggested by Fig. 5 and 6. According to theoretical models and simulations, both thin disks (e.g., Blandford & Payne 1982;Fendt 2006;Liska et al. 2019) and geometrically-thick hot disks (e.g., Blandford & Begelman 1999;McKinney 2006;Begelman 2012;Mościbrodzka & Falcke 2013;Mościbrodzka et al. 2016) can launch collimated outflows. Due to the higher mass loading and lower speed, the disk-driven jet is expected to dominate the emission in radio galaxies with respect to the de-boosted black hole-launched jet. This is confirmed in observations by the direct imaging of limb-brightened jet structures (e.g., Boccardi et al. 2016b;Mertens et al. 2016;Giovannini et al. 2018) as well by kinematic studies of radio galaxies, which generally show much lower intrinsic speeds than measured in blazars (Lister et al. 2019). As these properties are observed in high-luminosity and low-luminosity radio galaxies alike, a jet sheath must be produced from disks spanning different accretion regimes. Our results, however, indicate that the disk-driven jet in LERG originates at small disk radii (few R S , as measured in M 87), and indeed the expansion profiles of most of the LERG are well aligned with those of BL Lacs, which are expected to be dominated by the black hole-launched spine (see e.g., Ghisellini et al. 2014). This result is in broad agreement with models of jet launching from ADAFs (e.g., Cao 2011; Yuan & Narayan 2014, and references therein), which predict the formation of a thin and mildly-relativistic outer layer. ADAF models also predict the launch of a non-relativistic disk-wind component carrying the bulk of the disk mass outflow and spanning a large solid angle. There is no evidence for such a component based on the analyzed VLBI images, at least in the considered frequency regime. The jet profiles in HEGs, on the other hand, are all shifted upwards, and a back-extrapolation down to the jet base suggests that the jet sheath is launched at larger disk radii. Taking as a reference Cygnus A, which shows the thinnest jet among HEGs and for which an initial jet width of ∼200 R S was measured based on GMVA observations (Boccardi et al. 2016a), the present data suggest that thin disks could launch collimated winds with an initial outer radius 100 R S . This possible difference in the outer radius of the jet sheath is accompanied by a different extent of the collimation region in HEGs and LEGs (Fig. 6). Modeling of jet collimation by disk winds, presented by Globus & Levinson (2016), revealed a direct link between the wind outer radius and the collimation radius: for a given wind power, larger wind radii correspond to more extended collimation zones. A sufficiently high ratio (> 0.1) of wind power to jet power is required for this process to be efficient. When this condition is verified in reality, is a matter of debate. In recent simulations presented by Hervet et al. (2017), the diverse kinematic behavior of VLBI knots in blazars of different powers could be well explained by varying this ratio. Except for the least powerful class among BL Lacs (that of the High-frequency peaked BL Lacs, HBLs), whose properties could be reproduced by assuming an absent or very weak wind, ratios larger than 0.3 were suggested for blazars. A question remains concerning the portion of these winds which is actually detected in VLBI observations. When attempting to model the M 87 jet collimation profile, Globus & Levinson (2016) have suggested that the radio emission is produced in the shocked interface between the relativistic jet and the outer wind, which is undetected. Observational constraints on extended disk winds may be provided through other methods. For instance, we note that for all the HERGs in our sample (except PKS 1514+00) the detection of ultra-fast outflows was reported based on X-ray observations (Tombesi et al. 2010(Tombesi et al. , 2014Reynolds et al. 2015). These outflows, whose launching mechanism is unclear, are suggested to be characterized by mildly relativistic speeds, to originate at disk radii of 10 2 − 10 4 R S (in agreement with our findings), and to carry a significant fraction of the jet kinetic power. Thus collimation via the action of disk winds, where by disk winds we mean a mildly relativistic jet sheath plus possible broader outflows, appears to be a viable mechanism, especially for high-luminosity sources.

Implications for the FRI-FRII dichotomy
In this final section we discuss the possible impact that a more or less prominent disk-launched jet may have on the long-term jet stability and evolution. A surrounding jet sheath 7 , characterized by a smaller sound speed with respect to the central hot spine, provides further stability to relativistic outflows due to the increase in the jet inertia and the drop in instability growth rates (Hardee & Rosen 2002;Hardee & Hughes 2003;Perucho et al. 2005;Hardee 2007;Perucho & Lobanov 2007;Mizuno et al. 2007;Vega-García et al. 2019;Perucho et al. 2019). Furthermore, its presence mitigates the impact of surface perturbations that can be induced by the penetration of stars (which has been recently appointed as a possible triggering mechanism for jet deceleration of FRI jets, see Perucho 2020), since its smaller sound speed implies a slower propagation of the turbulent layer towards the jet axis. This provides the jet channel with time to reach the intergalactic medium without being decelerated (see the discussions in, e.g., Perucho 2012Perucho , 2016. In summary, a surrounding sheath or wind component shields the inner spine against entrainment from the ambient interstellar medium. According to this, we could expect a relation between the jet widths plotted in Fig. 5 and the large-scale morphologies of those jets. At zero order, this is the case, with the lower line populated by FRI jets alone, in addition to BL Lacs. Among the upper lines corresponding to thicker jets we observe the FRII radio galaxies Cygnus A, 3C 111, PKS 1514+00, and 3C390.3, plus other sources that do not develop clear FRI morphologies. NGC 1052 is a young source, and the large scale morphology may be consistent with that of a young FRII (see Wrobel 1984); 3C 84 is a recently reactivated source showing a bright hotspotlike feature (e.g., Nagai et al. 2017) which would approach it to a (temporary, at least) FRII morphology; 3C 120, while is more often classified as an FRI, presents a peculiar morphology, with a strongly bent jet and edge-brightened lobe (Walker et al. 1987).
The origin of the FRI/FRII dichotomy was widely discussed in the literature (e.g. Wold et al. 2007;Gendre et al. 2013;Tchekhovskoy & Bromberg 2016;Mingo et al. 2019), being attributed to a diversity in the nuclear properties or in the environmental conditions of the host. The different characteristics of disk winds ejected from HEGs and LEGs, suggested by our work, may provide a crucial link between these two scenarios, since the nature of the accretion disk has a direct impact on the properties of the environment, i.e. the wind, by which the jet is collimated and stabilized. We note that, while FRI and FRII morphologies usually develop in LERG and HERGs respectively, cross-populations formed by FRI-HERGs and FRII-LERG also exist, but differ significantly in size. FRI jets powered by powerful nuclei are highly infrequent (3C 120 is one of the few examples), while FRII morphologies are often found in low-luminosity galaxies (Heckman & Best 2014;Macconi et al. 2020, and references therein). This may suggest that when the jet is confined by a strong wind, the FRI morphology rarely develops. On the contrary, sources hosting inefficiently accreting disks that can still launch a relatively powerful wind may man-age to develop an FRII morphology as long as the jet reaches the intergalactic medium with a high degree of collimation.

Summary
In this paper we have presented a detailed study of the jet collimation in the low-luminosity radio galaxy NGC 315, based on a multi-frequency VLBI and VLA data set. At 86 GHz we have imaged the innermost jet base on scales of only ∼160 R S . We have then compared our findings to those obtained for other nearby sources, considering a sample of 27 objects classified as low-excitation (LEG) or high-excitation (HEG) galaxies. This classification reflects a different nature of the accretion, hot and radiatively inefficient in the first case, cold and radiatively efficient in the second case. The results are summarized in the following.
-The jet collimation in NGC 315 is completed on sub-parsec scales. A transition from a parabolic to conical jet shape is detected at a de-projected distance z t = 0.58 ± 0.28 pc (or ∼5 × 10 3 R S ), which is much smaller than the Bondi radius estimated based on X-ray data, r B ∼92 pc. While most of the jets analyzed in the literature collimate on larger scales, a similar behavior is observed in other low-luminosity galaxies (NGC 1052, NGC 4261, Cen A), whose jets are freely expanding at a distance of less than one parsec from the black hole. If the transition to a conical jet shape is induced by a change in the external pressure profile, this must occur in the nuclear regions. An initial confinement from a thick disk extending out to ∼10 3 -10 4 R S is possible for such objects. -We have compared the mass-scaled expansion profiles of the jets in our sample. Most jets in radio galaxies show a parabolic shape while most jets in BL Lacs, which are probed on large scales on average, show a conical shape at their base. At the same de-projected distance from the black hole, HEGs present "thicker" jets, while most of the jets in LEGs (including NGC 315 and M87) expand following profiles which are well aligned with those observed in BL Lacs at larger distances. We suggest that, while both hot and cold disks can launch collimated winds, jets in HEGs are surrounded by more prominent outer sheaths, with an outer launch radius > 100 R S . On the contrary, jet sheaths launched by hot disks, e.g., by ADAFs, are anchored in the innermost disk regions, as measured in M 87. -Jet collimation in HEGs tends to proceed over larger scales (>10 6 R S ) than in LEGs. This result matches the findings of Potter & Cotter (2015) obtained for blazars, since jets in FSRQs have been shown to accelerate over a more extended region (>10 5 R S ) with respect to jets in BL Lacs. This supports the idea that, both in high-luminosity and in lowluminosity sources, the jet acceleration and collimation processes are co-spatial, as expected based on theoretical models for magnetically-driven cold outflows (Komissarov et al. 2007;Tchekhovskoy et al. 2008;Lyubarsky 2009). -The possibility that relativistic jets are collimated by disk winds is discussed. The observation of more extended collimation zones in jets surrounded by thicker sheaths is in agreement with theoretical modeling describing the case of a Poynting-flux dominated jet confined by a wind layer (Globus & Levinson 2016). This mechanism requires the wind to carry a significant fraction (>10%) of the total jet power. According to recent simulations aimed at explaining the diverse kinematic behavior of VLBI knots in blazars of different powers (Hervet et al. 2017) this condition is verified in all cases, except for the weakest BL Lacs (HBLs). The Xray detection of ultra-fast outflows in most of the HEGs in our sample (Tombesi et al. 2010(Tombesi et al. , 2014Reynolds et al. 2015) supports the existence of disk winds originating at large radii and carrying significant kinetic power. Thus, particularly for high-luminosity sources, jet collimation by disk winds may be viable mechanism. -Motivated by the observation of mostly FRII morphologies among the sources presenting thicker jets, we have discussed the possible role of disk winds in the origin of the FRI/FRII dichotomy. A powerful sheath stabilizes the inner spine by shielding it against entrainment from the interstellar medium (Perucho 2012(Perucho , 2016, thus enabling the jet to reach the intergalactic medium with a high degree of collimation. This may explain the observed formation of FRII morphologies in HEGs but also in some LEGs -those producing sufficiently powerful winds -as well the rare occurrence of FRI morphologies in HEGs.