The Molecular Gas Distribution and Schmidt Law in M33

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.Comment: 19 pages + 5 figures. Accepted for publication in Ap

HI clouds in the proximity of M33

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}$.Comment: 16 pages, 19 figures. Accepted for publication in A&

Bar imprints on the inner gas kinematics of M33

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.Comment: 20 pages, 12 figures, ApJ in pres

Simulations of the flocculent spiral M33: what drives the spiral structure?

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.Comment: 17 pages, 17 figures, accepted for publication in MNRA

A warped disk model for M33 and the 21-cm line width in spiral galaxies

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.Comment: 27 pages, ps files only, ApJ in pres

A Comprehensive Statistical Analysis of the Gas Distribution in Lyman-limit and Damped Lyman-alpha Absorption Systems

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.Comment: 30 pages with 7 eps figures, LaTeX accepted for publication in ApJ main journa

Star formation in disk galaxies driven by primordial H_2

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.Comment: 16 pages (LaTeX), 2 Figures to be published in Astrophysical Journal Letter

Sharp HI edges at high z: the gas distribution from Damped Lyman-alpha to Lyman-limit absorption systems

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.Comment: 21 pages, 6 figures. Accepted for publication in Ap

Radial HI Profiles at the Periphery of Galactic Disks: The Role of Ionizing Background Radiation

Observations of neutral hydrogen in spiral galaxies reveal a sharp cutoff in the radial density profile at some distance from the center. Using 22 galaxies with known HI distributions as an example, we discuss the question of whether this effect can be associated exclusively with external ionizing radiation, as is commonly assumed. We show that before the surface density reaches $\sigma_{\textrm{HI}}\le 0.5 {\cal M}_\odot/{\textrm {pc}}^2$(the same for galaxies of different types), it is hard to expect the gas to be fully ionized by background radiation. For two of 13 galaxies with a sharp drop in the HI profile, the "steepening" can actually be caused by ionization. At the same time, for the remaining galaxies, the observed cutoff in the radial HI profile is closer to the center than if it was a consequence of ionization by background radiation and, therefore, it should be caused by other factors.Comment: 15 pages, 6 figure