MNRAS 000

, 1-17 (2024)

## A comparative study of radio signatures from winds and jets: Modelling synchrotron emission and polarization

Moun Meenakshi 1 , ★ Dipanjan Mukherjee 1 , † , Gianluigi Bodo 2 , Paola Rossi 2 and Chris M. Harrison 3

1 Inter-University Centre for Astronomy and Astrophysics, Pune- 411007, India

2 INAF, Osservatorio Astrofisico di Torino, Strada Osservatorio 20, I-10025 Pino Torinese, Italy

3 School of Mathematics, Statistics and Physics, Newcastle University, Newcastle upon Tyne NE1 7RU, UK

Accepted XXX. Received YYY; in original form ZZZ

## ABSTRACT

Outflows driven by active galactic nuclei (AGN) are seen in numerous compact sources; however, it has remained unclear how to distinguish between the driving mechanisms, such as winds and jets. Therefore, our study aims to offer observational insights from simulations to aid in this distinction. Specifically, in this paper, we investigate the evolution of wide-angled, moderately relativistic, magnetized winds and analyze their non-thermal radio emission and polarization properties. We find that the evolution of winds varies depending on factors such as power, density, and opening angle, which in turn influence their observable characteristics. Additionally, different viewing angles can lead to varying observations. Furthermore, we note distinctions in the evolution of winds compared to jets, resulting in disparities in their observable features. Jets typically exhibit a thin spine and hotspot(s). Winds manifest broader spines or an 'hourglass-shaped' bright emission in the cocoon, which are capped by bright arcs. Both display high polarization coinciding with the bright spine and hotspots/arcs, although these regions are relatively compact and localized in jets when compared to winds. We emphasize the importance of high resolution, as we demonstrate that emission features from both jets and winds can become indistinguishable at lower resolutions. The distribution of polarization is largely unaffected by resolution, though lower polarization becomes more noticeable when the resolution is decreased.

Key words: galaxies: active - ISM: jets and outflows - MHD - turbulence - methods: numerical

## 1 INTRODUCTION

Outflows driven by the Active Galactic Nuclei (AGN) have been observed through multi-wavelength studies in a plethora of sources (Tombesi et al. 2015; Nesvadba et al. 2021; Zovaro et al. 2019; Vayner et al. 2021; Girdhar et al. 2022; Shen et al. 2023). These outflows can arise due to either collimated jets or wide-angled winds resulting from the activity of the central supermassive black hole (Harrison &amp; Ramos Almeida 2024). Theoretical models and observations indicate that both jets and winds are launched by distinct mechanisms. The Blandford-Znajek mechanism is attributed to the generation of jets from the CSMBH (Blandford &amp; Znajek 1977), rendering them collimated and highly magnetized. The formation of disc winds is explained by three popular models, which are radiation driving (Murray et al. 1995; Proga et al. 2000), magnetic driving (Blandford &amp; Payne 1982; Yuan et al. 2015; Wang et al. 2022) and thermal driving (Begelman et al. 1983). These winds are wide-angled outflows, as presented in the models from Faucher-Giguère &amp; Quataert (2012) and Nims et al. (2015).

is detected as radio emission (Blandford &amp; Ostriker 1978; Melrose &amp; Pope 1993; Mukherjee et al. 2020, 2021). A recent survey by Sabater et al. (2019) has shown that almost all high mass galaxies ( 𝑀 &gt; 10 11 𝑀 ⊙ ) show radio emission beyond a threshold limit of 10 21 W Hz -1 . In radio-loud sources 1 , the jets are expected to dominate the radio emission, whereas the origin in radio-quiet AGN is widely debated (Panessa et al. 2019; Kharb &amp; Silpa 2023). In such sources, it can depend on various mechanisms, including synchrotron emission from electrons accelerated in low-power magnetized jets (Falcke et al. 2000; Gallimore et al. 2006), wide-angle winds (Zakamska et al. 2016; Peng et al. 2020), and star formation (de Vries et al. 2007; Condon et al. 2013).

The primary source of emission from both winds and jets is attributed to shocks that accelerate non-thermal electrons. These accelerated electrons generate incoherent synchrotron radiation, which

- ★ E-mail:mounmeenakshi@iucaa.in (MM)
- † E-mail:dipanjan@iucaa.in (DM)

Polarimetry studies (Sebastian et al. 2020) have found distinctions in polarization features that enable differentiation between radio emissions arising from star formation and AGN-driven outflows. Several recent studies indicate that traditional definitions of radioquiet may not always mean a lack of AGN winds or jets (see e.g. Jarvis et al. 2019, 2021; Silpa et al. 2020, 2023); nonetheless, in several instances, the nature of the emission remains ambiguous. A potential solution lies in the radio emission and polarization study

1 'radio-loud' sources are characterized by high radio luminosity (e.g. 𝑃 1 4 GHz . &gt;&gt; 10 24 WHz 1 ) and high ratio-to-optical flux density ratio (&gt;10), -and vice-versa for radio-quiet sources (Padovani 2017).

of resolved jets/winds, as it can aid in distinguishing between various radio emission mechanisms. The expectation stems from the understanding that jets and winds are launched through different mechanisms and are therefore expected to interact differently with their surroundings. Consequently, when these two outflows encounter a similar ambient medium, observable differences in their emission and polarization characteristics may arise, although this area remains largely unexplored. This serves as the motivation for our study (also discussed in Meenakshi et al. (2023), hereafter referred to as Paper I) of jets and winds spanning a few kpc in scale, where we aim to analyze their observable properties.

In Paper I, we focused on magnetized jets interacting with turbulent magnetic fields. The effect of magnetized Shocked Ambient Medium (SAM) on the net polarization from the jet was also explored. In this subsequent paper, we now focus on the above features from winds launched in a magnetized ambient medium. In our studies, the jets and winds are launched in a hot halo at temperatures of 10 7 K equipped with a turbulent magnetic field. Thus, our primary focus is on examining the evolution of jets/winds in an idealized scenario and understanding how this impacts their radio emission characteristics over time. This paper is organized as follows. The simulation setup for the study is briefly discussed in Sec. 2. Detailed results are presented in Sec. 3, where we specifically examine the dynamics, emission, and polarization characteristics in separate subsections. We discuss the distinctions in the emission and polarization features of jets from Paper I and winds from this study in Sec. 4. In Sec. 5, we discuss our findings and summarize our inferences regarding the jets and winds derived from our simulations.

## 2 SIMULATION SETUP

The setup for the ambient medium and the wind with mildly relativistic speeds in this study is kept similar to Paper I for the jets. Thus, initially, a hot halo is set up which is in the hydrostatic equilibrium in the presence of an external gravitational potential ( 𝜙 𝑟 ( ) ) as,

$$\frac { d p _ { a } } { d r } = - \rho _ { a } ( r ) \frac { d \phi ( r ) } { d r } ; \quad p _ { a } ( r ) = \frac { \rho _ { a } ( r ) } { \mu m _ { a } } k _ { B } T _ { h } \quad \quad ( 1 ) \quad \text{ in the 1} \\ \text{$\text{ Braito}$}.$$

from which we derive

$$p _ { a } ( r ) = ( n _ { 0 } k _ { B } T _ { h } ) \exp \left [ - \left ( \frac { \mu m _ { a } } { k _ { B } T _ { h } } \right ) \phi ( r ) \right ], T _ { h } = 1 0 ^ { 7 } \, K \quad \quad ( 2 ) \quad \text{analysis}$$

where 𝑝 𝑎 and 𝜌 𝑎 = 𝜇𝑚 𝑛 𝑎 ℎ are respectively the pressure and density of the halo gas at spherical radius 𝑟 , 𝜇 = 0 6 is the mean molecular . weight of the fully ionized gas, 𝑚 𝑎 is the atomic weight, and 𝑛 0 is the number density at 𝑟 = 0. The Boltzmann constant is shown using 𝑘 𝐵 , and 𝑇 ℎ is the temperature of the halo. The total gravitational potential ( 𝜙 𝑟 ( ) ) is the sum of the baryonic and dark-matter potential, as described in Paper I in detail. Then, the turbulent magnetic fields are introduced into the ambient medium, and the winds are launched from the lower 𝑍 side of the simulation domain.

However, there are several differences, which are explained below. The simulation domain is a large box with dimensions of 8 kpc × 8 kpc × 4 kpc, corresponding to a grid of 𝑋 × 𝑌 × 𝑍 , with a resolution of 800 × 800 × 400 cells. As the winds are expected to be broad, the wide box encompasses the cocoon and SAM for the winds by the time its head reaches the 𝑍 -altitude of 4 kpc. Similar to Paper I, we employ the RMHD module in pluto (Mignone et al. 2007) along with a five-wave HLLD Riemann solver (Mignone et al. 2009). We use a piece-wise parabolic scheme for reconstruction (Mignone et al. 2005) and a third-order Runge-Kutta method for time integration. In order to ensure numerical instability we use a more diffusive (HLL)

