Dust and hydrogen molecules in the extremely metal-poor dwarf galaxy SBS 0335–052

During the early stages of galaxy evolution, the metallicity is generally low and nearby metal-poor star-forming galaxies may provide templates for primordial star formation. In particular, the dust content of such objects is of great importance, because early molecular formation can take place on grains. To gain insight into primeval galaxies at high redshift, we examine the dust content of the nearby extremely low-metallicity galaxy SBS  0335–052  which hosts a very young starburst (≲10⁷ yr). In young galaxies, the dust formation rate in Type II supernovae governs the amount of dust, and by incorporating recent results on dust production in Type II supernovae we model the evolution of dust content. If the star-forming region is compact (≲100 pc), as suggested by observations of SBS  0335–052  , our models consistently explain the quantity of dust, far-infrared luminosity, and dust temperature in this low-metallicity object. We also discuss the H₂ abundance. The compactness of the region is important to H₂ formation, because the optical depth of dust for UV photons becomes large and H₂ dissociation is suppressed. We finally focus on implications for damped Ly α systems.     

Ram-pressure stripping of the interstellar medium in NGC 4485

This paper describes submm, ¹²CO (J = 2–1) observations of the interacting pair of galaxies NGC 4490 and 4485, and together with high resolution H  i and multifrequency radio continuum data we investigate the evolution of the ISM in this system. We find the following. (i) The smaller member of the pair, NGC 4485, has had the atomic, molecular and dust components of its ISM stripped via ram pressure during its recent passage through the extended H  i distribution of NGC 4490. A bow-shock is identified in the H  i ahead of the stripped gas. (ii) Within the disc of NGC 4490 we find a very low H₂-to-H  i ratio as well as a strong correlation between thermal emission and mass of H₂ suggesting that the star formation rate is limited in this case by the conversion of H  i to H₂. (iii) ¹²CO emission from an H  i and radio-continuum bridge between the two galaxies is detected.     

Molecular hydrogen in damped Lyα systems: clues to interstellar physics at high redshift

In order to interpret H₂ quasar absorption-line observations of damped Lyα systems (DLAs) and subDLAs, we model their H₂ abundance as a function of dust-to-gas ratio, including H₂ self-shielding and dust extinction against dissociating photons. Then, we constrain the physical state of the gas by using H₂ data. Using H₂ excitation data for DLAs with H₂ detections, we derive a gas density  1.5 ≲ log n (cm⁻³) ≲ 2.5  , temperature  1.5 ≲ log T (K) ≲ 3  , and an internal ultraviolet (UV) radiation field (in units of the Galactic value)  0.5 ≲ log χ≲ 1.5  . We then find that the observed relation between the molecular fraction and the dust-to-gas ratio of the sample is naturally explained by the above conditions. However, it is still possible that H₂ deficient DLAs and subDLAs with H₂ fractions less than  ∼10⁻⁶  are in a more diffuse and warmer state. The efficient photodissociation by the internal UV radiation field explains the extremely small H₂ fraction  (≲10⁻⁶)  observed for  κ≲ 1/30  (κ is the dust-to-gas ratio in units of the Galactic value); H₂ self-shielding causes a rapid increase in, and large variations of, H₂ abundance for  κ≳ 1/30  . We finally propose an independent method to estimate the star formation rates of DLAs from H₂ abundances; such rates are then critically compared with those derived from other proposed methods. The implications for the contribution of DLAs to the cosmic star formation history are briefly discussed.     

The fragmentation of pre-enriched primordial objects

