To explore the spatial variations of the regular (mean) magnetic field of the Andromeda galaxy (M31), we use Fourier analysis in azimuthal angle along four rings in the galaxy's plane. The Fourier coefficients give a quantitative measure of strength of the modes, enabling us to compare expectations from mean-field dynamo models of spiral galaxies. Earlier analyses indicated that the axisymmetric magnetic field (azimuthal Fourier mode m = 0) is sufficient to fit the observed polarization angles in a wide range of galactocentric distances. We apply a Bayesian inference approach to new, more sensitive radio continuum data at AA3.59, 6.18, and 11 . 33 cm and the earlier data at A20 . 46 cm to reveal subdominant contributions from the modes m = 1, 2, and 3 along with a dominant axisymmetric mode. Magnetic lines of the axisymmetric mode are close to trailing logarithmic spirals which are significantly more open than the spiral arms detectable in the interstellar dust and neutral hydrogen. The form of the m = 0 mode is consistent with galactic dynamo theory. Both the amplitudes and the pitch angles of the higher azimuthal modes ( m > 1) vary irregularly with r reflecting local variations in the magnetic field structure. The maximum strength of the mean magnetic field of 1 . 8-2 . 7 mu G (for the axisymmetric part of the field) occurs at 10-14 kpc but we find that its strength varies strongly along the azimuth; this variation gives rise to the m = 1 mode. We suggest a procedure of Bayesian inference which is independent of the specific nature of the depolarization and applies when the magneto-ionic layer observable in polarized emission is not symmetric along the line of sight because emission from its far side is completely depolarized.

We present the Very Large Array C-array S-band (2-4 GHz) radio continuum observations of the nearby edge-on spiral galaxy NGC 4565, a target from the Continuum Halos in Nearby Galaxies-an EVLA Survey. We conduct rotation measure (RM) synthesis to probe the magnetic field structure and analyze the vertical radio continuum intensity profiles using the 1D cosmic-ray (CR) transportation models. The radio continuum emission of NGC 4565 is vertically compact, with a vertical-to-radial extent ratio of similar to 1/6. Its vertical profile is optimally described by a two-component Gaussian distribution, yielding a mean Gaussian halo scale height of similar to 3.0 kpc. The magnetic field is weak, predominantly disk parallel, with an equipartition strength of less than or similar to 5 mu G and an RM profile indicative of an axisymmetric spiral structure. Nevertheless, we identify a localized, faint vertical magnetic field component in the northeastern region, hinting at an X-shaped structure that spatially coincides with extraplanar structures detected in H i and soft X-ray emission. The CR transport modeling favors a flux-tube advection scenario, with a slow initial velocity of v0 approximate to 60 km s-1, consistent with a limited energy input from star formation. Therefore, the absence of an extended radio halo can be explained by the low star formation rate, the weak magnetic field, and the inefficient CR transport. The localized X-shaped field may trace a weak, magnetically guided outflow or a tidal perturbation induced by the nearby companion. NGC 4565 is thus a key quiescent benchmark for understanding the physical conditions required to drive large-scale outflows and generate extended radio halos.

The intensities of total, linearly polarized, and unpolarized synchrotron emission are measures for the strengths of total, ordered, and isotropic turbulent fields in the sky plane. Faraday rotation measures (RMs) provide a model of the regular field. The quadratic difference between ordered and regular field strengths yields the strength of the anisotropic turbulent field. - Based on observations of M31 at 3.6 cm, 6.2 cm, and 20.5 cm wavelengths and assuming equipartition between the energy densities of total magnetic fields and total cosmic rays, we measured average equipartition strengths of the magnetic field in the emission torus of M31 of 6.3±.2 for the total, 5.4±0.2 for the isotropic turbulent, and 3.2±0.3 for the ordered field in the sky plane. The average strength of the axisymmetric regular field, Breg, is 2.0±0.5 and remains almost constant between 7 kpc and 12 kpc radius. Quadratic subtraction of the component Breg,perp in the sky plane from the ordered field Bord,perp yields the strength of the anisotropic turbulent field Ban,perp, which is 2.7±0.7 . - The average strength of the regular field is about 40