Riemann solver and MIN-MOD limiter for shocked cells which are identified using a pressure condition with 𝛿𝑝 / 𝑝 𝑚𝑖𝑛 &gt; 6 in regions near the injection zone ( | 𝑋,𝑌 | ≤ 0 5 and . 𝑍 ≤ 0 5 kpc). .

The turbulent magnetic field in the halo is incorporated by interpolating from a domain of size 4 kpc 3 , as detailed in Paper I for compact jets. Given the larger simulation box here, the field is extended along the 𝑋 and 𝑌 axes using periodic boundary conditions, as also done in Paper I for large-scale jet run. In this study, the maximum correlation length of the external magnetic field is 1000 pc. The investigation for varying the correlation lengths has already been performed in Paper I, and we anticipate similar outcomes in this regard. Consequently, our primary focus here is on examining the dynamics of winds with different powers launched in a similar turbulent magnetic field configuration. We use the constrained transport approach (Balsara &amp; Spicer 1999; Mignone &amp; Del Zanna 2021) to maintain the divergence-free nature of the magnetic field throughout the simulation. Additional specifics regarding the setup and numerical schemes can be found in Paper I, to which we direct the reader for further details.

The simulations performed in this study are listed in Table 1. The nomenclature of the simulations indicates the power of the injected wind and the correlation length of the ambient magnetic field, as explained in the footnote of Table 1. The half-cone opening angle for the wind ( 𝜃 𝑊 ) is set at 70 ◦ for a broader injection than the jets, which were launched with an angle of 7 ◦ . The radius of the injection zone is 0.2 kpc here; thus the diameter for the injected wind is resolved with 40 cells for a resolution of 10 pc in our runs. Simulations with highdensity contrast ( 𝜂 𝑤 ) between the wind and the ambient medium (W43\_L1000 D HD ) are also performed to explore the effect of a ( , ) dense wind on different physical measures. We also conduct runs for narrow winds, which are launched with a half-cone opening angle ( 𝜃 𝑤 ) of 35 ◦ , and explore this for both light and dense wind cases (W43\_L1000 N HDN ). It can be seen that the winds exhibit mildly ( , ) relativistic velocities in comparison to the jets in Paper I, which were launched with speeds reaching up to 0.99c. Outflows with velocities of around 0.25c, driven off by the accretion disc winds at distances of a few gravitational radii from the black hole, have been observed in the nuclear regions of several quasars (see e.g. Pounds et al. 2003; Braito et al. 2007; Reeves et al. 2009; Tombesi et al. 2015). These winds become mildly relativistic with distance and are likely to attain velocities similar to the values chosen in our study. For post-process analysis, we consider an electron population whose energy spectrum follows a power law with an index of 2.2, similar to the one assumed in Paper I. The synchrotron emissivities from radiation emitted by the acceleration of electrons in the presence of a magnetic field are estimated using the expressions given in Sec. 3 of Paper I. In the next section, we discuss the results of the dynamics and observable characteristics of winds from this study.

## 3 RESULTS

## 3.1 Dynamics of wind

In this section, we discuss the structure of the wind, and how its characteristics vary in response to various factors, including wind power ( 𝑃 𝑤 ), density contrast between the injected wind and the ambient medium ( 𝜂 𝑤 ), and the wind's opening angle ( 𝜃 𝑤 ). The general structure of light wind is illustrated through a schematic diagram and a 3D volume-rendered pressure map in Fig. 1 and 2,

Table 1. Parameters for simulations performed in this study.

| Simulation     | 𝜂 𝑤        | 𝑉 𝑤    | 𝑃 𝑤 [ erg s - 1 ]   |   𝐵 0 [ mG ] |   𝑀 𝑤 |
|----------------|------------|--------|---------------------|--------------|-------|
| W42_L1000      | 4 × 10 - 2 | 0.021c | 1 . 01 × 10 42      |        0.033 |  2.68 |
| W43_L1000      | 4 × 10 - 2 | 0.051c | 1 . 03 × 10 43      |        0.069 |  6.36 |
| W44_L1000      | 4 × 10 - 2 | 0.11c  | 1 . 04 × 10 44      |        0.148 | 14.02 |
| W43_L1000D a   | 4 × 10 - 1 | 0.024c | 1 . 04 × 10 43      |        0.1   |  9.5  |
| W43_L1000HD a  | 4          | 0.011c | 1 . 06 × 10 43      |        0.148 | 14.11 |
| W43_L1000N b   | 4 × 10 - 2 | 0.056c | 1 . 02 × 10 43      |        0.076 |  7.1  |
| W43_L1000HDN b | 4          | 0.011c | 0 . 8 × 10 43       |        0.15  | 14.11 |

The simulation label indicates both the wind power and the correlation length of the ambient magnetic field. For example, W42\_L1000 represents a simulation where a wind with a power of approximately 10 42 erg s -1 is launched into an ambient medium with a mean magnetic field of about 10 𝜇 G and a correlation length of 1000 pc.

- a Dense (D, 𝜂 𝑤 = 4 × 10 -1 ), and highly dense (HD, 𝜂 𝑤 = 4) wide-angled winds of power 10 43 erg s -1 .
- b Light and narrow (N, 𝜃 𝑤 = 35 ) and Highly dense and narrow (HDN, ◦ 𝜂 𝑤 = 4 and 𝜃 𝑤 = 35 ) winds of power 10 43 ◦ erg s -1 .
- 𝜂 𝑤 : Ratio of wind density to the ambient gas density at the radius of injection.
- 𝑉 𝑤 : Bulk velocity of the wind.
- 𝑃 𝑤 : Mechanical power of the wind.
- 𝐵 0 : Maximum strength of the toroidal magnetic field of the wind at the injection zone. Note that the magnetization of the wind ( 𝜎 ) as defined in Paper I (Eq. 7) is taken as 0.1 for all the simulations.
- 𝑀 𝑤 : Mach number of the wind.

Figure 1. Schematic diagram showing the structure of wind (left) and pressure map for W44\_L1000 (right) for illustration (similar to Fig. 3).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-fa93947e-30c4-4883-8341-32a7110651d3.jpg)

respectively. Slices depicting density, pressure, and wind tracer 2 in the 𝑌 -𝑍 plane from different simulations are shown in Figs. 3, 4, and 5. These snapshots correspond to the times when the front head of the forward shock reaches close to the top of the simulation domain; indicating that the spatial scale of the winds is similar across all simulations.

- · Structure of the wind: The general structure of the wind is expected to align with the model presented in Weaver et al. (1977) and Faucher-Giguère &amp; Quataert (2012). However, our study reveals
- 2 Tracers are passive scalar quantities between 0 and 1 in pluto, which gets advected with the fluid. The wind tracer (tr1) value at the injection zone is set to 1 at 𝑡 = 0, and zero elsewhere. See http://plutocode.ph.unito.it/ userguide.pdf for further clarification.

deviations in the wind morphology, which primarily depend on the initial density ratio between the wind and the ambient medium ( 𝜂 𝑤 ). Specifically, dense winds ( 𝜂 𝑤 ≳ 4 in this study) are stable and closely adhere to the wide-angle wind model, as depicted in above studies. Conversely, light winds ( 𝜂 𝑤 ≲ 0 4) are more susceptible to external . influences, resulting in instabilities and some deviations from this structure. Therefore, as the density of the wind decreases over time due to expansion, instabilities are likely to emerge. In this section, we focus primarily on the structure of the light wind, detailing the instabilities. A labeled schematic diagram for the structure of the light wind is shown in Fig. 1. Adjacent to it, the pressure map for W44\_L1000 is presented to illustrate the various elements. The main features inside the wind cocoon include a Mach disk illustrating the wind's termination shock, and lateral streams emerging from the

Figure 2. 3D volume rendering of logarithmic pressure and velocity streamlines showing the trajectories of fluid in the cocoon for W43\_L1000.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-b2b91953-f14a-4a55-a7f6-61a8c824543e.jpg)

Figure 3. Logarithmic density (left), pressure (middle), and wind tracer (right) for simulations of wide-angled ( 𝜃 𝑊 = 70 ) light winds ( ◦ 𝜂 𝑤 = 4 × 10 -2 ) of different powers (10 42 -10 44 erg s -1 ) in the 𝑌 -𝑍 plane. The region enclosed inside the cocoon and forward shock (see Appendix A) are represented by white and yellow contours in the right panels.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-0adcaddf-2971-4fb2-9301-7b5b5504d9ed.jpg)

edges of the barrel (or conical-shaped region), which are shown using red arrows. The wind plasma is abruptly stopped at the termination shock of the Mach disk. However, the oblique shock near the outer edges of the barrel deflects the gas moving downstream. This leads to high-velocity lateral streams emerging around the barrel, which can also be seen from the 3D volume rendering of the pressure along with fluid velocity streamlines in Fig. 2. The shock reflections can be seen in the two streams in the pressure map in Fig. 1. Similar features can also be seen for other cases in Fig. 3. This wind structure bears similarity to the under-expanded jets (e.g. Wilkes et al. 2012) and has also been observed in the previous simulations of wide-angled winds/jets (e.g. see Krause et al. 2012; Wagner et al. 2013; Mondal et al. 2022).