Recent theoretical investigations have suggested that the formation of the very first stars, forming out of metal-free gas, was fundamentally different from the present-day case. The question then arises which effect was responsible for this transition in the star formation properties. In this paper, we study the effect of metallicity on the evolution of the gas in a collapsing dark matter mini-halo. We model such a system as an isolated 3 σ peak of mass     that collapses at     , using smoothed particle hydrodynamics. The gas has a supposed level of pre-enrichment of either     or 10⁻³ Z_⊙. We assume that H₂ has been radiatively destroyed by the presence of a soft UV background. Metals therefore provide the only viable cooling at temperatures below 10⁴ K. We find that the evolution proceeds very differently for the two cases. The gas in the lower metallicity simulation fails to undergo continued collapse and fragmentation, whereas the gas in the higher metallicity case dissipatively settles into the centre of the dark matter halo. The central gas, characterized by densities     , and a temperature,     , that closely follows that of the cosmic microwave background, is gravitationally unstable and undergoes vigorous fragmentation. We discuss the physical reason for the existence of a critical metallicity,     , and its possible dependence on redshift. Compared with the pure H/He case, the fragmentation of the     gas leads to a larger relative number of low-mass clumps.     

The ortho- to para- ratio of H2 in the starburst of NGC 253

We present measurements of several near-infrared emission lines from the nearby galaxy NGC 253. We have been able to measure four H₂ lines across the circumnuclear starburst, from which we estimate the ortho- to para- ratio of excited H₂ to be ∼2. This indicates that the bulk of the H₂ emission arises from photodissociation regions (PDRs), rather than from shocks. This is the case across the entire region of active star formation.
As the H₂ emission arises from PDRs, it is likely that the ratio of H₂ to Brγ (the bright hydrogen recombination line) is a measure of the relative geometry of O and B stars and PDRs. Towards the nucleus of NGC 253 the geometry is deduced to be tightly clustered O and B stars in a few giant H  II regions that are encompassed by PDRs. Away from the nuclear region, the geometry becomes that of PDRs bathed in a relatively diffuse ultraviolet radiation field.
The rotation curves of 1–0 S(1) and Brγ suggest that the ionized gas is tracing a kinetic system different from that of the molecular gas in NGC 253, particularly away from the nucleus.     

Clues for the origin of the fundamental metallicity relations – I. The hierarchical building up of the structure

We analyse the evolutionary history of galaxies formed in a hierarchical scenario consistent with the concordance Lambda cold dark matter (ΛCDM) model focusing on the study of the relation between their chemical and dynamical properties. Our simulations consistently describe the formation of the structure and its chemical enrichment within a cosmological context. Our results indicate that the luminosity–metallicity and the stellar mass–metallicity (LZR and MZR) relations are naturally generated in a hierarchical scenario. Both relations are found to evolve with redshift. In the case of the MZR, the estimated evolution is weaker than that deduced from observational works by approximately 0.10 dex. We also determine a characteristic stellar mass, M _c≈ 3 × 10¹⁰ M_⊙, which segregates the simulated galaxy population into two distinctive groups and which remains unchanged since z ∼ 3, with a very weak evolution of its metallicity content. The value and role played by M _c is consistent with the characteristic mass estimated from the SDSS galaxy survey by Kauffmann et al. Our findings suggest that systems with stellar masses smaller than M _c are responsible for the evolution of this relation at least from z ≈ 3. Larger systems are stellar dominated and have formed more than 50 per cent of their stars at   z ≥ 2  , showing very weak evolution since this epoch. We also found bimodal metallicity and age distributions from z ∼ 3, which reflects the existence of two different galaxy populations. Although SN feedback may affect the properties of galaxies and help to shape the MZR, it is unlikely that it will significantly modify M _c since, from   z = 3  this stellar mass is found in systems with circular velocities larger than 100 km s⁻¹.     

Collapse of low-mass clouds in the presence of a UV radiation field

The collapse of marginally Jeans-unstable primordial gas clouds in the presence of a UV radiation field is discussed. Assuming that the dynamical collapse proceeds approximately in an isothermal self-similar fashion, we investigate the thermal evolution of the collapsing central core until H₂ cooling dominates photoheating and the temperature drops to below 10⁴ K. Consequently, the mass of the cooled core is evaluated as M _cool=3.6×10⁶ M_⊙ ( I ₂₁/1)^−0.32. This scale depends only on the incident UV intensity, and provides a lower limit to the mass of collapsed objects in the UV radiation field.     

