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.Comment: 38 pages, to appear in A&

Evolution of the ISM in Luminous IR Galaxies

Molecules that trace the high-density regions of the interstellar medium may be used to evaluate the changing physical and chemical environment during the ongoing nuclear activity in (Ultra-)Luminous Infrared Galaxies. The changing ratios of the HCN(1-0), HNC(1-0), HCO+(1-0), CN(1-0) and CN(2-1), and CS(3-2) transitions were compared with the HCN(1-0)/CO(1-0) ratio, which is proposed to represent the progression time scale of the starburst. These diagnostic diagrams were interpreted using the results of theoretical modeling using a large physical and chemical network to describe the state of the nuclear ISM in the evolving galaxies. Systematic changes are seen in the line ratios as the sources evolve from early stage for the nuclear starburst (ULIRGs) to later stages. These changes result from changing environmental conditions and particularly from the lowering of the average density of the medium. A temperature rise due to mechanical heating of the medium by feedback explains the lowering of the ratios at later evolutionary stages. Infrared pumping may affect the CN and HNC line ratios during early evolutionary stages. Molecular transitions display a behavior that relates to changes of the environment during an evolving nuclear starburst. Molecular properties may be used to designate the evolutionary stage of the nuclear starburst. The HCN(1-0)/CO(1-0) and HCO+(1-0)/HCN(1-0) ratios serve as indicators of the time evolution of the outburst.Comment: To be published in Astronomy and Astrophysics - 11 pages, 9 figures, 1 tabl

Molecular properties of (U)LIRGs: CO, HCN, HNC and HCO+

The observed molecular properties of a sample of FIR-luminous and OH megamaser (OH-MM) galaxies have been investigated. The ratio of high and low-density tracer lines is found to be determined by the progression of the star formation in the system. The HCO+/HCN and HCO+/HNC line ratios are good proxies for the density of the gas, and PDR and XDR sources can be distinguished using the HNC/HCN line ratio. The properties of the OH-MM sources in the sample can be explained by PDR chemistry in gas with densities higher than 10^5.5 cm^-3, confirming the classical OH-MM model of IR pumped amplification with (variable) low gains.Comment: 5 pages, 2 figures, to appear in: IAU Symposium 242 Astrophysical Masers and their Environment

FUV and X-ray irradiated protoplanetary disks: a grid of models I. The disk structure

Context. Planets are thought to eventually form from the mostly gaseous (~99% of the mass) disks around young stars. The density structure and chemical composition of protoplanetary disks are affected by the incident radiation field at optical, FUV, and X-ray wavelengths, as well as by the dust properties. Aims. The effect of FUV and X-rays on the disk structure and the gas chemical composition are investigated. This work forms the basis of a second paper, which discusses the impact on diagnostic lines of, e.g., C+, O, H2O, and Ne+ observed with facilities such as Spitzer and Herschel. Methods. A grid of 240 models is computed in which the X-ray and FUV luminosity, minimum grain size, dust size distribution, and surface density distribution are varied in a systematic way. The hydrostatic structure and the thermo-chemical structure are calculated using ProDiMo. Results. The abundance structure of neutral oxygen is stable to changes in the X-ray and FUV luminosity, and the emission lines will thus be useful tracers of the disk mass and temperature. The C+ abundance distribution is sensitive to both X-rays and FUV. The radial column density profile shows two peaks, one at the inner rim and a second one at a radius r=5-10 AU. Ne+ and other heavy elements have a very strong response to X-rays, and the column density in the inner disk increases by two orders of magnitude from the lowest (LX = 1e29 erg/s) to the highest considered X-ray flux (LX = 1e32 erg/s). FUV confines the Ne+ ionized region to areas closer to the star at low X-ray luminosities (LX = 1e29 erg/s). H2O abundances are enhanced by X-rays due to higher temperatures in the inner disk and higher ionization fractions in the outer disk. The line fluxes and profiles are affected by the effects on these species, thus providing diagnostic value in the study of FUV and X-ray irradiated disks around T Tauri stars. (abridged)Comment: 47 pages, accepted by Astronomy and Astrophysics, a high resolution version of the paper is located at http://www.astro.rug.nl/~meijerink/disk_paperI_xrays.pd