The backflows from the lateral streams are directed inwards i.e. in regions between the streams, as well as in the outward direction, as depicted using blue arrows in Fig. 1. Generally, the inward-directed backflows are weaker compared to the outward ones. As can be seen from the velocity streamlines shown in Fig. 2, these backflows lead to eddies/vortices in the cocoon. The eddies lead to velocity shears in the inner layers of the contact discontinuity, thereby facilitating some entrainment (or mixing) of the shocked ambient gas with the

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-874f4e88-d60f-4df6-9dac-de15f7556fad.jpg)

Figure 4. Logarithmic density (left), pressure (middle), and wind tracer (right) for the wide-angled dense ( 𝜂 𝑤 = 4 × 10 -1 ; top row) and highly dense wind ( 𝜂 𝑤 = 4; bottom row) with power 10 43 erg s -1 . The region enclosed inside the cocoon and forward shock (see Appendix A) are represented by white and yellow contours in the right panels.

Figure 5. Logarithmic density (left), pressure (middle), and wind tracer (right) for the light ( 𝜂 𝑤 = 4 × 10 -2 ; top row) and highly dense ( 𝜂 = 4; bottom row) winds with a narrow opening angle ( 𝜃 𝑤 = 35 ). The winds have a power of 10 43 ◦ erg s -1 . The region enclosed inside the cocoon and forward shock (see Appendix A) are represented by white and yellow contours in the right panels.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-57389146-eff1-4304-be08-7ef217f2645d.jpg)

Figure 6. Probability Distribution of magnetic field in the cocoon (tr1 &gt; 10 -3 ) and SAM (defined in Appendix A) for light winds W42, W43 and W44. The times of the snapshots are similar to those in Fig. 3, corresponding to when the head of the forward shock reaches close to the top of the simulation domain.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-6277d7fc-9bc2-45dc-97a3-d7cee6147837.jpg)

wind plasma. This causes wind tracer (tr1) values lower than 1 in these regions, which will be discussed in the next sections.

- · Dependenceonthewindpower: Thelogarithmic density, pressure, and wind tracer maps for simulations with different wind powers are shown in Fig. 3. The white and yellow contours in the right panels show the regions enclosed within the cocoon and forward shock, respectively. It can be seen that as the wind power decreases, the forward shock takes a more spherical shape. Conversely, high-power winds carry strong forward momentum, resulting in an elongated appearance along the 𝑍 direction. The SAM for W42 shows a doublelayered structure where the density is very high in regions closer to the cocoon, and decreases outwards. The inner SAM consists of the old swept-up ambient medium and hence depicts high gas density. The outer SAM also indicates that the forward shock from W42 is not as powerful as the high-power winds by the time it reaches a similar altitude. This can affect the acceleration of electrons at the forward shock and in turn, the polarized emission, which is discussed later in Sec. 3.2.2.

The backflows on the lateral sides of the cocoon generate eddies (as shown in Fig. 2), leading to a reduction in gas pressure in these regions. Simultaneously, a high-pressure region develops in the forward zone ahead of the Mach disk. The outward eddies act to confine the wind, resulting in collimation and instabilities in the fast streams from the edges of the barrel. Inward-directed backflows and highpressure zones ahead of the Mach disk, as observed in W43 and W44, counteract this collimation. However, this is not the case for W42 due to weak inward-directed backflows. This suggests that lowpower, light winds are likely to experience strong collimation from the surroundings, resulting in a considerably narrower cocoon, and conversely for the high-power ones.

The wind tracer maps in the right panels indicate that the cocoons undergo some mixing with the ambient gas. The tracer values are 1 inside the conical zone and along the fast lateral streams. High tracer values (tr1 &gt; 0 75) are also observed in the central and lateral . regions of the cocoon, where the backflowing wind plasma is present. Some mixing can be seen on the lateral sides, and is significant near the cocoon head, resulting in lower wind tracer values (tr1 ≲ 0 25). Initially, the formation of fast streams and strong backflows . in the lateral parts of the cocoon causes the ambient gas to fill most regions between the streams above the Mach disk. Inward-directed backflows cause shearing in the inner layers of SAM, leading to some entrainment of the external medium. This results in low wind tracer values near the top head of the cocoon. Stronger backflows in highpower winds cause more entrainment compared to low-power winds. The low-power wind W42 initially confines itself within a narrow cocoon, limiting the extent of its eddies and reducing mixing with the external medium.

- · Dependence on the density of the wind: Maps for different physical quantities for a dense and highly dense wind with a power of 10 43 erg s -1 are presented in Fig. 4 (referred to as W43\_L1000D and W43\_L1000HD in Table 1, respectively). It is evident that, as the wind density increases, the Mach disk in the winds is able to advance more along the 𝑍 -direction. This is because the dense winds carry a high momentum flux and are, thus, weakly affected by interaction with the surroundings. Consequently, the lateral streams manifest at later times when the wind has progressed further than the case of the light wind in Fig. 3, resulting in shorter streams. This implies that with increasing density, the winds become more stable, closely resembling the wide-angled wind model from Faucher-Giguère &amp; Quataert (2012).

It can also be seen that the cocoon of dense winds lacks the characteristic 'hump-like' feature at the head, as observed in the case of light winds in Fig. 3. While a small 'hump' is discernible in the dense wind, it is notably absent in the highly-dense case. In light winds, the lateral streams are unstable due to confinement force from eddies in the lateral regions, and thus the backflowing plasma is directed both inwards and outwards. The inward-directed backflows try to fill the volume ahead of the Mach disk, resulting in some plasma being pushed upward, creating a distinct 'hump-like' structure within the cocoon. In contrast, when the wind is dense, and thus, relatively stable, the lateral streams emanating from the edges of the barrel are comparatively stable. As a result, strong backflows from high-velocity regions occur primarily along the lateral sides and are weaker inwards, leading to the absence of a prominent hump at the cocoon's head.

The tracer maps displayed in the right panel indicate that some entrainment of the external gas (tr1 ≲ 0 25) is present near the cocoon . head of the dense wind. This results from the formation of fast lateral streams and inward-directed backflows at later stages. However, in the case of the highly dense wind, the lateral streams are still in their initial phase, and thus no entrainment of the shocked ambient gas is seen here.

- · Dependence on the opening angle of the wind: Maps for a light and highly dense wind with a half-cone opening angle ( 𝜃 𝑊 ) of 35 ◦ and power 10 43 erg s -1 are shown in Fig. 5. It can be seen that a narrow light wind experiences strong collimation from its surroundings when compared to the wide-angled wind of similar power in Fig. 3. Consequently, the wind becomes more constrained, resulting in the formation of a slender, elongated high-pressure zone ahead of the Mach disk. This also leads to high-pressure knots in these regions and a high-pressure cocoon head. Similar high pressure can also be seen near the cocoon head of the dense narrow wind, resulting from collimation from the surroundings. The dense wind is stable than the light wind here, and so fast lateral streams have not yet developed. The narrow winds undergo weak mixing with the external medium, as can be seen from the tracer map in the right panel.

Similar to the jets in Paper I, the bow shock of the wind compresses the external magnetic fields. Fig. 6 displays the probability density function (PDF) of the magnetic field within the cocoon and SAM across wind simulations with different powers. In all cases, the PDF for the SAM consistently lies on the left side of the cocoon, indicating comparatively weaker magnetic fields. However, it is noteworthy that the peaks in the PDF from the cocoon and SAM get close in magnitude as the wind power is decreased. This discrepancy arises due to the decay of the magnetic field within both the cocoon and the SAM as expansion occurs, as previously shown in Paper I. Consequently, the magnetic field within the cocoon of low-power winds, which are initiated with comparatively weaker magnetic fields (refer to Table 1), has significantly decayed to values comparable to those within the SAM. Conversely, the magnetic fields in the high-power winds, particularly for W44, are sustained at values higher than the SAMbythe time the cocoon's head reaches the top of the simulation domain. The distinct magnetic field strengths in the cocoon and SAM lead to varying effects on the emission and polarization originating from the wind, as elaborated in the subsequent sections.

## 3.2 Synchrotron emission and polarization from winds

In this section, we discuss the synchrotron emission and polarization fromthe wind simulations listed in Table 1. We present maps illustrating the logarithmic flux, polarization fractions ( 𝜋 𝑇 and 𝜋 𝑊 ), and the fractional change in polarization ( ( 𝜋 𝑊 -𝜋 𝑇 )/ 𝜋 𝑊 ) for winds with dif-

