Cosmic Ray Ionization Research Articles

X-ray radiative transfer in protoplanetary disks

Context. The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. Aims. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the chemical structure of disks. Further, we want to investigate the importance of the dust for X-ray radiative transfer in disks. Methods. We present a new X-ray radiative transfer module for the radiation thermo-chemical disk code PRODIMO (PROtoplanetary DIsk MOdel), which includes X-ray scattering and absorption by both the gas and dust component. The X-ray dust opacities can be calculated for various dust compositions and dust-size distributions. For the X-ray radiative transfer we consider irradiation by the star and by X-ray background fields. To study the impact of X-rays on the chemical structure of disks we use the well established disk ionization tracers N2H+ and HCO+. Results. For evolved dust populations (e.g. grain growth), X-ray opacities are mostly dominated by the gas; only for photon energies E ≳ 5−10 keV do dust opacities become relevant. Consequently the local disk X-ray radiation field is only affected in dense regions close to the disk midplane. X-ray background fields can dominate the local X-ray disk ionization rate for disk radii r ≳ 20 au. However, the N2H+ and HCO+ column densities are only significantly affected in cases of low cosmic-ray ionization rates (≲10-19 s-1), or if the background flux is at least a factor of ten higher than the flux level of ≈10-5 erg cm-2 s-1 expected for clusters typical for the solar vicinity. Conclusions. Observable signatures of X-ray background fields in low-mass star-formation regions, like Taurus, are only expected for cluster members experiencing a strong X-ray background field (e.g. due to their location within the cluster). For the majority of the cluster members, the X-ray background field has relatively little impact on the disk chemical structure.

CO line ratios in molecular clouds: the impact of environment

Line emission is strongly dependent on the local environmental conditions in which the emitting tracers reside. In this work, we focus on modelling the CO emission from simulated giant molecular clouds (GMCs), and study the variations in the resulting line ratios arising from the emission from the $J=1-0$, $J=2-1$ and $J=3-2$ transitions. We perform a set of smoothed particle hydrodynamics (SPH) simulations with time-dependent chemistry, in which environmental conditions -- including total cloud mass, density, size, velocity dispersion, metallicity, interstellar radiation field (ISRF) and the cosmic ray ionisation rate (CRIR) -- were systematically varied. The simulations were then post-processed using radiative transfer to produce synthetic emission maps in the 3 transitions quoted above. We find that the cloud-averaged values of the line ratios can vary by up to $\pm 0.3$ dex, triggered by changes in the environmental conditions. Changes in the ISRF and/or in the CRIR have the largest impact on line ratios since they directly affect the abundance, temperature and distribution of CO-rich gas within the clouds. We show that the standard methods used to convert CO emission to H$_2$ column density can underestimate the total H$_2$ molecular gas in GMCs by factors of 2 or 3, depending on the environmental conditions in the clouds.

The onset of energetic particle irradiation in Class 0 protostars

The early stages of low-mass star formation are likely to be subject to intense ionization by protostellar energetic MeV particles. As a result, the surrounding gas is enriched in molecular ions, such as HCO$^{+}$ and N$_{2}$H$^{+}$. Nonetheless, this phenomenon remains poorly understood for Class 0 objects. Recently, based on Herschel observations taken as part of the key program Chemical HErschel Surveys of Star forming regions (CHESS), a very low HCO$^{+}$/N$_{2}$H$^{+}$ abundance ratio of about 3-4, has been reported toward the protocluster OMC-2 FIR4. This finding suggests a cosmic-ray ionization rate in excess of 10$^{-14}$ s$^{-1}$, much higher than the canonical value of $\zeta$ = 3$\times$10$^{-17}$ s$^{-1}$ (value expected in quiescent dense clouds). To assess the specificity of OMC-2 FIR4, we have extended this study to a sample of sources in low- and intermediate mass. More specifically, we seek to measure the HCO$^{+}$/N$_2$H$^{+}$ abundance ratio from high energy lines (J $\ge$ 6) toward this source sample in order to infer the flux of energetic particles in the warm and dense gas surrounding the protostars. We use observations performed with the Heterodyne Instrument for the FarInfrared spectrometer on board the Herschel Space Observatory toward a sample of 9 protostars. We report HCO$^{+}$/N$_2$H$^{+}$ abundance ratios in the range of 5 up to 73 toward our source sample. The large error bars do not allow us to conclude whether OMC-2~FIR4 is a peculiar source. Nonetheless, an important result is that the measured HCO$^{+}$/N$_2$H$^{+}$ ratio does not vary with the source luminosity. At the present time, OMC-2 FIR4 remains the only source where a high flux of energetic particles is clearly evident. More sensitive and higher angular resolution observations are required to further investigate this process.

