The Chemistry of Interstellar OH+, H2O+, and H3O+: Inferring the Cosmic Ray Ionization Rates from Observations of Molecular Ions

We model the production of OH+, H2O+, and H3O+ in interstellar clouds, using a steady state photodissociation region code that treats the freeze-out of gas species, grain surface chemistry, and desorption of ices from grains. The code includes PAHs, which have important effects on the chemistry. All three ions generally have two peaks in abundance as a function of depth into the cloud, one at A_V<~1 and one at A_V~3-8, the exact values depending on the ratio of incident ultraviolet flux to gas density. For relatively low values of the incident far ultraviolet flux on the cloud ({\chi}<~ 1000; {\chi}= 1= local interstellar value), the columns of OH+ and H2O+ scale roughly as the cosmic ray primary ionization rate {\zeta}(crp) divided by the hydrogen nucleus density n. The H3O+ column is dominated by the second peak, and we show that if PAHs are present, N(H3O+) ~ 4x10^{13} cm^{-2} independent of {\zeta}(crp) or n. If there are no PAHs or very small grains at the second peak, N(H3O+) can attain such columns only if low ionization potential metals are heavily depleted. We also model diffuse and translucent clouds in the interstellar medium, and show how observations of N(OH+)/N(H) and N(OH+)/N(H2O+) can be used to estimate {\zeta}(crp)/n, {\chi}/n and A_V in them. We compare our models to Herschel observations of these two ions, and estimate {\zeta}(crp) ~ 4-6 x 10^-16 (n/100 cm^-3) s^-1 and \chi/n = 0.03 cm^3 for diffuse foreground clouds towards W49N

Spatially Resolved PAH Emission Features in Nearby, Low Metallicity, Star-Forming Galaxies

Low-resolution, mid-infrared Spitzer/IRS spectral maps are presented for three nearby, low-metallicity dwarf galaxies (NGC 55, NGC 3109 and IC 5152) for the purpose of examining the spatial distribution and variation of polycyclic aromatic hydrocarbon (PAH) emission. The sample straddles a metallicity of 12+log(O/H)~8.0, a transition point below which PAH intensity empirically drops and the character of the interstellar medium changes. We derive quantitative radiances of PAH features and atomic lines on both global and spatially-resolved scales. The Spitzer spectra, combined with extensive ancillary data from the UV through the mid-infrared, allow us to examine changes in the physical environments and in PAH feature radiances down to a physical scale of 50 pc. We discuss correlations between various PAH emission feature and atomic line radiances. The (6.2 micron)/(11.3 micron), (7.7 micron)/(11.3 micron), (8.6 micron)/(11.3 micron), (7.7 micron)/(6.2 micron), and (8.6 micron)/(6.2 micron) PAH radiance ratios are found to be independent of position across all three galaxies, although the ratios do vary from galaxy to galaxy. As seen in other galaxies, we find no variation in the grain size distribution as a function of local radiation field strength. Absolute PAH feature intensities as measured by a ratio of PAH/(24 micron) radiances are seen to vary both positionally within a given galaxy, and from one galaxy to another when integrated over the full observed extent of each system. We examine direct comparisons of CC mode PAH ratios (7.7 micron)/(6.2 micron) and (8.6 micron)/(6.2 micron) to the mixed (CC/CH) mode PAH ratio (7.7 micron)/(11.3 micron). We find little variation in either mode, and no difference in trends between modes. While the local conditions change markedly over the observed regions of these galaxies, the properties of PAH emission show a remarkable degree of uniformity.Comment: Astrophysical Journal, in pres

A general model for the identification of specific PAHs in the far-IR

