Search

We have carried out spatially resolved thermal/nonthermal separation on two edge-on galaxies, NGC~3044 and NGC~5775, using only radio data. Narrow-band imaging within a frequency band that is almost contiguous from 1.25 to 7.02 GHz (L-band, S-band and C-band) has allowed us to fit spectra and construct thermal, nonthermal, and nonthermal spectral index maps. This method does not require any ancillary H$\alpha$ and infrared data, or reliance on dust corrections that are challenging in edge-on galaxies. For NGC~3044, at 15 arcsec resolution, we find a median thermal fraction of $\sim\, 13$\% with an estimated uncertainty in this fraction of $\sim\, 50$\% at 4.13 GHz. This compares well with the H$\alpha$ mixture method results. We uncovered evidence for a vertical outflow feature reaching at least $z\,\sim\,3.5$ kpc in projection above the plane, reminiscent of M82's starburst wind. For the higher SFR galaxy, NGC~5775 at 12 arcsec resolution, we find a median thermal fraction of 44\% at 4.13 GHz with an estimated error on this fraction of 17\%. Both galaxies show a change of slope (flattening) in L-band. These results suggest that a radio-only method for separating thermal from nonthermal emission is not only feasible, but able to reveal new features that might otherwise be obscured in edge-on disks., Comment: Accepted to AJ on July 18, 2024

Aspera is a NASA Astrophysics Pioneers SmallSat mission designed to study diffuse OVI emission from the warm-hot phase gas in the halos of nearby galaxies. Its payload consists of two identical Rowland Circle-type long-slit spectrographs, sharing a single MicroChannel plate detector. Each spectrograph channel consists of an off-axis parabola primary mirror and a toroidal diffraction grating optimized for the 1013-1057 Angstroms bandpass. Despite the simple configuration, the optical alignment/integration process for Aspera is challenging due to tight optical alignment tolerances, driven by the compact form factor, and the contamination sensitivity of the Far-Ultraviolet optics and detectors. In this paper, we discuss implementing a novel multi-phase approach to meet these requirements using state-of-the-art optical metrology tools. For coarsely positioning the optics we use a blue-laser 3D scanner while the fine alignment is done with a Zygo interferometer and a custom computer-generated hologram. The detector focus requires iterative in-vacuum alignment using a Vacuum UV collimator. The alignment is done in a controlled cleanroom facility at the University of Arizona., Comment: Manuscript submitted for Proceedings of Astronomical Telescopes + Instrumentation 2024, Paper no. 13093-9

We report the discovery of the 10 kilo-parsec (kpc) scale radio lobes in the Sombrero galaxy (NGC 4594), using data from the Continuum Halos in Nearby Galaxies - an Expanded Very Large Array (VLA) Survey (CHANG-ES) 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 at least to r ~ 10 kpc. The feedback from the active galactic nucleus (AGN) 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., Comment: 8 pages, 3 figures, 2 tables, accepted for publication in ApJ

In 2013, the IceCube collaboration announced the detection of a diffuse high-energy astrophysical neutrino flux. The origin of this flux is still largely unknown. The most significant individual source is the close-by Seyfert galaxy NGC 1068 at 4.2-sigma level with a soft spectral index. To identify sources based on their counterpart, IceCube releases realtime alerts corresponding to neutrinos with a high probability of astrophysical origin. We report here the spatial coincidence of two neutrino alerts, IC220424A and IC230416A, with the Seyfert galaxy NGC 7469 at a distance of 70 Mpc. We evaluate, a-posteriori, the chance probability of such a coincidence and discuss this source as a neutrino emitter based on its multi-wavelength properties and in comparison to NGC 1068. To calculate the chance coincidence considering neutrino emission from a specific source population, we perform a Goodness-of-Fit test with a test statistic derived from a likelihood ratio that includes the neutrino angular uncertainty and the source distance. We apply this test first to a catalog of AGN sources and second to a catalog of Seyfert galaxies only. Our a-posteriori evaluation excludes the chance coincidence of the two neutrinos with the Seyfert galaxy NGC 7469 at 3.3-sigma level. To be compatible with non-detections of TeV neutrinos, the source would need to have a hard spectral index., Comment: 9 pages, 6 figures, submitted to Astronomy & Astrophysics

Aims. We analyze a MUSE optical integral field spectrum of the star-forming edge-on galaxy IC 1553 in order to study its extraplanar diffuse ionized gas (eDIG) and the processes shaping its disk-halo interface. Methods. We extracted the optical emission line properties from the integral field spectrum and generated the commonly used emission line diagnostic diagrams in order to analyze the ionization conditions and the distribution of the eDIG. Furthermore, we performed gravitational potential fitting to investigate the kinematics of a suspected galactic outflow. Results. We find that the eDIG scale height has a maximum value of approximately 1.0 kpc and decreases roughly linearly with the radial distance from the galactic center in projection. The ionization state of the eDIG is not consistent with a pure photoionization scenario and instead requires a significant contribution from shock ionization. This, in addition to the gas kinematics, strongly suggests the presence of a galactic scale outflow, the origin of which lies at least 1.4 kpc away from the galactic center. The inferred shock velocity in the eDIG of approximately 225 km s-1 is comparable to the escape velocity estimated from our potential modelling. The asymmetric distribution of currently star-forming clusters produces a range of different ionization conditions in the eDIG. As a result, the vertical emission line profiles vary quantitatively and qualitatively along the major axis of the galaxy. This analysis illustrates that it is crucial in studies of the eDIG to use observations that take the spatial and kinematical distributions into account, such as those done with integral field units, to form an accurate picture of the relevant physical properties., Comment: Accepted for publication by Astronomy & Astrophysics on August 1, 2023

