Source author record

F. P. Israel

F. P. Israel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2020arXiv

Central molecular zones in galaxies: 12CO-to-13CO ratios, carbon budget, and X factors

This paper presents ground-based 12CO and 13CO measurements of 126 nearby galaxy centers in various J transitions. More than 60 galaxies were measured in at least four lines. The average relative intensities of the first four 12CO transitions are 1.00 : 0.92 : 0.70 : 0.57. The average 12CO-to-13CO ratios are 13.0, 11.6, and 12.8 for the first three transitiions. The sizes of central CO concentrations are well defined in maps, but poorly determined by multi-aperture photometry. Using radiative transfer models (RADEX), we derived model gas parameters, where the assumed carbon elemental abundances and carbon gas depletion onto dust are the main causes of uncertainty. The new CO data and published [CI] and [CII] data imply that CO, Co, and C+ each represent about one-third of the gas-phase carbon in the molecular interstellar medium in galaxy centers Their mean beam-averaged central molecular hydrogen column density is N(H2) = (1.5+/-0.2) x 10^(21) cm^(2) and CO-to-H2 conversion factors are typically ten times lower than the `standard' Milky Way disk value, with a mean X(CO) = (1.9+/-0.2) x 10^(19) cm^(2)/K km/s. The corresponding [CI]-to-H2 factor is five times higher than X(CO), with X[CI] = (9+/-2) x 10^(19) cm^(2)/K km/s. No unique conversion factor can be determined for [CII]. The low molecular gas content of galaxy centers relative to their CO intensities is explained in roughly equal parts by higher central gas-phase carbon abundances, elevated gas temperatures, and larger gas velocity dispersions than found in galaxy disks.

preprint2019arXiv

Star Formation Efficiencies at Giant Molecular Cloud Scales in the Molecular Disk of the Elliptical Galaxy NGC 5128 (Centaurus A)

We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previously seen, up to a radius of 4 kpc. We find that the global star formation rate is $\sim$1 \Msol yr$^{-1}$, which yields a star formation efficiency (SFE) of 0.6 Gyr$^{-1}$ (depletion time $τ=$1.5 Gyr), similar to those in disk galaxies. We show the most detailed view to date (40\,pc resolution) of the relation between molecular gas and star formation within the stellar component of an elliptical galaxy, from several kpc scale to the circumnuclear region close to the powerful radio jet. Although on average the SFEs are similar to those of spiral galaxies, the circumnuclear disk (CND) presents SFEs of 0.3 Gyr$^{-1}$, lower by a factor of 4 than the outer disk. The low SFE in the CND is in contrast to the high SFEs found in the literature for the circumnuclear regions of some nearby disk galaxies with nuclear activity, probably as a result of larger shear motions and longer AGN feedback. The higher SFEs in the outer disk suggests that only central molecular gas or filaments with sufficient density and strong shear motions will remain in $\sim$1 Gyr, which will later result in the compact molecular distributions and low SFEs usually seen in other giant ellipticals with cold gas.

preprint2016arXiv

Carbon gas in SMC low-metallicity star-forming regions

This paper presents [CII], [CI] and CO emission line maps of the star-forming regions N66, N25+N26, and N88 in the metal-poor Local Group dwarf galaxy SMC. The spatial and velocity structure of the large HII region N66 reveals an expanding ring of shocked molecular gas centered on the exciting star cluster NGC346, whereas a more distant dense molecular cloud is being eroded by UV radiation from the same cluster. In the N25+N26 and N88 maps, diffuse [CII] emission at a relatively low surface brightness extends well beyond the compact boundaries of the bright emission associated with the [CII] regions. In all regions, the distribution of this bright [CII] emission and the less prominent [CI] emission closely follows the outline of the CO complexes, but the intensity of the [CII] and [CI] emission is generally anticorrelated, which can be understood by the action of photodissociation and photoionization processes. Notwithstanding the overall similarity of CO and [CII] maps, the intensity ratio of these lines varies significantly, mostly due to changes in CO brightness. [CII] emission line profiles are up to 50% wider in velocity than corresponding CO profiles. A radiative transfer analysis shows that the [CII] line is the dominant tracer of (CO-dark) molecular hydrogen in the SMC. CO emission traces only a minor fraction of the total amount of gas. The similarity of the spatial distribution and line profile shape, and the dominance of molecular gas associated with [[CII] rather than CO emission imply that in the low-metallicity environment of the SMC the small amount of dense molecular gas traced by CO is embedded in the much more extended molecular gas traced only by [CII] emission. The contribution from neutral atomic and ionized hydrogen zones is negligible in the star-forming regions observed.

preprint2016arXiv

Constraining cloud parameters using high density gas tracers in galaxies

Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic scales. We do this by re-sampling high density gas in a hydrodynamic model of a gas-rich disk galaxy, assuming that the density field of the interstellar medium of the model galaxy follows the probability density function (PDF) inferred from the resolved low density scales. We find that in a narrow gas density PDF, with a mean density of ~10 cm-3 and a dispersion σ= 2.1 in the log of the density, most of the emission of molecular lines, emanates from the 10-1000 cm-3 part of the PDF. We construct synthetic emission maps for the central 2 kpc of the galaxy and fit the line ratios of CO and 13CO up to J = 15-14, as well as HCN, HNC, and HCO+ up to J = 7-6, using one photo-dissociation region (PDR) model. We attribute the goodness of the one component fits for our model galaxy to the fact that the distribution of the luminosity, as a function of density, is peaked at gas densities between 10 and 1000 cm-3. We explore the impact of different log-normal density PDFs on the distribution of the line-luminosity as a function of density, and we show that it is necessary to have a broad dispersion, corresponding to Mach numbers >~ 30 in order to obtain significant emission from n > 10^4 cm-3 gas. Such Mach numbers are expected in star-forming galaxies, LIRGS, and ULIRGS. By fitting line ratios of HCN(1-0), HNC(1-0), and HCO+(1-0) for a sample of LIRGS and ULIRGS using mechanically heated PDRs, we constrain the Mach number of these galaxies to 29 < M < 77.

preprint2016arXiv

Continuum Observations of M51 and M83 at 1.1 mm with AzTEC

We observed the spiral galaxies M51 and M83 at 20" spatial resolution with the bolometer array AzTEC on the JCMT in the 1.1$\,$mm continuum, recovering the extended emission out to galactocentric radii of more than 12 kpc in both galaxies. The 1.1 mm-continuum fluxes are 5.6+/-0.7 and 9.9+/-1.4 Jy, with associated gas masses estimated at 9.4 X 10^9 Mo and 7.2 X 10^9 Mo for M51 and M83, respectively. In the interarm regions of both galaxies the N(H2)/I(CO) (or X-factor) ratios exceed those in the arms by factors of ~1.5-2. In the inner disks of both galaxies, the X-factor is about 1 X 10^20 cm^-2 / (K km s^-1). In the outer parts, the CO-dark molecular gas becomes more important. While the spiral density wave in M51 appears to influence the interstellar medium and stars in a similar way, the bar potential in M83 influences the interstellar medium and the stars differently. We confirm the result of Foyle et al. (2010) that the arms merely heighten the star formation rate and the gas surface density in the same proportion. Our maps reveal a threshold gas surface density for an SFR increase by two or more orders of magnitude. In both galaxy centers, the molecular gas depletion time is about 1 Gyr climbing to 10-20 Gyr at radii of 6-8 kpc. This is consistent with an inside-out depletion of the molecular gas in the disks of spiral galaxies.

preprint2016arXiv

Modelling Mechanical Heating in Star-Forming Galaxies: CO and 13CO Line Ratios as Sensitive Probes

We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3. We consider two galaxy models: a small disk galaxy of solar metallicity and a lighter dwarf galaxy with 0.2 \zsun metallicity. Elevated excitation temperatures for CO(1 - 0) correlate positively with mechanical feedback, that is enhanced towards the central region of both model galaxies. The emission maps of these model galaxies are used to compute line ratios of CO and 13CO transitions. These line ratios are used as diagnostics where we attempt to match them These line ratios are used as diagnostics where we attempt to match them to mechanically heated single component (i.e. uniform density, Far-UV flux, visual extinction and velocity gradient) equilibrium PDR models. We find that PDRs ignoring mechanical feedback in the heating budget over-estimate the gas density by a factor of 100 and the far-UV flux by factors of ~10 - 1000. In contrast, PDRs that take mechanical feedback into account are able to fit all the line ratios for the central < 2 kpc of the fiducial disk galaxy quite well. The mean mechanical heating rate per H atom that we recover from the line ratio fits of this region varies between $10^{-27}$ -- $10^{-26}$~erg s$^{-1}$.

preprint2016arXiv

Planck intermediate results. XLV. Radio spectra of northern extragalactic radio sources

Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The radio spectra peak at high frequencies and exhibit a variety of shapes. For a small set of low-z sources, we find a spectral upturn at high frequencies, indicating the presence of intrinsic cold dust. Variability can generally be approximated by achromatic variations, while sources with clear signatures of evolving shocks appear to be limited to the strongest outbursts.

preprint2015arXiv

Neutral carbon and CO in 76 (U)LIRGs and starburst galaxy centers A method to determine molecular gas properties in luminous galaxies

We present fluxes in both neutral carbon [CI] lines at the centers of 76 galaxies with FIR luminosities between 10^{9} and 10^{12} L(o) obtained with Herschel-SPIRE and with ground-based facilities, along with the J=7-6, J=4-3, J=2-1 12CO and J=2-1 13CO line fluxes. We investigate whether these lines can be used to characterize the molecular ISM of the parent galaxies in simple ways and how the molecular gas properties define the model results. In most starburst galaxies, the [CI]/13CO flux ratio is much higher than in Galactic star-forming regions, and it is correlated to the total FIR luminosity. The [CI](1-0)/CO(4-3), the [CI](2-1) (2-1)/CO(7-6), and the [CI] (2-1)/(1-0) flux ratios are also correlated, and trace the excitation of the molecular gas. In the most luminous infrared galaxies (LIRGs), the ISM is fully dominated by dense and moderately warm gas clouds that appear to have low [C]/[CO] and [13CO]/[12CO] abundances. In less luminous galaxies, emission from gas clouds at lower densities becomes progressively more important, and a multiple-phase analysis is required to determine consistent physical characteristics. Neither the CO nor the [CI] velocity-integrated line fluxes are good predictors of H2 column densities in individual galaxies, and X(CI) conversion factors are not superior to X(CO) factors. The methods and diagnostic diagrams outlined in this paper also provide a new and relatively straightforward means of deriving the physical characteristics of molecular gas in high-redshift galaxies up to z=5, which are otherwise hard to determine.

preprint2015arXiv

Planck intermediate results. XXV. The Andromeda Galaxy as seen by Planck

The Andromeda Galaxy (M31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiral arms and sub-features. We examine the morphology of this long-wavelength dust emission as seen by Planck, including a study of its outermost spiral arms, and investigate the dust heating mechanism across M31. We find that dust dominating the longer wavelength emission ($\gtrsim 0.3\,$mm) is heated by the diffuse stellar population (as traced by 3.6$\,μ$m emission), with the dust dominating the shorter wavelength emission heated by a mix of the old stellar population and star-forming regions (as traced by 24$\,μ$m emission). We also fit spectral energy distributions (SEDs) for individual 5' pixels and quantify the dust properties across the galaxy, taking into account these different heating mechanisms, finding that there is a linear decrease in temperature with galactocentric distance for dust heated by the old stellar population, as would be expected, with temperatures ranging from around 22$\,$K in the nucleus to 14$\,$K outside of the 10$\,$kpc ring. Finally, we measure the integrated spectrum of the whole galaxy, which we find to be well-fitted with a global dust temperature of ($18.2\pm1.0$)$\,$K with a spectral index of $1.62\pm0.11$ (assuming a single modified blackbody), and a significant amount of free-free emission at intermediate frequencies of 20-60$\,$GHz, which corresponds to a star formation rate of around $0.12$M$_\odot\,$yr$^{-1}$. We find a $2.3\,σ$ detection of the presence of spinning dust emission, with a 30$\,$GHz amplitude of $0.7\pm0.3\,$Jy, which is in line with expectations from our Galaxy.

preprint2015arXiv

The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $μ$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $μ$m, [OI] 63 $μ$m, and [OI] 145 $μ$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $μ$m, and [CI] 609 $μ$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $μ$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $μ$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

preprint2015arXiv

The JCMT Nearby Galaxies Legacy Survey X. Environmental Effects on the Molecular Gas and Star Formation Properties of Spiral Galaxies

We present a study of the molecular gas properties in a sample of 98 HI - flux selected spiral galaxies within $\sim25$ Mpc, using the CO $J=3-2$ line observed with the James Clerk Maxwell Telescope. We use the technique of survival analysis to incorporate galaxies with CO upper limits into our results. Comparing the group and Virgo samples, we find a larger mean H$_{2}$ mass in the Virgo galaxies, despite their lower mean HI mass. This leads to a significantly higher H$_{2}$ to HI ratio for Virgo galaxies. Combining our data with complementary H$α$ star formation rate measurements, Virgo galaxies have longer molecular gas depletion times compared to group galaxies, due to their higher H$_{2}$ masses and lower star formation rates. We suggest that the longer depletion times may be a result of heating processes in the cluster environment or differences in the turbulent pressure. From the full sample, we find that the molecular gas depletion time has a positive correlation with the stellar mass, indicative of differences in the star formation process between low and high mass galaxies, and a negative correlation between the molecular gas depletion time and the specific star formation rate.

preprint2014arXiv

Diagnostics of the molecular component of PDRs with mechanical heating. II: line intensities and ratios

CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an object and also study its significance on estimating the H2 mass. Equilibrium PDRs models were used to compute the thermal and chemical balance for the clouds. The equilibria were solved for numerically using the optimized version of the Leiden PDR-XDR code. Large velocity gradient calculations were done as post-processing on the output of the PDR models using RADEX. High-J CO line ratios are very sensitive to MH. Emission becomes at least one order of magnitude brighter in clouds with n~10^5~cm^-3 and a star formation rate of 1 Solar Mass per year (corresponding to a MH rate of 2 * 10^-19 erg cm^-3 s^-1). Emission of low-J CO lines is not as sensitive to MH, but they do become brighter in response to MH. Generally, for all of the lines we considered, MH increases excitation temperatures and decreases the optical depth at the line centre. Hence line ratios are also affected, strongly in some cases. Ratios involving HCN are a good diagnostic for MH, such as HCN(1-0)/CO(1-0) and HCN(1-0)/HCO^+(1-0). Both ratios increase by a factor 3 or more for a MH equivalent to > 5 percent of the surface heating, as opposed to pure PDRs. The first major conclusion is that low-J to high-J intensity ratios will yield a good estimate of the MH rate (as opposed to only low-J ratios). The second one is that the MH rate should be taken into account when determining A_V or equivalently N_H, and consequently the cloud mass. Ignoring MH will also lead to large errors in density and radiation field estimates.

preprint2014arXiv

Molecular gas heating in Arp 299

Understanding the heating and cooling mechanisms in nearby (Ultra) luminous infrared galaxies can give us insight into the driving mechanisms in their more distant counterparts. Molecular emission lines play a crucial role in cooling excited gas, and recently, with Herschel Space Observatory we have been able to observe the rich molecular spectrum. CO is the most abundant and one of the brightest molecules in the Herschel wavelength range. CO transitions are observed with Herschel, and together, these lines trace the excitation of CO. We study Arp 299, a colliding galaxy group, with one component harboring an AGN and two more undergoing intense star formation. For Arp 299 A, we present PACS spectrometer observations of high-J CO lines up to J=20-19 and JCMT observations of $^{13}$CO and HCN to discern between UV heating and alternative heating mechanisms. There is an immediately noticeable difference in the spectra of Arp 299 A and Arp 299 B+C, with source A having brighter high-J CO transitions. This is reflected in their respective spectral energy line distributions. We find that photon-dominated regions (PDRs) are unlikely to heat all the gas since a very extreme PDR is necessary to fit the high-J CO lines. In addition, this extreme PDR does not fit the HCN observations, and the dust spectral energy distribution shows that there is not enough hot dust to match the amount expected from such an extreme PDR. Therefore, we determine that the high-J CO and HCN transitions are heated by an additional mechanism, namely cosmic ray heating, mechanical heating, or X-ray heating. We find that mechanical heating, in combination with UV heating, is the only mechanism that fits all molecular transitions. We also constrain the molecular gas mass of Arp 299 A to 3e9 Msun and find that we need 4% of the total heating to be mechanical heating, with the rest UV heating.

preprint2014arXiv

Radiative and mechanical feedback into the molecular gas of NGC 253

Starburst galaxies are undergoing intense periods of star formation. Understanding the heating and cooling mechanisms in these galaxies can give us insight to the driving mechanisms that fuel the starburst. Molecular emission lines play a crucial role in the cooling of the excited gas. With SPIRE on the Herschel Space Observatory we have observed the rich molecular spectrum towards the central region of NGC 253. CO transitions from J=4-3 to 13-12 are observed and together with low-J line fluxes from ground based observations, these lines trace the excitation of CO. By studying the CO excitation ladder and comparing the intensities to models, we investigate whether the gas is excited by UV radiation, X-rays, cosmic rays, or turbulent heating. Comparing the $^{12}$CO and $^{13}$CO observations to large velocity gradient models and PDR models we find three main ISM phases. We estimate the density, temperature,and masses of these ISM phases. By adding $^{13}$CO, HCN, and HNC line intensities, we are able to constrain these degeneracies and determine the heating sources. The first ISM phase responsible for the low-J CO lines is excited by PDRs, but the second and third phases, responsible for the mid to high-J CO transitions, require an additional heating source. We find three possible combinations of models that can reproduce our observed molecular emission. Although we cannot determine which of these are preferable, we can conclude that mechanical heating is necessary to reproduce the observed molecular emission and cosmic ray heating is a negligible heating source. We then estimate the mass of each ISM phase; $6\times 10^7$ M$_\odot$ for phase 1 (low-J CO lines), $3\times 10^7$ M$_\odot$ for phase 2 (mid-J CO lines), and $9\times 10^6$ M$_\odot$ for phase 3 (high-J CO lines) for a total system mass of $1\times10^{8}$ M$_\odot$.

preprint2014arXiv

The molecular circumnuclear disk (CND) in Centaurus A: A multi-transition CO and [CI] survey with Herschel, APEX, JCMT, and SEST

We present new CO and C^o line measurements of the compact circumnuclear disk in the center of NGC 128 (Centaurus~A) obtained with the Herschel Space Observatory, as well as SEST, JCMT, and APEX. The Cen A center CO ladder is quite different from those of either star-burst galaxies or AGNs. In addition, the relative intensity of the central Cen A [CI] emission lines is much greate than that in any other galaxy. The CO surface brightness of the compact circumnuclear disk (CND) is significantly higher than that of the much more extended thin disk (ETD) in the same line of sight. Our LVG and PDR/XDR models suggest that much of the CND gas is relatively cool (25 - 80 K) and not very dense (~ 300 cm^{-3}) if the heating is by UV photons, although there is some gas in both the CND and the ETD with a much higher density of ~30 000 cm^{-3}. Finally, there is also high-excitation, high-density phase in the CND (but not in the ETD), either in the form of an extreme PDR but more likely in the form of an XDR. The total gas mass of the Cen A CND is 8.4 x 10^{7} M(sun), uncertain by a factor of two. The CO-H2 conversion factor is 4 x 10^{20} K km/s, also within a factor of two.

preprint2014arXiv

The nature of the low-frequency emission of M51: First observations of a nearby galaxy with LOFAR

The grand-design spiral galaxy M51 was observed with the LOFAR High Frequency Antennas (HBA) and imaged in total intensity and polarisation. This observation covered the frequencies between 115 MHz and 175 MHz. We produced an image of total emission of M51 at the mean frequency of 151 MHz with 20 arcsec resolution and 0.3 mJy rms noise, which is the most sensitive image of a galaxy at frequencies below 300 MHz so far. The integrated spectrum of total radio emission is described well by a power law, while flat spectral indices in the central region indicate thermal absorption. We observe that the disk extends out to 16 kpc and see a break in the radial profile near the optical radius of the disk. Our main results, the scale lengths of the inner and outer disks at 151 MHz and 1.4 GHz, arm--interarm contrast, and the break scales of the radio--far-infrared correlations, can be explained consistently by CRE diffusion, leading to a longer propagation length of CRE of lower energy. The distribution of CRE sources drops sharply at about 10 kpc radius, where the star formation rate also decreases sharply. We find evidence that thermal absorption is primarily caused by HII regions. The non-detection of polarisation from M51 at 151 MHz is consistent with the estimates of Faraday depolarisation. Future searches for polarised emission in this frequency range should concentrate on regions with low star formation rates.

preprint2013arXiv

A CO J=3-2 map of M51 with HARP-B: Radial properties of the spiral structure

We present the first complete CO J=3-2 map of the nearby grand-design spiral galaxy M51 (NGC 5194), at a spatial resolution of ~600 pc, obtained with the HARP-B instrument on the James Clerk Maxwell Telescope. The map covers the entire optical galaxy disk and out to the companion NGC 5195, with CO J=3-2 emission detected over an area of ~9'x6' (~21x14 kpc). We describe the CO J=3-2 integrated intensity map and combine our results with maps of CO J=2-1, CO J=1-0 and other data from the literature to investigate the variation of the molecular gas, atomic gas and polycyclic aromatic hydrocarbon (PAH) properties of M51 as a function of distance along the spiral structure on sub-kpc scales. We find that for the CO J=3-2 and CO J=2-1 transitions there is a clear difference between the variation of arm and inter-arm emission with galactocentric radius, with the inter-arm emission relatively constant with radius and the contrast between arm and inter-arm emission decreasing with radius. For CO J=1-0 and HI the variation with radius shows a similar trend for the arm and inter-arm regions, and the arm-inter-arm contrast appears relatively constant with radius. We investigate the variation of CO line ratios (J=3-2/2-1, J=2-1/1-0 and J=3-2/1-0) as a function of distance along the spiral structure. Line ratios are consistent with the range of typical values for other nearby galaxies in the literature. The highest CO J=3-2/2-1 line ratios are found in the central ~1 kpc and in the spiral arms and the lowest line ratios in the inter-arm regions.We find no clear evidence of a trend with radius for the spiral arms but for the inter-arm regions there appears to be a trend for all CO line ratios to increase with radius. We find a strong relationship between the ratio of CO J=3-2 intensity to stellar continuum-subtracted 8mu PAH surface brightness and the CO J=3-2 intensity that appears to vary with radius.

preprint2013arXiv

Dense Gas in M33 (HerM33es)

We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial distance of 3.4kpc. Optical, far-infrared, and submillimeter data from Herschel and other observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of photon-dominated regions, one for solar and one for M33-type subsolar metallicity. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The 12CO/13CO line ratio varies between 9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, while HCO+ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO+/HCN and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of moderate densities of n = 3x10^3 - 3x10^4 cm^-3, while the FUV field strength only has a small effect on the modeled line ratios. The line ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor role in influencing these line ratios.

preprint2013arXiv

Gas and dust cooling along the major axis of M33 (HerM33es): ISO/LWS CII observations

We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%.

preprint2013arXiv

The thermal dust emission in the N158-N159-N160 (LMC) star forming complex mapped by Spitzer, Herschel and LABOCA

We present a study of the infrared/submm emission of the LMC star forming complex N158-N159-N160. Combining observations from the Spitzer Space Telescope (3.6-70um), the Herschel Space Observatory (100-500um) and LABOCA (870um) allows us to work at the best angular resolution available now for an extragalactic source. We observe a remarkably good correlation between SPIRE and LABOCA emission and resolve the low surface brightnesses emission. We use the Spitzer and Herschel data to perform a resolved Spectral Energy Distribution (SED) modelling of the complex. Using MBB, we derive a global emissivity index beta_c of 1.47. If beta cold is fixed to 1.5, we find an average temperature of 27K. We also apply the Galliano et al. (2011) modelling technique (and amorphous carbon to model carbon dust) to derive maps of the star formation rate, the mean starlight intensity, the fraction of PAHs or the dust mass surface density of the region. We observe that the PAH fraction strongly decreases in the HII regions. This decrease coincides with peaks in the mean radiation field intensity map. The dust surface densities follow the FIR distribution, with a total dust mass of 2.1x10^4 Msolar (2.8 times less than when using graphite grains) in the resolved elements we model. We find a non-negligible amount of dust in the molecular cloud N159 South (showing no massive SF). We also investigate the drivers of the Herschel/PACS and SPIRE submm colours as well as the variations in the gas-to-dust mass ratio (G/D) and the XCO conversion factor in the region N159. We finally model individual regions to analyse variations in the SED shape across the complex and the 870um emission in more details. No measurable submm excess emission at 870um seems to be detected in these regions.

preprint2012arXiv

[FeII] as a tracer of Supernova Rate in Near-by Starburst Galaxies

Supernovae play an integral role in the feedback of processed material into the ISM of galaxies and are responsible for much of the chemical enrichment of the universe. The rate of supernovae can also reveal the star formation histories. Supernova rates are usually measured through the non-thermal radio continuum luminosity; however, a correlation between near-infrared [FeII] emission and supernova remnants has also been noted. We aim to find a quantitative relationship between the [FeII] at 1.26 um ([FeII]$_{1.26}$) luminosity and supernova rate in a sample of 11 near-by starburst galaxy centers. We perform a pixel-pixel analysis of this correlation on SINFONI data cubes. Using Br$γ$ equivalent width and luminosity as the only observational inputs into the Starburst 99 model, we derive the supernova rate at each pixel and thus create maps of supernova rates. We then compare these morphologically and quantitatively to the [FeII]$_{1.26}$ luminosity. We have found that a strong linear and morphological correlation exists between supernova rate and [FeII]$_{1.26}$ on a pixel-pixel basis: \[ log\frac{ν_{SNrate}}{yr^{-1}pc^{-2}} = 1.01 \pm 0.2\ast log\frac{[FeII]_{1.26}}{erg s^{-1}pc^{-2}} - 41.17 \pm 0.9\] This relation is valid for normal star forming galaxies but breaks down for extreme ultra luminous galaxies. The supernova rates derived from the Starburst 99 model are in good agreement with the radio-derived supernova rates, which underlines the strength of using [FeII] emission as a tracer of supernova rate. With the strong correlation found in this sample of galaxies, we conclude that [FeII]$_{1.26}$ emission can be generally used to derive accurate supernova rates on either a pixel-pixel or integrated galactic basis.

preprint2012arXiv

[FeII] as a tracer supernova rate

Supernovae play an integral role in the feedback of processed material into the ISMof galaxies and are responsible for most of the chemical enrichment of the universe. The rate of supernovae can also reveal the star formation histories. Supernova rate is usually measured through the non-thermal radio continuum luminosity, but in this paper we establish a quantitative relationship between the [FeII]1.26 luminosity and supernova rate in a sample of 11 near-by starburst galaxies. SINFONI data cubes are used to perform a pixel pixel analysis of this correlation. Using Br equivalent width and luminosity as the only observational inputs into Starburst 99, the supernova rate is derived at each pixel and a map of supernova rate is created. This is then compared morphologically and quantitatively to [FeII]1.26 luminosity map. We find a strong linear and morphological correlation between supernova rate and [FeII]1.26 on a pixel-pixel basis: log SNrate yr-1 pc-2 = (1.01 \pm 0.2) \ast log[FeII]1.26 ergs-1 pc-2 - 41.17 \pm 0.9 The Starburst 99 derived supernova rates are also in good agreement with the radio derived supernova rates, which further demonstrates the strength of [FeII] as a tracer of supernova rate. With the strong correlation found in this sample of galaxies, we now qualitatively use [FeII]1.26 to derive supernova rate on either a pixel-pixel or integrated galactic basis.

preprint2012arXiv

Diagnostics of the Molecular Component of PDRs with Mechanical Heating

Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy. Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. Methods. Clouds are modelled as 1D semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR-XDR code. Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and H2 O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease by a factor of at least two in high density regions with n \sim 105 cm-3, whereas that of HCN/HCO+ shows a strong dependence on mechanical heating for this same density range, with boosts of up to three orders of magnitude. Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the star formation rate to the gas density is close to, or exceeds, 7 \times 10-6 M yr-1 cm4.5 . If mechanical heating is not included, predicted column densities are underestimated, sometimes even by a few orders of magnitude. As a lower bound to its importance, we determined that it has non-negligible effects already when mechanical heating is as little as 1% of the UV heating in a PDR.

preprint2012arXiv

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

preprint2012arXiv

Spectrally resolved CII emission in M~33 ({\tt HerM33es}): Physical conditions and kinematics around BCLMP 691

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

preprint2012arXiv

The JCMT Nearby Galaxies Legacy Survey VIII. CO data and the L(CO3-2)-L(FIR) correlation in the SINGS sample

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

preprint2011arXiv

C+ Emission from the Magellanic Clouds II. [CII] maps of star-forming regions LMC-N 11, SMC-N 66, and several others

We study the 158 micron [CII] fine-structure line emission from star-forming regions as a function of metallicity. We have measured and mapped the [CII] emission from the very bright HII region complexes N 11 in the LMC and N 66 in the SMC, as well as the SMC HII regions N 25, N 27, N 83/N 84, and N 88, with the FIFI instrument on the Kuiper Airborne Observatory. In both the LMC and SMC, the ratio of the [CII] line to the CO line and to the far-infrared continuum emission is much higher than seen almost anywhere else, including Milky Way star-forming regions and whole galaxies. In the low metallicity, low dust-abundance environment of the LMC and the SMC, UV mean free path lengths are much greater than those in the higher-metallicity Milky Way. The increased photoelectric heating efficiencies cause significantly greater relative [CII] line emission strengths. At the same time, similar decreases in PAH abundances have the opposite effect, by diminishing photoelectric heating rates. Consequently, in low-metallicity environments the relative [CII] strengths are high but exhibit little further dependence on actual metallicity. Relative [CII] strengths are slightly higher in the LMC than in the SMC, which has both lower dust and lower PAH abundances.

preprint2010arXiv

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

preprint2010arXiv

Cold dust clumps in dynamically hot gas

We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.

preprint2010arXiv

Determining dust temperatures and masses in the Herschel era: the importance of observations longward of 200 micron

The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.

preprint2010arXiv

Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE): the Large Magellanic Cloud dust

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 μm and SPIRE 250, 350 and 500 μm. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 μm emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

preprint2010arXiv

Herschel observations of water vapour in Markarian 231

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

preprint2010arXiv

PACS and SPIRE photometer maps of M33: First results of the Herschel M33 extended survey (HERM33ES)

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

preprint2010arXiv

Sub-millimeter to centimeter excess emission from the Magellanic Clouds. I. Global spectral energy distribution

In order to reconstruct the global SEDs of the Magellanic Clouds over eight decades in spectral range, we combined literature flux densities representing the entire LMC and SMC respectively, and complemented these with maps extracted from the WMAP and COBE databases covering the missing the 23--90 GHz (13--3.2 mm) and the poorly sampled 1.25--250 THz (240--1.25 micron). We have discovered a pronounced excess of emission from both Magellanic Clouds, but especially the SMC, at millimeter and sub-millimeter wavelengths. We also determined accurate thermal radio fluxes and very low global extinctions for both LMC and SMC. Possible explanations are briefly considered but as long as the nature of the excess emission is unknown, the total dust masses and gas-to-dust ratios of the Magellanic Clouds cannot reliably be determined.

preprint2010arXiv

The Herschel revolution: unveiling the morphology of the high mass star formation sites N44 and N63 in the LMC

We study the structure of the medium surrounding sites of high-mass star formation to determine the interrelation between the HII regions and the environment from which they were formed. The density distribution of the surroundings is key in determining how the radiation of the newly formed stars interacts with the surrounds in a way that allows it to be used as a star formation tracer. We present new Herschel/SPIRE 250, 350 and 500 mum data of LHA 120-N44 and LHA 120-N63 in the LMC. We construct average spectral energy distributions (SEDs) for annuli centered on the IR bright part of the star formation sites. The annuli cover ~10-~100 pc. We use a phenomenological dust model to fit these SEDs to derive the dust column densities, characterise the incident radiation field and the abundance of polycyclic aromatic hydrocarbon molecules. We see a factor 5 decrease in the radiation field energy density as a function of radial distance around N63. N44 does not show a systematic trend. We construct a simple geometrical model to derive the 3-D density profile of the surroundings of these two regions. Herschel/SPIRE data have proven very efficient in deriving the dust mass distribution. We find that the radiation field in the two sources behaves very differently. N63 is more or less spherically symmetric and the average radiation field drops with distance. N44 shows no systematic decrease of the radiation intensity which is probably due to the inhomogeneity of the surrounding molecular material and to the complex distribution of several star forming clusters in the region.

preprint2010arXiv

The JCMT Nearby Galaxies Legacy Survey IV. Velocity Dispersions in the Molecular Interstellar Medium in Spiral Galaxies

An analysis of large-area CO J=3-2 maps from the James Clerk Maxwell Telescope for 12 nearby spiral galaxies reveals low velocity dispersions in the molecular component of the interstellar medium. The three lowest luminosity galaxies show a relatively flat velocity dispersion as a function of radius while the remaining nine galaxies show a central peak with a radial fall-off within 0.2-0.4 r(25). Correcting for the average contribution due to the internal velocitydispersions of a population of giant molecular clouds, the average cloud-cloud velocity dispersion across the galactic disks is 6.1 +/- 1.0 km/s (standard deviation 2.9 km/s), in reasonable agreement with previous measurements for the Galaxy andM33. The cloud-cloud velocity dispersion derived from the CO data is on average two times smaller than the HI velocity dispersion measured in the same galaxies. The low cloud-cloudvelocity dispersion implies that the molecular gas is the critical component determining the stability of the galactic disk against gravitational collapse, especially in those regions of the disk which are H2 dominated. The cloud-cloud velocity dispersion shows a significant positivecorrelation with both the far-infrared luminosity, which traces the star formation activity, and the K-band absolute magnitude, which traces the total stellar mass. For three galaxies in the Virgo cluster, smoothing the data to a resolution of 4.5 kpc (to match the typical resolution of high redshift CO observations) increases the measured velocity dispersion by roughly a factor of two, comparable to the dispersion measured recently in a normal galaxy at z=1. This comparison suggests that the mass and star formation rate surface densities may be similar in galaxies from z=0-1 and that the high star formation rates seen at z=1 may be partly due to the presence of physically larger molecular gas disks.

preprint2010arXiv

The JCMT Nearby Galaxies Legacy Survey VI: The distribution of gas and star formation in M81

We present a first complete 12 CO J=3-2 map of M81, observed as part of the Nearby Galaxies Legacy Survey being carried out at the James Clerk Maxwell Telescope. We detect 9 regions of significant CO emission located at different positions within the spiral arms, and confirm that the global CO emission in the galaxy is low. We combine these data with a new H-alpha map obtained using the Isaac Newton Telescope and archival HI, 24 microns and FUV images to uncover a correlation between the molecular gas and star forming regions in M81. For the nine regions detected in CO J=3-2, we combine our CO J=3-2 data with existing CO J=1-0 data to calculate line ratios. We find that the ratio J=(3-2)/(1-0) is in agreement with the range of typical values found in the literature (0.2-0.8). Making reasonable assumptions, this allows us to constrain the hydrogen density to the range (10^3-10^4) cm^{-3}. We also estimated the amount of hydrogen produced in photo-dissociation regions near the locations where CO J=3-2 was detected.

preprint2010arXiv

The Molecular Interstellar Medium of the Local Group Dwarf NGC6822

Do molecular clouds collapse to form stars at the same rate in all environments? In large spiral galaxies, the rate of transformation of H2 into stars (hereafter SFE) varies little. However, the SFE in distant objects (z~1) is much higher than in the large spiral disks that dominate the local universe. Some small local group galaxies share at least some of the characteristics of intermediate-redshift objects, such as size or color. Recent work has suggested that the Star Formation Efficiency (SFE, defined as the SFRate per unit H2) in local Dwarf galaxies may be as high as in the distant objects. A fundamental difficulty in these studies is the independent measure of the H2 mass in metal-deficient environments. At 490 kpc, NGC6822 is an excellent choice for this study; it has been mapped in the CO(2-1) line using the multibeam receiver HERA on the 30 meter IRAM telescope, yielding the largest sample of giant molecular clouds (GMCs) in this galaxy. Despite the much lower metallicity, we find no clear difference in the properties of the GMCs in NGC 6822 and those in the Milky Way except lower CO luminosities for a given mass. Several independent methods indicate that the total H2 mass in NGC 6822 is about 5 x 10^6 Msun in the area we mapped and less than 10^7 Msun in the whole galaxy. This corresponds to a NH2/ICO ~ 4 x 10^{21} cm^-2 /(Kkm/s) over large scales, such as would be observed in distant objects, and half that in individual GMCs. No evidence was found for H2 without CO emission. Our simulations of the radiative transfer in clouds are entirely compatible with these NH2/ICO values. The SFE implied is a factor 5 - 10 higher than what is observed in large local universe spirals.

preprint2009arXiv

CI and CO in nearby galaxy centers. The star-burst galaxies NGC 278, NGC 660, NGC 3628, NGC 4631, and NGC 4666

Maps and measurements of the J=1-0, J=2-1, J=3-2, J=4-3 12CO, the J=1-0, J=2-1 and J=3-2 13CO lines in the central arcminute squared of NGC 278, NGC 660, NGC 3628, NGC 4631, and NGC 4666, as well as 492 GHz [CI] maps in three of these are used to model the molecular gas. All five objects exhibit bright CO emission in the inner regions, with strong central concentrations in NGC 660, NGC 3628, and NGC 4666, but not in the weakest CO emitters NGC 278 and NGC 4631. In all cases, the observed lines could be modeled only with at least two distinct gas components. The physical condition of the molecular gas is found to differ from galaxy to galaxy. Relatively tenuous (density 100-1000 cm-3) and high kinetic temperature (100-150 K) gas occurs in all galaxies, except perhaps NGC 3628, and is mixed with cooler (10-30 K) and denser (3000-10000 cm-3) gas. In all galaxy centers, the CO-to-H2 conversion factor X is typically an order of magnitude smaller than the `standard' value for the Solar Neighborhood. The molecular gas is constrained within radii between 0.6 and 1.5 kpc from the nuclei. Within these radii, H2 masses are typically 0.6-1.5 x 10**8 M(O), which corresponds to no more than a few per cent of the dynamical mass in the same region.

