Large Velocity Gradient Model Research Articles

Molecular gas excitation in the circumgalactic medium of MACS1931–26

The evolution of galaxies is largely affected by exchanging material with their close environment, the circumgalactic medium (CGM). In this work, we investigate the CGM and the interstellar medium (ISM) of the bright central galaxy (BCG) of the galaxy cluster, MACS1931−26 at z ∼ 0.35. We detected [CI](2−1), 12CO(1−0), and 12CO(7−6) emission lines with the APEX 12-m and NRO 45-m telescopes. We complemented these single-dish observations with 12CO(1−0), 12CO(3−2), and 12CO(4−3) ALMA interferometric data and inferred the cold molecular hydrogen physical properties. Using a modified large velocity gradient (LVG) model, we modelled the CO and CI emission of the CGM and BCG to extract the gas thermodynamical properties, including the kinetic temperature, the density, and the virialisation factor. Our study shows that the gas in the BCG is highly excited, comparable to the gas in local ultra luminous infrared galaxies (ULIRGs), while the CGM is likely less excited, colder, less dense, and less bound compared to the ISM of the BCG. The molecular hydrogen mass of the whole system derived using [CI](2−1) is larger than the mass derived from 12CO(1−0) in literature, showing that part of the gas in this system is CO-poor. Additional spatially resolved CI observations in both transitions, [CI](1−0) and [CI](2−1), and the completion of the CO SLED with higher CO transitions are crucial to trace the different phases of the gas in such systems and constrain their properties.

Doubly substituted isotopologues of HCCCN in TMC-1: Detection of D13CCCN, DC13CCN, DCC13CN, DCCC15N, H13C13CCN, H13CC13CN, HC13C13CN, HCC13C15N, and HC13CC15N

We report the first detection in space of a complete sample of nine doubly substituted isotopologues of HCCCN towards the cyanopolyyne peak of TMC-1 using observations of the QUIJOTE1 line survey taken with the Yebes 40 m telescope. We detected D13CCCN, DC13CCN, DCC13CN, DCCC15N, H13C13CCN, H13CC13CN, HC13C13CN, HCC13C15N, and HC13CC15N through their J = 4 − 3 and J = 5 − 4 lines in the 7 mm window. In addition, we present an extensive analysis of the emission of HCCCN and its singly substituted isotopologues through a large velocity gradient model of the lines detected at 7 mm and 3 mm using the Yebes 40 m and the IRAM 30 m telescopes, respectively. The derived column densities for all the isotopologues are consistent in the two spectral bands for an H2 volume density of 1 × 104 cm−3 and a kinetic temperature of 10 K. Whereas we observed a 13C fractionation for HCC13CN and other double isotopologues with a 13C atom adjacent to the nitrogen atom, we derived similar C/13C abundance ratios for the three 13C substituted species of DCCCN. This suggests additional chemical discrimination for deuterated isotopologues of HCCCN. Finally, we present the spatial distribution of the J = 4 − 3 and J = 5 − 4 lines from the singly substituted species observed with the Yebes 40 m telescope. The emission peak of the spatial distribution of DCCCN appears to be displaced by ∼40″ with respect to that of HCCCN and the 13C and 15N isotopologues. In addition to a different formation route for the deuterated species, we could also expect that this differentiation owing to the deuterium fractionation is more efficient at low temperatures, and therefore, that deuterated species trace a colder region of the cloud.

Updated Inventory of Carbon Monoxide in the Taurus Molecular Cloud

The most extensive survey of carbon monoxide (CO) gas in the Taurus molecular cloud relied on 12CO and 13CO J = 1 → 0 emission only, distinguishing the region where 12CO is detected without 13CO (named mask 1 region) from the one where both are detected (mask 2 region) (Goldsmith et al. ; Pineda et al. ). We have taken advantage of recent 12CO J = 3 → 2 James Clerk Maxwell Telescope observations, where they include mask 1 regions to estimate density, temperature, and N(CO) with a large velocity gradient model. This represents 1395 pixels out of ∼1.2 million in the mark 1 region. Compared to Pineda et al. () results and assuming a T kin of 30 K, we find a higher volume density of molecular hydrogen of 3.3 × 103 cm−3, compared to their 250–700 cm−3, and a CO column density of 5.7 × 1015 cm−2, about a quarter of their value. The differences are important and show the necessity to observe several CO transitions to better describe the intermediate region between the dense cloud and the diffuse atomic medium. Future observations to extend the 12CO J = 3 → 2 mapping further away from the 13CO-detected region comprising mask 1 are needed to revisit our understanding of the diffuse portions of dark clouds.

