In this paper we focus on understanding the physical processes that lead to stable or unstable ionization fronts (I-fronts) observed in simulations of moving black holes (BHs). The front instability may trigger bursts of gas accretion, rendering the BH significantly more luminous than at steady-state. We perform a series of idealized three dimensional radiation hydrodynamics simulations resolving the I-fronts around BHs of mass $M_\mathrm{BH}$ and velocity $v_\infty$ accreting from a medium of density $n_\mathrm{H}$. The I-front, with radius $R_\mathrm{I}$, transitions from D-type to R-type as the BH velocity becomes larger than a critical value $v_\mathrm{R}\sim 40\,\mathrm{km/s}$. The D-type front is preceded by a bow-shock of thickness $\Delta R_\mathrm{I}$ that decreases as $v_\infty$ approaches $v_\mathrm{R}$. We find that both D-type and R-type fronts can be unstable given the following two conditions: i) for D-type fronts the shell thickness must be $\Delta R_\mathrm{I}/R_\mathrm{I}<0.05$ (i.e., $v_\infty \gtrsim 20\,\mathrm{km/s}$.), while no similar restriction holds for R-type fronts; ii) the temperature jump across the I-front must be $T_\mathrm{II}/T_\mathrm{I}>3$. This second condition is satisfied if $T_\mathrm{I}<5000\,\mathrm{K}$ or if $n_\mathrm{H}\,M_\mathrm{BH} \gtrsim 10^6\,M_\odot\,\mathrm{cm^{-3}}$. Due to X-ray pre-heating typically $T_\mathrm{I} \sim 10^4\,\mathrm{K}$, unless the D-type shell is optically thick to X-rays, which also happens when $n_\mathrm{H}\,M_\mathrm{BH}$ is greater than a metallicity-dependent critical value. We thus conclude that I-fronts around BHs are unstable only for relatively massive BHs moving trough very dense molecular clouds. We briefly discuss the observational consequences of the X-ray luminosity bursts likely associated with this instability.

Hyper-accretion occurs when the gas inflow rate onto a black hole (BH) is so high that the radiative feedback cannot reverse the accretion flow. This extreme process is a promising mechanism for the rapid growth of seed BHs in the early universe, which can explain high-redshift quasars powered by billion solar mass BHs. In theoretical models, spherical symmetry is commonly adopted for hyper-accretion flow, however its sustainability on timescales corresponding to the BH growth has not been addressed yet. Here we show that stochastic interactions between the ionizing radiation from the BH and non-uniform accretion flow can lead to a formation of rotating gas disk around the BH. Once the disk is formed, the supply of gas to the BH preferentially occurs via biconical dominated accretion flow (BDAF) perpendicular to the disk, avoiding the centrifugal barrier of the disk. BDAFs from opposite directions collide in the vicinity of the BH supplying high-density, low angular momentum gas to the BH, whereas most of the gas with non-negligible angular momentum is deflected to the rotationally supported outflowing decretion disk. The disk becomes reinforced progressively as more mass from the biconical flow transfers to the disk and some of the outflowing gas from the disk is redirected to the biconical accretion funnels through a meridional structure. This axisymmetric hydrodynamic structure of a BDAF and decretion disk continues to provide uninterrupted flow of high-density gas to the BH.

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 $100\times$ 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 $100\times$ denser than in today's Universe, which leads to a significant production of old globular clusters progenitors.

We present radiation-magneto-hydrodynamic simulations of star formation in self-gravitating, turbulent molecular clouds, modeling the formation of individual massive stars, including their UV radiation feedback. The set of simulations have cloud masses between $m_{\rm gas}=10^3$~M$_\odot$ to $3 \times 10^5$~M$_\odot$ and gas densities typical of clouds in the local universe ($\overline n_{\rm gas} \sim 1.8\times 10^2$~cm$^{-3}$) and 10$\times$ and 100$\times$ denser, expected to exist in high-redshift galaxies. The main results are: {\it i}) The observed Salpeter power-law slope and normalisation of the stellar initial mass function at the high-mass end can be reproduced if we assume that each star-forming gas clump (sink particle) fragments into stars producing on average a maximum stellar mass about $40\%$ of the mass of the sink particle, while the remaining $60\%$ is distributed into smaller mass stars. Assuming that the sinks fragment according to a power-law mass function flatter than Salpeter, with log-slope $0.8$, satisfy this empirical prescription. {\it ii}) The star formation law that best describes our set of simulation is $d\rho_*/dt \propto \rho_{gas}^{1.5}$ if $\overline n_{gas}<n_{cri}\approx 10^3$~cm$^{-3}$, and $d\rho_*/dt \propto \rho_{\rm gas}^{2.5}$ otherwise. The duration of the star formation episode is roughly $6$ cloud's sound crossing times (with $c_s=10$~km/s). {\it iii}) The total star formation efficiency in the cloud is $f_*=2\% (m_{\rm gas}/10^4~M_\odot)^{0.4}(1+\overline n_{\rm gas}/n_{\rm cri})^{0.91}$, for gas at solar metallicity, while for metallicity $Z<0.1$~Z$_\odot$, based on our limited sample, $f_*$ is reduced by a factor of $\sim 5$. {\it iv)} The most compact and massive clouds appear to form globular cluster progenitors, in the sense that star clusters remain gravitationally bound after the gas has been expelled.

If primordial black holes of $\mathcal{O}(1-100) M_{\odot}$ constitute a significant portion of the dark matter in the Universe, they should be very abundant in our Galaxy. We present here a detailed analysis of the radio and X-ray emission that these objects are expected to produce due to the accretion of gas from the interstellar medium. With respect to previous studies, we relax the assumption of a monochromatic mass function, and introduce an improved treatment of the physics of gas accretion onto isolated, moving compact objects, based on a set of state-of-the-art numerical simulations. By comparing our predictions with known radio and X-ray sources in the Galactic center region, we show that the maximum relic density of primordial black holes in the mass range of interest is $\sim 10^{-3}$ smaller than that of dark matter. The new upper bound is two orders of magnitude stronger with respect to previous results, based on a conservative phenomenological treatment of the accretion physics. We also provide a comprehensive critical discussion on the reliability of this bound, and on possible future developments in the field. We argue in particular that future multi-wavelength searches will soon start to probe the galactic population of astrophysical black holes.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes).

Stellar halos observed around normal galaxies are extended and faint stellar structures formed by debris of tidally disrupted dwarf galaxies accreted over time by the host galaxy. Around dwarf galaxies, these stellar halos may not exist if all the accreted satellites are dark halos without stars. However, if a stellar halo is found in sufficiently small mass dwarfs, the whole stellar halo is composed of tidal debris of fossil galaxies, and we refer to it as ghostly halo. Fossil galaxies are called so because they formed most of their stars before the epoch of reionization, and have been identified as the ultra-faint dwarf galaxies found around the Milky Way and M31. In this paper we carry out semi-analytical simulations to characterize the sizes and stellar masses of ghostly stellar halos in dwarf galaxies as a function of their dark matter halo mass. By comparing the models to observations of six isolated dwarf galaxies in the Local Group showing evidence of extended stellar halos, we are able to constrain the star formation efficiency in fossil galaxies. We find that at redshift z~6, dark matter halos in the mass range 10^7-10^9 M_sun have a mean star formation efficiency f_* = M_*/M_{dm} ~ 0.05%-0.1% nearly constant as a function of the dark matter halo mass. There is a tentative indication of a sharp increase of f_* at M_{dm} ~ 10^6-10^7 M_sun, but it is very uncertain and based on only one dwarf galaxy (Leo~T).

We present the methodology of a photon-conserving, spatially-adaptive, ray-tracing radiative transfer algorithm, designed to run on multiple parallel Graphic Processing Units (GPUs). Each GPU has thousands computing cores, making them ideally suited to the task of tracing independent rays. This ray-tracing implementation has speed competitive with approximate momentum methods, even with thousands of ionization sources, without sacrificing accuracy and resolution. Here, we validate our implementation with the selection of tests presented in the "cosmological radiative transfer codes comparison project," to demonstrate the correct behavior of the code. We also present a selection of benchmarks to demonstrate the performance and computational scaling of the code. As expected, our method scales linearly with the number of sources and with the square of the dimension of the 3D computational grid. Our current implementation is scalable to an arbitrary number of nodes possessing GPUs, but is limited to a uniform resolution 3D grid. Cosmological simulations of reionization with tens of thousands of radiation sources and intergalactic volumes sampled with 1024$^3$ grid points take about 30 days on 64 GPUs to reach complete reionization.