Molecular transitions as probes of the physical conditions of extragalactic environments

Ab initio grids of time dependent chemical models, varying in gas density, temperature, cosmic ray ionization rate, and radiation field, are used as input to RADEX calculations. Tables of abundances, column densities, theoretical line intensities, and line ratios for some of the most used dense gas tracers are provided. The degree of correlation as well as degeneracy inherent in molecular ratios is discussed. Comparisons of the theoretical intensities with example observations are also provided. We find that, within the parameters space explored, chemical abundances can be constrained by a well defined set of gas density-gas temperature-cosmic ray ionization rate for the species we investigate here. However, line intensities, as well as, more importantly, line ratios, from different chemical models can be very similar leading to a clear degeneracy. We also find that the gas subjected to a galactic cosmic ray ionization rate will not necessarily have reached steady state by 1 Myr. The species most affected by time dependency effects are HCN and CS, both high density tracers. We use our method to fit an example set of data from two galaxies. We find that (i) molecular line ratios can be easily matched even with erroneous individual line intensities; (ii) no set of species can be matched by a one-component ISM; (iii) a species may be a good tracer of an energetic process but only under specific density and temperature conditions. We show that by taking into consideration the chemistry behind each species and the individual line intensities, many degeneracies that arise by just using molecular line ratios can be avoided. Finally we show that using a species or a ratio as a tracer of an individual energetic process (e.g. cosmic rays, UV) ought to be done with caution.

Heart in space: effect of the extraterrestrial environment on the cardiovascular system.

National space agencies and private corporations aim at an extended presence of humans in space in the medium to long term. Together with currently suboptimal technology, microgravity and cosmic rays raise health concerns about deep-space exploration missions. Both of these physical factors affect the cardiovascular system, whose gravity-dependence is pronounced. Heart and vascular function are, therefore, susceptible to substantial changes in weightlessness. The altered cardiovascular function in space causes physiological problems in the postflight period. A compromised cardiovascular system can be excessively vulnerable to space radiation, synergistically resulting in increased damage. The space radiation dose is significantly lower than in patients undergoing radiotherapy, in whom cardiac damage is well-documented following cancer therapy in the thoracic region. Nevertheless, epidemiological findings suggest an increased risk of late cardiovascular disease even with low doses of radiation. Moreover, the peculiar biological effectiveness of heavy ions in cosmic rays might increase this risk substantially. However, whether radiation-induced cardiovascular effects have a threshold at low doses is still unclear. The main countermeasures to mitigate the effect of the space environment on cardiac function are physical exercise, antioxidants, nutraceuticals, and radiation shielding.

IceCube Constraints on the Fermi Bubbles

Abstract We analyze the IceCube four-year neutrino data in search of a signal from the Fermi bubbles. No signal is found from the bubbles or from their dense shell, even when taking into account the softer background. This imposes a conservative ξ i < 8% upper limit on the cosmic-ray ion (CRI) acceleration efficiency, and an η ≡ ξ e /ξ i ≳ 0.006 lower limit on the electron-to-ion ratio of acceleration efficiencies (at the 2σ confidence level). For typical ξ i , a signal should surface once the number of IceCube neutrinos increases by ∼an order of magnitude, unless there is a <PeV cutoff on the CRI spectrum.

Modelling the ArH+ emission from the Crab nebula

