The grand-design face-on spiral galaxy M51 is an excellent laboratory for studying magnetic fields in galaxies. We present new observations of M51 using the VLA at the frequency range of S-band (2-4GHz), to shed new light on the transition region between the disk and halo. We present images of the distributions of the total intensity, polarized intensity, degree of polarization, and rotation measure (RM). The RM distribution in S-band shows a fluctuating pattern without any apparent large-scale structure. We discuss a model of the depolarization of synchrotron radiation in a multi-layer magneto-ionic medium and compare the model predictions to the polarization data of M51 between 1-8GHz. Since the model predictions strongly differ within the wavelength range of the S-band, the new data are essential. The parameters of the model are adjusted to fit to the data of polarization fractions in a few selected regions. In three spiral arm regions, the turbulent field in the disk dominates with strengths between 18muG and 24muG, while the regular field strengths are 8-16muG. In one inter-arm region, the regular field strength of 18muG exceeds that of the turbulent field of 11muG. The regular field strengths in the halo are 3-5muG. The observed RMs in the disk-halo transition region are probably dominated by tangled regular fields, as predicted from models of evolving dynamos, and/or vertical fields, as predicted from numerical simulations of Parker instabilities or galactic winds. Both types of magnetic fields have frequent reversals on scales similar to or larger than the beam size (550pc) that contribute to an increase of the RM dispersion and to distortions of any large-scale pattern of the regular field. Our study devises new ways of analyzing and interpreting broadband multi-frequency polarization data that will be applicable to future data from, for example, the Square Kilometre Array.

M33 is a gas rich spiral galaxy of the Local Group. We investigate the relationship between the two major gas cooling lines and the total infrared (TIR) dust continuum. We mapped the emission of gas and dust in M33 using the far-infrared lines of [CII] and [OI](63um) and the TIR. The line maps were observed with Herschel/PACS. These maps have 50pc resolution and form a ~370pc wide stripe along its major axis covering the sites of bright HII regions, but also more quiescent arm and inter-arm regions from the southern arm at 2kpc galacto-centric distance to the south out to 5.7kpc distance to the north. Full-galaxy maps of the continuum emission at 24um from Spitzer/MIPS, and at 70um, 100um, and 160um from PACS were combined to obtain a map of the TIR. TIR and [CII] intensities are correlated over more than two orders of magnitude. The range of TIR translates to a range of far ultraviolet (FUV) emission of G0,obs~2 to 200 in units of the average Galactic radiation field. The binned [CII]/TIR ratio drops with rising TIR, with large, but decreasing scatter. Fits of modified black bodies (MBBs) to the continuum emission were used to estimate dust mass surface densities and total gas column densities. A correction for possible foreground absorption by cold gas was applied to the [OI] data before comparing it with models of photon dominated regions (PDRs). Most of the ratios of [CII]/[OI] and ([CII]+[OI])/TIR are consistent with two model solutions. The median ratios are consistent with one solution at n~2x10^2 cm-3, G0~60, and and a second low-FUV solution at n~10^4 cm-3, G0~1.5. The bulk of the gas along the lines-of-sight is represented by a low-density, high-FUV phase with low beam filling factors ~1. A fraction of the gas may, however, be represented by the second solution.

We investigate thermal and non-thermal radio continuum associated with the early formation and evolution of Young Stellar Clusters (YSCs) selected by their MIR emission in M33. For the first time in an external galaxy it has been possible to identify radio counterparts to more than 300 star forming regions. We proof the nature of candidate YSCs fully embedded in molecular clouds, by recovering their associated faint radio continuum luminosities. Using the Halpha line to identify free-free radio emission at 5 GHz in the more evolved, partially exposed YSCs, we retrieve information on the relevance of magnetic fields and cosmic rays across the M33 disk at 25 pc spatial scales. A cross-correlation of MIR and radio continuum luminosities is established from bright to very faint YSCs, with MIR-to-radio emission ratio showing a gradual decline towards the outer disk, while the magnetic field is pervasive at all radii. We establish and discuss the tight relation between radio continuum and other star formation indicators, such as Halpha. This relation holds for individual YSCs over four orders of magnitude as well as for molecular clouds hosting YSCs. On average about half of radio emission at 5 GHz in YSCs is non-thermal. For exposed but compact YSCs the non-thermal radio fraction increases with source brightness, while for large HII regions the fraction is lower and shows no clear trend. This has been found for YSCs with and without identified SNRs and underlines the possible role of massive stars in triggering particle acceleration through winds and shocks: these particles diffuse throughout the native molecular cloud prior to cloud dispersal.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

Magnetic fields constitute an energetic component of the interstellar medium in galaxies and hence can affect the formation of galactic structures. Sensitive resolved radio continuum observations together with statistical studies in galaxy samples are performed to investigate the origin and the impact of the magnetic fields. The JVLA cloud-scale survey of the local group galaxy, M33, unveils strong tangled magnetic field along the spiral arms and in the galaxy center indicating amplification due to compression and local shear in molecular clouds. Studying a sample of non-cluster, non-interacting galaxies, we find that the large-scale ordered magnetic field scales with the rotation speed of galaxies. The total and disordered magnetic fields scales with the star formation rate in normal star forming galaxies. On the other hand, a strong magnetic field in the center of NGC 1097 -- a massive galaxy undergoing quenching -- is found responsible for decelerating its massive star formation. Putting these results together, it is deduced that 1) the Universe was highly magnetized short after the peak of the massive star birth at about 1 Gyr after the Big Bang and 2) the strong magnetic field has possibly acted as back reaction after this time, quenching the massive star formation and stimulating the creation of low-mass stars in massive galaxies.

Correlations between the radio continuum, infrared and CO emission are known to exist for several types of galaxies and across several orders of magnitude. However, the low-mass, low-luminosity and low-metallicity regime of these correlations is not well known. A sample of metal-rich and metal-poor dwarf galaxies from the literature has been assembled to explore this extreme regime. The results demonstrate that the properties of dwarf galaxies are not simple extensions of those of more massive galaxies; the different correlations reflect different star-forming conditions and different coupling between the star formation and the various quantities. It is found that dwarfs show increasingly weaker CO and infrared emission for their luminosity, as expected for galaxies with a low dust content, slower reaction rates, and a hard ionizing radiation field. In the higher-luminosity dwarf regime (L_1.4GHz > 10^27 W, where L_1.4GHz ~ 10^29 W for a Milky Way star formation rate of ~1 M_sun yr^-1), the total and non-thermal radio continuum emission appear to adequately trace the star formation rate. A breakdown of the dependence of the (Halpha-based) thermal, non-thermal, and, hence, total radio continuum emission on star formation rate occurs below L_1.4GHz ~ 10^27 W, resulting in a steepening or downturn of the relations at extreme low luminosity. Below L_FIR ~ 10^36 W ~ 3 x 10^9 L_sun, the infrared emission ceases to adequately trace the star formation rate. A lack of a correlation between the magnetic field strength and the star formation rate in low star formation rate dwarfs suggests a breakdown of the equipartition assumption. As extremely metal-poor dwarfs mostly populate the low star formation rate and low luminosity regime, they stand out in their infrared, radio continuum and CO properties.

Radio continuum observations have proven to be a workhorse in our understanding of the star formation process (i.e., stellar birth and death) from galaxies both in the nearby universe and out to the highest redshifts. In this article we focus on how the ngVLA will transform our understanding of star formation by enabling one to map and decompose the radio continuum emission from large, heterogeneous samples of nearby galaxies on $\gtrsim 10$\,pc scales to conduct a proper accounting of the energetic processes powering it. At the discussed sensitivity and angular resolution, the ngVLA will simultaneously be able to create maps of current star formation activity at $\sim$100\,pc scales, as well as detect and characterize (e.g., size, spectral shape, density, etc.) discrete H{\sc ii} regions and supernova remnants on 10\,pc scales in galaxies out to the distance of the Virgo cluster. Their properties can then be used to see how they relate to the local and global ISM and star formation conditions. Such investigations are essential for understanding the astrophysics of high-$z$ measurements of galaxies, allowing for proper modeling of galaxy formation and evolution.