Fossil H ii regions: self-limiting star formation at high redshift

Recent results from the Wilkinson Microwave Anisotropy Probe ( WMAP ) satellite suggest that the intergalactic medium (IGM) was significantly reionized at redshifts as high as   z ∼ 17  . At this early epoch, the first ionizing sources probably appeared in the shallow potential wells of mini-haloes with virial temperatures   T _vir < 10⁴ K  . Once such an ionizing source turns off, its surrounding H ii region Compton cools and recombines. None the less, we show that the 'fossil' H ii regions left behind remain at high adiabats, prohibiting gas accretion and cooling in subsequent generations of mini-haloes. This greatly amplifies feedback effects explored in previous studies, and early star formation is self-limiting. We quantify this effect to show that star formation in mini-haloes cannot account for the bulk of the electron scattering opacity measured by WMAP , which must be due to more massive objects. We argue that gas entropy, rather than IGM metallicity, regulates the evolution of the global ionizing emissivity and impedes full reionization until lower redshifts. We discuss several important consequences of this early entropy floor for reionization. It reduces gas clumping, curtailing the required photon budget for reionization. An entropy floor also prevents H₂ formation and cooling, due to reduced gas densities: it greatly enhances feedback from ultraviolet photodissociation of H₂. An early X-ray background would also furnish an entropy floor to the entire IGM; thus, X-rays impede rather than enhance H₂ formation. Future 21-cm observations may probe the topology of fossil H ii regions.     

On the possibility of observing H2 emission from primordial molecular cloud kernels

We study the prospects for observing H₂ emission during the assembly of primordial molecular cloud kernels. The primordial molecular cloud cores, which resemble those at the present epoch, can emerge around  1+ z ∼20  according to recent numerical simulations. The kernels form inside the cores, and the first stars will appear inside the kernels. A kernel typically contracts to form one of the first generation stars with an accretion rate that is as large as ∼0.01 M_⊙ yr⁻¹. This occurs owing to the primordial abundances, which result in a kernel temperature of order 1000 K, and the collapsing kernel emits H₂ line radiation at a rate ∼10³⁵ erg s⁻¹. Predominantly   J =5-3   ( v =0)  rotational emission of H₂ is expected. At redshift  1+ z ∼20  , the expected flux is ∼0.01 μJy for a single kernel. While an individual object is not observable by any facilities available in the near future, the expected assembly of primordial star clusters on subgalactic scales can result in fluxes at the sub-mJy level. This is marginally observable with ASTRO-F and ALMA. We also examine the rotational   J =2-0   ( v =0)  and vibrational   δv =1  emission lines. The former may possibly be detectable with ALMA.     

Dust formation and evolution in a Ca–Fe–SiO–H2–O2 vapour phase condensation experiment and astronomical implications

Here, we report on a kinetically controlled vapour phase condensation experiment using a low-calcium Ca–Fe–SiO–H₂–O₂ vapour. Under these conditions of extreme disequilibrium, the condensate properties become predictable. They are amorphous solids with (predictable) deep metastable eutectic compositions. This study also shows how chemical evolution of the condensate grains will lead to chemically complex amorphous solids. The highly disordered structure of the deep metastable eutectic condensates is the very key to this predictable chemical evolution to grains with a silicate mineral composition, yet being amorphous. We compare our results with astronomical observations of dust around young stellar objects.     

A statistically selected Chandra sample of 20 galaxy clusters – II. Gas properties and cool core/non-cool core bimodality