We present the serendipitous discovery of (1) a large double radio relic associated with the galaxy cluster PSZ2 G277.93+12.34 and (2) a new odd radio circle, ORC J1027-4422, both found in the same deep MeerKAT 1.3 GHz wide-band radio continuum image. The angular separation of the two arc-shaped cluster relics is ~16 arcmin or ~2.6 Mpc for a cluster redshift of z ~ 0.158. The thin southern relic, which shows several ridges/shocks including one possibly moving inwards, has a linear extent of ~1.64 Mpc. In contrast, the northern relic is about twice as wide, twice as bright, but only has a largest linear size of ~0.66 Mpc. Complementary SRG/eROSITA X-ray images reveal extended emission from hot intracluster gas between the two relics and around the narrow-angle tail (NAT) radio galaxy PMN J1033-4335 (z ~ 0.153) located just east of the northern relic. The radio morphologies of the NAT galaxy and the northern relic, which are also detected with the Australian Square Kilometer Array Pathfinder (ASKAP) at 888 MHz, suggest both are moving in the same outward direction. The discovery of ORC J1027-4422 in a different part of the same MeerKAT image makes it the 4th known single ORC. It has a diameter of ~90 arcsec corresponding to 400 kpc at a tentative redshift of z ~ 0.3 and remains undetected in X-ray emission. Supported by simulations, we discuss similarities between outward moving galaxy and cluster merger shocks as the formation mechanisms for ORCs and radio relics, respectively., Comment: 16 pages, 11 figures, MNRAS, in press

Interstellar scattering (ISS) of radio pulsar emission can be used as a probe of the ionised interstellar medium (IISM) and causes corruptions in pulsar timing experiments. Two types of ISS phenomena (intensity scintillation and pulse broadening) are caused by electron density fluctuations on small scales (< 0.01 AU). Theory predicts that these are related, and both have been widely employed to study the properties of the IISM. Larger scales ($\sim$1-100\,AU) cause measurable changes in dispersion and these can be correlated with ISS observations to estimate the fluctuation spectrum over a very wide scale range. IISM measurements can often be modeled by a homogeneous power-law spatial spectrum of electron density with the Kolmogorov ($-11/3$) spectral exponent. Here we aim to test the validity of using the Kolmogorov exponent with PSR~J0826+2637. We do so using observations of intensity scintillation, pulse broadening and dispersion variations across a wide fractional bandwidth (20 -- 180\,MHz). We present that the frequency dependence of the intensity scintillation in the high frequency band matches the expectations of a Kolmogorov spectral exponent but the pulse broadening in the low frequency band does not change as rapidly as predicted with this assumption. We show that this behavior is due to an inhomogeneity in the scattering region, specifically that the scattering is dominated by a region of transverse size $\sim$40\,AU. The power spectrum of the electron density, however, maintains the Kolmogorov spectral exponent from spatial scales of 5$\times10^{-6}$\,AU to $\sim$100\,AU., Comment: Accepted for publication in MNRAS, typo fixed

We consider the dynamics of and emission from growing superbubbles in a stratified interstellar gaseous disc driven by energy release from supernovae explosions in stellar clusters with {masses $M_{cl}= 10^5-1.6\times 10^6~M_\odot$}. Supernovae are spread randomly within a sphere of $r_c=60$ pc, and inject energy episodically with a specific rate $1/130~M_\odot^{-1}$ proportional to the star formation rate (SFR) in the cluster. Models are run for several values of SFR in the range $0.01$ to $0.1~M_\odot$ yr$^{-1}$, with the corresponding average surface energy input rate $\sim 0.04-0.4$ erg cm$^{-2}$ s$^{-1}$. We find that the discrete energy injection by isolated SNe are more efficient in blowing superbubbles: asymptotically they reach heights of up to 3 to 16 kpc for $M_{cl}=10^5-1.6\times 10^5~M_\odot$, correspondingly, and stay filled with a hot and dilute plasma for at least 30 Myr. During this time they emit X-ray, H$\alpha$ and dust infrared emission. X-ray liminosities $L_X\propto {\rm SFR}^{3/5}$ that we derive here are consistent with observations in star-forming galaxies. Even though dust particles of small sizes $a\leq 0.03~\mu$m are sputtered in the interior of bubbles, larger grains still contribute considerably ensuring the bubble luminosity $L_{\rm IR}/{\rm SFR}\sim 5\times 10^7 L_\odot M_\odot^{-1} ~{\rm yr}$. It is shown that the origin of the North Polar Spur in the Milky Way can be connected with activity of a cluster with the stellar mass of $\sim 10^5~M_\odot$ and the ${\rm SFR}\sim 0.1~M_\odot$ yr$^{-1}$ some 25--30 Myr ago. Extended luminous haloes observed in edge-on galaxies (NGC 891 as an example) can be maintained by disc spread stellar clusters of smaller masses $M_\ast \simlt 10^5~M_\odot$., Comment: 13 pages, 18 figures, accepted to MNRAS

AT 2000ch is a highly variable massive star and supernova imposter in NGC 3432 first detected in 2000. It is similar and often compared to SN 2009ip, and it is therefore expected to undergo a core-collapse supernova -- a SN imposter of similar brightness -- in the near future. We characterize the long-term variability of AT 2000ch in the radio and optical regimes with archival data reaching back to the year 1984. We use these newly reduced observations in addition to observations in the literature to restrict the mass-loss rates of AT 2000ch at multiple epochs based on different approaches, and to infer the general properties of its circumstellar nebula with respect to the detected radio brightness. We extend the known optical light curve of AT 2000ch up to the beginning of 2022 by performing point spread function photometry on archival data from the Palomar Transient Factory and the Zwicky Transient Facility. We reduced archival radio continuum observations obtained with the Very Large Array using standard calibration and imaging methods and complemented these with pre-reduced \changes observations as well as observations obtained with the Westerbork Synthesis Radio Telescope and LOw Frequency ARray. For the analysis of AT 2000ch, we consider the optical light curve and color evolution, its radio continuum brightness at different frequencies and times, and the corresponding spectral indices. We estimated mass-loss rates and optical depths based on radio continuum brightnesses and Ha fluxes. We report two newly detected outbursts of AT 2000ch similar to those found in the 2000s and 13 re-brightening events, of which at least four are not conclusively detected because of insufficient sampling of the light curve. The dates of all outbursts and significant, well-sampled re-brightening events are consistent with a period of $\sim 201 \pm 12\,$days over a total time-span of two decades. Such a behavior has never been found for any SN imposter, especially not for candidate SN~2009ip analogs. During 2010 to 2012 and 2014 to 2018, we only have a few detections, which is insufficient to come to any conclusion as to a possible less eruptive phase of the transient. We find steady dimming after the most recent re-brightening events and possible evidence of porosity in the circumstellar envelope, suggesting AT~2000ch may currently be in transition to a state of relative calm. We identified a second, unrelated source at a projected distance of $\sim 23\,$pc ($\sim0.5^{\prime\prime}$) that has contaminated the optical measurements of AT~2000ch at its minimum luminosity over the last two decades probably on a $5\%-10\,\%$ level, but this does not affect our overall findings and is negligible during re-brightening. We are able to restrict the mass-loss rate of AT~2000ch to range between several $10^{-6}\,\textrm{M}_{\odot}/\textrm{yr}$ and several $10^{-5}\,\textrm{M}_{\odot}/\textrm{yr}$. The fresh ejecta appear to be optically thick to radio continuum emission at least within the first $\sim 25\,$days after significant re-brightening., Comment: Accepted for publication on A&A on 23/10/2022

