A transition from parabolic to conical shape as a common effect in nearby AGN jets

TL;DR

This study identifies a common transition from parabolic to conical shapes in nearby AGN jets, suggesting a universal effect linked to jet energetics and ambient medium conditions, based on VLBA observations of 367 AGN.

Contribution

It provides the first systematic search for shape transitions in a large AGN sample, revealing their prevalence and proposing a physical model for the transition mechanism.

Findings

10 out of 29 nearby AGN show shape transition

Break occurs at 10^5-10^6 gravitational radii

Jet shape transition linked to energy flux balance

Abstract

Observational studies of collimation in jets in active galactic nuclei (AGN) are a key to understanding their formation and acceleration processes. We have performed an automated search for jet shape transitions in a sample of 367 AGN using VLBA data at 15 GHz and 1.4 GHz. This search has found ten out of 29 nearby jets at redshifts z<0.07 with a transition from a parabolic to conical shape, while the full analyzed sample is dominated by distant AGN with a typical z about 1. The ten AGN are UGC00773, NGC1052, 3C111, 3C120, TXS0815-094, Mrk180, PKS1514+00, NGC6251, 3C371, and BL Lac. We conclude that the geometry transition may be a common effect in AGN jets. It can be observed only when sufficient linear resolution is obtained. Supplementing these results with previously reported shape breaks in the nearby AGN 1H0323+342 and M87, we estimate that the break occurs at 10^5-10^6 gravitational radii from the nucleus. We suggest that the jet shape transition happens when the bulk plasma kinetic energy flux becomes equal to the Poynting energy flux, while the ambient medium pressure is assumed to be governed by Bondi accretion. In general, the break point may not coincide with the Bondi radius. The observational data supports our model predictions on the jet acceleration and properties of the break point.

Tables

Table 1. Table 2: Derived best-fit parameters of the two fitted dependencies d = a 1 ​ ( r + r 0 ) k 1 𝑑 subscript 𝑎 1 superscript 𝑟 subscript 𝑟 0 subscript 𝑘 1 d=a_{1}(r+r_{0})^{k_{1}} and d = a 2 ​ ( r + r 1 ) k 2 𝑑 subscript 𝑎 2 superscript 𝑟 subscript 𝑟 1 subscript 𝑘 2 d=a_{2}(r+r_{1})^{k_{2}} before and after the jet break, respectively ( Figure 1 ). We used the VLBA data at 15 GHz only (band ‘U’) or 15 GHz and 1.4 GHz (band ‘UL’) between r min subscript 𝑟 min r_{\mathrm{min}} and r max subscript 𝑟 max r_{\mathrm{max}} distance from the apparent core. Note that all the values of r 𝑟 r are projected on the plane of the sky.

Source	Band	$r_{\min }$	$r_{\max }$	$a_1$		$r_0$		$k_1$
		(mas)	(mas)	(pc${}^{1-k_1}$)	(mas${}^{1-k_1}$)	(pc)	(mas)
0111$+$021	U	0.2	1.5	$0.179\pm 0.010$	$0.188\pm 0.011$	$0.143\pm 0.085$	$0.157\pm 0.093$	$0.497\pm 0.077$
0238$-$084	U	0.3	2.5	$0.078\pm 0.007$	$0.319\pm 0.046$	$0.074\pm 0.038$	$0.740\pm 0.380$	$0.391\pm 0.048$
0415$+$379	U	0.2	6.0	$0.305\pm 0.011$	$0.313\pm 0.011$	$0.042\pm 0.020$	$0.044\pm 0.021$	$0.468\pm 0.026$
0430$+$052	U	0.2	2.0	$0.202\pm 0.015$	$0.245\pm 0.020$	$0.122\pm 0.071$	$0.188\pm 0.109$	$0.556\pm 0.070$
0815$-$094	U	0.2	1.0	…	$0.294\pm 0.015$	…	$0.163\pm 0.048$	$0.527\pm 0.044$
1133$+$704	U	0.3	1.5	$0.437\pm 0.013$	$0.464\pm 0.014$	$0.061\pm 0.046$	$0.069\pm 0.052$	$0.528\pm 0.040$
1514$+$004	U	0.2	3.5	$0.171\pm 0.011$	$0.171\pm 0.011$	$0.189\pm 0.088$	$0.189\pm 0.088$	$0.564\pm 0.048$
1637$+$826	U	0.2	3.0	$0.155\pm 0.005$	$0.223\pm 0.010$	$0.098\pm 0.044$	$0.204\pm 0.092$	$0.506\pm 0.041$
1807$+$698	U	0.2	1.4	$0.207\pm 0.016$	$0.210\pm 0.016$	$0.130\pm 0.089$	$0.133\pm 0.091$	$0.388\pm 0.087$
2200$+$420	U	0.9	2.0	$0.505\pm 0.029$	$0.449\pm 0.027$	$0.087\pm 0.096$	$0.067\pm 0.074$	$0.537\pm 0.057$
Source	Band	$r_{\min }$	$r_{\max }$	$a_2$		$r_1$		$k_2$
		(mas)	(mas)	(pc${}^{1-k_2}$)	(mas${}^{1-k_2}$)	(pc)	(mas)
0111$+$021	U	1.5	14.0	$0.252\pm 0.031$	$0.254\pm 0.031$	$-1.126\pm 0.181$	$-1.237\pm 0.199$	$0.934\pm 0.054$
0238$-$084	U	2.5	14.1	$0.252\pm 0.021$	$0.228\pm 0.047$	$-0.069\pm 0.067$	$-0.690\pm 0.670$	$1.052\pm 0.081$
0415$+$379	UL	6.0	60.8	$0.123\pm 0.019$	$0.122\pm 0.019$	$-2.043\pm 0.100$	$-2.151\pm 0.105$	$1.175\pm 0.046$
0430$+$052	UL	2.0	121.9	$0.229\pm 0.022$	$0.216\pm 0.021$	$-0.500\pm 0.188$	$-0.769\pm 0.289$	$1.131\pm 0.027$
0815$-$094	U	1.0	14.3	…	$0.282\pm 0.033$	…	$-0.085\pm 0.097$	$1.032\pm 0.049$
1133$+$704	U	1.5	5.1	$0.921\pm 0.057$	$0.941\pm 0.058$	$-0.753\pm 0.083$	$-0.857\pm 0.094$	$0.828\pm 0.047$
1514$+$004	U	3.5	14.3	$0.185\pm 0.010$	$0.185\pm 0.010$	$-1.167\pm 0.102$	$-1.167\pm 0.102$	$0.886\pm 0.022$
1637$+$826	U	3.0	8.9	$0.175\pm 0.007$	$0.213\pm 0.010$	$-0.089\pm 0.048$	$-0.185\pm 0.100$	$0.730\pm 0.029$
1807$+$698	UL	1.4	85.2	$0.179\pm 0.018$	$0.179\pm 0.018$	$-0.120\pm 0.188$	$-0.122\pm 0.192$	$1.023\pm 0.025$
2200$+$420	UL	2.0	49.0	$0.433\pm 0.016$	$0.447\pm 0.016$	$-1.142\pm 0.093$	$-0.885\pm 0.072$	$1.124\pm 0.009$

Table 2. Table 3: Derived best-fit parameters of a single fit dependence d = a ​ ( r + r 0 ) k 𝑑 𝑎 superscript 𝑟 subscript 𝑟 0 𝑘 d=a(r+r_{0})^{k} for 319 AGN. Their k 𝑘 k values are presented in Figure 2 . We used the VLBA data at 15 GHz only (band ‘U’) or 15 GHz and 1.4 GHz (band ‘UL’) between r min subscript 𝑟 min r_{\mathrm{min}} and r max subscript 𝑟 max r_{\mathrm{max}} distance from the apparent core. Full table is available online only, the first seven rows are shown here for guidance.

Source	Band	$r_{\min }$	$r_{\max }$	$a$	$r_0$	$k$
		(mas)	(mas)	(mas1-k)	(mas)
0003$+$380	U	0.5	2.5	$0.213\pm 0.059$	$0.043\pm 0.285$	$1.064\pm 0.218$
0006$+$061	U	0.5	2.8	$0.245\pm 0.070$	$0.205\pm 0.246$	$1.378\pm 0.187$
0007$+$106	U	0.5	1.6	$0.356\pm 0.162$	$0.092\pm 0.259$	$2.110\pm 0.489$
0010$+$405	U	0.5	7.0	$0.220\pm 0.051$	$-0.087\pm 0.268$	$0.804\pm 0.190$
0011$+$189	U	0.5	11.7	$0.266\pm 0.059$	$0.204\pm 0.323$	$0.744\pm 0.177$
0015$-$054	U	0.5	12.4	$0.200\pm 0.049$	$0.333\pm 0.340$	$0.814\pm 0.193$
0016$+$731	UL	0.5	38.2	$0.675\pm 0.065$	$-0.230\pm 0.080$	$0.867\pm 0.030$

Table 3. Table 4: Derived parameters of the jet shape break for 10 AGN with addition of 0321+340 adopted from Hada et al. ( 2018 ) and 1228+126 from Nokhrina et al. ( 2019 ) . Columns are as follows: (1) source name (B1950); (2) jet width at the break in mas; (3) same as (2) but in pc; (4) projected distance of the break from the apparent core along the jet in mas; (5) projected distance of the break from the BH along the jet in mas; (6) same as (5) but in pc; (7) deprojected distance of the break from the BH along the jet in pc, the parameter uses the estimated viewing angle; (8) presence of a bright low pattern speed feature (Lister et al., 2019 ) .