We investigate the thermodynamic and chemical structure of the intracluster medium (ICM) across a statistical sample of 20 galaxy clusters analysed with the Chandra X-ray satellite. In particular, we focus on the scaling properties of the gas density, metallicity and entropy and the comparison between clusters with and without cool cores (CCs). We find marked differences between the two categories except for the gas metallicity, which declines strongly with radius for all clusters  ( Z ∝ r ^−0.31)  , outside  ∼0.02 r ₅₀₀  . The scaling of gas entropy is non-self-similar and we find clear evidence of bimodality in the distribution of logarithmic slopes of the entropy profiles. With only one exception, the steeper sloped entropy profiles are found in CC clusters whereas the flatter slope population are all non-CC clusters. We explore the role of thermal conduction in stabilizing the ICM and conclude that this mechanism alone is sufficient to balance cooling in non-CC clusters. However, CC clusters appear to form a distinct population in which heating from feedback is required in addition to conduction. Under the assumption that non-CC clusters are thermally stabilized by conduction alone, we find the distribution of Spitzer conduction suppression factors, f _c, to be lognormal, with a log (base 10) mean of  −1.50 ± 0.03  (i.e.   f _c= 0.032  ) and log standard deviation  0.39 ± 0.02  .     

Dust formation in primordial Type II supernovae

We have investigated the formation of dust in the ejecta of Type II supernovae (SNe), mostly of primordial composition, to answer the question of where the first solid particles are formed in the Universe. However, we have also considered non-zero progenitor metallicity values up to Z = Z_⊙ . The calculations are based on standard nucleation theory, and the scheme has been tested for the first time on the well-studied case of SN1987A, yielding results that are in agreement with the available data. We find that: (i) the first dust grains are predominantly made of silicates, amorphous carbon (AC), magnetite and corundum; and (ii) the largest grains are the AC ones, with sizes around 300 Å, whereas the other grain types have smaller radii, around 10–20 Å . The grain size distribution depends somewhat on the thermodynamics of the ejecta expansion, and variations in the results by a factor ≈2 might occur within reasonable estimates of the relevant parameters. Also, and for the same reason, the grain size distribution is essentially unaffected by metallicity changes. The predictions on the amount of dust formed are very robust: for Z =0 , we find that SNe with masses in the range (12–35) M_⊙ produce about 0.08 M_⊙≲ M _d≲0.3 M_⊙ of dust per supernova. The above range increases by roughly three times as the metallicity is increased to solar values. We discuss the implications and the cosmological consequences of the results.     

Molecular hydrogen formation on porous dust grains

Recent laboratory experiments on interstellar dust analogues have shown that H₂ formation on dust-grain surfaces is efficient in a range of grain temperatures below 20 K. These results indicate that surface processes may account for the observed H₂ abundance in cold diffuse and dense clouds. However, high abundances of H₂ have also been observed in warmer clouds, including photon-dominated regions (PDRs), where grain temperatures may reach 50 K, making the surface processes extremely inefficient. It was suggested that this apparent discrepancy can be resolved by chemisorption sites. However, recent experiments indicate that chemisorption processes may not be efficient at PDR temperatures. Here we consider the effect of grain porosity on H₂ formation, and analyse it using a rate-equation model. It is found that porosity extends the efficiency of the recombination process to higher temperatures. This is because H atoms that desorb from the internal surfaces of the pores may re-adsorb many times and thus stay longer on the surface. However, this porosity-driven extension may enable efficient H₂ formation in PDRs only if porosity also contributes to significant cooling of the grains, compared to non-porous grains.     

A possible theoretical explanation of metallicity gradients in elliptical galaxies

Models of chemical evolution of elliptical galaxies taking into account different escape velocities at different galactocentric radii are presented. As a consequence of this, the chemical evolution develops differently in different galactic regions; in particular, we find that the galactic wind, powered by supernovae (of Type II and I) starts, under suitable conditions, in the outer regions and successively develops in the central ones. The star formation is assumed to stop after the onset of the galactic wind in each region. The main result found in the present work is that this mechanism is able to reproduce metallicity gradients, namely the gradients in the Mg₂ index, in good agreement with observational data. We also find that in order to honour the constant [〈Mg/Fe〉] ratio with galactocentric distance, as inferred from metallicity indices, a variable initial mass function as a function of galactocentric distance is required. This is only a suggestion, as trends on abundances inferred purely from metallicity indices are still uncertain.     