preprint2009arXiv

Dense Star-forming Gas and Dust in the Magellanic Clouds

The early stages of star formation are closely related to the ambient conditions in the interstellar medium. Important questions such as dust abundance, size distribution, temperature distribution, fraction of molecular gas, fraction of dense gas, gas surface density and total amount of gas and dust require separation of metallicity and radiation effects. The Magellanic Clouds provide an ideal laboratory to carry out such studies. They are prominent targets for space observatories (Spitzer, Herschel), but an important role remains for large groundbased facilities, such as a 25 m class sub-millimeter telescope on Dome C.

preprint2009arXiv

Disentangling the circumnuclear environs of Centaurus A: I. High resolution molecular gas imaging

We present high resolution images of the 12CO(2-1) emission in the central 1' (1 kpc) of NGC 5128 (Centaurus A), observed using the SMA. We elucidate for the first time the distribution and kinematics of the molecular gas in this region with a resolution of 6'.0 x 2'.4 (100 pc x 40 pc). We spatially resolve the circumnuclear molecular gas in the inner 24'' x 12'' (400 pc x 200 pc), which is elongated along a position angle P.A. = 155 deg and perpendicular to the radio/X-ray jet. The SE and NW components of the circumnuclear gas are connected to molecular gas found at larger radii. This gas appears as two parallel filaments at P.A. = 120 deg, which are coextensive with the long sides of the 3 kiloparsec parallelogram shape of the previously observed dust continuum, as well as ionized and pure rotational H2 lines. Spatial and kinematical asymmetries are apparent in both the circumnuclear and outer gas, suggesting non-coplanar and/or non-circular motions. We extend to inner radii (r < 200 pc) previously studied warped disk models built to reproduce the central parallelogram-shaped structure. Adopting the warped disk model we would confirm a gap in emission between the radii r = 200 - 800 pc (12'' - 50''), as has been suggested previously. Although this model explains this prominent feature, however, our 12CO(2-1) observations show relevant deviations from this model. Namely, the physical connection between the circumnuclear gas and that at larger radii, brighter SE and NW sides on the parallelogram-shaped feature, and an outer curvature of its long sides. Overall it resembles more closely an S-shaped morphology, a trend that is also found in other molecular species. Hence, we explore qualitatively the possible contribution of a weak bi-symmetric potential which would naturally explain these peculiarities.