Evidence for dust around supermassive black holes (SMBHs) in the early Universe is strongly suggested by recent observations. However, the accretion mechanism of SMBHs in dusty gas is not well understood yet. We investigate the growth of intermediate-mass black-holes (IMBHs) of $\sim 10^{5}~M_{\odot}$ in dusty clouds by using one-dimensional radiative-hydrodynamics simulations. We find that the accretion of dusty gas onto IMBHs proceeds gently with small fluctuations of the accretion rate, whereas that of pristine gas causes more violent periodic bursts. At dust-to-gas mass ratios similar to the solar neighborhood, the time averaged luminosity becomes smaller than that for primordial gas by one order of magnitude and the time-averaged Eddington ratio ranges from $\sim 10^{-4}$ to $\sim 10^{-2}$ in clouds with initial gas densities of $n_{\rm H} = 10 - 1000~\rm cm^{-3}$. Our calculations show that the effect of dust opacity alone is secondary compared to the radiation pressure on dust in regulating the BH growth. We also derive spectral energy distributions at IR bands by calculating dust thermal emission and show that the flux ratio between $\lambda \lesssim 20~\rm \mu m$ and $\gtrsim 100~\rm \mu m$ is closely related to the Eddington ratio. Thermal emission from hot dust near the BH dominates only during the high accretion phase, producing higher flux density at $\lesssim 20~\rm \mu m$. Therefore, we suggest that the combinations of MIR observations by JWST and FIR observation by ALMA or Spitzer can be used to estimate the Eddington ratio of massive BHs.

We present theoretical constraints for the formation of the newly discovered dark star clusters (DSCs) with high mass-to-light (M/L) ratios, from Taylor et al (2015). These compact stellar systems photometrically resemble globular clusters (GCs) but have dynamical M/L ratios of ~ 10 - 100, closer to the expectations for dwarf galaxies. The baryonic properties of the dark star clusters (DSCs) suggest their host dark matter halos likely virialized at high redshift with M > 10^8 M_sun. We use a new set of high-resolution N-body simulations of Centaurus A to determine if there is a set of z=0 subhalos whose properties are in line with these observations. While we find such a set of subhalos, when we extrapolate the dark matter density profiles into the inner 20 pc, no dark matter halo associated with Centaurus A in our simulations, at any redshift, can replicate the extremely high central mass densities of the DSCs. Among the most likely options for explaining 10^5 - 10^7 M_sun within 10 pc diameter subhalos is the presence of a central massive black hole. We, therefore, propose that the DSCs are remnant cusps of stellar systems surrounding the central black holes of dwarf galaxies which have been almost completely destroyed by interactions with Centaurus A.

Pop III stars are typically massive stars of primordial composition forming at the centers of the first collapsed dark matter structures. Here we estimate the optimal X-ray emission in the early universe for promoting the formation of Pop III stars. This is important in determining the number of dwarf galaxies formed before reionization and their fossils in the local universe, as well as the number of intermediate-mass seed black holes. A mean X-ray emission per source above the optimal level reduces the number of Pop III stars because of the increased Jeans mass of the intergalactic medium (IGM), while a lower emission suppresses the formation rate of H2 preventing or delaying star formation in dark matter minihalos above the Jeans mass. The build up of the H2 dissociating background is slower than the X-ray background due to the shielding effect of resonant hydrogen Lyman lines. Hence, the nearly unavoidable X-ray emission from supernova remnants of Pop III stars is sufficient to boost their number to few tens per comoving Mpc^3 by redshift z~15. We find that there is a critical X-ray to UV energy ratio emitted per source that produces a universe where the number of Pop III stars is largest: 400 per comoving- Mpc^3. This critical ratio is very close to the one provided by 20-40 M_sun Pop III stars exploding as hypernovae. High mass X-ray binaries in dwarf galaxies are far less effective at increasing the number of Pop III stars than normal supernova remnants, we thus conclude that supernovae drove the formation of Pop III stars.

In this paper, the first in a series on galaxy formation before reionization, we focus on understanding what determines the size and morphology of stellar objects in the first low mass galaxies, using parsec- scale cosmological simulations performed with an adaptive mesh hydrodynamics code. Although the dense gas in which stars are formed tends to have a disk structure, stars are found in spheroids with little rotation. Halos with masses between 10^6 M_sun and 5x10^8 M_sun form stars stochastically, with stellar masses in the range 10^4 M_sun to 2x10^6 M_sun. Nearly independent of stellar mass, we observe a large range of half-light radii for the stars, from a few parsecs to a few hundred parsecs and surface brightnesses and mass-to-light ratios ranging from those typical of globular clusters to ultra-faint dwarfs. In our simulations, stars form in dense stellar clusters with high gas-to-star conversion efficiencies and rather uniform metallicities. A fraction of these clusters remain bound after the gas is removed by feedback, but others are destroyed, and their stars, which typically have velocity dispersions of 20 to 40 km/s, expand until they become bound by the dark matter halo. We thus speculate that the stars in ultra-faint dwarf galaxies may show kinematic and chemical signatures consistent with their origin in a few distinct stellar clusters. On the other hand, some globular clusters may form at the center of primordial dwarf galaxies and may contain dark matter, perhaps detectable in the outer parts.

We present semi-analytic models of the epoch of reionization focusing on the differences between continuous and bursty star formation (SF). Our model utilizes physically motivated analytic fits to 1D radiative transfer simulations of HII regions around dark matter halos in a representative cosmic volume. Constraining our simulations with observed and extrapolated UV luminosity functions of high redshift galaxies, we find that for a fixed halo mass, stellar populations forming in bursty models produce larger HII regions which leave behind long-lived relic HII regions which are able to maintain partial ionization in the intergalactic medium (IGM) in a manner similar to an early X-ray background. The overall effect is a significant increase in the optical depth of the IGM, $\tau_e$, and a milder increase of the redshift of reionization. To produce $\tau_e =0.066$ observed by Planck and complete reionization by redshift $z_{\rm re}\sim 6$, models with bursty SF require an escape fraction $f_{\rm esc}\sim 2\%-10\%$ that is 2-10 times lower than $f_{\rm esc}\sim 17\%$ found assuming continuous SF and is consistent with upper limits on $f_{\rm esc}$ from observations at $z=0$ and $z\sim 1.3-6$. The ionizing photon budget needed to reproduce the observed $\tau_e$ and $z_{\rm re}$ depends on the period and duty cycle of the bursts of SF and the temperature of the neutral IGM. These results suggest that the tension between observed and predicted ionizing photon budget for reionization can be alleviated if reionization is driven by short bursts of SF, perhaps relating to the formation of Population~III stars and compact star clusters such as proto-globular clusters.

We examine radiation-regulated accretion onto intermediate-mass and massive black holes (BHs) embedded in a bulge component. Using spherically symmetric one-dimensional radiation-hydrodynamics simulations, we track the growth of BHs accreting from a cold, neutral gas reservoir with temperature T=10^4 K. We find that the accretion rate of BHs embedded in bulges is proportional to r_{B,eff}/r_B, where r_{B,eff} is the increased effective Bondi radius that includes the gravitational potential of the bulge, and r_B is the Bondi radius of the BH. The radiative feedback from the BH suppresses the cold accretion rate to ~1 percent of the Bondi rate when a bulge is not considered. However, we find that the BH fueling rate increases rapidly when the bulge mass M_bulge is greater than the critical value of 10^6 M_sun and is proportional to M_bulge. Since the critical bulge mass is independent of the central BH mass M_{BH}, the growth rate of BHs with masses of 10^2, 10^4, and 10^6 M_sun exhibits distinct dependencies on the bulge-to-BH mass ratio. Our results imply that light seed BHs (<= 10^2 M_sun) which might be the remnants of the Pop III stars, cannot grow through accretion coevally with the early assembly of the bulge of the host galaxies until the bulge reaches the critical mass. However, massive BH seeds (>= 10^5 M_sun) that may form via direct collapse, are more likely to be embedded in a supercritical bulge and thus can grow efficiently coupling to the host galaxies and driving the early evolution of the M_{BH}-$\sigma$ relationship.

The study of the Cosmic Near-Infrared Background (CIB) light after subtraction of resolved sources can push the limits of current observations and infer the level of galaxy and black hole activity in the early universe. However, disentangling the relative contribution from low- and high-redshift sources is not trivial. Spatial fluctuations of the CIB exhibit a clustering excess at angular scales $\sim 1^\circ$ whose origin has not been conclusively identified. We explore the likelihood that this signal is dominated by emission from galaxies and accreting black holes in the early Universe. We find that, if the first small mass galaxies have a normal IMF, the light of their ageing stars (fossils) integrated over cosmic time contributes a comparable amount to the CIB as their pre-reionization progenitors. However, the measured fluctuation signal is too large to be produced by galaxies at redshifts $z>8$ unless their star formation efficiencies are much larger than those inferred from the observed Lyman-dropout population. In order to produce the observed level of CIB fluctuation without violating constraints from galaxy counts and the electron optical depth of the IGM, minihalos at $z>12$ must form stars with efficiency $f_\star \gtrsim 0.1$ and, although a top-heavy IMF is preferred, have a very low escape fraction of ionizing radiation, $f_{\rm esc}<0.01$. If instead the CIB fluctuations are produced by high-$z$ black holes, one requires vigorous accretion in the early universe reaching $\rho_{\rm acc} \gtrsim 10^5M_\odot{\rm Mpc^{-3}}$ by $z\simeq 10$. This growth must stop by $z \sim 6$ and be significantly obscured not to overproduce the soft cosmic X-ray background (CXB) and its observed coherence with the CIB. We therefore find the range of suitable possibilities at high-$z$ to be narrow, but could possibly be widened by including additional physics and evolution at those epochs.