Molecular gas properties of Q1700-MD94: A massive main-sequence galaxy at z ≈ 2

We use a combination of new NOrthern Extended Millimeter Array (NOEMA) observations of the pair of [CI] transitions, the CO(7-6) line, and the dust continuum, in addition to ancillary CO(1-0) and CO(3-2) data, to study the molecular gas properties of Q1700-MD94. This is a massive, main-sequence galaxy at z ≈ 2. We find that for a reasonable set of assumptions for a typical massive star-forming galaxy, the CO(1-0), the [CI](1-0) and the dust continuum yield molecular gas masses that are consistent within a factor of ∼2. The global excitation properties of the molecular gas as traced by the [CI] and CO transitions are similar to those observed in other massive star-forming galaxies at z ∼ 2. Our large velocity gradient modeling using RADEX of the CO and [CI] spectral line energy distributions suggests the presence of relatively warm (Tkin = 41 K), dense (nH2 = 8 × 103 cm−3) molecular gas, comparable to the high-excitation molecular gas component observed in main-sequence star-forming galaxies at z ∼ 1. The galaxy size in the CO(1-0) and CO(7-6) line emission is comparable, which suggests that the highly excited molecular gas is distributed throughout the disk, powered by intense star formation activity. A confirmation of this scenario will require spatially resolved observations of the CO and [CI] lines, which can now be obtained with NOEMA upgraded capabilities.

Fast outflows in protoplanetary nebulae and young planetary nebulae observed by Herschel/HIFI

Context. Fast outflows and their interaction with slow shells (generally known as the fossil circumstellar envelope of asymptotic giant branch stars) play an important role in the structure and kinematics of protoplanetary and planetary nebulae (pPNe, PNe). To properly study their effects within these objects, we also need to observe the intermediate-temperature gas, which is only detectable in the far-infrared and submillimetre (submm) transitions. Aims. We study the physical conditions of the outflows presented in a number of pPNe and PNe, with a focus on their temperature and excitation states. Methods. We carried out Herschel/HIFI observations in the submm lines of 12CO in nine pPNe and nine PNe and complemented them with low-J CO spectra obtained with the IRAM 30m telescope and taken from the literature. The spectral resolution of HIFI allows us to identify and measure the different nebular components in the line profiles. The comparison with large velocity gradient model predictions was used to estimate the physical conditions of the warm gas in the nebulae, such as excitation conditions, temperature, and density. Results. We found high kinetic temperatures for the fast winds of pPNe, typically reaching between 75 K and 200 K. In contrast, the high-velocity gas in the sampled PNe is colder, with characteristic temperatures between 25 K and 75 K, and it is found in a lower excitation state. We interpret this correlation of the kinetic temperature and excitation state of fast outflows with the amount of time elapsed since their acceleration (probably driven by shocks) as a consequence of the cooling that occurred during the pPN phase.

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: CO Excitation and Atomic Carbon in Star-forming Galaxies at z = 1–3

Abstract We investigate the CO excitation and interstellar medium (ISM) conditions in a cold gas mass-selected sample of 22 star-forming galaxies at z = 0.46–3.60, observed as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). Combined with Very Large Array follow-up observations, we detect a total of 34 CO transitions with J = 1 up to 8 (and an additional 21 upper limits, up to J = 10) and 6 and transitions (and 12 upper limits). The CO(2–1) and CO(3–2)-selected galaxies, at and 2.5, respectively, exhibit a range in excitation in their mid-J = 4, 5 and high-J = 7, 8 lines, on average lower than ( -brighter) BzK-color- and submillimeter-selected galaxies at similar redshifts. The former implies that a warm ISM component is not necessarily prevalent in gas mass-selected galaxies at . We use stacking and Large Velocity Gradient models to measure and predict the average CO ladders at z < 2 and z ≥ 2, finding and , respectively. From the models, we infer that the galaxies at z ≥ 2 have intrinsically higher excitation than those at z < 2. This fits a picture in which the global excitation is driven by an increase in the star formation rate surface density of galaxies with redshift. We derive a neutral atomic carbon abundance of , comparable to the Milky Way and main-sequence galaxies at similar redshifts, and fairly high densities (≥104 cm−3), consistent with the low-J CO excitation. Our results imply a decrease in the cosmic molecular gas mass density at z ≥ 2 compared to previous ASPECS measurements.