Abstract We have performed combined photoionization and photodissociation region (PDR) modelling of a Crab nebula filament subjected to the synchrotron radiation from the central pulsar wind nebula, and to a high flux of charged particles; a greatly enhanced cosmic-ray ionization rate over the standard interstellar value, ζ0, is required to account for the lack of detected [C i] emission in published Herschel SPIRE FTS observations of the Crab nebula. The observed line surface brightness ratios of the OH+ and ArH+ transitions seen in the SPIRE FTS frequency range can only be explained with both a high cosmic-ray ionization rate and a reduced ArH+ dissociative recombination rate compared to that used by previous authors, although consistent with experimental upper limits. We find that the ArH+/OH+ line strengths and the observed H2 vibration–rotation emission can be reproduced by model filaments with nH = 2 × 104 cm−3, ζ = 107ζ0 and visual extinctions within the range found for dusty globules in the Crab nebula, although far-infrared emission from [O i] and [C ii] is higher than the observational constraints. Models with nH = 1900 cm−3 underpredict the H2 surface brightness, but agree with the ArH+ and OH+ surface brightnesses and predict [O i] and [C ii] line ratios consistent with observations. These models predict HeH+ rotational emission above detection thresholds, but consideration of the formation time-scale suggests that the abundance of this molecule in the Crab nebula should be lower than the equilibrium values obtained in our analysis.

Seeds of Life in Space (SOLIS)

The interstellar delivery of carbon atoms locked into molecules might be one of the key ingredients for the emergence of life. Cyanopolyynes are carbon chains delimited at their two extremities by an atom of hydrogen and a cyano group, so that they might be excellent reservoirs of carbon. The simplest member, HC3N, is ubiquitous in the galactic interstellar medium and found also in external galaxies. Thus, understanding the growth of cyanopolyynes in regions forming stars similar to our Sun, and what affects it, is particularly relevant. In the framework of the IRAM/NOEMA Large Program SOLIS (Seeds Of Life In Space), we have obtained a map of two cyanopolyynes, HC3N and HC5N, in the protocluster OMC2-FIR4. Because our Sun is thought to be born in a rich cluster, OMC2-FIR4 is one of the closest and best known representatives of the environment in which the Sun may have been born. We find a HC3N/HC5N abundance ratio across the source in the range ~ 1 - 30, with the smallest values (< 10) in FIR5 and in the Eastern region of FIR4. The ratios < 10 can be reproduced by chemical models only if: (1) the cosmic-ray ionisation rate $\zeta$ is ~ $4 \times 10^{-14}$ s$^{-1}$; (2) the gaseous elemental ratio C/O is close to unity; (3) O and C are largely depleted. The large $\zeta$ is comparable to that measured in FIR4 by previous works and was interpreted as due to a flux of energetic (> 10 MeV) particles from embedded sources. We suggest that these sources could lie East of FIR4 and FIR5. A temperature gradient across FIR4, with T decreasing by about 10 K, could also explain the observed change in the HC3N/HC5N line ratio, without the need of a cosmic ray ionisation rate gradient. However, even in this case, a high constant cosmic-ray ionisation rate (of the order of $10^{-14}$ s$^{-1}$) is necessary to reproduce the observations.

The Cosmic-Ray Ionization Rate in the Galactic Disk, as Determined from Observations of Molecular Ions

Abstract We have obtained estimates for the cosmic-ray ionization rate (CRIR) in the Galactic disk, using a detailed model for the physics and chemistry of diffuse interstellar gas clouds to interpret previously published measurements of the abundance of four molecular ions: ArH+, OH+, , and . For diffuse atomic clouds at Galactocentric distances in the range , observations of ArH+, OH+, and imply a mean primary CRIR of per hydrogen atom, where . Within diffuse molecular clouds observed toward stars in the solar neighborhood, measurements of and H2 imply a primary CRIR of per H atom, corresponding to a total ionization rate per H2 molecule of in good accord with previous estimates. These estimates are also in good agreement with a rederivation, presented here, of the CRIR implied by recent observations of carbon and hydrogen radio recombination lines along the sight line to Cas A. Here, our best-fit estimate for the primary CRIR is per H atom. Our results show marginal evidence that the CRIR in diffuse molecular clouds decreases with cloud extinction, , with a best-fit dependence for .

