In the second work of this series, we analyse the connection between the availability of gas and the position of a region with respect to the spatially resolved main sequence (MS) relation. Following the procedure presented in Paper I we obtain 500pc scales estimates of stellar mass and star formation rate surface densities ($\Sigma_{\star}$ and $\Sigma_{\rm{SFR}}$). Our sample consists of five face-on, grand design spiral galaxies located on the MS. Thanks to HI 21cm and $^{12}$CO(2-1) maps, we connect the gas surface densities and gas fractions to the observed star formation properties of each region. We find that the spatially resolved MS ($\sigma=0.23$ dex) is the combination of two relations: the Kennicutt-Schmidt law ($\sigma=0.19$ dex) and the molecular gas MS (MGMS, $\sigma=0.22$ dex); $\Sigma_{\star}$, $\Sigma_{\rm{SFR}}$ and the surface density of the molecular gas, $\Sigma_{\rm{H_2}}$, define a 3D relation as proposed by \citet{2019ApJ...884L..33L}. We find that $\Sigma_{\rm{H_2}}$ steadily increases along the MS relation, varies little towards higher $\Sigma_{\rm{SFR}}$ at fixed stellar surface densities (not enough to sustain the change in SFR), and it is almost constant perpendicular to the relation. The surface density of neutral gas ($\Sigma_{\rm{HI}}$) is constant along the MS, and increases in its upper envelop. $\Sigma_{\rm{SFR}}$ can be expressed as a function of $\Sigma_{\star}$ and $\Sigma_{\rm{HI}}$, following the Equation: $\log\Sigma_{\rm{SFR}}$ = 0.97$\log\Sigma_{\star}$ + 1.99$\log\Sigma_{\rm{HI}}$ - 11.11. Finally, we show that f$_{\rm{gas}}$ increases significantly towards the starburst region in the $\log\Sigma_{\star}$ - $\log\Sigma_{\rm{SFR}}$ plane, accompanied by a slight increase in SFE.

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.

The purpose of this work is the characterization of the main scaling relations between all the ISM components (dust, atomic/molecular/total gas), gas-phase metallicity, and other galaxy properties, such as Mstar and galaxy morphology, for late-type galaxies in the Local Universe. This study is performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage from T = 1 to 10, Mstar from 6 x 10^7 to 3 x 10^11 Msun, SFR from 6 x 10^(-4) to 60 Msun/yr, and 12 + log(O/H) from 8 to 9.5. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent XCO. The analysis has been performed by means of the survival analysis technique. We confirm that the dust mass correlates very well with the total gas mass, and find -- for the first time -- that the dust mass correlates better with the atomic gas mass than the molecular one. We characterize important mass ratios such as gas fraction, molecular-to-atomic gas mass ratio, dust-to-total gas mass ratio (DGR), and dust-to-stellar mass ratio. Only the assumption of a metallicity-dependent XCO reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. DGR, gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to the galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for molecular gas rich galaxies this trend breaks down. This trend has never been found previously, to our knowledge. It provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies.

Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made.

We compare the properties of clouds in simulated M33 galaxies to those observed in the real M33. We apply a friends of friends algorithm and CPROPS to identify clouds, as well as a pixel by pixel analysis. We obtain very good agreement between the number of clouds, and maximum mass of clouds. Both are lower than occurs for a Milky Way-type galaxy and thus are a function of the surface density, size and galactic potential of M33. We reproduce the observed dependence of molecular cloud properties on radius in the simulations, and find this is due to the variation in gas surface density with radius. The cloud spectra also show good agreement between the simulations and observations, but the exact slope and shape of the spectra depends on the algorithm used to find clouds, and the range of cloud masses included when fitting the slope. Properties such as cloud angular momentum, velocity dispersions and virial relation are also in good agreement between the simulations and observations, but do not necessarily distinguish between simulations of M33 and other galaxy simulations. Our results are not strongly dependent on the level of feedback used here (10 and 20%) although they suggest that 15% feedback efficiency may be optimal. Overall our results suggest that the molecular cloud properties are primarily dependent on the gas and mass surface density, and less dependent on the localised physics such as the details of stellar feedback, or the numerical code used.

We carried out deep searches for CO line emission in the outer disk of M33, at R>7 kpc, and examined the dynamical conditions that can explain variations in the mass distribution of the molecular cloud throughout the disk of M33. We used the IRAM-30~m telescope to search for CO lines in the outer disk toward 12 faint mid-infrared (MIR) selected sources and in an area of the southern outer disk hosting MA1, a bright HII region. We detect narrow CO lines at the location of two MIR sources at galactocentric distances of about 8 kpc that are associated with low-mass young stellar clusters, and at four locations in the proximity of MA1. The paucity of CO lines at the location of weak MIR-selected sources probably arises because most of them are not star-forming sites in M33, but background sources. Although very uncertain, the total molecular mass of the detected clouds around MA1 is lower than expected given the stellar mass of the cluster, because dispersal of the molecular gas is taking place as the HII region expands. The mean mass of the giant molecular clouds (GMCs) in M33 decreases radially by a factor 2 from the center out to 4 kpc, then it stays constant until it drops at R>7 kpc. We suggest that GMCs become more massive toward the center because of the fast rotation of the disk, which drives mass growth by coalescence of smaller condensations as they cross the arms. The analysis of both HI and CO spectral data gives the consistent result that corotation of the two main arms in this galaxy is at a radius of 4.7+-0.3 kpc, and spiral shock waves become subsonic beyond 3.9 kpc. Perturbations are quenched beyond 6.5 kpc, where CO lines have been detected only around sporadic condensations associated with UV and MIR emission.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

The probability distribution functions (PDFs) for atomic, molecular, and total gas surface densities of M33 are determined at a resolution of about 50~pc over regions that share coherent morphological properties to unveil fingerprints of self-gravity across the star-forming disk. Most of the total gas PDFs from the central region to the edge of the star-forming disk are well-fitted by log-normal functions whose width decreases radially outwards. Because the HI velocity dispersion is approximately constant across the disk, the decrease of the PDF width is consistent with a lower Mach number for the turbulent ISM at large galactocentric radii where a higher fraction of HI is in the warm phase. The atomic gas is found mostly at face-on column densities below N$_{H}^{lim}$=2.5 10$^{21}$~cm$^{-2}$, with small radial variations of N$_{H}^{lim}$. The molecular gas PDFs do not show strong deviations from log-normal functions in the central region where molecular fractions are high. Here the high pressure and rate of star formation shapes the PDF as a log-normal function dispersing self-gravitating complexes with intense feedback at all column densities that are spatially resolved. Power law PDFs for the molecules are found near and above N$_H^{lim}$, in the well defined southern spiral arm and in a continuous dense filament extending at larger galactocentric radii; this is evident in cloud samples at different evolutionary stages along the star formation cycle. In the filament nearly half of the molecular gas departs from a log-normal PDF and power laws are also observed in pre-star forming molecular complexes. The slope of the power law is between -1 and -2. This slope, combined with maps showing where the different parts of the power law PDFs come from, suggest a power-law stratification of density within molecular cloud complexes, which is consistent with the dominance of self-gravity.

We perform simulations of isolated galaxies in order to investigate the likely origin of the spiral structure in M33. In our models, we find that gravitational instabilities in the stars and gas are able to reproduce the observed spiral pattern and velocity field of M33, as seen in HI, and no interaction is required. We also find that the optimum models have high levels of stellar feedback which create large holes similar to those observed in M33, whilst lower levels of feedback tend to produce a large amount of small scale structure, and undisturbed long filaments of high surface density gas, hardly detected in the M33 disc. The gas component appears to have a significant role in producing the structure, so if there is little feedback, both the gas and stars organise into clear spiral arms, likely due to a lower combined $Q$ (using gas and stars), and the ready ability of cold gas to undergo spiral shocks. By contrast models with higher feedback have weaker spiral structure, especially in the stellar component, compared to grand design galaxies. We did not see a large difference in the behaviour of $Q_{stars}$ with most of these models, however, because $Q_{stars}$ stayed relatively constant unless the disc was more strongly unstable. Our models suggest that although the stars produce some underlying spiral structure, this is relatively weak, and the gas physics has a considerable role in producing the large scale structure of the ISM in flocculent spirals.

The sample of 566 molecular clouds identified in the CO(2--1) IRAM survey covering the disk of M~33 is explored in detail.The clouds were found using CPROPS and were subsequently catalogued in terms of their star-forming properties as non-star-forming (A), with embedded star formation (B), or with exposed star formation C.We find that the size-linewidth relation among the M~33 clouds is quite weak but, when comparing with clouds in other nearby galaxies, the linewidth scales with average metallicity.The linewidth and particularly the line brightness decrease with galactocentric distance.The large number of clouds makes it possible to calculate well-sampled cloud mass spectra and mass spectra of subsamples.As noted earlier, but considerably better defined here, the mass spectrum steepens (i.e. higher fraction of small clouds) with galactocentric distance.A new finding is that the mass spectrum of A clouds is much steeper than that of the star-forming clouds.Further dividing the sample, this difference is strong at both large and small galactocentric distances and the A vs C difference is a stronger effect than the inner/outer disk difference in mass spectra.Velocity gradients are identified in the clouds using standard techniques.The gradients are weak and are dominated by prograde rotation; the effect is stronger for the high signal-to-noise clouds.A discussion of the uncertainties is presented.The angular momenta are low but compatible with at least some simulations.The cloud and galactic gradients are similar; the cloud rotation periods are much longer than cloud lifetimes and comparable to the galactic rotation period.The rotational kinetic energy is 1-2\% of the gravitational potential energy and the cloud edge velocity is well below the escape velocity, such that cloud-scale rotation probably has little influence on the evolution of molecular clouds.