An excellent laboratory for studying large scale magnetic fields is the grand de- sign face-on spiral galaxy M51. Due to wavelength-dependent Faraday depolarization, linearly polarized synchrotron emission at different radio frequencies gives a picture of the galaxy at different depths: Observations at L-band (1-2 GHz) probe the halo region while at C- and X- band (4-8 GHz) the linearly polarized emission probe the disk region of M51. We present new observations of M51 using the Karl G. Jansky Very Large Array (VLA) at S-band (2-4 GHz), where previously no polarization observations existed, to shed new light on the transition region between the disk and the halo. We discuss a model of the depolarization of synchrotron radiation in a multilayer magneto-ionic medium and compare the model predictions to the multi-frequency polarization data of M51 between 1-8GHz. The new S-band data are essential to distinguish between different models. Our study shows that the initial model parameters, i.e. the total reg- ular and turbulent magnetic field strengths in the disk and halo of M51, need to be adjusted to successfully fit the models to the data.

Low-frequency radio continuum observations of edge-on galaxies are ideal to study cosmic-ray electrons (CREs) in halos via radio synchrotron emission and to measure magnetic field strengths. We obtained new observations of the edge-on spiral galaxy NGC 891 at 129-163 MHz with the LOw Frequency ARray (LOFAR) and at 13-18 GHz with the Arcminute Microkelvin Imager (AMI) and combine them with recent high-resolution Very Large Array (VLA) observations at 1-2 GHz, enabling us to study the radio continuum emission over two orders of magnitude in frequency. The spectrum of the integrated nonthermal flux density can be fitted by a power law with a spectral steepening towards higher frequencies or by a curved polynomial. Spectral flattening at low frequencies due to free-free absorption is detected in star-forming regions of the disk. The mean magnetic field strength in the halo is 7 +- 2 $\mu$G. The scale heights of the nonthermal halo emission at 146 MHz are larger than those at 1.5 GHz everywhere, with a mean ratio of 1.7 +- 0.3, indicating that spectral ageing of CREs is important and that diffusive propagation dominates. The halo scale heights at 146 MHz decrease with increasing magnetic field strengths which is a signature of dominating synchrotron losses of CREs. On the other hand, the spectral index between 146 MHz and 1.5 GHz linearly steepens from the disk to the halo, indicating that advection rather than diffusion is the dominating CRE transport process. This issue calls for refined modelling of CRE propagation.

Low-mass galaxies are subject to strong galactic outflows, in which cosmic rays may play an important role, they can be best traced with low-frequency radio continuum observations, which are less affected by spectral ageing. We present a study of the nearby star burst dwarf irregular galaxy IC 10 using observations at 140 MHz with the LOw-Frequency ARray (LOFAR), at 1580 MHz with the Very Large Array (VLA) and at 6200 MHz with the VLA and the 100-m Effelsberg telescope. We find that IC 10 has a low-frequency radio halo, which manifests itself as a second component (thick disc) in the minor axis profiles of the non-thermal radio continuum emission at 140 and 1580 MHz. These profiles are then fitted with 1D cosmic-ray transport models for pure diffusion and advection. We find that a diffusion model fits best, with a diffusion coefficient of $D=(0.4$-$0.8) \times 10^{26}(E/{\rm GeV})^{0.5}~{\rm cm^2\,s^{-1}}$, which is at least an order of magnitude smaller than estimates both from anisotropic diffusion and the diffusion length. In contrast, advection models, which cannot be ruled out due to the mild inclination, while providing poorer fits, result in advection speeds close to the escape velocity of $\approx$$50~\rm km\,s^{-1}$, as expected for a cosmic-ray driven wind. Our favoured model with an accelerating wind provides a self-consistent solution, where the magnetic field is in energy equipartition with both the warm neutral and warm ionized medium with an important contribution from cosmic rays. Consequently, cosmic rays can play a vital role for the launching of galactic winds in the disc--halo interface.

The relative abundance of the dust grain types in the interstellar medium (ISM) is directly linked to physical quantities that trace the evolution of galaxies. We study the dust properties of the whole disc of M33 at spatial scales of ~170 pc. This analysis allows us to infer how the relative dust grain abundance changes with the conditions of the ISM, study the existence of a submillimetre excess and look for trends of the gas-to-dust mass ratio (GDR) with other physical properties of the galaxy. For each pixel in the disc of M33 we fit the infrared SED using a physically motivated dust model that assumes an emissivity index beta close to 2. We derive the relative amount of the different dust grains in the model, the total dust mass, and the strength of the interstellar radiation field (ISRF) heating the dust at each spatial location. The relative abundance of very small grains tends to increase, and for big grains to decrease, at high values of Halpha luminosity. This shows that the dust grains are modified inside the star-forming regions, in agreement with a theoretical framework of dust evolution under different physical conditions. The radial dependence of the GDR is consistent with the shallow metallicity gradient observed in this galaxy. The strength of the ISRF derived in our model correlates with the star formation rate in the galaxy in a pixel by pixel basis. Although this is expected it is the first time that a correlation between both quantities is reported. We produce a map of submillimetre excess in the 500 microns SPIRE band for the disc of M33. The excess can be as high as 50% and increases at large galactocentric distances. We further study the relation of the excess with other physical properties of the galaxy and find that the excess is prominent in zones of diffuse ISM outside the main star-forming regions, where the molecular gas and dust surface density are low.

Observations show that massive star formation quenches first at centers of galaxies. To understand quenching mechanisms, we investigate the thermal and nonthermal energy balance in the central kpc of NGC1097- a prototypical galaxy undergoing quenching- and present a systematic study of the nuclear star formation efficiency and its dependencies. This region is dominated by the nonthermal pressure from the magnetic field, cosmic rays, and turbulence. A comparison of the mass-to-magnetic flux ratio of the molecular clouds shows that most of them are magnetically critical or supported against gravitational collapse needed to form cores of massive stars. Moreover, the star formation efficiency of the clouds drops with the magnetic field strength. Such an anti-correlation holds with neither the turbulent nor the thermal pressure. Hence, a progressive built up of the magnetic field results in high-mass stars forming inefficiently, and it may be the cause of the low-mass stellar population in the bulges of galaxies.

We study the spectral energy distribution (SED) of the radio continuum emission from the KINGFISH sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4GHz, 4.8GHz, 8.5GHz, and 10.5GHz combined with archive data allow us, for the first time, to determine the mid-radio continuum (1-10 GHz, MRC) bolometric luminosities and further present calibration relations vs. the monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo (MCMC) technique leading to measurements for the nonthermal spectral index and the thermal fraction f_th with mean values of alpha_nt=0.97+-0.16 (0.79+-0.15 for the total spectral index) and f_th= 10% +- 9% at 1.4 GHz. The MRC luminosity changes over ~3 orders of magnitude in the sample. The thermal emission is responsible for ~23% of the MRC on average. We also compare the extinction-corrected diagnostics of star formation rate with the thermal and nonthermal radio tracers and derive the first star formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing star formation rate surface density, indicating the effect of the star formation feedback on the cosmic ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with star formation rate, due to the amplification of the magnetic fields in star forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star forming galaxies are detected.

Recent studies show the importance of the star formation feedback in changing the energetic and structure of galaxies. Dissecting the physics of the feedback is hence crucial to understand the evolution of galaxies. Full polarization radio continuum surveys can be ideally performed to trace not only star formation but also the energetic components of the interstellar medium (ISM), the magnetic fields and cosmic ray electrons. Using the SKA precursors, we investigate the effect of the massive star formation on the ISM energy balance in nearby galaxies. Our multi-scale and multi-frequency surveys show that cosmic rays are injected in star forming regions and lose energy propagating away from their birth place. Due to the star formation feedback, cosmic ray electron population becomes younger and more energetic. Star formation also amplifies the turbulent magnetic field inserting a high pressure which is important in energy balance in the ISM and structure formation in the host galaxy.

