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.     

Population III stars: hidden or disappeared?

A Population III/Population II transition from massive to normal stars is predicted to occur when the metallicity of the star-forming gas crosses the critical range   Z _cr= 10^−5±1 Z_⊙  . To investigate the cosmic implications of such a process, we use numerical simulations which follow the evolution, metal enrichment and energy deposition of both Population II and Population III stars. We find that: (i) due to inefficient heavy element transport by outflows and slow 'genetic' transmission during hierarchical growth, large fluctuations around the average metallicity arise; as a result, Population III star formation continues down to   z = 2.5  , but at a low peak rate of  10⁻⁵ M_⊙ yr⁻¹ Mpc⁻³  occurring at   z ≈ 6  (about 10⁻⁴ of the Population II one); and (ii) Population III star formation proceeds in an 'inside–out' mode in which formation sites are progressively confined to the periphery of collapsed structures, where the low gas density and correspondingly long free-fall time-scales result in a very inefficient astration. These conclusions strongly encourage deep searches for pristine star formation sites at moderate  (2 < z < 5)  redshifts where metal-free stars are likely to be hidden.     

The cooling of shock-compressed primordial gas

We find that at redshifts   z ≳ 10, HD  line cooling allows strongly shocked primordial gas to cool to the temperature of the cosmic microwave background (CMB). This temperature is the minimum value attainable via radiative cooling. Provided that the abundance of HD, normalized to the total number density, exceeds a critical level of  ∼10⁻⁸  , the CMB temperature floor is reached in a time which is short in comparison to the Hubble time. We estimate the characteristic masses of stars formed out of shocked primordial gas in the wake of the first supernovae, and resulting from the virialization of dark matter haloes during hierarchical structure formation to be  ∼10 M_⊙  . In addition, we show that cooling by HD enables the primordial gas in relic H  ii regions to cool to temperatures considerably lower than those reached via H₂ cooling alone. We confirm that HD cooling is unimportant in cases where the primordial gas does not go through an ionized phase, as in the formation process of the very first stars in   z ≳ 20  minihaloes of mass  ∼10⁶ M_⊙  .     

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⁻¹.     

Dust coagulation in star formation with different metallicities

Dust grains coagulate into larger aggregates in dense gas. This changes their size distribution and possibly affects the thermal evolution of star-forming clouds. We here investigate dust coagulation in collapsing pre-stellar cores with different metallicities by considering the thermal motions of grains. We show that coagulation does occur even at low metallicity  ∼10⁻⁶ Z_⊙  . However, we also find (i) that the H₂ formation rate on dust grains is reduced only after the majority of H₂ is formed and (ii) that the dust opacity is modified only after the core becomes optically thick. Therefore, we conclude that the effects of dust coagulation can safely be neglected in discussing the temperature evolution of the pre-stellar cores for any metallicity as long as the grain motions are thermal.     

On the first generation of stars

We argue that the first stars may have spanned the conventional mass range rather than be identified with the very massive objects  (∼100–10³ M_⊙)  favoured by numerical simulations. Specifically, we find that magnetic field generation processes acting in the first protostellar systems suffice to produce fields that exceed the threshold for magneto-rotational instability (MRI) to operate, and thereby allow the MRI dynamo to generate equipartition-amplitude magnetic fields on protostellar mass scales below  ∼50 M_⊙  . Such fields allow primordial star formation to occur at essentially any metallicity by regulating angular momentum transfer, fragmentation, accretion and feedback in much the same way as occurs in conventional molecular clouds.     

Binary formation with different metallicities: dependence on initial conditions

The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for  <10⁻⁴ Z_⊙  . With initial mass  <0.1 M_⊙  , if fragments are ejected in embryo from the host clouds by multibody interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.     

Universal gas density and temperature profile

We explore possibilities of collapse and star formation in Population III objects exposed to the external ultraviolet background (UVB) radiation. Assuming spherical symmetry, we solve self-consistently radiative transfer of photons, non-equilibrium H₂ chemistry and gas hydrodynamics. Although the UVB does suppress the formation of low-mass objects, the negative feedback turns out to be weaker than previously suggested. In particular, the cut-off scale of collapse drops significantly below the virial temperature T _vir∼10⁴ K at weak UV intensities ( J ₂₁≲10⁻²) , owing to both self-shielding of the gas and H₂ cooling. Clouds above this cut-off tend to contract highly dynamically, further promoting self-shielding and H₂ formation. For plausible radiation intensities and spectra, the collapsing gas can cool efficiently to temperatures well below 10⁴ K before rotationally supported and the final H₂ fraction reaches ∼ 10⁻³.
Our results imply that star formation can take place in low-mass objects collapsing in the UVB. The threshold baryon mass for star formation is ∼ 10⁹ M_⊙ for clouds collapsing at redshifts z ≲3 , but drops significantly at higher redshifts. In a conventional cold dark matter universe, the latter coincides roughly with that of the 1 σ density fluctuations. Objects near and above this threshold can thus constitute 'building blocks' of luminous structures, and we discuss their links to dwarf spheroidal/elliptical galaxies and faint blue objects. These results suggest that the UVB can play a key role in regulating the star formation history of the Universe.     