Multi-wavelength observations indicate that some starburst galaxies show a dominant non-thermal contribution from their central region. These active galactic nuclei (AGN)-starburst composites are of special interest, as both phenomena on their own are potential sources of highly-energetic cosmic rays and associated gamma-ray and neutrino emission. In this work, a homogeneous, steady-state two-zone multi-messenger model of the non-thermal emission from the AGN corona as well as the circumnuclear starburst region is developed and subsequently applied to the case of NGC 1068, which has recently shown some first indications of high-energy neutrino emission. Here, we show that the entire spectrum of multi-messenger data - from radio to gamma-rays including the neutrino constraint - can be described very well if both, starburst and AGN corona, are taken into account. Using only a single emission region is not sufficient.

Observations of the neutral atomic hydrogen (HI) in the nuclear starburst galaxy NGC 4945 with MeerKAT are presented. We find a large amount of halo gas, previously missed by HI observations, accounting for 6.8% of the total HI mass. This is most likely gas blown into the halo by star formation. Our maps go down to a $3\sigma$ column density level of $5\times10^{18} cm^{-2}$ . We model the HI distribution using tilted-ring fitting techniques and find a warp on the galaxy's approaching and receding sides. The HI in the northern side of the galaxy appears to be suppressed. This may be the result of ionisation by the starburst activity in the galaxy, as suggested by a previous study. The origin of the warp is unclear but could be due to past interactions or ram pressure stripping. Broad, asymmetric HI absorption lines extending beyond the HI emission velocity channels are present towards the nuclear region of NGC 4945. Such broad lines suggest the existence of a nuclear ring moving at a high circular velocity. This is supported by the clear rotation patterns in the HI absorption velocity field. The asymmetry of the absorption spectra can be caused by outflows or inflows of gas in the nuclear region of NGC 4945. The continuum map shows small extensions on both sides of the galaxy's major axis that might be signs of outflows resulting from the starburst activity., Comment: 20 pages, 15 figures, 3 tables, accepted for publication in MNRAS

Recent MeerKAT neutral hydrogen (HI) observations of Fornax A reveal tidal material intersecting in projection the western lobe of this radio galaxy. We found a spatial coincidence between the northern HI tail and a depolarized structure observed for the first time with the Australian Square Kilometre Array Pathfinder (ASKAP) at 1.2 GHz. We analyzed the properties of the rotation measure (RM) image obtained with ASKAP data at the location of the HI tail and in its neighborhood. We modeled the observed RM structure function to investigate the magnetic field power spectrum at the location of the HI tail and in a nearby control region. We found that the observed RM, in the control region and in a region enclosing the HI tail, cannot be due to the intracluster Faraday screen caused by the Fornax cluster. An intragroup magnetized medium with a central magnetic field strength of 18.5 $\rm\mu$G can explain the control region RM, but it is clear that there is an excess in correspondence with the HI tail region. We evaluated several scenarios in which the HI tail is either in the lobe foreground or embedded in the lobe. We determined a magnetic field strength on the order of $\sim$9.5$-$11 $\mu$G in the HI tail, a value consistent with constraints derived from narrowband H$\alpha$ imaging of the ionized gas. The spatial coincidence between HI tail and depolarization analyzed in this paper could be the first observed evidence of a magnetic field that either has passed through a radio galaxy lobe or has survived the lobe expansion., Comment: 11 pages, 9 figures, 1 table, accepted by A&A

Ram-pressure stripping is a crucial evolutionary driver for cluster galaxies and jellyfish galaxies, characterized by very extended tails of stripped gas, are the most striking examples of it in action. Recently, those extended tails are found to show on-going star formation raising the question how the stripped, cold gas can survive long enough to form new stars outside the stellar disk. In this work, we summarize the most recent results achieved within the GASP collaboration to provide a holistic explanation for this phenomenon. We focus on two textbook examples of jellyfish galaxies, JO206 and JW100 for which, via multi-wavelength observations from radio to X-ray and numerical simulations, we have explored the different gas phases (neutral, molecular, diffuse-ionized, and hot). Based on additional multi-phase gas studies we now propose a scenario of stripped tail evolution including all phases that is driven by a magnetic draping sheath, where the intracluster turbulent magnetized plasma condenses onto the galaxy disk and tail and produces a magnetized interface that protects the stripped galaxy tail gas from evaporation. In such a scenario, the accreted environmental plasma can cool down and eventually join the tail gas, hence providing additional gas to form stars. The implications of our findings can shed light on the more general scenario of draping, condensation, and cooling of hot gas surrounding cold clouds that is fundamental in many astrophysical phenomena.

We present deep JVLA observations at 1.4 GHz and 2.7 GHz (full polarization), as well as optical OmegaWINGS/WINGS and X-ray observations of two extended radio galaxies in the IIZW108 galaxy cluster at z = 0.04889. They show a bent tail morphology in agreement with a radio lobed galaxy falling into the cluster potential. Both galaxies are found to possess properties comparable with {narrow-angle} tail galaxies in the literature even though they are part of a low mass cluster. We find a spectral index steepening and an increase in fractional polarization through the galaxy jets and an ordered magnetic field component mostly aligned with the jet direction. This is likely caused by either shear due to the velocity difference of the intracluster medium and the jet fluid and/or magnetic draping of the intracluster medium across the galaxy jets. We find clear evidence that one source is showing two active galactic nuclei (AGN) outbursts from which we expect the AGN has never turned off completely. We show that pure standard electron cooling cannot explain the jet length. We demonstrate therefore that these galaxies can be used as a laboratory to study gentle re-acceleration of relativistic electrons in galaxy jets via transition from laminar to turbulent motion., Comment: accepted for publication in MNRAS