Stellar energetic particle ionization in protoplanetary disks around T Tauri stars

Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a high amount of energetic particles (E$\gtrsim$10 MeV). The particle flux of the contemporary Sun cannot explain these anomalies. However, similar to T Tauri stars the young Sun was more active and probably produced enough high energy particles to explain those anomalies. We want to study the interaction of stellar energetic particles with the gas component of the disk and identify possible observational tracers of this interaction. We use a 2D radiation thermo-chemical protoplanetary disk code to model a disk representative for T Tauri stars. We use a particle energy distribution derived from solar flare observations and an enhanced stellar particle flux proposed for T Tauri stars. For this particle spectrum we calculate the stellar particle ionization rate throughout the disk with an accurate particle transport model. We study the impact of stellar particles for models with varying X-ray and cosmic-ray ionization rates. We find that stellar particle ionization has a significant impact on the abundances of the common disk ionization tracers HCO$^+$ and N$_2$H$^+$, especially in models with low cosmic-ray ionization rates. In contrast to cosmic rays and X-rays, stellar particles cannot reach the midplane of the disk. Therefore molecular ions residing in the disk surface layers are more affected by stellar particle ionization than molecular ions tracing the cold layers/midplane of the disk. Spatially resolved observations of molecular ions tracing different vertical layers of the disk allow to disentangle the contribution of stellar particle ionization from other competing ionization sources. Modeling such observations with a model like the one presented here allows to constrain the stellar particle flux in disks around T Tauri stars.

A Simple and Accurate Network for Hydrogen and Carbon Chemistry in the Interstellar Medium

Abstract Chemistry plays an important role in the interstellar medium (ISM), regulating the heating and cooling of the gas and determining abundances of molecular species that trace gas properties in observations. Although solving the time-dependent equations is necessary for accurate abundances and temperature in the dynamic ISM, a full chemical network is too computationally expensive to incorporate into numerical simulations. In this paper, we propose a new simplified chemical network for hydrogen and carbon chemistry in the atomic and molecular ISM. We compare results from our chemical network in detail with results from a full photodissociation region (PDR) code, and also with the Nelson &amp; Langer (NL99) network previously adopted in the simulation literature. We show that our chemical network gives similar results to the PDR code in the equilibrium abundances of all species over a wide range of densities, temperature, and metallicities, whereas the NL99 network shows significant disagreement. Applying our network to 1D models, we find that the CO-dominated regime delimits the coldest gas and that the corresponding temperature tracks the cosmic-ray ionization rate in molecular clouds. We provide a simple fit for the locus of CO-dominated regions as a function of gas density and column. We also compare with observations of diffuse and translucent clouds. We find that the CO, , and abundances are consistent with equilibrium predictions for densities , but the predicted equilibrium C abundance is higher than that seen in observations, signaling the potential importance of non-equilibrium/dynamical effects.

Star formation rates and efficiencies in the Galactic Centre

The inner few hundred parsecs of the Milky Way harbours gas densities, pressures, velocity dispersions, an interstellar radiation field and a cosmic ray ionisation rate orders of magnitude higher than the disc; akin to the environment found in star-forming galaxies at high-redshift. Previous studies have shown that this region is forming stars at a rate per unit mass of dense gas which is at least an order of magnitude lower than in the disc, potentially violating theoretical predictions. We show that all observational star formation rate diagnostics - both direct counting of young stellar objects and integrated light measurements - are in agreement within a factor two, hence the low star formation rate is not the result of the systematic uncertainties that affect any one method. As these methods trace the star formation over different timescales, from $0.1 - 5$ Myr, we conclude that the star formation rate has been constant to within a factor of a few within this time period. We investigate the progression of star formation within gravitationally bound clouds on $\sim$ parsec scales and find $1 - 4$ per cent of the cloud masses are converted into stars per free-fall time, consistent with a subset of the considered "volumetric" star formation models. However, discriminating between these models is obstructed by the current uncertainties on the input observables and, most importantly and urgently, by their dependence on ill-constrained free parameters. The lack of empirical constraints on these parameters therefore represents a key challenge in the further verification or falsification of current star formation theories.

Cosmic-ray Induced Destruction of CO in Star-forming Galaxies