To shed light on the time evolution of local star formation episodes in M33, we study the association between 566 Giant Molecular Clouds (GMCs), identified through the CO (J=2-1) IRAM-all-disk survey, and 630 Young Stellar Cluster Candidates (YSCCs), selected via Spitzer-24~$\mu$m emission. The spatial correlation between YSCCs and GMCs is extremely strong, with a typical separation of 17~pc, less than half the CO(2--1) beamsize, illustrating the remarkable physical link between the two populations. GMCs and YSCCs follow the HI filaments, except in the outermost regions where the survey finds fewer GMCs than YSCCs, likely due to undetected, low CO-luminosity clouds. The GMCs have masses between 2$\times 10^4$ and 2$\times 10^6$ M$_\odot$ and are classified according to different cloud evolutionary stages: inactive clouds are 32$\%$ of the total, classified clouds with embedded and exposed star formation are 16$\%$ and 52$\%$ of the total respectively. Across the regular southern spiral arm, inactive clouds are preferentially located in the inner part of the arm, possibly suggesting a triggering of star formation as the cloud crosses the arm. Some YSCCs are embedded star-forming sites while the majority have GALEX-UV and H$\alpha$ counterparts with estimated cluster masses and ages. The distribution of the non-embedded YSCC ages peaks around 5~Myrs with only a few being as old as 8--10~Myrs. These age estimates together with the number of GMCs in the various evolutionary stages lead us to conclude that 14~Myrs is a typical lifetime of a GMC in M33, prior to cloud dispersal. The inactive and embedded phases are short, lasting about 4 and 2~Myrs respectively. This underlines that embedded YSCCs rapidly break out from the clouds and become partially visible in H$\alpha$ or UV long before cloud dispersal.

The stellar and gaseous mass distributions, as well as the extended rotation curve in the nearby galaxy M33 are used to derive the radial distribution of dark matter density in the halo and to test cosmological models of galaxy formation and evolution. Two methods are examined to constrain dark mass density profiles. The first method deals directly with fitting the rotation curve data in the range of galactocentric distances $0.24\,\text{kpc}\leq r\leq22.72\,\text{kpc}.$ As found in a previous paper by \citet{Corbelli:2014lga}, and using the results of collisionless $\Lambda-$Cold Dark Matter numerical simulations, we confirm that the Navarro-Frenkel-White (hereafter NFW) dark matter profile provides a better fit to the rotation curve data than the cored Burkert profile (hereafter BRK) profile. The second method relies on the local equation of centrifugal equilibrium and on the rotation curve slope. In the aforementioned range of distances we fit the observed velocity profile, using a function which has a rational dependence on the radius, and derive the slope of the rotation curve. Following \citet{Salucci:2010qr} we then infer the effective matter densities. In the radial range $9.53\,\text{kpc}\leq r\leq22.72\,\text{kpc}$ the uncertainties induced by the luminous matter (stars and gas) becomes negligible, because the dark matter density dominates, and we can determine locally the radial distribution of dark matter. With this second method we tested the NFW and the BRK dark matter profiles and confirm that both profiles are compatible with the data even though in this case the cored BRK density profile provides a more reasonable value for the baryonic-to-dark matter ratio.

We present $^{12}$CO(1-0) and $^{12}$CO(2-1) observations of a sample of 20 star-forming dwarfs selected from the Herschel Virgo Cluster Survey, with oxygen abundances ranging from 12 + log(O/H) ~ 8.1 to 8.8. CO emission is observed in ten galaxies and marginally detected in another one. CO fluxes correlate with the FIR 250 $\mu$m emission, and the dwarfs follow the same linear relation that holds for more massive spiral galaxies extended to a wider dynamical range. We compare different methods to estimate H2 molecular masses, namely a metallicity-dependent CO-to-H2 conversion factor and one dependent on H-band luminosity. The molecular-to-stellar mass ratio remains nearly constant at stellar masses <~ 10$^9$ M$_{\odot}$, contrary to the atomic hydrogen fraction, M$_{HI}$/M$_*$, which increases inversely with M$_*$. The flattening of the M$_{H_2}$/M$_*$ ratio at low stellar masses does not seem to be related to the effects of the cluster environment because it occurs for both HI-deficient and HI-normal dwarfs. The molecular-to-atomic ratio is more tightly correlated with stellar surface density than metallicity, confirming that the interstellar gas pressure plays a key role in determining the balance between the two gaseous components of the interstellar medium. Virgo dwarfs follow the same linear trend between molecular gas mass and star formation rate as more massive spirals, but gas depletion timescales, $\tau_{dep}$, are not constant and range between 100 Myr and 6 Gyr. The interaction with the Virgo cluster environment is removing the atomic gas and dust components of the dwarfs, but the molecular gas appears to be less affected at the current stage of evolution within the cluster. However, the correlation between HI deficiency and the molecular gas depletion time suggests that the lack of gas replenishment from the outer regions of the disc is lowering the star formation activity.

To assess the effects of the cluster environment on the different components of the interstellar medium, we analyse the FIR-submm properties of a sample of star-forming dwarf (SFD) galaxies detected by the Herschel Virgo Cluster Survey (HeViCS). We determine dust masses and dust temperatures by fitting a modified black body (MBB) function to the spectral energy distributions (SEDs). Stellar and gas masses, star formation rates (SFRs), and metallicities are obtained from the analysis of a set of ancillary data. Dust is detected in 49 out of 140 optically identified dwarfs covered by the HeViCS field; considering only dwarfs brighter than $m_B$ = 18 mag, this gives a detection rate of 43%. After evaluating different emissivity indices, we find that the FIR-submm SEDs are best-fit by $\beta$=1.5, with a median dust temperature $T_d$ = 22.4 K. Assuming $\beta$=1.5, 67% of the 23 galaxies detected in all five Herschel bands show emission at 500 $\mu$m in excess of the MBB model. The excess is inversely correlated with SFR and stellar masses. To study the variations in the global properties of our sample due to environmental effects, we compare the Virgo SFDs to other Herschel surveys, such as KINGFISH, the Dwarf Galaxy Survey (DGS), and the HeViCS bright galaxy catalogue (BGC). We explore the relations between stellar mass and HI fraction, specific SFR, dust fraction, gas-to-dust ratio over a wide range of stellar masses. Highly HI-deficient Virgo dwarf galaxies are mostly characterised by quenched star formation activity and lower dust fractions giving hints for dust stripping in cluster dwarfs. However, we find that the fraction of dust removed has to be less than that of the HI component. Since the Virgo SFDs are likely to be crossing the cluster for the first time, a longer timescale might be necessary to strip the more centrally concentrated dust distribution.

We present three independent catalogs of point-sources extracted from SPIRE images at 250, 350, and 500 micron as a part of the Herschel Virgo Cluster Survey (HeViCS). The source positions are determined by estimating the likelihood to be a real source for each peak on the maps and the flux densities are estimated using the sourceExtractorTimeline, a timeline-based point source fitter. Afterwards, each source is subtracted from the maps, removing a Gaussian function in every position with the full width half maximum equal to that estimated in sourceExtractorTimeline. This procedure improves the robustness of our algorithm in terms of source identification. The HeViCS catalogs contain about 52000, 42200, and 18700 sources selected at 250, 350, and 500 micron above 3sigma and are ~ 75%, 62%, and 50% complete at flux densities of 20 mJy at 250, 350, 500 micron, respectively. We then measured source number counts at 250, 350, and 500 micron and we also cross-correlated the catalogs with the Sloan Digital Sky Survey to investigate the redshift distribution of the nearby sources. From this cross-correlation, we select ~2000 sources with reliable fluxes and a high signal-to-noise ratio, finding an average redshift z~0.3+/-0.22. The number counts at 250, 350, and 500 micron show an increase in the slope below 200 mJy, indicating a strong evolution in number of density for galaxies at these fluxes. In general, models tend to overpredict the counts at brighter flux densities, underlying the importance of studying the Rayleigh-Jeans part of the spectral energy distribution to refine the theoretical recipes of the models. Our iterative method for source identification allowed the detection of a family of 500 micron sources that are not foreground objects belonging to Virgo and not found in other catalogs.

We map the neutral atomic gas content of M33 using high resolution VLA and GBT observations and fit a tilted ring model to determine the orientation of the extended gaseous disk and its rotation curve. The disk of M33 warps from 8 kpc outwards without substantial change of its inclination with respect to the line of sight. Rotational velocities rise steeply with radius in the inner disk, reaching 100 km/s in 4 kpc, then the rotation curve becomes more perturbed and flatter with velocities as high as 120-130 km/s out to 23 kpc. We derive the stellar mass surface density map of M33's optical disk, via pixel -SED fitting methods based on population synthesis models, which highlights variations in the mass-to-light ratio. The stellar mass surface further out is estimated from deep images of outer disk fields. Stellar and gas maps are then used in the dynamical analysis of the rotation curve to constrain the dark matter distribution which is relevant at all radii. A dark matter halo with a Navarro-Frenk-White density profile in a LCDM cosmology, provides the best fit to the rotation curve for a dark halo concentration C=10 and a total halo mass of 4.3 10^{11}Msun. This imples a baryonic fraction of order 0.02 and the evolutionary history of this galaxy should account for loss of a large fraction of its original baryonic content.

In order to study the ISM and the interplay between the atomic and molecular components in a low-metallicity environment, we present a complete high angular and spectral resolution map and data cube of the 12CO(2-1) emission from the Local Group galaxy M33. Its metallicity is roughly half-solar, such that we can compare its ISM with that of the Milky Way with the main changes being the metallicity and the gas mass fraction. The data have a 12" angular resolution (50pc) with a spectral resolution of 2.6 km/s and a mean noise level of 20 mK per channel in antenna temperature. A radial cut along the major axis was also observed in the 12CO(1-0) line. The CO data cube and integrated intensity map are optimal when using HI data to define the baseline window and the velocities over which the CO emission is integrated. Great care was taken when building these maps, testing different windowing and baseline options and investigating the effect of error beam pickup. The total CO(2-1) luminosity is 2.8e7 K km/s pc2, following the spiral arms in the inner disk. There is no clear variation in the CO(2-1/1-0) intensity ratio with radius and the average value is roughly 0.8. The total molecular gas mass is estimated, using a N(H2)/Ico(1-0)=4e20 cm-2/(K km/s) conversion factor, to be 3.1e8 Msol. The CO spectra in the cube were shifted to zero velocity by subtracting the velocity of the HI peak from the CO spectra. Hence, the velocity dispersion between the atomic and molecular components is extremely low, independently justifying the use of the HI line in building our maps. Stacking the spectra in concentric rings shows that the CO linewidth and possibly the CO-HI velocity dispersion decrease in the outer disk. Using the CO(2-1) emission to trace the molecular gas, the probability distribution function of the H2 column density shows an excess at high column density above a log normal distribution.