Source	$d_{\mathrm{break}}$		$r_{\mathrm{break},\mathrm{app}}^{\mathrm{proj}}$	$r_{\mathrm{break}}^{\mathrm{proj}}$	$r_{\mathrm{break}}^{\mathrm{proj}}$	$r_{\mathrm{break}}^{\mathrm{deproj}}$	Stationary
	(mas)	(pc)	(mas)	(mas)	(pc)	(pc)	jet feature
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)
0111$+$021	$0.30\pm 0.03$	$0.28\pm 0.03$	$2.46\pm 0.27$	$2.62\pm 0.29$	$2.38\pm 0.26$	$27.31$	Y
0238$-$084	$0.53\pm 0.05$	$0.05\pm 0.01$	$2.93\pm 0.57$	$3.73\pm 0.65$	$0.37\pm 0.06$	$0.49$	Y
0321$+$340	1	1.16	10	10.04	11.64	$106.07$	Ya
0415$+$379	$0.78\pm 0.03$	$0.74\pm 0.03$	$7.03\pm 0.50$	$7.07\pm 0.50$	$6.72\pm 0.47$	$29.00$	…
0430$+$052	$0.45\pm 0.06$	$0.29\pm 0.04$	$2.67\pm 0.40$	$2.85\pm 0.41$	$1.85\pm 0.27$	$5.77$	…
0815$-$094	$0.37\pm 0.05$	…	$1.37\pm 0.30$	$1.54\pm 0.30$	…	…	…
1133$+$704	$0.57\pm 0.02$	$0.50\pm 0.02$	$1.39\pm 0.09$	$1.46\pm 0.10$	$1.29\pm 0.09$	$14.80$	…
1228$+$126	$13.00\pm 0.50$	$1.20\pm 0.04$	…	$131\pm 6$	$10.50\pm 0.46$	43.41	Yb
1514$+$004	$0.34\pm 0.02$	$0.34\pm 0.02$	$3.10\pm 0.22$	$3.39\pm 0.30$	$3.39\pm 0.30$	$13.10$	Y
1637$+$826	$0.32\pm 0.02$	$0.16\pm 0.01$	$1.92\pm 0.26$	$2.13\pm 0.28$	$1.02\pm 0.13$	$3.30$	Y
1807$+$698	$0.26\pm 0.04$	$0.25\pm 0.04$	$1.53\pm 0.33$	$1.67\pm 0.34$	$1.63\pm 0.33$	$12.83$	Y
2200$+$420	$0.74\pm 0.03$	$0.95\pm 0.04$	$2.45\pm 0.10$	$2.52\pm 0.13$	$3.25\pm 0.16$	$24.57$	…

Table 4. Table 5: Angular size of the VLBA core at 15 GHz, d c MC superscript subscript 𝑑 c MC d_{\mathrm{c}}^{\mathrm{MC}} , and its offset from the true jet origin, r 0 MC superscript subscript 𝑟 0 MC r_{0}^{\text{MC}} , derived from our Monte Carlo modeling of the jet width compared with independent MOJAVE core size measurements in the visibility plane, d c u ​ v superscript subscript 𝑑 c 𝑢 𝑣 d_{\mathrm{c}}^{uv} (Lister et al., 2019 ) , and the core offset, r 0 cs superscript subscript 𝑟 0 cs r_{0}^{\text{cs}} , estimated from the multi-frequency core shift measurements (Pushkarev et al., 2012a ) . The shown d c u ​ v superscript subscript 𝑑 c 𝑢 𝑣 d_{\mathrm{c}}^{uv} values are medians over all epochs from Lister et al. ( 2019 ) .

Source	$d_c^{\mathrm{MC}}$	$d_c^{uv}$	$r_0^{\text{MC}}$	$r_0^{\text{cs}}$
	(mas)	(mas)	(mas)	(mas)
(1)	(2)	(3)	(4)	(5)
0111$+$021	$0.075\pm 0.025$	$0.079\pm 0.030$	$0.157\pm 0.093$	$0.159\pm 0.050$
0238$-$084	$0.282\pm 0.067$	$0.284\pm 0.042$	$0.740\pm 0.380$	…
0415$+$379	$0.073\pm 0.017$	$0.075\pm 0.008$	$0.044\pm 0.021$	$0.275\pm 0.050$
0430$+$052	$0.096\pm 0.034$	$0.182\pm 0.014$	$0.188\pm 0.109$	$0.051\pm 0.050$
0815$-$094	$0.113\pm 0.020$	$0.062\pm 0.041$	$0.163\pm 0.048$	…
1133$+$704	$0.116\pm 0.049$	$0.089\pm 0.048$	$0.072\pm 0.056$	…
1514$+$004	$0.067\pm 0.018$	$0.043\pm 0.027$	$0.189\pm 0.088$	…
1637$+$826	$0.100\pm 0.025$	$0.069\pm 0.004$	$0.204\pm 0.092$	$0.198\pm 0.050$
1807$+$698	$0.096\pm 0.031$	$0.067\pm 0.007$	$0.133\pm 0.091$	$0.240\pm 0.050$
2200$+$420	$0.105\pm 0.064$	$0.044\pm 0.003$	$0.067\pm 0.074$	$0.090\pm 0.050$

Table 5. Table 1: Properties for 12 sources with a detected jet shape break from this study ( Figure 1 ) as well as Hada et al. ( 2018 , 1H 0323+342) and Asada & Nakamura ( 2012 , M 87) . They are supplemented by the MOJAVE-1 sources for which redshift values, Doppler factor estimates, and robust jet shape fits ( Table 3 ) are available. Columns are as follows: (1) B1950 name; (2) alias; (3) optical class, where Q = quasar, B = BL Lac, G = radio galaxy, N = Narrow Line Seyfert 1 (NLSy1); (4) redshift; (5) literature reference for the data in column (4); (6) maximum apparent radial speed from Lister et al. ( 2019 ) , (7) variability Doppler factor from Hovatta et al. ( 2009 ) ; (8) viewing angle; (9) black hole mass estimated basing on assumption of virialized broad lines region (BLR) movement and correlation between the size of BLR and UV/optical luminosity; (10) literature reference for the data in column (9); (11) black hole mass estimated by a stellar velocity dispersion method and associated fundamental plane method (for 2200 + + 420); (12) literature reference for the data in column (11). Names of the sources with the shape break are highlighted by the boldface font.