ferent powers at 𝜃 𝐼 = 90 ◦ (angle from the + 𝑍 -axis) and 𝜙 𝐼 = 0 (measured counter-clockwise from the X-Z plane) in Fig. 7. The polarization vectors are overlaid on the polarization map with black arrows. Similar to Paper I, 𝜋 𝑇 maps correspond to the total polarization obtained from the cocoon + SAM, and 𝜋 𝑊 shows the polarization from the cocoon only. The conditions for identifying the regions inside the forward shock and SAM are mentioned in Appendix A. The fractional change in polarization due to the contribution of the SAM along the LOS is shown using Δ 𝜋 = ( 𝜋 𝑊 -𝜋 𝑇 )/ 𝜋 𝑊 , i.e., Δ 𝜋 = 1 implies very strong depolarization due to the SAM, and 0 means no depolarization. Fig. 8 showcase maps from dense (W43\_L1000D) and highly dense (W43\_L1000HD) winds. Additionally, in Fig. 9, maps for a narrow light and dense wind (W43\_L1000N and W43\_L1000HDN) are presented. In the paper, we assume 𝜖 𝑊 and 𝜖 𝑆𝐴𝑀 = 0 1, implying . that the electrons carry exactly 10% of the fluid energy density in the cocoon and SAM for almost all the cases. However, the acceleration efficiencies for non-thermal electrons in the cocoon and SAM can vary, which is explored in Fig. 10. Figs. 11 and 12 showcase maps depicting winds with different powers and the dense wind case, respectively, for an image plane at 𝜃 𝐼 = 60 . Below, we discuss the ◦ common characteristics of synchrotron emission and polarization, focusing on maps from the light winds at 𝜃 𝐼 = 90 ◦ illustrated in Fig. 7. How different physical parameters (wind power ( 𝑃 𝑤 ), density contrast ( 𝜂 𝑤 ), opening angle ( 𝜃 𝑤 ), viewing orientation of the observer ( 𝜃 𝐼 ), etc.) can affect these features in the winds is discussed in the subsequent section (Sec. 3.2.2). The depolarization effect from the turbulent field in the SAM is discussed in Sec. 3.2.3.

## 3.2.1 General characteristics of emission morphology and polarization

## Morphology of emission:

As can be seen from Fig. 7, the general characteristics of wind emission at 𝜃 𝐼 = 90 ◦ are diffuse spherical extended emission with a bright cylindrical structure that is capped by a horizontal bright band. However, the shapes of these substructures can vary depending on different parameters and also the observer's orientation, which we discuss in subsequent sections. The high emission in the central parts originates from the high-pressure zones situated between the lateral streams. Notably, nearly horizontal arcs of high emission are visible near the cocoon head (refer to H1, H2, and H3). Such a feature in the light winds arises due to elevated pressure and high magnetic fields, particularly the 𝐵 𝑦 component of the wind's magnetic field in the 𝑌 -𝑍 planes 3 .

Additionally, shearing of the fields (particularly 𝐵 𝑦 and 𝐵 𝑧 ) can also cause bright filamentary structures (F1 and F2) in the regions between the lateral streams. The vertical filaments are primarily influenced by the poloidal component (i.e., 𝐵 𝑧 ) of the magnetic field, as can be seen from the middle row in Fig. B1. The lateral regions, characterized by backflows and eddies (refer to Fig. 1), depict low pressure, consequently exhibiting weak emission.

Distribution of polarization: The polarization distribution ( 𝜋 𝑇 ) shows high polarization primarily from the forward shock, the central regions of the wind, and the horizontal arc, which coincides with a bright band in the synchrotron maps. High polarization at the forward shock arises from the compression of the external magnetic fields parallel to the shock. This outcome is similar to the findings observed

3 It is important to note that only the synchrotron emission resulting from the field components perpendicular to the line of sight contributes to the total emission at the designated image plane.

for jets in Paper I. The polarization vectors align perpendicularly to the shock surface, indicating the compression of magnetic fields parallel to it. Within the SAM, there is a gradual decay of turbulent fields, resulting in low polarization levels due to the small-scale coherence of the fields. Regions above and below the high-polarized bright horizontal arc also display low-polarization values. This is due to turbulent fields in the SAM, and backflows in regions between lateral streams respectively in upper and lower zones. The poloidal fields between the lateral streams lead to low polarization in the cocoon, which is caused by the vector cancellation of orthogonally polarized components along LOS.

In 𝜋 𝑊 (also 𝜋 𝑇 ) maps, there is a prevalence of high polarization, attributed to the injected toroidal magnetic field within the wind, as indicated by the vertical polarization vectors. However, certain regions near the cocoon head and lateral sides exhibit low polarization values, resulting from the generation of poloidal fields. The poloidal fields are introduced due to velocity shears caused by backflows in these regions (see Fig. 1). Additionally, the mixing with the surrounding SAM, facilitated by eddies in the backflows, introduces further poloidal fields, which are subsequently amplified by turbulent motions within the cocoon (see 𝐵 𝑧 map in Fig. B1). Such shearing of the fields near the contact surface leads to high polarization, particularly due to the small integration length in these regions. Consequently, the polarization vectors align perpendicularly to the projected surface of the cocoon in 𝜋 𝑊 maps.

## 3.2.2 Effect of different parameters on the observable features

· Dependence on the wind power: The emission maps in Fig. 7 indicate that the bright region in the central parts gets wider with increasing wind power. This, as shown in Fig. 3, is attributed to the weak collimation of the lateral streams, resulting in wide highpressure zones above the Mach disk for the high-power winds. As a result, the bright horizontal arc is more extended as the wind power is increased. These arcs near the cocoon's head display a double-layered appearance for light low-power winds here, which is clearly visible for W42andW43andislessprominentforW44.Theupperarcisformed due to the shearing of the 𝐵 𝑦 component near the 𝑌 -𝑍 midplane, resulting from inward-directed backflows. This is evident from the evolution of the 𝐵 𝑦 field displayed in the top row of Fig. B1, where fields are sheared to form an elongated, nearly horizontal structure. The lower arc is a consequence of the 𝐵 𝑦 component of the toroidal magnetic field of the wind, coinciding with high-pressure zones at the top edges of the lateral streams. We confirm that such features may not be as distinctly observable from other lines of sight (e.g., in the 𝑋 -𝑍 image plane), as it depends on the instantaneous flows in the wind cocoon. As indicated in Table 1, the same magnetization value of 0.1 across all winds results in relatively weak magnetic fields for low-power winds in comparison to high-power ones. Consequently, as the cocoon expands, the local magnetic field gradually decreases, reaching a level comparable to that in the SAM (Fig. 6). This results in nearly uniform emissions from the W42 wind. Conversely, highpower wind cocoons exhibit stronger magnetic fields than the SAM. Consequently, the emissions from the central parts of such cocoons are considerably higher compared to those emanating solely from the SAM.

From the polarization maps ( 𝜋 𝑇 ) in Fig. 7, it can be seen that the high polarization at the forward shock is particularly evident in the W43 and W44 cases, but not in W42 due to weak forward shock (see Fig. 3). The polarization map in W42 shows larger areas of low polarization when compared to cases with higher power. This occurs as the SAM in W42 is broader, which enhances the decay and

Figure 7. Logarithmic synchrotron flux (first row), total and cocoon polarization fraction (second and third row), and the fractional change in polarization (Fourth row) at 𝜃 𝐼 = 90 ◦ image plane for different simulations. H1, H2, and H3 show horizontal bright arcs near the cocoon's head. F1 and F2 show filamentary-like structures resulting from the shearing of magnetic fields.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-6c7772eb-5fe1-45ed-9b75-9be295a47d19.jpg)

Figure 8. Top: Logarithmic synchrotron flux (left), total and cocoon polarization fraction middle and right) at 𝜃 𝐼 = 90 ◦ image plane for dense wind case (W43\_L1000D). H4 shows a horizontal bright arc near the cocoon's head. Bottom: Same as above for highly dense wind (W43\_L1000HD). H5 shows a high-emission extended region near the cocoon's head.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-17d9e1b0-2c3d-4d17-b879-adc2af3901c9.jpg)

Figure 9. Top: Logarithmic synchrotron flux (left), total and cocoon polarization fraction (middle and right) at 𝜃 𝐼 = 90 ◦ image plane for narrow light wind case (W43\_L1000N). H6 shows a high-emission region near the cocoon's head, and N1 are the high-pressure knots that form due to the strong collimation of the wind. Bottom: Same as above for highly dense narrow wind (W43\_L1000HDN). A bright compact region (N2) can be seen above the horizontal bright arc (H7).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-f8021d2e-5fd8-4335-aeef-e5552c689c12.jpg)

Figure 10. Top: Logarithmic synchrotron flux (left), total polarization fraction (middle) and fractional change in polarization (right) at 𝜃 𝐼 = 90 ◦ image plane for W43\_L1000. The figures in the top row correspond to a stronger acceleration of electrons in the cocoon than the SAM ( 𝜖 𝑊 &gt; 𝜖 𝑆𝐴𝑀 ), and conversely, in the bottom row.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-b27226a9-d3ab-42e5-a3ec-5a789b331e12.jpg)