Aims. We study the dust content of a large optical input sample of 910 early-type galaxies (ETG) in the Virgo cluster, extending also to the dwarf ETG, and examine the results in relation with those on the other cold ISM components. Methods. We searched for far-infrared emission in all galaxies of the input sample using the 250 micron image of the Herschel Virgo Cluster Survey (HeViCS). This image covers a large fraction of the cluster. For the detected ETG we measured fluxes in 5 bands from 100 to 500 micron, and estimated the dust mass and temperature with modified black-body fits. Results. Dust is detected above the completeness limit of 25.4 mJy at 250 micron in 46 ETG, 43 of which are in the optically complete part of the input sample. In addition dust is present at fainter levels in another 6 ETG. We detect dust in the 4 ETG with synchrotron emission, including M 87. Dust appears to be much more concentrated than stars and more luminous ETG have higher dust temperatures. Dust detection rates down to the 25.4 mJy limit are 17% for ellipticals, about 40% for lenticulars (S0 + S0a) and around 3% for dwarf ETG. Dust mass does not correlate clearly with stellar mass and is often much more than that expected for a passive galaxy in a closed-box model. The dust-to-stars mass ratio anticorrelates with galaxy luminosity, and for some dwarf ETG reaches values as high as for dusty late-type galaxies. In the Virgo cluster slow rotators appear more likely to contain dust than fast ones. Comparing the dust results with those on HI from ALFALFA, there are only 8 ETG detected both in dust and in HI in the HeViCS area; 39 have dust but only an upper limit on HI, and 8 have HI but only an upper limit on dust. The locations of these galaxies in the cluster are different, with the dusty ETG concentrated in the densest regions, while the HI rich ETG are at the periphery.

The Herschel Virgo Cluster Survey (HeViCS) is the deepest, confusion-limited survey of the Virgo Cluster at far-infrared (FIR) wavelengths. The entire survey at full depth covers $\sim$55 sq. deg. in 5 bands (100-500 \micron), encompassing the areas around the central dominant elliptical galaxies (M87, M86 & M49) and extends as far as the NW cloud, the W cloud and the Southern extension. The survey extends beyond this region with lower sensitivity so that the total area covered is 84 sq. deg. In this paper we describe the data, the data acquisition techniques and present the detection rates of the optically selected Virgo Cluster Catalogue (VCC). We detect 254 (34%) of 750 VCC galaxies found within the survey boundary in at least one band and 171 galaxies are detected in all five bands. For the remainder of the galaxies we have measured strict upper limits for their FIR emission. The population of detected galaxies contains early- as well as late-types although the latter dominate the detection statistics. We have modelled 168 galaxies, showing no evidence of a strong synchrotron component in their FIR spectra, using a single-temperature modified blackbody spectrum with a fixed emissivity index ($\beta = 2$). A study of the $\chi^2$ distribution indicates that this model is not appropriate in all cases, and this is supported by the FIR colours which indicate a spread in $\beta$=1--2. Statistical comparison of the dust mass and temperature distributions from 140 galaxies with $\chi^2_{dof=3} < 7.8$ (95% confidence level) shows that late-types have typically colder, more massive dust reservoirs; the early-type dust masses have a mean of ${\rm log}(<M> / M_{\sun}) = 6.3 \pm 0.3 $, while for late-types ${\rm log}(<M> / M_{\sun}) =7.1 \pm 0.1$... (abridged)

We investigate the dust-to-gas mass ratio and the environmental effects on the various components of the interstellar medium for a spatially resolved sample of Virgo spirals. We have used the IRAM-30m telescope to map over their full extent NGC 4189, NGC 4298, NGC 4388, and NGC 4299 in the 12CO(1-0) and the 12CO(2-1) lines. We observed the same lines in selected regions of NGC 4351, NGC 4294, and NGC 4424. The CO observations are combined with Herschel maps in 5 bands between 100-500 {\mu}m from the HeViCS survey, and with HI data from the VIVA survey, to obtain spatially resolved dust and gas distributions. We studied the environmental dependencies by adding to our sample eight galaxies with 12CO(1-0) maps from the literature. We estimate the integrated mass of molecular hydrogen for the galaxies observed in the CO lines. We find molecular-to-total gas mass fractions between 0.04 \leq fmol \leq 0.65, with the lowest values for the dimmest galaxy in the B-band. The integrated dust-to-gas ratio ranges between 0.011 and 0.004. For the 12 mapped galaxies we derive the radial distributions of the atomic gas, molecular gas, and dust. We also study the effect of different CO-to-H2 conversion factors. Both the molecular gas and the dust distributions show steeper radial profiles for HI-deficient galaxies and the average dust-to-gas ratio for these galaxies increases or stays radially constant. On scales of \sim 3 kpc, we find a strong correlation between the molecular gas and the 250 micron surface brightness that is tighter than average for non-deficient galaxies. The correlation becomes linear if we consider the total gas surface mass density. However, the inclusion of atomic hydrogen does not improve the statistical significance of the correlation. The environment can modify the distributions of molecules and dust within a galaxy, although these components are more tightly bound than the atomic gas.

Using the far-infrared emission, as observed by the Herschel Virgo Cluster Survey (HeViCS), and the integrated HI and CO brightness, we infer the dust and total gas mass for a magnitude limited sample of 35 metal rich spiral galaxies in Virgo. The CO flux correlates tightly and linearly with far-infrared fluxes observed by Herschel. Molecules in these galaxies are more closely related to cold dust rather than to dust heated by star formation or to optical/NIR brightness. We show that dust mass establishes a stronger correlation with the total gas mass than with the atomic or molecular component alone. The dust-to-gas ratio increases as the HI deficiency increases, but in highly HI deficient galaxies it stays constant. Dust is in fact less affected than atomic gas by weak cluster interactions, which remove most of the HI gas from outer and high latitudes regions. Highly disturbed galaxies, in a dense cluster environment, can instead loose a considerable fraction of gas and dust from the inner regions of the disk keeping constant the dust-to-gas ratio. There is evidence that the molecular phase is also quenched. This quencing becomes evident by considering the molecular gas mass per unit stellar mass. Its amplitude, if confirmed by future studies, highlights that molecules are missing in Virgo HI deficient spirals, but to a somewhat lesser extent than dust.

We combine new Herschel/SPIRE sub-millimeter observations with existing multiwavelength data to investigate the dust scaling relations of the Herschel Reference Survey, a magnitude-, volume-limited sample of ~300 nearby galaxies in different environments. We show that the dust-to-stellar mass ratio anti-correlates with stellar mass, stellar mass surface density and NUV-r colour across the whole range of parameters covered by our sample. Moreover, the dust-to-stellar mass ratio decreases significantly when moving from late- to early-type galaxies. These scaling relations are similar to those observed for the HI gas-fraction, supporting the idea that the cold dust is tightly coupled to the cold atomic gas component in the interstellar medium. We also find a weak increase of the dust-to-HI mass ratio with stellar mass and colour but no trend is seen with stellar mass surface density. By comparing galaxies in different environments we show that, although these scaling relations are followed by both cluster and field galaxies, HI-deficient systems have, at fixed stellar mass, stellar mass surface density and morphological type systematically lower dust-to-stellar mass and higher dust-to-HI mass ratios than HI-normal/field galaxies. This provides clear evidence that dust is removed from the star-forming disk of cluster galaxies but the effect of the environment is less strong than what is observed in the case of the HI disk. Such effects naturally arise if the dust disk is less extended than the HI and follows more closely the distribution of the molecular gas phase, i.e., if the dust-to-atomic gas ratio monotonically decreases with distance from the galactic center.

We present an analysis of the systematic CO(2-1) survey at 12" resolution covering most of the local group spiral M 33 which, at a distance of 840 kpc, is close enough that individual giant molecular clouds (GMCs) can be identified. The goal of this work is to study the properties of the GMCs in this subsolar metallicity galaxy. The CPROPS (Cloud Properties) algorithm (Rosolowsky & Leroy 2006) was used to identify 337 GMCs in M 33, the largest sample to date in an external galaxy. The sample is used to study the GMC luminosity function, or mass spectrum under the assumption of a constant N(H2)/ICO ratio. We find that n(L)dL = K*L^(-2.0\pm0.1) for the entire sample. However, when the sample is divided into inner and outer disk samples, the exponent changes from 1.6 \pm 0.2 for the centre 2 kpc to 2.3 \pm 0.2 for galactocentric distances larger than 2 kpc. Based on the emission in the FUV, Halpha, 8mu, and 24mu bands, each cloud was classified in terms of its star forming activity - no star formation, embedded, or exposed star formation (visible in FUV and Halpha). At least one sixth of the clouds had no (massive) star formation, suggesting that the average time required for star formation to start is about one sixth of the total time for which the object is identifiable as a GMC. The clouds without star formation have significantly lower CO luminosities than those with star formation, whether embedded or exposed, presumably related to the lack of heating sources. Taking the cloud sample as a whole, the main non-trivial correlation was the decrease in cloud CO brightness (or luminosity) with galactocentric radius. The complete cloud catalog, including CO and HI spectra and the CO contours on the FUV, Halpha, 8mu, and 24mu images is presented in the appendix.