preprint2009arXiv

The JCMT Nearby Galaxies Legacy Survey III: Comparisons of cold dust, polycyclic aromatic hydrocarbons, molecular gas, and atomic gas in NGC 2403

We used 3.6, 8.0, 70, 160 micron Spitzer Space Telescope data, James Clerk Maxwell Telescope HARP-B CO J=(3-2) data, National Radio Astronomy Observatory 12 meter telescope CO J=(1-0) data, and Very Large Array HI data to investigate the relations among PAHs, cold (~20 K) dust, molecular gas, and atomic gas within NGC 2403, an SABcd galaxy at a distance of 3.13 Mpc. The dust surface density is mainly a function of the total (atomic and molecular) gas surface density and galactocentric radius. The gas-to-dust ratio monotonically increases with radius, varying from ~100 in the nucleus to ~400 at 5.5 kpc. The slope of the gas-to-dust ratio is close to that of the oxygen abundance, suggesting that metallicity strongly affects the gas-to-dust ratio within this galaxy. The exponential scale length of the radial profile for the CO J=(3-2) emission is statistically identical to the scale length for the stellar continuum-subtracted 8 micron (PAH 8 micron) emission. However, CO J=(3-2) and PAH 8 micron surface brightnesses appear uncorrelated when examining sub-kpc sized regions.