Abstract We explore the effects of the expected higher cosmic ray (CR) ionization rates on the abundances of carbon monoxide (CO), atomic carbon (C), and ionized carbon (C+) in the H2 clouds of star-forming galaxies. The study of Bisbas et al. is expanded by (a) using realistic inhomogeneous giant molecular cloud (GMC) structures, (b) a detailed chemical analysis behind the CR-induced destruction of CO, and (c) exploring the thermal state of CR-irradiated molecular gas. CRs permeating the interstellar medium with are found to significantly reduce the [CO]/[H2] abundance ratios throughout the mass of a GMC. CO rotational line imaging will then show much clumpier structures than the actual ones. For (Galactic) this bias becomes severe, limiting the usefulness of CO lines for recovering structural and dynamical characteristics of H2-rich galaxies throughout the universe, including many of the so-called main-sequence galaxies where the bulk of cosmic star formation occurs. Both C+ and C abundances increase with rising , with C remaining the most abundant of the two throughout H2 clouds, when (Galactic). C+ starts to dominate for (Galactic). The thermal state of the gas in the inner and denser regions of GMCs is invariant with for (Galactic). For (Galactic) this is no longer the case and are reached. Finally, we identify OH as the key species whose T gas-sensitive abundance could mitigate the destruction of CO at high temperatures.

On the chemistry of the young massive protostellar core NGC 2264 CMM3

We present the first gas-grain astrochemical model of the NGC 2264 CMM3 protostellar core. The chemical evolution of the core is affected by changing its physical parameters such as the total density and the amount of gas-depletion onto grain surfaces as well as the cosmic ray ionisation rate, $\zeta$ . We estimated $\zeta_{\text{CMM3}} = 1.6 \times 10^{-17}~\mbox{s}^{-1}$ . This value is 1.3 times higher than the standard CR ionisation rate, $\zeta _{\text{ISM}} = 1.3 \times 10^{-17}~\mbox{s}^{-1}$ . Species response differently to changes into the core physical conditions, but they are more sensitive to changes in the depletion percentage and CR ionisation rate than to variations in the core density. Gas-phase models highlighted the importance of surface reactions as factories of large molecules and showed that for sulphur bearing species depletion is important to reproduce observations. Comparing the results of the reference model with the most recent millimeter observations of the NGC 2264 CMM3 core showed that our model is capable of reproducing the observed abundances of most of the species during early stages ( $\le 3\times10^{4}~\mbox{yrs}$ ) of their chemical evolution. Models with variations in the core density between $1\mbox{--}20 \times 10^{6}~\mbox{cm}^{-3}$ are also in good agreement with observations during the early time interval $1 \times 10^{4} < t\ (\mbox{yr}) < 5 \times 10^{4}$ . In addition, models with higher CR ionisation rates $(5\mbox{--}10) \times\zeta_{\text{ISM}}$ are often overestimating the fractional abundances of the species. However, models with $\zeta_{\text{CMM3}} = 5 \zeta_{\text{ISM}}$ may best fit observations at times $\sim 2 \times 10^{4}~\mbox{yrs}$ . Our results suggest that CMM3 is $(1\mbox{--}5) \times 10^{4}~\mbox{yrs}$ old. Therefore, the core is chemically young and it may host a Class 0 object as suggested by previous studies.

Robustness of N2H+as tracer of the CO snowline

[Abridged] Snowlines in protoplanetary disks play an important role in planet formation and composition. Since the CO snowline is difficult to observe directly with CO emission, its location has been inferred in several disks from spatially resolved ALMA observations of DCO+ and N2H+. N2H+ is considered to be a good tracer of the CO snowline based on astrochemical considerations predicting an anti-correlation between N2H+ and gas-phase CO. In this work, the robustness of N2H+ as a tracer of the CO snowline is investigated. A simple chemical network is used in combination with the radiative transfer code LIME to model the N2H+ distribution and corresponding emission in the disk around TW Hya. The assumed CO and N2 abundances, corresponding binding energies, cosmic ray ionization rate, and degree of large-grain settling are varied to determine the effects on the N2H+ emission and its relation to the CO snowline. For the adopted physical structure of the TW Hya disk and molecular binding energies for pure ices, the balance between freeze-out and thermal desorption predicts a CO snowline at 19 AU, corresponding to a CO midplane freeze-out temperature of 20 K. A model with a total, i.e. gas plus ice, CO abundance of 3e-6 with respect to H2 fits the position of the emission peak observed by Qi et al. 2013 for the TW Hya disk. However, the relationship between N2H+ and the CO snowline is more complicated than generally assumed: for the investigated parameters, the N2H+ column density peaks at least 5 AU outside the CO snowline. Moreover, the N2H+ emission can peak much further out, as far as ~50 AU beyond the snowline. Hence, chemical modeling, as done here, is necessary to derive a CO snowline location from N2H+ observations.