Resolved studies of the correlation between the radio and far-infrared (FIR) emission from galaxies at different frequencies can unveil the interplay between star formation and the relativistic interstellar medium (ISM). Thanks to the LOFAR LoTSS observations combined with VLA, Herschel, and WISE data, we study the role of cosmic rays and magnetic fields in the radio-FIR correlation on scales of greater than or similar to 200 pc in the nearby galaxy IC 342. The thermal emission traced by the 22 mu m emission, constitutes about 6%, 13%, and 30% of the observed radio emission at 0.14, 1.4, and 4.8 GHz, respectively, in star-forming regions and less in other parts. The nonthermal spectral index becomes flatter at frequencies lower than 1.4 GHz (alpha(n) = -0.51 +/- 0.09, S-nu proportional to nu(alpha)(n)) than between 1.4 and 4.8 GHz (alpha(n) = -1.06 +/- 0.19) on average, and this flattening occurs not only in star-forming regions but also in the diffuse ISM. The radio-FIR correlation holds at all radio frequencies; however, it is tighter at higher radio frequencies. A multi-scale analysis shows that this correlation cannot be maintained on small scales due to diffusion of cosmic ray electrons (CREs). The correlation breaks at a larger scale (similar or equal to 320 pc) at 0.14 GHz than at 1.4 GHz (similar or equal to 200 pc), indicating that the CREs traced at lower frequencies have diffused a longer path in the ISM. We find that the energy index of CREs becomes flatter in star-forming regions, in agreement with previous studies. Cooling of CREs due to the magnetic field is evident globally only after compensating for the effect of star formation activity that both accelerates CREs and amplifies magnetic fields. Compared with other nearby galaxies, we show that the smallest scale of the radio-FIR correlation is proportional to the propagation length of the CREs on which the ordered magnetic field has an important effect.

