CHANG-ES: XVIII---The CHANG-ES Survey and Selected Results
Judith Irwin Ancor Damas-Segovia Marita Krause Arpad Miskolczi Jiangtao Li Yelena Stein Jayanne English Richard Henriksen Rainer Beck Theresa Wiegert Ralf-Juergen Dettmar

 













1
Department of Physics, Engineering Physics & Astronomy, Queen's University, Kingston K7L 3N6, ON, Canada; irwinja@queensu.ca; henriksn@astro.queensu.ca; theresa.wiegert@gmail.com
2Departamento de Astronomía Extragaláctica, Instituto de Astrofísica de Andalucía, Glorieta de la Astronomía sn, 18008, Granada, Spain; adamas@iaa.es
3Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany; mkrause@mpifr-bonn.mpg.de; rbeck@mpifr-bonn.mpg.de
4Faculty of Physics & Astronomy, Astronomical Institute, Ruhr-University Bochum, 44780 Bochum, Germany; miskolczi@astro.rub.de; ystein@astro.ruhr-uni-bochum.de, dettmar@astro.rub.de
5Department of Astronomy, University of Michigan, 311 West Hall, 1085 S. University Ave., Ann Arbor, 48109 MI, USA; jiangtal@umich.edu
6Observatoire astronomique de Strasbourg, Université de Strasbourg, CNRS, UMR 7550, 11 rue de l'Université, F-67000 Strasbourg, France
7Deptartment of Physics & Astronomy, University of Manitoba, Winnipeg, MB, R3T 2N2, Canada; jayanne_english@umanitoba.ca














The CHANG-ES (Continuum Halos in Nearby Galaxies) survey of 35 nearby edge-on galaxies is revealing new and sometimes unexpected and startling results in their radio continuum emission. The observations were in wide bandwidths centred at 1.6 and 6.0 GHz. Unique to this survey is full polarization data showing magnetic field structures in unprecedented detail, resolution and sensitivity for such a large sample.  A wide range of new results are reported here, some never before seen in any galaxy. We see circular polarization and variability in active galactic nuclei (AGNs), in-disk discrete features, disk-halo structures sometimes only seen in polarization, and broad-scale halos with reversing magnetic fields, among others. This paper summarizes some of the CHANG-ES results seen thus far. Released images can be found at <https://www.queensu.ca/changes>
  .









































§ INTRODUCTION TO THE CHANG-ES SURVEY




 §.§ Motivation