Recent multi-wavelength observations indicate that some starburst galaxies show a dominant nonthermal contribution from its central region. These active galactic nuclei (AGN)-starburst composites are of special interest, as both phenomena on their own are potential sources of the high-energetic cosmic rays and their accompanied gamma-ray and neutrino emission. Here, we will focus on NGC 1068, which is known since several years from its atypical radio-gamma-ray correlation. Recently this source has also shown strong indications of high energy neutrino emission. A first semi-analytical, two-component multi-messenger model is presented that already gives some constraints on the AGN-starburst composite characteristics of NGC 1068 and exposes the need to include both starburst \emph{and} AGN corona to describe the multi-messenger data.

ESO 149-G003 is a close-by, isolated dwarf irregular galaxy. Previous observations with the ATCA indicated the presence of anomalous neutral hydrogen (HI) deviating from the kinematics of a regularly rotating disc. We conducted follow-up observations with the MeerKAT radio telescope during the 16-dish Early Science programme as well as with the MeerLICHT optical telescope. Our more sensitive radio observations confirm the presence of anomalous gas in ESO 149-G003, and further confirm the formerly tentative detection of an extraplanar HI component in the galaxy. Employing a simple tilted-ring model, in which the kinematics is determined with only four parameters but including morphological asymmetries, we reproduce the galaxy's morphology, which shows a high degree of asymmetry. By comparing our model with the observed HI, we find that in our model we cannot account for a significant (but not dominant) fraction of the gas. From the differences between our model and the observed data cube we estimate that at least 7%-8% of the HI in the galaxy exhibits anomalous kinematics, while we estimate a minimum mass fraction of less than 1% for the morphologically confirmed extraplanar component. We investigate a number of global scaling relations and find that, besides being gas-dominated with a neutral gas-to-stellar mass ratio of 1.7, the galaxy does not show any obvious global peculiarities. Given its isolation, as confirmed by optical observations, we conclude that the galaxy is likely currently acquiring neutral gas. It is either re-accreting gas expelled from the galaxy or accreting pristine intergalactic material., Comment: 16 pages, 14 figures, MNRAS, Accepted 2020 December 2, in original form 2019 September 18

Jellyfish galaxies have long tails of gas that is stripped from the disc by ram pressure due to the motion of galaxies in the intracluster medium in galaxy clusters. We present the first measurement of the magnetic field strength and orientation within the disc and the (90$\,$kpc-long) $\rm H\alpha$-emitting tail of the jellyfish galaxy JO206. The tail has a large-scale magnetic field ($>4.1\,\mu$G), a steep radio spectral index ($\alpha \sim -2.0$), indicating an aging of the electrons propagating away from the star-forming regions, and extremely high fractional polarisation ($>50\,$%), indicating low turbulent motions. The magnetic field vectors are aligned with (parallel to) the direction of the ionised-gas tail and stripping direction. High-resolution simulations of a large, cold gas cloud that is exposed to a hot, magnetised turbulent wind show that the high fractional polarisation and the ordered magnetic field can be explained by accretion of draped magnetised plasma from the hot wind that condenses onto the external layers of the tail, where it is adiabatically compressed and sheared. The ordered magnetic field, preventing heat and momentum exchange, may be a key factor in allowing in-situ star formation in the tail., Comment: accepted for publication in Nature Astronomy, Sharedlt: https://rdcu.be/b9c0L (free read only)

We report the detection of small scale bending waves, also known as corrugations, in the dust lanes of five nearby edge-on disc galaxies. This phenomenon, where the disc mid-plane bends to become wavy, just as in warps but on a smaller scale, is seen here for the first time, in the dust lanes running across the discs. Because they are seen in absorption, this feature must be present in the dust disc in the outskirts of these galaxies. We enhance the visibility of these features using unsharp masking, trace the dust mid-plane across the disc, measure the corrugation amplitude by eye and the corrugation wavelength using Fourier analysis. The corrugation amplitude is found to be in the range of 70 - 200 pc and the wavelengths lie between 1 - 5 kpc. In this limited sample, we find that the amplitude of the corrugations tends to be larger for lower mass galaxies, whereas the wavelength of corrugation does not seem to depend on the mass of host galaxies. Linear stability analysis is performed to find out the dynamical state of these dust discs. Based on WKB analysis, we find that the dust corrugations in about half of our sample are stable. Further analysis, on a larger sample would be useful to strengthen the above results., Comment: 11 pages, 10 figures, accepted for publication by MNRAS

Context. The magnetic field in spiral galaxies is known to have a large-scale spiral structure along the galactic disk and is observed as X-shaped in the halo of some galaxies. While the disk field can be well explained by dynamo action, the 3-dimensional structure of the halo field and its physical nature is still unclear. Aims. As first steps towards understanding the halo fields, we want to clarify whether the observed X-shaped field is a wide-spread pattern in the halos of spiral galaxies and whether these halo fields are just turbulent fields ordered by compression or shear (anisotropic turbulent fields), or have a large-scale regular structure. Methods. The analysis of the Faraday rotation in the halo is the tool to discern anisotropic turbulent fields from large-scale magnetic fields. This, however, has been challenging until recently because of the faint halo emission in linear polarization. Our sensitive VLA broadband observations C-band and L-band of 35 spiral galaxies seen edge-on (called CHANG-ES) allowed us to perform RM-synthesis in their halos and to analyze the results. We further accomplished a stacking of the observed polarization maps of 28 CHANG-ES galaxies at C-band. Results. Though the stacked edge-on galaxies were of different Hubble types, star formation and interaction activities, the stacked image clearly reveals an X-shaped structure of the apparent magnetic field. We detected a large-scale (coherent) halo field in all 16 galaxies that have extended polarized intensity in their halos. We detected large-scale field reversals in all of their halos. In six galaxies they are along lines about vertical to the galactic midplane (vertical RMTL) with about 2 kpc separation. Only in NGC 3044 and possibly in NGC 3448 we observed vertical giant magnetic ropes (GMRs) similar to those detected recently in NGC 4631., Comment: 40 pages, 29 figures, 5 tables, accepted for publication by A&A. Conclusions in abstract abridged to meet ArXiv size limit