Impact of cosmic rays on Population III star formation

We explore the implications of a possible cosmic-ray (CR) background generated during the first supernova explosions that end the brief lives of massive Population III stars. We show that such a CR background could have significantly influenced the cooling and collapse of primordial gas clouds in minihaloes around redshifts of   z ∼ 15–20  , provided the CR flux was sufficient to yield an ionization rate greater than about 10⁻¹⁹ s⁻¹ near the centre of the minihalo. The presence of CRs with energies  ≲10⁷  eV would indirectly enhance the molecular cooling in these regions, and we estimate that the resulting lower temperatures in these minihaloes would yield a characteristic stellar mass as low as  ∼10 M_⊙  . CRs have a less-pronounced effect on the cooling and collapse of primordial gas clouds inside more massive dark matter haloes with virial masses  ≳10⁸ M_⊙  at the later stages of cosmological structure formation around   z ∼ 10–15  . In these clouds, even without CR flux the molecular abundance is already sufficient to allow cooling to the floor set by the temperature of the cosmic microwave background.     

Star formation and chemical evolution of damped Lyman α systems

In this paper, we investigate the star formation and chemical evolution of damped Lyman α systems (DLAs) based on the disc galaxy formation model developed by Mo, Mao & White. We propose that the DLAs are the central galaxies of less-massive dark haloes present at redshifts z ∼3, and they should inhabit haloes of moderately low circular velocity. The empirical Schmidt law of star formation rates, and closed box model of chemical evolution that an approximation known as instantaneous recycling is assumed, are adopted. In our models, when the predicted distribution of metallicity for DLAs is calculated, two cases are considered. One is that, using the closed-box model, empirical Schmidt law and star formation time, the distribution of metallicity can be directly calculated. The other is that, when the simple gravitational instability of a thin isothermal gas disc as first discussed by Toomre is considered, the star formation occurs only in the region where the surface density of gas satisfies the critical value, not everywhere of a gas disc. In this case, we first obtain the region where the star formation can occur by assuming that the disc has a flat rotation curve and rotational velocity is equal to the circular velocity of the surrounding dark matter halo, and then calculate the metallicity distribution as in case one. We assume that star formation in each DLA lasts for a period of 1 Gyr from redshifts z =3. There is only one output parameter in our models, i.e. the stellar yield, which relates to the time of star formation history and is obtained by normalizing the predicted distribution of metallicity to the mean value of 1/13 Z_⊙ as presented by Pettini et al.. The predicted metallicity distribution is consistent with the current (rather limited) observational data. A random distribution of galactic discs is taken into account.     

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.     

Semi-analytic simulations of galactic winds: volume filling factor, ejection of metals and parameter study

We present a semi-analytic treatment of galactic winds within high-resolution, large-scale cosmological N -body simulations of a Λ cold dark matter (ΛCDM) universe. The evolution of winds is investigated by following the expansion of supernova-driven superbubbles around the several hundred thousand galaxies that form in an approximately spherical region of space with diameter 52  h ⁻¹ Mpc and mean density close to the mean density of the universe. We focus our attention on the impact of winds on the diffuse intergalactic medium. Initial conditions for mass loss at the base of winds are taken from Shu, Mo & Mao. Results are presented for the volume filling factor and the mass fraction of the intergalactic medium (IGM) affected by winds, and their dependence on the model parameters is carefully investigated. The mass-loading efficiency of bubbles is a key factor to determine the evolution of winds and their global impact on the IGM: the higher the mass loading, the later the IGM is enriched with metals. Galaxies with 10⁹ < M _★ < 10¹⁰ M_⊙ are responsible for most of the metals ejected into the IGM at   z = 3  , while galaxies with   M _★ < 10⁹ M_⊙   give a non-negligible contribution only at higher redshifts, when larger galaxies have not yet assembled. We find a higher mean IGM metallicity than Lyα forest observations suggest, and we argue that the discrepancy may be explained by the high temperatures of a large fraction of the metals in winds, which may not leave detectable imprints in absorption in the Lyα forest.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

Wide triplets of galaxies: collapse on spatial scales of ∼1 Mpc

We study triple systems of galaxies with mean projected harmonic separation ≃0.6  h ⁻¹ Mpc     We call the systems 'wide triplets', in contrast to compact triplets with mean projected harmonic separation ≃0.04  h ⁻¹ Mpc, studied by Karachentsev et al. Data are taken for 108 wide triplets from a list compiled by Trofimov & Chernin; at least one-third of them are considered to be probably isolated physical systems. With typical crossing times of about the Hubble time, the wide triplets seem to be in a state of ongoing collapse. This is confirmed by a set of computer models which simulate well the observational characteristics of the ensemble of wide triplets. The simulations also give a statistical estimate of the total mass of a typical wide triplet: it proves to be ≃10¹³ M_⊙. This figure indicates that the dark matter mass is 15–30 times the mass of baryonic matter in the systems. The dynamics of wide triplets, as well as their dark matter content, provide new direct cosmological constraints by establishing that hierarchical evolution is occurring on a mass scale of ∼10¹³ M_⊙ and a spatial scale of ∼1 Mpc.     