The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk

Context.Cloud-scale surveys of molecular gas reveal the link between giant molecular cloud properties and star formation across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the normal star formation process, while galaxy centers tend to harbor denser gas that exhibits more extreme star formation. At high resolution, however, molecular clouds with exceptional gas properties and star formation activity may also be observed in normal disk environments. In this paper we study the brightest cloud traced in CO(2–1) emission in the disk of nearby spiral galaxy NGC 628.Aims.We characterize the properties of the molecular and ionized gas that is spatially coincident with an extremely bright H IIregion in the context of the NGC 628 galactic environment. We investigate how feedback and large-scale processes influence the properties of the molecular gas in this region.Methods.High-resolution ALMA observations of CO(2–1) and CO(1−0) emission were used to characterize the mass and dynamical state of the “headlight” molecular cloud. The characteristics of this cloud are compared to the typical properties of molecular clouds in NGC 628. A simple large velocity gradient (LVG) analysis incorporating additional ALMA observations of13CO(1−0), HCO+(1−0), and HCN(1−0) emission was used to constrain the beam-diluted density and temperature of the molecular gas. We analyzed the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the H IIregion.Results.The unusually bright headlight cloud is massive (1 − 2 × 107 M⊙), with a beam-diluted density ofnH2 = 5 × 104cm−3based on LVG modeling. It has a low virial parameter, suggesting that the CO emission associated with this cloud may be overluminous due to heating by the H IIregion. A young (2 − 4 Myr) stellar population with mass 3 × 105 M⊙is associated.Conclusions.We argue that the headlight cloud is currently being destroyed by feedback from young massive stars. Due to the large mass of the cloud, this phase of the its evolution is long enough for the impact of feedback on the excitation of the gas to be observed. The high mass of the headlight cloud may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk and receives a sustained inflow of gas that can promote the mass growth of the cloud.

C i](1–0) and [C i](2–1) in Resolved Local Galaxies*