Reliable identification of the origin of the high-energy non-thermal emission from the Galactic Center (GC) is not achievable without adequate consideration of the ambient conditions such as the magnetic field configuration or gas distribution. In a first step, we present a model that can explain the diffuse gamma-ray emission as measured by H.E.S.S. for small longitudes in the Galactic Center region but comes to grief with higher longitudes. The model is given via the solution of a transport equation that allows for a radial dependency of the mass distribution. In order to move from this semi-analytical approximation toward a full understanding of the PeVatron signature, we present a new 3D analytical model of the gas distribution in the Central Molecular Zone (CMZ). Furthermore, we derive for the first time a 3D model of the magnetic field configuration and strength in the CMZ, which is analytical and divergence-free. The model is built via a combination of a model for the diffuse inter-cloud medium, local molecular clouds and non-thermal filaments at which local information are based on investigations from previous works and the molecular gas density. It can be shown that without an efficient longitudinal CR entrapment, a single source at the center does not facilely suffice the diffuse gamma-ray detection. Further, we show that the new magnetic field model GBFD19 is compatible with recent polarization data and has a significant impact on the longitudinal profiles of CR propagation., Comment: Proceeding for the ICRC 2019 (in Madison)

Context: Cosmic-ray propagation is strongly dependent on the large-scale configuration of the Galactic magnetic field. In particular, the Galactic center (GC) region provides highly interesting cosmic-ray data from gamma-ray maps and it is clear that a large fraction of the cosmic rays detected at Earth originate in this region of the Galaxy. Yet because of confusion from line-of-sight integration, the magnetic field structure in the Galactic center is not well known and no large-scale magnetic field model exists at present. Aims: In this paper, we develop a magnetic field model, derived from observational data on the diffuse gas, nonthermal radio filaments, and molecular clouds. Methods. We derive an analytical description of the magnetic field structure in the central molecular zone by combining observational data with the theoretical modeling of the basic properties of magnetic fields. Results. We provide a first description of the large-scale magnetic field in the Galactic center region. We present first test simulations of cosmic-ray propagation and the impact of the magnetic field structure on the cosmic-ray distribution in the three dimensions. Conclusions: Our magnetic field model is able to describe the main features of polarization maps; it is particularly important to note that they are significantly better than standard global Galactic magnetic field models. It can also be used to model cosmic-ray propagation in the Galactic Center region more accurately., Comment: 21 pages, 29 figures

The CHANG-ES galaxy sample consists of 35 nearby edge-on galaxies that have been observed using the VLA at 1.6 GHz and 6.0 GHz. Here we present the 3rd data release of our sample, namely the B-configuration 1.6 GHz sample. In addition, we make available the {\it band-to-band} spectral index maps between 1.6 GHz and 6.0 GHz, the latter taken in the matching resolution C-configuration. The images can be downloaded from https://www.queensu.ca/changes. These are our highest resolution images ($\approx$ 3 arcsec) and we examine the possible presence of low luminosity active galactic nuclei in the sample as well as some in-disk structure. New features can be seen in the spectral index maps that are masked in the total intensity emission, including hidden spiral arms in NGC~3448 and two previously unknown radio lobes on either side of the nucleus of NGC~3628. Our AGN detection rate, using only radio criteria, is 55\% which we take as a lower limit because some weaker embedded AGNs are likely present which could be revealed at higher resolution. Archival XMM-Newton data were used to search for further fingerprints of the AGNs in the studied sample. In galaxy disks, discrete regions of flat spectral index are seen, likely due to a thermal emission fraction that is higher than the global average., Comment: 36 pages, 7 figures plus a series of Panel images for Appendix B. Accepted to AJ, May 3, 2019 Version 2 changes the reference to previous observations of N4438 on p. 19

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 centered 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., Comment: 16 pages, accepted for publication in Galaxies: "New Perspectives on Galactic Magnetism", eds. Ann Mao and Andrew Fletcher

Galactic all-sky maps at very disparate frequencies, like in the radio and $\gamma$-ray regime, show similar morphological structures. This mutual information reflects the imprint of the various physical components of the interstellar medium. We want to use multifrequency all-sky observations to test resolution improvement and restoration of unobserved areas for maps in certain frequency ranges. For this we aim to reconstruct or predict from sets of other maps all-sky maps that, in their original form, lack a high resolution compared to other available all-sky surveys or are incomplete in their spatial coverage. Additionally, we want to investigate the commonalities and differences that the ISM components exhibit over the electromagnetic spectrum. We build an $n$-dimensional representation of the joint pixel-brightness distribution of $n$ maps using a Gaussian mixture model and see how predictive it is: How well can one map be reproduced based on subsets of other maps? Tests with mock data show that reconstructing the map of a certain frequency from other frequency regimes works astonishingly well, predicting reliably small-scale details well below the spatial resolution of the initially learned map. Applied to the observed multifrequency data sets of the Milky Way this technique is able to improve the resolution of, e.g., the low-resolution Fermi LAT maps as well as to recover the sky from artifact-contaminated data like the ROSAT 0.855 keV map. The predicted maps generally show less imaging artifacts compared to the original ones. A comparison of predicted and original maps highlights surprising structures, imaging artifacts (fortunately not reproduced in the prediction), and features genuine to the respective frequency range that are not present at other frequency bands. We discuss limitations of this machine learning approach and ideas how to overcome them.