We use N-body simulations of dark matter haloes in cold dark matter (CDM) and a large set of different warm dark matter (WDM) cosmologies to demonstrate that the spherically averaged density profile of dark matter haloes has a shape that depends on the power spectrum of matter perturbations. Density profiles are steeper in WDM but become shallower at scales less than one percent of the virial radius. Virialization isotropizes the velocity dispersion in the inner regions of the halo but does not erase the memory of the initial conditions in phase space. The location of the observed deviations from CDM in the density profile and in phase space can be directly related to the ratio between the halo mass and the filtering mass and are most evident in small mass haloes, even for a 34 keV thermal relic WDM. The rearrangement of mass within the haloes supports analytic models of halo structure that include angular momentum. We also find evidence of a dependence of the slope of the inner density profile in CDM cosmologies on the halo mass with more massive haloes exhibiting steeper profiles, in agreement with the model predictions and with previous simulation results. Our work complements recent studies of microhaloes near the filtering scale in CDM and strongly argue against a universal shape for the density profile.

We investigate the role of Compton heating in radiation-regulated accretion on to black holes from a neutral dense medium using 1D radiation-hydrodynamic simulations. We focus on the relative effects of Compton-heating and photo-heating as a function of the spectral slope {\alpha}, assuming a power-law spectrum in the energy range of 13.6 eV--100 keV. While Compton heating is dominant only close to the black hole, it can reduce the accretion rate to 0.1 % ($l \propto \dot{m}^2$ model)--0.01 % ($l \propto \dot{m}$ model) of the Bondi accretion rate when the BH radiation is hard ({\alpha} ~ 1), where $l$ and $\dot{m}$ are the luminosity and accretion rate normalised by Eddington rates, respectively. The oscillatory behaviour otherwise typically seen in simulations with {\alpha} > 1, become suppressed when {\alpha} ~ 1 only for the $l \propto \dot{m}$ model. The relative importance of the Compton heating over photo-heating decreases and the oscillatory behaviour becomes stronger as the spectrum softens. When the spectrum is soft ({\alpha} > 1.5), photo-heating prevails regardless of models making the effect of Compton heating negligible. On the scale of the ionization front, where the gas supply into the Str\"omgren sphere from large scale is regulated, photo-heating dominates. Our simulations show consistent results with the advection-dominated accretion flow ($l \propto \dot{m}^2$) where the accretion is inefficient and the spectrum is hard ({\alpha} ~ 1).

Globular clusters are unique tracers of ancient star formation. We determine the formation efficiencies of globular clusters across cosmic time by modeling the formation and dynamical evolution of the globular cluster population of a Milky Way type galaxy in hierarchical cosmology, using the merger tree from the Via Lactea II simulation. All of the models are constrained to reproduce the observed specific frequency and initial mass function of globular clusters in isolated dwarfs. Globular cluster orbits are then computed in a time varying gravitational potential after they are either accreted from a satellite halo or formed in situ, within the Milky Way halo. We find that the Galactocentric distances and metallicity distribution of globular clusters are very sensitive to the formation efficiencies of globular clusters as a function of redshift and halo mass. Our most accurate models reveal two distinct peaks in the globular cluster formation efficiency at z~2 and z~7-12 and prefer a formation efficiency that is mildly increasing with decreasing halo mass, the opposite of what expected for feedback-regulated star formation. This model accurately reproduces the positions, velocities, mass function, metallicity distribution, and age distribution of globular clusters in the Milky Way and predicts that ~ 40% formed in situ, within the Milky Way halo, while the other ~ 60% were accreted from about 20 satellite dwarf galaxies with Vc > 30 km/s, and about 29% or all globular clusters formed at redshifts z > 7. These results further strengthen the notion that globular cluster formation was an important mode of star formation in high-redshift galaxies and likely played a significant role in the reionization of the intergalactic medium

A spatial clustering signal has been established in Spitzer/IRAC measurements of the unresolved Cosmic near-Infrared Background (CIB) out to large angular scales, ~1 deg. This CIB signal, while significantly exceeding the contribution from the remaining known galaxies, was further found to be coherent at a highly statistically significant level with the unresolved soft Cosmic X-ray Background (CXB). This measurement probes the unresolved CXB to very faint source levels using deep near-IR source subtraction. We study contributions from extragalactic populations at low to intermediate redshifts to the measured positive cross-power signal of the CIB fluctuations with the CXB. We model the X-ray emission from AGN, normal galaxies and hot gas residing in virialized structures, calculating their CXB contribution including their spatial coherence with all infrared emitting counterparts. We use a halo model framework to calculate the auto and cross-power spectra of the unresolved fluctuations based on the latest constraints of the halo occupation distribution and the biasing of AGN, galaxies and diffuse emission. At small angular scales (<1'), the 4.5mic vs 0.5-2 keV coherence can be explained by shot noise from galaxies and AGN. However, at large angular scales (~10 arcmin) we find that the net contribution from the modeled populations is only able to account for ~3% of the measured CIBxCXB cross-power. The discrepancy suggests that the CIBxCXB signal originates from the same unknown source population producing the CIB clustering signal out to ~1 deg.

We examine the role of ionizing radiation emitted from black holes (BHs) in suppressing the growth of the Rayleigh-Taylor instability (RTI) across the ionization front (I-front) that forms when the gas fueling the BH is neutral. We use radiation-hydrodynamic simulations to show that the RTI is suppressed for non-accelerating fronts on all scales resolved in our simulations. A necessary condition for the stability of the I-front is that the radius of the Str\"omgren sphere is larger than the Bondi radius. When this condition is violated the I-front collapses producing an accretion luminosity burst. Transient growth of the RTI occurs only during the accretion burst when the effective acceleration in the frame of reference of the I-front increases significantly due to the rapid expansion of the Str\"omgren sphere.

The linear growth rate of the Rayleigh-Taylor instability (RTI) at ionization fronts is investigated via perturbation analysis in the limit of incompressible fluids. In agreement with previous numerical studies is found that absorption of ionizing radiation inside the HII region due to hydrogen recombinations suppresses the growth of instabilities. In the limit of a large density contrast at the ionization front the RTI growth rate has the simple analytical solution n=-nur+(nur^2+gk)^(1/2), where nur is the hydrogen recombination rate inside the HII region, k is the perturbation's wavenumber and g is the effective acceleration in the frame of reference of the front. Therefore, the growth of surface perturbations with wavelengths lambda >> lambda_{cr} = 2\pi g/nur^2 is suppressed by a factor (lambda_{cr}/4lambda)^(1/2) with respect to the non-radiative incompressible RTI. Implications on stellar and black hole feedback are briefly discussed.

We investigate the claim that the largest subhaloes in high resolution dissipationless cold dark matter (CDM) simulations of the Milky Way are dynamically inconsistent with observations of its most luminous satellites. We find that the inconsistency is largely attributable to the large values of \sigma_8 and n_s adopted in the discrepant simulations producing satellites that form too early and therefore are too dense. We find the tension between observations and simulations adopting parameters consistent with WMAP9 is greatly diminished making the satellites a sensitive test of CDM. We find the Via Lactea II halo to be atypical for haloes in a WMAP3 cosmology, a discrepancy that we attribute to its earlier formation epoch than the mean for its mass. We also explore warm dark matter (WDM) cosmologies for 1--4 keV thermal relics. In 1 keV cosmologies subhaloes have circular velocities at kpc scales ~ 60% lower than their CDM counterparts, but are reduced by only 10% in 4 keV cosmologies. Since relic masses < 2-3 keV are ruled out by constraints from the number of Milky Way satellites and Lyman-\alpha\ forest, WDM has a minor effect in reducing the densities of massive satellites. Given the uncertainties on the mass and formation epoch of the Milky Way, the need for reducing the satellite densities with baryonic effects or WDM is alleviated.

In this third paper of a series, we study the growth and luminosity of black holes (BHs) in motion with respect to their surrounding medium. We run a large set of two-dimensional (2D) axis-symmetric simulations to explore a large parameter space of initial conditions and formulate an analytical model for the accretion. Contrary to the case without radiation feedback, the accretion rate increases with increasing BH velocity $\vbh$ reaching a maximum value at $\vbh=2c_{\rm s,in}\sim 50$ km/s, where $c_{\rm s,in}$ is the sound speed inside the "cometary-shaped" \hii region around the BH, before decreasing as $\vbh^{-3}$ when the I-front becomes {\it R}-type (rarefied) and the accretion rate approaches the classical Bondi-Hoyle-Lyttleton solution. The increase of the accretion rate with $\vbh$ is produced by the formation of a {\it D}-type (dense) ionization front (I-front) preceded by a standing bow-shock that reduces the downstream gas velocity to transonic values. There is a range of densities and velocities where the dense shell is unstable producing periodic accretion rate peaks which can significantly increase the detectability of intermediate mass BHs. We find that the mean accretion rate for a moving BH is larger than that of a stationary BH of the same mass if the medium temperature is $T_\infty<10^4$ K. This result could be important for the growth of seed BHs in the multi-phase medium of the first galaxies and for building and early X-ray background that may affect the formation of the first galaxies and the reionization process.

The ages of globular clusters in our own Milky Way are known with precision of about \pm 1 Gyr, hence their formation at redshifts z>~3 and their role in hierarchical cosmology and the reionization of the intergalactic medium remain relatively undetermined. Here we analyze the effect of globular cluster formation on the observed rest-frame UV luminosity functions (LFs) and UV continuum slopes of high redshift galaxies in the Hubble Ultra Deep Fields. We find that the majority of present day globular clusters have formed during two distinct epochs: at redshifts z ~ 2-3 and at redshifts z>~6. The birth of proto-GC systems produce the steep, faint-end slopes of the galaxy LFs, and because the brightness of proto-GCs fades 5 Myrs after their formation, their blue colors are in excellent agreement with observations. Our results suggest that: i) the bulk of the old globular cluster population with estimated ages >~12 Gyr (about 50% of the total population), formed in the relatively massive dwarf galaxies at redshifts z>~6; ii) proto-GC formation was an important mode of star formation in those dwarf galaxies, and likely dominated the reionization process. Another consequence of this scenario is that some of the most massive Milky Way satellites may be faint and yet undiscovered because tidal stripping of a dominant GC population precedes significant stripping of the dark matter halos of these satellites. This scenario may alleviate some remaining tensions between CDM simulations and observations.

