Mapping CS in Starburst Galaxies: Disentangling and Characterising Dense Gas

Aims. We observe the dense gas tracer CS in two nearby starburst galaxies to determine how the conditions of the dense gas varies across the circumnuclear regions in starburst galaxies. Methods. Using the IRAM-30m telescope, we mapped the distribution of the CS(2-1) and CS(3-2) lines in the circumnuclear regions of the nearby starburst galaxies NGC 3079 and NGC 6946. We also detected the formaldehyde (H2CO) and methanol (CH3OH) in both galaxies. We marginally detect the isotopologue C34S. Results. We calculate column densities under LTE conditions for CS and CH3OH. Using the detections accumulated here to guide our inputs, we link a time and depth dependent chemical model with a molecular line radiative transfer model; we reproduce the observations, showing how conditions where CS is present are likely to vary away from the galactic centres. Conclusions. Using the rotational diagram method for CH3OH, we obtain a lower limit temperature of 14 K. In addition to this, by comparing the chemical and radiative transfer models to observations, we determine the properties of the dense gas as traced by CS (and CH3OH). We also estimate the quantity of the dense gas. We find that, provided that there are a between 10^5 and 10^6 dense cores in our beam, for both target galaxies, emission of CS from warm (T = 100 - 400 K), dense (n(H2) = 10^5-6 cm-3) cores, possibly with a high cosmic ray ionisation rate (zeta = 100 zeta0) best describes conditions for our central pointing. In NGC 6946, conditions are generally cooler and/or less dense further from the centre, whereas in NGC 3079, conditions are more uniform. The inclusion of shocks allows for more efficient CS formation, leading to an order of magnitude less dense gas being required to replicate observations in some cases.Comment: 14 pages, 10 figures, accepted to A&

Are 12^{12}CO lines good indicators of the star formation rate in galaxies?

In this paper, we investigate the relevance of using the $^{12}$CO line emissions as indicators of star formation rates (SFR). For the first time, we present this study for a relatively large number of $^{12}$CO transitions (12) as well as over a large interval in redshift (from z$\sim$0 to z$\sim$6). For the nearby sources (D$\leq$10 Mpc), we have used homogeneous sample of $^{12}$CO data provided by Bayet et al. (2004, 2006), mixing observational and modelled line intensities. For higher-z sources (z $\geq$ 1), we have collected $^{12}$CO observations from various papers and have completed the data set of line intensities with model predictions which we also present in this paper. Finally, for increasing the statistics, we have included recent $^{12}$CO(1-0) and $^{12}$CO(3-2) observations of intermediate-z sources. Linear regressions have been calculated for identifying the tightest SFR-$^{12}$CO line luminosity relationships. We show that the \emph{total} $^{12}$CO, the $^{12}$CO(5-4), the $^{12}$CO(6-5) and the $^{12}$CO(7-6) luminosities are the best indicators of SFR (as measured by the far-infrared luminosity). Comparisons with theoretical approaches from Krumholz and Thompson (2007) and Narayanan et al. (2008) are also performed in this paper. Although in general agreement, the predictions made by these authors and the observational results we present here show small and interesting discrepancies. In particular, the slope of the linear regressions, for J$_{upper}\geq$ 4 $^{12}$CO lines are not similar between theoretical studies and observations. On one hand, a larger high-J $^{12}$CO data set of observations might help to better agree with models, increasing the statistics. On the other hand, theoretical studies extended to high redshift sources might also reduce such discrepancies.Comment: 10 pages, 3 figures and 4 tables, Accepted in MNRA

Dense molecular gas toward W49A: A template for extragalactic starbursts?