CHANG-ES (Continuum Halos in Nearby Galaxies—an EVLA Survey) is a project to  observe 35 edge-on galaxies in the radio continuum with wide bandwidths centred at 1.6 GHz (L-band) and 6.0 GHz (C-band,  using the Karl J. Jansky Very Large Array (hereafter, the VLA, formerly the Expanded Very Large Array or EVLA).  The wide bandwidths, Δ, of the VLA (Table <ref>) have made such a survey feasible since the theoretical signal-to-noise (S/N) improves as 1/√(Δ), all else being equal.  For CHANG-ES, this has meant an order of magnitude improvement, on average, compared to any similar survey.   For example, <cit.> achieved sensitivities ranging from 80 to 130 Jy beam^-1 rms and <cit.> achieved 45 to 90 Jy beam^-1  rms at 5 GHz.  By contrast, typical rms values of the CHANG-ES D-configuration images are ≈ 6Jy beam^-1 at the same frequency (range of values for all configurations and frequencies are summarized in Table <ref>).
Additionally, no previous survey has included polarization. This has opened up the possibility of exploring the magnetic fields in galaxy disks and halos for a well-defined galaxy sample.


lcc
 
Overview of galaxy and observing parameters over frequencies and VLA 
  array configurations.
0pt

Parameter	               Value	          Ref

Galaxy Distance range                      4.4→42  Mpc            <cit.>
	
Galaxy SFR range                           0.02→7.29 M_⊙/yr          <cit.>

Galaxy optical diameter range               3.9→15.8 arcmin      <cit.>

Array configurations (L-band, C-band)      B C D, C D                  <cit.>

ν_0 (L-band, C-band)	       1.58, 6.00 GHz			             <cit.>

Δ ν (L-band, C-band)     512 MHz, 2.0 GHz		             <cit.>

No. of spectral channels  (L-band, C-band)     2048, 1024              <cit.>



Approx. Spatial resolution range     3→60 arcsec      <cit.> 



Typical resolutions in arcsec         <cit.>

(high-res, low-res uv weightings):   

     B configuration L-band      3, 6

     C configuration C-band      3, 6

     C configuration L-band      10, 15    

     D configuration C-band      10, 15   

     D configuration L-band      35, 45    

Approx. rms range       

     (all configurations and bands)      3→ 95 Jy/beam       









Wide bandwidths, when broken up into many individual channels, also provide advantages other than improved S/N: (a) radio frequency interference (RFI) can be more readily identified and excised, (b) in-band spectral index maps can be formed (e.g., Section <ref>), (c) Rotation Measure (RM) analysis can be carried out (Sections <ref> and <ref>), and (d) imaging can make use of the power of multi-frequency synthesis.







CHANG-ES has a variety of goals. We wanted to understand the incidence of radio halos in `normal' spiral galaxies, determine their scale heights and correlate these results with other properties such as the underlying star formation rate (SFR) or SFR `surface density' (SFR per unit area, SFR_ SA). Measurements of outflow speeds can help to determine whether galaxies are more likely to experience winds or `fountain' flow.  Outflows in nearby galaxies could also represent low-energy analogues of the `feedback' that appears to be required in galaxy formation scenarios in order to control the star formation rate <cit.>.  Since outflows can be observed with spectacular resolution in nearby galaxies, the physics of such systems may provide meaningful constraints on galaxy formation in the early universe. A consistent galaxy sample, with complementary data at other wavelengths, can also help us explore pressure balance in the ISM.  For more information on our science goals, see <cit.>. 


As indicated above, unique to the survey is that all polarization products were measured, namely Stokes I, Q, U, and V.  With analysis of linear polarization (P = √(Q^2+U^2)), along with RMs, an in-depth understanding of the magnetic fields in these galaxies can be obtained, and several examples (Section <ref>) have been very illuminating in this regard. Some of our galaxies show lagging halos (e.g., see others <cit.>) either in optical emission lines or HI.  The reason for such lags is not yet understood, but an association with magnetic fields is possible <cit.>.

A survey size that samples a sufficient number of galaxies in a reasonable (but large) total observing time resulted in 35 galaxies in total.  Clearly, rms values that improve upon previous measurements were desired. Since vertical features in galaxies show structures on virtually all scales, we required observations over a range of VLA array configurations. Thus, we have observed at B, C, and D-configurations in L-band and in C and D-configurations in C-band.  This also provides us with matching resolution for spectral index maps which are essential in understanding the origin of the emission (e.g., thermal or non-thermal fractions) as well as the propagation of cosmic ray electrons (CREs). 


Our galaxy selection criteria were (1) an inclination > 75^∘ as given in the Nearby Galaxies Catalog <cit.> in order to easily see a halo,  (2) a declination  > -23^∘ for detectability at the VLA, (3) an optical size between 4 and 15 arcmin for a good match to the array configurations, and (4) flux densities at L-band of at least 23 mJy in order for a detection to be more likely.  Two more galaxies were added that did not meet these criteria but were included  because of other interesting known halo properties (NGC 5775 which was smaller than the lower angular size limit, and NGC 4244 which had a flux density below the minimum cut-off). 

Notice that selection effects should be minor, the main one being the possibility that, by adopting a radio flux density limit, we could inadvertently be choosing galaxies with high SFRs.  However, Table <ref> reveals that our SFR range is quite typical of `normal' spirals. We do not choose our galaxies based on whether or not they have halos (except for the two mentioned above) nor do we specifically choose any active galactic nuclei (AGNs) although several were found to be present once the sample was defined.
Table <ref> provides a summary of our data.  In all, over 405 hours of VLA time were granted for this large project.



 §.§ Techniques and Data Products



Data reduction was carried out using the Common Astronomy Software Applications (CASA) package[<http://casa.nrao.edu>].  This allowed us to take advantage of a number of sophisticated imaging algorithms such as multi-scale multi-frequency synthesis (ms-mfs, <cit.>).  `Multi-scale' assumes that any intensity feature can be represented by spatial scales of varying width prior to convolution with the synthesized beam, as opposed to the traditional assumption that an intensity feature is represented as a series of points.  `Multi-frequency' employs the technique of sampling differing spatial scales as the frequency varies across the band, thereby more thoroughly filling in the uv plane.  In addition, these techniques, as employed in CASA, allow for the fitting of in-band spectral indices across any individual band. In principle, one can solve for curvature, β, in the spectral index as well  (i.e., I_ν ∝ ν^α+β logν but our tests showed that the S/N for any given data set was not high enough to take advantage of the curvature option <cit.>.  Hence, a straight line in log space  (β = 0) was fitted across the band while imaging.  Thus, a solution for the  in-band spectral index, α, was obtained point-by-point, in each galaxy.

At the time of writing, D-configuration images <cit.> have been released at our data release website[<https://www.queensu.ca/changes>], the B-configuration images have also been released <cit.> and the C-configuration <cit.> data release is imminent.  Release data include total intensity images at two different uv weightings,  i.e., `robust 0' and `robust 0 + uv taper', specified as high-res and low-res, respectively, in Table <ref>. Both non-primary-beam-corrected and primary-beam-corrected versions are included. At D-configuration, linear polarization and polarization angle maps, and in-band spectral index maps are also available.  At high resolution, where emission is weaker, we opted to include band-to-band spectral index maps, rather than in-band spectral index maps, i.e.,  C-configuration C-band to B-configuration L-band α maps.

Because the largest angular size detectable at the VLA is 16 arcmin at L-band and 4 arcmin at C-band, it was necessary to supplement the VLA data with single-dish data for the largest galaxies at L-band and all galaxies at C-band (Table <ref>).  This has been accomplished with over 200 h of supplementary Greenbank Telescope (GBT) observations.  Because the GBT wide-band instrumentation was delayed and because data reduction algorithms had to be developed, all GBT observations are completed but the data reduction is currently in progress.  For several galaxies (NGC 891, NGC 4565, and NGC 4631), Effelsberg data have been used for the large-scale total intensity emission.

Complementary to our radio data, we have obtained
Hα images using the Apache Point Observatory 3.5-m telescope <cit.> for 21 galaxies, supplemented with images from the literature for the remaining galaxies.  This data set was obtained to assist with the thermal/non-thermal separation analysis (Section <ref>) and will also soon be publicly released.  We also have WISE enhanced resolution images of all galaxies and XMM-Newton and Chandra data for a large fraction of the sample.











§ CHANG-ES HIGHLIGHTS


In the following sections, we review CHANG-ES results.  Rather than presenting the various papers in chronological order, we present the results topically.







 §.§ AGNs



As indicated in the introduction, we did not select AGNs in the CHANG-ES sample nor was the sample designed specifically for AGN detection.  However, galaxies with AGNs have turned out to provide some of the more interesting results to date, as follows.



  §.§.§ A Remarkable Nearby TDE: NGC 4845



Observations of NGC 4845 resulted in a serendipitous radio detection <cit.> of a Tidal Disruption event (TDE) approximately one year after its hard X-ray outburst <cit.>.  These authors attribute the X-ray variability to a super-Jupiter mass object being devoured by a 10^5 M_⊙ black hole.  CHANG-ES detected the radio emission as well as the variability of the radio emission from its five independent data sets spanning approximately 6 months. Remarkably, this is a very close TDE since the galaxy is in the Virgo Cluster.

The power of in-band spectral fitting paid off for this strong source.  Because we had 3 total intensity data points at L-band, 2 total intensity points at C-band, as well as in-band spectral indices for each one of these data points, it was possible to fit a spectrum to the AGN.  As shown in Figure <ref>a, the AGN is a GHz-peaked source (GPS) which is an order of magnitude closer than any GPS source seen before.  The spectral peak both declines with time as well as shifts to lower frequencies.



The low frequency turnover has been explained by synchrotron self-absorption and the behaviour of the emission modelled by adiabatic expansion of conical or jet-like outflow <cit.>. Our model further predicts the rate of decline of the radio emission with time which is S_ν ∝ t^-5/3 for cone-like geometry (Equation (22) of <cit.>).   Follow-up observations 3 and 4 years later
show a good fit to the predicted t^-5/3 decline (<cit.>, and additional as yet unpublished data).


The fact that such a nearby galaxy harbours a TDE-related outburst means that such systems can be monitored with unprecedented resolution. The relationship between classical AGNs and TDEs may simply be a matter of how much mass is involved in the infall and how sporadically the related outburst occurs.  A second TDE has now also been identified in NGC 660 in the X-ray.













  §.§.§ Circular Polarization


Circular Polarization (CP) is a difficult quantity to measure in extragalactic sources since it is typically very weak (V/I ≈ 0.2%) and is technically challenging.  Nevertheless, the TDE source, NGC 4845, shows a very strong CP signal of ≈ 2% at L-band near the turnover to optical thickness in the spectrum <cit.>. This unusual result prompted us to search through all of the galaxies at high resolution (B configuration L-band), resulting in three more detections (NGC 4388, NGC 660, NGC 3628) and a marginal detection in NGC 3079 <cit.>. In all cases, the detection was right at the core of the galaxy. Weaker signals (less than approx. 0.2%) could potentially be present in the sample but could not have been detected with our technical set-up.  The explanation that is most consistent with the data is Faraday conversion which converts linear to circular polarization along a line of sight (e.g., <cit.>).  The facts that no linear polarization is observed, that the Stokes V spectrum is very steep (α_V ≈ -3), and that CP is observed close to the optically thick limit, are in agreement with this interpretation.  We provide practical expressions for the Stokes V flux density for a variety of physical conditions in galaxy cores (<cit.>, Appendix A).  We also note that the two galaxies with the strongest 2% signals are also those which have the most luminous X-ray cores.

Unlike previous measurements (with the exception of the Galactic Center), CHANG-ES has now revealed nearby galaxies with a CP signal, again providing clearer laboratories than before for studying such phenomena.   Such a study is only in its infancy but has the potential to reveal activity deep into the core where linear polarization is undetectable.  Never-before-seen-related phenomena have now been observed, for example, what we have called the `conversion disk' in NGC 4845 as seen in Figure <ref>b.  Here, Stokes V contours are shown over the total intensity emission at L-band, revealing a disk that is tilted with respect to the total intensity emission.  Presumably, the conversion disk has something to do with the shape of a linearly polarized signal prior to conversion.



  §.§.§ Hidden AGN-Related Activity Revealed


The polarization capability of CHANG-ES has revealed previously hidden AGN activity which had either been masked or was not as clear in total intensity.  A good example is NGC 2992, whose bipolar radio structure is only revealed in polarization <cit.>.  NGC 3079 is another galaxy whose radio lobes are more clearly seen in polarization.

Hidden AGN-related activity can also sometimes be revealed via the spectral index maps.  An example is NGC 3628 which shows bipolar radio emission perpendicular to the major axis via flatter spectral indices in comparison to the surrounding emission <cit.>.

These results show the importance of both polarization and spectral index mapping in disentangling AGN-related emission from emission related to other processes.



  §.§.§ Frequency of AGNs


For the CHANG-ES galaxies, like other spirals, it has historically been a challenge to understand whether outflows might be due to stellar winds and supernovae (SNe) (SF-related activity), or whether the outflow is due to AGNs.  This is because there could very well be both types of activity in galaxy nuclei.  Theoretically, we expect supermassive black holes to be present in all galaxies (e.g., <cit.>), although this does not necessarily imply AGN activity.  Additionally, low-luminosity AGNs (LLAGNs), when present, can be `buried' in other types of emission.

Since the CHANG-ES galaxies are edge-on, obscuration can be severe at other wavelengths so adopting radio criteria for AGNs is useful.  We use (1) the presence of point-like cores, (2) a luminosity that is too high for a collection of SNe, (3) flat or positive spectral indices at the centre, (4) the presence of bipolar or lobe-like structures, (5) time variability, and (6) CP.  Using these criteria where possible, we find an occurrence of 55% of the sample <cit.>. This is a lower limit because weak LLAGNs could still be present in some cases and missed by our criteria.  Higher resolution observations in the future would help to further distinguish AGNs from central starbursts.










 §.§ Edge-on Disks and Thermal/Nonthermal Separation



Since we would like to understand the propagation of CREs as a function of position in our galaxies, we need the non-thermal spectral index, α_NT, as well as the spatially resolved thermal fraction in the CHANG-ES galaxies.  To assist with this, we have obtained Hα maps of most of our galaxies, with the remainder taken from previously published data.  Since both thermal radio emission and Hα emission are proportional to the emission measure, EM, the Hα maps, corrected for dust, can be used to predict the thermal fraction in the radio. The non-thermal emission can then be found as a function of position in the galaxy.

This issue has been addressed by others in the past
(e.g., <cit.>), but not specifically for edge-on galaxies for which dust extinction can be severe.  The first indication of this comes from <cit.> in a plot of L-band luminosity against SFR for the CHANG-ES galaxies, with the SFR determined from the 22 m luminosity.  This plot is essentially the mid-IR radio continuum relation for our galaxies and shows a normalization that differs from face-on systems.  That is, for a given radio continuum luminosity, the SFR is low by a factor of two to three, suggesting that the 22 m flux may itself be suffering from some extinction.  Thus, standard corrections for dust need to be modified for edge-on or very dusty galaxies.

A variety of tests <cit.> have shown that, for resolutions of typically 10 arcsec (≈ 1 kpc at 20 Mpc),

    L( Hα_corr) = L( Hα_obs) + 0.042 ν L_ν(24 m)

where L is the luminosity, corr and obs refer to corrected and observed, respectively, and the 24  m luminosity (approximately the same as at 22m) corrects for dust.  The constant is about 1.4 to 2 × larger than previous authors and takes the higher extinction into account.



It should be pointed out that in the highest resolution images (3 arcsec, or ≈ 300 pc at 20 Mpc), in which HII region complexes in disks start to be resolved, the thermal fraction can be higher than Equation (<ref>) predicts when applied to those discrete HII regions <cit.>.





 §.§ Vertical Scale Heights, Winds, and Some Unexpected Results





An important goal of the CHANG-ES survey is to measure vertical scale heights in our galaxies, an improvement over simply measuring the radio halo extent to some contour level since scale heights will not be strongly dependent on the sensitivity of the image.  Such information (e.g., the ratio of scale heights at two frequencies) can tell us something about the way in which cosmic rays propagate away from the disk, e.g., via diffusion or winds, as outlined via the SPINNAKER code of <cit.>. Careful measurements are required, ensuring that only the highest inclination galaxies and galaxies for which there are no missing large spatial scales are included in the sample. We have therefore measured scale heights in 13 galaxies <cit.> finding sample averages of 1.4 ± 0.7 kpc at L-band and 1.1 ± 0.3 kpc at C-band.

An unexpected result is that, when scale heights are correlated with a variety of parameters (e.g., SFR, B-field strength, etc.), the strongest correlation is actually with galaxy size, i.e., larger galaxies have larger radio halos (see Figure <ref>a).  Further analysis suggests a good correlation with global SFR <cit.>.
The physical reason for this correlation and its relation to star formation is discussed in Krause (2019, this volume).


If we remove the dependence on size by forming a normalized scale height scale, we now see an inverse dependence of the normalized scale height on underlying mass surface density (Figure <ref>b). That is, as the underlying surface mass increases, the normalized scale heights decrease.  This is the first evidence for gravitational deceleration in vertical outflows. Subsequent observations (e.g., NGC 4666 <cit.>) agree with this result.




Once the normalized scale heights are considered, SPINNAKER results show that outflow velocities are consistent with galactic winds with speeds exceeding the escape velocity of the galaxy.  Again, subsequent observations <cit.> agree with this.



A beautiful illustration of global halos can also be seen in
Figure <ref> in which we show the `median edge-on galaxy' in total intensity modified from <cit.>. Thirty galaxies have been scaled to the angular size of the largest galaxy, the median formed, and then superimposed on an optical scaled image of the galaxy, NGC 5775.




 §.§ Magnetic Fields and their Orientation—Reversing Fields







One of the most remarkable results from the CHANG-ES survey is the discovery of reversing magnetic fields in our galaxy disks and halos.  Such a careful and detailed look at halo magnetic fields has never before been carried out and is assisted by the application of the Rotation Measure Synthesis (RMS) technique <cit.>. A number of CHANG-ES galaxies have now had their magnetic fields probed using this technique (e.g., Section <ref> below).



The most stunning result can be seen in NGC 4631 (Figure <ref>) which has a well-known strong halo.  Here, we see systematic reversals  (positive/negative) in the RM (implying systematic reversals in the magnetic field direction) in the halo of the galaxy on kpc scales.  This is the first time that such a regular pattern has been seen in any galaxy halo.  
Follow-up modeling of NGC 4631 by <cit.> and earlier papers <cit.>  have shown that such fields follow from scale-invariant or self-similar solutions of the classical dynamo equations.
The amplitude of these fields is constant on cones in the halo and disc
but the field lines are not confined to these cones. Rather, the field lines  have quite general poloidal components that frequently loop over the spiral arm projections for the asymmetric modes. The axially symmetric mode also has spiral and poloidal projections.  This structure remains a prediction to be tested. However, such field geometry is compatible with the existing observations in CHANG-ES galaxies.

 Sign reversals in edge-on disks have now also been observed <cit.> as well as in other CHANG-ES galaxies (not yet published).




 §.§ Individual Galaxies
 


In this section, we focus on some individual galaxies that have received attention from the CHANG-ES consortium so far.  They were studied in some detail for a variety of reasons, some of which are related to the PhD programs of our student members.



  §.§.§ UGC 10288






UGC 10288 (Figure <ref>) was studied in some detail early in the project <cit.>.  For this galaxy, we combined all five total intensity data sets to obtain a single image at a frequency intermediate between L-band and C-band (4.1 GHz).  The fact that in-band spectral index fitting occurred at every position in the galaxy (Section <ref>) allowed such a combination to be carried out.




Surprisingly, most of the emission from the location of this galaxy is due to a bright background double-lobed radio source (`CHANG-ES A') whose central AGN is situated just above the disk of UGC 10288 and whose southern radio lobe shines right through the disk.  The radio emission from UGC 10288 itself would have fallen below our flux density limit for CHANG-ES, had the background galaxy flux density been subtracted off.  This result warns that care is required in interpreting discrete features that appear to be in the galaxy but may in fact be background sources, especially for such sensitive observations (see also <cit.>).

In spite of the corresponding downwards revision of the SFR for UGC 10288 (≈0.4 M_⊙/yr), vertical discrete radio continuum features can still be seen extending out of the disk of this galaxy.  A preliminary RM analysis of the light from the background source as it passes through the disk and halo region of UGC 10288 gave hints of RM reversals, though not with the stunning regularity shown for NGC 4631 in Figure <ref>.




  §.§.§ NGC 4666


The CHANG-ES radio continuum data of the superwind and starburst galaxy NGC 4666 (Figure <ref> and <cit.>) show a boxy radio halo in total intensity as well as in the X-ray.  The wind is likely driven by supernovae across almost the entire disk.
For the first time, the central point source was observed at radio
wavelengths. In our X-ray data, clear AGN fingerprints were detected, which supports previous findings from X-ray observations that NGC 4666 harbors an AGN.




The high resolution combined C-band data of NGC 4666 (Figure <ref>) show remarkable filamentary vertical structures reaching into the halo-like fingers.
A scale height analysis reveals

a thick disk with a mean scale height of 1.57 ± 0.21  kpc (C-band) and 2.16 ± 0.36 kpc (L-band). The CR transport was modelled with the SPINNAKER code and found to have an advective wind 
of

310 km s^-1. This is different from the CHANG-ES galaxy NGC 4013 <cit.>, where diffusion seems to be the dominant CR transport process.

NGC 4666 is characterized by an extended X-shaped structure with magnetic field vectors reaching far into the halo.

The mean disk magnetic field strength is 12.3 G, which is in good agreement with the field strength of other spiral galaxies <cit.>.
The analysis of the C-band RM-synthesis map along the major axis of NGC 4666 shows indication of one disk magnetic field reversal at a radius of 30" (4 kpc). The field orientation is likely axisymmetric and points inwards and then outwards with increasing radial distance
from the centre (please see <cit.> for more details on the analysis). This is the first time that an indication of
a radial field reversal within the disk is observed in an external
galaxy. As we saw for the magnetic field reversals in the halo of NGC 4631 (Section <ref>), these findings are consistent with mean-field dynamo theory (e.g., <cit.>). 




  §.§.§ NGC 4388


In the case of this Virgo galaxy, the sensitivity provided by the use of the RM-synthesis technique on our C band-C configuration data <cit.> was key to revealing a large number of new faint features in the halo (Figure <ref>). With an exceptionally low rms noise of 2.3 Jy/beam and a resolution of 5.33 arcsec in polarized intensity, we are also able to resolve some interesting features of the disk never seen before.


This is a Seyfert 2 galaxy with many processes occurring at the same time. There are several outflow events detected at different wavelengths (i.e., Hα, X-ray, HI) whose kinematics and morphology are strongly modified by the strong ram-pressure action of the intracluster medium (ICM). The brightest emission in radio polarization in our VLA map is detected at the northern part of the nuclear outflow, where the magnetic field vectors are oriented along the outflow structure. Our data also show for the first time a well defined southern nuclear outflow where the magnetic field vectors are distributed in an S-shaped structure, which is an indication of a connection between the bright radio spots detected above the disk and the nucleus of the galaxy. The halo also shows two polarization filaments to the north and south of the nuclear outflow that extend up to ∼ 5 kpc in the case of the northern filament. These structures seem to also be connected to the nuclear activity due to its symmetry with respect to the centre. It is still unclear whether their origin is AGN- or star formation-related.

NGC 4388 shows two extended blobs of polarized signal towards the southeast and the northwest parts of the halo and their magnetic field vectors are oriented perpendicular to the disk of the galaxy. By computing the half-life time of the CREs in these regions, we can estimate how far into the halo these particles can reach before losing all their energy. We found that in order to reach heights between 2.9 and ∼ 3.7 kpc, they would need to travel at a speed of 270 km s^-1, which is typical for a galactic wind
.




  §.§.§ NGC 3556 and LOFAR Data


Total power results from CHANG-ES observations of NGC 3556 were combined with results from 144 MHz observations carried out using the Low Frequency Array (LOFAR) <cit.>. These data were taken as part of the LOFAR Two-metre Sky Survey (LoTSS) <cit.>. The combined results were used to calculate the spectral index between L-band and 144 MHz and the magnetic field strength in the galaxy using the equipartition approach by <cit.>. This showed that the disk of the galaxy has an average magnetic field strength of 9 G, decreasing to around 5 G towards the halo. The magnetic field strength map is presented in Figure <ref>.



The total power values and the magnetic field strength were used to model the cosmic ray electron propagation using SPINNAKER <cit.> and its interactive wrapper SPINTERACTIVE [<https://github.com/vheesen/Spinnaker>]. Multiple models of diffusive and advective propagation were tested. The best fitting result was an advective transport of the CREs with a linearly increasing wind speed.



  §.§.§ NGC 891 and NGC 4565

The two nearby spiral galaxies, NGC 891 and NGC 4565, differ largely in their detectable halo extent and their star formation rates. The question is how these differences are related to the (advective and/or diffusive) CRE transport in the disk and halo. In C-band, the total power emission of both galaxies has been corrected for missing zero spacing flux with the single-dish 100-m Effelsberg telescope. After subtracting off the thermal emission, maps of the total magnetic field strength and the nonthermal spectral indices as well as vertical nonthermal scale heights at different radii could be determined <cit.>.

These quantities were used to model the CRE transport by SPINNAKER. <cit.> found that the CRE transport in the halo of NGC 891 is probably dominated by advection with an accelerated galactic wind reaching escape velocity at a height of 9 to 17 kpc. This wind is likely to coexist with diffusion-dominated regions. For NGC 4565, however, the halo is diffusion-dominated up to a height of at least 1 kpc with a diffusion coefficient being only weakly energy dependent.

The high-resolution total intensity radio map of NGC 4565 in C-band has a linear resolution of about 200 pc that allows the possibility of separating both sides of the large radio ring (which could also represent tightly wound inner spiral arms) in the disk, as shown in Figure <ref>. This ring-like feature must be very thin with a vertical thickness of only 260 ± 100 pc and horizontal width of 2.2 ± 2.2 kpc.















§ CONCLUSIONS


The CHANG-ES sample of 35 edge-on galaxies has been observed in wide bands centred at 1.6 and 6.0 GHz over a variety of VLA array configurations.  This survey, which includes spatially resolved polarization and spectral index information, is providing new insight into the halo and disk-halo activity in spiral galaxies.

With thermal/non-thermal separation of the emission, measurements of scale-heights, and the sophisticated modeling of CRE outflows, the physical conditions in such regions are being probed in considerably more detail than ever before.  The results of such modeling suggest that most of our galaxies are experiencing outflowing winds.

AGNs have been revealed in the CHANG-ES sample in a variety of ways, including in circular polarization. Bipolar features can be seen in polarization or spectral index structures which otherwise have been masked in total intensity.

The polarization data together with a rotation measure synthesis analysis have now revealed structures that have never before been seen in any galaxy.  These are regularly reversing magnetic fields on kpc scales both in the disk and halo, likely because of magnetic dynamo action.

We summarize CHANG-ES papers to date in Table <ref>. CHANG-ES FITS images are available at the data release web site:
<https://www.queensu.ca/changes>.



lc
 
 List of CHANG-ES papers to date.
 0pt


  
 
   Paper Title          Ref
 

CHANG-ES I: Continuum Halos in Nearby Galaxies: An EVLA Survey: Introduction to the Survey     <cit.>


CHANG-ES II: Continuum Halos in Nearby Galaxies: An EVLA Survey: First Results on NGC 4631     <cit.>


CHANG-ES III: UGC 10288—An Edge-on Galaxy with a Background Double-lobed Radio Source     <cit.>


CHANG-ES IV: Radio continuum emission of 35 edge-on galaxies observed with the 

                             Karl G. Jansky Very Large Array in D-configuration — Data Release 1     <cit.>


CHANG-ES V: Nuclear Outflow in a Virgo Cluster Spiral after a Tidal Disruption Event     <cit.>


CHANG-ES VI: Probing Supernova energy deposition in spiral galaxies through Multiwavelength

                             relationships     <cit.>


CHANG-ES VII: Magnetic Outflows from the Virgo Cluster Galaxy NGC 4388   <cit.>


CHANG-ES VIII: Uncovering Hidden AGN activity in radio polarization     <cit.>


CHANG-ES IX: Radio scale heights and scale lengths of a consistent sample of 13 spiral galaxies

                             seen edge-on and their correlations     <cit.>


CHANG-ES X: Spatially Resolved Separation of Thermal Contribution from Radio Continuum

                             Emission in Edge-on Galaxies     <cit.>


CHANG-ES XI: Circular polarization in the cores of nearby galaxies     <cit.>


CHANG-ES XII: A LOFAR and VLA View of the Edge-on Star Forming Galaxy, NGC 3556     <cit.>


CHANG-ES XIII: Transport processes and the Magnetic Fields of NGC 4666 — Indication of a

                             Reversing Disk Magnetic Field     <cit.>


CHANG-ES XIV: Cosmic-ray Propagation and Magnetic Field Strengths in the Radio Halo of NGC 4631     <cit.>


CHANG-ES XV:  Large-scale Magnetic Field Reversals in the Radio Halo of NGC 4631    <cit.>


CHANG-ES XVI: An In-Depth View of the Cosmic-ray Transport in the Edge-on Spiral Galaxies

                             NGC 891 and NGC 4565     <cit.>


CHANG-ES XVII: H-alpha Imaging of Nearby Edge-on Galaxies, New SFRs, and an Extreme

                             Star Formation Region — Data Release 2    <cit.>


CHANG-ES XVIII: The CHANG-ES Survey and Selected Results     <cit.>


CHANG-ES XIX: The Galaxy NGC 4013 — A diffusion-dominated radio halo with plane-parallel   

                             disk and vertical halo magnetic fields     <cit.>


CHANG-ES XX: High Resolution Radio Continuum Images of Edge-on Galaxies and their AGNs

                             — Data Release 3   <cit.>


CHANG-ES XXI: Radio continuum emission of 35 edge-on galaxies observed with the

                             Karl G. Jansky Very Large Array in C-configuration — Data Release 4   <cit.>
































The first author wishes to thank the Natural Sciences and Engineering Research Council of Canada for a Discovery Grant.




































999






[Beck & Krause (2005)]bec05
Beck, R.; Krause, M.  Astron. Nachrichten 2005, 326, 414–427
[Beck (2015)]bec16
Beck, R.  Astron. Astrophys. Rev. 2016, 24, 4
[Brentjens & de Bruyn (2005)]bre05
Brentjens, M.A.; de Bruyn, A.G.  Astron. Astrophys. 2005, 441, 1217–1228
[Damas-Segovia (2016)]VII
Damas-Segovia, A.; Beck, R.; Vollmer, B.; et al.  Astrophys. J. 2016, 824, 30 
[Heesen et al. (2016)]hee16
Heesen, V.; Dettmar, R.-J.; Krause, M.; et al. Mon. Not. R. Astron. Soc. 2016, 458, 332–353 
[Heesen et al. (2018)]hee18
Heesen, V.; Krause, M.; Beck, R.; et al.  Mon. Not. R. Astron. Soc. 2018, 476, 158–183
[Henriksen & Irwin (2016)]hen16
Henriksen, R.N.; Irwin, J.A.  Mon. Not. R. Astron. Soc. 2016, 458, 4210–4221
[Henriksen et al. (2018a)]hen18a
Henriksen, R.N.; Woodfinden, A.; Irwin, J.A. arXiv 2018, arXiv:1802.07689. 
[Henriksen et al. (2018b)]hen18b
Henriksen, R.N.; Woodfinden, A.; Irwin, J. A.  Mon. Not. R. Astron. Soc. 2018, 476, 635–645
[Hummel et al. (1991)]hum91
Hummel, E.; Beck, R.; Dettmar, R.-J.  Astron. Astrophys. Suppl. Ser. 1991, 87, 309–317.
[Irwin et al. (1999)]irw99
Irwin, J.; English, J.; Sorathia, B.  Astron. J. 1999, 117, 2102–2140
[Irwin et al. (2012a)]I
Irwin, J.; Beck, R.; Benjamin, R.; et al.  Astron. J. 2012, 144, 43 
[Irwin et al. (2012b)]II
Irwin, J.; Beck, R.; Benjamin, R.A.; et al.  Astron. J. 2012, 144, 44 
[Irwin et al. (2013)]III
Irwin, J.; Krause, M.; English, J.; et al.  Astron. J. 2013, 146, 164 
[Irwin et al. (2015)]V
Irwin, J.A.; Henriksen, R.N.; Krause, M.; et al.  Astrophys. J. 2015, 809, 172 
[Irwin et al. (2017)] VIII
  Irwin, J.A.; Schmidt, P.; Damas-Segovia, A.; et al.  Mon. Not. R. Astron. Soc. 2017, 464, 1333–1346  
[Irwin et al. (2018)]XI
Irwin, J.A.; Henriksen, R.N.; Weżgowiec, M.; et al.  Mon. Not. R. Astron. Soc. 2018, 476, 5057–5074 
[Irwin et al. (2019)]XVIII
Irwin, J.; Damas-Segovia, A.; Krause, M.; et al.  Galaxies 2019, this paper (accepted)

[Irwin et al. (2019)]XX
Irwin, J.; Wiegert, T.; Merritt, A.; et al.  2019, in prep. 
[Kennett & Melrose (1998)]ken98 Kennett, M.; Melrose, D.  Publ. Astron. Soc. Aust. 1998, 15, 211–216
[King (2016)]kin16
King, A.  Mon. Not. R. Astron. Soc. 2016, 456, L109–L112
[Krause et al. (2018)]IX
Krause, M.; Irwin, J.; Wiegert, T.; et al.  Astron. Astrophys. 2018, 611, A72 
[Li et al. (2016)]VI
Li, J.-T.; Beck, R.; Dettmar, R.-J.; et al.  Mon. Not. R. Astron. Soc. 2016, 456, 1723–1738 
[Miskolczi et al. (2019)]XII
Miskolczi, A.; Heesen, V.; Horellou, C.; et al.  Astron. Astrophys. 2019, 622, A9 
[Mora-Partiarroyo et al. (2019a)]XIV
Mora-Partiarroyo, S.C.; Krause, M.; Basu, A.; et al.  Astron. Astrophys., 2019, submitted. 
[Mora-Partiarroyo et al. (2019b)]XV
Mora-Partiarroyo, S.C.; Krause, M.; Basu, A.; et al.  Astron. Astrophys. 2019, submitted. 
[Moss et al. (2012)]moss2012
Moss, D.; Stepanov, R.; Arshakian, T.G.; et al.  Astron. Astrophys. 2012, 537, A68
[Nikolajuk & Walter (2013)]nik13
Nikolajuk, M.; Walter, R.  Astron. Astrophys. 2013, 552, A75 
[Perlman et al. (2017)]per17
Perlman, E.S.; Meyer, E.T.; Wang, Q.D.; et al.  Astrophys. J. 2017 842, 126
[Rau & Cornwell (2011)]rau11
Rau, U.; Cornwell, T.J.  Astron. Astrophys. 2011, 532, A71
[Schmidt et al. (2019)]XVI
Schmidt, P.; Krause, M.; Heesen, V.; et al. Astron. Astrophys. 2019, submitted. 
[Shimwell et al. (2019)]LoTSS
Shimwell, T.W.; Tasse, C.; Hardcastle, M.J.; et al.  Astron. Astrophys. 2019, 622, A1


[Silk & Mamon (2012)]sil12Silk, J.; Mamon, G. A.  Res. Astron. Astrophys. 2012, 12, 917–946

[Stein et al. (2019a)]XIII
Stein, Y.; Dettmar, R.-J.; Irwin, J.; et al.  Astron. Astrophys. 2019, 623, A33 

[Stein et al. (2019b)]XIX
Stein, Y.; Dettmar, R.-J.; Weżgowiec, M.; et al.  Astron. Astrophys. 2019, submitted. 

[Tabatabaei et al. (2007)]tab07
Tabatabaei, F.S.; Beck, R.; Krügel, E.; et al. Astron. Astrophys. 2007, 475, 133–142

[Tabatabaei et al. (2017)]tab17
Tabatabaei, F.S.; Schinnerer, E.; Krause, M.; et al.  Astrophys. J. 2017, 836, 185

[Tully (1988)]tul88
Tully, R.B. Nearby Galaxies Catalog; Cambridge University Press: Cambridge, UK; New York, NY, USA, 1988.

[Vargas et al. (2018)]X
Vargas, C.J.; Mora-Partiarroyo, S.C.; Schmidt, P.; et al. Astrophys. J. 2018, 853, 128 

[Vargas et al. (2019)]XVII
Vargas, C.; Walterbos, R.A.M.; Rand, R.J.; et al.  Astrophys. J. 2019, submitted. 



[Walterbos et al. (2019)]XXI
Walterbos, R.; Irwin, J.; Wiegert, T.; et al.  2019, in prep.

[Wiegert et al. (2015)]IV
Wiegert, T.; Irwin, J.; Miskolczi, A.; et al.  Astron. J. 2015, 150, 81 

[Woodfinden et al. (2019)]woo19
Woodfinden, A.; Henriksen, R.N.; Irwin, J.  Mon. Not. R. Astron. Soc. 2019, submitted. 

[Zschaechner et al. (2015)]zsc15
Zschaechner, L.; Rand, R.J.  Astrophys. J. 2015, 808, 153

