We have conducted a near-infrared monitoring campaign at the UK InfraRed Telescope (UKIRT), of the Local Group spiral galaxy M33 (Triangulum). On the basis of their variability, we have identified stars in the very final stage of their evolution, and for which the luminosity is more directly related to the birth mass than the more numerous less-evolved giant stars that continue to increase in luminosity. In this fifth paper of the series, we construct the birth mass function and hence derive the star formation history across the galactic disc of M33. The star formation rate has varied between $\sim0.010\pm0.001$ ($\sim0.012\pm0.007$) and 0.060$\pm0.005$ (0.052$\pm0.009$)M$_\odot$ yr$^{-1}$ kpc$^{-2}$ statistically (systematically) in the central square kiloparsec of M33, comparable with the values derived previously with another camera. The total star formation rate in M33 within a galactocentric radius of 14 kpc has varied between $\sim0.110\pm0.005$ ($\sim0.174\pm0.060$) and $\sim0.560\pm0.028$ ($\sim0.503\pm0.100$) M$_\odot$ yr$^{-1}$ statistically (systematically). We find evidence of two epochs during which the star formation rate was enhanced by a factor of a few -- one that started $\sim6$ Gyr ago and lasted $\sim3$ Gyr and produced $\geq71$% of the total mass in stars, and one $\sim250$ Myr ago that lasted $\sim200$ Myr and formed $\leq13$% of the mass in stars. Radial star formation history profiles suggest that the inner disc of M33 was formed in an inside-out formation scenario. The outskirts of the disc are dominated by the old population, which may be the result of dynamical effects over many Gyr. We find correspondence to spiral structure for all stars, but enhanced only for stars younger than $\sim100$ Myr; this suggests that the spiral arms are transient features and not part of a global density wave potential.

[Abridged] Do some environments favor efficient conversion of molecular gas into stars? To answer this, we need to be able to estimate the H2 mass. Traditionally, this is done using CO and a few assumptions but the Herschel observations in the FIR make it possible to estimate the molecular gas mass independently of CO. Previous attempts to derive gas masses from dust emission suffered from biases. Generally, dust surface densities, HI column densities, and CO intensities are used to derive a gas-to-dust ratio (GDR) and the local CO intensity to H2 column density ratio (XCO), sometimes allowing for an additional CO-dark gas component (Kdark). We tested earlier methods, revealing degeneracies among the parameters, and then used a Bayesian formalism to derive the most likely values for each of the parameters mentioned above as a function of position in the nearby low metallicity spiral galaxy M33. The data are from the IRAM 30m CO(2-1) line, high-resolution HI and Herschel dust continuum observations. Solving for GDR, XCO, and Kdark in macro pixels 500 pc in size, we find that (i) allowing for CO-dark gas significantly improves fits; (ii) Kdark decreases with galactocentric distance; (iii) GDR is slightly higher than initially expected and increases with galactocentric distance; (iv) the total amount of dark gas closely follows the radially decreasing CO emission, as might be expected if the dark gas is H2 where CO is photodissociated. The total amount of H2, including dark gas, yields an average XCO of twice the galactic value of 2e20 cm^-2/(K km/s), 55% of this traced directly through CO. The rather constant fraction of dark gas suggests that there is no large population of diffuse H2 clouds (unrelated to GMCs) without CO emission. Unlike in large spirals, we detect no systematic radial trend in XCO, possibly linked to the absence of a radial decrease in CO line ratios.

The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code.

Recent studies show the importance of feedback in the evolution of the star formation rate in the Universe. However, the nature and physics of the feedback are still pressing questions. Radio continuum observations can provide unique dust-unbiased tracers of massive star formation and of the interstellar medium (ISM) and hence are ideal to address the regulation of star formation in galaxies. Our multi-frequency and multi-resolution radio surveys in nearby galaxies enable us to trace various phases of star formation and dissect the thermal and nonthermal ISM in galaxies. Mapping the cosmic ray electron energy index and magnetic field strength, we have found observational evidence that massive star formation significantly affects the energy balance in the ISM through the injection and acceleration of cosmic rays and the amplification of magnetic fields. How the next generation of stars could form in such a magnetized and turbulent ISM will be addressed in our 'EVLA cloud-scale survey of the local group galaxy M33' and in forthcoming surveys with the SKA.

To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quantifying what physical processes drive the variations between different hybrid estimators. Doing so, we aim at deriving new universal UV+IR hybrid estimators to correct the UV for dust attenuation, taking into account the intrinsic physical properties of galaxies. We use the CIGALE code to model the spatially-resolved FUV to FIR SED of eight nearby star-forming galaxies drawn from the KINGFISH sample. This allows us to determine their local physical properties, and in particular their UV attenuation, average SFR, average specific SFR (sSFR), and their stellar mass. We then examine how hybrid estimators depend on said properties. We find that hybrid UV+IR estimators strongly depend on the stellar mass surface density (in particular at 70 and 100 micron) and on the sSFR (in particular at 24 micron and the TIR). Consequently, the IR scaling coefficients for UV obscuration can vary by almost an order of magnitude. This result contrasts with other groups who found relatively constant coefficients with small deviations. We exploit these variations to construct a new class of hybrid estimators based on observed UV to near-IR colours and near-IR luminosity densities per unit area. We find that they can reliably be extended to entire galaxies. The new estimators provide better estimates of attenuation-corrected UV emission than classical hybrid estimators. Naturally taking into account the variable impact of dust heated by old stellar populations, they constitute a step towards universal estimators.

The origin and evolution of cosmic magnetic fields as well as the influence of the magnetic fields on the evolution of galaxies are unknown. Though not without challenges, the dynamo theory can explain the large-scale coherent magnetic fields which govern galaxies, but observational evidence for the theory is so far very scarce. Putting together the available data of non-interacting, non-cluster galaxies with known large-scale magnetic fields, we find a tight correlation between the integrated polarized flux density, S(PI), and the rotation speed, v(rot), of galaxies. This leads to an almost linear correlation between the large-scale magnetic field B and v(rot), assuming that the number of cosmic ray electrons is proportional to the star formation rate, and a super-linear correlation assuming equipartition between magnetic fields and cosmic rays. This correlation cannot be attributed to an active linear alpha-Omega dynamo, as no correlation holds with global shear or angular speed. It indicates instead a coupling between the large-scale magnetic field and the dynamical mass of the galaxies, B ~ M^(0.25-0.4). Hence, faster rotating and/or more massive galaxies have stronger large-scale magnetic fields. The observed B-v(rot) correlation shows that the anisotropic turbulent magnetic field dominates B in fast rotating galaxies as the turbulent magnetic field, coupled with gas, is enhanced and ordered due to the strong gas compression and/or local shear in these systems. This study supports an stationary condition for the large-scale magnetic field as long as the dynamical mass of galaxies is constant.

Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.

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.

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.

In this chapter, we highlight a number of science investigations that are enabled by the inclusion of Band~5 ($4.6-13.8$ GHz) for SKA1-MID science operations, while focusing on the astrophysics of star formation over cosmic time. For studying the detailed astrophysics of star formation at high-redshift, surveys at frequencies $\gtrsim$10 GHz have the distinct advantage over traditional $\sim$1.4 GHz surveys as they are able to yield higher angular resolution imaging while probing higher rest frame frequencies of galaxies with increasing redshift, where emission of star-forming galaxies becomes dominated by thermal (free-free) radiation. In doing so, surveys carried out at $\gtrsim$10 GHz provide a robust, dust-unbiased measurement of the massive star formation rate by being highly sensitive to the number of ionizing photons that are produced. To access this powerful star formation rate diagnostic requires that Band~5 be available for SKA1-MID. We additionally present a detailed science case for frequency coverage extending up to 30 GHz during full SKA2 operations, as this allows for highly diverse science while additionally providing contiguous frequency coverage between the SKA and ALMA, which will likely be the two most powerful interferometers for the coming decades. To enable this synergy, it is crucial that the dish design of the SKA be flexible enough to include the possibility of being fit with receivers operating up to 30 GHz.