Source	Alias	Opt.	$z$	Reference	${\beta }_{\mathrm{app}}$	$\delta$	$\theta$	$M_{BH}$	Reference	$M_{BH}$	Reference
		ID			(c)		(deg)	($\log \left(M_{☉}\right)$)		($\log \left(M_{☉}\right)$)
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)	(11)	(12)
0016$+$731	S5 0016$+$73	Q	1.781	Lawrence et al. (1986)	8.22	7.8	7.4	9.18	Torrealba et al. (2012)	…	…
0059$+$581	TXS 0059$+$581	Q	0.644	Sowards-Emmerd et al. (2005)	6.84	21.9a	1.5	…	…	…	…
0106$+$013	4C $+$01.02	Q	2.099	Hewett et al. (1995)	10.91	18.2	2.8	…	…	…	…
0111$+$021	UGC 00773	B	0.047	Wills & Wills (1976)	0.17	…	5.0b	…	…	…	…
0133$+$476	DA 55	Q	0.859	Lawrence et al. (1986)	16.53	20.5	2.7	9.07	Torrealba et al. (2012)	…	…
0215$+$015	OD 026	Q	1.715	Boisse & Bergeron (1988)	15.83	27.1a	1.8	…	…	…	…
0234$+$285	4C $+$28.07	Q	1.206	Shaw et al. (2012)	24.58	16.0	3.3	9.22	Shaw et al. (2012)	…	…
0238$-$084	NGC 1052	G	0.005	Denicoló et al. (2005)	0.42	0.3a	49.0	5.51	Torrealba et al. (2012)	8.19	Woo & Urry (2002)
0321$+$340	1H 0323$+$342	N	0.061	Marcha et al. (1996)	9.02	…	6.3	7.30	Landt et al. (2017)	…	…
0336$-$019	CTA 26	Q	0.852	Wills & Lynds (1978)	24.45	17.2	3.1	8.74	Torrealba et al. (2012)	…	…
0415$+$379	3C 111	G	0.049	Eracleous & Halpern (2004)	8.15	2.0a	13.4	8.21	Torrealba et al. (2012)	…	…
0420$-$014	PKS 0420$-$01	Q	0.916	Jones et al. (2009)	2.36	19.7	0.7	8.98	Torrealba et al. (2012)	…	…
0430$+$052	3C 120	G	0.033	Michel & Huchra (1988)	5.27	2.1a	18.7	7.52	Torrealba et al. (2012)	8.13	Woo & Urry (2002)
0458$-$020	S3 0458$-$02	Q	2.286	Strittmatter et al. (1974)	6.07	15.7	2.5	…	…	…	…
0528$+$134	PKS 0528$+$134	Q	2.070	Hunter et al. (1993)	10.91	30.9	1.2	9.03	Palma et al. (2011)	…	…
0605$-$085	OC $-$010	Q	0.870	Shaw et al. (2012)	31.98	7.5	3.4	8.63	Shaw et al. (2012)	…	…
0642$+$449	OH 471	Q	3.396	Osmer et al. (1994)	8.53	10.6	5.3	9.12	Torrealba et al. (2012)	…	…
0735$+$178	OI 158	B	0.450	Nilsson et al. (2012)	5.04	4.5a	12.7	…	…	…	…
0736$+$017	OI 061	Q	0.189	Ho & Kim (2009)	11.89	8.5	6.4	8.10	Torrealba et al. (2012)	…	…
0754$+$100	PKS 0754$+$100	B	0.266	Carangelo et al. (2003)	1.06	5.5	4.0	…	…	…	…
0805$-$077	PKS 0805$-$07	Q	1.837	White et al. (1988)	39.71	14.9a	2.5	…	…	…	…
0815$-$094	TXS 0815$-$094	B	…	…	…	…	5.0b	…	…	…	…
0827$+$243	OJ 248	Q	0.942	Shaw et al. (2012)	19.81	13.0	4.0	8.77	Torrealba et al. (2012)	…	…
0851$+$202	OJ 287	B	0.306	Stickel et al. (1989)	6.59	16.8	2.3	…	…	…	…
0923$+$392	4C $+$39.25	Q	0.695	Abazajian et al. (2004)	2.76	4.3	12.4	8.75	Torrealba et al. (2012)	…	…
0945$+$408	4C $+$40.24	Q	1.249	Abazajian et al. (2004)	20.20	6.3	5.2	8.66	Torrealba et al. (2012)	…	…
1055$+$018	4C $+$01.28	Q	0.888	Shaw et al. (2012)	7.00	12.1	4.1	9.14	Torrealba et al. (2012)	…	…
1127$-$145	PKS 1127$-$14	Q	1.184	Wilkes (1986)	18.93	21.9a	2.6	9.03	Torrealba et al. (2012)	…	…
1133$+$704	Mrk 180	B	0.045	Falco et al. (1999)	…	…	5.0b	…	…	8.21	Woo & Urry (2002)
1156$+$295	4C $+$29.45	Q	0.725	Shaw et al. (2012)	24.59	28.2	2.0	8.35	Torrealba et al. (2012)	…	…
1228$+$126	M 87	G	0.004	Smith et al. (2000)	…	…	14.0	…	…	9.82	Gebhardt et al. (2011)
1253$-$055	3C 279	Q	0.536	Marziani et al. (1996)	20.58	23.8	2.4	9.52	Torrealba et al. (2012)	…	…
1308$+$326	OP 313	Q	0.997	Shaw et al. (2012)	21.30	15.3	3.5	8.76	Torrealba et al. (2012)	…	…
1413$+$135	PKS 1413$+$135	B	0.247	Stocke et al. (1992)	1.20	12.1	0.9	…	…	…	…
1502$+$106	OR 103	Q	1.839	Adelman-McCarthy et al. (2008)	18.25	11.9	4.4	9.44	Torrealba et al. (2012)	…	…
1510$-$089	PKS 1510$-$08	Q	0.360	Thompson et al. (1990)	21.56	16.5	3.4	8.36	Torrealba et al. (2012)	…	…
1514$+$004	PKS 1514$+$00	G	0.052	Healey et al. (2008)	0.45	…	15.0c	…	…	…	…
1538$+$149	4C $+$14.60	B	0.606	Shaw et al. (2013)	8.75	4.3	10.5	…	…	…	…
1611$+$343	DA 406	Q	1.400	Shaw et al. (2012)	31.08	13.6	3.1	9.00	Torrealba et al. (2012)	…	…
1633$+$382	4C $+$38.41	Q	1.813	Shaw et al. (2012)	29.22	21.3	2.6	9.00	Torrealba et al. (2012)	…	…
1637$+$574	OS 562	Q	0.751	Marziani et al. (1996)	13.59	13.9	4.1	8.72	Torrealba et al. (2012)	…	…
1637$+$826	NGC 6251	G	0.024	Wegner et al. (2003)	…	…	18.0	…	…	8.78	Ferrarese & Ford (1999)
1641$+$399	3C 345	Q	0.593	Marziani et al. (1996)	19.28	7.7	5.1	9.00	Torrealba et al. (2012)	…	…
1730$-$130	NRAO 530	Q	0.902	Junkkarinen (1984)	13.73	10.6	5.2	…	…	…	…
1749$+$096	OT 081	B	0.322	Stickel et al. (1988)	7.91	11.9	4.5	…	…	…	…
1803$+$784	S5 1803$+$784	B	0.680	Lawrence et al. (1996)	1.16	12.1	0.9	9.20	Torrealba et al. (2012)	…	…
1807$+$698	3C 371	B	0.051	de Grijp et al. (1992)	0.01	1.1	7.3	7.14	Torrealba et al. (2012)	8.51	Woo & Urry (2002)
1828$+$487	3C 380	Q	0.692	Lawrence et al. (1996)	13.03	5.6	7.4	8.59	Torrealba et al. (2012)	…	…
1849$+$670	S4 1849$+$67	Q	0.657	Stickel & Kuehr (1993)	21.43	8.1a	4.7	8.83	Torrealba et al. (2012)	…	…
1928$+$738	4C $+$73.18	Q	0.302	Lawrence et al. (1986)	7.55	1.9	14.2	8.62	Torrealba et al. (2012)	…	…
2136$+$141	OX 161	Q	2.427	Wills & Wills (1974)	2.58	8.2	4.0	…	…	…	…
2145$+$067	4C $+$06.69	Q	0.999	Steidel & Sargent (1991)	3.09	15.5	1.4	9.28	Torrealba et al. (2012)	…	…
2200$+$420	BL Lac	B	0.069	Vermeulen et al. (1995)	9.96	7.2	7.6	…	…	8.23	Woo & Urry (2002)
2201$+$171	PKS 2201$+$171	Q	1.076	Smith et al. (1977)	17.56	10.0a	4.9	…	…	…	…
2201$+$315	4C $+$31.63	Q	0.295	Marziani et al. (1996)	7.99	6.6	8.5	8.78	Torrealba et al. (2012)	…	…
2223$-$052	3C 446	Q	1.404	Wright et al. (1983)	17.74	15.9	3.6	…	…	…	…
2227$-$088	PHL 5225	Q	1.559	Abazajian et al. (2004)	3.07	15.8	1.4	8.90	Torrealba et al. (2012)	…	…
2230$+$114	CTA 102	Q	1.037	Falomo et al. (1994)	17.73	15.5	3.7	…	…	…	…
2251$+$158	3C 454.3	Q	0.859	Jackson & Browne (1991)	3.46	32.9	0.4	9.07	Torrealba et al. (2012)	…	…
a Doppler factor value is from Liodakis et al. (2017).
b Assumed $\theta$ value as typical for BL Lacs.
c Assumed $\theta$ value as typical for radio galaxies in the list which do not show a strong counter-jet.

Equations

