No Arabic abstract
The ionization fraction plays a key role in the chemistry and dynamics of molecular clouds. We study the H13CO+, DCO+ and HOC+ line emission towards the Horsehead, from the shielded core to the UV irradiated cloud edge, i.e., the Photodissociation Region (PDR), as a template to investigate the ionization fraction gradient in molecular clouds. We analyze a PdBI map of the H13CO+ J=1-0 line, complemented with IRAM-30m H13CO+ and DCO+ higher-J line maps and new HOC+ and CO+ observations. We compare self-consistently the observed spatial distribution and line intensities with detailed depth-dependent predictions of a PDR model coupled with a nonlocal radiative transfer calculation. The chemical network includes deuterated species, 13C fractionation reactions and HCO+/HOC+ isomerization reactions. The role of neutral and charged PAHs in the cloud chemistry and ionization balance is investigated. The detection of HOC+ reactive ion towards the Horsehead PDR proves the high ionization fraction of the outer UV irradiated regions, where we derive a low [HCO+]/[HOC+]~75-200 abundance ratio. In the absence of PAHs, we reproduce the observations with gas-phase metal abundances, [Fe+Mg+...], lower than 4x10(-9) (with respect to H) and a cosmic-rays ionization rate of zeta=(5+/-3)x10(-17) s(-1). The inclusion of PAHs modifies the ionization fraction gradient and increases the required metal abundance. The ionization fraction in the Horsehead edge follows a steep gradient, with a scale length of ~0.05 pc (or ~25), from [e-]~10(-4) (or n_e ~ 1-5 cm(-3)) in the PDR to a few times ~10(-9) in the core. PAH^- anions play a role in the charge balance of the cold and neutral gas if substantial amounts of free PAHs are present ([PAH] >10(-8)).
Micro-physical processes on interstellar dust surfaces are tightly connected to dust properties (i.e. dust composition, size and shape) and play a key role in numerous phenomena in the interstellar medium (ISM). The large disparity in physical conditions (i.e. density, gas temperature) in the ISM triggers an evolution of dust properties. The analysis of how dust evolves with the physical conditions is a stepping-stone towards a more thorough understanding of interstellar dust. The aim of this paper is to highlight dust evolution in the Horsehead Nebula PDR region. We use Spitzer/IRAC (3.6, 4.5, 5.8 and 8 {mu}m), Spitzer/MIPS (24 {mu}m) together with Herschel/PACS (70 and 160 {mu}m) and Herschel/SPIRE (250, 350 and 500 {mu}m) to map the spatial distribution of dust in the Horsehead over the entire emission spectral range. We model dust emission and scattering using the THEMIS interstellar dust model together with the 3D radiative transfer code SOC. We find that the nano-grains dust-to-gas ratio in the irradiated outer part of the Horsehead is 6 to 10 times lower than in the diffuse ISM. Their minimum size is 2 to 2.25 times larger than in the diffuse ISM and the power-law exponent of their size distribution, 1.1 to 1.4 times lower than in the diffuse ISM. Regarding the denser part of the Horsehead, it is necessary to use evolved grains (i.e. aggregates, with or without an ice mantle). It is not possible to explain the observations using grains from the diffuse medium. We therefore propose the following scenario to explain our results. In the outer part of the Horsehead, all the nano-grains have not yet had time to re-form completely through photo-fragmentation of aggregates and the smallest of the nano-grains that are sensitive to the radiation field are photo-destroyed. In the inner part of the Horsehead, grains most likely consist of multi-compositional, mantled aggregates.
We calculate the hydrogen and helium-ionizing radiation escaping star-forming molecular clouds, as a function of the star cluster mass and compactness, using a set of high-resolution radiation-magneto-hydrodynamic simulations of star formation in self-gravitating, turbulent molecular clouds. In these simulations, presented in He, Ricotti and Geen (2019), the formation of individual massive stars are well resolved, and their UV radiation feedback and lifetime on the main sequence are modeled self-consistently. We find that the escape fraction of ionizing radiation from molecular clouds, $langle f_{rm esc}^{scriptscriptstyle rm MC}rangle$, decreases with increasing mass of the star cluster and with decreasing compactness. Molecular clouds with densities typically found in the local Universe have negligible $langle f_{rm esc}^{scriptscriptstyle rm MC}rangle$, ranging between $0.5%$ to $5%$. Ten times denser molecular clouds have $langle f_{rm esc}^{scriptscriptstyle rm MC}rangle approx 10%-20%$, while $100times$ denser clouds, which produce globular cluster progenitors, have $langle f_{rm esc}^{scriptscriptstyle rm MC}rangle approx 20%-60%$. We find that $langle f_{rm esc}^{scriptscriptstyle rm MC}rangle$ increases with decreasing gas metallicity, even when ignoring dust extinction, due to stronger radiation feedback. However, the total number of escaping ionizing photons decreases with decreasing metallicity because the star formation efficiency is reduced. We conclude that the sources of reionization at $z>6$ must have been very compact star clusters forming in molecular clouds about $100times$ denser than in todays Universe, which leads to a significant production of old globular clusters progenitors.
We carried out deep searches for CO line emission in the outer disk of M33, at R>7 kpc, and examined the dynamical conditions that can explain variations in the mass distribution of the molecular cloud throughout the disk of M33. We used the IRAM-30~m telescope to search for CO lines in the outer disk toward 12 faint mid-infrared (MIR) selected sources and in an area of the southern outer disk hosting MA1, a bright HII region. We detect narrow CO lines at the location of two MIR sources at galactocentric distances of about 8 kpc that are associated with low-mass young stellar clusters, and at four locations in the proximity of MA1. The paucity of CO lines at the location of weak MIR-selected sources probably arises because most of them are not star-forming sites in M33, but background sources. Although very uncertain, the total molecular mass of the detected clouds around MA1 is lower than expected given the stellar mass of the cluster, because dispersal of the molecular gas is taking place as the HII region expands. The mean mass of the giant molecular clouds (GMCs) in M33 decreases radially by a factor 2 from the center out to 4 kpc, then it stays constant until it drops at R>7 kpc. We suggest that GMCs become more massive toward the center because of the fast rotation of the disk, which drives mass growth by coalescence of smaller condensations as they cross the arms. The analysis of both HI and CO spectral data gives the consistent result that corotation of the two main arms in this galaxy is at a radius of 4.7+-0.3 kpc, and spiral shock waves become subsonic beyond 3.9 kpc. Perturbations are quenched beyond 6.5 kpc, where CO lines have been detected only around sporadic condensations associated with UV and MIR emission.
Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic center occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties.
We study the behavior of eight diffuse interstellar bands (DIBs) in different interstellar environments, as characterized by the fraction of hydrogen in molecular form [$f$(H$_2$)], with comparisons to the corresponding behavior of various known atomic and molecular species. The equivalent widths of the five normal DIBs ($lambdalambda$5780.5, 5797.1, 6196.0, 6283.8, and 6613.6), normalized to $E(B-V)$, show a Lambda-shaped behavior: they increase at low $f$(H$_2$), peak at $f$(H$_2$) ~ 0.3, and then decrease. The similarly normalized column densities of Ca, Ca$^+$, Ti$^+$, and CH$^+$ also decline for $f$(H$_2$) > 0.3. In contrast, the normalized column densities of Na, K, CH, CN, and CO increase monotonically with $f$(H$_2$), and the trends exhibited by the three C$_2$ DIBs ($lambdalambda$4726.8, 4963.9, and 4984.8) lie between those two general behaviors. These trends with $f$(H$_2$) are accompanied by cosmic scatter, the dispersion at any given $f$(H$_2$) being significantly larger than the individual errors of measurement. The Lambda-shaped trends suggest the balance between creation and destruction of the DIB carriers differs dramatically between diffuse atomic and diffuse molecular clouds; additional processes besides ionization and shielding are needed to explain those observed trends. Except for several special cases, the highest $W$(5780)/$W$(5797) ratios, characterizing the so-called sigma-zeta effect, occur only at $f$(H$_2$) < 0.2. We propose a sequence of DIBs based on trends in their pair-wise strength ratios with increasing $f$(H$_2$). In order of increasing environmental density, we find the $lambda$6283.8 and $lambda$5780.5 DIBs, the $lambda$6196.0 DIB, the $lambda$6613.6 DIB, the $lambda$5797.1 DIB, and the C$_2$ DIBs.