We present a study of the magneto-ionic medium in the Whirlpool galaxy (M51) using new wide-band multi-configuration polarization data at L band (1-2 GHz) obtained at the Karl G. Jansky Very Large Array. By fitting the observed diffuse complex polarization $Q$+$iU$ as a function of wavelength directly to various depolarization models, we find that polarized emission from M51 at 1-2 GHz originates from the top of the synchrotron disk and then experiences Faraday rotation in the near-side thermal halo of the galaxy. Thus, the scale height of the thermal gas must exceed that of the synchrotron emitting gas at L band. The observed Faraday depth distribution at L band is consistent with a halo field that comprises of a plane-parallel bisymmetric component and a vertical component which produces a Faraday rotation of $\sim$ $-$9 rad m$^{-2}$. The derived rotation measure structure functions indicate a characteristic scale of rotation measure fluctuations of less than 560 pc in the disk and approximately 1 kpc in the halo. The outer scale of turbulence of 1 kpc found in the halo of M51 is consistent with superbubbles and the Parker instability being the main energy injection mechanisms in galactic halos.

The [CII] 157.74 $\mu$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $\Sigma_{\rm [CII]}-\Sigma_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $\Sigma_{\rm[CII]}-\Sigma_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $\Sigma_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.

We investigate the far-infrared (IR) dust emission for 20 local star forming galaxies from the Key Insights on Nearby Galaxies: A Far-IR Survey with Herschel (KINGFISH) sample. We model the far-IR/submillimeter spectral energy distribution (SED) using images from Spitzer Space Telescope and Herschel Space Observatory. We calculate the cold dust temperature (T(cold)) and emissivity (beta) on a pixel by pixel basis (where each pixel ranges from 0.1-3 kpc^2) using a two temperature modified blackbody fitting routine. Our fitting method allows us to investigate the resolved nature of temperature and emissivity variations by modeling from the galaxy centers to the outskirts (physical scales of ~15-50 kpc, depending on the size of the galaxy). We fit each SED in two ways: (1) fit T(cold) and beta simultaneously, (2) hold beta constant and fit T(cold). We compare T(cold) and beta with star formation rates (calculated from L(Halpha) and L(24)), the luminosity of the old stellar population (traced through L(3.6), and the dust mass surface density (traced by 500 micron luminosity, L(500)). We find a significant trend between SFR/L(500) and T(cold), implying that the flux of hard UV photons relative to the amount of dust is significantly contributing to the heating of the cold, or diffuse, dust component. We also see a trend between L(3.6)/L(500) and beta, indicating that the old stellar population contributes to the heating at far-IR/submillimeter wavelengths. Finally, we find that when beta is held constant, T(cold) exhibits a strongly decreasing radial trend, illustrating that the shape of the far-IR SED is changing radially through a galaxy, thus confirming on a sample almost double in size the trends observed in Galametz et al. (2012).

NGC 2146, a nearby luminous infrared galaxy (LIRG), presents evidence for outflows along the disk minor axis in all gas phases (ionized, neutral atomic and molecular). We present an analysis of the multi-phase starburst driven superwind in the central 5 kpc as traced in spatially resolved spectral line observations, using far-IR Herschel PACS spectroscopy, to probe the effects on the atomic and ionized gas, and optical integral field spectroscopy to examine the ionized gas through diagnostic line ratios. We observe an increased ~250 km/s velocity dispersion in the [OI] 63 micron, [OIII] 88 micron, [NII] 122 micron and [CII] 158 micron fine-structure lines that is spatially coincident with high excitation gas above and below the disk. We model this with a slow ~200 km/s shock and trace the superwind to the edge of our field of view 2.5 kpc above the disk. We present new SOFIA 37 micron observations to explore the warm dust distribution, and detect no clear dust entrainment in the outflow. The stellar kinematics appear decoupled from the regular disk rotation seen in all gas phases, consistent with a recent merger event disrupting the system. We consider the role of the superwind in the evolution of NGC 2146 and speculate on the evolutionary future of the system. Our observations of NGC 2146 in the far-IR allow an unobscured view of the wind, crucial for tracing the superwind to the launching region at the disk center, and provide a local analog for future ALMA observations of outflows in high redshift systems.

We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temperature. Both beta and the dust temperature are higher in the inner disk than in the outer disk, contrary to reported beta-T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Halpha) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate fashion. The flocculent spiral arms and the dust lanes are evident in the map of the cold component. Both cold and warm dust column densities are high in star forming regions and reach their maxima toward the giant star forming complexes NGC604 and NGC595. beta declines from close to 2 in the center to about 1.3 in the outer disk. beta is positively correlated with star formation and with molecular gas column, as traced by Halpha and CO emission. The lower dust emissivity index in the outer parts of M33 is likely related to the reduced metallicity (different grain composition) and possibly different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the FIR-bright regions in the outer disk also have a low beta. Like most spirals, M33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Halpha or CO tend to trace the inner disk, making it difficult to distinguish between their effects on the dust.

We derive the magnitude of fluctuations in total synchrotron intensity in the Milky Way and M33, from both observations and theory under various assumption about the relation between cosmic rays and interstellar magnetic fields. Given the relative magnitude of the fluctuations in the Galactic magnetic field (the ratio of the rms fluctuations to the mean magnetic field strength) suggested by Faraday rotation and synchrotron polarization, the observations are inconsistent with local energy equipartition between cosmic rays and magnetic fields. Our analysis of relative synchrotron intensity fluctuations indicates that the distribution of cosmic rays is nearly uniform at the scales of the order of and exceeding $100\p$, in contrast to strong fluctuations in the interstellar magnetic field at those scales. A conservative upper limit on the ratio of the the fluctuation magnitude in the cosmic ray number density to its mean value is 0.2--0.4 at scales of order 100\,pc. Our results are consistent with a mild anticorrelation between cosmic-ray and magnetic energy densities at these scales, in both the Milky Way and M33. Energy equipartition between cosmic rays and magnetic fields may still hold, but at scales exceeding 1\,kpc. Therefore, we suggest that equipartition estimates be applied to the observed synchrotron intensity smoothed to a linear scale of kiloparsec order (in spiral galaxies) to obtain the cosmic ray distribution and a large-scale magnetic field. Then the resulting cosmic ray distribution can be used to derive the fluctuating magnetic field strength from the data at the original resolution. The resulting random magnetic field is likely to be significantly stronger than existing estimates.

We investigate the effect of propagation of cosmic-ray electrons (CRE) on the nonthermal (synchrotron) - far-infrared correlations in M31 and M33. The thermal (TH) and nonthermal (NTH) emission components of the radio continuum emission at 1.4 GHz and one higher frequency are compared with dust emission from M31 and M33 using Spitzer data. In both galaxies the TH emission is linearly correlated with the emission from warm dust (24 \mu m, 70 \mu m), but the power laws of the NTH-FIR correlations have exponents b < 1 that increase with increasing frequency. Furthermore, the values of b for M33 are significantly smaller (b ~ 0.4) than those for M31 (b ~ 0.6). We interpret the differences in b as differences in the diffusion length of the CRE. We estimate the diffusion length in two ways: (1) by smoothing the NTH emission at the higher frequency until the correlation with NTH emission at 1.4 GHz has b = 1, and (2) by smoothing the TH emission until the correlation with the NTH emission at the same frequency has b = 1, assuming that the TH emission represents the source distribution of the CRE. Our smoothing experiments show that M31 only has a thin NTH disk with a scale height of h = 0.3-0.4 kpc at 1.4 GHz, whereas M33 has a similar thin disk as well as a thick disk with scale height h_thick ~ 2 kpc. In the thin disks, the (deprojected) diffusion length at 1.4 GHz is ~ 1.5 kpc, yielding a diffusion coefficient of ~ 2 10^28 cm^2/s. The structure, strength and regularity of the magnetic field in a galaxy as well as the existence of a thick disk determine the diffusion of the CRE, and hence, the power-law exponent of the NTH-FIR correlations.

