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.     

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.     

Occurrence of metal-free galaxies in the early Universe

The character of the first galaxies at redshifts z ≳ 10 strongly depends on their level of pre-enrichment, which is in turn determined by the rate of primordial star formation prior to their assembly. In order for the first galaxies to remain metal-free, star formation in minihaloes must be highly suppressed, most likely by H₂-dissociating Lyman–Werner (LW) radiation. We show that the build-up of such a strong LW background is hindered by two effects. First, the level of the LW background is self-regulated, being produced by the Population III (Pop III) star formation which it, in turn, suppresses. Secondly, the high opacity to LW photons which is built up in the relic H  ii regions left by the first stars acts to diminish the global LW background. Accounting for a self-regulated LW background, we estimate a lower limit for the rate of Pop III star formation in minihaloes at z ≳ 15. Further, we simulate the formation of a 'first galaxy' with virial temperature   T _vir≳ 10⁴ K  and total mass  ≳10⁸ M_⊙  at z ≳ 10, and find that complete suppression of previous Pop III star formation is unlikely, with stars of  ≳100 M_⊙  (Pop III.1) and  ≳10 M_⊙  (Pop III.2) likely forming. Finally, we discuss the implications of these results for the nature of the first galaxies, which may be observed by future missions such as the James Webb Space Telescope .     

Two populations of metal-free stars in the early Universe

We construct star formation histories at redshifts z ≳ 5 for two physically distinct populations of primordial, metal-free stars, motivated by theoretical and observational arguments that have hinted towards the existence of an intermediate stellar generation between Population III and Population I/II. Taking into account the cosmological parameters as recently revised by the Wilkinson Microwave Anisotropy Probe after three years of operation, we determine self-consistent reionization histories and discuss the resulting chemical enrichment from these early stellar generations. We find that the bulk of ionizing photons and heavy elements produced at high redshifts must have originated in Population II.5 stars, which formed out of primordial gas in haloes with virial temperatures ≳10⁴ K, and had typical masses ≳10 M_⊙. Classical Population III stars, formed in minihaloes and having masses ≳100 M_⊙, on the other hand, had only a minor impact on reionization and early metal enrichment. Specifically, we conclude that only ≃10 per cent by mass of metal-free star formation went into Population III.     

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.     

Non-linear evolution of suppressed dark matter primordial power spectra

We address the degree and rapidity of generation of small-scale power over the course of structure formation in cosmologies where the primordial power spectrum is strongly suppressed beyond a given wavenumber. We first summarize the situations where one expects such suppressed power spectra and point out their diversity. We then employ an exponential cut-off, which characterizes warm dark matter (WDM) models, as a template for the shape of the cut-off and focus on damping scales ranging from 10⁶ to  10⁹  h ⁻¹ M_⊙  . Using high-resolution simulations, we show that the suppressed part of the power spectrum is quickly (re)generated and catches up with both the linear and the non-linear evolution of the unsuppressed power spectrum. From   z = 2  onwards, a power spectrum with a primordial cut-off at  10⁹  h ⁻¹ M_⊙  becomes virtually indistinguishable from an evolved cold dark matter (CDM) power spectrum. An attractor such as that described in Zaldarriaga, Scoccimarro & Hui for power spectra with different spectral indices also emerges in the case of truncated power spectra. Measurements of   z ∼ 0  non-linear power spectra at  ∼100  h ⁻¹ kpc  cannot rule out the possibility of linear power spectra damped below  ∼10⁹  h ⁻¹ M_⊙  . Therefore, WDM or scenarios with similar features should be difficult to exclude in this way.     

Formation of HD molecules in merging dark matter haloes

HD molecules can be an important cooling agent of the primordial gas behind the shock waves originated through merging of the dark matter haloes at epochs when the first luminous objects formed. We study the necessary conditions for the HD cooling to switch on in the low-temperature range   T < 200 K  . We show that these conditions are fulfilled in merging haloes with total (dark matter and baryon) mass in excess of   M _cr∼ 10⁷[(1 + z )/20]⁻² M_⊙  . Haloes with masses   M > M _cr  may be the sites of low-mass star formation.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.     

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.     

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.     

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 census of nuclear activity of late-type galaxies in the Virgo cluster