Modelling the chemistry of star-forming filaments – II. Testing filament characteristics with synthetic observations

We present synthetic continuum and $^{13}$CO and C$^{18}$O line emission observations of dense and cold filaments. The filaments are dynamically evolved using 3D-MHD simulations that include one of the largest on-the-fly chemical networks used to date, which models the detailed evolution of H$_2$ and CO. We investigate the reliability of observable properties, in particular filament mass and width, under different simulation conditions like magnetic field orientation and cosmic ray ionisation rate. We find that filament widths of $\sim$0.1 pc can be probed with both line and continuum emission observations with a high accuracy (deviations $\leq$ 20%). However, the width of more narrow filaments can be significantly overestimated by up to a factor of a few. Masses obtained via the dust emission are accurate within a few percent whereas the masses inferred from molecular line emission observations deviate from the actual mass by up to a factor of 10 and show large differences for different $J$ transitions. The inaccurate estimate of filament masses and widths of narrow filaments using molecular line observations can be attributed to (i) the non-isothermal state of the filaments, (ii) optical depth effects, and (iii) the subthermally excited state of CO, while inclination effects and opacity correction only influence the obtained masses and widths by less than 50%. Both, mass and width estimates, can be improved by using two isotopes to correct for the optical depth. Since gas and dust temperature generally differ (by up to 25 K), the filaments appear more gravitationally unstable if the (too low) dust temperature is used for the stability analysis.

A study of methanol and silicon monoxide production through episodic explosions of grain mantles in the Central Molecular Zone

Methanol (CH$_3$OH) is found to be abundant and widespread towards the Central Molecular Zone, the inner few hundred parsecs of our Galaxy. Its origin is, however, not fully understood. It was proposed that the high cosmic ray ionisation rate in this region could lead to a more efficient non-thermal desorption of this species formed on grain surfaces, but it would also mean that this species is destroyed in a relatively short timescale. In a first step, we run chemical models with a high cosmic ray ionisation rate and find that this scenario can only reproduce the lowest abundances of methanol derived in this region ($\sim$10$^{-9}$-10$^{-8}$). In a second step, we investigate another scenario based on episodic explosions of grain mantles. We find a good agreement between the predicted abundances of methanol and the observations. We find that the dominant route for the formation of methanol is through hydrogenation of CO on the grains followed by the desorption due to the grain mantle explosion. The cyclic aspect of this model can explain the widespread presence of methanol without requiring any additional mechanism. We also model silicon monoxide (SiO), another species detected in several molecular clouds of the Galactic Centre. An agreement is found with observations for a high depletion of Si (Si/H $\sim$ 10$^{-8}$) with respect to the solar abundance.

Discovery of a Perseus-like cloud in the early Universe

We present the discovery of a molecular cloud at zabs ≈ 2.5255 along the line of sight to the quasar SDSS J 000015.17+004833.3. We use a high-resolution spectrum obtained with the Ultraviolet and Visual Echelle Spectrograph together with a deep multi-wavelength medium-resolution spectrum obtained with X-shooter (both on the Very Large Telescope) to perform a detailed analysis of the absorption lines from ionic, neutral atomic and molecular species in different excitation levels, as well as the broad-band dust extinction. We find that the absorber classifies as a Damped Lyman-α system (DLA) with log N(H i) (cm-2) = 20.8 ± 0.1. The DLA has super-solar metallicity (Z ~ 2.5 Z⊙, albeit to within a factor of two to three) with a depletion pattern typical of cold gas and an overall molecular fraction f = 2N(H2)/(2N(H2) + N(H i)) ~ 50%. This is the highest f-value observed to date in a high-z intervening system. Most of the molecular hydrogen arises from a clearly identified narrow (b ~ 0.7km s-1), cold component in which carbon monoxide molecules are also found, with log N(CO)≈ 15. With the help of the spectral synthesis code Cloudy, we study the chemical and physical conditions in the cold gas. We find that the line of sight probes the gas deep after the H i-to-H2 transition in a ~4−5 pc-size cloud with volumic density nH ~ 80 cm-3 and temperature of only 50 K. Our model suggests that the presence of small dust grains (down to about 0.001 μm) and high cosmic ray ionisation rate (ζH ~ a few times 10-15 s-1) are needed to explain the observed atomic and molecular abundances. The presence of small grains is also in agreement with the observed steep extinction curve that also features a 2175 Å bump. Interestingly, the chemical and physical properties of this cloud are very similar to what is seen in diffuse molecular regions of the nearby Perseus complex, despite the former being observed when the Universe was only 2.5 Gyr old. The high excitation temperature of CO rotational levels towards J0000+0048 betrays however the higher temperature of the cosmic microwave background. Using the derived physical conditions, we correct for a small contribution (0.3 K) of collisional excitation and obtain TCMB(z = 2.53) ≈ 9.6 K, in perfect agreement with the predicted adiabatic cooling of the Universe.