The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of h _I radio ∼ h _PI radio = 0.23 ± 0.03 kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane ( h _I FIR ∼ h _PI FIR = 0.11 ± 0.01 kpc). Radio emission traces magnetic fields at scales of ≳300 pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution (<300 pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism.

We report the discovery of the 10 kpc scale radio lobes in the Sombrero galaxy (NGC 4594), using data from the Continuum Halos in Nearby Galaxies—an Expanded Very Large Array Survey project. We further examine the balance between the magnetic pressure inside the lobes and the thermal pressure of the ambient hot gas. At the radii r of ∼(1–10) kpc, the magnetic pressure inside the lobes and the thermal pressure of the ambient hot gas are generally in balance. This implies that the jets could expand into the surroundings to at least r ∼ 10 kpc. The feedback from the active galactic nucleus jet responsible for the large-scale lobes may help to explain the unusually high X-ray luminosity of this massive quiescent isolated disk galaxy, although more theoretical work is needed to further examine this possibility.

CHANG-ES (Continuum Halos in Nearby Galaxies—an EVLA Survey) is an ambitious project to target 35 nearby disk galaxies that are edge-on to the line of sight. The orientation permits both the disk and halo regions to be studied. The observations were initially at 1.5 GHz (L-band) and 6.0 GHz (C-band) in a variety of VLA array configurations, and in all four Stokes parameters, which allowed for spatially resolved images in total intensity plus polarization. The inclusion of polarization is unique to an edge-on galaxy survey and reveals the galaxies’ halo magnetic fields. This paper will summarize the results to date, some of which are new phenomena, never seen prior to CHANG-ES. For example, we see that ‘X-type’ fields, as well as rotation measure reversals, are common features of spiral galaxies. Further observations at 3.0 GHz (S-band) as well as future scientific opportunities will also be described.

Context. The existence of magnetic fields in the circumgalactic medium (CGM) is largely unconstrained. Their detection is important as magnetic fields can have a significant impact on the evolution of the CGM and, in turn, the fields can serve as tracers for dynamical processes in the CGM. Aims. With Faraday rotation of polarised background sources, we aim to detect a possible excess of the rotation measure in the surrounding area of nearby galaxies. Methods. We use 2,461 residual rotation measures (RRMs) observed with the LOw Frequency ARray (LOFAR), where the foreground contribution from the Milky Way is subtracted. The RRMs are then studied around a subset of 183 nearby galaxies that was selected by apparent $B$-band magnitude. Results. We find that, in general, the RRMs show no significant excess for small impact parameters (i.e. the perpendicular distance to the line of sight). However, if we only consider galaxies at higher inclination angles and sight lines that pass close to the minor axis of the galaxies, we find significant excess at impact parameters of less than 100 kpc. The excess in |RRM| is 3.7 $\rm rad\,m^{-2}$ with an uncertainty between $\pm 0.9~\rm rad\,m^{-2}$ and $\pm 1.3~\rm rad\,m^{-2}$ depending on the statistical properties of the background (2.8$\sigma$-4.1$\sigma$). With electron densities of ~$10^{-4}~\rm cm^{-3}$ this suggests magnetic field strengths of a few tenths of a micro Gauss. Conclusions. Our results suggest a slow decrease of the magnetic field strength with distance from the galactic disc such as expected if the CGM is magnetised by galactic winds and outflows.

The extraplanar diffuse ionized gas (eDIG) represents the cool/warm ionized gas reservoir around galaxies. We present spatial analysis of the H alpha images of 22 nearby edge-on spiral galaxies taken with the Apache Point Observatory 3.5-m telescope (eDIG-CHANGES). We conduct an exponential fit to the vertical H alpha intensity profiles of the galaxies, of which 16 can be decomposed into thin + thick disk components. The median value of the H alpha scale height of the thick disk is (h(H alpha)) =1.13 +/- 0.14 kpc. We further examine the dependence of h(H alpha) on the stellar mass, SFR, and SFR surface density (SFRSD) of the galaxies. We find a tight sublinear correlation between h(H alpha) and SFR, expressed in h(H alpha) proportional to SFR alpha, where asymptotic to 0.29. Moreover, the offset of individual galaxies from the best-fit SFR-h(H alpha) relation, expressed in h(H alpha)/SFR alpha, shows significant anti-correlation with SFRSD. We further compare the vertical extension of the eDIG to multi-wavelength measurements of other CGM phases. We find the eDIG slightly more extended than the neutral gas. This indicates the existence of some extended ionizing sources, in addition to the leaking photons from the disk star formation regions. Most galaxies have an X-ray scale height smaller than H alpha, suggesting the majority of the X-ray photons are actually from the thick disk instead of the extended CGM. h(H alpha) is comparable to the L-band radio continuum scale height. This indicates that the thermal and non-thermal electrons have similar spatial distributions, a natural result if both are transported outwards by a galactic wind. This further indicates the thermal gas, cosmic rays, and magnetic field may be close to energy equipartition.

Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multiphase ISM of the closest merger of two spiral galaxies, the Antennae galaxies. We compare the inferred B-fields using 154 μm thermal dust and 11 cm radio synchrotron emission polarimetric observations. We find that the 154 μm B-fields are more ordered across the Antennae galaxies than the 11 cm B-fields. The turbulent-to-ordered 154 μm B-field increases at the galaxy cores and star-forming regions. The relic spiral arm has an ordered spiral 154 μm B-field, while the 11 cm B-field is radial. The 154 μm B-field may be dominated by turbulent dynamos with high 12CO(1–0) velocity dispersion driven by star-forming regions, while the 11 cm B-field is cospatial with high H i velocity dispersion driven by galaxy interaction. This result shows the dissociation between the warm gas mainly disturbed by the merger, and the dense gas still following the dynamics of the relic spiral arm. We find a ∼8.9 kpc scale ordered B-field connecting the two galaxies. The base of the tidal tail is cospatial with the H i and 12CO(1–0) emission and has compressed and/or sheared 154 μm and 11 cm B-fields driven by the merger. We suggest that amplified B-fields, with respect to the rest of the system and other spiral galaxies, may be supporting the gas flow between both galaxies and the tidal tail.

Recent detection of polarized thermal emission from dust grains in a high-redshift, rapidly star-forming galaxy can give us an insight into the formation and evolution of magnetic fields in large-scale structures of the early Universe.

A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg(2) of sky (corresponding to a polarized source areal number density of similar to 0.43 deg(-2)). The linear polarization and RM properties were derived using RM synthesis from the Stokes Q and U channel images at an angular resolution of 20 arcsec across a frequency range of 120 to 168 MHz with a channel bandwidth of 97.6 kHz. The fraction of total intensity sources (>1 mJy beam(-1)) found to be polarized was similar to 0.2 per cent. The median detection threshold was 0.6 mJy beam(-1) (8 sigma(QU)), with a median RM uncertainty of 0.06 rad m(-2) (although a systematic uncertainty of up to 0.3 rad m(-2) is possible, after the ionosphere RM correction). The median degree of polarization of the detected sources is 1.8 per cent, with a range of 0.05 per cent to 31 per cent. Comparisons with cm-wavelength RMs indicate minimal amounts of Faraday complexity in the LoTSS detections, making them ideal sources for RM Grid studies. Host galaxy identifications were obtained for 88 per cent of the sources, along with redshifts for 79 per cent (both photometric and spectroscopic), with the median redshift being 0.6. The focus of the current catalogue was on reliability rather than completeness, and we expect future versions of the LoTSS RM Grid to have a higher areal number density. In addition, 25 pulsars were identified, mainly through their high degrees of linear polarization.

Galactic winds play a key role in regulating the evolution of galaxies over cosmic time. In recent years, the role of cosmic rays (CR) in the formation of the galactic wind has increasingly gained attention. Therefore, we use radio continuum data to analyse the cosmic ray transport in edge-on galaxies. Data from the LOFAR Two-metre Sky Survey (LoTSS) data release 2 at 144 MHz (HBA) and reprocessed VLA data at 1.6 GHz (L-band) from the Continuum Halos in Nearby Galaxies - an EVLA Survey (CHANG-ES) enable us to increase the extent of the analysed radio continuum profile significantly (up to a factor of 2) compared to previous studies. We compute thermal emission maps using a mixture approach of H-alpha and near infrared data, which is then subtracted to yield radio synchrotron emission maps. Then we compile non-thermal spectral index maps and compute intensity profiles using a box integration approach. Lastly, we perform 1D cosmic ray transport modelling. The non-thermal spectral index maps show evidence that the LoTSS maps are affected by thermal absorption, in star forming regions. The scale height analysis reveals that most of the galaxies are equally well fitted with an one-component instead of a twocomponent exponential profile. We find a bi-modality within our sample. While NGC 3432 and NGC 4013 have similar scale heights in the L-band and HBA, the low-frequency scale heights of NGC 891, NGC 4157, and NGC 4631 exceed their high-frequency counterpart significantly. The 1D CR transport modelling shows agreement of the predicted magnetic field strength and the magnetic field strength estimates of equipartition measurements. Additionally we find an increasing difference of wind velocities (with increasing height over the galactic disk) between central and outer regions of the analysed galaxies.

Context. The details of cosmic-ray transport have a strong impact on galaxy evolution. The peak of the cosmic-ray energy distribution is observable in the radio continuum using the electrons as proxy. Aims. We measure the length that the cosmic-ray electrons (CRE) are transported during their lifetime in the nearby galaxy M 51 across one order of magnitude in cosmic-ray energy (approximately 1-10 GeV). To this end we use new ultra-low frequency observations from the LOw Frequency ARay (LOFAR) at 54 MHz and ancillary data between 144 and 8350 MHz. Methods. As the the CRE originate from supernova remnants, the radio maps are smoothed in comparison to the distribution of the star formation. By convolving the map of the star-formation rate (SFR) surface density with a Gaussian kernel, we can linearise the radio-SFR relation. The best-fitting convolution kernel is then our estimate of the CRE transport length. Results. We find that the CRE transport length increases at low frequencies, as expected since the CRE have longer lifetimes. The CRE transport length is $l_{\rm CRE} = \sqrt{4Dt_{\rm syn}}$, where $D$ is the isotropic diffusion coefficient and $t_{\rm syn}$ is the CRE lifetime as given by synchrotron and inverse Compton losses. We find that the data can be well fitted by diffusion, where $D=(2.14\pm 0.13) \times 10^{28}~\rm cm^2\,s^{-1}$. With $D\propto E^{0.001\pm 0.185}$, the diffusion coefficient is independent of the CRE energy $E$ in the range considered. Conclusions. Our results suggest that the transport of GeV-cosmic ray electrons in the star-forming discs of galaxies is governed by energy-independent diffusion.

The relation between the resolved star formation rate (SFR) per unit area and the nonthermal radio continuum emission is studied in 21 Virgo cluster galaxies and the two nearby spiral galaxies, NGC 6946 and M 51. For the interpretation and understanding of our results, we used a 3D model where star formation, 2D cosmic-ray (CR) propagation, and the physics of synchrotron emission are included. Based on the linear correlation between the SFR per unit area and the synchrotron emission and its scatter, radio-bright and radio-dim regions can be robustly defined for our sample of spiral galaxies. We identified CR diffusion or streaming as the physical causes of radio-bright regions of unperturbed symmetric spiral galaxies as NGC 6946. The enhanced magnetic field in the region of interstellar medium (ISM) compression via ram pressure is responsible for the southwestern radio-bright region in NGC 4501. We identified the probable causes of radio-bright regions in several galaxies as CR transport, via either gravitational tides (M 51) or galactic winds (NGC 4532) or ram pressure stripping (NGC 4330 and NGC 4522). Three galaxies are overall radio dim: NGC 4298, NGC 4535, and NGC 4567. Based on our model of synchrotron-emitting disks, we suggest that the overall radio-dim galaxies have a significantly lower magnetic field than expected by equipartition between the magnetic and turbulent energy densities. We suggest that this is linked to the difference between the timescales of the variation in the SFR and the small-scale dynamo. In NGC 4535, shear motions increase the total magnetic field strength via the induction equation, which leads to enhanced synchrotron emission with respect to the SFR in an otherwise radio-dim galactic disk. Radio-bright regions frequently coincide with asymmetric ridges of polarized radio continuum emission, and we found a clear albeit moderate correlation between the polarized radio continuum emission and the radio/SFR ratio. When compression or shear motions of the ISM are present in the galactic disk, the radio-bright regions are linked to the commonly observed asymmetric ridges of polarized radio continuum emission and represent a useful tool for the interaction diagnostics. The magnetic field is enhanced (as observed in NGC 4535 and NGC 4501) and ordered by these ISM compression and shear motions. Whereas the enhancement of the magnetic field is rather modest and does not significantly influence the radio-SFR correlation, the main effect of ISM compression and shear motions is the ordering of the magnetic field, which significantly affects the CR transport. Cosmic-ray energy losses and transport also affect the spectral index, which we measured between 4.85 and 1.4 GHz. The influence of CR losses and transport on the spectral index distribution with respect to the synchrotron/SFR ratio is discussed with the help of model calculations. Based on our results, we propose a scenario for the interplay between star formation, CR electrons, and magnetic fields in spiral galaxies.

We present the analysis of the magnetic field (B-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey of extragALactic magnetiSm with SOFIA of 14 nearby (<20 Mpc) galaxies with resolved (5 '' -18 ''; 90 pc-1 kpc) imaging polarimetric observations using SOFIA/HAWC+ from 53 to 214 mu m. We compute the magnetic pitch-angle (psi(B)) profiles as a function of the galactocentric radius. We introduce a new magnetic alignment parameter (.) to estimate the disordered-to-ordered ratio of spiral B-fields. We find FIR and radio wavelengths to not generally trace the same B-field morphology in galaxies. The psi(B) profiles tend to be more ordered across all galactocentric radii in radio (zeta(6cm) = 0.93 +/- 0.03) than in FIR (zeta(154 mu m) = 0.84 +/- 0.14). For spiral galaxies, FIR B-fields are 2%-75% more turbulent than the radio B-fields. For starburst galaxies, we find that FIR polarization is a better tracer of the B-fields along the galactic outflows than radio polarization. Our results suggest that the B-fields associated with dense, dusty, turbulent star-forming regions (those traced at FIR) are less ordered than warmer, less dense regions (those traced at radio) of the interstellar medium. The FIR B-fields seem to be more sensitive to the activity of the star-forming regions and molecular clouds within a vertical height of a few hundred parsecs in the disk of spiral galaxies than the radio B-fields.

ABSTRACT Studying the interplay between massive star formation and the interstellar medium (ISM) is paramount to understand the evolution of galaxies. Radio continuum (RC) emission serves as an extinction-free tracer of both massive star formation and the energetic components of the ISM. We present a multiband RC survey of the Local Group galaxy M 33 down to ≃30 pc linear resolution observed with the Karl G. Jansky Very Large Array (VLA). We calibrate the star formation rate surface density and investigate the impact of diffuse emission on this calibration using a structural decomposition. Separating the thermal and non-thermal emission components, the correlation between different phases of the ISM, and the impact of massive star formation are being investigated. Radio sources with sizes ≲200 pc constitute about 36 per cent (46 per cent) of the total RC emission at 1.5 GHz (6.3 GHz) in the inner 18 × 18 arcmin2 (or 4 kpc × 4 kpc) disc of M 33. The non-thermal spectral index becomes flatter with increasing star formation rate surface density, indicating the escape of cosmic ray electrons from their birth places. The magnetic field strength also increases with star formation rate following a bi-modal relation, indicating that the small-scale turbulent dynamo acts more efficiently at higher luminosities and star formation rates. Although the correlations are tighter in star-forming regions, the non-thermal emission is also correlated with the more quiescent molecular gas in the ISM. An almost linear molecular star formation law exists in M 33 when excluding diffuse structures. Massive star formation amplifies the magnetic field and increases the number of high-energy cosmic ray electrons, which can help the onset of winds and outflows.

We present the HI distribution of galaxies from the Continuum Halos in Nearby Galaxies - an EVLA Survey (CHANG-ES). Though the observational mode was not optimized for detecting HI, we successfully produce HI cubes for 19 galaxies. The moment-0 maps from this work are available on CHANG-ES data release website, i.e., https://www.queensu.ca/changes. Our sample is dominated by star-forming, HI-rich galaxies at distances from 6.27 to 34.1 Mpc. HI interferometric images on two of these galaxies (NGC 5792 and UGC 10288) are presented here for the first time, while 12 of our remaining sample galaxies now have better HI spatial resolutions and/or sensitivities of intensity maps than those in existing publications. We characterize the average scale heights of the HI distributions for a subset of most inclined galaxies (inclination > 80 deg), and compare them to the radio continuum intensity scale heights, which have been derived in a similar way. The two types of scale heights are well correlated, with similar dependence on disk radial extension and star formation rate surface density but different dependence on mass surface density. This result indicates that the vertical distribution of the two components may be governed by similar fundamental physics but with subtle differences.

We study the vertical distribution of the highly inclined galaxies from the Continuum Halos in Nearby Galaxies—an EVLA Survey (CHANG-ES). We explore the feasibility of photometrically deriving the H I disk scale heights from the moment-0 images of the relatively edge-on galaxies with inclination>80°, by quantifying the systematic broadening effects and thus deriving correction equations for direct measurements. The corrected H I disk scale heights of the relatively edge-on galaxies from the CHANG-ES sample show trends consistent with the quasi-equilibrium model of the vertical structure of gas disks. The procedure provides a convenient way to derive the scale heights and can easily be applied to statistical samples in the future.

ABSTRACT Active galactic nucleus (AGN) bubbles could play an important role in accelerating high-energy cosmic rays (CRs) and galactic feedback. Only in nearby galaxies could we have high enough angular resolution in multiwavelengths to study the sub-kpc environment of the AGN, where the bubbles are produced and strongly interact with the surrounding interstellar medium. In this paper, we present the latest Chandra observations of the Virgo cluster galaxy NGC 4438, which hosts multiscale bubbles detected in various bands. The galaxy also has low current star formation activity, so these bubbles are evidently produced by the AGN rather than a starburst. We present spatially resolved spectral analysis of the Chandra data of the ∼3 arcsec × 5 arcsec (${\sim} 200{\rm ~pc}\times 350\rm ~pc$) nuclear bubble of NGC 4438. The power-law tail in the X-ray spectra can be most naturally explained as synchrotron emission from high-energy CR leptons. The hot gas temperature increases, while the overall contribution of the non-thermal X-ray emission decreases with the vertical distance from the galactic plane. We calculate the synchrotron cooling time-scale of the CR leptons responsible for the non-thermal hard X-ray emission to be only a few tens to a few hundreds of years. The thermal pressure of the hot gas is about three times the magnetic pressure, but the current data cannot rule out the possibility that they are still in pressure balance. The spatially resolved spectroscopy presented in this paper may have important constraints on how the AGN accelerates CRs and drives outflows. We also discover a transient X-ray source only ∼5 arcsec from the nucleus of NGC 4438. The source was not detected in 2002 and 2008, but became quite X-ray bright in 2020 March, with an average 0.5–7 keV luminosity of ${\sim} 10^{39}\rm ~erg~s^{-1}$.