Far-Infrared and Sub-Millimeter Observations and Physical Models of the Reflection Nebula Ced 201

ISO [C II] 158 micron, [O I] 63 micron, and H_2 9 and 17 micron observations are presented of the reflection nebula Ced 201, which is a photon-dominated region illuminated by a B9.5 star with a color temperature of 10,000 K (a cool PDR). In combination with ground based [C I] 609 micron, CO, 13CO, CS and HCO+ data, the carbon budget and physical structure of the reflection nebula are constrained. The obtained data set is the first one to contain all important cooling lines of a cool PDR, and allows a comparison to be made with classical PDRs. To this effect one- and three-dimensional PDR models are presented which incorporate the physical characteristics of the source, and are aimed at understanding the dominant heating processes of the cloud. The contribution of very small grains to the photo-electric heating rate is estimated from these models and used to constrain the total abundance of PAHs and small grains. Observations of the pure rotational H_2 lines with ISO, in particular the S(3) line, indicate the presence of a small amount of very warm, approximately 330 K, molecular gas. This gas cannot be accommodated by the presented models.Comment: 32 pages, 7 figures, in LaTeX. To be published in Ap

The irradiated ISM of ULIRGs

The nuclei of ULIRGs harbor massive young stars, an accreting central black hole, or both. Results are presented for molecular gas that is exposed to X-rays (1-100 keV, XDRs) and far-ultraviolet radiation (6-13.6 eV, PDRs). Attention is paid to species like HCO+, HCN, HNC, OH, H2O and CO. Line ratios of HCN/HCO+ and HNC/HCN discriminate between PDRs and XDRs. Very high J (>10) CO lines, observable with HIFI/Herschel, discriminate very well between XDRs and PDRs. In XDRs, it is easy to produce large abundances of warm (T>100 K) H2O and OH. In PDRs, only OH is produced similarly well.Comment: 5 pages, 6 figures, to appear in: IAU Symposium 242 Astrophysical Masers and their Environment

Search for Interstellar Water in the Translucent Molecular Cloud toward HD 154368

We report an upper limit of 9 x 10^{12} cm-2 on the column density of water in the translucent cloud along the line of sight toward HD 154368. This result is based upon a search for the C-X band of water near 1240 \AA carried out using the Goddard High Resolution Spectrograph of the Hubble Space Telescope. Our observational limit on the water abundance together with detailed chemical models of translucent clouds and previous measurements of OH along the line of sight constrain the branching ratio in the dissociative recombination of H_3O+ to form water. We find at the $3\sigma$ level that no more than 30% of dissociative recombinations of H_3O+ can lead to H_2O. The observed spectrum also yielded high-resolution observations of the Mg II doublet at 1239.9 \AA and 1240.4 \AA, allowing the velocity structure of the dominant ionization state of magnesium to be studied along the line of sight. The Mg II spectrum is consistent with GHRS observations at lower spectral resolution that were obtained previously but allow an additional velocity component to be identified.Comment: Accepted by ApJ, uses aasp

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.Comment: accepted by Astrophysical Journal Letter

Search for grain growth towards the center of L1544

In dense and cold molecular clouds dust grains are surrounded by thick icy mantles. It is however not clear if dust growth and coagulation take place before the switch-on of a protostar. This is an important issue, as the presence of large grains may affect the chemical structure of dense cloud cores, including the dynamically important ionization fraction, and the future evolution of solids in protoplanetary disks. To study this further, we focus on L1544, one of the most centrally concentrated pre-stellar cores on the verge of star formation, and with a well-known physical structure. We observed L1544 at 1.2 and 2 mm using NIKA, a new receiver at the IRAM 30 m telescope, and we used data from the Herschel Space Observatory archive. We find no evidence of grain growth towards the center of L1544 at the available angular resolution. Therefore, we conclude that single dish observations do not allow us to investigate grain growth toward the pre-stellar core L1544 and high sensitivity interferometer observations are needed. We predict that dust grains can grow to 200 $\mu$m in size toward the central ~300 au of L1544. This will imply a dust opacity change by a factor of ~2.5 at 1.2 mm, which can be detected using the Atacama Large Millimeter and submillimeter Array (ALMA) at different wavelengths and with an angular resolution of 2".Comment: 12 pages, 14 figures. Accepted for publication in A&