We describe the Herschel Virgo Cluster Survey (HeViCS) and the first data that cover the complete survey area (four 4 x 4 deg2 regions). We use these data to measure and compare the global far infrared properties of 78 optically bright galaxies that are selected at 500 \mum and detected in all five far-infrared bands. We show that our measurements and calibration are broadly consistent with previous data obtained by IRAS, ISO, Spitzer and Planck. We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 \mum cluster luminosity distributions. These luminosity distributions are not power laws, but peaked, with small numbers of both faint and bright galaxies. We measure a cluster 100-500 micron far-infrared luminosity density of 1.6(7.0) \pm 0.2 x 10^9 Lsun/Mpc3. This compares to a cluster 0.4-2.5 \mum optical luminosity density of 5.0(20.0) x 10^9 Lsun/Mpc3, some 3.2(2.9) times larger than the far-infrared. A typical photon originates from an optical depth of 0.4\pm0.1. Most of our sample galaxies are well fitted by a single modified blackbody (beta=2), leading to a mean dust mass of log Mdust = 7.31 Msun and temperature of 20.0 K. We also derive both stellar and atomic hydrogen masses from which we calculate mean values for the stars:gas(atomic) and gas(atomic): dust mass ratios of 15.1 and 58.2 respectively. Using our derived dust, atomic gas and stellar masses we estimate cluster mass densities of 8.6(27.8) x 10^6, 4.6(13.9) x 10^8, 7.8(29.7) x 10^9 Msun/Mpc3, respectively for dust, atomic gas and stars. These values are higher than those derived for field galaxies by factors of 39(126), 6(18) and 34(129) respectively. In the above luminosity/mass densities are given using the whole sample with values in brackets using just those galaxies that lie between 17 and 23 Mpc. We provide a data table of flux densities in all the Herschel bands for all 78 bright Virgo cluster galaxies.

The properties of young stellar clusters (YSCs) in M33, identified from the center out to about twice the size of the bright star-forming disk,are investigated. We find 915 discrete MIR sources as far as the extent of the warped HI disk, i.e. 16 kpc from the galaxy center. Their surface density has a steep radial decline beyond 4.5 kpc, and flattens out beyond the optical radius at 8.5 kpc. We are able to identify YSCs out to 12 kpc. At large galactocentric radii, the paucity of very luminous clusters and the relevance of hot dust emission become evident from the analysis of the bolometric and MIR luminosity functions. The YSC mass and size are correlated with a log-log slope of 2.09, similar to that measured for giant molecular clouds in M33 and the Milky Way, which represent the protocluster environment. Most of the YSCs in our sample have low extinction and ages between 3 and 10 Myr. In the inner regions of M33 the clusters span a wide range of mass (10^2<M<3 10^5 msun) and luminosity 10^38<L{bol}<3 10^{41}erg/s, while at galactocentric radii larger than 4 kpc we find a deficiency of massive clusters. Beyond 7 kpc, where the Halpha surface brightness drops significantly, the dominant YSC population has M<10^3 msun and a slightly older age (10 Myrs). This implies the occurrence of star formation events about 10 Myr ago as far as 10-12 kpc from the center of M33. The cluster L{FUV}--L{Halpha} relation is non-linear for L{FUV}<10^{39}erg/s, in agreement with randomly sampled models of the IMF which, furthermore, shows no appreciable variation throughout the M33 disk.

We present Subaru/Suprime-Cam deep V and I imaging of seven fields in the outer regions of M33. Our aim is to search for stellar structures corresponding to extended HI clouds found in a recent 21-cm survey of the galaxy. Three fields probe a large HI complex to the southeastern (SE) side of the galaxy. An additional three fields cover the northwestern (NW) side of the galaxy along the HI warp. A final target field was chosen further north, at a projected distance of approximately 25 kpc, to study part of the large stellar plume recently discovered around M33. We analyse the stellar population at R > 10 kpc by means of V, I colour magnitude diagrams reaching the red clump. Evolved stellar populations are found in all fields out to 120' (~ 30 kpc), while a diffuse population of young stars (~ 200 Myr) is detected out to a galactocentric radius of 15 kpc. The mean metallicity in the southern fields remains approximately constant at [M/H] = -0.7 beyond the edge of the optical disc, from 40' out to 80'. Along the northern fields probing the outer \hi disc, we also find a metallicity of [M/H] = -0.7 between 35' and 70' from the centre, which decreases to [M/H] = -1.0 at larger angular radii out to 120'. In the northernmost field, outside the disc extent, the stellar population of the large stellar feature possibly related to a M33-M31 interaction is on average more metal-poor ([M/H] = -1.3) and older (> 6 Gyr). An exponential disc with a large scale-length (~ 7 kpc) fits well the average distribution of stars detected in both the SE and NW regions from a galactocentric distance of 11 kpc out to 30 kpc. The stellar distribution at large radii is disturbed and, although there is no clear correlation between the stellar substructures and the location of the HI clouds, this gives evidence for tidal interaction or accretion events.

We analyse multi-wavelength observations of 32 young star clusters and associations in M33 with known oxygen abundance (8 < 12 + log(O/H) < 8.7), using ultraviolet (UV), optical, mid-infrared (MIR), CO (1-0) and 21-cm line (HI) observations. We derive their spectral energy distribution, and we determine age, bolometric luminosities, masses and the extinction, by comparing the multi-band integrated photometry to single-age stellar population models. The stellar system ages range between 2 and 15 Myr, masses are between 3 x 10^2 and 4 x 10^4 M_sun, and the intrinsic extinction, A_V, varies from 0.3 to 1 mag. We find a correlation between age and extinction, and between the cluster mass and size. The MIR emission shows the presence of a dust component around the clusters whose fractional luminosity at 24 um, L_{24}/L_{Bol}, decreases with the galactocentric distance. However, the total IR luminosity inferred from L_{24} is smaller than what we derive from the extinction corrections. The Halpha luminosity predicted by population synthesis models is larger than the observed one, especially for low-mass systems (M < 10^4 M_sun). Such a difference is reduced, but not erased, when the incomplete sampling of the initial mass function (IMF) at the high-mass end is taken into account. Our results suggest that a non-negligible fraction of UV ionising and non-ionising radiation is leaking into the ISM outside the HII regions. This would be in agreement with the large UV and Halpha diffuse fractions observed in M33, but it implies that stellar systems younger than 3 Myr retain, on average, only 30% of their Lyman continuum photons. However, the uncertainties on cluster ages and the stochastic fluctuations of the IMF do not allow to accurately quantify this issue. We also consider the possibility that this discrepancy is the consequence of a suppressed or delayed formation of the most massive stars.

Using Herschel data from the Open Time Key Project the Herschel Virgo Cluster Survey (HeViCS), we investigated the relationship between the metallicity gradients expressed by metal abundances in the gas phase as traced by the chemical composition of HII regions, and in the solid phase, as traced by the dust-to-gas mass ratio. We derived the radial gradient of the dust-to-gas mass ratio for all galaxies observed by HeViCS whose metallicity gradients are available in the literature. They are all late type Sbc galaxies, namely NGC4254, NGC4303, NGC4321, and NGC4501. We examined different dependencies on metallicity of the CO-to-H$_2$ conversion factor (\xco), used to transform the $^{12}$CO observations into the amount of molecular hydrogen. We found that in these galaxies the dust-to-gas mass ratio radial profile is extremely sensitive to choice of the \xco\ value, since the molecular gas is the dominant component in the inner parts. We found that for three galaxies of our sample, namely NGC4254, NGC4321, and NGC4501, the slopes of the oxygen and of the dust-to-gas radial gradients agree up to $\sim$0.6-0.7R$_{25}$ using \xco\ values in the range 1/3-1/2 Galactic \xco. For NGC4303 a lower value of \xco$\sim0.1\times$ 10$^{20}$ is necessary. We suggest that such low \xco\ values might be due to a metallicity dependence of \xco (from close to linear for NGC4254, NGC4321, and NGC4501 to superlinear for NGC4303), especially in the radial regions R$_G<$0.6-0.7R$_{25}$ where the molecular gas dominates. On the other hand, the outer regions, where the atomic gas component is dominant, are less affected by the choice of \xco, and thus we cannot put constraints on its value.

We estimate the conversion factor relating CO emission to H2 mass, alpha_CO, in five Local Group galaxies that span approximately an order of magnitude in metallicity - M31, M 33, the Large Magellanic Cloud (LMC), NGC 6822, and the Small Magellanic Cloud (SMC). We model the dust mass along the line of sight from infrared (IR) emission and then solve for the alpha_CO that best allows a single gas-to-dust ratio (delta_GDR) to describe each system. This approach remains sensitive to CO-dark envelopes of H2 surrounding molecular clouds. In M 31, M 33, and the LMC we find alpha_CO \approx 3-9 M_sun pc^-2 (K km s^-1)^-1, consistent with the Milky Way value within the uncertainties. The two lowest metallicity galaxies in our sample, NGC 6822 and the SMC (12 + log(O/H) \approx 8.2 and 8.0), exhibit a much higher alpha_CO. Our best estimates are \alpha_NGC6822 \approx 30 M_sun/pc^-2 (K km s^-1)^-1 and \alpha_SMC \approx 70 M_sun/pc^-2 (K km s-1)-1. These results are consistent with the conversion factor becoming CO a strong function of metallicity around 12 + log(O/H) \sim 8.4 - 8.2. We favor an interpretation where decreased dust-shielding leads to the dominance of CO-free envelopes around molecular clouds below this metallicity.

We investigate the nature of 24micron sources in M33 which have weak or no associated Halpha emission. Both bright evolved stars and embedded star forming regions are visible as compact infrared sources in the 8 and 24micron maps of M33 and contribute to the more diffuse and faint emission in these bands. Can we distinguish the two populations? We carry out deep CO J=2-1 and J=1-0 line searches at the location of compact mid-IR sources to unveil an ongoing star formation process. We use different assumptions to estimate cloud masses from pointed observations and analyze if SED and mid-IR colours can be used to discriminate between evolved stars and star forming regions. Molecular emission has been detected at the location of several sources at the level of 0.3 K km/s or higher in at least one of the CO rotational lines. Even though there are no giant molecular clouds beyond 4kpc in M33, our deep observations have revealed that clouds of smaller mass are very common. Sources which are known to be evolved variable stars show weaker or undetectable CO lines. Evolved stars occupy a well defined region of the IRAC color-color diagrams. Star forming regions are scattered throughout a larger area even though the bulk of the distribution has different IRAC colors than evolved variable stars. We estimate that about half of the 24 micron sources without an Halpha counterpart are genuine embedded star forming regions. Sources with faint but compact Halpha emission have an incomplete Initial Mass Function (IMF) at the high-mass end and are compatible with a population of young clusters with a stochastically sampled, universal IMF.