Interstellar magnetic fields and the propagation of cosmic ray electrons have an important impact on the radio-infrared (IR) correlation in galaxies. This becomes evident when studying different spatial scales within galaxies. We investigate the correlation between the infrared (IR) and free-free/synchrotron radio continuum emission at 20 cm from the two local group galaxies M31 and M33 on spatial scales between 0.4 and 10 kpc. The multi-scale radio-IR correlations have been carried out using a wavelet analysis. The free-free and IR emission are correlated on all scales, but on some scales the synchrotron emission is only marginally correlated with the IR emission. The synchrotron-IR correlation is stronger in M33 than in M31 on small scales (<1 kpc), but it is weaker than in M31 on larger scales. Taking the smallest scale on which the synchrotron-IR correlation exists as the propagation length of cosmic ray electrons, we show that the difference on small scales can be explained by the smaller propagation length in M33 than in M31. On large scales, the difference is due to the thick disk/halo in M33, which is absent in M31. A comparison of our data with data on NGC6946, the LMC and M51 suggests that the propagation length is determined by the ratio of ordered-to-turbulent magnetic field strength, which is consistent with diffusion of CR electrons in the ISM. As the diffusion length of CR electrons influences the radio-IR correlation, this dependence is a direct observational evidence of the importance of magnetic fields for the radio-IR correlation within galaxies. The star formation rate per surface area only indirectly influences the diffusion length as it increases the strength of the turbulent magnetic field.

Dust has long been identified as a barrier to measuring inherent galaxy properties. However, the link between dust and attenuation is not straightforward and depends on both the amount of dust and its distribution. Herschel imaging of nearby galaxies undertaken as part of the KINGFISH project allows us to map the dust as seen in emission with unprecedented sensitivity and ~1 kpc resolution. We present here new optical integral field unit spectroscopy for eight of these galaxies that provides complementary 100-200 pc scale maps of the dust attenuation through observation of the reddening in both the Balmer decrement and the stellar continuum. The stellar continuum reddening, which is systematically less than that observed in the Balmer decrement, shows no clear correlation with the dust, suggesting that the distribution of stellar reddening acts as a poor tracer of the overall dust content. The brightest HII regions are observed to be preferentially located in dusty regions, and we do find a correlation between the Balmer line reddening and the dust mass surface density for which we provide an empirical relation. Some of the high-inclination systems in our sample exhibit high extinction, but we also find evidence that unresolved variations in the dust distribution on scales smaller than 500 pc may contribute to the scatter in this relation. We caution against the use of integrated A_V measures to infer global dust properties.

Observations of galaxies at sub-millimeter wavelengths, where the emission is mainly due to cold dust, are required to constrain the dust physical properties and provide important insight on the gas content of galaxies. We mapped NGC1365 at 870um with LABOCA, the Large APEX Bolometer Camera, allowing us to probe the central mass concentration as well as the rate at which the gas flows to the center. We obtained the dust physical properties both globally and locally for different locations in the galaxy. A 20 K modified black body represents about 98% of the total dust content of the galaxy, the rest can be represented by a warmer dust component of 40 K. The bar exhibits an east-west asymmetry in the dust distribution: The eastern bar is heavier than the western bar by more than a factor of 4. Integrating the dust SED, we derive a total infrared (IR) luminosity of 9.8 x 10^{10} L_{\odot} leading to a dust-enshrouded star formation rate of ~16.7 M_{\odot} per year in NGC1365. We derive the gas mass from the measurements of the dust emission resulting in a CO-to-H2 conversion factor of X_{CO}~1.2x10^{20} mol cm^{-2} (K km s^{-1})^{-1} in the central disk including the bar. Taking into account the metallicity variation, the central gas mass concentration is only ~20% at R<40" (3.6 kpc). On the other hand, the time-scale with which the gas flows into the center, ~300 Myr, is rather short. This indicates that the current central mass in NGC1365 is evolving fast due to the strong bar.

We present new empirical calibrations to estimate resolved and integrated total infrared luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed with Herschel. We perform a resolved SED modelling of these objects using the Draine and Li (2007) dust models and investigate the influence of the addition of SPIRE measurements in the estimation of LTIR. We find that using data up to 250 um leads to local LTIR values consistent with those obtained with a complete coverage (up to 500 um) within 10 per cent for most of our resolved elements. We then study the distribution of energy in the resolved SEDs of our galaxies. The bulk of energy (30-50 per cent) is contained in the (70-160 um) band. The (24-70 um) fraction decreases with increasing metallicity. The (160-1100 um) submillimeter band can account for up to 25 per cent of the LTIR in metal-rich galaxies. We investigate the correlation between TIR surface brightnesses/luminosities and monochromatic Spitzer and Herschel surface brightnesses/luminosities. The three PACS bands can be used as reliable monochromatic estimators of the LTIR, the 100 um band being the most reliable monochromatic tracer. There is also a strong correlation between the SPIRE 250um and LTIR, although with more scatter than for the PACS relations. We also study the ability of our monochromatic relations to reproduce integrated LTIR of nearby galaxies as well as LTIR of z=1-3 sources. Finally, we provide calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated uncertainties.

We derive the distribution of the synchrotron spectral index across NGC6946 and investigate the correlation between the radio continuum (synchrotron) and far-infrared (FIR) emission using the KINGFISH Herschel PACS and SPIRE data. The radio--FIR correlation is studied as a function of star formation rate, magnetic field strength, radiation field strength, and the total gas surface brightness. The synchrotron emission follows both star-forming regions and the so-called magnetic arms present in the inter-arm regions. The synchrotron spectral index is steepest along the magnetic arms ($\alpha_n \sim 1$), while it is flat in places of giant H{\sc ii} regions and in the center of the galaxy ($\alpha_n \sim 0.6-0.7$). The map of $\alpha_n$ provides an observational evidence for aging and energy loss of cosmic ray electrons propagating in the disk of the galaxy. Variations in the synchrotron--FIR correlation across the galaxy are shown to be a function of both star formation and magnetic fields. We find that the synchrotron emission correlates better with cold rather than with warm dust emission, when the interstellar radiation field is the main heating source of dust. The synchrotron--FIR correlation suggests a coupling between the magnetic field and the gas density. NGC6946 shows a power-law behavior between the total (turbulent) magnetic field strength B and the star formation rate surface density $\Sigma_{\rm SFR}$ with an index of 0.14\,(0.16)$\pm$0.01. This indicates an efficient production of the turbulent magnetic field with the increasing gas turbulence expected in actively star forming regions. The scale-by-scale analysis of the synchrotron--FIR correlation indicates that the ISM affects the propagation of old/diffused cosmic ray electrons, resulting in a diffusion coefficient of $D_0=4.6\times 10^{28}$\,cm$^2$\,s$^{-1}$ for 2.2\,GeV CREs.

Within the framework of the Herschel M 33 extended survey HerM33es we study the Spectral Energy Distribution (SED) of a set of HII regions in M 33 as a function of the morphology. We present a catalogue of 119 HII regions morphologically classified: 9 filled, 47 mixed, 36 shell, and 27 clear shell HII regions. For each object we extract the photometry at twelve available wavelength bands (from FUV-1516A to IR-250mi) and obtain the SED. We also obtain emission line profiles across the regions to study the location of the stellar, ionised gas, and dust components. We find trends for the SEDs related to the morphology, showing that the star and gas-dust configuration affects the ratios of the emission in different bands. The mixed and filled regions show higher emission at 24mi than the shells and clear shells, which could be due to the proximity of the dust to the stellar clusters in the case of mixed and filled regions. The FIR peak for shells and clear shells seems to be located towards longer wavelengths, indicating that the dust is colder for this type of objects.The logarithmic 100/70mi ratio for filled and mixed regions remains constant over one order of magnitude in Halpha and FUV surface brightness, while the shells and clear shells exhibit a wider range of values of almost two orders of magnitude. We derive dust masses and temperatures fitting the individual SEDs with dust models proposed in the literature. The derived dust mass range is between 10^2-10^4 Msun and the cold dust temperature spans T(cold)~12-27 K. The spherical geometrical model proposed for the Halpha clear shells is confirmed by the emission profile obtained from the observations and is used to infer the electron density within the envelope: the typical electron density is 0.7+-0.3 cm^-3, while filled regions can reach values two to five times higher.