Abstract We present resolved [C i] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from J up = 1 to 7 to interpret what phase of the interstellar medium the [C i] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4–3) and [C i](2–1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [C i] and CO line emission using large-velocity gradient models combined with an empirical template. According to this modeling, the [C i](1–0) line is clearly dominated by the low-excitation component. We determine [C i] to molecular mass conversion factors for both the [C i](1–0) and [C i](2–1) lines, with mean values of α [C i](1−0) = 7.3 M ⊙ K−1 km−1 s pc−2 and α [C i](2−1) = 34 M ⊙ K−1 km−1 s pc−2 with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of to (derived using the CO(2–1) line) suggests that [C i](1–0) may be just as good a tracer of cold molecular gas as CO(2–1) in galaxies of this type. On the other hand, the wider spread of α [C i](2−1) and the tight relation found between [C i](2–1) and CO(4–3) suggest that much of the [C i](2–1) emission may originate in warmer molecular gas.

The water line emission and ortho-to-para ratio in the Orion Bar photon-dominated region

Context. The ortho-to-para ratio (OPR) of water in the interstellar medium (ISM) is often assumed to be related to the formation temperature of water molecules, making it a potentially interesting tracer of the thermal history of interstellar gas. Aims. A very low OPR of 0.1–0.5 was previously reported in the Orion Bar photon-dominated region (PDR), based on observations of two optically thin H218O lines which were analyzed by using a single-slab large velocity gradient (LVG) model. The corresponding spin temperature does not coincide with the kinetic temperature of the molecular gas in this UV-illuminated region. This was interpreted as an indication of water molecules being formed on cold icy grains which were subsequently released by UV photodesorption. Methods. A more complete set of water observations in the Orion Bar, including seven H216O lines and one H218O line, carried out using Herschel/HIFI instrument, was reanalyzed using the Meudon PDR code to derive gas-phase water abundance and the OPR. The model takes into account the steep density and temperature gradients present in the region. Results. The model line intensities are in good agreement with the observations assuming that water molecules formed with an OPR corresponding to thermal equilibrium conditions at the local kinetic temperature of the gas and when solely considering gas-phase chemistry and water gas-grain exchanges through adsorption and desorption. Gas-phase water is predicted to arise from a region deep into the cloud, corresponding to a visual extinction of AV ~ 9, with a H216O fractional abundance of ~2 × 10−7 and column density of (1.4 ± 0.8) × 1015 cm−2 for a total cloud depth of AV = 15. A line-of-sight average OPR of 2.8 ± 0.2 is derived. Conclusions. The observational data are consistent with a nuclear spin isomer repartition corresponding to the thermal equilibrium at a temperature of 36 ± 2 K, much higher than the spin temperature previously reported for this region and close to the gas kinetic temperature in the water-emitting gas.

Discovery of a New Class I Methanol Maser Transition at 266.8 GHz

Abstract We report the detection of a new class I methanol maser candidate from the 52−41 E transition (266.8 GHz). This methanol transition has been detected toward a nearby high-mass star-forming region G352.630-1.067 (distance ∼0.7 kpc), in Submillimeter Array (SMA) observations. The new candidate transition has a similar spatial distribution as the 42−31 E (218.4 GHz) and 8−1−70 E (229.7 GHz) transitions, which are known class I maser transitions. Thermal methanol emission in this source is confined to a central hot core, while the three class I maser transitions are detected in two additional regions. These two maser-only emission regions are clearly associated with shocked gas traced by 2 μm Ks-band and thermal v = 0, J = 5−4 SiO molecular emission. In contrast to the thermal methanol emission from the hot core, the three class I maser transitions show an positive trend in the rotation diagram for the two maser regions. Large velocity gradient modeling of the 266.8, 218.4, and 229.7 GHz transitions shows that the 266.8 GHz transition can be a maser for a wide range of conditions. The intensity ratios for the three methanol transitions detected in maser regions can be reproduced under conditions that are typical for class I methanol maser sites. These facts all support the hypothesis that the detected emission from the 266.8 GHz methanol (52−41 E) transition is masing.

Source Structure and Molecular Gas Properties from High-resolution CO Imaging of SPT-selected Dusty Star-forming Galaxies

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of high-J CO lines ($J_\mathrm{up}=6$, 7, 8) and associated dust continuum towards five strongly lensed, dusty, star-forming galaxies (DSFGs) at redshift $z = 2.7$-5.7. These galaxies, discovered in the South Pole Telescope survey, are observed at $0.2''$-$0.4''$ resolution with ALMA. Our high-resolution imaging coupled with the lensing magnification provides a measurement of the structure and kinematics of molecular gas in the background galaxies with spatial resolutions down to kiloparsec scales. We derive visibility-based lens models for each galaxy, accurately reproducing observations of four of the galaxies. Of these four targets, three show clear velocity gradients, of which two are likely rotating disks. We find that the reconstructed region of CO emission is less concentrated than the region emitting dust continuum even for the moderate-excitation CO lines, similar to what has been seen in the literature for lower-excitation transitions. We find that the lensing magnification of a given source can vary by 20-50% across the line profile, between the continuum and line, and between different CO transitions. We apply Large Velocity Gradient (LVG) modeling using apparent and intrinsic line ratios between lower-J and high-J CO lines. Ignoring these magnification variations can bias the estimate of physical properties of interstellar medium of the galaxies. The magnitude of the bias varies from galaxy to galaxy and is not necessarily predictable without high resolution observations.

HOW DENSE IS YOUR GAS? ON THE RECOVERABILITY OF LVG MODEL PARAMETERS

ABSTRACT We explore the recoverability of gas physical conditions with the Large Velocity Gradient (LVG) model, using the public code RADEX and the molecules HCN and CO. Examining a wide parameter range with a series of models of increasing complexity, we use both grid and Monte Carlo Markov Chain methods to recover the input conditions, and quantify the inherent and noise-induced uncertainties in the model results. We find that even with the benefit of generous assumptions, the LVG models struggle to recover any parameter better than to within half a dex, although we find no evidence of systemic offsets. Examining isotopologue lines, we demonstrate that it is always preferable to model the isotopologue abundance ratio as a free parameter due to large biases introduced in all other parameters when an incorrect ratio is assumed. Finally, we explore the effects of the background radiation temperature on CO and HCN line ratios, with an emphasis on the effect of the CMB at z > 4 ?> , and show that while the effect on the line ratios is minor, the effect on the spectral line energy distribution peak is significant and that the CO(1–0) line luminosity to H2 mass conversion factor (α CO) needs to be altered to account for the loss of contrast against the hotter CMB as redshift increases.

MODELING THE MOLECULAR GAS IN NGC 6240

We present the first observations of H$^{13}$CN$(1-0)$, H$^{13}$CO$^+(1-0)$ and SiO$(2-1)$ in NGC\,6240, obtained with the IRAM PdBI. Combining a Markov Chain Monte Carlo (MCMC) code with Large Velocity Gradient (LVG) modelling, and with additional data from the literature, we simultaneously fit three gas phases and six molecular species to constrain the physical condition of the molecular gas, including mass$-$luminosity conversion factors. We find $\sim10^{10}M_\odot$ of dense molecular gas in cold, dense clouds ($T_{\rm k}\sim10$\,K, $n_{{\rm H}_2}\sim10^6$\,cm$^{-3}$) with a volume filling factor $<0.002$, embedded in a shock heated molecular medium ($T_{\rm k}\sim2000$\,K, $n_{{\rm H}_2}\sim10^{3.6}$\,cm$^{-3}$), both surrounded by an extended diffuse phase ($T_{\rm k}\sim200$\,K, $n_{{\rm H}_2}\sim10^{2.5}$\,cm$^{-3}$). We derive a global $\alpha_{\rm CO}=1.5^{7.1}_{1.1}$ with gas masses $\log_{10}\left(M / [M_\odot]\right)=10.1_{10.0}^{10.8}$, dominated by the dense gas. We also find $\alpha_{\rm HCN} = 32^{89}_{13}$, which traces the cold, dense gas. The [$^{12}$C]/[$^{13}$C] ratio is only slightly elevated ($98^{230}_{65}$), contrary to the very high [CO]/[$^{13}$CO] ratio (300-500) reported in the literature. However, we find very high [HCN]/[H$^{13}$CN] and [HCO$^+$]/[H$^{13}$CO$^+$] abundance ratios $(300^{500}_{200})$ which we attribute to isotope fractionation in the cold, dense clouds.

Predicting the intensity mapping signal for multi-J CO lines

We present a novel approach to estimating the intensity mapping signal of any CO rotational line emitted during the Epoch of Reionization (EoR). Our approach is based on large velocity gradient (LVG) modeling, a radiative transfer modeling technique that generates the full CO spectral line energy distribution (SLED) for a specified gas kinetic temperature, volume density, velocity gradient, molecular abundance, and column density. These parameters, which drive the physics of CO transitions and ultimately dictate the shape and amplitude of the CO SLED, can be linked to the global properties of the host galaxy, mainly the star formation rate (SFR) and the SFR surface density. By further employing an empirically derived SFR−M relation for high redshift galaxies, we can express the LVG parameters, and thus the specific intensity of any CO rotational transition, as functions of the host halo mass M and redshift z. Integrating over the range of halo masses expected to host CO-luminous galaxies, we predict a mean CO(1-0) brightness temperature ranging from ∼ 0.6 μK at z = 6 to ∼ 0.03 μK at z = 10 with brightness temperature fluctuations of ΔCO2 ∼ 0.1 and 0.005 μK respectively, at k = 0.1 Mpc−1. In this model, the CO emission signal remains strong for higher rotational levels at z = 6, with ⟨ TCO ⟩ ∼ 0.3 and 0.05 μK for the CO J = 6arrow5 and CO J = 10arrow9 transitions respectively. Including the effects of CO photodissociation in these molecular clouds, especially at low metallicities, results in the overall reduction in the amplitude of the CO signal, with the low- and high-J lines weakening by 2–20% and 10–45%, respectively, over the redshift range 4 < z < 10.

SIMULTANEOUS OBSERVATION OF WATER AND CLASS I METHANOL MASERS TOWARD CLASS II METHANOL MASER SOURCES

We present a simultaneous single-dish survey of 22 GHz water maser and 44 GHz and 95 GHz class I methanol masers toward 77 6.7 GHz class II methanol maser sources, which were selected from the Arecibo methanol maser Galactic plane survey (AMGPS) catalog.Water maser emission is detected in 39 (51%) sources, of which 15 are new detections. Methanol maser emission at 44 GHz and 95 GHz is found in 25 (32%) and 19 (25%) sources, of which 21 and 13 sources are newly detected, respectively. We find 4 high-velocity (> 30 km/s) water maser sources, including 3 dominant blue- or redshifted outflows.The 95 GHz masers always appear with the 44 GHz maser emission. They are strongly correlated with 44 GHz masers in velocity, flux density, and luminosity, while they are not correlated with either water or 6.7 GHz class II methanol masers. The average peak flux density ratio of 95 GHz to 44 GHz masers is close to unity, which is two times higher than previous estimates. The flux densities of class I methanol masers are more closely correlated with the associated BGPS core mass than those of water or class II methanol masers. Using the large velocity gradient (LVG) model and assuming unsaturated class I methanol maser emission, we derive the fractional abundance of methanol to be in a range of 4.2*10^-8 to 2.3*10^-6, with a median value of 3.3\pm2.7*10^-7.

The Central Molecular Zone in H13CO+ and thermal SiO emission lines

Abstract We present results from a high-resolution wide-field imaging observation of the Central Molecular Zone (CMZ) in H13CO+J = 1–0 and SiO v = 0, J = 2–1 emission lines using the 45-m telescope at the Nobeyama Radio Observatory in order to depict the high-density molecular gas mass distribution and explore molecular gas affected by interstellar shocks. By using the clumpfind algorithm we identify 114 molecular clumps in the H13CO+ data cube. The clump mass function (CMF) can be described by a power law (dN/dM ∝ M−γ, with γ ≃ 2) in the range of 5 × 103 &amp;lt; Mclump &amp;lt; 2 × 104 M⊙. The brightness temperature ratio RT = T(SiO)/T(H13CO+) indicates the strength of the shock. The Sgr B2 complex has several expanding shell-like structures with a large velocity-width (ΔV &amp;gt; 50 km s−1). They correspond to the high brightness temperature ratio area (RT ≲ 8). A Large Velocity Gradient model analysis with the RADEX program confirms the enhancement of SiO molecules in the shells. On the other hand, there is no such large shell with a high ratio in the Sgr A complex. A hollow hemispherical structure is found in the l–b–v data cube of the Sgr B2 complex. The central part of the structure has a high ratio (RT ≲ 8). The ratio decreases with proximity to the outer area. The properties are consistent with those expected from simulations of cloud–cloud collisions. The structure is probably a site of ongoing cloud–cloud collision.

SI-BEARING MOLECULES TOWARD IRC+10216: ALMA UNVEILS THE MOLECULAR ENVELOPE OF CWLEO.

We report the detection of SiS rotational lines in high-vibrational states as well as SiO and SiC2 lines in their ground vibrational state toward IRC+10216 during the Atacama Large Millimeter Array Cycle 0. The spatial distribution of these molecules shows compact emission for SiS and a more extended emission for SiO and SiC2, and also proves the existence of an increase in the SiC2 emission at the outer shells of the circumstellar envelope. We analyze the excitation conditions of the vibrationally excited SiS using the population diagram technique, and we use a large velocity gradient model to compare with the observations. We found moderate discrepancies between the observations and the models that could be explained if SiS lines detected are optically thick. Additionally, the line profiles of the detected rotational lines in the high energy vibrational states show a decreasing linewidth with increasing energy levels. This may be evidence that these lines could be excited only in the inner shells, i.e., the densest and hottest, of the circumstellar envelope of IRC+10216.

HIGH-JCO SLEDs IN NEARBY INFRARED BRIGHT GALAXIES OBSERVED BYHERSCHEL/PACS

We report the detection of far-infrared (FIR) CO rotational emission from nearby active galactic nuclei (AGN) and starburst galaxies, as well as several merging systems and Ultra-Luminous Infrared Galaxies (ULIRGs). Using Herschel-PACS, we have detected transitions in the J$_{upp}$ = 14 - 20 range ($\lambda \sim$ 130 - 185 $\mu$m, $\nu \sim$ 1612 - 2300 GHz) with upper limits on (and in two cases, detections of) CO line fluxes up to J$_{upp}$ = 30. The PACS CO data obtained here provide the first well-sampled FIR extragalactic CO SLEDs for this range, and will be an essential reference for future high redshift studies. We find a large range in the overall SLED shape, even amongst galaxies of similar type, demonstrating the uncertainties in relying solely on high-J CO diagnostics to characterize the excitation source of a galaxy. Combining our data with low-J line intensities taken from the literature, we present a CO ratio-ratio diagram and discuss its potential diagnostic value in distinguishing excitation sources and physical properties of the molecular gas. The position of a galaxy on such a diagram is less a signature of its excitation mechanism, than an indicator of the presence (or absence) of warm, dense molecular gas. We then quantitatively analyze the CO emission from a subset of the detected sources with Large Velocity Gradient (LVG) radiative transfer models to fit the CO SLEDs. Using both single-component and two-component LVG models to fit the kinetic temperature, velocity gradient, number density and column density of the gas, we derive the molecular gas mass and the corresponding CO-to-H$_2$ conversion factor, $\alpha_{CO}$, for each respective source. For the ULIRGs we find $\alpha$ values in the canonical range 0.4 - 5 M$_\odot$/(K kms$^{-1}$pc$^2$), while for the other objects, $\alpha$ varies between 0.2 and 14.} Finally, we compare our best-fit LVG model ..