We present a new method to analyze the IMF at its high mass end in young stellar clusters, which rely on two integrated observables: the cluster bolometric and Halpha luminosity. Using several cluster samples selected in M33 we show that a stochastically sampled universal IMF is in better agreement with the data than a truncated IMF whose maximum stellar mass depends on cluster mass. We also discuss the possibility that a delayed formation of massive stars is taking place in low density star forming regions as an alternative to a strong leakage of ionizing photons from HII regions of young luminous clusters.

The origin of the far-infrared emission from the nearby radio galaxy M87 remains a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust emission, whereas others propose that the far-infrared emission can be explained by synchrotron emission without the need for an additional dust emission component. We observed M87 with PACS and SPIRE as part of the Herschel Virgo Cluster Survey (HeViCS). We compare the new Herschel data with a synchrotron model based on infrared, submm and radio data to investigate the origin of the far-infrared emission. We find that both the integrated SED and the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds with previous claims, we find no evidence of a diffuse dust component in M87.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

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.

Passive early-type galaxies (ETGs) provide an ideal laboratory for studying the interplay between dust formation around evolved stars and its subsequent destruction in a hot gas. Using Spitzer-IRS and Herschel data we compare the dust production rate in the envelopes of evolved AGB stars with a constraint on the total dust mass. Early-type galaxies which appear to be truly passively evolving are not detected by Herschel. We thus derive a distance independent upper limit to the dust grain survival time in the hostile environment of ETGs of < 46 +/- 25 Myr for amorphous silicate grains. This implies that ETGs which are detected at far-infrared wavelengths have acquired a cool dusty medium via interaction. Given likely time-scales for ram-pressure stripping, this also implies that only galaxies with dust in a cool (atomic) medium can release dust into the intra-cluster medium.

We use the Science Demonstration Phase data of the Herschel Virgo Cluster Survey to search for dust emission of early-type dwarf galaxies in the central regions of the Virgo Cluster as an alternative way of identifying the interstellar medium.We present the first possible far-infrared detection of cluster early-type dwarf galaxies: VCC781 and VCC951 are detected at the 10 sigma level in the SPIRE 250 micron image. Both detected galaxies have dust masses of the order of 10^5 Msun and average dust temperatures ~20K. The detection rate (less than 1%) is quite high compared to the 1.7% detection rate for Hi emission, considering that dwarfs in the central regions are more Hi deficient. We conclude that the removal of interstellar dust from dwarf galaxies resulting from ram pressure stripping, harassment, or tidal effects must be as efficient as the removal of interstellar gas.

The origin of the far-infrared emission from the nearby radio galaxy M87 remains a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust emission, whereas others propose that the far-infrared emission can be explained by synchrotron emission without the need for an additional dust emission component. We present Herschel PACS and SPIRE observations of M87, taken as part of the science demonstration phase observations of the Herschel Virgo Cluster Survey. We compare these data with a synchrotron model based on mid-infrared, far-infrared, submm and radio data from the literature to investigate the origin of the far-infrared emission. Both the integrated SED and the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds with previous claims, we find no evidence of a diffuse dust component in M87, which is not unexpected in the harsh X-ray environment of this radio galaxy sitting at the core of the Virgo Cluster.

We present the dust properties of a small sample of Virgo cluster dwarf galaxies drawn from the science demonstration phase data set of the Herschel Virgo Cluster Survey. These galaxies have low metallicities (7.8 < 12 + log(O/H) < 8.3) and star-formation rates < 10^{-1} M_{sun}/yr. We measure the spectral energy distribution (SED) from 100 to 500 um and derive dust temperatures and dust masses. The SEDs are fitted by a cool component of temperature T < 20 K, implying dust masses around 10^{5} M_{sun} and dust-to-gas ratios D within the range 10^{-3}-10^{-2}. The completion of the full survey will yield a larger set of galaxies, which will provide more stringent constraints on the dust content of star-forming dwarf galaxies.

We present a resolved dust analysis of three of the largest angular size spiral galaxies, NGC 4501 and NGC 4567/8, in the Herschel Virgo Cluster Survey (HeViCS) Science Demonstration field. Herschel has unprecedented spatial resolution at far-infrared wavelengths and with the PACS and SPIRE instruments samples both sides of the peak in the far infrared spectral energy distribution (SED).We present maps of dust temperature, dust mass, and gas-to-dust ratio, produced by fitting modified black bodies to the SED for each pixel. We find that the distribution of dust temperature in both systems is in the range ~19 - 22 K and peaks away from the centres of the galaxies. The distribution of dust mass in both systems is symmetrical and exhibits a single peak coincident with the galaxy centres. This Letter provides a first insight into the future analysis possible with a large sample of resolved galaxies to be observed by Herschel.

We describe the Herschel Virgo Cluster Survey (HeViCS) and the first data obtained as part of the Science Demonstration Phase (SDP). The data cover a central 4x4 sq deg region of the cluster. We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 micron luminosity functions (LFs) for optically bright galaxies that are selected at 500 micron and detected in all bands. We compare these LFs with those previously derived using IRAS, BLAST and Herschel-ATLAS data. The Virgo Cluster LFs do not have the large numbers of faint galaxies or examples of very luminous galaxies seen previously in surveys covering less dense environments.

By combining Herschel-SPIRE observations obtained as part of the Herschel Virgo Cluster Survey with 21 cm HI data from the literature, we investigate the role of the cluster environment on the dust content of Virgo spiral galaxies.We show for the first time that the extent of the dust disk is significantly reduced in HI-deficient galaxies, following remarkably well the observed 'truncation' of the HI disk. The ratio of the submillimetre-to- optical diameter correlates with the HI-deficiency, suggesting that the cluster environment is able to strip dust as well as gas. These results provide important insights not only into the evolution of cluster galaxies but also into the metal enrichment of the intra-cluster medium.

We present infrared colours (in the 25-500 mic spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is dominated bynthe synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-S0a than in star-forming systems probably because of the different nature of their dust heating sources (evolved stellar populations, X-ray, fast electrons) and dust grain properties. In contrast to the colour temperature of the warm dust, the f350/f500 index sensitive to the cold dust decreases with star formation and increases with metallicity, suggesting an overabundance of cold dust or an emissivity parameter beta<2 in low metallicity, active systems.

Nearby galaxies are ideal places to study in detail metallicity gradients and their time evolution. We consider chemical abundances of a new sample of \hii\ regions complemented with previous literature data-sets. We compare \hii\ region and PN abundances obtained with a common set of observations taken at MMT. With an updated theoretical model, we follow the time evolution of the baryonic components and chemical abundances in the disk of M33, assuming that the galaxy is accreting gas from an external reservoir. Supported by a uniform sample of nebular spectroscopic observations, we conclude that: {\em i}) the metallicity distribution in M33 is very complex, showing a central depression in metallicity probably due to observational bias; {\em ii}) the metallicity gradient in the disk of M33 has a slope of -0.037$\pm$ 0.009 dex kpc$^{-1}$ in the whole radial range up to $\sim$8 kpc, and -0.044$\pm$ 0.009 dex kpc$^{-1}$ excluding the central kpc; {\em iii}) there is a small evolution of the slope with time from the epoch of PN progenitor formation to the present-time.}

We test cosmological models of structure formation using the rotation curve of the nearest spiral galaxy, M31, determined using a recent deep, full-disk 21-cm imaging survey smoothed to 466 pc resolution. We fit a tilted ring model to the HI data from 8 to 37 kpc. The disk of M31 warps from 25 kpc outwards and becomes more inclined with respect to our line of sight. Newtonian dynamics without a dark matter halo provide a very poor fit to the rotation curve derived using the warp model. In the framework of modified Newtonian dynamic however the 21-cm rotation curve is well fitted by the gravitational potential traced by the baryonic matter density alone. The inclusion of a dark matter halo with a density profile as predicted by structure formation in a hierarchical clustering LambdaCDM cosmology makes the mass model in newtonian dynamic compatible with the rotation curve data. The dark halo concentration for the best fit is C=12 and its total mass is 1.2 10^{12} Msun. If a dark halo model with a constant density core is considered, the core radius has to be larger than 20 kpc in order for the model to fit to the data. We extrapolate the best-fit LambdaCDM and constant-density core mass models to very large galactocentric radii, comparable to the size of the dark matter halo. A comparison of the predicted mass with the M31 mass determined at such large radii using other dynamical tracers, confirms the validity of our results. In particular the LambdaCDM dark halo model which best fits the 21-cm data well reproduces the M31 mass traced out to 560 kpc. Our estimated total mass of M31 is 1.3 10^{12} Msun, with 12% baryonic fraction and only 6% of the baryons in neutral gas.

In the Local Group spiral galaxy M33, we investigate the correlation between the star formation rate (SFR) surface density, Sigma_SFR, and the gas density Sigma_gas (molecular, atomic, and total). We also explore whether there are other physical quantities, such as the hydrostatic pressure and dust optical depth, which establish a good correlation with Sigma_SFR. We use the Ha, far-ultraviolet (FUV), and bolometric emission maps to infer the SFR locally at different spatial scales, and in radial bins using azimuthally averaged values. Most of the local analysis is done using the highest spatial resolution allowed by gas surveys, 180 pc. The Kennicutt-Schmidt (KS) law, Sigma_SFR \propto (Sigma_gas)^n is analyzed by three statistical methods. At all spatial scales, with Ha emission as a SFR tracer, the KS indices n are always steeper than those derived with the FUV and bolometric emissions. We attribute this to the lack of Ha emission in low luminosity regions where most stars form in small clusters with an incomplete initial mass function at their high mass end. For azimuthally averaged values the depletion timescale for the molecular gas is constant, and the KS index is n_H2 = 1.1 +- 0.1. Locally, at a spatial resolution of 180 pc, the correlation between Sigma_SFR and Sigma_gas is generally poor, even though it is tighter with the molecular and total gas than with the atomic gas alone. Considering only positions where the CO J=1-0 line is above the 2-sigma detection threshold and taking into account uncertainties in Sigma_H2 and Sigma_SFR, we obtain a steeper KS index than obtained with radial averages: n_H2 = 2.22 +- 0.07 (for FUV and bolometric SFR tracers), flatter than that relative to the total gas (n_Htot = 2.59 +- 0.05). [abridged]