Spatial Fluctuations in the Cosmic Infrared Background have now been measured out to sub-degree scales showing a strong clustering signal from unresolved sources. We attempt to explain these measurement by considering faint galaxy populations at z<6 as the underlying sources for this signal using 233 measured UV, optical and NIR luminosity functions (LF) from a variety of surveys covering a wide range of redshifts. We populate the lightcone and calculate the total emission redshifted into the near-IR bands in the observer frame and recover the observed optical and near-IR galaxy counts to a good accuracy. Using a halo model for the clustering of galaxies with an underlying LCDM density field, we find that fluctuations from known galaxy populations are unable to account for the large scale CIB clustering signal seen by HST/NICMOS, Spitzer/IRAC and AKARI/IRC and continue to diverge out to larger angular scales. Our purely empirical reconstruction shows that known galaxy populations are not responsible for the bulk of the fluctuation signal seen in the measurements and suggests an unknown population of very faint and highly clustered sources dominating the signal.

We model fluctuations in the Cosmic Infrared Background (CIB) arising from known galaxy populations using 233 measured UV, optical and NIR luminosity functions (LF) from a variety of surveys spanning a wide range of redshifts. We compare best-fit Schechter parameters across the literature and find clear indication of evolution with redshift. Providing fitting formulae for the multi-band evolution of the LFs out to z~5, we calculate the total emission redshifted into the near-IR bands in the observer frame and recover the observed optical and near-IR galaxy counts to a good accuracy. Our empirical approach, in conjunction with a halo model describing the clustering of galaxies, allows us to compute the fluctuations of the unresolved CIB and compare the models to current measurements. We find that fluctuations from known galaxy populations are unable to account for the large scale CIB clustering signal seen by Spitzer/IRAC and AKARI/IRC and continue to diverge out to larger angular scales. This holds true even if the LFs are extrapolated out to faint magnitudes with a steep faint-end slope all the way to z=8. We also show that removing resolved sources to progressively fainter magnitude limits, isolates CIB fluctuations to increasingly higher redshifts. Our empirical approach suggests that known galaxy populations are not responsible for the bulk of the fluctuation signal seen in the measurements and favors a very faint population of highly clustered sources.

This paper, the second of a series on radiation-regulated accretion onto black holes (BHs) from galactic scales, focuses on the effects of radiation pressure and angular momentum of the accreting gas. We simulate accretion onto intermediate-mass black holes, but we derive general scaling relationships that are solutions of the Bondi problem with radiation feedback valid for any mass of the BH $M_{\rm bh}$. Thermal pressure of the ionized sphere around the BH regulates the accretion rate producing periodic and short-lived luminosity bursts. We find that for ambient gas densities exceeding $n^{\rm cr}_{\rm H,\infty} \propto M_{\rm bh}^{-1}$, the period of the oscillations decreases rapidly and the duty cycle increases from 6%, in agreement with observations of the fraction of active galactic nuclei at $z\sim 3$, to 50%. The mean accretion rate becomes Eddington limited for $n_{\rm H,\infty}>n_{\rm H,\infty}^{\rm Edd} \simeq n_{\rm H,\infty}^{\rm cr} T_{\infty,4}^{-1}$ where $T_{\infty,4}$ is the gas temperature in units of $10^4$ K. In the sub-Eddington regime, the mean accretion rate onto BHs is about $1% T_{\infty,4}^{2.5}$ of the Bondi rate, thus is proportional to the thermal pressure of the ambient medium. The period of the oscillations coincides with depletion timescale of the gas inside the ionized bubble surrounding the BH. Gas depletion is dominated by a pressure gradient pushing the gas outward if $n_{\rm H,\infty}<n^{\rm cr}_{\rm H,\infty}$ and by accretion onto the BH otherwise. Generally, for $n_{\rm H,\infty}<n^{\rm cr}_{\rm H,\infty}$ angular momentum does not affect significantly the accretion rate and period of the oscillations.