We detect 5 galaxies in the CHANG-ES (Continuum Halos in Nearby Galaxies -- an EVLA Survey) sample that show circular polarization (CP) at L-band in our high resolution data sets. Two of the galaxies (NGC~4388 and NGC~4845) show strong Stokes $V/I\,\equiv\,m_C\,\sim\,2$\%, two (NGC~660 and NGC~3628) have values of $m_C\sim \,0.3$\%, and NGC~3079 is a marginal detection at $m_C\sim \,0.2$\%. The two strongest $m_C$ galaxies also have the most luminous X-ray cores and the strongest internal absorption in X-rays. We have expanded on our previous Faraday conversion interpretation and analysis and provide analytical expressions for the expected $V$ signal for a general case in which the cosmic ray electron energy spectral index can take on any value. We provide examples as to how such expressions could be used to estimate magnetic field strengths and the lower energy cutoff for CR electrons. Four out of our detections are {\it resolved}, showing unique structures, including a {\it jet} in NGC~4388 and a CP `conversion disk' in NGC~4845. The conversion disk is inclined to the galactic disk but is perpendicular to a possible outflow direction. Such CP structures have never before been seen in any galaxy to our knowledge. None of the galaxy cores show linear polarization at L-band. Thus CP may provide a unique probe of physical conditions deep into radio AGNs., Comment: 30 pages, 4 figures, accepted to MNRAS

The vertical halo scale height is a crucial parameter to understand the transport of cosmic-ray electrons (CRE) and their energy loss mechanisms in spiral galaxies. Until now, the radio scale height could only be determined for a few edge-on galaxies because of missing sensitivity at high resolution. We developed a sophisticated method for the scale height determination of edge-on galaxies. With this we determined the scale heights and radial scale lengths for a sample of 13 galaxies from the CHANG-ES radio continuum survey in two frequency bands. The sample average value for the radio scale heights of the halo are 1.1 +/- 0.3kpc in C-band and 1.4 +/- 0.7kpc in L-band. From the frequency dependence analysis of the halo scale heights we found that the wind velocities (estimated using the adiabatic loss time) are above the escape velocity. We found that the halo scale heights increase linearly with the radio diameters. In order to exclude the diameter dependence, we defined a normalized scale height which is quite similar for all sample galaxies at both frequency bands and does not depend on the star formation rate or the magnetic field strength. However, the normalized scale height shows a tight anticorrelation with the mass surface density. The sample galaxies with smaller scale lengths are more spherical in the radio emission, while those with larger scale lengths are flatter. The radio scale height depends mainly on the radio diameter of the galaxy. The sample galaxies are consistent with an escape-dominated radio halo with convective cosmic ray propagation, indicating that galactic winds are a widespread phenomenon in spiral galaxies. While a higher star formation rate or star formation surface density does not lead to a higher wind velocity, we deceleration of CRE outflow, e.g. a lowering of the wind velocity from the galactic disk., Comment: accepted for publication in Astronomy & Astrophysics

The luminosity ratio of electrons to protons as it is produced in stochastic acceleration processes in cosmic ray sources is an important quantity relevant for several aspects of the modeling of the sources themselves. It is usually assumed to be around 1:100 in the case of Galactic sources, while a value of 1:10 is typically assumed when describing extragalactic sources. It is supported by observations that the average ratios should be close to these values. At this point, however, there is no possibility to investigate how each individual source behaves. When looking at the physics aspects, a 1:100 ratio is well supported in theory when making the following assumptions: (1) the total number of electrons and protons that is accelerated are the same; (2) the spectral index of both populations after acceleration is $a_e=a_p\approx 2.2$. In this paper, we reinvestigate these assumptions. In particular, assumption (2) is not supported by observational data of the sources and PIC simulation yield different spectral indices as well. We present the detailed calculation of the electron-to-proton ratio, dropping the assumption of equal spectral indices. We distinguish between the ratio of luminosities and the ratio of the differential spectral behavior, which becomes necessary for cases where the spectral indices of the two particle populations are not the same. We discuss the possible range of values when allowing for different spectral indices concerning the spectral behavior of electrons and protons. Additionally, it is shown that the minimum energy of the accelerated population can have a large influence on the results., Comment: 30 pages, 12 figures, accepted for publication in Astroparticle Physics (02/24/2017)

We present 1D cosmic ray transport models, numerically solving equations of pure advection and diffusion for the electrons and calculating synchrotron emission spectra. We find that for exponential halo magnetic field distributions advection leads to approximately exponential radio continuum intensity profiles, whereas diffusion leads to profiles that can be better approximated by a Gaussian function. Accordingly, the vertical radio spectral profiles for advection are approximately linear, whereas for diffusion they are of `parabolic' shape. We compare our models with deep ATCA observations of two edge-on galaxies, NGC 7090 and 7462, at $\lambda\lambda$ 22 and 6 cm. Our result is that the cosmic ray transport in NGC 7090 is advection dominated with $V=150^{+80}_{-30}~\rm km\,s^{-1}$, and that the one in NGC 7462 is diffusion dominated with $D=3.0\pm 1.0 \times 10^{28}E_{\rm GeV}^{0.5}~\rm cm^2\,s^{-1}$. NGC 7090 has both a thin and thick radio disc with respective magnetic field scale heights of $h_{\rm B1}=0.8\pm 0.1$ kpc and $h_{\rm B2}=4.7\pm 1.0$ kpc. NGC 7462 has only a thick radio disc with $h_{\rm B2}=3.8\pm 1.0$ kpc. In both galaxies, the magnetic field scale heights are significantly smaller than what estimates from energy equipartition would suggest. A non-negligible fraction of cosmic ray electrons can escape from NGC 7090, so that this galaxy is not an electron calorimeter., Comment: MNRAS accepted. 19 pages, 10 figures, 5 tables plus 4 pages of appendix