disorder of the magnetic fields within it. Consequently, the bright arc near the cocoon head is notably visible in the W42 case.

- · Dependence on the density of the wind: In order to compare how the wind density contrast can affect the observable features, we focus on the light, dense and highly dense wind cases with power 10 43 erg s -1 . The maps from the light wind case are displayed in the middle column in Fig. 7. The emission and polarization maps for the dense and highly dense winds at 𝜃 𝐼 = 90 ◦ image plane are illustrated in Fig. 8. It can be seen that the bright and highly polarized arc at the cocoon's head appears more extended for the dense winds when compared to the light wind case. This occurs due to more stability of the dense winds (as discussed in Sec. 3.1) when compared to the light winds. Also, the absence of a double-layered structure in the horizontal bright band, as observed in lighter wind cases, is evident here. This absence, as illustrated in the bottom row of Fig. B1, is attributed to the lack of enhancement in the 𝐵 𝑦 component within the 𝑌 -𝑍 midplane, a phenomenon resulting from inward-directed backflows in lighter wind cases. Consequently, the 𝐵 𝑦 component of the toroidal field, coinciding with high-pressure zones at the top edges of the lateral streams (as depicted in the 𝐵 𝑥 map in the bottom row), contributes to the formation of the bright horizontal arc when integrated along the line of sight in this scenario. Some emission can

be seen between the Mach disk and horizontal arc for dense wind, giving rise to a 'mushroom-shaped' morphology. This, however, is not seen for the highly dense wind, as the whole regions above the Mach disk contribute to the wide bright band here.

Also, the horizontal arcs are quite prominent in the dense wind cases when compared to the light wind case. It can be seen that the emission from the central parts of the light wind is almost uniform, making the arcs less prominent. Inward-directed backflows in the light winds introduce toroidal as well as poloidal components between the lateral streams (see the middle row in Fig. B1), leading to almost homogeneous emission from regions above the Mach disk, which however is not the case for the dense winds (bottom row in Fig. B1). Contrarily, backflows in the dense wind mainly occur in the lateral sides of the high-velocity streams, which also generate poloidal fields. This leads to low polarization values in these regions, as can be seen from the 𝜋 𝑊 maps.

- · Dependence on the opening angle of the wind: Fig. 9 shows the emission and polarization maps for narrow light and highly dense winds, respectively, observed at an angle of 𝜃 𝐼 = 90 . The collima-◦ tion of lateral streams in the light wind case leads to bright knots (labeled as N1). Both cases exhibit intense compact emissions labeled as H6 and N2, stemming from high-pressure zones near the

Figure 11. Logarithmic synchrotron flux (top), total polarization fraction (bottom) at 𝜃 𝐼 = 60 ◦ image plane for simulations with different wind powers.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-92a264ef-f279-4019-b2b8-f40e54460c37.jpg)

Figure 12. Logarithmic synchrotron flux (left) and total polarization fraction (right) for W43\_L1000D at 𝜃 𝐼 = 60 ◦ image plane.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-e85b0dfe-b200-4cd7-bcee-9889a2936ecb.jpg)

cocoon head, which also coincide with high-polarization values. The highly dense wind case displays a horizontal bright and highly polarized arc from regions above the Mach disk (see H7 labeled region). However, in this instance, the width of the arc is relatively narrower, as the regions above are more constrained, resulting in a bright compact spot (N2) ahead.

- · Different particle acceleration scenarios: In Fig. 10, the emission and polarization maps correspond to varying values of 𝜖 for the cocoon and SAM, i.e. implying different acceleration efficiency of the electrons in these regions. The top row shows the maps for higher 𝜖 𝑊 (0.1) than 𝜖 𝑆𝐴𝑀 (0.01); implying higher acceleration in the cocoon than SAM, and conversely in the bottom row. The electrons within the cocoon experience acceleration primarily at the Mach disk, while in

the SAM, acceleration can be expected at the forward shock. When the efficiency of acceleration in the cocoon is notably high in the presence of magnetic fields (for magnetically driven winds, see e.g. Fukumura et al. 2015), it results in stronger radio emissions originating from the cocoon compared to that from the SAM (as shown in the top row of Fig. 10). Conversely, if the SAM exhibits substantially higher acceleration efficiency (for radiation-pressure-driven winds, see e.g. Nomura et al. 2013), the observable radio emission characteristics from the wind cocoon become less discernible (see bottom row).

The polarization maps ( 𝜋 𝑇 ) indicate that when the wind is strong in accelerating electrons in the cocoon, most of the regions from the wind display high polarization. Conversely, the region exhibiting low

polarization from the wind becomes more pronounced in the bottom row. This observation suggests that for radiation-driven winds, where the magnetic fields inside the winds may not be very strong, it is improbable to obtain significant radio polarization from the wind. However, high polarization can be observed from the forward shock region in SAM only if the emission is strong enough to be detected.

Moreover, it should be noted that the acceleration efficiency of the electrons at the forward shock can also vary in time. For example, as can be seen from density and pressure maps in Fig. 3, the forward shock in W42 has considerably weakened compared to the other cases, developing a double-layered structure. This weakening diminishes the impact of the outer shock on the compression of external magnetic fields, consequently lowering the acceleration efficiency of electrons in this region. This leads to low polarization values at the forward shock for W42, as can be seen from Fig. 7.

- · Dependence on the viewing angle: At an inclined view of 𝜃 𝐼 = 60 ◦ in Fig. 11 and 12 the emission from the central parts of wind appears 'hourglass-shaped'. The majority of this emission originates from regions situated above the Mach disk, with additional contributions coming from the areas surrounding the barrel zone, as illustrated at 𝜃 𝐼 = 90 ◦ in Fig. 7. Consequently, in light winds (see e.g. W42), these structures are narrower due to the strong collimation and relatively broader for the high-power light winds. The injected toroidal magnetic field results in high polarization from these regions in all the cases. In the surrounding areas with weak flux, low polarization is attributed to vector cancellation arising from orthogonally polarized components, along with contributions from the SAM. It should be noted that an inclined LOS, the toroidal field components are not uniformly perpendicular or parallel to the line of sight. This also add on to the low-polarized regions seen around the central spines. The distinctiveness of low polarization in W42 is emphasized due to its broader SAM, leading to rapid decay and disorder of the turbulent fields. The broad annular rings at the top of the 'hourglass' structure (particularly visible on the approaching side for W43 and W44) represent high-pressure annular regions above the edges of the lateral streams coinciding with high toroidal field components. Due to weak inward-directed backflows in the dense wind when compared to the light winds, any emission from the regions between the arcs is considerably weaker. Nonetheless, the bright outer arcs display high polarization in all cases. These arcs are distinctly visible for W42, whereas they seem merged with the high-polarization regions from the SAM in others.

## 3.2.3 Depolarization effect of the SAM

- · Dependence on the wind power: The fourth row in Fig. 7 depicts the fractional change in the polarization of the cocoon ( Δ 𝜋 ) after the contribution from the SAM. More reddening for the W42 clearly indicates that the presence of a turbulent field in the SAM leads to more lowering in the polarization from the cocoon for the lowpower winds. This is because the magnetic fields in the SAM are quite comparable to the cocoon for the low-power winds, as can be seen from the PDFs of the magnetic fields in Fig. 6. Contrarily, the highpower winds carry stronger magnetic fields than the SAM, and thus the depolarization effect is weaker. As found in Paper I, the magnetic fields in the cocoon and SAM decay with the expansion. When the fields in both regions become comparable, the depolarization from the SAM gets strengthened.
- · Dependence on 𝜖 𝑊 : It can be seen from Δ 𝜋 maps in Fig. 10, the depolarization from the SAM becomes prominent for the case with 𝜖 𝑆𝐴𝑀 = 0 1 and . 𝜖 𝑊 = 0 01. In the context of a radiation-driven . wind, the non-thermal radio emission from the wind cocoon may not

be very strong. Under such circumstances, the SAM's contribution takes precedence in the total non-thermal emission, as illustrated in the left panels of the figure. Consequently, the diminished polarization of the wind primarily results from the contribution of turbulent magnetic fields within the SAM.

## 4 COMPARING SYNCHROTRON AND POLARIZATION CHARACTERISTICS FROM JETS AND WINDS