The first spectroscopic census of active galactic nuclei (AGNs) associated with late-type galaxies in the Virgo cluster was carried out by observing 213 out of a complete set of 237 galaxies more massive than   M _dyn > 10^8.5 M_⊙  . Among them, 77 are classified as AGNs [including 21 transition objects, 47 low-ionization nuclear emission regions (LINERs) and nine Seyferts] and comprise 32 per cent of the late-type galaxies in Virgo. Due to spectroscopic incompleteness, at most 21 AGNs are missed in the survey, so that the fraction would increase up to 41 per cent. Using corollary near-infrared observations that enable us to estimate galaxy dynamical masses, it is found that AGNs are hosted exclusively in massive galaxies, i.e.   M _dyn≳ 10¹⁰ M_⊙  . Their frequency increases steeply with the dynamical mass from zero at   M _dyn≈ 10^9.5 M_⊙  to virtually 1 at   M _dyn > 10^11.5 M_⊙  . These frequencies are consistent with those of low-luminosity AGNs found in the general field by the Sloan Digital Sky Survey. Massive galaxies that harbour AGNs commonly show conspicuous r -band star-like nuclear enhancements. Conversely, they often, but not necessarily, contain massive bulges. A few well-known AGNs (e.g. M61, M100, NGC 4535) are found in massive Sc galaxies with little or no bulge. The AGN fraction seems to be only marginally sensitive to galaxy environment. We infer the black hole masses using the known scaling relations of quiescent black holes. No black holes lighter than  ∼10⁶ M_⊙  are found active in our sample.     

Unresolved emission and ionized gas in the bulge of M31

We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 10²⁷ erg s⁻¹ L⁻¹_⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 10⁶ M_⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L _X/SFR  (star formation rate) ratio equals  ∼9 × 10³⁸ (erg s⁻¹)(M_⊙ yr⁻¹)⁻¹  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.     

On the origin of high-velocity runaway stars

We explore the hypothesis that some high-velocity runaway stars attain their peculiar velocities in the course of exchange encounters between hard massive binaries and a very massive star (either an ordinary  50–100 M_⊙  star or a more massive one, formed through runaway mergers of ordinary stars in the core of a young massive star cluster). In this process, one of the binary components becomes gravitationally bound to the very massive star, while the second one is ejected, sometimes with a high speed. We performed three-body scattering experiments and found that early B-type stars (the progenitors of the majority of neutron stars) can be ejected with velocities of  ≳200–400 km s⁻¹  (typical of pulsars), while  3–4 M_⊙  stars can attain velocities of  ≳300–400 km s⁻¹  (typical of the bound population of halo late B-type stars). We also found that the ejected stars can occasionally attain velocities exceeding the Milky Ways's escape velocity.     

Cosmological origin of the lowest metallicity halo stars

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars that are likely to be found within our Galaxy today. By PopIII we refer to stars formed at large redshift ( z >4), with low metallicities ([ Z /Z_⊙]<−2.5) and in small systems (total mass ≲ 2×10⁸ M_⊙) that are extremely sensitive to stellar feedback, and which through a prescribed merging history end up becoming part of the Milky Way today. An analytic, extended Press–Schechter formalism is used to obtain the mass functions of haloes which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modelled as a mini-galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disc instabilities, all of which determine the fraction of the baryons that are subject to star formation. The use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, and give accurate estimates of their expected numbers today and their location in L /L_⊙ versus T /K Hertzsprung–Russell (HR) diagrams. A first comparison with observational data suggests that the initial mass function (IMF) of the first stars was increasingly high-mass weighted towards high redshifts, levelling off at z ≳9 at a characteristic stellar mass scale m _s=10–15 M_⊙.     

Angular momentum and clustering properties of early dark matter haloes

In this paper, we study the angular momentum properties of simulated dark matter haloes at high redshifts that likely host the first stars in the Universe. Calculating the spin distributions of these  10⁶– 10⁷ M_⊙  haloes in redshift slices from   z = 15  to 6, we find that they are well fit by a lognormal distribution as is found for lower redshift and more massive haloes in earlier work. We find that both the mean value of the spin and dispersion are largely unchanged with redshift for all haloes. Our key result is that subsamples of low- and high-spin, 10⁶ and  10⁷ M_⊙  , haloes show difference in clustering strength. In both mass bins, higher spin haloes are more strongly clustered in concordance with a tidal torquing picture for the growth of angular momentum in dark matter haloes in the cold dark matter paradigm.     

A search for the progenitors of two Type Ia Supernovae in NGC 1316

Recent evidence of a young progenitor population for many Type Ia supernovae (SNe Ia) raises the possibility that evolved intermediate-mass progenitor stars may be detected in pre-explosion images. NGC 1316, a radio galaxy in the Fornax cluster, is a prolific producer of SNe Ia, with four detected since 1980. We analyse Hubble Space Telescope ( HST ) pre-explosion images of the sites of two of the SNe Ia that exploded in this galaxy, SN2006dd (a normal Type Ia) and SN2006mr (likely a subluminous, 1991bg-like, SN Ia). Astrometric positions are obtained from optical and near-infrared ground-based images of the events. We find no candidate point sources at either location, and set upper limits on the flux in B, V and I from any such progenitors. We also estimate the amount of extinction that could be present, based on analysis of the surface-brightness inhomogeneities in the HST images themselves. At the distance of NGC 1316, the limits correspond to absolute magnitudes of  ∼−5.5, −5.4  and −6.0 mag in   M _B , M _V   and   M _I   , respectively. Comparison to stellar evolution models argues against the presence at the supernova sites, 3 yr prior to the explosion, of normal stars with initial masses  ≳6 M_⊙  at the tip of their asymptotic-giant branch (AGB) evolution, young post-AGB stars that had initial masses  ≳4 M_⊙  and post-red giant stars of initial masses  ≳9 M_⊙  .     