Colour–magnitude relations and spectral line strengths in the Coma cluster

We use the C₂4668, Fe4383, H γ _A and H δ _A spectral absorption line indices, together with U - and V -band photometry of 101 galaxies in the Coma cluster, to investigate how mean age and metal abundance correlate with galaxy luminosity. In particular, we use the line index measurements to address the origin of the colour–magnitude relation (CMR). We find that the CMR in Coma is driven primarily by a luminosity–metallicity correlation. We additionally show evidence for a relation between age and luminosity, in the direction predicted by the semi-analytic hierarchical clustering models of Kauffmann & Charlot, but this is only present in the C₂4668 index models, and could be an effect of the lack of non-solar abundance ratios in the Worthey models used.
By comparing deviations from the CMR with deviations in absorption index from analogous 'index–magnitude' relations, we find that colour deviations bluewards of the mean relation are strongly correlated with the hydrogen Balmer line series absorption. We show that the properties of these blue galaxies are consistent with the presence of a young stellar population in the galaxies, rather than with a reduced metallicity.     

Star formation triggered by supernova explosions in young galaxies

We study the evolution of supernova remnants in a low-metallicity medium   Z /Z_⊙= 10⁻⁴ to 10⁻²  in the early universe, using one-dimensional hydrodynamics with non-equilibrium chemistry. Once a post-shock layer is able to cool radiatively, a dense shell forms behind the shock. If this shell becomes gravitationally unstable and fragments into pieces, next-generation stars are expected to form from these fragments. To explore the possibility of this triggered star formation, we apply a linear perturbation analysis of an expanding shell to our results and constrain the parameter range of ambient density, explosion energy and metallicity where fragmentation of the shell occurs. For the explosion energy of  10⁵¹ erg (10⁵² erg)  , the shell fragmentation occurs for ambient densities higher than  ≳10² cm⁻³ (10 cm⁻³  ), respectively. This condition depends little on the metallicity in the ranges we examined. We find that the mode of star formation triggered occurs only in massive  (≳10⁸ M_⊙)  haloes.     

Shattering and coagulation of dust grains in interstellar turbulence

We investigate shattering and coagulation of dust grains in turbulent interstellar medium (ISM). The typical velocity of dust grain as a function of grain size has been calculated for various ISM phases based on a theory of grain dynamics in compressible magnetohydrodynamic turbulence. In this paper, we develop a scheme of grain shattering and coagulation and apply it to turbulent ISM by using the grain velocities predicted by the above turbulence theory. Since large grains tend to acquire large velocity dispersions as shown by earlier studies, large grains tend to be shattered. Large shattering effects are indeed seen in warm ionized medium within a few Myr for grains with radius   a ≳ 10⁻⁶  cm. We also show that shattering in warm neutral medium can limit the largest grain size in ISM  ( a ∼ 2 × 10⁻⁵ cm)  . On the other hand, coagulation tends to modify small grains since it only occurs when the grain velocity is small enough. Coagulation significantly modifies the grain size distribution in dense clouds (DC), where a large fraction of the grains with   a < 10⁻⁶ cm  coagulate in 10 Myr. In fact, the correlation among   R_V   , the carbon bump strength and the ultraviolet slope in the observed Milky Way extinction curves can be explained by the coagulation in DC. It is possible that the grain size distribution in the Milky Way is determined by a combination of all the above effects of shattering and coagulation. Considering that shattering and coagulation in turbulence are effective if dust-to-gas ratio is typically more than ∼1/10 of the Galactic value, the regulation mechanism of grain size distribution should be different between metal-poor and metal-rich environments.     