We have undertaken a deep, wide-field HI imaging survey of M31, reaching a maximum resolution of about 50 pc and 2 km/s across a 95x48 kpc region. The HI mass and brightness sensitivity at 100 pc resolution for a 25 km/s wide spectral feature is 1500 M_Sun and 0.28 K. Our study reveals ubiquitous HI self-opacity features, discernible in the first instance as filamentary local minima in images of the peak HI brightness temperature. Local minima are organized into complexes of more than kpc length and are particularly associated with the leading edge of spiral arm features. Just as in the Galaxy, there is only patchy correspondence of self-opaque features with CO(1-0) emission. Localized opacity corrections to the column density exceed an order of magnitude in many cases and add globally to a 30% increase in the atomic gas mass over that inferred from the integrated brightness under the usual assumption of negligible self-opacity. Opaque atomic gas first increases from 20 to 60 K in spin temperature with radius to 12 kpc but then declines again to 20 K beyond 25 kpc. We have extended the resolved star formation law down to physical scales more than an order of magnitude smaller in area and mass than has been possible previously. The relation between total-gas-mass- and star-formation-rate-density is significantly tighter than that with molecular-mass and is fully consistent in both slope and normalization with the power law index of 1.56 found in the molecule-dominated disk of M51 at 500 pc resolution. Below a gas-mass-density of about 5 M_Sun/pc^2, there is a down-turn in star-formation-rate-density which may represent a real local threshold for massive star formation at a cloud mass of about 5x10^4 M_Sun.

We test the reliability of infrared (IR) emission to trace star formation in individual star-forming sites of M33, and outline a new method for testing the distribution function of massive stars in newly formed clusters. We select IR sources from the Spitzer survey of M33 and show that the IR and Halpha luminosities are not correlated. Complementing the infrared photometry with GALEX-UV data, we estimate the source bolometric luminosities. For a given stellar IMF we simulate a theoretical curve for the expected bolometric-to-Halpha luminosity ratio, along which stellar clusters are born. We call this the cluster birthline in the Lbol--Lbol/LHal plane. The birthline is flat for Lbol>3x10^{39}erg/s because all clusters fully sample the IMF and it increases toward lower luminosities as the upper end of the IMF becomes incompletely sampled. The observations of M33 show that young isolated clusters lie close to the theoretical birthline for a wide range of Lbol. The luminosity is not proportional to Halpha emission for low mass clusters and aging moves clusters above the birthline. The best fit to the birthline is for a randomly sampled IMF, in which the mass of most massive star in a cluster is not strictly limited by the cluster's mass. We also find that the IR luminosity of young stellar clusters in M33 is not proportional to their bolometric luminosity. This irregularity could be the result of low and patchy dust abundance in M33.

Aims. We use different tracers, such as Ha, ultraviolet (UV), and infrared (IR) emissions at various wavelengths, to study the dust and star-formation (SF) conditions throughout the disk of M33. Methods. We derive the radial distribution of dust, of the old and young stellar population using Spitzer and GALEX data, complemented by ground-based optical data and available surveys of atomic and molecular gas. We separate the contribution of discrete sources to the IR brightness from the diffuse emission. Results. At 8 and 24 um, discrete sources account for 40% of the IR emission in the innermost 3 kpc, and for 20% further out. We find that stochastic emission from very small grains in the diffuse interstellar medium accounts for only 10% of the diffuse 24 um emission, and that dusty circumstellar shells of unresolved, evolved AGB stars (carbon stars) are a viable alternative. The 8 um profile suggests that PAH emission declines faster with radius than the dust continuum. In annular regions, 0.24 kpc wide, we find a mean extinction value for stellar continuum AV = 0.25 mag with a weak dependence on radius, consistent with the shallow metallicity gradient observed. Dust opacity derived from the 160 um emission decreases instead by a factor 10 from the center to edge of the SF disk. Conclusions. Using extinction corrected UV and Ha maps we find the global SF rate in M33, over the last 100 Myr, to be 0.45 +- 0.10 M yr-1. FIR (far-IR) and TIR (total-IR) luminosities can trace SF even though a high conversion factor is required to recover the effective rate. If carbon stars are powering the diffuse 24 um emission in M33 this can trace star formation 1 Gyr ago and provide a more complete view of the SF history of the galaxy. [abridged]

Neutral hydrogen clouds are found in the Milky Way and Andromeda halo both as large complexes and smaller isolated clouds. Here we present a search for Hi clouds in the halo of M33, the third spiral galaxy of the Local Group. We have used two complementary data sets: a 3^o x 3^o map of the area provided by the Arecibo Legacy Fast ALFA (ALFALFA) survey and deeper pointed observations carried out with the Arecibo telescope in two fields that permit sampling of the north eastern and south-western edges of the HI disc. The total amount of Hi around M33 detected by our survey is $\sim 10^7$ M$_{\odot}$. At least 50% of this mass is made of HI clouds that are related both in space and velocity to the galaxy. We discuss several scenarios for the origin of these clouds focusing on the two most interesting ones: $(a)$ dark-matter dominated gaseous satellites, $(b)$ debris from filaments flowing into M33 from the intergalactic medium or generated by a previous interaction with M31. Both scenarios seem to fit with the observed cloud properties. Some structures are found at anomalous velocities, particularly an extended HI complex previously detected by Thilker et al. (2002). Even though the ALFALFA observations seem to indicate that this cloud is possibly connected to M33 by a faint gas bridge, we cannot firmly establish its extragalactic nature or its relation to M33. Taking into account that the clouds associated with M33 are likely to be highly ionised by the extragalactic UV radiation, we predict that the total gas mass associated with them is > 5 x 10^7 M$_{\odot}$. If the gas is steadily falling towards the M33 disc it can provide the fuel needed to sustain a current star formation rate of 0.5 M$_{\odot}$ yr$^{-1}$.

Spiral galaxies appear to be dynamical systems whose disks are still forming at the current epoch and which continue to accrete mass. The presence of extraplanar gas in spirals indicates that galactic halos can contain at least part of the material needed to fuel the star formation activity in their disks. Here we present the analysis of the ALFALFA survey data in the region of M33 aimed at searching high velocity clouds around this Local Group galaxy. We find a varied population of HI clouds with masses ranging between 4 x 10^4 and few times 10^6 M_sun. We also detect an extended HI complex at anomalous velocities, whose extragalactic nature cannot be firmly established. We estimate that the total amount of neutral hydrogen mass associated to these clouds is around 10^7 M_sun.

The Local Group member M33 is a pure disk galaxy bearing no prominent bulge or stellar halo. It constitutes a challenge for any hierarchical galaxy formation theory and an ideal laboratory for studying quiescent star formation. Using multiwavelength observations of the gas and stellar component in this nearby galaxy we are able to constrain the gas accretion and star formation history. In the centermost region we find kinematical evidence of a weak bar, which explains the central light excess and the enhanced metallicity. In the more extended disk the lack of strong gradients of metal and dust abundances supports the picture that the slow radial decline of the star formation rate is due to a change in the large scale disk perturbations: bright HII regions and giant molecular clouds being born only in the inner disk. The analysis of the infrared Spitzer maps has however revealed hundreds of low luminosity star forming sites in places with a variety of dust content. These are essential ingredients for understanding the overall gas to star formation process in M33 and in more distant late type galaxies.

Combining the relative vicinity of the Local Group spiral galaxy M33 with the Spitzer images, we investigate the properties of infrared (IR) emission sites and assess the reliability of the IR emission as a star formation tracer. The mid- and far-IR emission of M33 was obtained from IRAC and MIPS images from the Spitzer archive. We compared the photometric results for several samples of three known types of discrete sources (HII regions, supernovae remnants and planetary nebulae) with theoretical diagnostic diagrams, and derived the spectral energy distribution (from 3.6 to 24 micron) of each type of object. Moreover, we generated a catalogue of 24 micron sources and inferred their nature from the observed and theoretical colours of the known type sources. We estimated the star formation rate in M33 both globally and locally, from the IR emission and from the Halpha emission line. The colours of the typical IR emissions of HII regions, supernovae remnants and planetary nebulae are continuous among the different samples, with overlapping regions in the diagnostic diagrams. The comparison between the model results and the colours of HII regions indicates a dusty envelope at relatively high temperatures ~600 K, and moderate extinction Av < 10. The 24 micron sources IR colours follow the regions observationally defined by the three classes of known objects but the majority of them represent HII regions. The derived total IR luminosity function is in fact very similar to the HII luminosity function observed in the Milky Way and in other late type spirals. Even though our completeness limit is 5x10^37 ergs s-1, in low density regions we are able to detect sources five times fainter than this, corresponding to the faintest possible HII region. [abridged]

Combining the relative vicinity of the Local Group spiral galaxy M33 with the Spitzer images, we investigate the properties of infrared (IR) emission sites and assess the reliability of the IR emission as a star formation tracer. We compared the photometric results for several samples of three known types of discrete sources (HII regions, supernovae remnants and planetary nebulae) with theoretical diagnostic diagrams, and derived the spectral energy distribution (from 3.6 to 24 microns) of each type of object. Moreover, we generated a catalogue of 24 microns sources and inferred their nature from the observed and theoretical colours of the known type sources. We estimated the star formation rate in M33 both globally and locally, from the IR emission and from the Halpha emission line.

We present measurements of the stellar and gaseous velocities in the central 5' of the Local Group spiral M33. The data were obtained with the ARC 3.5m telescope. Blue and red spectra with resolutions from 2 to 4\AA covering the principal gaseous emission and stellar absorption lines were obtained along the major and minor axes and six other position angles. The observed radial velocities of the ionized gas along the photometric major axis of M33 remain flat at ~22 km s^{-1} all the way into the center, while the stellar velocities show a gradual rise from zero to 22 km s^{-1} over that same region. The central star cluster is at or very close to the dynamical center, with a velocity that is in accordance with M33's systemic velocity to within our uncertainties. Velocities on the minor axis are non-zero out to about 1' from the center in both the stars and gas. Together with the major axis velocities, they point at significant deviations from circular rotation. The most likely explanation for the bulk of the velocity patterns are streaming motions along a weak inner bar with a PA close to that of the minor axis, as suggested by previously published IR photometric images. The presence of bar imprints in M33 implies that all major Local Group galaxies are barred. The non-circular motions over the inner 200 pc make it difficult to constrain the shape of M33's inner dark matter halo profile. If the non-circular motions we find in this nearby Sc galaxy are present in other more distant late-type galaxies, they might be difficult to recognize.