On the formation of massive galaxies: a simultaneous study of number density, size and intrinsic colour evolution in GOODS

The evolution of number density, size and intrinsic colour is determined for a volume-limited sample of visually classified early-type galaxies selected from the Hubble Space Telescope /Advanced Camera for Surveys images of the Great Observatories Origins Deep Survey (GOODS) North and South fields (version 2). The sample comprises 457 galaxies over 320 arcmin² with stellar masses above  3 × 10¹⁰ M_⊙  in the redshift range  0.4 < z < 1.2  . Our data allow a simultaneous study of number density, intrinsic colour distribution and size. We find that the most massive systems  (≳3 × 10¹¹ M_⊙)  do not show any appreciable change in comoving number density or size in our data. Furthermore, when including the results from 2dF galaxy redshift survey, we find that the number density of massive early-type galaxies is consistent with no evolution between   z = 1.2  and 0, i.e. over an epoch spanning more than half of the current age of the Universe. We find large discrepancies between the predictions of semi-analytic models. Massive galaxies show very homogeneous intrinsic colour distributions, with nearly flat radial colour gradients, but with a significant negative correlation between stellar mass and colour gradient, such that red cores appear predominantly in massive galaxies. The distribution of half-light radii – when compared to   z ∼ 0  and   z > 1  samples – is compatible with the predictions of semi-analytic models relating size evolution to the amount of dissipation during major mergers.     

Neutralino annihilations and the gas temperature in the dark ages

Assuming that the dark matter is entirely made up of neutralinos, we re-visit the role of their annihilation on the temperature of diffuse gas in the high-redshift universe  ( z > 10)  , before the formation of luminous structures. We consider neutralinos of particle mass 36 and 100 GeV. The former is able to produce  ∼7  e ⁻ e ⁺  particles per annihilation through the fremionic channel, and the latter ∼53 particles assuming a purely bosonic channel. High-energy   e ⁻ e ⁺  particles up-scatter the cosmic microwave background (CMB) photons into higher energies via the inverse-Compton scattering. The process produces a power-law   e ⁻ e ⁺  energy spectrum of index −1 in the energy range of interest, independent of the initial energy distribution. The corresponding energy spectrum of the up-scattered photons is a power law of index −1/2, if absorption by the gas is not included. The scattered photons photoheat the gas by releasing electrons which deposit a fraction (14 per cent) of their energy as heat into the ambient medium. For uniformly distributed neutralinos, the heating is insignificant. The effect is greatly enhanced by the clumping of neutralinos into dense haloes. We use a time-dependent clumping model which takes into account the damping of density fluctuations on mass-scales smaller than  ∼10⁻⁶ M_⊙  . With this clumping model, the heating mechanism boosts the gas temperature above that of the CMB after a redshift of   z ∼ 30  . By   z ≈ 10  , the gas temperature is nearly 100 times its temperature when no heating is invoked. Similar increase is obtained for the two neutralino masses considered.     

X-ray analysis of the cluster 3C 129

The cluster 3C 129 is classified as a rich cluster. An analysis of the properties of the cluster 3C 129 from ROSAT PSPC and HRI, Einstein IPC, and EXOSAT ME observations is presented. The mean temperature from a joint fit of the ROSAT PSPC and EXOSAT ME data is 5.5(±0.2) keV. The luminosity is 0.6×10⁴⁴ erg s⁻¹ in 0.2–2.4 keV and 2.7×10⁴⁴ erg s⁻¹ in 0.2–10 keV. We find a cooling flow with a rate of ∼84 M_⊙ yr⁻¹. The central gas density is 6×10⁻³ cm⁻³, and the ICM mass is 3.6×10¹³ M_⊙. The total cluster mass is ∼5×10¹⁴ M_⊙. The X-ray morphology shows an east–west elongation, which is evidence for a recent merger event. The radio source 3C 129.1 is located near the X-ray centre. Another cluster member galaxy (the radio galaxy 3C 129) is a prototype of head-tailed radio galaxies, and is located in the west part of the cluster. The tail points along the gradient of intracluster gas pressure. There are no significant point X-ray sources associated with the AGNs of the two radio galaxies.