The primary motivation of this section is to identify some noticeable observational features for the non-thermal synchrotron emission and polarization to distinguish between the jets and winds. Emission and polarization maps depicting jets and winds with varying powers and viewing angles are presented in Figs. 13-15 and Figs. C1-C2. To compare with observed results, these maps are convolved with a Gaussian filter with a different standard deviation of 100 pc and 300 pc and truncated at twice these scales 4 . The corresponding FWHM gives an angular resolution of 125 and 375 milliarcsec, respectively, for sources at a redshift ( 𝑧 ) = 0.1 (assuming 𝐻 0 = 67 kms -1 Mpc -1 , and Ω 𝑀 = 0 3). The latter values are comparable to the resolution . scales reached with previous studies (see e.g. Jarvis et al. 2019, 2021; Pierce et al. 2020) and can be achieved in the local universe (Baldi et al. 2018). For a larger redshift of 2, they correspond to a resolution of around 27.4 and 82 mas, which is typically not achieved by the studies looking at this scale (Alexandroff et al. 2016; Rosario et al. 2021); although can be possible using high-resolution data (Njeri et al. 2024). In these figures, we mainly focus on one order in magnitude for the flux dynamic range in the fiducial maps presented here and also show maps for the two-order-of-magnitude flux range for W42andW44.Thepolarization vectors are depicted as black arrows. These snapshots correspond to times when the jet or wind head approaches near the top of the 𝑍 axis, and the forward shock is confined within the box. The correlation length of the external magnetic field is 1000 pc in all the simulations, except for J42, where, we ran the simulation with only a 500 pc coherence scale in Paper I.

## 4.1 Characteristics in emission morphology

The winds and jets display some common traits as well as some differences in emission morphology, which are listed below:

- · In the emission morphologies, it is evident that bright features emerge from the central regions in both jets and winds, which are surrounded by weak extended emission (see Figs. 13 and 15). A jet has a systematically narrow spine with defined hotspot, or hotspot complexes for kink-unstable jets (see h1, h2, and h3). Winds, on the other hand, display broader spines or 'hourglass-shaped' bright emission. The bright annular rings/arcs observed at an inclined view of 60 ◦ (R1, R2, and R3) manifest into bright horizontal arcs near the edges of the cocoon at 90 ◦ (H2 and H5). Such rings/arcs are not seen for jets and are a particular feature in the winds.
- · Winds exhibit a diverse range of morphologies depending on their power and density. Light and low-power winds tend to become collimated, unlike high-power winds. As a result, the presence of a hotspot may not reliably distinguish jets from winds, as pronounced

4 This implies FWHM of 240 and 720 pc, and truncated at scales of 400 and 1200 pc, respectively.

Figure 13. Logarithmic synchrotron flux from jet simulations in Paper I and wind simulations in this study, respectively, is presented at the 𝜃 𝐼 = 60 ◦ image plane. The color bar range is maintained at one order of magnitude in logarithm in all figures except for the top and bottom right ones, where a two-orders-of-flux range is used for comparison. Regions with flux values lower than the minimum in the color bar are removed, and the flux is convolved with a Gaussian filter with an FWHM of 240 pc. The hotspots in J42, J43 and J44 are indicated as h1, h2 ('multiple') and h3, respectively. F3 represents a bright filamentary-like structure in J43. R1, R2, and R3 represent the bright annular rings near the cocoon head in the winds.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-4aafb8ab-fd7a-4f7f-87c1-ceba45662638.jpg)

collimation in low-power winds can produce similar features (as observed for W42 in Fig. 13). We have confirmed that such strong collimation and hotspot-like regions do not form in highly dense ( 𝜂 𝑤 = 4) wide-angled winds with a power of 10 42 erg s -1 . Additionally, the emission closely resembles that of the dense wind case illustrated in Fig. 12. This suggests that dense winds are more stable and thus, exhibit prominent bright annular rings (and bright distinguished hor- izontal arcs at 𝜃 𝐼 = 90 , as inferred from Fig. 15) compared to lighter ◦ winds. However, it is worth noting that the hotspot-like features can still form in dense narrow-angled winds, as depicted in Fig. 9.

- · Effect of dynamic range: To illustrate the impact of dynamic range, we present the emission maps for W42 and W44 using a twoorders-of-magnitude flux range in Fig. 13 (see top and bottom right plots). Notably, the extended emission is narrower for jets than for

## Polarization fraction (Tr)

Figure 14. Convolved total polarization fraction ( 𝜋 𝑇 ) corresponding to the emission maps displayed in Fig. 13. The Gaussian filter used for convolution has a FWHMof 240 pc.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-66169757-016b-478b-b908-98bdfba3b346.jpg)

winds. It can be observed that the emission from regions near the forward shock in the SAM is lower compared to that from the central parts of the wind. Consequently, while the central bright regions are unlikely to be missed when observations are limited to one order of magnitude in dynamic range, outer regions, which are part of the SAM, can be overlooked.

- · Effect of resolution: The extent of bright emission from the central parts varies in both jets and winds with resolution. To facilitate comparison, we calculated the average width of the regions surrounding the central vertical axis, encompassing 80% of the maximum flux

at the given 𝑍 ′ value, within a range of 𝑍 ′ = ± 1 kpc in Fig. 13. At high resolutions, these widths measure 420 pc and 310 pc for J42 and J44, respectively, and approximately 1.2 kpc and 1.3 kpc for W42 and W44. In low-resolution maps shown in Fig. C1, the widths of the bright spines have increased to quite comparable values of 920 pc and 1.2 kpc for J42 and W42, while for J44 and W44, these are 630 pc and 1.36 kpc, respectively. Furthermore, it is apparent that at low resolution, several features, such as the spine, hotspots, filaments, and arcs, have become less prominent. This leads to some resemblance in appearance between jets and some wind models. For example, in

Figure 15. Logarithmic synchrotron flux (top) and total polarization fraction (bottom) for J43 and W43 at 𝜃 𝐼 = 90 ◦ image plane. Both emission and polarization maps correspond to a one-order-of-magnitude flux range and are convolved with a Gaussian filter with an FWHM of 240 pc. The multiple hotspots in J43 are shown as h2, and the bright horizontal arcs in W43(HD) are presented as H2 and H5.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-75584b2f-b77f-4bd9-8fc4-6e8bcad7e39a.jpg)

Table 2. Selection criteria to distinguish between jets and winds. These traits are discussed in detail in Sec. 4.

| Characteristic features   | Jets                                              | Winds                                                                                    | Comment                                                                                                                                  | Criteria   |
|---------------------------|---------------------------------------------------|------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------|------------|
| Regions of highest flux   | Spine and hotspots(s)                             | Similar, but broader spine and arcs                                                      | Resemblance for some models with low resolution                                                                                          | Maybe      |
| Width of bright emission  | Narrow                                            | Comparatively broader                                                                    | Can be comparable for low resolution                                                                                                     | Maybe      |
| High polarization         | Localized: hotspots, spine                        | Wide arcs and spine                                                                      | Persist with lowering resolution, but low-values dominate                                                                                | Yes        |
| Depolarization in cocoon  | More                                              | Less                                                                                     | Same as above                                                                                                                            | Yes        |
| Collimation               | Mostly collimated, cocoon width varies with power | Dense winds are broader; light winds (especially low power) can undergo some collimation | Low-power jets have broader cocoons, while low-power winds, due to collimation, have narrower cocoons than their high-power counterparts | Maybe      |

light winds, bright 'hourglass-shaped' regions appear close to the spine of the jets. Interestingly, bright regions in W42 closely match the emission morphology from the jet in J44. Although a hotspot in W42 might resemble that of a jet, it is important to recognize the absence of a bright, thin spine typically associated with jets, which becomes evident at higher resolutions. Since the emission morphology is often used to distinguish between jets and winds, our study shows that some features can become blended when poor resolution is used. This underscores the importance of higher-resolution radio observations; however, noting that 'only' utilizing the very highest resolution images may resolve away important diffuse structures that also contain important diagnostic information in radio morphology and polarization.

## 4.2 Distribution of polarization

Below, we list some common characteristics and differences in the polarization structure from the jets and winds.

- · Some notable features can be clearly seen from the maps depicted in Figs. 14 and 15. High polarization is mainly detected along the central spine and hotspots/arcs, which are localized in jets and

are more widespread in the winds. Additionally, high polarization is observed along the filament (F3) within the jet, although it is reduced due to smoothing.

- · In jets, particularly those with a higher power, the hotspots exhibit noticeable high-polarization values, typically found near the upper edges of the cocoon. These hotspot(s) exhibit compactness, resulting in polarization originating from them being concentrated at specific locations (e.g. refer to h1, h2 h3). However, in winds, high polarization mainly occurs from the central parts and top edges of the cocoon (refer to R1, R2, R3, H2, and H5 from emission maps). These arcs are surrounded by comparatively low-polarization values (e.g. see maps for W42 and W43 in Fig. 14). The lateral edges of the maps can also display high polarization values, as can be seen particularly in the middle row left panel in Fig. 14.
- · In Fig. 14, both winds and jets display low polarization values in their extended regions, which mainly arises because of the vector cancellation of orthogonally polarized components, and also some contributions from the SAM. It is clear that in the jets' cocoon, where stronger poloidal fields are generated due to pronounced backflows, experiences more internal depolarization compared to winds (e.g. see 𝜋 𝐽 maps in Fig. 9 of Paper I and 𝜋 𝑊 in Figs. 7-8). The higher prevalence of low polarization in W42 is linked to its wider SAM and

relatively narrower cocoon, as depicted in Fig. 3, resulting in more decay of turbulent fields.