Modeling CO emission from hydrodynamic simulations of nearby spirals, starbursting mergers, and high-redshift galaxies

We model the intensity of emission lines from the CO molecule, based on hydrodynamic simulations of spirals, mergers, and high-redshift galaxies with very high resolutions (3pc and 10^3 Msun) and detailed models for the phase-space structure of the interstellar gas including shock heating, stellar feedback processes and galactic winds. The simulations are analyzed with a Large Velocity Gradient (LVG) model to compute the local emission in various molecular lines in each resolution element, radiation transfer and opacity effects, and the intensity emerging from galaxies, to generate synthetic spectra for various transitions of the CO molecule. This model reproduces the known properties of CO spectra and CO-to-H2 conversion factors in nearby spirals and starbursting major mergers. The high excitation of CO lines in mergers is dominated by an excess of high-density gas, and the high turbulent velocities and compression that create this dense gas excess result in broad linewidths and low CO intensity-to-H2 mass ratios. When applied to high-redshift gas-rich disks galaxies, the same model predicts that their CO-to-H2 conversion factor is almost as high as in nearby spirals, and much higher than in starbursting mergers. High-redshift disk galaxies contain giant star-forming clumps that host a high-excitation component associated to gas warmed by the spatially-concentrated stellar feedback sources, although CO(1-0) to CO(3-2) emission is overall dominated by low-excitation gas around the densest clumps. These results overall highlight a strong dependence of CO excitation and the CO-to-H2 conversion factor on galaxy type, even at similar star formation rates or densities. The underlying processes are driven by the interstellar medium structure and turbulence and its response to stellar feedback, which depend on global galaxy structure and in turn impact the CO emission properties.