The evolution of radial gradients of metallicity in disk galaxies and its relation with the disk formation are not well understood. Theoretical models of galactic chemical evolution make contrasting predictions about the time evolution of metallicity gradients. To test chemical evolution models and trace the star formation and accretion history of low luminosity disk galaxies we focus on the Local Group galaxy M33. We analyze O/H and S/H abundances in planetary nebulae, H{\sc ii} regions, and young stars, together with known [Fe/H] abundances in the old stellar population of M33. With a theoretical model, we follow the time evolution of gas (diffuse and condensed in clouds), stars, and chemical abundances in the disk of M33, assuming that the galaxy is accreting gas from an external reservoir. Our model is able to reproduce the available observational constraints on the distribution of gas and stars in M33 and to predict the time evolution of several chemical abundances. In particular, we find that a model characterized by a continuous infall of gas on the disk, at a rate of $\dot M_{\rm inf}\approx 1$ $M_\odot$ yr$^{-1}$, almost constant with time, can also account for the relatively high rate of star formation and for the shallow chemical gradients. Supported by a large sample of high resolution observations for this nearby galaxy, we conclude that the metallicity in the disk of M33 has increased with time at all radii, with a continuous flattening of the gradient over the last $\sim 8$ Gyr.

The rotation curves and the relative mass distributions of the two nearby Local Group spiral galaxies, M31 and M33, show discrepancies with Modified Newtonian dynamic (MOND) predictions. In M33 the discrepancy lies in the kinematics of the outermost regions. It can be alleviated by adopting tilted ring models compatible with the 21-cm datacube but different from the one that best fits the data. In M31 MOND fails to fit the falling part of the rotation curve at intermediate radii, before the curve flattens out in the outermost regions. Newtonian dynamics in a framework of a stellar disc embedded in a dark halo can explain the complex rotation curve profiles of these two galaxies, while MOND has some difficulties. However, given the present uncertainties in the kinematics of these nearby galaxies, we cannot address the success or failure of MOND theory in a definite way. More sensitive and extended observations around the critical regions, suggested by MOND fits discussed in this paper, may lead to a definite conclusion.

Murphy et al. (2003a, MNRAS, 345, 609) claim to find evidence of cosmological variations of the fine structure constant $\alpha$ in the spectra of intervening QSO absorption line systems. We find that this result is affected by systematic effects. The $\alpha$ values estimated in individual line systems depend on the set of atomic transitions used and therefore the quoted dependence on the cosmic age may reflect the fact that different sets of transitions are used at different redshifts. A correlation between line shifts and relative optical depths of the atomic transitions is also present. This correlation is very tight for a high-redshift subsample and accounts for the anomalous dispersion of the $\alpha$ values found by Murphy et al. (2003a) in this subsample. The above correlations are consistent with a scenario in which gravitational redshift, caused by the gravity of the dark halo, contributes to the shift in frequency of individual components. Gravitational redshift causes differential spectral shifts of the same order as magnitude of those measured. In the presence of line misidentification, these shifts can be interpreted in terms of a variable $\alpha$. In order to verify the gravitational redshift hypothesis, a direct access to Murphy et al. (2003a) data, or to a large amount of new high resolution data, is necessary.

The SKA is a unique instrument to open a window on many aspects of galaxy formation and evolution which can be examined in our Local Universe. Here I will focus on the outermost regions of galaxies which can be observed with sufficient sensitivity and resolution at 21-cm as to enlighten the interplay between galaxies and the intergalactic gas, the competing race between the local dark matter gravity and the external ionizing radiation field. Tracing the gas distribution out to large galactocentric radii will be complementary to QSO's Ly$\alpha$ absorption studies for understanding the evolution of the dark and visible matter of the Universe.

Green Bank Telescope (GBT) 21cm observations have revealed a faint, yet extensive HI cloud population surrounding the Andromeda Galaxy (M31). The newfound objects are likely analogs to the high-velocity HI clouds (HVCs) seen around the Milky Way. At least 20 discrete features are detected within 50 kpc of the M31 disk, with radial velocities that are comparable to those of outer disk rotation. In addition, a filamentary ``halo'' component of at least 30 kpc extent is concentrated at the M31 systemic velocity. Some of the discrete features are organized into elongated systems with velocity continuity, suggestive of tidal streams. The discrete population can be characterized by a steep power-law distribution of number versus HI mass in the range between 10^5 and 10^7 M_sun. The velocity line-width of discrete clouds is correlated with the cloud HI mass: such that if the clouds are gravitationally bound this implies a dark- to HI mass ratio of ~ 100:1. Possible origins for the discrete and ``halo'' M31 features include: a Local Group ``cooling flow'', tidal debris from recent mergers or interactions, and the gaseous counterparts of low mass dark-matter halos.

The relationship between the star formation rate and surface density of neutral gas within the disk of M33 is examined with new imaging observations of CO J=1-0 emission gathered with the FCRAO 14m telescope and IRAS HiRes images of the 60 micron and 100 micron emission. The Schmidt law, Sigma_SFR ~ Sigma_gas^n, is constructed using radial profiles of the HI 21cm, CO, and far infrared emission. A strong correlation is identified between the star formation rate and molecular gas surface density. This suggests that the condensation of giant molecular clouds is the limiting step to star formation within the M33 disk. The corresponding molecular Schmidt index, n_{mol}, is 1.36 +/- 0.08. The star formation rate has a steep dependence on total mass gas surface density, (Sigma_{HI}+Sigma_{H_2}), owing to the shallow radial profile of the atomic gas which dominates the total gas surface density for most radii. The disk pressure of the gas is shown to play a prominent role in regulating the molecular gas fraction in M33.

In this paper we present new measurements of the gas kinematics in M33 using the CO J=1-0 line. The resulting rotational velocities complement previous 21-cm line data for a very accurate and extended rotation curve. The implied dark matter mass, within the total gaseous extent, is a factor 5 higher than the visible baryonic mass. Dark matter density profiles with an inner cusp as steep as R^{-1}, are compatible with M33 data. The dark matter concentrations required for fitting the M33 rotation curve are very low but marginally consistent with halos forming in a standard Cold Dark Matter cosmology. The M33 virialized dark halo is at least 50 times more massive than the visible baryons and its size is comparable with the M33-M31 separation. Inner cusps as steep as R^{-1.5} are ruled out, while halo models with a large size core of constant density are consistent with the M33 data. A central spheroid of stars is needed and we evaluate its dynamical mass range. Using accurate rotational velocity gradients and the azimuthally averaged baryonic surface densities, we show that a disk instability can regulate the star formation activity in M33. Considering the gaseous surface density alone, the predicted outer star formation threshold radius is consistent with the observed drop of the H-alpha surface brightness if a shear rate criterion is used. The classical Toomre criterion predicts correctly the size of the unstable region only when the stellar or dark halo gravity, derived in this paper, is added to that of the gaseous disk.

We present a catalogue of the 322 damped Lyman alpha absorbers taken from the literature. All damped Lyman alpha absorbers are included, with no selection on redshift or quasar magnitude. Of these, 123 are candidates and await confirmation using high resolution spectroscopy. For all 322 objects we catalogue the radio properties of the background quasars, where known. Around 60 quasars have radio flux densities above 0.1 Jy and approximately half of these have optical magnitudes brighter than V = 18. This compilation should prove useful in several areas of extragalactic/cosmological research.

Using data on the HI column density distribution in the local Universe, f(N_{HI}), we show how to determine g(N_{H}), the distribution of the total gas (HI+HII) column density. A simple power law fit to f(N_{HI}) fails due to bendings in the distributions when N_{HI}<10^{20} cm^{-2} and H is no longer fully neutral. If an ultraviolet background is responsible for the gas ionization, and g(N_{H}) is proportional N_{H}^{-alpha}, we find the values of alpha and of the intensity of the background radiation which are compatible with the present data. These best fitting values, however, depend upon the scaling law of the the gas volume densities with N_{H} and cannot be determined unambiguously. We examine in detail two models: one in which the average gas volume density decreases steadily with N_H, while in the other it stays constant at low column densities. The former model leads to a steep power law fit for g(N_{H}), with alpha of order 3.3 +- 0.4 and requires an ultraviolet flux larger than what the QSOs alone produce at z=0. For the latter alpha is 1.5 +- 0.1 and a lower ionizing flux is required. The ambiguities about the modelling and the resulting steep or shallow N_{H} distribution can be resolved only if new 21-cm observations and QSOs Lyman limit absorbers searches will provide more data in the HI-HII transition region at low redshifts. Using the best fit obtained for higher redshift data we outline two possible scenarios for the evolution of gaseous structures, compatible with the available data at z=0.

In this paper we show how to use data on Lyman-limit and Damped Lyman-alpha absorption systems to derive the hydrogen ionization fractions and the distribution of the face-on total gas column density. We consider axially symmetric, randomly oriented absorbers, ionized by an external background radiation field in order to relate the face-on total gas distribution to that of the neutral hydrogen observed along the line of sight. We devise a statistical procedure based on the Maximum Likelihood criterion, that is able to treat simultaneously data coming from different surveys and statistically recovers the "true" column densities in the presence of large uncertainties: this is especially important for Lyman-limit systems which leave an unmeasurable residual flux at wavelengths shorter than the Lyman break. We make use of simulated data to look for possible observational biases and extensively test our procedure. For a large statistical sample of real data in the redshift range [1.75,3.25] (collected from all published surveys) our Maximum Likelihood procedure gives a power-law slope for the total hydrogen distribution of -2.7. All together Lyman-limit systems therefore contain more gas than Damped Lyman-alpha systems. Analysis of data at other redshifts shows that more observations are needed to reach a compelling evidence for a cosmological evolution of the slope of the gas distribution.