- · Effect of dynamic range: The top and bottom right plots for total polarization (see Fig. 14) indicate that a two-order-of-magnitude flux range enables the observation of polarized emission from the outer parts of the SAM. This remains invisible with a lower dynamic range, as shown in the plots on the left side. These regions may exhibit either high or low polarization values in different wind power simulations. In the case of W42, the weak strength of the forward shock leads to low polarization, whereas relatively higher polarization is seen in W44. The lateral regions in W44 (which corresponds to SAM) also display high-polarization values, similar to those seen for the jets (see left panels), which were not observed at the lower dynamic range. At 𝜃 𝐼 = 90 , high values like these are not observed in 1-dex ◦ range maps (see Fig. 15). However, we confirm that they are seen, primarily near the bow shock's head, in 2-dex range maps.
- · Effect of resolution: In Fig. C2, it is evident that decreasing the resolution has not altered the qualitative distribution of polarization characteristics. However, low-polarization values have become more prevalent. Polarization values near the spine and hotspots in jets appear lower compared to the high-resolution maps. In winds, the highly polarized arcs near the cocoon head seem to blend with regions from the SAM, and the resulting values are diminished after smoothing. This suggests that polarization can be an important factor in distinguishing between the jets and winds.

## 5 SUMMARY AND DISCUSSION

AGN-driven jets and winds are believed to be launched by distinct mechanisms (Königl 2006). Numerous simulations have been performed for both (see e.g. Proga &amp; Kallman 2004; Perucho &amp; Martí 2007; Fukumura et al. 2015; Perucho et al. 2022), shedding light on various aspects relevant to real systems. It is found that the evolution and stability of jets are influenced by factors such as Mach number, Lorentz factor, and density contrast with respect to the ambient medium (Bodo et al. 1994; Perucho et al. 2005). Conversely, the understanding of winds remains less clear, with their origin being a subject of ongoing debate. It is plausible that a single mechanism may not suffice to explain the observed ultrafast outflows (UFOs) in quasars, suggesting the contribution of multiple processes (Wang et al. 2022). This underscores the necessity for incorporating more sophisticated physics and a deep understanding of physical conditions in AGNs to comprehensively account for the winds.

Observations have identified clear indications of AGN-driven winds, whether through direct or indirect means. For instance, an inverse relationship between the radio loudness and the column density of ionized wind is found in radio-loud AGN by Mehdipour &amp; Costantini (2019). This connection suggests a shared magnetic mechanism responsible for initiating both wind and jet formations (see e.g. King et al. 2013). Excess radio emissions within several quasars also appear to be predominantly influenced by the presence of either winds or small-scale jets (Zakamska &amp; Greene 2014; Jarvis et al. 2019; Petley et al. 2022; Andonie et al. 2022).

Theoretical simulations can help establish a connection with observational findings, bridging the gap in our understanding of AGNdriven jets and winds. In our work (Paper I and this paper), we have attempted to gain insights into the evolution of jets and winds within a medium enriched with turbulent magnetic fields. We have examined the behavior of these at different energies, all of which are initiated with a toroidal magnetic field in our studies. Our analysis reveals that the evolution of jets and winds follows different behaviors, leading to observable distinctions in their features, which we summarize in the following section.

## 5.1 Dynamical and observational perspectives on jets and winds

Asindicated in Sec. 3.1 of this study and Sec. 4.1 in Paper I, the evolution of jets and winds exhibits distinct evolution characteristics. Jets, characterized by high beaming and relativistic speeds, advance more rapidly in their surroundings compared to winds. In contrast, winds exhibit mild relativistic speeds with a wide-angle spread, providing a larger surface area for interaction with the surrounding medium. This wide-angle nature results in instabilities at the periphery of the barrel zone, particularly evident in light winds, leading to lateral streams from the sides of the barrel (see Fig. 1). These streams influence wind evolution by inducing backflows ahead of the Mach disk and contributing to mixing with the SAM. Dense winds, in contrast, demonstrate stability and propagate without being much affected by their surroundings. The emission and polarization features from the winds with varying parameters (such as power, density, and opening angle) are discussed in detail in Sec. 3.2.

We have analyzed polarization and emission from jets and winds with different models and compared them in Sec. 4. The insights from the observable distinguishing traits between jets and winds are outlined in Table 2. We find that several characteristic features, such as hotspots and spines, may not be good indicators for poorly resolved sources, highlighting the need for highly resolved observations. The polarization distribution appears to be a reliable indicator for distinguishing between jets and winds, as it is not significantly affected by resolution. Jets exhibit compact, bright emissions from narrow spines and hotspots, aligning with high polarization values. Winds display similar widespread features, with the polarized spine nearly hourglass-shaped and capped with polarized bright arcs. Additionally, our studies suggest that depolarization, due to poloidal fields generated in the cocoon, is more dominant in jets compared to winds. It should be noted that both jets and winds are initially launched with a toroidal magnetic field in our study; nevertheless, a poloidal component at the base region is expected to exist, giving rise to helical magnetic fields. Consequently, the amplification of the poloidal field is likely to be much more pronounced in jets compared to winds. Contrarily, a weak poloidal component is seen in winds, with a dominance of toroidal field elsewhere in the cocoon. In several radio sources, such transverse or toroidal components at the base of the jets (with poloidal fields) are associated with AGN wind (Silpa et al. 2021a,b). The jets remain collimated over longer spatial scales, although the cocoons widen as jet power decreases. However, denser winds consistently remain broad. Lighter, low-power winds can become unstable at earlier stages, leading to collimation by the surroundings and resulting in narrower cocoons than those of their high-power counterparts.

## 5.2 Limitations and future aspects

In our studies (Paper I and this work), we have focused on the evolution of magnetized compact collimated jets and wide-angled winds in a hot halo carrying turbulent magnetic fields. However, in realistic scenarios, such compact jets/winds are expected to interact with the multiphase ISM of the host galaxy. This interaction, in turn, may also influence their synchrotron and polarization characteristics, which we reserve for future exploration. In addition to investigating the morphology of emission and polarization, simulations can also provide valuable quantitative parameters, such as spectra originating

from the acceleration of non-thermal electrons within the jets and winds. Particle acceleration processes are likely to differ between jets and winds. For example, significant acceleration in jets may occur along the spine and at the hotspot, while in winds, the Mach disk can accelerate particles. Furthermore, differences in the strength of forward shocks and SAM properties in both jets and winds may have varying effects on particle spectra originating from the SAM. This aspect is reserved for future research, wherein non-thermal electrons (Lagrangian Particle in pluto: Vaidya et al. 2018) will be introduced into the simulation domain and their evolution will be followed. Theoretical results can subsequently be juxtaposed with direct predictions from recent observational missions like LOw-Frequency ARray (LOFAR) 5 , eMERLIN 6 , Event Horizon Telescope (EHT) 7 , and Very Large Array (VLA) 8 , thereby facilitating a more clear understanding of jets and winds.

## ACKNOWLEDGEMENTS

We express our gratitude to the Accordo Quadro INAF-CINECA 2017 and 'Pegasus' HPC facility at IUCAA, Pune, India, for high-performance computing resources. We also acknowledge the assistance of the INDO-ITALIAN mobility grant (INT/Italy/P37/2022(ER)(G)). GB and PR acknowledge support by Next Generation EU through PRIN MUR 2022 (grant n. 2022C9TNNX) and by the INAF Theory Grant 'Multi-scale simulations of relativistic jets'. C.M.H acknowledges funding from a United Kingdom Research and Innovation grant (code: MR/V022830/1). MM expresses gratitude to Ankush Mandal for the valuable discussions.

## DATA AVAILABILITY

The simulations generated for this study will be shared upon a reasonable request to the corresponding authors.

## REFERENCES

| Alexandroff R. M., Zakamska N. L., van Velzen S., Greene J. E., Strauss M. A., 2016, MNRAS, 463, 3056 Andonie C., et al., 2022, MNRAS, 517, 2577   |
|----------------------------------------------------------------------------------------------------------------------------------------------------|
| Baldi R. D., et al., 2018, MNRAS, 476, 3478                                                                                                        |
| Balsara D. S., Spicer D. S., 1999, Journal of Computational Physics, 149, 270 M. C., McKee C. F., Shields G. A., 1983, ApJ, 271, 70                |
| Begelman Blandford R. D., Ostriker J. P., 1978, ApJ, 221, L29                                                                                      |
| Blandford R. D., Payne D. G., 1982, MNRAS, 199, 883                                                                                                |
| Blandford R. D., Znajek R. L., 1977, MNRAS, 179, 433                                                                                               |
| Bodo G., Massaglia S., Ferrari A., Trussoni E., 1994, A&A, 283, 655                                                                                |
| Braito V., et al., 2007, ApJ, 670, 978                                                                                                             |
| Condon J. J., Kellermann K. I., Kimball A. E., Ivezić Ž., Perley R. A., 2013, ApJ, 768, 37                                                         |
| Falcke H., Nagar N. M., Wilson A. S., Ulvestad J. S., 2000, ApJ, 542, 197                                                                          |
| Faucher-Giguère C.-A., Quataert E., 2012, MNRAS, 425, 605                                                                                          |
| Fukumura K., Tombesi F., Kazanas D., Shrader C., Behar E., Contopoulos I., 2015, ApJ, 805, 17                                                      |
| Gallimore J. F., Axon D. J., O'Dea C. P., Baum S. A., Pedlar A., 2006, AJ, 132, 546                                                                |
| Girdhar A., et al., 2022, MNRAS, 512, 1608                                                                                                         |