This work presents high spectral resolution observations of the \CII\ line at 158 \micron, one of the major cooling lines of the interstellar medium, taken with the HIFI heterodyne spectrometer on the Herschel satellite. In BCLMP 691, an \HII\ region far north (3.3 kpc) in the disk of M 33, the \CII\ and CO line profiles show similar velocities within $0.5 \kms$, while the \HI\ line velocities are systematically shifted towards lower rotation velocities by $\sim 5\kms$. Observed at the same $12"$ angular resolution, the \CII\ lines are broader than those of CO by about 50% but narrower than the \HI\ lines. The \CII\ line to far-infrared continuum ratio suggests a photoelectric heating efficiency of 1.1%. The data, together with published models indicate a UV field $G_0 \sim 100$ in units of the solar neighborhood value, a gas density $n_H \sim 1000 \cc$, and a gas temperature $T\sim 200$ K. Adopting these values, we estimate the C$^+$ column density to be $N_{C^+} \approx 1.3 \times 10^{17} \cmt$. The \CII\ emission comes predominantly from the warm neutral region between the \HII\ region and the cool molecular cloud behind it. From published abundances, the inferred C$^+$ column corresponds to a hydrogen column density of $N_H \sim 2 \times 10^{21} \cmt$. The CO observations suggest that $N_H = 2 N_{H_2} \sim 3.2 \times 10^{21} \cmt$ and 21cm measurements, also at $12"$ resolution, yield $N_\HI \approx 1.2 \times 10^{21} \cmt$ within the \CII\ velocity range. Thus, some H$_2$ not detected in CO must be present, in agreement with earlier findings based on the SPIRE 250 -- 500 $\mu$m emission.

We examine Herschel Space Observatory images of one nearby prototypical outer ring galaxy, NGC 1291, and show that the ring becomes more prominent at wavelengths longer than 160um. The mass of cool dust in the ring dominates the total dust mass of the galaxy, accounting for at least 70% of it. The temperature of the emitting dust in the ring (T=19.5+/-0.3K) is cooler than that of the inner galaxy (T=25.7+/-0.7K). We discuss several explanations for the difference in dust temperature, including age and density differences in the stellar populations of the ring versus the bulge.

Taking advantage of the sensitivity and angular resolution of the Herschel Space Observatory at far-infrared and submm wavelengths, we aim to characterize the physical properties of cold dust within nearby galaxies and study the robustness of the parameters we derive using different modified blackbody models. For a pilot subsample of the KINGFISH program, we perform 2 temperature fits of the Spitzer and Herschel photometric data (24 to 500um), with a warm and a cold component, globally and in each resolution element.At global scales, we observe ranges of values for beta_c(0.8 to 2.5) and Tc(19.1 to 25.1K).We compute maps of our parameters with beta fixed or free to test the robustness of the temperature and dust surface density maps we deduce. When the emissivity is fixed, we observe temperature gradients as a function of radius.When the emissivity is fitted as a free parameter, barred galaxies tend to have uniform fitted emissivities.Gathering resolved elements in a Tc-beta_c diagram underlines an anti-correlation between the two parameters.It remains difficult to assess whether the dominant effect is the physics of dust grains, noise, or mixing along the line of sight and in the beam. We finally observe in both cases that the dust column density peaks in central regions of galaxies and bar ends (coinciding with molecular gas density enhancements usually found in these locations).We also quantify how the total dust mass varies with our assumptions about the emissivity index as well as the influence of the wavelength coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (beta_c < 2.0) lead to significantly lower dust masses compared to the beta_c < 2.0 case, with dust masses lower by up to 50% if beta_c=1.5 for instance.The working resolution affects our total dust mass estimates: masses increase from global fits to spatially-resolved fits.

We use new Herschel multi-band imaging of the Andromeda galaxy to analyze how dust heating occurs in the central regions of galaxy spheroids that are essentially devoid of young stars. We construct a dust temperature map of M31 through fitting modified blackbody SEDs to the Herschel data, and find that the temperature within 2 kpc rises strongly from the mean value in the disk of 17 pm 1K to \sim35K at the centre. UV to near-IR imaging of the central few kpc shows directly the absence of young stellar populations, delineates the radial profile of the stellar density, and demonstrates that even the near-UV dust extinction is optically thin in M31's bulge. This allows the direct calculation of the stellar radiation heating in the bulge, U\ast(r), as a function of radius. The increasing temperature profile in the centre matches that expected from the stellar heating, i.e. that the dust heating and cooling rates track each other over nearly two orders of magnitude in U\ast. The modelled dust heating is in excess of the observed dust temperatures, suggesting that it is more than sufficient to explain the observed IR emission. Together with the wavelength dependent absorption cross section of the dust, this demonstrates directly that it is the optical, not UV, radiation that sets the heating rate. This analysis shows that neither young stellar populations nor stellar near-UV radiation are necessary to heat dust to warm temperatures in galaxy spheroids. Rather, it is the high densities of Gyr-old stellar populations that provide a sufficiently strong diffuse radiation field to heat the dust. To the extent which these results pertain to the tenuous dust found in the centres of early-type galaxies remains yet to be explored.

We study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy).

Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics.

We investigate the correlation between the far-infrared (FIR) and radio continuum emission from NGC6946 on spatial scales between 0.9 and 17 kpc. We use the Herschel PACS (70, 100, 160$\mu$m) and SPIRE (250$\mu$m) data from the KINGFISH project. Separating the free-free and synchrotron components of the radio continuum emission, we find that FIR is better correlated with the free-free than the synchrotron emission. Compared to a similar study in M33 and M31, we find that the scale dependence of the synchrotron--FIR correlation in NGC6946 is more similar to M31 than M33. The scale dependence of the synchrotron--FIR correlation can be explained by the turbulent-to-ordered magnetic field ratio or, equivalently, the diffusion length of the cosmic ray electrons in these galaxies.

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

Using the Herschel PACS and SPIRE FIR/submm data, we investigate variations in the dust spectral index $\beta$ in the nearby spiral galaxy M33 at a linear resolution of 160 pc. We use an iteration method in two different approaches, single and two-component modified black body models. In both approaches, $\beta$ is higher in the central disk than in the outer disk similar to the dust temperature. There is a positive correlation between $\beta$ and Halpha as well as with the molecular gas traced by CO(2-1). A Monte-Carlo simulation shows that the physical parameters are better constrained when using the two-component model.

Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region.

We observed M51 at three frequencies, 1.4GHz (20cm), 4.9GHz (6cm) and 8.4GHz (3.6cm), with the VLA and the Effelsberg 100m telescope to obtain the highest quality radio continuum images of a nearby spiral galaxy. These radio data were combined with deconvolved Spitzer IRAC 8mum and MIPS 24mum images to search for and investigate local changes in the radio-IR correlation. Utilizing wavelet decomposition, we compare the distribution of the radio and IR emission on spatial scales between 200pc and 30kpc. We show that the radio-IR correlation is not uniform across the galactic disk. It presents a complex behavior with local extrema corresponding to various galactic structures, such as complexes of HII regions, spiral arms and interarm filaments, indicating that the contribution of the thermal and non-thermal radio emission is a strong function of environment. In particular, the relation of the 24mum and 20cm emission presents a linear relation within the spiral arms and globally over the galaxy, while it deviates from linearity in the interarm and outer regions as well in the inner region, with two different behaviours: it is sublinear in the interarm and outer region and over-linear in the central 3.5kpc. Our analysis suggests that the changes in the radio/IR correlation reflect variations of ISM properties between spiral arms and interarm region. The good correlation in the spiral arms implies that 24mum and 20cm are tracing recent star formation, while a change in the dust opacity, 'Cirrus' contribution to the IR emission and/or the relation between the magnetic field strength and the gas density can explain the different relations found in the interarm, outer and inner regions.