Context. In the framework of the interstellar PAH hypothesis, far-IR skeletal bands are expected to be a fingerprint of single species in this class. Aims. A detailed model of the photophysics of interstellar PAHs is required for such single-molecule identification of their far-IR features in the presently available Infrared Space Observatory data and in those of the forthcoming Herschel Space Observatory mission. Methods. We modelled the detailed photophysics of a vast sample of species in different radiation fields, using a compendium of Monte-Carlo techniques and quantum-chemical calculations. This enabled us to validate the use of purely theoretical data and assess the expected accuracy and reliability of the resulting synthetic far-IR emission spectra. Results. We produce positions and intensities of the expected far-IR features which ought to be emitted by each species in the sample in the considered radiation fields. A composite emission spectrum for our sample is computed for one of the most favourable sources for detection, namely the Red Rectangle nebula. The resulting spectrum is compared with the estimated dust emission in the same source, to assess the dependence of detectability on key molecular parameters. Conclusions. Identifying specific PAHs from their far-IR features is going to be a difficult feat in general, still it may well be possible under favourable conditions.Comment: 14 pages, 9 figures + 18 pages of online appendix. Accepted for publication in A&A (09/06/2006

Properties of dust in the high-latitude translucent cloud L1780 I: Spatially distinct dust populations and increased dust emissivity from ISO observations

We have analyzed the properties of dust in the high galactic latitude translucent cloud Lynds 1780 using ISOPHOT maps at 100 and 200 micrometers and raster scans at 60, 80, 100, 120, 150 and 200 micrometers. In far-infrared (FIR) emission, the cloud has a single core that coincides with the maxima of visual extinction and 200um optical depth. At the resolution of 3.0 arcmin, the maximum visual extinction is 4.0 mag. At the cloud core, the minimum temperature and the maximum 200um optical depth are 14.9+/-0.4 K and 2.0+/-0.2x10^{-3}, respectively, at the resolution of 1.5 arcmin. The cloud mass is estimated to be 18M_{SUN}. The FIR observations, combined with IRAS observations, suggest the presence of different, spatially distinct dust grain populations in the cloud: the FIR core region is the realm of the "classical" large grains, whereas the very small grains and the PAHs have separate maxima on the Eastern side of the cold core, towards the "tail" of this cometary-shaped cloud. The color ratios indicate an overabundance of PAHs and VSGs in L1780. Our FIR observations combined with the optical extinction data indicate an increase of the emissivity of the big grain dust component in the cold core, suggesting grain coagulation or some other change in the properties of the large grains. Based on our observations, we also address the question, to what extent the 80um emission and even the 100um and the 120um emission contain a contribution from the small-grain component.Comment: 12 pages, 9 figures, minor changes, one table adde

Gas cooling within the diffuse ISM of late-type galaxies

We combine observations of spiral galaxies in the [CII] line at 158 micron, made with the Long Wavelength Spectrometer aboard ISO, with previous data from the Kuiper Airborne Observatory to study the origin of this line, which is the main coolant of the interstellar medium at relatively low temperatures. We also use HI and CO(1-0) observations of these galaxies and estimate the respective line fluxes in the same beam as the [CII] observations. We confirm the existence of a linear relation between the [CII] line intensity and the CO(1-0) line intensity, that we extend to intrinsically fainter galaxies. The dispersion around this relation is significant and due to variations in the far-UV flux, thus in the star formation rate. We find that for the least active galaxies of our sample, in terms of star formation, the rate of [CII] line emission per interstellar hydrogen atom is similar to that in the Solar neighbourhood. For those galaxies, most of the [CII] line emission comes probably from the diffuse cold atomic medium. In more active galaxies, considered globally, the average [CII] line emission is dominated by dense photodissociation regions and to some extent by the warm ionized diffuse medium. This is true in the central region of many spiral galaxies, and probably even in the interarm regions of the most actively star-forming ones.Comment: 8 pages, 6 Postscript figures, 2 Postscript tables, accepted to A&

Influence of UV radiation from a massive YSO on the chemistry of its envelope

We have studied the influence of far ultraviolet (UV) radiation from a massive young stellar object (YSO) on the chemistry of its own envelope by extending the models of Doty et al. (2002) to include a central source of UV radiation. The models are applied to the massive star-forming region AFGL 2591 for different inner UV field strengths. Depth-dependent abundance profiles for several molecules are presented and discussed. We predict enhanced column densities for more than 30 species, especially radicals and ions. Comparison between observations and models is improved with a moderate UV field incident on the inner envelope, corresponding to an enhancement factor G0~10-100 at 200 AU from the star with an optical depth tau~15-17. Subtle differences are found compared with traditional models of Photon Dominated Regions (PDRs) because of the higher temperatures and higher gas-phase H2O abundance caused by evaporation of ices in the inner region. In particular, the CN/HCN ratio is not a sensitive tracer of the inner UV field, in contrast with the situation for normal PDRs: for low UV fields, the extra CN reacts with H2 in the inner dense and warm region and produces more HCN. It is found that the CH+ abundance is strongly enhanced and grows steadily with increasing UV field. High-J lines of molecules like CN and HCN are most sensitive to the inner dense region where UV radiation plays a role. Thus, even though the total column density affected by UV photons is small, comparison of high-J and low-J lines can selectively trace and distinguish the inner UV field from the outer one. In addition, future Herschel-HIFI observations of hydrides can sensitively probe the inner UV field.Comment: Accepted for publication in A&A. 13 pages, 10 figure

ISM Properties in Low-Metallicity Environments II. The Dust Spectral Energy Distribution of NGC 1569

We present new 450 and 850 microns SCUBA data of the dwarf galaxy NGC 1569. We construct the mid-infrared to millimeter SED of NGC 1569, using ISOCAM, ISOPHOT, IRAS, KAO, SCUBA and MAMBO data, and model the SED in order to explore the nature of the dust in low metallicity environments. The detailed modeling is performed in a self-consistent way, synthesizing the global ISRF of the galaxy using an evolutionary synthesis model with further constraints provided by the observed MIR ionic lines and a photoionisation model. Our results show that the dust properties are different in this low metallicity galaxy compared to other more metal rich galaxies. The results indicate a paucity of PAHs probably due to the destructive effects of the ISRF penetrating a clumpy environment and a size-segregation of grains where the emission is dominated by small grains of size ~3 nm, consistent with the idea of shocks having a dramatic effect on the dust properties in NGC 1569. A significant millimetre excess is present in the dust SED which can be explained by the presence of ubiquitous very cold dust (T = 5-7 K). This dust component accounts for 40 to 70 % of the total dust mass in the galaxy (1.6 - 3.4 10^5 Msol) and could be distributed in small clumps (size a few pc) throughout the galaxy. We find a gas-to-dust mass ratio of 740 - 1600, larger than that of the Galaxy and a dust-to-metals ratio of 1/4 to 1/7. We generate an extinction curve for NGC 1569, consistent with the modeled dust size distribution. This extinction curve has relatively steep FUV rise and smaller 2175 Angstroms bump, resembling the observed extinction curve of some regions in the Large Magellanic Cloud.Comment: 20 pages, 20 figures, accepted by A&

Modeling of Protostellar Clouds and their Observational Properties

A physical model and two-dimensional numerical method for computing the evolution and spectra of protostellar clouds are described. The physical model is based on a system of magneto-gasdynamical equations, including ohmic and ambipolar diffusion, and a scheme for calculating the thermal and ionization structure of a cloud. The dust and gas temperatures are determined during the calculations of the thermal structure of the cloud. The results of computing the dynamical and thermal structure of the cloud are used to model the radiative transfer in continuum and in molecular lines. We presented the results for clouds in hydrostatic and thermal equilibrium. The evolution of a rotating magnetic protostellar cloud starting from a quasi-static state is also considered. Spectral maps for optically thick lines of linear molecules are analyzed. We have shown that the influence of the magnetic field and rotation can lead to a redistribution of angular momentum in the cloud and the formation of a characteristic rotational velocity structure. As a result, the distribution of the velocity centroid of the molecular lines can acquire an hourglass shape. We plan to use the developed program package together with a model for the chemical evolution to interpret and model observed starless and protostellar cores.Comment: Accepted to Astronomy Report

Water destruction by X-rays in young stellar objects

We study the H2O chemistry in star-forming environments under the influence of a central X-ray source and a central far ultraviolet (FUV) radiation field. The gas-phase water chemistry is modeled as a function of time, hydrogen density and X-ray flux. To cover a wide range of physical environments, densities between n_H = 10^4-10^9 cm^-3 and temperatures between T = 10-1000 K are studied. Three different regimes are found: For T < 100 K, the water abundance is of order 10^-7-10^-6 and can be somewhat enhanced or reduced due to X-rays, depending on time and density. For 100 K < T < 250 K, H2O is reduced from initial x(H2O) ~ 10^-4 following ice evaporation to x(H2O) ~ 10^-6 for F_X > 10^-3 ergs s-1 cm^-2 (t = 10^4 yrs) and for F_X > 10^-4 ergs s^-1 cm^-2 (t = 10^5 yrs). At higher temperatures (T > 250 K) and hydrogen densities, water can persist with x(H2O) ~ 10^-4 even for high X-ray fluxes. The X-ray and FUV models are applied to envelopes around low-mass Class 0 and I young stellar objects (YSOs). Water is destroyed in both Class 0 and I envelopes on relatively short timescales (t ~ 5000 yrs) for realistic X-ray fluxes, although the effect is less prominent in Class 0 envelopes due to the higher X-ray absorbing densities there. FUV photons from the central source are not effective in destroying water. The average water abundance in Class I sources for L_X > 10^27 ergs s^-1 is predicted to be x(H2O) < 10^-6.Comment: 12 pages, 14 figures, Accepted for publication in A&

Kinetics of the urea–urease clock reaction with urease immobilized in hydrogel beads

Feedback driven by enzyme catalyzed reactions occurs widely in biology and has been well characterized in single celled organisms such as yeast. There are still few examples of robust enzyme oscillators in vitro that might be used to study nonlinear dynamical behavior. One of the simplest is the urea–urease reaction that displays autocatalysis driven by the increase in pH accompanying the production of ammonia. A clock reaction was obtained from low to high pH in batch reactor and bistability and oscillations were reported in a continuous flow rector. However, the oscillations were found to be irreproducible and one contributing factor may be the lack of stability of the enzyme in solution at room temperature. Here, we investigated the effect of immobilizing urease in thiol-poly(ethylene glycol) acrylate (PEGDA) hydrogel beads, prepared using emulsion polymerization, on the urea–urease reaction. The resultant mm-sized beads were found to reproduce the pH clock and, under the conditions employed here, the stability of the enzyme was increased from hours to days