5 https://lofar-surveys.org/

6 https://www.e-merlin.ac.uk/

7 https://eventhorizontelescope.org/

8 http://www.vla.nrao.edu/

| Harrison C. M., Ramos Almeida C., arXiv:2404.08050                                                                                                                           | 2024, arXiv                                                                                                                                                                  |
|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Jarvis M. E., et al., 2019, MNRAS, 485, 2710 Jarvis M. E., et al., 2021, MNRAS, 503, 1780 Kharb P., Silpa S., 2023, Galaxies, 11, 27 King A. L., et al., 2013, ApJ, 762, 103 | Jarvis M. E., et al., 2019, MNRAS, 485, 2710 Jarvis M. E., et al., 2021, MNRAS, 503, 1780 Kharb P., Silpa S., 2023, Galaxies, 11, 27 King A. L., et al., 2013, ApJ, 762, 103 |

```
                                                                                                                                                                                                       
                                                                                                                                                                                                        (
                                                                                                                                                                                                       (                                                                                                                                                                                                        with a                                                                                                                                                                                                        and C:                                                                                                                                                                                                        one-or                                                                                                                                                                                                        or                                                                                                                                                                                                        out                                                                                                                                                                                                        ine-or                                                                                                                                                                                                     out  
                                                                                                                                                                                                     out  

                                                                                                                                                                                                        +                                                                                                                                                                                                        -                                                                                                                                                                                                        end                                                                                                                                                                                                        int                                                                                                                                                                                                        for                                                                                                                                                                                                        if                                                                                                                                                                                                        all
```

## APPENDIX A: IDENTIFICATION OF FORWARD SHOCK AND SAM

The SAM region in different simulations is identified using the pressure ( 𝑃 ), density ( 𝜌 ), and jet-tracer (tr1) as:

```
W { \text{42} } \colon ( P > 2 \times P _ { 0 } \text{ or} r 1 > 0 \text{ or} \rho > 0. 0 0 7 \text{ cm} ^ { - 3 } ) \ & \text{tr} 1 \leq 1 0 ^ { - 3 } \\ W { \text{43} } ( N ) \colon ( P > 2 \times P _ { 0 } \text{ or} r 1 > 1 0 ^ { - 2 0 } ) \ & \text{tr} 1 \leq 1 0 ^ { - 3 } \\ W { \text{43} } \text{HD} \colon ( P > 1. 6 \times P _ { 0 } \text{ or} r 1 > 1 0 ^ { - 2 0 } ) \ & \text{tr} 1 \leq 1 0 ^ { - 1 } \\ W { \text{43} } \text{HD} \colon ( P > 2 \times P _ { 0 } \text{ or} r 1 > 0 ) \ & \text{tr} 1 \leq 1 0 ^ { - 4 } \\ W { \text{44} } \colon ( P > 4 \times P _ { 0 } \text{ or} r 1 > 1 0 ^ { - 2 0 } ) \ & \text{tr} 1 \leq 1 0 ^ { - 3 } \\ \text{ Here, } P _ { 0 } \text{ is the initial pressure in the simulation domain. The}
```

Here, 𝑃 0 is the initial pressure in the simulation domain. The conditions in the brackets select the regions inside the forward shock. The SAM is then identified by the wind tracer limit at the end, which is also used to identify the regions inside the wind cocoon.

## APPENDIX B: EVOLUTION OF MAGNETIC FIELDS IN THE LIGHT AND HIGHLY DENSE WINDS

Fig. B1 illustrates the time evolution of the 𝐵 𝑦 component of the magnetic field for the light wind W43 in the top row. It can be seen that the inward-directed backflows cause shearing and amplification of 𝐵 𝑦 fields near the cocoon's head. The substructures of fields from opposite sides can combine, forming an elongated amplified magnetic component. Such phenomenon gives rise to the upper bright horizontal arcs in the light wind cases, depicted in Fig. 7.

In the middle and bottom rows, we have shown the different magnetic field components in the 𝑌 -𝑍 midplane for the light and highly dense winds, respectively. It is clear that the evolution of wind differs between these cases, as elaborated in Sec. 3.1, resulting in varying magnetic field structures. The 𝐵 𝑥 component represents the injected toroidal field and exhibits enhanced values in both cases. As shown in the top row, the inward-directed backflows enhance the 𝐵 𝑦 component of the field near the cocoon head in light wind. Moreover, the 𝐵 𝑧 component is also introduced above the Mach disk. On the contrary, a highly dense wind does not exhibit a significant enhancement in the magnetic field ( 𝐵 𝑦 and 𝐵 𝑧 ) above the Mach disk. The dominating toroidal component is notably higher in this case, particularly near the top edges of the lateral streams, which have partially emerged in this case.

## APPENDIX C: JETS VS WINDS: SYNCHROTRON EMISSION AND POLARIZATION

Themapsforthelogarithmic synchrotron flux and polarization for the jets from Paper I and the winds from this study, observed at a viewing angle of 𝜃 𝐼 = 60 ◦ are presented in this section. Specifically, the maps correspond to low resolution, where we employ a Gaussian filter with a FWHM of 720 pc for convolution, and are shown in Figs. C1 and C2. Our focus primarily lies on displaying flux levels within a one-order dynamic range across all figures. However, for comparative purposes, we also include maps for W42 and W44 spanning a twoorders-of-flux range, which are shown in the top and bottom right plots. Regions below the lower limit of the flux color bar are omitted from both emission and polarization images.

This paper has been typeset from a T E X/L A T E X file prepared by the author.

5

20

~20

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-3fb08f3c-6f39-4876-8461-b9b2f218463a.jpg)

3

1

2

3

Y [kpc]

Y [kpc]

Figure B1. Top: Time evolution of the 𝐵 𝑦 component of the magnetic field in the 𝑌 -𝑍 midplane for light wind run of W43 from Table 1. An elongated horizontal structure with high 𝐵 𝑦 values forms near the cocoon head with time. This high-field component gives rise to the upper bright horizontal arc observed in the middle top panel of Fig. 7. Middle and Bottom: 𝑋,𝑌 and 𝑍 -components of the magnetic field in the 𝑌 -𝑍 plane for light (W43\_L1000, top row) and highly dense wind (W43\_L1000HD; bottom row) simulations of W43. 𝐵 𝑦 component of the toroidal field gives rise to a horizontal bright arc in the emission maps (as discussed in Sec. 3.2.1) at 𝜃 𝐼 = 90 ◦ , 𝜙 𝐼 = 0 ◦ image plane.

20

10

~10

20

10

~10

## Synchrotron flux

Figure C1. Same as Fig. 13, but convolved with a Gaussian filter with a FWHM of 720 pc.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-043bda04-50fd-4509-b7bf-fcaa2e1d51a7.jpg)

## Polarization fraction (Tr)

Figure C2. Convolved polarization fraction maps corresponding to the emission maps displayed in Fig. C1. The Gaussian filter used for convolution has a FWHMof 720 pc.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00099v1-figure-ddd4fa73-d9d0-4faa-a528-b7eb02fd89c6.jpg)
Moun Meenakshi Dipanjan Mukherjee Gianluigi Bodo Paola Rossi Chris M. Harrison
A comparative study of radio signatures from winds and jets: Modelling synchrotron emission and polarization
Outflows driven by active galactic nuclei (AGN) are seen in numerous compact
sources; however, it has remained unclear how to distinguish between the
driving mechanisms, such as winds and jets. Therefore, our study aims to offer
observational insights from simulations to aid in this distinction.
Specifically, in this paper, we investigate the evolution of wide-angled,
moderately relativistic, magnetized winds and analyze their non-thermal radio
emission and polarization properties. We find that the evolution of winds
varies depending on factors such as power, density, and opening angle, which in
turn influence their observable characteristics. Additionally, different
viewing angles can lead to varying observations. Furthermore, we note
distinctions in the evolution of winds compared to jets, resulting in
disparities in their observable features. Jets typically exhibit a thin spine
and hotspot(s). Winds manifest broader spines or an "hourglass-shaped" bright
emission in the cocoon, which are capped by bright arcs. Both display high
polarization coinciding with the bright spine and hotspots/arcs, although these
regions are relatively compact and localized in jets when compared to winds. We
emphasize the importance of high resolution, as we demonstrate that emission
features from both jets and winds can become indistinguishable at lower
resolutions. The distribution of polarization is largely unaffected by
resolution, though lower polarization becomes more noticeable when the
resolution is decreased.