We obtained Spitzer/IRAC 3.6-8 micron images of the nearby spiral galaxy NGC 4258 to study possible interactions between dust and the radio jet. In our analysis we also included high-resolution radio continuum, H-alpha, CO, and X-ray data. Our data reveal that the 8 micron emission, believed to originate largely from PAH molecules and hot dust, is an excellent tracer of the normal spiral structure in NGC 4258, and hence it originates from the galactic plane. We investigated the possibility of dust destruction by the radio jet by calculating correlation coefficients between the 8 micron and radio continuum emissions along the jet in two independent ways, namely (i) from wavelet-transformed maps of the original images at different spatial scales, and (ii) from one-dimensional intensity cuts perpendicular to the projected path of the radio jet on the sky. No definitive sign of a correlation (or anticorrelation) was detected on relevant spatial scales with either approach, implying that any dust destruction must take place at spatial scales that are not resolved by our observations.

Does star formation proceed in the same way in large spirals such as the Milky Way and in smaller chemically younger galaxies? Earlier work suggests a more rapid transformation of H$_2$ into stars in these objects but (1) a doubt remains about the validity of the H$_2$ mass estimates and (2) there is currently no explanation for why star formation should be more efficient. M~33, a local group spiral with a mass $\sim 10$\% and a metallicity half that of the Galaxy, represents a first step towards the metal poor Dwarf Galaxies. We have searched for molecular clouds in the outer disk of M~33 and present here a set of detections of both $^{12}$CO and $^{13}$CO, including the only detections (for both lines) beyond the R$_{25}$ radius in a subsolar metallicity galaxy. The spatial resolution enables mass estimates for the clouds and thus a measure of the $N({\rm H}_2) / I_{\rm CO}$ ratio, which in turn enables a more reliable calculation of the H$_2$ mass. Our estimate for the outer disk of M~33 is $N({\rm H}_2) / I_{\rm CO(1-0)} \sim 5 \times 10^{20} \,{\rm cm^{-2}/(K{\rm \ km \ s^{-1}})}$ with an estimated uncertainty of a factor $\le 2$. While the $^{12/13}$CO line ratios do not provide a reliable measure of $N({\rm H}_2) / I_{\rm CO}$, the values we find are slightly greater than Galactic and corroborate a somewhat higher $N({\rm H}_2) / I_{\rm CO}$ value. Comparing the CO observations with other tracers of the interstellar medium, no reliable means of predicting where CO would be detected was identified. In particular, CO detections were often not directly on local HI or FIR or H$\alpha$ peaks, although generally in regions with FIR emission and high HI column density. The results presented here provide support for the quicker transformation of H$_2$ into stars in M~33 than in large local universe spirals.

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $\sigma$ using SPIRE and recent CO and \HI\ observations; a variation in $\sigma$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250um in order to investigate the nature of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also focus on a more precise study of some striking Halpha shells in the northern part of the galaxy. We present a catalogue of 159 compact emission sources in M33 identified by SExtractor in the 250um SPIRE band that is the one that provides the best spatial resolution. We also measured fluxes at 24um and Halpha for those 159 extracted sources. The morphological study of the shells also benefits from a multiwavelength approach including Halpha, far-UV from GALEX, and infrared from both Spitzer IRAC 8um and MIPS 24um in order to make comparisons. For the 159 compact sources selected at 250um, we find a very strong Pearson correlation coefficient with the MIPS 24um emission (r24 = 0.94) and a rather strong correlation with the Halpha emission, although with more scatter (rHa = 0.83). The morphological study of the Halpha shells shows a displacement between far-ultraviolet, Halpha, and the SPIRE bands. The cool dust emission from SPIRE clearly delineates the Halpha shell structures. The very strong link between the 250um compact emission and the 24um and Halpha emissions, by recovering the star formation rate from standard recipes for HII regions, allows us to provide star formation rate calibrations based on the 250um compact emission alone. The different locations of the Halpha and far-ultraviolet emissions with respect to the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Halpha shells and the associated cool dust.

Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies.

We derive distributions of dust temperature and dust opacity across M31 at 45" resolution using the Spitzer data. With the opacity map and a standard dust model we de-redden the Ha emission yielding the first de-reddened Ha map of M31. We compare the emissions from dust, Ha, HI and H2 by means of radial distributions, pixel-to-pixel correlations and wavelet cross-correlations. The dust temperature steeply decreases from 30K near the center to 15K at large radii. The mean dust optical depth at the Ha wavelength along the line of sight is about 0.7. The radial decrease of the dust-to-gas ratio is similar to that of the oxygen abundance. On scales<2kpc, cold dust emission is best correlated with that of neutral gas and warm dust emission with that of ionized gas. Ha emission is slightly better correlated with emission at 70um than at 24um. In the area 6kpc<R< 17kpc, the total SFR is ~0.3Msun/yr. The Kennicutt-Schmidt law between SFR and total gas has a power-law index of 1.30+-0.05 in the radial range of R=7-11kpc increasing by about 0.3 for R=11-13kpc. The lack of H2 in the central region could be related to the lack of HI and the low opacity/high temperature of the dust. Since neither SFR nor SFE is well correlated with the surface density of H2 or total gas, other factors than gas density must play an important role in the formation of massive stars in M31. The molecular depletion time scale of 1.1 Gyr indicates that M31 is about three times less efficient in forming young massive stars than M33.

Among various approaches for indirect detection of dark matter, synchrotron emission due to secondary electrons/positrons produced in galactic WIMPs annihilation is raising an increasing interest. In this paper we propose a new method to derive bounds in the mchi - <sigmaA*v> plane by using radio continuum observations of Messier 33, paying particular attention to a low emitting Radio Cavity. The comparison of the expected radio emission due to the galactic dark matter distribution with the observed one provides bounds which are comparable to those obtained from a similar analysis of the Milky Way. Remarkably, the present results are simply based on archival data and thus largely improvable by means of specifically tailored observations. The potentiality of the method compared with more standard searches is discussed by considering the optimistic situation of a vanishing flux (within the experimental sensitivity) measured inside the cavity by a high resolution radio telescope like ALMA. Under the best conditions our technique is able to produce bounds which are comparable to the ones expected after five years of Fermi LAT data taking for an hadronic annihilation channel. Furthermore, it allows to test the hypothesis that space telescopes like Pamela and Fermi LAT are actually observing electrons and positrons due to galactic dark matter annihilation into leptons.

The radio-infrared correlation was explained as a direct and linear relationship between star formation and IR emission. However, one fact making the IR-star formation linkage less obvious is that the IR emission consists of at least two emission components, cold dust and warm dust. The cold dust emission may not be directly linked to the young stellar population. Furthermore, understanding the origin of the radio-IR correlation requires to discriminate between the two main components of the radio continuum emission, free-free and synchrotron emission. Here, we present a multi-scale study of the correlation of IR with both the thermal and non-thermal (synchrotron) components of the radio continuum emission from the nearby galaxies M33 and M31.