CHEMICALLY DISTINCT NUCLEI AND OUTFLOWING SHOCKED MOLECULAR GAS IN Arp 220

We present the results of interferometric spectral line observations of Arp 220 at 3.5mm and 1.2mm from the Plateau de Bure Interferometer (PdBI), imaging the two nuclear disks in H$^{13}$CN$(1 - 0)$ and $(3 - 2)$, H$^{13}$CO$^+(1 - 0)$ and $(3 - 2)$, and HN$^{13}$C$(3 - 2)$ as well as SiO$(2 - 1)$ and $(6 - 5)$, HC$^{15}$N$(3 - 2)$, and SO$(6_6 - 5_5)$. The gas traced by SiO$(6 - 5)$ has a complex and extended kinematic signature including a prominent P Cygni profile, almost identical to previous observations of HCO$^+(3 - 2)$. Spatial offsets $0.1''$ north and south of the continuum centre in the emission and absorption of the SiO$(6 - 5)$ P Cygni profile in the western nucleus (WN) imply a bipolar outflow, delineating the northern and southern edges of its disk and suggesting a disk radius of $\sim40$ pc, consistent with that found by ALMA observations of Arp 220. We address the blending of SiO$(6 - 5)$ and H$^{13}$CO$^+(3 - 2)$ by considering two limiting cases with regards to the H$^{13}$CO$^+$ emission throughout our analysis. Large velocity gradient (LVG) modelling is used to constrain the physical conditions of the gas and to infer abundance ratios in the two nuclei. Our most conservative lower limit on the [H$^{13}$CN]/[H$^{13}$CO$^+$] abundance ratio is 11 in the WN, cf. 0.10 in the eastern nucleus (EN). Comparing these ratios to the literature we argue on chemical grounds for an energetically significant AGN in the WN driving either X-ray or shock chemistry, and a dominant starburst in the EN.