Chemically driven desorption of CO from icy grains in dark clouds

The chemical desorption of an adsorbed CO molecule in the vicinity of H₂-forming sites on cosmic dust grains in cold dense clouds is investigated theoretically, mainly using a model based on a classical molecular dynamics computational simulation. As a model surface for icy mantles of dust grains, an amorphous water ice slab is generated at 10 K, and the first and the second H atoms are thrown on to the model surface to reproduce the recombination process of the two H atoms, H+H→H₂. Then, the time and space dependence of the local temperature increase of icy mantles caused by the release of H₂ formation energy in the vicinity of H₂-forming sites is examined. It is found that icy mantles are heated locally up to about 30 K in the surface region at R 4 Å and about 20 K at 4 R 6 Å, where R is the distance from the H₂-forming site. The critical temperature of CO desorption is estimated to be about 20–30 K under conditions in typical dense clouds, which might be seen to be comparable to the above result. However, the lifetime of local heating of icy mantles is found to be too short, compared with the time-scale of CO desorption (10¹³ s) and that for H₂ forming in the vicinity of an adsorbed CO molecule (more than 2×10¹³ s). Thus, it is found that the efficiency of chemical desorption of CO on a large dust grain is negligible. On the other hand, chemical desorption can occur on a small dust grain with size less than 20 Å.     

Shocked gas around Cepheus A: evidence for multiple outflows from H2S and SO2 observations

The Cepheus A star-forming region has been investigated through a multiline H₂S and SO₂ survey at millimetre wavelengths. Large-scale maps and high-resolution line profiles reveal the occurrence of several outflows. Cep A East is associated with multiple mass-loss processes: in particular, we detect a 0.6-pc jet-like structure which shows for the first time that the Cep A East young stellar objects are driving a collimated outflow moving towards the south.
The observed outflows show different clumps associated with definitely different H₂S/SO₂ integrated emission ratios, indicating that the gas chemistry in Cepheus A has been altered by the passage of shocks. H₂S appears to be more abundant than SO₂ in high-velocity clumps, in agreement with chemical models. However, we also find quite small H₂S linewidths, suggestive of regions where the evaporated H₂S molecules had enough time to slow down but not to freeze out on to dust grains. Finally, comparison between the line profiles indicates that the excitation conditions increase with the velocity, as expected for a propagation of collimated bow shocks.     

Molecular hydrogen emission in Cygnus A

We present J , H and K -band spectroscopy of Cygnus A, spanning 1.0–2.4 μm in the rest-frame and hence several rovibrational H₂, H recombination and [Fe  ii ] emission lines. The lines are spatially extended by up to 6 kpc from the nucleus, but their distinct kinematics indicate that the three groups (H, H₂ and [Fe  ii ]) are not wholly produced in the same gas. The broadest line, [Fe  ii ] λ 1.644, exhibits a non-Gaussian profile with a broad base (FWHM≃1040 km s⁻¹), perhaps because of the interaction with the radio source. Extinctions to the line-emitting regions substantially exceed earlier measurements based on optical H recombination lines.
Hard X-rays from the quasar nucleus are likely to dominate the excitation of the H₂ emission. The results of Maloney, Hollenbach & Tielens are thus used to infer the total mass of gas in H₂ v=1–0 S(1)-emitting clouds as a function of radius, for gas densities of 10³ and 10⁵ cm⁻³, and stopping column densities N _H=10²²–10²⁴ cm⁻². Assuming azimuthal symmetry, at least 2.3×10⁸ M_⊙ of such material is present within 5 kpc of the nucleus, if the line-emitting clouds see an unobscured quasar spectrum. Alternatively, if the bulk of the X-ray absorption to the nucleus inferred by Ueno et al. actually arises in a circumnuclear torus, the implied gas mass rises to ∼10¹⁰ M_⊙. The latter plausibly accounts for 10⁹ yr of mass deposition from the cluster cooling flow, for which within this radius.