The HCN, HCO+, and HNC molecules are commonly used as tracers of dense star-forming gas in external galaxies, but such observations are spatially unresolved. Reliably inferring the properties of galactic nuclei and disks requires detailed studies of sources whose structure is spatially resolved. We compare the spatial distributions and abundance ratios of HCN, HCO+, and HNC in W49A, the most massive and luminous star-forming region in the Galactic disk, based on maps of a 2' (6.6 pc) field at 14" (0.83 pc) resolution of the J=4-3 transitions of HCN, H13CN, HC15N, HCO+, H13CO+, HC18O+ and HNC. The kinematics of the molecular gas in W49A appears complex, with a mixture of infall and outflow motions. Both the line profiles and comparison of the main and rarer species show that the main species are optically thick. Two 'clumps' of infalling gas appear to be at ~40 K, compared to ~100 K at the source centre, and may be ~10x denser than the rest of the outer cloud. Chemical modelling suggests that the HCN/HNC ratio probes the current gas temperature, while the HCN/HCO+ ratio and the deuterium fractionation were set during an earlier, colder phase of evolution. The data suggest that W49A is an appropriate analogue of an extragalactic star forming region. Our data show that the use of HCN/HNC/HCO+ line ratios as proxies for the abundance ratios is incorrect for W49A, suggesting the same for galactic nuclei. Our observed isotopic line ratios such as H13CN/H13CO+ approach our modeled abundance ratios quite well in W49A. The 4-3 lines of HCN and HCO+ are much better tracers of the dense star-forming gas in W49A than the 1-0 lines. Our observed HCN/HNC and HCN/HCO+ ratios in W49A are inconsistent with homogeneous PDR or XDR models, indicating that irradiation hardly affects the gas chemistry in W49A. Overall, the W49A region appears to be a useful template for starburst galaxies.Comment: Accepted by A&A; 17 pages, 15 figure

The ATLAS3D project - XXV: Two-dimensional kinematic analysis of simulated galaxies and the cosmological origin of fast and slow rotators

We present a detailed two-dimensional stellar dynamical analysis of as ample of 44 cosmological hydrodynamical simulations of individual central galaxies with stellar masses of 2 x 1010Msun ∼≤ Mstar ∼≤ 6x 1011Msun. Kinematic maps of the stellar line-of-sight velocity, velocity dispersion, and higher-order Gauss-Hermite moments h3 and h4 are constructed for each central galaxy and for the most massive satellites. The amount of rotation is quantified using the λR-parameter. The velocity, velocity dispersion, h3, and h4 fields of the simulated galaxies show a diversity similar to observed kinematic maps of early-type galaxies in the ATLAS3D survey. This includes fast (regular), slow, and misaligned rotation, hot spheroids with embedded cold disk components as well as galaxies with counter-rotating cores or central depressions in the velocity dispersion. We link the present-day kinematic properties to the individual cosmological formation histories of the galaxies. In general, major galaxy mergers have a significant influence on the rotation properties resulting in both a spin-down as well as a spin-up of the merger remnant. Lower mass galaxies with significant in-situ formation of stars, or with additional gas-rich major mergers - resulting in a spin-up - in their formation history, form elongated fast rotators with a clear anti-correlation of h3 and v/σ. An additional formation path for fast rotators includes gas-poor major mergers leading to a spin-up of the remnants. This formation path does not result in anti-correlated h3 and v/σ. The galaxies most consistent with the rare class of non-rotating round early-type galaxies grow by gas-poor minor mergers alone. In general, more massive galaxies have less in-situ star formation since z ∼ 2, rotate slower and have older stellar populations. (shortened)PostprintPeer reviewe

Extended warm and dense gas towards W49A: starburst conditions in our Galaxy?

The star formation rates in starburst galaxies are orders of magnitude higher than in local star-forming regions, and the origin of this difference is not well understood. We use sub-mm spectral line maps to characterize the physical conditions of the molecular gas in the luminous Galactic star-forming region W49A and compare them with the conditions in starburst galaxies. We probe the temperature and density structure of W49A using H_2CO and HCN line ratios over a 2'x2' (6.6x6.6 pc) field with an angular resolution of 15" (~0.8 pc) provided by the JCMT Spectral Legacy Survey. We analyze the rotation diagrams of lines with multiple transitions with corrections for optical depth and beam dilution, and estimate excitation temperatures and column densities. Comparing the observed line intensity ratios with non-LTE radiative transfer models, our results reveal an extended region (about 1'x1', equivalent to ~3x3 pc at the distance of W49A) of warm (> 100 K) and dense (>10^5 cm^-3) molecular gas, with a mass of 2x10^4 - 2x10^5 M_Sun (by applying abundances derived for other regions of massive star-formation). These temperatures and densities in W49A are comparable to those found in clouds near the center of the Milky Way and in starburst galaxies. The highly excited gas is likely to be heated via shocks from the stellar winds of embedded, O-type stars or alternatively due to UV irradiation, or possibly a combination of these two processes. Cosmic rays, X-ray irradiation and gas-grain collisional heating are less likely to be the source of the heating in the case of W49A.Comment: Accepted for publication in A&A; 11 pages, 9 figure