$$d=\{\begin{array}{ll}a(r+r_0{)}^k+\mathcal{N}(0,{\sigma }_1^2),& \text{mbox}ifr<R\\ b+\mathcal{N}(0,{\sigma }_2^2),& \text{mbox}ifr>R\end{array}$$

$$\theta =\arctan \frac{2{\beta }_{\text{app}}}{{\beta }_{\text{app}}^2+{\delta }_{\text{var}}^2-1}$$

$${\sigma }_{\mathrm{M}}=\frac{{\Omega }_0^2{\Psi }_0}{8{\pi }^2\mu \eta c^2}$$

$$\gamma (r_{\perp })=\frac{r_{\perp }}{R_{\mathrm{L}}},$$

$$w=\frac{c_{\mathrm{jet}}^2}{(\Gamma -1)}{\left(\frac{n}{n_{\mathrm{jet}}}\right)}^{\Gamma -1},$$

$${\mathcal{M}}^{2}n=4\pi {\eta }^2{m}_{\mathrm{p}}{c}^2\left[1+\frac{1}{\Gamma -1}\frac{c_{\mathrm{jet}}^2}{c^2}{\left(\frac{n}{n_{\mathrm{jet}}}\right)}^{\Gamma -1}\right].$$

$${P|}_{r_{\perp }=d/2-0}=P_{\mathrm{ext}}.$$

$$\frac{d}{d{r}_{\perp }}\left(P+\frac{B^2}{8\pi }\right)=0.$$

$$p_0={\left(\frac{{\Psi }_0}{2\pi R_{\mathrm{L}}^2{\sigma }_{\mathrm{M}}}\right)}^2,$$

$$P\propto d^{-3.7}$$

$$P\propto d^{-2.4}$$

$$P_{\mathrm{ext}}=P_0{\left(\frac{r}{r_0}\right)}^{-b}.$$

$$d\propto r^{0.54}.$$

$$d\propto r^{0.83}.$$

$$\mathbf{S}=\frac{c}{4\pi }\mathbf{E}\times \mathbf{B}$$

$$\mathbf{K}=\gamma m{c}^{2}n{\mathbf{u}}_{\mathrm{p}},$$

$$\sigma =\frac{|\mathbf{S}|}{|\mathbf{K}|}=\frac{{\Omega }_{\mathrm{F}}I}{2\pi c\gamma \mu \eta }.$$

$$\sigma ={\Omega }_{\mathrm{F}}\frac{L-{\Omega }_{\mathrm{F}}{r}_{\perp }^{2}E/c^2}{E-{\Omega }_{\mathrm{F}}L-{\mathcal{M}}^{2}E}.$$

$$|\mathbf{S}|\propto I=\pi j{r}_{\perp }^2+o(r_{\perp }^2)$$

Full text

††pubyear: 2020††pagerange: A transition from parabolic to conical shape as a common effect in nearby AGN jets–1

1 Introduction

Understanding the physical processes that determine the formation, acceleration and collimation of relativistic jets in active galactic nuclei (AGN) continues to be among the most challenging problems of modern astrophysics. There are a wide variety of analytical and numerical models for jet acceleration and its confinement (e.g., Vlahakis & Königl, 2003; Beskin & Nokhrina, 2006; McKinney, 2006; Komissarov et al., 2007; Tchekhovskoy et al., 2011; McKinney et al., 2012; Potter & Cotter, 2015) that consider different solutions for jet shapes, such as cylindrical, conical and parabolic. General relativistic magnetohydrodynamic simulations (e.g., McKinney et al., 2012) predict that a jet starting from its apex has a parabolic streamline within the magnetically dominated acceleration zone. At other scales it transitions to a conical geometry associated with equipartition between energy densities of the magnetic field and the radiating particle populations. It has been shown for cold jets that acceleration should not occur in a conical jet. This requires something akin to a parabolic jet shape closer to the jet base to allow differential expansion (Vlahakis & Königl, 2004; Komissarov, 2012).

In order to investigate these theories it is important to collect observational data on jet profile shapes for a large enough sample of AGN whose properties are well understood. The first observational evidence for a transition from parabolic to conical jet shape was detected in M87 (Asada & Nakamura, 2012) at a distance of about 900 mas near the feature HST-1, about 70 pc in projection, corresponding to $10^{5}$ Schwarzschild radii. A few more studies of nearby AGNs to probe their innermost jet regions were performed recently: Mkn 501 (Giroletti et al., 2008), Centaurus A (Müller et al., 2014), Cygnus A (Boccardi et al., 2016; Nakahara et al., 2019), NGC 6251 (Tseng et al., 2016), 1H 0323+342 (Hada et al., 2018), 3C 273 (Akiyama et al., 2018), NGC 4261 (Nakahara et al., 2018), 3C 84 (Giovannini et al., 2018), 3C 264 (Boccardi et al., 2019), NGC 1052 (Nakahara et al., 2020). Hovatta et al. (2019) have indirectly addressed this question for the 3C 273 jet close to the central engine on the basis of a model analysis of ALMA rotation measure data. Larger survey studies (Pushkarev et al., 2009) have typically probed regions farther away from the central nucleus, although Algaba et al. (2017) have used apparent parsec-scale jet base parameters closer in.

In a previous work (Pushkarev et al., 2017), we analyzed parsec-scale radio VLBI images of jets in 362 active galaxies from the MOJAVE program (Lister et al., 2018). This sample is dominated by compact radio bright blazars with a jet at a small angle to the line of sight and a typical redshift $z\approx 1$. However, some low luminosity nearby radio galaxies were also included. Pushkarev et al. (2017) show that while the majority of resolved jets have a shape close to conical, a significant fraction of the sample has observed deviations. A systematic change in jet width profile has been noted by Hervet et al. (2017), who explain it by using a stratified jet model with a fast spine and slow but relatively powerful outer layer. In this paper, we investigate if this outcome is partly affected by the typical finite angular resolution of VLBI observations. We probe a possible dependence of the jet shape on the distance $r$ from the nucleus. Furthermore, we perform a systematic search for a possible transition from one jet shape to another on the basis of 15 GHz and 1.4 GHz VLBA images.

The observation of jets with a change from parabolic to conical shape may provide an instrument to probe the MHD acceleration mechanism models as well as the ambient medium conditions. The change in jet shape in M87 (Asada & Nakamura, 2012) is coincident with the stationary bright feature HST-1, which can be associated with the change in ambient pressure profile and appearance of a recollimation shock due to pressure drop and abrupt expansion. This interpretation is supported by the measurements of external medium pressure by Russell et al. (2015) almost down to the Bondi radius $r_{\mathrm{B}}=2GM/c_{\mathrm{s}}^{2}$ (sphere of influence), with an observed mass density profile $\rho\propto r^{-1}$ (here $c_{\mathrm{s}}$ is a sound speed). The recently observed jet shape in 1H 0323+342 (Hada et al., 2018) demonstrates a similar behavior. On the other hand, there are models predicting a jet shape transition for a single power law pressure profile. The analytical model by Lyubarsky (2009) predicts the transition from parabolic to conical form for certain regimes, as well as quasi-oscillations in jet shape in the conical domain. This solution has been applied to the reconstruction of the recollimation shock properties of M87 by Levinson (2017), with a predicted total jet power on the order of $10^{43}$ erg/s. The recent semi-analytical results for the warm jet matching the ambient medium with a total electric current closed inside a jet by Beskin et al. (2017) predicts a change in a jet shape from parabolic to conical for the Bondi pressure profile $P\propto r^{-2}$. In this work we follow the latter model and consider the results for a warm outflow in more detail.

The structure of the paper is the following: section 2 presents our results of a search for the jet profile change from parabolic to conical in a large sample of AGN jets, we suggest a model and interpret our findings in section 3, a discussion is presented in section 4. We summarize our work in section 5. Throughout this paper we will use the term “core” as the apparent origin of AGN jets that commonly appears as the brightest feature in VLBI images of blazars (e.g., Lobanov, 1998; Marscher, 2008). We adopt a cosmology with $\Omega_{m}=0.27$, $\Omega_{\Lambda}=0.73$ and $H_{0}=71$ km s*-1* Mpc*-1* (Komatsu et al., 2009).

2 A discovery of shape transition as a common effect in AGN jets

2.1 Automated search of candidates with a change in jet geometry

For the purposes of our study, we made use of data at 15 GHz from the MOJAVE program, the 2 cm VLBA Survey, and the National Radio Astronomy Observatory (NRAO) data archive (Lister et al., 2018) for those sources that have at least five VLBA observing epochs at 15 GHz between 1994 August 31 and 2016 December 26 inclusive. We used the 15 GHz VLBA total intensity MOJAVE stacked epoch images supplemented by single epoch 1.4 GHz VLBA images to derive apparent jet widths, $d$, as a function of projected distance $r$ from the jet core, and determined jet shapes similar to Pushkarev et al. (2017). In that work we fitted the $d$–$r$ dependence with a single power law $d\propto r^{k}$. The index is expected to be $k\approx 0.5$ for a quasi-parabolic shape and 1.0 for a conical jet. We note that even single-epoch observations at 1.4 GHz adequately reproduce source morphology, i.e., effectively fill jet cross-section due to a steep spectrum of synchrotron emission of the outflow, with a typical spectral index $-0.7$ measured between 2 and 8 GHz (Pushkarev & Kovalev, 2012) and $-1.0$ between 8 and 15 GHz (Hovatta et al., 2014), making the low-frequency observations sensitive enough to probe jet morphology at larger scales. In our analysis we use the jet width measurements made at 1.4 GHz only on large scales, not covered by the 15 GHz data. These scales are typically beyond 10 mas. This allows us to neglect the core shift effect (e.g., Kovalev et al., 2008; Sokolovsky et al., 2011; Plavin et al., 2019a), which is expected to be about 1 mas between 1.4 and 15 GHz. Its value can not be easily derived since it requires simultaneous observations at different frequencies. As a result, the jet widths estimated at 15 GHz smoothly transition to those at 1.4 GHz.

We have carried out a similar analysis allowing for a change in the jet shape. Using all available data (15 GHz only or combined data set at 15 GHz and 1.4 GHz) for each source, we performed a double power law fit of the jet width as a function of distance, dividing the jet path length in a logarithmic scale by two parts in proportion of 1:1, 1:2 and 1:3 to search for cases when the fitted $k$-index at inner scales was $0.5\pm 0.2$, while at outer scales it was $1.0\pm 0.2$. After such cases were identified automatically, we tuned the fits by setting the distance of the transition region by eye.

We ended up dropping 36 AGN jets from the original samle of 367 objects as having unsatisfactory fits caused by either (i) non-optimal ridge line reconstruction for jets with strong bending, (ii) numerous large gaps in jet emission, (iii) too short a jet length (iv) low intensity regions not captured well by our jet width fitting. This resulted in a sample comprising 331 AGN jets.

As a result of this analysis, we found a shape transition in ten jets (Figure 1, LABEL:t:sourprop) out of 367 analyzed. We emphasize that all the AGNs with detected transition of the jet shape turned out to have low redshifts $z<0.07$, i.e., have a high linear resolution of 15 GHz VLBA observations — better than 1 pc. This is highly unlikely to occur by chance and provides additional strong evidence that this result is not an observational artifact but a real effect. See discussion of the rest of analyzed low redshift AGN in the sample in subsection 2.4. Among the ten sources, there is one, the radio galaxy 0238$-$084 (NGC 1052), that shows a two-sided jet morphology. For this object, we analyzed the approaching, brighter outflow propagating to north-east direction, determining the position of a virtual VLBI core using a kinematic-based minimization method described in Vermeulen et al. (2003).

Following our discovery of the shape transition preferentially occurring in nearby AGN, we supplemented our initial AGN sample of 362 targets from Pushkarev et al. (2017) with stacked images of five more low-$z$ AGN which had five or more 15 GHz VLBA observing epochs after the Pushkarev et al. analysis was finished. These were: 0615$-$172, 1133+704, 1200+608, 1216+061, 1741+196. All the stacked images are available from the MOJAVE database111http://www.physics.purdue.edu/MOJAVE/.

2.2 Rigorous fitting of the jet shape

For each of the 10 sources found to have a jet geometry transition, we fit the data with the following dependencies: $d=a_{1}(r+r_{0})^{k_{1}}$ and $d=a_{2}(r+r_{1})^{k_{2}}$, describing a jet shape before and after the break. Here $r_{0}$ is understood as the separation of the 15 GHz apparent core from the true jet origin due to the synchrotron opacity (e.g., Lobanov, 1998; Pushkarev et al., 2012a), while $r_{1}$ shows how much one underestimates the jet length if it is derived from the data only beyond the geometrical transition of the jet. We note that this approach is more accurate but more computationally intensive than that used by Pushkarev et al. (2017) and applied in the original selection of jet break candidates. It is needed in order to better fit for jet shape close to the apex.

We fit these dependencies with Bayesian modeling using the NUTS Markov Chain Monte Carlo sampler based on the gradient of the log posterior density. It was implemented in PYMC3 (Salvatier et al., 2016), which automatically accounts for uncertainties of all the parameters in further inferences. The best fit parameters are listed in Table 2, showing that initially the jets are quasi-parabolic with $k_{1}$ close to 0.5, while beyond the break point region the outflow manifests a streamline close to conical, with $k_{2}\approx 1$. The location of the jet shape break given in Table 4 is estimated as the intersection point of these two $d-r$ dependencies. Note that Table 4 includes results on the jet shape transition region for two more sources, 1H 0323+342 and M87 taken from Hada et al. (2018) and Nokhrina et al. (2019), respectively. We also note that the shown error of the deprojected position of the break is propagated from the fitting procedure, it does not include uncertainties on the viewing angle and black hole mass.

For other sources without a detected shape break, we fit a single power-law $d=a(r+r_{0})^{k}$ for consistency. We excluded objects with unreliable ridge line detection or patchy structure in images (15 sources) and those with nonphysical $d-r$ dependence (24 sources) after visual inspection. They constitute only about one tenth of the dataset and thus the exclusion should not bias our estimates. To account for increased uncertainties of jet width measurements further from the core, the power law model is complemented as following:

[TABLE]

Here all of $a,r_{0},k,R,\sigma_{1},b,\sigma_{2}$ are treated as unknown parameters and inferred simultaneously using a Nested Sampling algorithm as implemented in PolyChord (Handley et al., 2015). As expected, $\sigma_{2}$ is typically significantly larger than $\sigma_{1}$. We find that this model generally captures the $d-r$ dependence and its uncertainty well. Fitting results are given in Table 3 and the source distribution of exponents $k$ is shown in Figure 2. Even though the estimates for individual sources have a large spread, the median exponent is very close to $1$. This indicates a conical average outflow shape, and agrees with previous results using slightly different estimation method (Pushkarev et al., 2012b). We note the peak in the histogram bin at $k=0.5$ which corresponds to the parabolic jet shape; the number of objects with $k\approx 0.5$ is not high enough in the sample to make it significant (Table 3).

2.3 Checking consistency of the fits and analyzing for possible biases

By setting $r=0$ we can estimate the apparent core size $d_{\mathrm{c}}^{\mathrm{MC}}$ at 15 GHz from the Monte-Carlo fit of the jet width as $a_{1}r_{0}^{k_{1}}$ and compare it with a median value of the core size $d_{\mathrm{c}}^{uv}$ derived from structure modelfit in visibility plane taken from Lister et al. (2019), see Table 5. Two sources, 0111+021 and 0415+379, show a good agreement between $d_{\mathrm{c}}^{\mathrm{MC}}$ and $d_{\mathrm{c}}^{uv}$, while for the other seven objects $d_{\mathrm{c}}^{\mathrm{MC}}$ is somewhat larger than $d_{\mathrm{c}}^{uv}$. This is likely due to a non-ideal determination of the core position throughout the epochs, which is used to align single-epoch maps to produce stacked images. The radio galaxy 0430+052 is the only source having $d_{\mathrm{c}}^{\mathrm{MC}}<d_{\mathrm{c}}^{uv}$ for reasons that are unclear.

A bias related to this effect might affect the results. Statistically analyzing jet shapes for the whole sample of 331 sources with stacked VLBA images by introducing different ridge line path length limits we have found the following. A near-parabolic streamline for quasars and BL Lacs can be derived if the innermost jet, only up to $\sim$1 mas from the apparent core, is considered. This is not a real effect. The bias is found to be the most pronounced for curved jets or jets with features emerging at different position angles over time (Lister et al., 2013). This is confirmed by an apparent artificial correlation of median jet width with the number of epochs in a stacked image for such AGN. Uncertainties in the core position also contribute to this effect due to the imperfect alignment of images while performing the stacking. Variability of opacity conditions and apparent position of the core (Plavin et al., 2019a) affect this partially even though the alignment of the stacked single epoch images is done on the core position. Together, it causes an additional artificial widening of the jet near the core region up to distances $r\approx 0.3$ mas. The effect quickly vanishes at larger scales. Thus, if we exclude jet width measurements at distances $\lesssim 0.4$ mas, the effect becomes much weaker and disappears completely if we rule out the measurements within 0.5 mas from the core. We also note that radio galaxies, being at low redshift and thus having apparently wider outflows, are much less subject to this effect. The same is true for the sources with a jet shape break shown in Figure 1, as these are low-redshift objects. Only for BL Lac, as the most remote source among them and also having a bright quasi-stationary component near the core (Cohen et al., 2014), we put a conservative limit of 0.9 mas. For the other sources we used the non-cut intervals listed in Table 2, because dropping measurements at $r<0.5$ mas did not significantly change the fit parameters. For the remaining sources, we have dropped all measurement for $r<0.5$ mas while analyzing the data (Table 3).

Another possible problem might be related to cases where the jet width is completely unresolved. Indeed, this was found for some AGN targets at some epochs from the visibility model fitting of the core (e.g., Kovalev et al., 2005; Lister et al., 2019). We have addressed this issue by dropping all measurement for $r<0.5$ mas. Interestingly, the rest of the measured deconvolved jet width values are always positive. If we assume that this is some sort of a positive bias overestimating the width, it should not depend on $r$ for unresolved jets and will result in $k$ values close to zero. This behavior was not seen in our fitting results.

We have also compared the fitted parameter $r_{0}$ with the core offset from the jet base estimated from the core shift measured between 15 GHz and 8 GHz (Pushkarev et al., 2012a) assuming an inverse frequency dependence $r\propto\nu^{-1}$. These quantities, also listed in Table 5, agree well within the errors in four out of six sources having measured core shifts. The large discrepancy for two sources can be explained by the recently recently established phenomenon of significant core shift variability (Plavin et al., 2019a) or the difference between the true jet shape derived by us and the assumed conical jet shape in (Pushkarev et al., 2012a). We note that this result opens a new way to estimate the distance to the true jet origin which does not require an assumption regarding the jet geometry.

We checked and complemented our analysis using 43 GHz data from the Boston University (BU) AGN group222https://www.bu.edu/blazars/VLBAproject.html for 0430+052 and 2200+420 (Figure 1), which are present in both the MOJAVE and BU samples. For each of these sources we (i) produced stacked total intensity 43 GHz maps, aligning single epoch-images by the position of the VLBA core derived from structure modelfitting of the visibility data, (ii) determined the reconstructed jet ridge line, and (iii) fitted the transverse jet width as a function of distance from the core (Figure 3). It resulted in the same $k$-values before and after the break as in our original analysis within the errors (compare Figure 3 and Table 2). The jet shape transition region is found at core separations comparable to those from the 15 GHz data fits but has shifted slightly. We note however that 7 mm jet width estimates are systematically lower than those found from the 15 GHz data due to the weak high frequency synchrotron emission coming from the jet edges. Robust estimates of jet geometry and particularly of the jet width require high dynamic range images which are better sampled at intermediate radio frequencies. A good agreement between 15 GHz and 1.4 GHz width measurements increases the robustness of our results.

We warn readers about deriving jet shapes from structure model fitting of single-epoch data (e.g., Hervet et al., 2017), as the jet may appear quasi-parabolic ($k<1$) up to a certain (typically short) distance from the core and then change its shape to conical ($k\approx 1$). This effect occurs in the sources that show variations in their inner jet position angle. Lister et al. (2013) established this as a common, decade-timescale phenomenon for the most heavily monitored AGNs in the MOJAVE sample. Thus, single-epoch VLBI maps may not reveal the whole jet cross-section, but rather a portion of it, especially in the inner jet regions where images are dynamic range limited. Therefore, the conclusions regarding jet geometry based strictly on a modelfit approach should be treated with caution.

2.4 Jet shape transition: a common effect in AGN jets, its consequences and prospects

We have found evidence for geometry transition in many jets for which sufficient linear resolution was achieved. This means that a change in jet shape is a common phenomenon which has significant consequences for many high angular resolution astrophysical and astrometric studies. It is difficult to conclude if the geometry transition with measured properties is specific to only nearby radio galaxies and BL Lacs, or can be extended to the AGN class in general. The radio luminosities of the nearby ($z<0.07$) jets are much lower than the rest of the sample and this might affect the geometry and transition zone. We note that Figure 4 presents a consistent picture of the power index dependence on the downstream distance for nearby and distant jets.

In total, indications of the transition from parabolic to conical shape are found in 10 out of 29 nearby ($z<0.07$) jets observed as part of the MOJAVE program or by other investigators. VLBA archival data from the latter were processed by the MOJAVE team. The reasons for non-detection of a geometry transition in nearby AGN jets are varied. Some jets, e.g., 0007+106 and 1959+650, have too compact structure to study their shapes. Some others, e.g., 0241+622, 0316+413, 1216+061, show purely parabolic streamlines (Table 3), and their transition regions are expected at larger angular scales than those probed by our observations. E-MERLIN or low-frequency VLBA observations are needed. For example, the nearby radio galaxy 1216+061 ($z=0.0075$, scale factor 0.15 pc mas*-1*; not shown in Figure 4) has a parabolic streamline with $k=0.64\pm 0.05$ out to 7 mas at 15 GHz, corresponding to a deprojected distance of only $\approx 1$ pc. We are studying the remaining 12 low-redshift jets that show no sign of a profile break in a followup approved VLBA program.

The other jets in the sample (Table 3), namely 97 %, do not show a clear significant change in jet geometry. We explain this by (i) a large scale factor of the order of 8 pc mas*-1* for a typical source in the sample at a redshift of $z\sim 1$ and (ii) a small viewing angle typically about several degrees (Pushkarev et al., 2017). Jet power may also play a role, since the MOJAVE sample is flux-density limited and the AGN with $z>0.1$ typically have jet luminosities $\sim 2$ orders of magnitude higher than the lower-redshift ones. The jets with a detected shape change have an average scaling factor of 0.7 pc mas*-1* and, on average, larger viewing angle since 6 out of 12 are radio galaxies. Thus, if a transition region is located at a distance of a few tens of pc, it corresponds to a projected angular separation of $\la$1 mas from the apparent jet base at 15 GHz, which is comparable to the typical interferometric restoring beam size. VLBI observations at higher frequencies may be more effective in registering the jet shape transition, since they provide a better angular resolution and are less subject to opacity effects. This would probe scales closer to the jet apex and possible dependencies between acceleration zone extension and the maximum bulk Lorentz factor or jet power, as predicted by Potter & Cotter (2015). On the other hand, the steep spectrum of the optically thin jet emission hinders the tracing of the jet for long distances. The small viewing angles of the bright AGN jets set another limit on any jet shape investigation in the innermost parts. The streamline of an outflow can be studied down to distances at which the jet half-opening angle is still smaller than viewing angle. As shown by Pushkarev et al. (2017), the intrinsic jet opening angle reaches values of a few degrees at scales of the order of 10 pc. This suggests that the jet shape transition phenomenon might be more effectively studied for nearby AGNs that are oriented at larger angles to the line of sight. After considering all the points discussed above, we have begun a dedicated VLBA program in 2019 to search for geometry transitions in 61 AGN jets with $z<0.07$ from observations at 15 GHz and 1.4 GHz.

It is a challenging problem to estimate the consequences of this result on astrometry and astrophysics of AGN. VLBI astrometry delivers the position of the true jet apex only if the opacity driven core shift is proportional to the frequency as $r\propto\nu^{-1}$ (Porcas, 2009). However, this is expected only for conical jets and synchrotron opacity (Lobanov, 1998). A non-conical jet base results in an extension of the true jet length between the apex and the observed opaque core. This also produces somewhat larger VLBI-Gaia offsets for AGN positions (Kovalev et al., 2017; Plavin et al., 2019b) than predicted by Kovalev et al. (2008).

2.5 Deprojected position of the jet break

We chose the MOJAVE-1 sample of 135 AGN (Lister et al., 2009) to perform a direct comparison with the 12 jets showing the breaks. Our reasoning is as follows. Most of MOJAVE-1 targets were observed by VLBA not only at multiple 15 GHz epochs but also in a single epoch at 1.4 GHz, which increases the jet distance probed by our analysis. In addition, VLBI measurements of the apparent kinematics $\beta_{\mathrm{app}}$ (Lister et al., 2019) and variability Doppler factor estimates $\delta$ (Hovatta et al., 2009; Liodakis et al., 2017) are available for a large fraction of the sample. We need this information to derive deprojected distance values. These requirements result in a sample of 65 sources (LABEL:t:sourprop) described in Pushkarev et al. (2017).

We derived viewing angle estimates through the relation

[TABLE]

to convert the jet distance from angular projected to linear deprojected. Note that this assumes the same beaming parameters for the flux density variability and jet kinematics. For $\beta_{\mathrm{app}}$ we used the fastest non-accelerating apparent jet speeds from the MOJAVE kinematic analysis. For 1H 0323+342 we use $\theta=6.3\degr$, based on the observed superluminal motion (Lister et al., 2016) assuming $\theta=(1+\beta_{\mathrm{app}})^{-0.5}=\gamma^{-1}$, which minimizes the required bulk Lorentz factor $\gamma$. The other possible viewing angle value for this target $\theta=4\degr$ is based on the variability time scale (Hada et al., 2018). For the BL Lac objects 0111+021 and 1133+704 we assumed a viewing angle of $5\degr$, typical for this class of AGN (Hovatta et al., 2009; Savolainen et al., 2010; Pushkarev et al., 2017; Liodakis et al., 2017). For the radio galaxy 1514+004 we assumed a viewing angle of $15\degr$ which is typical for this class of AGN in our sample.

In Figure 4, we plot the corresponding single power-law $k$-index values derived from the 15 GHz and 1.4 GHz VLBA data (Pushkarev et al., 2017) versus deprojected distance from the 15 GHz VLBA core for 62 sources. There are eleven sources with known deprojected linear jet distance that have a jet shape transition (Figure 1, LABEL:t:sourprop). They are shown by a pair of points each from the double power-law fits. The BL Lac object 0815$-$094 is not shown in Fig. 4, as it does not have a measured spectroscopic redshift. Our results on jet shape transition (Table 2, Table 4, Figure 4) are supplemented by multi-frequency data for M87 from Nakamura et al. (2018), with $k_{1}=0.57$, $k_{2}=0.90$, and break point position obtained by Nokhrina et al. (2019). For M87 we adopt $\theta=14\degr$ (Wang & Zhou, 2009), consistent with more recent results by Mertens et al. (2016). For NLSy1 1H 0323+342 we use 1.4–2.3 GHz measurements from VLBA observations (Hada et al., 2018), with $k_{1}=0.6$ and $k_{2}=1.41$, for which the jet shape break point position is estimated.

Horizontal lines represent the scales at which $k$-indices were derived, starting from several tens of mas distance from the 15 GHz VLBA core (see subsection 2.3) and up to distances limited by the sensitivity of our observations. The nearby jets, for which we are probing closer to the central engine, have low $k$ values and show a transition from quasi-parabolic values at small scales to quasi-conical at larger scales (Figure 4). It is possible that at scales greater than $\sim 100$ kpc, where jets become diffuse and disruptive, their geometry further changes from conical to hyperbolic, characterized by more rapid expansion (Owen et al., 2000).

In order to plot the observed $k$-index values as a function of the deprojected distance along jets in gravitational radius $r_{\mathrm{g}}=GM/c^{2}$ units, we use the black hole masses estimated assuming virialized broad lines region (BLR) motion and correlation between BLR size and UV/optical luminosity (Torrealba et al., 2012; McLure & Jarvis, 2002; Vestergaard & Peterson, 2006; Landt et al., 2017; Palma et al., 2011; Shaw et al., 2012; Liu et al., 2006). We also use mass values inferred by stellar or gas kinematics methods (e.g., Woo & Urry, 2002) for the closest sources. The mass values and references can be found in LABEL:t:sourprop. We plot the data in Figure 5. It turns out that the sources with BH masses obtained by stellar velocity dispersion method or stellar/gas kinematics measurements are the subset of the sources with the detected jet shape break (i.e., the closest ones).

Since estimating the black hole mass is a complicated and strongly model-dependent method, some of the values might be significantly in error. By dropping the highest and lowest values as possible outliers of the derived jet break position $r_{\mathrm{break}}$ measured in $r_{\mathrm{g}}$ we are able to bound its values in the narrower range $r_{\mathrm{break}}\in(10^{5},\,10^{6})r_{\mathrm{g}}$. This is an important result, especially when taken together with our finding that the jet shape transition may be a common phenomenon in nearby or even most of the AGN.

We note the following. The black hole mass of 1H 0323+342 is suspected to be underestimated (León Tavares et al., 2014; Hada et al., 2018). If we use for this source the mass $M=10^{8.6}M_{\odot}$, obtained using the relation between black hole mass and bulge luminosity (León Tavares et al., 2014), 1H 0323+342 yields $r_{\mathrm{break}}=5.6\times 10^{6}r_{\mathrm{g}}$, falling much closer to the discussed above range of $r_{\mathrm{break}}/r_{\mathrm{g}}$ distances. This may provide an additional argument favoring a higher black hole mass for this source.

We have compared our results for the radio galaxy NGC 6251 with those obtained earlier for this source by Tseng et al. (2016). We have found that the jet shape transition region in this source is at $(1.6\pm 0.2)\times 10^{5}\,r_{\mathrm{g}}$, assuming viewing angle of $18\degr$ and black hole mass of $6\times 10^{8}\,M_{\odot}$ (see LABEL:t:sourprop, Table 4). This is slightly smaller compared to $(1-2)\times 10^{5}$ Schwarzschild radius estimated by Tseng et al. (2016), who assumed the same black hole mass and a viewing angle of $19\degr$. The small difference might be caused by different techniques used to derive it. First, we measured transverse jet widths from the stacked image of the source, using 14 epochs at 15 GHz from the MOJAVE program and archival VLBA data. Second, we have taken into account the synchrotron opacity of the jet base by introducing the parameter $r_{0}$ that reflects an offset of the apparent 15 GHz core from the true jet apex.

Of 12 sources with observed change in a jet boundary shape 6 are FR I type, 2 are FR II type, and 4 have uncertain classification based on published radio images. This may mean that different environments expected in these two different types of sources on large scales are either the same on the smaller scales, or affect the jet shape in the same way up to $10^{6}r_{\mathrm{g}}$.

3 Modeling relativistic jet with a shape break

3.1 Qualitative consideration

Both analytical (see below) and phenomenological (Potter & Cotter, 2013, 2015) considerations as well as numerical simulations (Komissarov et al., 2009; Tchekhovskoy et al., 2009; Porth et al., 2011) show that for moderate initial magnetization of a jet $\sigma_{\rm M}\sim 10$–$10^{2}$, where

[TABLE]

is the Michel magnetization parameter, the flow transits from a magnetically dominated regime at small distances $r$ from the origin to a particle dominated regime at larger distances. Here $\Psi_{0}$ and $\Omega_{0}$ are the total magnetic flux and characteristic angular velocity of the “central engine” respectively. Accordingly, $\mu=m_{\rm p}c^{2}+m_{\rm p}w$ is the relativistic enthalpy, where $w$ is the non relativistic enthalpy, and $m_{\rm p}$ is a particle mass. Here we assume a leptonic jet, so $m_{\rm p}$ is the electron mass. Below for simplicity we consider not so large temperatures, so that $w\ll c^{2}$. Finally, $\eta$ is the particle-to-magnetic flux ratio.

Indeed, the physical meaning of the Michel magnetization parameter is the maximum Lorentz factor $\gamma$ of the hydrodynamical flow when all the electromagnetic energy flux is transferred to particles. On the other hand, for quasi-cylindrical jets the following asymptotic solution for magnetically dominated flow exists (see e.g., Beskin, 2009)

[TABLE]

where $R_{\rm L}=c/\Omega_{0}$ is the light cylinder radius, and $r_{\perp}$ is the distance from the jet axis. For the black hole spin $a_{*}=0.5$, $R_{\rm L}\approx 14.9\,r_{\rm g}\approx 2.2\times 10^{15}\;(M_{\rm BH}/10^{9}\,M_{\odot})\;{\rm cm}\approx 7.1\times 10^{-4}\;(M_{\rm BH}/10^{9}\,M_{\odot})\;{\rm pc}$. Here and below we use the maximum BH energy extraction rate condition $\Omega_{\rm F}=\Omega_{\rm H}/2$ (Blandford & Znajek, 1977). For observed pc scale jets, the jet width $d$ at the jet shape break point reaches 1 pc. This means that at the transition point $d/2R_{\rm L}>\sigma_{\rm M}$, and the flow cannot be still magnetically dominated. As was shown by Nokhrina et al. (2015) who have analysed about 100 AGN jets, $\sigma_{\rm M}\sim 10-50$ is a reasonable value constrained by the observations. The observed median value of 1.02 for the $k$-index also clearly points to a ballistic plasma motion. This suggests that the jet is dominated by the plasma bulk motion kinetic energy at the deprojected distance longer than $\sim 100$ pc or $\sim 10^{7}r_{\mathrm{g}}$ rather than by the Poynting flux, as expected close to the launching region.

For this reason we aim to explain the break in the $d(r)$ dependence as a consequence of a transition from the magnetically dominated to the particle dominated regime. Below we present the main results of our semi-analytical consideration. Our goal is in evaluating the dependence of the jet width $d$ on an ambient pressure profile $P_{\rm ext}(r)$. The results for the cold jet are presented in Beskin et al. (2017), while here we consider the semi-analytical results for a warm outflow.

3.2 Semi-analytical model

Basic equations describing the internal structure of relativistic and non relativistic jets within the Grad-Shafranov approach are now well-established (Heyvaerts & Norman, 1989; Pelletier & Pudritz, 1992; Lery et al., 1998; Beskin & Malyshkin, 2000; Beskin, 2009; Lyubarsky, 2009). This approach allows us to formulate the problem of finding a stationary axisymmetric magnetohydrodynamic outflow structure (a jet solution) using a set of two differential equations on a magnetic flux function $\Psi$ and an Alfvénic Mach number ${\cal M}$. These equations are Bernoulli equation and Grad–Shafranov equation of a force balance perpendicular to magnetic surfaces. The approach allows us to determine the internal structure of axisymmetric stationary jets knowing in general case five “integrals of motion”, i.e., energy $E(\Psi)$ and angular momentum $L(\Psi)$ flux, electric potential which connects with angular velocity $\Omega_{\rm F}(\Psi)$, entropy $s(\Psi)$, and the particle-to-magnetic flux ratio $\eta(\Psi)$. All these values are to be constant along magnetic surfaces $\Psi=\mathrm{const}$. Once the Grad–Shafranov and Bernoulli equations are solved for the given integrals, all the other flow properties, such as particle number density, four-velocity, electric current, and Lorentz factor, can be determined from algebraic equations (e.g., Beskin, 2009). In particular, it was shown that a jet with total zero electric current can exist only in the presence of an external medium with non-negligible pressure $P_{\rm ext}$. Thus, it is the ambient pressure $P_{\rm ext}$ that is expected to determine the transverse dimension of astrophysical jets. In general, it is a complicated problem to solve the set of Bernoulli and Grad–Shafranov equations. An additional complication is connected with the change of a system type from elliptical to hyperbolic. So, to tackle the problem different simplifications are introduced. Here we simplify the problem, assuming the flow is highly collimated and can be described within the cylindrical geometry, in which case it can be solved numerically (Beskin & Malyshkin, 2000).

On the other hand, careful matching of a solution inside the jet with the external medium has not been achieved up to now. The difficulty arises with having a very low energy density of the external medium in comparison with the energy density inside the relativistic jet. For this reason, in most cases an infinitely thin current sheet was introduced. Moreover, an ambient pressure was often modelled by homogeneous magnetic field $B_{\rm ext}^{2}/8\pi=P_{\rm ext}$.

Below we use the approach developed by Beskin et al. (2017). This paper is later referred to as B17. We propose a flow with an electric current closing fully inside a jet. This is achieved by a natural assumption that the integrals $L$ and $\Omega_{\rm F}$ vanish at the jet boundary. The second assumption of the model is a vanishing flow velocity at the jet boundary, which leads to vanishing of a poloidal magnetic field component along with a toroidal due to current closure. As a consequence, only a thermal pressure, defined by a sound velocity $c_{\rm jet}$ and particle number density $n_{\rm jet}$, is left at the jet boundary to balance the ambient medium pressure without a current sheet. We solve Grad–Shafranov and Bernoulli equations for the flux function $\Psi(r_{\perp})$ and the square of an Alfvénic Mach number ${\cal M}^{2}(r_{\perp})$. The local non-relativistic enthalpy $w$ for a polytropic equation of state with politropic index $\Gamma=5/3$ can be written as

[TABLE]

where the local particle number density $n$ is obtained from the equation

[TABLE]

We solve the system of MHD equations (B17) for the boundary conditions $\Psi(0)=0$ and

[TABLE]

We should note that due to vanishing of the integrals $L(\Psi)$ and $\Omega_{\rm F}(\Psi)$ at the jet boundary, the thickness of the final current closure domain tends to zero and in B17 it was is not resolved. However, as it was shown, that the total pressure in this region is strictly conserved:

[TABLE]

This means that the solution we obtain up to the boundary does contain the residual current and, thus, the toroidal magnetic field $B_{\varphi}$.

The main difference between the result presented here and the result by B17 is in more accurate account for the thermal terms, which can be seen in Equation 4. To obtain the solution, we employ the following iterative procedure. For each fixed fast magnetosonic Mach number at the axis ${\cal{M}}_{0}^{2}$ we initially set $P_{\rm ext}$ at the jet boundary. It defines the particle number density at the boundary $n_{\rm jet}$, and together with ${\cal{M}}_{0}^{2}$ — the particle number density at the axis $n_{0}$. Having set the latter, we solve MHD equations across a jet from the axis outwards and calculate the jet pressure at the boundary provided by the solution $P^{\mathrm{(solution)}}$. By iterations we find self-consistently such $P_{\mathrm{ext}}$ that is equal to one, provided by the solution: $P^{\mathrm{(solution)}}=P_{\mathrm{ext}}$. Thus, we obtain the dependence of a jet pressure at the boundary as a function of a local jet width $d$.

This procedure fully determines the solution of our problem. For each magnitude of the external pressure the obtained solution is a crosscut at $r={\rm const}$. Piling of these different crosscuts is a solution for an outflow in which one may neglect by the derivatives over $r$ in comparison with the derivatives over $r_{\perp}$ in the two-dimensional Grad–Shafranov and Bernoulli equations. This can be done for highly collimated, at least as a parabola, outflows (Nokhrina et al., 2015) and flows with small opening angles (Tchekhovskoy et al., 2009).

We find that for the chosen sound velocity at the boundary $c_{0}^{2}=0.001c^{2}$ the thermal effects may be neglected in the outflow volume, playing an important role only at the outflow boundary. It turns out that the resultant dependence of pressure at the jet boundary as a function of jet radius obtained by B17 and here start to differ somewhat only for large ${\cal{M}}_{0}^{2}$ (this value is of an order of 10, but depends on the initial magnetization), affecting the flow boundary shape downstream of the equipartition transition, and the effect on $k_{2}$-index is of the order of a few per cent. We will address the particular effects of higher temperature in the future work.

The proposed jet model with an electric current enclosed inside the jet has a natural sheath structure, observed, for example, in the M87 jet (Mertens et al., 2016). Due to choice of integrals, the outer parts of a jet have slower velocities, tending to non-relativistic with $\gamma(d/2)=1$. Such a sheath may be produced by different mechanisms: it may be a slower disk wind or an outer jet disturbed and slowed down by the pinch instability (Chatterjee et al., 2019). In our model it appears naturally as a consequence of a jet transiting into the ambient medium with the hydrodynamical discontinuity only (B17).

3.3 Transition from magnetically dominated to particle dominated flow

It is necessary to stress that this system of equations can describe both magnetically and particle dominated flow, with the physical answer (including the jet boundary radius $d/2$) depending on one external parameter only, namely, on the ambient pressure $P_{\rm ext}$. In Figure 6 we show the dependence of the dimensionless ambient pressure $\tilde{p}$ on a dimensionless jet width $\tilde{d}$ obtained by solving numerically the system of Grad–Shafranov and Bernoulli equations B17. The pressure is plotted in units of

[TABLE]

so that $P_{\rm ext}=\tilde{p}\,p_{0}$, and the jet width in units of light cylinder radius is $d=\tilde{d}\,R_{\rm L}$. We observe (see Figure 6) that the pressure has a different power law dependence on the jet radius for small and large $d$. For each magnetization $\sigma_{\rm M}$, this behavior holds, with the change between two profiles occurring at different jet widths. For $\sigma_{\rm M}=50$ the pressure changes its dependence on $d$ from

[TABLE]

closer to the jet base to

[TABLE]

further downstream. The particular exponents of the power laws depend weakly on $\sigma_{\rm M}$.

We assume the equilibrium between jet and ambient medium pressure. In order to model a jet shape break position along the jet, we need to introduce the exerted pressure dependence on $r$, which we choose in the power law form

[TABLE]

Such a pressure profile is consistent with Bondi flow (Quataert & Narayan, 2000; Shcherbakov, 2008; Narayan & Fabian, 2011) having $b\in(1.5;\;2.5)$ for different models, with the limiting value 2.5 for classical supersonic Bondi flow. This power law with $b\approx 2.0$ allows us to reproduce well both the parabolic jet form upstream the break and conical downstream. Using power laws Equation 8, Equation 9, and Equation 10, we obtain for small distances $r$ (magnetically dominated regime)

[TABLE]

Accordingly, for large distances (saturation regime)

[TABLE]

As we see, qualitatively, the power indices are in good agreement with the observational data. Thus, we are able to reproduce the jet boundary shape behaviour without introducing two different pressure profiles, as was done in (Asada & Nakamura, 2012). Having the reasonable pressure dependence on a distance, we reproduce both power laws in a jet shape. For example, for a central mass $M=10^{9}\,M_{\odot}$ and black hole spin $a_{*}=0.5$ the light cylinder radius is $R_{\rm L}\approx 7\times 10^{-4}$ pc. We also set the total magnetic flux in an outflow $\Psi_{0}=10^{32}\;{\rm G\,cm^{2}}$ (Zamaninasab et al., 2014; Zdziarski et al., 2015; Nokhrina, 2017), which gives the value $B(r_{\rm g})\approx 1400$ G. Thus, for these test parameters the jet width at the break, designated by a star in Figure 6, has typical values $0.2-1.0$ pc in agreement with the observational results in Table 4.

In dimensionless units the point of transition from one power law for pressure as a function of a jet width to the other is defined by one parameter only: the jet initial magnetization. In the equipartition regime the jet bulk Lorentz factor is $\gamma=\sigma_{\rm M}/2$. The observed kinematics in parsec-scale jets constrains the initial magnetization to a value $\lesssim 100$ (Lister et al., 2016), while estimates for $\sigma_{\rm M}$ based on core-shift effect measurements provide the preferred value $\lesssim 20$ (Nokhrina et al., 2015). In dimensional units the jet width at the break depends also on BH mass and spin. The distance to a shape transition along the jet is determined by the total magnetic flux in a jet and the ambient medium pressure. We address the question of bounding these parameters in the next paper (in preparation).

3.4 Magnetization

In this subsection we check whether the break in a jet shape corresponds to the transition from the magnetically-dominated into the equipartition regime. The jet magnetization is defined as the ratio of Poynting flux

[TABLE]

to particle kinetic energy flux

[TABLE]

where $n$ is particle number density in the jet proper frame. Using the standard expressions for ideal MHD velocities and electric and magnetic fields, one obtains the following expression for the magnetization:

[TABLE]

Using the definitions of bulk Lorentz factor $\gamma$ and total current $I$, we rewrite it as

[TABLE]

In order to check $\sigma$ along the jet, we calculate the maximal magnetization across the jet for each given distance $r$. The magnetization is always much less than the unity at the jet axis and at the jet boundary. The first holds everywhere, since the Poynting flux behaves at the jet axis as

[TABLE]

if the current density $j$ has no singular behavior at $r_{\perp}=0$. Thus, $\sigma\rightarrow 0$ at the axis. The same holds for the boundary in a case of the full electric current closure. Due to specific choice of integrals $E(\Psi)$, $L(\Psi)$, and $\Omega_{\rm F}(\Psi)$ (B17), the Poynting flux together with the magnetization reach their maximum values at $\Psi=\Psi_{0}/2$. It is at this magnetic field line the flow attains its highest Lorentz factor across the jet for the given distance from the central source. Thus, we choose the maximal magnetization reaching approximately unity as a criteria of a flow attaining the ideal MHD equipartition regime.

In Figure 7 we present the maximal magnetization and the break in a jet form. We plot the modelled jet boundary shape for $\sigma_{\rm M}=50$, BH and jet parameters the same as in subsection 3.3. The position of a jet shape break along a jet depends on an ambient pressure profile (Equation 7 and Equation 10), and we use here, as an example, $P_{0}=10^{-6}\;{\rm dyn/cm^{2}}$ at $r_{0}=10$ pc. We see that the break in jet shape occurs roughly at the distance from the BH, where the flow magnetization becomes equal to unity. For the higher initial magnetization it takes the larger transverse jet dimension in $R_{\rm L}$ to accelerate the flow up to equipartition, according to Equation 2.

4 Discussion

4.1 Role of a Bondi sphere

In this paper we propose that the jet form change, observed in a dozen of nearby sources, may be explained by an internal flow transition from magnetically dominated to particle dominated regime with the smooth external pressure profile $P\propto r^{-2}$. There are indications, however, that the ambient pressure may have different profiles at different scales. The measurements of particle number density in ISM by Russell et al. (2015) suggest $\rho\propto r^{-1}$ from about 400 pc down to expected Bondi radius $r_{\rm B}\sim 100-250$ pc. The temperature profile on scales $100-1000$ pc is roughly constant. This means that just outside, or even inside, the Bondi radius, pressure profile is $P\propto r^{-1}$, with no information on it inside a sphere $\sim 150$ pc. The position of a sphere of influence is expected to be at a distance $10^{5}-10^{6}\;{r_{\mathrm{g}}}$ (Blandford et al., 2019). The position of a transition point $r_{\mathrm{break}}$ from magnetically dominated to particle dominated regimes predicted by our model for reasonable parameters lay in general in the same interval or inside $r_{\mathrm{B}}$. For example, in the case of M87 we observe $r_{\mathrm{break}}\approx 40$ pc (Nokhrina et al., 2019) smaller than $r_{\mathrm{B}}$. The same phenomenon has been noted by Nakahara et al. (2018) for NGC 4261, where the structural transition lies well inside the expected sphere of influence. In Figure 8 we present a cartoon for a jet shape with different ambient pressure profile. Inside the Bondi sphere the jet is accelerating effectively up to the distance $r_{\mathrm{break}}$, with predicted parabolic boundary shape described by Equation 11. This is the acceleration and collimation zone (ACZ) discussed by Blandford et al. (2019). For $r_{\mathrm{break}}<r<r_{\mathrm{B}}$ the jet assumes a close to conical form Equation 12. Up to $r_{\mathrm{B}}$ the jet stays in equilibrium with the ambient pressure $P_{\mathrm{jet}}=P_{\mathrm{ext}}$. If for $r>r_{\mathrm{B}}$ the ambient pressure has a more shallow profile, the conical particle-dominated jet may become overpressured with a possible appearance of a standing shock. Thus, we predict the presence of a standing bright feature, associated with a shock, outside the Bondi sphere and downstream the break in jet shape. At this shock we may expect plasma heating, with the flow continuing a conical expansion (Blandford et al., 2019). The position of HST-1 in M87 jet in a close vicinity of expected $r_{\mathrm{B}}$ and downstream the $r_{\mathrm{break}}$ supports this picture.

4.2 Additional observational evidence of the break point and predicted evolution of plasma acceleration

For each of the 10 sources with a jet geometry transition detected (Table 4), we checked for slow pattern ($\beta_{\textrm{app}}<0.2c)$ jet features in Lister et al. (2019). We examined if their median locations with respect to the core are positionally associated, i.e., they match within the errors with the position of the derived jet shape break. We found that five sources have a quasi-stationary feature in the region where jet changes its shape, as expected (see discussion in subsection 4.1). This is a factor of 1.5 larger compared to a ratio from the overall statistics of jet kinematics analysis performed at 15 GHz, which reveals a fraction of quasi-stationary jet features to be about 30% (Lister et al., 2019), applying the criterion $\beta_{\textrm{app}}<0.2c$. We underline that the MOJAVE kinematic analysis uses conservative criteria in cross-identifying components between epochs and selecting robust ones (Lister et al., 2019). This means that the 50% fraction of sources which show a standing feature in the break point region should be considered as a lower limit. This analysis is also conservative because of the requirement of the feature to be coincident with the detected break point. As discussed above, the shock may be located downstream the jet in the vicinity of $r_{\rm B}$, which position is usually not known. We note that two sources included from other studies, 1H 0323+342 and M87, have the jet shape transition at distances larger than maximum angular scales probed by the MOJAVE 15 GHz observations.

We plot in Figure 9 the maximum Lorentz factor of a bulk plasma motion along a jet, which we obtain within our semi-analytical model. The predicted pattern of a bulk Lorentz factor acceleration in magnetically dominated domain is $\gamma\propto r_{\perp}$, which provides for a parabolic jet $\gamma\propto r^{0.5}$. After the flow reaches equipartition, the acceleration continues slower than any power-law (logarithmically slow) (e.g., Beskin & Nokhrina, 2006). There is also a transitional zone between the two regimes. Thus, we would expect for the sources with the detected jet shape break and superluminal motion the following kinematics pattern: efficient Lorentz factor growth before the break point, and cessation of it in the conical region. This expected Lorentz factor behaviour was reported by Hada et al. (2018). The observed in radio band velocity map in M87 (Mertens et al., 2016) shows the acceleration saturation much earlier than the jet shape break. However, observations in the optical band (Biretta et al., 1999) support the acceleration of plasma continuing further, with reported $\gamma=6$ at HST-1, situated downstream the jet shape break. This may point to non detection of fast components in radio.

This prediction is consistent with observations by the MOJAVE program that acceleration is a common property of jet features (e.g., Homan et al., 2015; Lister et al., 2019), reflecting a tendency for increasing Lorentz factors near the base of the jet, with decreasing or constant speeds being more common at projected distances $\gtrsim 10-20$ parsecs (Homan et al., 2015). While decreasing speeds are not a prediction of this model for a change in jet shape, they could naturally occur if the reduction in positive acceleration is also accompanied by entrainment of external material into the jet.

5 Summary

Pushkarev et al. (2017) studied AGN jet shapes by measuring the power low index $k$ assuming a $d\propto r^{k}$ dependence of the observed deconvolved jet width $d$ on the apparent distance from its core $r$. Most of the jets exhibited $k$ values in the range from 0.5 to 1.5. As it was clearly demonstrated by Pushkarev et al. (2017), high-quality, high-dynamic-range stacked images are needed for an analysis of this kind in order to trace the full jet channel. In view of a few recent exciting reports on jet shape transitions from parabolic to conical (e.g., Asada & Nakamura, 2012; Giroletti et al., 2008; Tseng et al., 2016; Hervet et al., 2017; Hada et al., 2018; Akiyama et al., 2018; Nakahara et al., 2018, 2019), we have performed a systematic search of such transition using MOJAVE 15 GHz stacked images, supplementing some of them with available single epoch 1.4 GHz VLBA images to trace larger scales.

Using an automated analysis approach, we have found 10 jets with such transition out of 367 analyzed: 0111$+$021, 0238$-$084, 0415$+$379, 0430$+$052, 0815$-$094, 1133$+$704, 1514$+$004, 1637$+$826, 1807$+$698, 2200$+$420. Their redshifts lie in the range $z<0.07$ except for 0815$-$094, whose redshift is unknown. For the full analyzed sample the redshift values cover the range from 0.004 to 3.6 with the typical value being about 1. This low-$z$ coincidence is unlikely to have occurred by chance. Taken together with an analysis of possible biases, we conclude that a genuine effect is present in the data for which VLBA reaches the linear resolution better than 1 pc. We would also predict that the BL Lac object 0815$-$094 is a nearby AGN.

This finding leads to the following important conclusion. A transition from parabolic to conical shape may be a general property of AGN jets. At the same time, we note that AGN observed at higher redshifts typically have higher luminosities and kinetic power, which can affect the collimation properties. This conclusion has important implications for jet models, astrophysics and astrometry of AGN. Measuring this phenomenon requires a search within nearby AGN which is the subject of our current followup study, or increasing the resolution by using Space VLBI (e.g., Giovannini et al., 2018) or high dynamic range high frequency VLBI imaging.

The deprojected distance $r_{\mathrm{break}}$ from the nucleus to the break zone is found to be typically 10 pc. Even more interesting due to its relation to jet formation and acceleration models is this value measured in gravitational radius units. We find the range to be $r_{\mathrm{break}}\in(10^{5},\,10^{6})r_{\mathrm{g}}$ which corresponds to the typical Bondi radius.

We have developed the following model to explain the observed jet shape break. The accurate matching of a jet outflow with an ambient medium B17 predicts a change in jet shape from parabolic to conical if the ambient medium pressure is assumed to be governed by Bondi accretion. Within the model, a smaller external pressure is needed to support a jet than in earlier models. The transition of predicted jet shape from parabolic to conical occurs in the domain where the bulk plasma kinetic energy flux becomes equal to the Poynting energy flux, i.e., where the bulk flow acceleration reaches saturation (Beskin & Nokhrina, 2006). From studying the break properties we can estimate black hole spin and/or mass, jet total magnetic flux, and ambient medium properties as discussed by Nokhrina et al. (in prep.).

The following two model predictions are supported observationally. The break point, where jets start to be plasma dominated energetically, might be a preferable domain for shocks. We detect standing jet features in this region from MOJAVE analysis (Lister et al., 2019) in at least a half of the AGN targets. The plasma acceleration is predicted to decrease significantly at the transition region, which is consistent with MOJAVE acceleration results (Homan et al., 2015; Lister et al., 2019).

Our finding also implies the following (see also discussion in Algaba et al., 2017). The well-known effect of the apparent shift of the core position with frequency due to synchrotron self-absorption does not follow the $r_{\mathrm{core}}\propto\nu^{-1}$ law all the way up to the true jet base, since a $-1$ power low index is expected only for a conical jet (Blandford & Königl, 1979; Lobanov, 1998). Geometrical and physical estimates made on the basis of core shift measurements will need to take this into account while VLBI and VLBI-Gaia astrometry applications will need to correct for it (Porcas, 2009) in cases where very high accuracy is required.

Acknowledgments

We thank Dan Homan, Eduardo Ros, Luigi Foschini, Dave Meier as well as the anonymous referee for valuable comments and suggestions. This study has been supported in part by the Russian Science Foundation: project 16-12-10481. TS was supported by the Academy of Finland projects 274477, 284495, and 31249. This research has made use of data from the MOJAVE database that is maintained by the MOJAVE team (Lister et al., 2018). The MOJAVE program is supported under NASA-Fermi grant NNX15AU76G. This study makes use of 43 GHz VLBA data from the VLBA-BU Blazar Monitoring Program, funded by NASA through the Fermi Guest Investigator Program. The VLBA is an instrument of the National Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated by Associated Universities, Inc. This research has made use of NASA’s Astrophysics Data System.