We derive the distribution of neutral and ionized gas in high redshift clouds which are optically thick to hydrogen ionizing radiation, using published data on Lyman-limit and Damped Lyman-alpha absorption systems in the redshift range 1.75 < z < 3.25. We assume that the distribution of the hydrogen total (HI+HII) column density in the absorbers follows a power law K N_H^{-alpha}, whereas the observed HI column density distribution deviates from a pure power law as a result of ionization from a background radiation field. Comparison of the models and observations give Maximum Likelihood solutions for the exponent alpha and for X, the value of log(N_H/N_HI) when the Lyman-limit optical depth is unity: alpha=2.7^{+1.0}_{-0.7} and X=2.75\pm0.35. X is much lower than what would be obtained for a gaseous distribution in equilibrium under its own gravity but the ratio of dark matter to gas density is not well constrained being log(eta_0)=1.1\pm 0.8. An extrapolation of our derived power law distribution towards systems of lower column density, the Lyman-alpha forest, favours models with log(eta_0) < 1.1 and alpha=2.7-3.3. With alpha appreciably larger than 2, Lyman-limit systems contain more gas than Damped Lyman-alpha systems and Lyman-alpha forest clouds even more. Estimates of the cosmological gas and dark matter density due to absorbers of different column density around z=2.5 are also given.

We present the 21-cm rotation curve of the nearby galaxy M33 out to a galactocentric distance of 16 kpc (13 disk scale-lengths). The rotation curve keeps rising out to the last measured point and implies a dark halo mass larger than 5 10^{10} solar masses. The stellar and gaseous disks provide virtually equal contributions to the galaxy gravitational potential at large galactocentric radii but no obvious correlation is found between the radial distribution of dark matter and the distribution of stars or gas. Results of the best fit to the mass distribution in M33 picture a dark halo which controls the gravitational potential from 3 kpc outward, with a matter density which decreases radially as R^{-1.3}. The density profile is consistent with the theoretical predictions for structure formation in hierarchical clustering cold dark matter models but mass concentrations are lower than those expected in the standard cosmogony.

We show that gaseous \HI disks of primordial composition irradiated by an external radiation field can develop a multiphase medium with temperatures between 10^2 and 10^4 K due to the formation of molecular hydrogen. For a given \HI column density there is a critical value of the radiation field below which only the cold \HI phase can exist. Due to a time decreasing quasar background, the gas starts cooling slowly after recombination until the lowest stable temperature in the warm phase is reached at a critical redshift $z=z_{cr}$. Below this redshift the formation of molecular hydrogen promotes a rapid transition towards the cold \HI phase. We find that disks of protogalaxies with $10^{20}\simlt N_{HI}\simlt 10^{21} cm^{-2}$ are gravitationally stable at $T\sim 10^4$ K and can start their star formation history only at $z \simlt z_{cr}\sim 2$, after the gas in the central portion of the disk has cooled to temperatures $T\simlt 300$ K. Such a delayed starbust phase in galaxies of low gas surface density and low dynamical mass can disrupt the disks and cause them to fade away. These objects could contribute significantly to the faint blue galaxy population.

To determine the actual HI distribution and the velocity field in the outermost disk of the spiral galaxy M33, a tilted-ring model is fitted to 21-cm line data taken with the Arecibo Telescope. Since M33 is one of the main calibrators for the extragalactic distance scale derived through the Tully-Fisher relation, the outer disk warping is of interest for a correct determination and deprojection of the galaxy's line width. Even though our best model predicts small effects on the observed line width of M33, we show that similar outer disk warping in galaxies oriented differently along our line of sight could affect the widths considerably. Therefore there may be systematic effects in the determination of the rotation velocities and dynamic masses of spiral galaxies, whose exact value depends also on which method is used for measuring the galaxy's total line width.

We discuss a possible delay experienced by protogalaxies with low column density of gas in forming stars over large scales. After the hydrogen has recombined, as the external ionizing UV flux decreases and the metal abundance $Z$ increases, the HI, initially in the warm phase ($T\simgt 5000$ K), makes a transition to the cool phase ($T\simlt 100$ K). The minimum abundance $Z_{min}$ for which this phase transition takes place in a small fraction of the Hubble time decreases rapidly with increasing gas column density. Therefore in the ``anemic'' disk galaxies, where $N_{HI}$ is up to ten times smaller than for normal large spirals, the onset of the cool HI phase is delayed. The onset of gravitational instability is also delayed, since these objects are more likely to be gravitationally stable in the warm phase than progenitors of today's large spiral galaxies. The first substantial burst of star formation may occur only as late as at redshifts $z \sim 0.5$ and give a temporary high peak luminosity, which may be related to the ``faint blue objects". Galaxy disks of lower column density tend to have lower escape velocities and a starburst/galactic fountain instability which decreases the gas content of the inner disk drastically.

We present a linear analysis of acoustic and thermo-reactive instabilities in a diffuse gas, photoionized and heated by a radiation field, cooled by collisional excitation of hydrogen and metal lines. The hydrogen recombination reaction has a stabilizing effect on the thermal mode found by Field (1965) since the condensation instability occurs in a narrower region of the parameter space and grows on longer time scales due to its oscillatory character. This effect is stronger when the mean photon energy is not much larger than the hydrogen ionization energy. Moreover, for fixed values of thermal pressure and photoionization rate, there are thermo- reactive unstable equilibria for which no transition to a stable phase is possible. By extending our analysis of the thermo-reactive modes to the nonlinear regime we show that when no phase transition is possible the medium evolves through a series of nonequilibrium states characterized by large amplitude, nonlinear periodic oscillations of temperature, density and hydrogen ionization fraction. We find also unstable acoustic waves which, for solar metal abundances, are the fastest growing modes in two temperature intervals: around $T\sim 100$~K and $T\sim 8000$~K (i.e. cold and warm phase, respectively), independent of the mean photon energy. Possible implications for the interstellar medium and intergalactic medium are briefly outlined.

We study thermoreactive and acoustic instabilities in a diffuse gas, photoionized and heated by a radiation field. The analysis of the thermal instability by Field (1965) is extended to include the effects of the hydrogen recombination reaction, which, in general, is found to be a stabilizing agent for the condensation mode. This effect is stronger when the mean photon energy is not much larger than the hydrogen ionization energy. In addition, there are thermoreactive unstable equilibria for which an isobaric transition to a stable phase is not possible and the system evolves toward a dynamic state characterized by large amplitude, nonlinear periodic oscillations of temperature, density and hydrogen ionization fraction. Acoustic waves are found to be unstable for some temperatures of both the cold and the warm phase of the ISM, also in regions where the gas is {\it thermally stable}, independently of the mean photon energy.

ABSRACT:We show how to use 21-cm emission and absorption studies to estimate the heat inputs to the neutral gas in low pressure environments, such as in outer disks of spiral galaxies, in galactic halos or in intergalactic space. For a range of model parameters we calculate the gas kinetic and spin temperatures ($T_K$ and $T_S$) and the relation between $T_S$ and the heat input to the gas. We outline the conditions for a ``two phase medium'' to exist. We find that although $T_S$ can be much smaller than $T_K$, $T_S$ is always $ \gg 3$ K for column densities greater that $5 \times 10^{18}$ cm$^{-2}$. This excludes the possibility that relevant HI masses at the periphery of galaxies are invisible at 21-cm in emission. The outermost interstellar gas in a disk galaxy is more directly affected by external processes and in this paper we estimate the intensity of the extragalactic background at energies close to 0.1 keV by comparing our theoretical results with HI emission/absorption studies. We take into account the possibility that some energy produced in the inner regions affects the energy balance in outer regions. We find that in the absence of any other local heat source QSO dominated background models are still compatible with the spin temperature limits derived for the two best documented HI emission/absorption studies in outer regions.

Observations indicate that some extended outer disks have a sharp cut off in the surface density of neutral hydrogen when this approaches the value of $\sim 2\times 10^{19}$ cm$^{-2}$. In this paper we model these HI edges as places where the ratio of neutral to ionized hydrogen drops rapidly due to ionizing radiation. We use two different models for the vertical distribution of gas above the outer galactic plane: in the first model we derive the density from the ideal gas law while in the second model we insert a macroscopic pressure term and derive the density as for an isothermal slab. We consider two different sources of ionizing photons: external fluxes of different intensity and spectral index due to quasars, and a monochromatic UV flux due to neutrino decays. We find that galaxies which have a smaller gas scale height should show outer HI disks to a lower column density and smoother HI edges. The sharpness of the HI-HII transition and the total column density at which the medium is 50$\%$ ionized, are strongly correlated. We present several model fits to the HI sharp edge observed in M33. If today's UV background is dominated by attenuated quasar light, which gives $\sim 10^{-14}$ H ionizations s$^{-1}$, a large gas scale height or equivalently a nearly spherical halo is preferred. If ionizing photons from decaying neutrinos are responsible for the sharp edge, a flat dark matter halo, is required.

3C275.1 and \NGC4651 make a particularly interesting quasar/galaxy pairing because of the alignment of such a strong radio emitter behind the outer H\I disk of a relatively undisturbed spiral galaxy. This provides an opportunity to study the spin-temperature characteristics of atomic hydrogen at low column densities, in an apparently star-free environment. We previously reported a tentative detection of absorption against the quasar based on VLA C--array observations; we have now made more sensitive maps of the H\I emission from \NGC4651 with the VLA D--array, and we have attempted to confirm the weak H\I absorption against the quasar at higher spatial and spectral resolution in VLA B--array. The possible absorption feature against this quasar appears to be weaker than we previously suspected, even though it seems fairly clear that H\I emission is present close to the line of sight to the quasar.The weakness of the possible absorption seems also to confirm, conversely, the trend found in previous observations that where strong absorption lines are seen, the galaxies show evidence of disturbance. The possible detection of (or limits on) absorption suggest that the neutral gas in the outer disk is quite warm. We use the absorption and emission measurements to set lower limits on the combination of heating inputs outside the star-forming regions of a disk galaxy and/or the intensity of the cosmic background radiation around 100 eV.