CONSTRAINING THE ISM PROPERTIES OF THE CLOVERLEAF QUASAR HOST GALAXY WITH HERSCHEL SPECTROSCOPY

ABSTRACT We present Herschel observations of the far-infrared (FIR) fine-structure (FS) lines [C ii]158 μm, [O i]63 μm, [O iii]52 μm, and [Si ii]35 μm in the z = 2.56 Cloverleaf quasar, and combine them with published data in an analysis of the dense interstellar medium (ISM) in this system. Observed [C ii]158 μm, [O i]63 μm, and FIR continuum flux ratios are reproduced with photodissociation region (PDR) models characterized by moderate far-ultraviolet (FUV) radiation fields with 0.3–1 × 103 and atomic gas densities 3–5 × 103 cm−3, depending on contributions to [C ii]158 μm from ionized gas. We assess the contribution to the [C ii]158 μm flux from an active galactic nucleus (AGN) narrow line region (NLR) using ground-based measurements of the [N ii]122 μm transition, finding that the NLR can contribute at most 20%–30% of the observed [C ii]158 μm flux. The PDR density and far-UV radiation fields inferred from the atomic lines are not consistent with the CO emission, indicating that the molecular gas excitation is not solely provided via UV heating from local star formation (SF), but requires an additional heating source. X-ray heating from the AGN is explored, and we find that X-ray-dominated region (XDR) models, in combination with PDR models, can match the CO cooling without overproducing the observed FS line emission. While this XDR/PDR solution is favored given the evidence for both X-rays and SF in the Cloverleaf, we also investigate alternatives for the warm molecular gas, finding that either mechanical heating via low-velocity shocks or an enhanced cosmic-ray ionization rate may also contribute. Finally, we include upper limits on two other measurements attempted in the Herschel program: [C ii]158 μm in FSC 10214 and [O i]63 μm in APM 08279+5255.

Carbon and hydrogen radio recombination lines from the cold clouds towards Cassiopeia A

We use the Low Frequency Array to perform a systematic high spectral resolution investigation of the low-frequency 33-78 MHz spectrum along the line of sight to Cassiopeia A. We complement this with a 304-386 MHz Westerbork Synthesis Radio telescope observation. In this first paper we focus on the carbon radio recombination lines. We detect Cn$\alpha$ lines at -47 and -38 km s$^{-1}$ in absorption for quantum numbers n=438-584 and in emission for n=257-278 with high signal to noise. These lines are associated with cold clouds in the Perseus spiral arm component. Hn$\alpha$ lines are detected in emission for n=257-278. In addition, we also detect Cn$\alpha$ lines at 0 km s$^{-1}$ associated with the Orion arm. We analyze the optical depth of these transitions and their line width. Our models show that the carbon line components in the Perseus arm are best fit with an electron temperature 85 K and an electron density 0.04 cm$^{-3}$ and can be constrained to within 15\%. The electron pressure is constrained to within 20\%. We argue that much of these carbon radio recombination lines arise in the CO-dark surface layers of molecular clouds where most of the carbon is ionized but hydrogen has made the transition from atomic to molecular. The hydrogen lines are clearly associated with the carbon line emitting clouds, but the low-frequency upperlimits indicate that they likely do not trace the same gas. Combining the hydrogen and carbon results we arrive at a firm lower limit to the cosmic ray ionization rate of 2.5$\times$10$^{-18}$ s$^{-1}$, but the actual value is likely much larger.