The molecular gas in Luminous Infrared Galaxies I: CO lines, extreme physical conditions, and their drivers

We report results from a large molecular line survey of Luminous Infrared Galaxies (L_{IR} >= 10^{11} L_sol) in the local Universe (z<=0.1), conducted during the last decade with the James Clerk Maxwell Telescope (JCMT) and the IRAM 30-m telescope. This work presents the CO and {13}CO line data for 36 galaxies, further augmented by multi-J total CO luminosities available for other IR-bright galaxies from the literature. This yields a sample of N=70 galaxies with the star-formation (SF) powered fraction of their IR luminosities spanning L_{IR} (10^{10}-2x10^{12}) L_sol and a wide range of morphologies. Simple comparisons of their available CO Spectral Line Energy Distributions (SLEDs) with local ones, as well as radiative transfer models discern a surprisingly wide range of average ISM conditions, with most of the surprises found in the high-excitation regime. These take the form of global CO SLEDs dominated by a very warm (T_{kin}>=100 K) and dense (n>=10^4 cm^{-3}) gas phase, involving galaxy-sized (~(few)x10^9 M_sol) gas mass reservoirs under conditions that would otherwise amount only ~1% of mass per typical SF molecular cloud in the Galaxy. Some of the highest excitation CO SLEDs are found in the so-called Ultra Luminous Infrared Galaxies and seem irreducible to ensembles of ordinary SF-powered regions. Highly supersonic turbulence and high cosmic ray (CR) energy densities rather than far-UV/optical photons or SNR-induced shocks from individual SF sites can globally warm the large amounts of dense gas found in these merger-driven starbursts and easily power their extraordinary CO line excitation.....Comment: 29 pages, 12 Figures, 8 Tables, originally submitted and now accepted for publication in the Monthly Notices of the Royal Astronomical Society (only minor modifications with respect to the first version

Parameterizing the interstellar dust temperature

The temperature of interstellar dust particles is of great importance to astronomers. It plays a crucial role in the thermodynamics of interstellar clouds, because of the gas-dust collisional coupling. It is also a key parameter in astrochemical studies that governs the rate at which molecules form on dust. In 3D (magneto)hydrodynamic simulations often a simple expression for the dust temperature is adopted, because of computational constraints, while astrochemical modelers tend to keep the dust temperature constant over a large range of parameter space. Our aim is to provide an easy-to-use parametric expression for the dust temperature as a function of visual extinction ($A_{\rm V}$) and to shed light on the critical dependencies of the dust temperature on the grain composition. We obtain an expression for the dust temperature by semi-analytically solving the dust thermal balance for different types of grains and compare to a collection of recent observational measurements. We also explore the effect of ices on the dust temperature. Our results show that a mixed carbonaceous-silicate type dust with a high carbon volume fraction matches the observations best. We find that ice formation allows the dust to be warmer by up to 15% at high optical depths ($A_{\rm V}> 20$ mag) in the interstellar medium. Our parametric expression for the dust temperature is presented as $T_{\rm d} = \left[ 11 + 5.7\times \tanh\bigl( 0.61 - \log_{10}(A_{\rm V})\bigr) \right] \, \chi_{\rm uv}^{1/5.9}$, where $\chi_{\rm uv}$ is in units of the Draine (1978) UV fieldComment: 16 pages, 17 figures, 4 tables. Accepted for publication in A&A. Version 2: the omission of factor 0.921 in equation 4 is correcte

Exceptional AGN-driven turbulence inhibits star formation in the 3C 326N radio galaxy

We detect bright [CII]158μm line emission from the radio galaxy 3C 326N at z=0.09, which shows weak star formation (SFR⊙~yr−1) despite having strong H2 line emission and 2×109M⊙ of molecular gas. The [CII] line is twice as strong as the 0-0S(1) 17μm H2 line, and both lines are much in excess what is expected from UV heating. We combine infrared Spitzer and Herschel data with gas and dust modeling to infer the gas physical conditions. The [CII] line traces 30 to 50% of the molecular gas mass, which is warm (70−3. The [CII] line is broad with a blue-shifted wing, and likely to be shaped by a combination of rotation, outflowing gas, and turbulence. It matches the near-infrared H2 and the Na D optical absorption lines. If the wing is interpreted as an outflow, the mass loss rate would be larger than 20M⊙/yr, and the depletion timescale shorter than the orbital timescale (108yr). These outflow rates may be over-estimated because the stochastic injection of turbulence on galactic scales can contribute to the skewness of the line profile and mimic outflowing gas. We argue that the dissipation of turbulence is the main heating process of this gas. Cosmic rays can also contribute to the heating but they require an average gas density larger than the observational constraints. We show that strong turbulent support maintains a high gas vertical scale height (0.3-4kpc) in the disk and can inhibit the formation of gravitationally-bound structures at all scales, offering a natural explanation for the weakness of star formation in 3C 326N. To conclude, the bright [CII] line indicates that strong AGN jet-driven turbulence may play a key role in enhancing the amount of molecular gas (positive feedback) but yet can prevent star formation on galactic scales (negative feedback)

Photon dominated regions in the spiral arms of M83 and M51

We present CI 3P1-3P0 spectra at four spiral arm positions and the nuclei of the nearby galaxies M83 and M51 obtained at the JCMT. This data is complemented with maps of CO 1-0, 2-1, and 3-2, and ISO/LWS far-infrared data of CII (158 micron), OI (63 micron), and NII (122 micron) allowing for the investigation of a complete set of all major gas cooling lines. From the intensity of the NII line, we estimate that between 15% and 30% of the observed CII emission originate from the dense ionized phase of the ISM. The analysis indicates that emission from the diffuse ionized medium is negligible. In combination with the FIR dust continuum, we find gas heating efficiencies below ~0.21% in the nuclei, and between 0.25 and 0.36% at the outer positions. Comparison with models of photon-dominated regions (PDRs) of Kaufman et al. (1999) with the standard ratios OI(63)/CII_PDR and (OI(63)+CII_PDR) vs. TIR, the total infrared intensity, yields two solutions. The physically most plausible solution exhibits slightly lower densities and higher FUV fields than found when using a full set of line ratios, CII_PDR/CI(1-0), CI(1-0)/CO(3-2), CO(3-2)/CO(1-0), CII/CO(3-2), and, OI(63)/CII_PDR. The best fits to the latter ratios yield densities of 10^4 cm^-3 and FUV fields of ~G_0=20-30 times the average interstellar field without much variation. At the outer positions, the observed total infrared intensities are in perfect agreement with the derived best fitting FUV intensities. The ratio of the two intensities lies at 4-5 at the nuclei, indicating the presence of other mechanisms heating the dust

The initial conditions of star formation: cosmic rays as the fundamental regulators

Cosmic rays (CRs) control the thermal, ionization and chemical state of the dense H_2 gas regions that otherwise remain shielded from far-UV and optical stellar radiation propagating through the dusty ISM of galaxies. It is in such CR-dominated regions (CRDRs) rather than Photon-dominated regions (PDRs) of H_2 clouds where the star formation initial conditions are set, making CRs the ultimate star-formation feedback factor in galaxies, able to operate even in their most deeply dust-enshrouded environments. CR-controlled star formation initial conditions naturally set the stage for a near-invariant stellar Initial Mass Function (IMF) in galaxies as long as their average CR energy density U_{CR} permeating their molecular ISM remains within a factor of ~10 of its Galactic value. Nevertheless, in the extreme environments of the compact starbursts found in merging galaxies, where U_{CR}\sim(few)x10^{3}U_{CR,Gal}, CRs dramatically alter the initial conditions of star formation. In the resulting extreme CRDRs H_2 cloud fragmentation will produce far fewer low mass (<8 M_{sol}) stars, yielding a top-heavy stellar IMF. This will be a generic feature of CR-controlled star-formation initial conditions, lending a physical base for a bimodal IMF during galaxy formation, with a top-heavy one for compact merger-induced starbursts, and an ordinary IMF preserved for star formation in isolated gas-rich disks. In this scheme the integrated galactic IMFs (IGIMF) are expected to be strong functions of the star formation history of galaxies.Comment: 21 pages, 8 figures, Invited contribution to the Sant Cugat Forum on Astrophysics: Workshop on Cosmic-ray induced phenomenology in star-forming environments (April 16-19, 2012) (Version 3, with final, minor, corrections