The end of the dark ages in modified Newtonian dynamics

We study the evolution of a spherically symmetric density perturbation in the Modified Newtonian Dynamics model applied to the net acceleration over Hubble flow. The background cosmological model is a Λ-dominated, low-Ω_b Friedmann model with no cold dark matter. We include thermal processes and non-equilibrium chemical evolution of the collapsing gas. We find that under these assumptions the first low-mass objects  ( M ≤ 3 × 10⁴ M_⊙)  may collapse already for   z ∼ 30  , which is in quite good agreement with the recent Wilkinson Microwave Anisotropy Probe results. A lower value of a ₀ would lead to much slower collapse of such objects.     

Baryonic conversion tree: the global assembly of stars and dark matter in galaxies from the Sloan Digital Sky Survey

Using the spectroscopic sample of the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed approximately 80 per cent of their stars at   z > 1  , while for the less massive ones the value is only approximately 20 per cent. By comparing the BCT with the dark matter merger tree, we find indications that star formation efficiency at   z > 1  had to be approximately a factor of two higher than today (∼10 per cent) in galaxies with present-day stellar mass larger than  2 × 10¹¹ M_⊙  , if this early star formation occurred in the main progenitor. Therefore, the λ cold dark matter (LCDM) paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than  10¹¹ M_⊙  , efficient star formation seems to have been triggered at   z ∼ 0.2  . We show that there is a characteristic mass  ( M _*∼ 10¹⁰ M_⊙)  for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical models of galaxy formation.     

Expected planets in globular clusters

We argue that all transient searches for planets in globular clusters have a very low detection probability. Planets of low-metallicity stars typically do not reside at small orbital separations. The dependence of planetary system properties on metallicity is clearly seen when the quantity   I _e ≡ M _p[ a (1 − e )]²  is considered;   M _p, a   and e are the planet mass, semimajor axis and eccentricity, respectively. In high-metallicity systems, there is a concentration of systems at high and low values of I _e , with a low-populated gap near   I _e ∼ 0.3 M _J au²  , where M _J is Jupiter's mass. In low-metallicity systems, the concentration is only at the higher range of I _e , with a tail to low values of I _e . Therefore, it is still possible that planets exist around main-sequence stars in globular clusters, although at small numbers because of the low metallicity, and at orbital periods of ≳10 d. We discuss the implications of our conclusions on the role that companions can play in the evolution of their parent stars in globular clusters, for example, influencing the distribution of horizontal branch stars on the Hertzsprung–Russell diagram of some globular clusters, and in forming low-mass white dwarfs.     

The orbit and mass of the Sagittarius dwarf galaxy

Possible orbital histories of the Sgr dwarf galaxy are explored. A special-purpose N -body code is used to construct the first models of the Milky Way–Sgr dwarf system in which both the Milky Way and the Sgr dwarf are represented by full N -body systems and followed for a Hubble time. These models are used to calibrate a semi-analytic model of the Sgr dwarf's orbit that enables us to explore a wider parameter space than is accessible to the N -body models. We conclude that the extant data on the Sgr dwarf are compatible with a wide range of orbital histories. At one extreme the Sgr dwarf initially possesses ∼10¹¹ M_⊙ and starts from a Galactocentric distance R _D(0)≳200 kpc. At the other extreme the Sgr dwarf starts with ∼10⁹ M_⊙ and R _D(0)∼60 kpc, similar to its present apocentric distance. In all cases the Sgr dwarf is initially dark matter dominated and the current velocity dispersion of the Sgr dwarf's dark matter is tightly constrained to be 21±2 km s⁻¹. This number is probably compatible with the smaller measured dispersion of the Sgr dwarf's stars because of (i) the dynamical difference between dark and luminous matter, and (ii) velocity anisotropy.     

Dry mergers: a crucial test for galaxy formation

We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a nearly constant dry merger rate of  ∼6 × 10⁻⁵ Mpc⁻³ Gyr⁻¹  at   z ≤ 1  and a steep decline at larger z . Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early- and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of   M _{*, c} ∼ 6.3 × 10¹⁰ M_⊙  by mergers only. As a consequence of this transition, we find that around   M _{*, c}   , the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggests that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z . We however find that only 10–20 per cent of galaxies more massive than   M _{*, c}   experience dry major mergers within their last Gyr at any given redshift   z ≤ 1  .     

Dark matter subhaloes in numerical simulations

We use cosmological Λ cold dark matter (CDM) numerical simulations to model the evolution of the substructure population in 16 dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10¹¹ to  10¹⁴ M_⊙  at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation  d n /d v = (1/8)( v _cmax/ v _cmax,host)⁻⁴  . Halo to halo variance in the VDF is a factor of roughly 2 to 4. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for   r ≳ 0.3 r _vir  , where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease towards the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias towards more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power-law cosmologies of   P ∝ kⁿ   , and demonstrate that a steeper spectral index, n , results in significantly less substructure.