Using high-resolution data of the linearly polarized intensity and polarization angle at 3.6, 6.2, and 20 cm together with a 3-D model of the regular magnetic field, we study variations of the structure, strength, and energy density of the magnetic field in the Scd galaxy M33. The regular magnetic field consists of a horizontal component (represented by an axisymmetric mode from 1 to 3 kpc radius and a superposition of axisymmetric and bisymmetric modes from 3 to 5 kpc radius) and a vertical component. However, the inferred `vertical field' may be partly due to a galactic warp. We estimate the average total and regular magnetic field strengths as ~ 6.4 and 2.5 $\mu$G, respectively. Generation of interstellar magnetic fields by turbulent gas motion in M33 is indicated as the turbulent and magnetic energy densities are about equal.

Using the linearly polarized intensity and polarization angle data at 3.6, 6.2 and 20 cm, we determine variations of Faraday rotation and depolarization across the nearby galaxy M33. A 3-D model of the regular magnetic field is fitted to the observed azimuthal distribution of polarization angles. Faraday rotation, measured between 3.6 and 6.2 cm at a linear resolution of 0.7 kpc, shows more variation in the south than in the north of the galaxy. About 10% of the nonthermal emission from M33 at 3.6 cm is polarized. We estimate the average total and regular magnetic field strengths in M33 as ~ 6.4 and 2.5 $\mu$G, respectively. Under the assumption that the disk of M33 is flat, the regular magnetic field consists of horizontal and vertical components: however the inferred vertical field may be partly due to a galactic warp. The horizontal field is represented by an axisymmetric (m=0) mode from 1 to 3 kpc radius and a superposition of axisymmetric and bisymmetric (m=0+1) modes from 3 to 5 kpc radius. An excess of differential Faraday rotation in the southern half together with strong Faraday dispersion in the southern spiral arms seem to be responsible for the north-south asymmetry in the observed wavelength dependent depolarization. The presence of an axisymmetric m=0 mode of the regular magnetic field in each ring suggests that a galactic dynamo is operating in M33. The pitch angles of the spiral regular magnetic field are generally smaller than the pitch angles of the optical spiral arms but are twice as big as simple estimates based on the mean-field dynamo theory and M33's rotation curve. Generation of interstellar magnetic fields from turbulent gas motions in M33 is indicated by the equipartition of turbulent and magnetic energy densities.

We determine the variation in the nonthermal radio spectral index in the nearby spiral galaxy M33. We separate the thermal and nonthermal components of the radio continuum emission without the assumption of a constant nonthermal spectral index. Using the Spitzer FIR data at 70 and 160 $\mu$m and a standard dust model, we de-redden the H$\alpha$ emission. The extinction-corrected H$\alpha$ emission serves as a template for the thermal free-free radio emission. Subtracting from the observed 3.6 and 20cm emission (Effelsberg and the VLA) this free-free emission, we obtain maps of the nonthermal intensity and spectral index.

The origin of the tight radio--IR correlation in galaxies has not been fully understood. One reason is the uncertainty about which heating sources (stars or diffuse interstellar radiation field)provide the energy that is absorbed by dust and re-radiated in IR. Another problem is caused by comparing the IR emission with the thermal and nonthermal components of the radio continuum emission separated by simplistically assuming a constant nonthermal spectral index. We use the data at the $Spitzer$ MIPS wavelengths of 24, 70, and 160 $\mu$m, as well as recent radio continuum map at 3.6 cm observed with the 100--m Effelsberg telescope. Using the wavelet transformation, we separate diffuse emission components from compact sources and study the radio-IR correlation at various scales. We also investigate the IR correlations with the thermal and nonthermal radio emissions separated by our developed method. A H$\alpha$ map serves as a tracer of star forming regions.

We determine the variation in the nonthermal radio spectral index in the nearby spiral galaxy M33 at a linear resolution of 360 pc. We separate the thermal and nonthermal components of the radio continuum emission without the assumption of a constant nonthermal spectral index. Using the Spitzer FIR data at 70 and 160 $\mu$m and a standard dust model, we deredden the H$\alpha$ emission. The extinction corrected H$\alpha$ emission serves as a template for the thermal free-free radio emission. Subtracting from the observed 3.6 cm and 20 cm emission (Effelsberg and the VLA) this free-free emission, we obtain the nonthermal maps. A constant electron temperature used to obtain the thermal radio intensity seems appropriate for M~33 which, unlike the Milky Way, has a shallow metallicity gradient. For the first time, we derive the distribution of the nonthermal spectral index across a galaxy, M33. We detect strong nonthermal emission from the spiral arms and star-forming regions. Wavelet analysis shows that at 3.6 cm the nonthermal emission is dominated by contributions from star-forming regions, while it is smoothly distributed at 20 cm. For the whole galaxy, we obtain thermal fractions of 51% and 18% at 3.6 cm and 20 cm, respectively. The thermal emission is slightly stronger in the southern than in the northern half of the galaxy. We find a clear radial gradient of mean extinction in the galactic plane. The nonthermal spectral index map indicates that the relativistic electrons suffer energy-loss when diffusing from their origin in star-forming regions towards interarm regions and the outer parts of the galaxy. We also conclude that the radio emission is mostly nonthermal at R $>$ 5 kpc in M33.

We study the exponential scale length of total radio emission, the spectral index distribution, and the linear radio polarization in the Scd galaxy M33. Observations were carried out using the 3.6 cm dual channel and the 6.2 cm four channel receivers of the 100--m Effelsberg telescope along with the L-band VLA D--array at 20 cm. High spatial resolution and sensitivity in both total and linearly polarized radio continuum emission from M33 were achieved. We found considerable extended emission, not only from the main arms I S and I N, but also from the weaker arms. The large--scale magnetic field exhibits well--ordered spiral structure with almost the same orientation as that of the optical spiral arms, however, it does not show a clear structural correlation or anti--correlation with the optical arms. There is a north-south asymmetry in polarization that is frequency-dependent. We found that the ring mean spectral index versus radius increases faster beyond $R$ = 4 kpc. At each wavelength, the exponential scale length is larger inside than outside $R$ = 4 kpc. From the larger scales lengths at $R$ $<$ 4 kpc, we conclude that star forming regions are mainly spread over the region $R$ $<$ 4 kpc without a dominant nuclear concentration. Furthermore, at $R$ $<$ 4 kpc, a spatial correlation between cosmic rays and star forming regions may exist. From the behaviour of the mean spectral indices obtained from different pairs of the radio continuum data at 3.6, 6.2, and 20 cm, we confirm that a decrease in the thermal fraction causes an increase in the spectral index. The frequency-dependent asymmetry in the polarization hints to an asymmetry in Faraday depolarization.

We investigate the energy sources of the infrared (IR) emission and their relation to the radio continuum emission at various spatial scales within the Scd galaxy M33. We use the wavelet transform to analyze IR data at the Spitzer wavelengths of 24, 70, and 160$\mu$m, as well as recent radio continuum data at 3.6cm and 20cm. An H$\alpha$ map serves as a tracer of the star forming regions and as an indicator of the thermal radio emission. We find that the dominant scale of the 70$\mu$m emission is larger than that of the 24$\mu$m emission, while the 160$\mu$m emission shows a smooth wavelet spectrum. The radio and H$\alpha$ maps are well correlated with all 3 MIPS maps, although their correlations with the 160$\mu$m map are weaker. After subtracting the bright HII regions, the 24 and 70$\mu$m maps show weaker correlations with the 20cm map than with the 3.6cm map at most scales. We also find a strong correlation between the 3.6cm and H$\alpha$ emission at all scales. Comparing the results with and without the bright HII regions, we conclude that the IR emission is influenced by young, massive stars increasingly with decreasing wavelength from 160 to 24$\mu$m. The radio-IR correlations indicate that the warm dust-thermal radio correlation is stronger than the cold dust-nonthermal radio correlation at scales smaller than 4kpc. A perfect 3.6cm-H$\alpha$ correlation implies that extinction has no significant effect on H$\alpha$ emitting structures.

Using LBQS data, it seems that the predictions of the model is confirmed by observations. The behavior of comoving density versus redshift is consistent with LDDE for QSQs. Considering the cosmic evolution of SFR, a unified aspect for origin and evolution of QSOs and ordinary galaxies is implied by this model.

The contraction model of Field and Colgate for proto-galaxies, first proposed to describe the observed properties of quasars, is generalized and used to investigate the evolution of galaxies. The LEDA data base for elliptical, spiral, compact and diffuse galaxies is employed and it is shown that the above model is consistent with observational evidences regarding their dynamical evolution, star formation rate and different morphologies.