We present the first part of the observations made for the Continuum Halos in Nearby Galaxies, an EVLA Survey (CHANG-ES) project. The aim of the CHANG-ES project is to study and characterize the nature of radio halos, their prevalence as well as their magnetic fields, and the cosmic rays illuminating these fields. This paper reports observations with the compact D configuration of the Karl G. Jansky Very Large Array (VLA) for the sample of 35 nearby edge-on galaxies of CHANG-ES. With the new wide bandwidth capabilities of the VLA, an unprecedented sensitivity was achieved for all polarization products. The beam resolution is an average of 9.6" and 36" with noise levels reaching approximately 6 and 30 microJy per beam for C- and L-bands, respectively (robust weighting). We present intensity maps in these two frequency bands (C and L), with different weightings, as well as spectral index maps, polarization maps, and new measurements of star formation rates (SFRs). The data products described herein are available to the public in the CHANG-ES data release available at www.queensu.ca/changes. We also present evidence of a trend among galaxies with larger halos having higher SFR surface density, and we show, for the first time, a radio continuum image of the median galaxy, taking advantage of the collective signal-to-noise ratio of 30 of our galaxies. This image shows clearly that a typical spiral galaxy is surrounded by a halo of magnetic fields and cosmic rays., Comment: 70 pages, of which 35 pages present the data of each galaxy