preprint2006arXiv

Irradiated ISM: Discriminating between Cosmic Rays and X-rays

The ISM of active galaxy centers is exposed to a combination of cosmic ray, FUV and X-ray radiation. We apply PDR models to this ISM with both `normal' and highly elevated (5\times 10^{-15}s^-1) cosmic-ray rates and compare the results to those obtained for XDRs. Our existing PDR-XDR code is used to construct models over a 10^3-10^5 cm^-3 density range and for 0.16-160 erg s^-1 cm^-2 impingent fluxes. We obtain larger high J (J>10) CO ratios in PDRs when we use the highly elevated cosmic ray rate, but these are always exceeded by the corresponding XDR ratios. The [CI] 609 mum/13CO(2-1) line ratio is boosted by a factor of a few in PDRs with n~10^3 cm^-3 exposed to a high cosmic ray rate. At higher densities ratios become identical irrespective of cosmic ray flux, while XDRs always show elevated [CI] emission per CO column. The HCN/CO and HCN/HCO+ line ratios, combined with high J CO emission lines, are good diagnostics to distinguish between PDRs under either low or high cosmic ray irradiation conditions, and XDRs. Hence, the HIFI instrument on Herschel, which can detect these CO lines, will be crucial in the study of active galaxies.

preprint2003arXiv

Results of the ESO-SEST Key Programme on CO in the Magellanic Clouds X. CO emission from star formation regions in LMC and SMC

We present J=1-0 and J=2-1 12CO maps of several star-forming regions in both the Large and the Small Magellanic Cloud, and briefly discuss their structure. Many of the detected molecular clouds are relatively isolated and quite small with dimensions of typically 20 pc. Some larger complexes have been detected, but in all cases the extent of the molecular clouds sampled by CO emission is significantly less than the extent of the ionized gas of the star-formation region. Very little diffuse extended CO emission was seen; diffuse CO in between or surrounding the detected discrete clouds is either very weak or absent. The majority of all LMC lines of sight detected in 13CO has an isotopic emission ratio I(12CO)/I(13CO) of about 10, i.e. twice higher than found in Galactic star-forming complexes. At the lowest 12CO intensities, the spread of isotopic emission ratios rapidly increases, low ratios representing relatively dense and cold molecular gas and high ratios marking CO photo-dissociation at cloud edges.

preprint2002arXiv

Results of the ESO-SEST Key Programme on CO in the Magellanic Clouds. IX. The giant LMC HII region complex N11

We present maps and a catalogue containing the J=1-0 12CO parameters of 29 individual molecular clouds in the second-brightest LMC star formation complex, N11. In the southwestern part of N11, molecular clouds occur in a ring or shell surrounding the major OB star association LH9. In the northeastern part, a chain of molecular clouds delineates the rim of one of the so-called supergiant shells in the LMC. The well-defined clouds have dimensions close to those of the survey beam (diameters of 25 pc or less). Some of the clouds were also observed in J=2-1 12CO, and in the lower two transitions of 13CO. Clouds mapped with a twice higher angular resolution in J=2-1 12CO show substructure with dimensions once again comparable to those of the mapping beam. The few clouds for which we could model physical parameters have fairly warm (T(kin) = 60 - 150 K) and moderately dense (n(H2) = 3000 cm-3) gas. The northeastern chain of CO clouds, although lacking in diffuse intercloud emission, is characteristic of the more quiescent regions of the LMC and appears to have been subject to relatively little photo-processing. The clouds forming part of the southwestern shell or ring, however, are almost devoid of diffuse intercloud emission and also exhibit other characteristics of an extreme photon-dominated region (PDR).

preprint2002arXiv

The Peculiar Nebula Simeis 57: I. Ionized Gas and Dust Extinction

We present high resolution radio continuum maps of the Galactic nebula Simeis 57 (= HS 191 = G 80.3+4.7) made with the WSRT and the DRAObservatory at frequencies of 609, 1412 and 1420 MHz. At optical and at radio wavelengths, the nebula has a peculiar ``S'' shape, crossed by long, thin and straight filaments. The radio maps, combined with other maps from existing databases, show essentially all radio emission from the peculiar and complex nebula to be thermal and optically thin. Although neither the distance nor the source of excitation of Simeis 57 are known, the nebula can only have a moderate electron density of typically n(e) = 100 cm-3. Its mass is also low, not exceeding some tens of solar masses. Peak emission measures are 5000 pc cm-6.Obscuring dust is closely associated with the nebula, but seems to occur mostly in front of it. Extinctions vary from A(V) = 1.0 mag to A(V) = 2.8 mag with a mean of about 2 mag. The extinction and the far-infrared emission at 100 microns are well-correlated.