We use a new set of cold dark matter simulations of the local universe to investigate the distribution of fossils of primordial dwarf galaxies within, and around the Milky Way. Throughout, we build upon previous results showing agreement between the observed stellar properties of a subset of the ultra-faint dwarfs and our simulated fossils. Here, we show that fossils of the first galaxies have galactocentric distributions and cumulative luminosity functions consistent with observations. In our model there are ~ 300 luminous satellites orbiting the Milky Way, ~50-70% of which are well preserved fossils, with this fraction decreasing with galactocentric distance. Within the Milky Way virial radius, the majority of these fossils have luminosities L_V<10^5 L_solar. This work produces an overabundance of bright dwarf satellites (L_V > 10^4 L_solar) with respect to observations where observations are nearly complete. The "bright satellite problem" is most evident in the outer parts of the Milky Way. We estimate that, although relatively bright, the primordial stellar populations are very diffuse, producing a population with surface brightnesses below surveys` detection limits and are easily stripped by tidal forces. Although we cannot yet present unmistakable evidence for the existence of the fossils of first galaxies in the Local Group, the results of our studies suggest observational strategies that may demonstrate their existence. Primarily, the detection of "ghost halos" of primordial stars around isolated dwarfs would prove that stars formed in minihalos (M<10^8 M_solar) before reionization, and strongly suggest that at least a fraction of the ultra-faint dwarfs are fossils of the first galaxies.

We present a new method for generating initial conditions for LCDM N-body simulations which provides the dynamical range necessary to follow the evolution and distribution of the fossils of the first galaxies on Local Volume, 5-10 Mpc, scales. The initial distribution of particles represents the position, velocity and mass distribution of the dark and luminous halos extracted from pre-reionization simulations. We confirm previous results that ultra-faint dwarfs have properties compatible with being well preserved fossils of the first galaxies. However, because the brightest pre-reionization dwarfs form preferentially in biased regions, they most likely merge into non-fossil halos with circular velocities >20-30 km/s. Hence, we find that the maximum luminosity of true-fossils in the Milky Way is L_V<10^5 L_solar, casting doubts on the interpretation that some classical dSphs are true-fossils. In addition, we argue that most ultra-faints at small galactocentric distance, R<50 kpc, had their stellar properties modified by tides, while a large population of fossils is still undetected due to their extremely low surface brightness log(Sigma_V) < -1.4. We estimate that the region outside R_50 (~ 400 kpc) up to 1 Mpc from the Milky Way contains about a hundred true fossils of the first galaxies with V-band luminosities 10^3 - 10^5 L_solar and half-light radii, r_hl ~ 100-1000 pc.

We present a comprehensive analysis of the halo model of cosmological large to small-scale structure statistics in the case of warm dark matter (WDM) structure formation scenarios. We include the effects of WDM on the linear matter power spectrum, halo density profile, halo concentration relation, halo mass function, subhalo density profile, subhalo mass function and biasing of the smooth dark matter component. As expected, we find large differences at the smallest physical scales in the nonlinear matter power spectrum predicted in the halo model between WDM and cold dark matter even for reasonably high-scale WDM particle masses. We find that significant effects are contributed from the alteration of the halo density profile and concentration, as well as the halo mass function. We further find that the effects of WDM on the subhalo population are important but sub-dominant. Clustering effects of the biasing of the smooth component in WDM is not largely significant.

We study the effect of radiative feedback on accretion onto intermediate mass black holes (IMBHs) using the hydrodynamical code ZEUS-MP with a radiative transfer algorithm. In this paper, the first of a series, we assume accretion from a uniformly dense gas with zero angular momentum and extremely low metallicity. Our 1D and 2D simulations explore how X-ray and UV radiation emitted near the black hole regulates the gas supply from large scales. Both 1D and 2D simulations show similar accretion rate and period between peaks in accretion, meaning that the hydro-instabilities that develop in 2D simulations do not affect the mean flow properties. We present a suite of simulations exploring accretion across a large parameter space, including different radiative efficiencies and radiation spectra, black hole masses, density and temperature, $T_\infty$, of the neighboring gas. In agreement with previous studies we find regular oscillatory behavior of the accretion rate, with duty cycle $\sim 6%$, mean accretion rate 3% $(T_{\infty}/10^4 {\rm K})^{2.5}$ of the Bondi rate and peak accretion $\sim 10$ times the mean for $T_{\infty}$ ranging between 3000K and 15000K. We derive parametric formulas for the period between bursts, the mean accretion rate and the peak luminosity of the bursts and thus provide a formulation of how feedback regulated accretion operates. The temperature profile of the hot ionized gas is crucial in determining the accretion rate, while the period of the bursts is proportional to the mean size of the Str\"{o}mgren sphere and we find qualitatively different modes of accretion in the high vs. low density regimes. We also find that softer spectrum of radiation produces higher mean accretion rate.

We have conducted N-body simulations of the growth of Milky Way-sized halos in cold and warm dark matter cosmologies. The number of dark matter satellites in our simulated Milky Ways decreases with decreasing mass of the dark matter particle. Assuming that the number of dark matter satellites exceeds or equals the number of observed satellites of the Milky Way we derive lower limits on the dark matter particle mass. We find with 95% confidence m_s > 13.3 keV for a sterile neutrino produced by the Dodelson and Widrow mechanism, m_s > 8.9 keV for the Shi and Fuller mechanism, m_s > 3.0 keV for the Higgs decay mechanism, and m_{WDM} > 2.3 keV for a thermal dark matter particle. The recent discovery of many new dark matter dominated satellites of the Milky Way in the Sloan Digital Sky Survey allows us to set lower limits comparable to constraints from the complementary methods of Lyman-alpha forest modeling and X-ray observations of the unresolved cosmic X-ray background and of dark matter halos from dwarf galaxy to cluster scales. Future surveys like LSST, DES, PanSTARRS, and SkyMapper have the potential to discover many more satellites and further improve constraints on the dark matter particle mass.

In cold dark matter cosmologies, small mass halos outnumber larger mass halos at any redshift. However, the lower bound for the mass of a galaxy is unknown, as are the typical luminosity of the smallest galaxies and their numbers in the universe. The answers depend on the extent to which star formation in the first population of small mass halos may be suppressed by radiative feedback loops operating over cosmological distance scales. If early populations of dwarf galaxies did form in significant number, their relics should be found today in the Local Group. These relics have been termed "fossils of the first galaxies". This paper is a review that summarizes our ongoing efforts to simulate and identify these fossils around the Milky Way and Andromeda. It is widely believed that reionization of the intergalactic medium would have stopped star formation in the fossils of the first galaxies. Thus, they should be among the oldest objects in the Universe. However, here we dispute this idea and discuss a physical mechanism whereby relatively recent episodes of gas accretion and star formation would be produced in some fossils of the first galaxies. We argue that fossils may be characterized either by a single old population of stars or by a bimodal star formation history. We also propose that the same mechanism could turn small mass dark halos formed before reionization into gas-rich but starless "dark galaxies". We believe that current observational data support the thesis that a fraction of the new ultra-faint dwarfs recently discovered in the Local Group are in fact fossils of the first galaxies.

We propose the existence of ultracompact minihalos as a new type of massive compact halo object (MACHO) and suggest an observational test to discover them. These new MACHOs are a powerful probe into the nature of dark matter and physics in the high energy Universe. Non-Gaussian energy-density fluctuations produced at phase transitions (e.g., QCD) or by features in the inflaton potential can trigger primordial black hole (PBH) formation if their amplitudes are delta > 30%. We show that a PBH accumulates over time a sufficiently massive and compact minihalo to be able to modify or dominate its microlensing magnification light curve. Perturbations of amplitude 0.03% < delta < 30% are too small to form PBHs, but can nonetheless seed the growth of ultracompact minihalos. Thus, the likelihood of ultracompact minihalos as MACHOs is greater than that of PBHs. In addition, depending on their mass, they may be sites of formation of the first PopIII stars. Ultracompact minihalos and PBHs produce a microlensing light curve that can be distinguished from that of a "point-like" object if high-quality photometric data are taken for a sufficiently long time after the peak of the magnification event. This enables them to be detected below the stellar-lensing "background" toward both the Magellanic Clouds and the Galactic bulge.

In this letter we revisit the idea of reionization feedback on dwarf galaxy formation. We show that primordial minihalos with v_cir<20 km/s stop accreting gas after reionization, as it is usually assumed, but in virtue of their increasing concentration and the decreasing temperature of the intergalactic medium as redshift decreases below z=3, they have a late phase of gas accretion and possibly star formation. We expect that pre-reionization fossils that evolved on the outskirts of the Milky Way or in isolation show a bimodal star formation history with 12 Gyr old and <10 Gyr old population of stars. Leo T fits with this scenario. Another prediction of the model is the possible existence of a population of gas rich minihalos that never formed stars. More work is needed to understand whether a subset of compact high-velocity clouds can be identified as such objects or whether an undiscovered population exists in the voids between galaxies.

Probing the growth of structure from the epoch of hydrogen recombination to the formation of the first stars and galaxies is one of the most important uncharted areas of observational cosmology. Far-IR spectroscopy covering $\lambda$ 100-500 microns from space, and narrow partial transmission atmospheric bands available from the ground, opens up the possibility of probing the molecular hydrogen and metal fine-structure lines from primordial clouds from which the first stars and galaxies formed at 6 < z $<$ 15. Building on Spitzer observations of unexpectedly powerful H2 emission from shocks, we argue that next-generation far-IR space telescopes may open a new window into the main cloud cooling processes and feedback effects which characterized this vital, but unexplored epoch. Without this window, we are essential blind to the dominant cloud cooling which inevitably led to star formation and cosmic reionization.

We argue that, at least a fraction of the newly discovered population of ultra-faint dwarf spheroidal galaxies in the Local Group constitute the fossil relic of a once ubiquitous population of dwarf galaxies formed before reionization with circular velocities smaller than $v_{c}^{cr} \sim 20$ km/s. We present several arguments in support of this model. The number of luminous Milky Way satellites inferred from observations is larger than the estimated number of dark halos in the Galaxy that have, or had in the past, circular velocity $>v_{c}^{cr}$, as predicted by the "Via Lactea" simulation. This implies that some ultra-faint dwarfs are fossils. However, this argument is weakened by recent results from the "Aquarius" simulations showing that the number of Galactic dark matter satellites is 2.5 larger than previously believed. Secondly, the existence of a population of ultra-faint dwarfs was predicted by cosmological simulations in which star formation in the first minihalos is reduced -- but not suppressed -- by radiative feedback. Here, we show the statistical properties of the fossil galaxies in those simulations are consistent with observations of the new dwarf population and with the number and radial distribution of Milky Way satellites as a function of their luminosity. Finally, the observed Galactocentric distribution of dwarfs is consistent with a fraction of dSphs being fossils. To make our case more compelling, future work should determine whether stellar chemical abundances of simulated "fossils" can reproduce observations and whether the tidal scenarios for the formation of Local Group dwarf spheroidals are equally consistent with all available observations.

The fact that dark matter (DM), thus far, has revealed itself only on scales of galaxies and larger, again thrusts onto astrophysics the opportunity and the responsibility to confront the age old mystery "What is the nature of matter?" By deriving basic data on the nature of DM - e.g., mass of its particle(s), present mean temperature, distribution in galaxies and other structures in the universe, and capacity for dissipational collapse - we will be uncovering the properties of the dominant species of matter in the universe and significantly extending the standard models of particle physics. Determining the mass of the DM particle to an order of magnitude would help to sort out the particle family to which it (or they) belongs. Beyond mass, there are issues of stability. The DM particle may be unstable with a measurable half-life, or it may become unstable after absorbing a certain amount of energy from collisions. In both cases it would contribute to the present hot dark matter component. Some key parameters of DM can most accurately be measured in the very nearby universe because DM dominates the mass in the outer Milky Way (MW), in other galaxies in the Local Group, and in the Local Group in its entirety. The presence and amount of DM can be quantified by study of dynamical processes observable in fine detail within these entities. Precise measurements of 3-D velocities for stars, coherent star streams, and stars in satellite stellar systems out to the edge of the Galaxy can reveal "what is the shape, orientation, density law, and lumpiness of the dark matter halo" as well as "what is the total mass of the Galaxy?"

In two previous papers, we presented simulations of the first galaxies in a representative volume of the Universe. The simulations are unique because we model feedback-regulated galaxy formation, using time-dependent, spatially-inhomogeneous radiative transfer coupled to hydrodynamics. Here, we study the properties of simulated primordial dwarf galaxies with masses <2x10^8 Msolar and investigate their impact on the intergalactic medium. While many primordial galaxies are dark, about 100--500 per comoving Mpc^3 are luminous but relatively faint. They form preferentially in chain structures, and have low surface brightness stellar spheroids extending to 20% of the virial radius. Their interstellar medium has mean density n_H~10--100 cm^-3, metallicity Z~ 0.01--0.1 Zsolar and can sustain a multi-phase structure. With large scatter, the mean efficiency of star formation scales with halo mass, <f_*> \propto M_dm^2, independent of redshift. Because of feedback, halos smaller than a critical mass, M_crit(z), are devoid of most of their baryons. More interestingly, we find that dark halos have always a smaller M_crit(z) than luminous ones. Metal enrichment of the intergalactic medium is inhomogeneous, with only a 1%--10% volume filling factor of enriched gas with [Z/H]>-3.0 and 10%--50% with [Z/H]>-5.0. At z=10, the fraction of stars with metallicity Z<10^-3 Zsolar is 10^-6 of the total stellar mass. Although detections of high-redshift dwarf galaxies with the James Webb Space Telescope will be a challenge, studies of their fossil records in the local Universe are promising because of their large spatial density.

We investigate the effect of non-evaporating primordial black holes (PBHs) on the ionization and thermal history of the universe. X-rays emitted by gas accretion onto PBHs modify the cosmic recombination history, producing measurable effects on the spectrum and anisotropies of the Cosmic Microwave Background (CMB). Using the third-year WMAP data and FIRAS data we improve existing upper limits on the abundance of PBHs with masses >0.1 Msun by several orders of magnitude. Fitting WMAP3 data with cosmological models that do not allow for non-standard recombination histories, as produced by PBHs or other early energy sources, may lead to an underestimate of the best-fit values of the amplitude of linear density fluctuations (sigma_8) and the scalar spectral index (n_s). Cosmological parameter estimates are affected because models with PBHs allow for larger values of the Thomson scattering optical depth, whose correlation with other parameters may not be correctly taken into account when PBHs are ignored. Values of tau_e=0.2, n_s=1 and sigma_8=0.9 are allowed at 95% CF. This result that may relieve recent tension between WMAP3 data and clusters data on the value of sigma_8. PBHs may increase the primordial molecular hydrogen abundance by up to two orders of magnitude, this promoting cooling and star formation. The suppression of galaxy formation due to X-ray heating is negligible for models consistent with the CMB data. Thus, the formation rate of the first galaxies and stars would be enhanced by a population of PBHs.

Several recent studies suggest a correlation between dark matter halo mass and the shape of the density profile. We re-analyze simulations from Ricotti (2003) in which such a correlation was proposed. We use a standard analysis of the halo density profiles and compare the old simulations to new ones performed with Gadget2, including higher resolution runs. We confirm Ricotti's result that, at virialization, the central log slopes alpha, at 5%-10% of the virial radius are correlated with the halo mass and that the halo concentration is a universal constant. Our results do not contradict the majority of published papers: when using a split power law to fit the density profiles, due to the alpha-concentration degeneracy, the fits are consistent with halos having a universal shape with alpha=1 or 1.5 and concentrations that depend on the mass, in agreement with results published elsewhere. Recently, several groups have found no evidence for convergence of the inner halo profile to a constant power law. The choice of a split power law parameterization used in this letter is motivated by the need to compare our results to previous ones and is formally valid because we are not able to resolve regions where the slope of the fitting function reaches its asymptotic constant value. Using a non-parameterized technique, we also show that the density profiles of dwarf galaxies at z ~ 10 have a log slope shallower than 0.5 within 5% of the virial radius.

This paper presents a study of quasi-steady spherical accretion in the early Universe, before the formation of the first stars and galaxies. The main motivation is to derive the basic formulas that will be used in a companion paper to calculate the accretion luminosity of primordial black holes and their effect on the cosmic ionization history. The following cosmological effects are investigated: the coupling of the gas to the CMB photon fluid (i.e., Compton drag), Hubble expansion, and the growth of the dark matter halo seeded by the gravitational potential of the central point mass. The gas equations of motion are solved assuming either a polytropic or an isothermal equation of state. We consider the cases in which the accreting object is a point mass or a spherical dark matter halo with power-law density profile, as predicted by the theory of "secondary infall''. Analytical solutions for the sonic radius and fitting formulas for the accretion rate are provided. Different accretion regimes exist depending on the mass of the accreting object. If the black hole mass is smaller than 50-100 Msun, gas accretion is unaffected by Compton drag. A point mass and an extended dark halo of equal mass accrete at the same rate if M>5000 Msun, while smaller mass dark halos accrete less efficiently than the equivalent point mass. For masses M>3 x 10^4 Msun, the viscous term due to the Hubble expansion becomes important and the assumption of quasi-steady flow fails. Hence, the steady Bondi solutions transition to the time-dependent self-similar solutions for "cold cosmological infall".

We discuss the possibilities for primordial black holes (PBHs) to grow via the accretion of dark matter. In agreement with previous works, we find that accretion during the radiation-dominated era does not lead to a significant mass increase. However, during matter-domination, PBHs may grow by up to two orders of magnitude in mass through the acquisition of large dark matter halos. We discuss the possibility of PBHs being an important component in dark matter halos of galaxies as well as their potential to explain the ultra-luminous x-ray sources (ULXs) observed in nearby galactic disks. We point out that although PBHs are ruled out as the dominant component of dark matter, there is still a great deal of parameter space open to them playing a role in the modern-day universe. For example, a primordial halo population of PBHs each at $10^{2.5} M_\odot$ making up 0.1% of the dark matter grow to $10^{4.5} M_\odot$ via the accumulation of dark matter halos to account for $\sim 10%$ of the dark matter mass by a redshift of $z \approx 30$. These intermediate mass black holes may then ``light up'' when passing through molecular clouds, becoming visible as ULXs at the present day, or they may form the seeds for supermassive black holes at the centers of galaxies.

A model for gas outflows is proposed which simultaneously explains the correlations between the (i) equivalent width of low ionization and Ly-alpha lines, (ii) outflow velocity, and (iii) star formation rate observed in Lyman Break Galaxies (LBGs). Our interpretation implies that LBGs host short-lived (30 +/- 5 Myr) starburst episodes observed at different evolutionary phases. Initially, the starburst powers a hot wind bound by a denser cold shell, which after about 5 Myr becomes dynamically unstable and fragments; afterwards the fragment evolution is approximately ballistic while the hot bubble continues to expand. As the fragments are gravitationally decelerated, their screening ability of the starlight decreases as the UV starburst luminosity progressively dims. LBG observations sample all these evolutionary phases. Finally, the fragments fall back onto the galaxy after approximately 60 Myr. This phase cannot be easily probed as it occurs when the starburst UV luminosity has already largely faded; however, galaxies dimmer in the UV than LBGs should show infalling gas.

We compare the properties of dwarf galaxies in the Local Group with the simulated galaxies formed before reionization in a cosmological simulation of unprecedented spatial and mass resolution,including radiative feedback effects. We find that a subset of the Local Group dwarfs are remarkably similar to the simulated dwarf galaxies in all their properties already before reionization. Based on this similarity, we propose the hypothesis that Local Group dwarfs form in a variety of ways: some of them are ``true fossils'' of the pre-reionization era, some of them form most of their stars later, after reionization (we call them ``survivors'' of the reionization era), and the rest of them form an intermediate group of ``polluted fossils''. We also identify a simple observational test that is able to falsify our hypothesis.

We use cosmological simulations to study the X-ray ionisation and heating of the intergalactic medium by an early population of accreting black holes. By considering observational constraints from the X-ray background, we find an upper limit for the optical depth to Thompson scattering tau_e~0.17. The redshifted soft X-ray background from these early sources produces:(i) fully ionised atomic hydrogen in the low density intergalactic medium before redshift z~7 (consequently stellar reionisation is characterised by an instantaneous overlap phase of HII regions),(ii) a second HeII reionisation at z~3 and (iii) heats the intergalactic medium to near 10000 K at low redshifts. The typical luminosity in the soft X-ray band of the galaxies hosting the black holes is about one order of magnitude below the sensitivity limit of the Chandra deep field. About a thousand of these sources may be present per square arcmin of the sky, producing detectable fluctuations. Few rarer objects could be luminous enough to be visible in the Chandra deep field. XEUS and Con-X satellites will be able to detect more of these sources that, if radio loud, could be used to study the 21cm forest in absorption.A signature of an early X-ray preionisation is the production of secondary CMB anisotropies on small angular scales. We find that in these models the power spectrum of temperature fluctuations increases with decreasing angular scale (dT~16 muK at 1arcsec scales), while for stellar reionisation scenarios the power decreases on smaller scales. We also show that the redshifted 21 cm radiation from neutral hydrogen can be marginally detected in emission at redshifts 7<z<12. At a redshift of about z~30 a stronger and narrower (in redshift space) signal in absorption against the CMB, that is peculiar to these models,could be detectable.[abridged]

We investigate the possibility that there is a first phase of partial ionisation due to X-rays produced by black hole accretion in small-mass galaxies at redshifts 7<z<20. This is followed by complete reionisation by stellar sources at z~7. This scenario is motivated by the large optical depth to Thompson scattering, tau_e=0.17, measured by WMAP. But it is also consistent with the observed Gunn-Peterson trough in the spectra of quasars at z~5-6. We use a semianalytic code to explore models with different black hole accretion histories and cosmological parameters.We find that ``preionisation'' by X-rays can increase the intergalactic medium (IGM) optical depth from tau_e~0.06 given by stellar sources only, to 0.1<tau_e<0.2, if a fraction of baryons 5x10^-5 is accreted onto seed black holes produced in the collapse of low metallicity, high mass stars before z~15. To be effective, preionisation requires a non-negligible star formation in the first small-mass galaxies in hich seed black holes are formed.The model predicts that dwarf spheroidal galaxies may host a mass in black holes that is 5-40% of their stellar mass. The redshifted X-ray background produced by this early epoch of black hole accretion constitutes about 5-10% of the X-ray background in the 2-50 keV bands and roughtly half of the currently estimated black hole mass density was formed at early times.Moreover, in most models, the photons from the redshifted background are sufficient to fully reionise HeII at redshift z~3 without any additional contribution from quasars at lower redshifts and the temperature of the mean density intergalactic medium remains close to 10^4 K down to redshift z~1.[abridged]

In a recent paper Pell'o et al. have reported a candidate z=10 galaxy, A1835#1916, which was found in a near-infrared survey of the central regions of the gravitational lensing cluster A1835. If this detection is confirmed and the detection rate turns out to be typical, then the volume averaged ultraviolet emissivity must be rising rapidly with increasing redshift. For a magnification due to gravitational lensing by a factor M>25 estimated by Pell'o et al., the inferred star formation rate density at z=10 would be about one order of magnitude higher than estimates of the star formation rate density at z=6. Objects at z=10 would contribute substantially to the total source counts at 1.6 micronm and the estimated space density of sources may exceed the space density of dark matter haloes in a LCDM model. We therefore argue that if A1835#1916 is indeed at z=10 then either the magnification factor may have been overestimated or the galaxy has a top-heavy initial mass function. Sources with the ultraviolet flux and space density of A1835#1916 may be sufficient to reionise most of the diffuse hydrogen in the Universe at z=10. We further use a correlation between the equivalent width and the redshift of the Ly-alpha emission line with respect to the systemic redshift observed in Lyman break galaxies to obtain constraints on the ionisation state of the surrounding intergalactic medium (IGM) from the Gunn-Peterson absorption. These constraints also argue in favour of the surrounding IGM being partially or fully ionised.

In this paper, we study the dynamical stability and time evolution of the central dark matter cores of low-mass (about 10^8-10^9 M_solar) galactic haloes found in recent cold dark matter simulations at high redshift. From those simulations we extract three haloes, assembled by hierarchical merging, that at redshift z>10 display a core and we evolve them without further merging to low redshift using direct N-body integration. The central core in the dark matter profile is found to be dynamically stable: it survives for many crossing times without evolution into a cusp. This result supports the claim that the mass dependence of the central dark matter profiles of simulated haloes is a direct consequence of the power spectrum of primordial density fluctuations. In addition, we show that the simulated dark matter profiles, if they evolved in isolation, are consistent with the observed velocity dispersion profile of stars in the inner parts of the Draco dwarf spheroidal galaxy. Simple scaling arguments are reviewed which explain the evolution of the concentration parameter with redshift. We also review some arguments used to derive the logarithmic slope of the inner and outer density profile.

In this talk I discuss the role of proto-globular clusters as the dominant sources of radiation that reionised hydrogen in the intergalactic medium (IGM) at redshift z~6. Observations at lower redshifts indicate that only a small fraction, f_esc, of hydrogen ionising radiation emitted from massive stars can escape unabsorbed by the galaxy into the IGM. High redshift galaxies are expected to be more compact and gas rich than present day galaxies, consequently f_esc from their disks or spheroids might have been very small. But if the sites of star formation in the galaxies are off-centre and if the star formation efficiency of the proto-clusters is high, then the mean f_esc calculated for these objects only, is expected to be close to unity. Here I argue that this mode of star formation is consistent with several models for globular clusters formation. Using simple arguments based on the observed number of globular cluster systems in the local universe and assuming that the oldest globular clusters formed before reionisation and had f_esc~1, I show that they produced enough ionising photons to reionise the IGM at z~6. I also enphasize that globular cluster formation might have been the dominant mode of star formation at redshifts from 6 to 12.

We investigate the effects of a top-heavy stellar initial mass function on the reionisation history of the intergalactic medium (IGM). We use cosmological simulations that include self-consistently the feedback from ionising radiation, H_2 dissociating radiation and supernova (SN) explosions. We find that it is difficult to reionise the IGM at z_rei>10 with stellar sources even after making extreme assumptions. If star formation in 10^9 M_\odot galaxies is not suppressed by SN explosions, the optical depth to Thomson scattering is tau_e< 0.13. If we allow for the normal energy input from SNe or if pair-instability SNe are dominant, we find taue<0.09. Assuming normal yields for the first stars (popIII), the mean metallicity of the IGM is already Z/Z_\odot=2x10^-3 (10^-3<Z/Z_\odot<1 in overdense regions) when the IGM mean ionisation fraction is less than 10%. For these reasons popIII stars cannot contribute significantly to reionisation unless the mechanical energy input from SNe is greatly reduced and either the metal yield or the mixing efficiency is reduced by a factor of 1000. Both problems have a solution if popIII stars collapse to black holes. This can happen if, having masses M_*< 130 M_\odot, they are characterised by heavy element fall-back or if, having masses M_*>260 M_\odot, they collapse directly onto black holes without exploding as SNe.If metal-poor stars are initially important and collapse to black holes is the typical outcome, then the secondary emission of ionising radiation from accretion on SN induced seed black holes, might be more important than the primary emission. Surprisingly, including feedback effects, we estimate that a warm dark matter scenario (with particle mass of 1.25 keV) reduces tau_e by only approximately 10%.[abridged]

I compare the density profile of dark matter (DM) halos in cold dark matter (CDM) N-body simulations with 1 Mpc, 32 Mpc, 256 Mpc and 1024 Mpc box sizes. In dimensionless units the simulations differ only for the initial power spectrum of density perturbations. I compare the profiles when the most massive halos are composed of about 10^5 DM particles. The DM density profiles of the halos in the 1 Mpc box show systematically shallower cores with respect to the corresponding halos in the 32 Mpc simulation that have masses, M_{dm}, typical of the Milky Way and are fitted by a NFW profile. The DM density profiles of the halos in the 256 Mpc box are consistent with having steeper cores than the corresponding halos in the 32 Mpc simulation, but higher mass resolution simulations are needed to strengthen this result. Combined, these results indicate that the density profile of DM halos is not universal, presenting shallower cores in dwarf galaxies and steeper cores in clusters. Physically the result sustains the hypothesis that the mass function of the accreting satellites determines the inner slope of the DM profile. In comoving coordinates, r, the profile \rho_{dm} \propto 1/(X^\alpha(1+X)^{3-\alpha}), with X=c_\Delta r/r_\Delta, r_\Delta is the virial radius and \alpha =\alpha(M_{dm}), provides a good fit to all the DM halos from dwarf galaxies to clusters at any redshift with the same concentration parameter c_\Delta ~ 7. The slope, \gamma, of the outer parts of the halo appears to depend on the acceleration of the universe: when the scale parameter is a=(1+z)^{-1} < 1, the slope is \gamma ~ 3 as in the NFW profile, but \gamma ~ 4 at a > 1 when \Omega_\Lambda ~ 1 and the universe is inflating.[abridged]

I analyze the properties of the first galaxies in cosmological simulations with radiative feedback. Preliminary results indicate similarities with the observed properties of the bulk of dwarf spheroidal galaxies (dSphs) in the Local Group and Andromeda. I briefly discuss observational tests that could help in understanding the impact of a population of small primordial objects on the cosmic evolution.

A problem still unsolved in cosmology is the identification of the sources of radiation able to reionize HI in the intergalactic medium (IGM) by z~6. Theoretical works and observations seem to indicate that the fraction, f_esc, of HI ionizing radiation emitted from galaxies that escapes into the IGM is small in the local universe (f_esc < 10%). At high redshift galaxies are more compact and probably gas rich implying smaller values of f_esc from their disks or spheroids. But if the sites of star formation are displaced from the disk or spheroid and the star formation efficiency of the proto-clusters is high, then f_esc should be about one. This star formation scenario is consistent with several models for globular clusters formation. Using simple arguments based on the observed number of globular cluster systems in the local universe and assuming that the oldest globular clusters formed before reionization and had f_esc~1, I show that they produced enough ionizing photons to reionize the IGM at z=6.

[abridged]We use 3D cosmological simulations with radiative transfer to study the formation and evolution of the first galaxies in a LCDM cosmology. We find that the first luminous objects ("small-halos") are characterized by a bursting star formation (SF) that is self-regulated by a feedback process acting on cosmological instead of galactic scales. The global star formation history is regulated by the mean number of ionizing photons that escape from each source, \epsilon_{UV}\fesc. It is almost independent of the assumed star formation efficiency parameter, \epsilon_*, and the intensity of the dissociating background. The main feedback process that regulates the SF is the re-formation of H_2 in front of HII regions and inside relic HII regions. The HII regions remain confined inside filaments, maximizing the production of H_2 in overdense regions through cyclic destruction/reformation of H_2. If \epsilon_{UV}\fesc > 10^{-8}/\epsilon_* the SF is self-regulated, photo-evaporation of "small-halo" objects dominate the metal pollution of the low density IGM, and the mass of produced metals depends only on \fesc. If \epsilon_{UV}\fesc \simlt 10^{-8}/\epsilon_*, positive feedback dominates, and "small-halo" objects constitute the bulk of the mass in stars and metals at least until redshift z \sim 10. "Small-halo" objects cannot reionize the universe because the feedback mechanism confines the HII regions inside the large scale structure filaments. In contrast to massive objects, which can reionize voids, "small-halo" objects partially ionize only the dense filaments while leaving the voids mostly neutral.

In cold dark matter (CDM) cosmogonies, low-mass objects play an important role in the evolution of the universe. Not only are they the first luminous objects to shed light in a previously dark universe, but, if their formation is not inhibited by their own feedback, they dominate the galaxy mass function until redshift z \sim 5. In this paper we present and discuss the implementation of a 3D cosmological code that includes most of the needed physics to simulate the formation and evolution of the first galaxies with a self-consistent treatment of radiative feedback. The simulation includes continuum radiative transfer using the ``Optically Thin Variable Eddington Tensor'' (OTVET) approximation and line-radiative transfer in the H_2 Lyman-Werner bands of the background radiation. We include detailed chemistry for H_2 formation/destruction, molecular and atomic cooling/heating processes, ionization by secondary electrons, and heating by Ly\alpha resonant scattering. We find that the first galaxies ("small-halos") are characterized by a bursting star formation, self-regulated by a feedback process that acts on cosmological scales. Their formation is not suppressed by feedback processes; therefore, their impact on cosmic evolution cannot be neglected. The main focus of this paper is on the methodology of the simulations, and we only briefly introduce some of the results. An extensive discussion of the results and the nature of the feedback mechanism are the focus of a companion paper.

Kinetic energy stored in ISM bulk/turbulent motions is a crucial ingredient to properly describe most properties of observed galaxies. By using Monte Carlo simulations, we investigate how this energy is injected by supernovae and dissipated via cloud collisions and derive the corresponding ISM velocity probability distribution function (PDF). The functional form of the PDF for the modulus of the velocity dispersion is $$p(v) \propto v^2 \exp[-(v/\sigma)^\beta].$$ The power-law index of the PDF depends only on the value of the average cloud collision elasticity < \epsilon > as \beta = 2\exp(<\epsilon > -1). If \beta and the gas velocity dispersion \sigma are known, the specific kinetic energy dissipated by collisions is found to be \propto \sigma^2 \ln (2 / \beta)/(\beta-0.947); in steady state, this is equal to the energy input from SNe. We predict that in a multiphase, low metallicity (Z \approx 5 \times 10^{-3} Z_\odot) ISM the PDF should be close to a Maxwellian (\beta = 2) with velocity dispersion \sigma \simgt 11$ km s\m; in more metal rich systems (Z \simgt 5 \times 10^{-2} Z_\odot), instead, we expect to observe almost exponential PDFs. This is in good agreement with a number of observations that we review and might explain the different star formation modes seen in dwarfs and spiral galaxies.

We use one-dimensional radiative transfer simulations to study the evolution of H_2 gas-phase (H^- catalyzed) formation and photo-dissociation regions in the primordial universe. We find a new positive feedback mechanism capable of producing shells of H_2 in the intergalactic medium, which are optically thick in some Lyman-Werner bands. While these shells exist, this feedback effect is important in reducing the H_2 dissociating background flux and the size of photo-dissociation spheres around each luminous object. The maximum background opacity of the IGM in the H_2 Lyman-Werner bands is \tau_{H_2} ~ 1-2 for a relic molecular fraction x_{H_2}=2 x 10^{-6}, about 6 times greater than found by Haiman, Abel & Rees. Therefore, the relic molecular hydrogen can decrease the photo-dissociation rate by about an order of magnitude. The problem is relevant to the formation of small primordial galaxies with masses M_{DM} < 10^8 M_\odot, that rely on molecular hydrogen cooling to collapse. Alternatively, the universe may have remained dark for several hundred million years after the birth of the first stars, until galaxies with virial temperature T_{vir} > 10^4 K formed.

Several observational and theoretical arguments suggest that starburst galaxies may rival quasars as sources of metagalactic ionizing radiation at redshifts z > 3. Reionization of the intergalactic medium (IGM) at z > 5 may arise, in part, from the first luminous massive stars. To be an important source of radiative feedback from star formation, a substantial fraction (\sim 1-10%) of the ionizing photons must escape the gas layers of the galaxies. Using models of smoothly distributed gas confined by dark-matter (DM) potentials, we estimate the fraction, f_esc, of Lyc flux that escapes the halos of spherical galaxies as a function of their mass and virialization redshift. The gas density profile is found by solving the equation of hydrostatic equilibrium for the baryonic matter in the potential well of a DM halo with the density profile found by Navarro et. al. (1996). We then perform a parametric study to understand the dependence of fesc on redshift, mass, star-formation efficiency (SFE), stellar density distribution, and OB association luminosity function. We give useful analytical formulas for f_esc(z,M_{DM}). Using the Press-Schechter formalism, we find that stellar reionization at z \sim 7 is probably dominated by small galactic sub-units, with M_{DM} \simlt 10^7 M_\odot and SFE \sim 10 times that in the Milky Way. This may affect the distribution of heavy elements throughout the intergalactic medium.

We develop a method to extract the "effective equation of state" of the intergalactic medium from the doppler b parameter distribution of the low-density Lyman-alpha forest. We test the method on numerical simulations and then apply it to published observations of the Lyman-alpha forest at redshifts z from 0 to 4. We find that the effective equation of state is close to isothermal at redshift z=3, indicating that a second reheating of the IGM took place at z=3. This reheating can plausibly be identified with the reionization of HeII observed to occur at z about 3.

A new class of disk MHD equilibrium solutions is described, which is valid within the standard local (``shearing sheet'') approximation scheme. These solutions have the following remarkable property: velocity streamlines and magnetic lines of force rotate rigidly, even in the presence of differential rotation. This situation comes about because the Lorentz forces acting upon modified epicycles compel fluid elements to follow magnetic lines of force. Field line (and streamline) configurations may be elliptical or hyperbolic, prograde or retrograde. These structures have previously known hydrodynamical analogs: the ``planet'' solutions described by Goodman, Narayan, & Goldreich. The primary focus of this investigation is configurations in the disk plane. A related family of solutions lying in a vertical plane is briefly discussed; other families of solutions may exist. Whether these MHD structures are stable is not yet known, but could readily be determined by three-dimensional simulations. If stable or quasi-stable, these simple structures may find important applications in both accretion and galactic disks.

We present a study of the kinetic energy dissipation in interstellar cloud collisions. The main aim is to understand the dependence of the elasticity (defined as the ratio of the final to the initial kinetic energy of the clouds) on the velocity and mass ratio of the colliding clouds, magnetic field strength, and gas metallicity for head-on collisions. The problem has been studied both analytically and via numerical simulations. We have derived handy analytical relationships that well approximate the analogous numerical results. The main findings of this work are: (i) the kinetic energy dissipation in cloud collisions is minimum (i.e. the collision elasticity is maximum) for a cloud relative velocity $v_r \simeq 30 km s^{-1}$; (ii) the above minimum value is proportional $Z L_c^2$, where $Z$ is the metallicity and $L_c$ is the cloud size: the larger is $Z L_c^2$ the more dissipative (i.e. inelastic) the collision will be; (iii) in general, we find that the energy dissipation decreases when the magnetic field strength, and mass ratio of the clouds are increased and the metallicity is decreased, respectively. We briefly discuss the relevance of this study to the global structure of the interstellar medium and to galaxy formation and evolution.