Magnetic fields are an important ingredient of the interstellar medium (ISM). Besides their importance for star formation, they govern the transport of cosmic rays, relevant to the launch and regulation of galactic outflows and winds, which in turn are pivotal in shaping the structure of halo magnetic fields. Mapping the small-scale structure of interstellar magnetic fields in many nearby galaxies is crucial to understand the interaction between gas and magnetic fields, in particular how gas flows are affected. Elucidation of the magnetic role in, e.g., triggering star formation, forming and stabilising spiral arms, driving outflows, gas heating by reconnection and magnetising the intergalactic medium has the potential to revolutionise our physical picture of the ISM and galaxy evolution in general. Radio polarisation observations in the very nearest galaxies at high frequencies (>= 3 GHz) and with high spatial resolution (<= 5") hold the key here. The galaxy survey with SKA1 that we propose will also be a major step to understand the galactic dynamo, which is important for models of galaxy evolution and for astrophysical magnetohydrodynamics in general. Field amplification by turbulent gas motions, which is crucial for efficient dynamo action, has been investigated so far only in simulations, while compelling evidence of turbulent fields from observations is still lacking., Comment: 18 pages, 3 figures; to appear as part of 'Cosmic Magnetism' in Proceedings 'Advancing Astrophysics with the SKA' (AASKA14) as article PoS(AASKA14)094

Magnetic fields play an important role in shaping the structure and evolution of the interstellar medium (ISM) of galaxies, but the details of this relationship remain unclear. With SKA1, the 3D structure of galactic magnetic fields and its connection to star formation will be revealed. A highly sensitive probe of the internal structure of the magnetoionized ISM is the partial depolarization of synchrotron radiation from inside the volume. Different configurations of magnetic field and ionized gas within the resolution element of the telescope lead to frequency-dependent changes in the observed degree of polarization. The results of spectro-polarimetric observations are tied to physical structure in the ISM through comparison with detailed modeling, supplemented with the use of new analysis techniques that are being actively developed and studied within the community such as Rotation Measure Synthesis. The SKA will enable this field to come into its own and begin the study of the detailed structure of the magnetized ISM in a sample of nearby galaxies, thanks to its extraordinary wideband capabilities coupled with the combination of excellent surface brightness sensitivity and angular resolution., Comment: 11 pages, 1 figure; to appear as part of 'Cosmic Magnetism' in Proceedings 'Advancing Astrophysics with the SKA (AASKA14)', PoS(AASKA14)106

Simulations of galaxy growth need to invoke strong negative feedback from active galactic nuclei (AGN) to suppress the formation of stars and thus prevent the over-production of very massive systems. While some observations provide evidence for such negative feedback, other studies find either no feedback, or even positive feedback, with increased star formation associated with higher AGN luminosities. Here we report an analysis of several hundred AGN and their host galaxies in the Chandra Deep Field South using X-ray and radio data for sample selection. Combined with archival far infrared data as a reliable tracer of star formation activity in the AGN host galaxies, we find that AGN with pronounced radio jets exhibit a much higher star formation rate than the purely X-ray selected ones, even at the same X-ray luminosities. This difference implies that positive AGN feedback plays an important role, too, and therefore has to be accounted for in all future simulation work. We interpret this to indicate that the enhanced star formation rate of radio selected AGN arises because of jet-induced star formation, as is hinted by the different jet powers among our AGN samples, while the suppressed star formation rate of X-ray selected AGN is caused by heating and photo-dissociation of molecular gas by the hot AGN accretion disc., Comment: 21 pages, 9 figures, 3 tables, accepted for publication by ApJ

The global properties of the interstellar medium with processes such as infall and outflow of gas and a large scale circulation of matter and its consequences for star formation and chemical enrichment are important for the understanding of galaxy evolution. In this paper we studied the kinematics and morphology of the diffuse ionized gas (DIG) in the disk and in the halo of the star forming spiral galaxy NGC 4666 to derive information about its kinematical properties. Especially, we searched for infalling and outflowing ionized gas. We determined surface brightness, radial velocity, and velocity dispersion of the warm ionized gas via high spectral resolution (R ~ 9000) Fabry-P\'erot interferometry. This allows the determination of the global velocity field and the detection of local deviations from this verlocity field. We calculated models of the DIG distribution and its kinematics for comparison with the measured data. In this way we determined fundamental parameters such as the inclination and the scale height of NGC 4666, and established the need for an additional gas component to fit our observed data. We found individual areas, especially along the minor axis, with gas components reaching into the halo which we interpret as an outflowing component of the diffuse ionized gas. As the main result of our study, we were able to determine that the vertical structure of the DIG distribution in NGC 4666 is best modeled with two components of ionized gas, a thick and a thin disk with 0.8 kpc and 0.2 kpc scale height, respectively. Therefore, the enhanced star formation in NGC 4666 drives an outflow and also maintains a thick ionized gas layer reminiscent of the Reynold's layer in the Milky Way., Comment: 12 pages, 10 figures, 3 tables

We present the first results from the CHANG-ES survey, a new survey of 35 edge-on galaxies to search for both in-disk as well as extra-planar radio continuum emission. The motivation and science case for the survey are presented in a companion paper (Paper I). In this paper (Paper II), we outline the observations and data reduction steps required for wide-band calibration and mapping of EVLA data, including polarization, based on C-array test observations of NGC 4631. With modest on-source observing times (30 minutes at 1.5 GHz and 75 minutes at 6 GHz for the test data) we have achieved best rms noise levels of 22 and 3.5 $\mu$Jy beam$^{-1}$ at 1.5 GHz and 6 GHz, respectively. New disk-halo features have been detected, among them two at 1.5 GHz that appear as loops in projection. We present the first 1.5 GHz spectral index map of NGC 4631 to be formed from a single wide-band observation in a single array configuration. This map represents tangent slopes to the intensities within the band centered at 1.5 GHz, rather than fits across widely separated frequencies as has been done in the past and is also the highest spatial resolution spectral index map yet presented for this galaxy. The average spectral index in the disk is $\bar\alpha_{1.5 GHz}\,=\,-0.84\,\pm\,0.05$ indicating that the emission is largely non-thermal, but a small global thermal contribution is sufficient to explain a positive curvature term in the spectral index over the band. Two specific star forming regions have spectral indices that are consistent with thermal emission. Polarization results (uncorrected for internal Faraday rotation) are consistent with previous observations and also reveal some new features. On broad scales, we find strong support for the notion that magnetic fields constrain the X-ray emitting hot gas., Comment: Accepted to the Astronomical Journal, Version 2 changes: Added acknowledgement to NRAO

We introduce a new survey to map the radio continuum halos of a sample of 35 edge-on spiral galaxies at 1.5 GHz and 6 GHz in all polarization products. The survey is exploiting the new wide bandwidth capabilities of the Karl G. Jansky Very Large Array (i.e. the Expanded Very Large Array, or EVLA) in a variety of array configurations (B, C, and D) in order to compile the most comprehensive data set yet obtained for the study of radio halo properties. This is the first survey of radio halos to include all polarization products. In this first paper, we outline the scientific motivation of the survey, the specific science goals, and the expected improvements in noise levels and spatial coverage from the survey. Our goals include investigating the physical conditions and origin of halos, characterizing cosmic ray transport and wind speed, measuring Faraday rotation and mapping the magnetic field, probing the in-disk and extraplanar far-infrared - radio continuum relation, and reconciling non-thermal radio emission with high-energy gamma-ray models. The sample size allows us to search for correlations between radio halos and other properties, including environment, star formation rate, and the presence of AGNs. In a companion paper (Paper II) we outline the data reduction steps and present the first results of the survey for the galaxy, NGC 4631., Comment: 17 pages, 1 figure, accepted to the Astronomical Journal, Version 2 changes: added acknowledgement to NRAO

Magnetic fields are a good tracer for gas compression by shock waves, which can be caused by interaction of star-formation driven outflows from individual star formation sites as described in the chimney model. We study the magnetic field structure in the central part of the nuclear starburst galaxy NGC 253 with spatial resolutions between 40 and 150 pc to detect any filamentary emission associated with the nuclear outflow. New VLA observations at 3 cm with 7.5" resolution were combined with archive data at 20 and 6 cm. We find filamentary radio continuum emission in a geometrical distribution that we interpret as the boundary of the northwestern nuclear outflow cone. The scaleheight of the continuum emission is 150+/-20 pc, regardless of the observing frequency. The equipartition magnetic field strength is 46+/-10 microG for the total field and 21+/-5 microG for the regular field in the filaments. The ordered magnetic field is aligned along the filaments, in agreement with amplification due to compression. The perpendicular diffusion coefficient across the filaments is kappa_perp = 1.5 x 10^28 cm^2 s^-1 E(GeV)^(0.5+/-0.7). In the SE part of the nuclear outflow cone the magnetic field is pointing away from the disc in form of a helix, with an azimuthal component increasing up to at least 1200 pc height, where it is about equal to the total component. The ordered magnetic field in the disc is anisotropic within a radius of 2.2 kpc. At larger radii, the large-scale field is regular and of even parity. The magnetic field is able to collimate the outflow, which can explain the observed small opening angle of ~26 degree. Due to angular momentum conservation, the field lines are frozen into the plasma and are wound up into a helix. Strong adiabatic losses of the cosmic-ray electrons can partly explain why the radio luminosity of the nucleus lies below the radio-FIR correlation., Comment: 18 pages, 17 figures, replaced with version published in Astronomy and Astrophysics

We present a classification-based approach to identify quasi-stellar radio sources (quasars) in the Sloan Digital Sky Survey and evaluate its performance on a manually labeled training set. While reasonable results can already be obtained via approaches working only on photometric data, our experiments indicate that simple but problem-specific features extracted from spectroscopic data can significantly improve the classification performance. Since our approach works orthogonal to existing classification schemes used for building the spectroscopic catalogs, our classification results are well suited for a mutual assessment of the approaches' accuracies., Comment: 6 pages, 8 figures, published in proceedings of 2010 Ninth International Conference on Machine Learning and Applications (ICMLA) of the IEEE

In order to better determine the large-scale 3D-structure of magnetic fields in spiral galaxies we present a Faraday rotation analysis of the edge-on spiral galaxy NGC 5775. Deep radio-continuum observations in total power and linear polarization were performed at 8.46 GHz with the VLA and the 100-m Effelsberg telescope. They were analyzed together with archival 4.86 and 1.49 GHz VLA-data. We thus can derive rotation measures from a comparison of three frequencies and determine the intrinsic magnetic field structure. A very extended halo is detected in NGC 5775, with magnetic field lines forming an X-shaped structure. Close to the galactic disk the magnetic field is plane-parallel. The scaleheights of the radio emission esimated for NGC 5775 are comaprable with other galaxies. The rotation measure distribution varies smoothly on both sides along the major axis from positive to negative values. From the derived distribution of rotation measures and the plane-parallel intrinsic magnetic field orientation along the galactic midplane we conclude that NGC 5775 has an 'even axisymmetric' large-scale magnetic field configuration in the disk as generated by an \alpha \Omega -dynamo which is accompanied by a quadrupolar poloidal field. The magnetic field lines of the plane-parallel component are pointing 'outwards'. The observed X-shaped halo magnetic field, however, cannot be explained by the action of the disk's mean-field dynamo alone. It is probably due to the influence of the galactic wind together with the dynamo action., Comment: 11 pages, 13 figures, A&A accepted