Allowing for enhanced Ly$\alpha$ photon line emission from Population III dominated stellar systems in the first forming galaxies, we show the 21-cm cosmic dawn signal at $10<z<30$ may substantially differ from standard scenarios. Energy transfer by Ly$\alpha$ photons emerging from galaxies may heat intergalactic gas if HII regions within galaxies are recombination bound, or cool the gas faster than by adiabatic expansion if reddened by winds internal to the haloes. In some cases, differential 21-cm antenna temperatures near $-500$ mK may be achieved at $15<z<25$, similar to the signature detected by the EDGES 21-cm cosmic dawn experiment.

The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory.

We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos in which cooling and condensation of gas with primordial composition is possible at very high redshift ($z > 20$), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies.

Using a 1D Lagrangian code specifically designed to assess the impact of multiple, time-resolved supernovae (SNe) from a single star cluster on the surrounding medium, we test three commonly used feedback recipes: delayed cooling (e.g., used in the GASOLINE-2 code), momentum-energy injection (a resolution-dependent transition between momentum-dominated feedback and energy-dominated feedback used, e.g., in the FIRE-2 code), and simultaneous energy injection (e.g., used in the EAGLE simulations). Our work provides an intermediary test for these recipes: we analyse a setting that is more complex than the simplified scenarios for which many were designed, but one more controlled than a full galactic simulation. In particular, we test how well these models reproduce the enhanced momentum efficiency seen for an 11 SN cluster simulated at high resolution (0.6 pc; a factor of 12 enhancement relative to the isolated SN case) when these subgrid recipes are implemented in low resolution (20 pc) runs. We find that: 1) the delayed cooling model performs well -- resulting in 9 times the momentum efficiency of the fiducial isolated SN value -- when SNe are clustered and $10^{51}$ erg are injected per SN, while clearly over-predicting the momentum efficiency in the single SN test case; 2) the momentum-energy model always achieves good results, with a factor of 5 boost in momentum efficiency; and 3) injecting the energy from all SNe simultaneously does little to prevent over-cooling and greatly under-produces the momentum deposited by clustered SNe, resulting in a factor of 3 decrease in momentum efficiency on the average.

We investigate Lyman-$\alpha$ transmission spikes at $z > 5$ in synthetic quasar spectra, and discuss their connection to the properties of the intergalactic medium and their ability to constrain reionization models. We use state-of-the-art radiation-hydrodynamic simulations from the Cosmic Reionization On Computers series to predict the number of transmission spikes as a function of redshift, both in the ideal case of infinite spectral resolution and in a realistic observational setting. Transmission spikes are produced in highly-ionized underdense regions located in the vicinity of UV sources. We find that most of the predicted spikes are unresolved by current observations, and show that our mock spectra are consistent with observations of the quasar ULAS J1120+0641 in about 15% of the realizations. The spike height correlates with both the gas density and the ionized fraction, but the former link is erased when synthetic spectra are smoothed to realistically achievable spectral resolutions. There exists a linear relationship between spike width and the extent of the associated underdense region, with a slope that is redshift-dependent. In agreement with observations, the spike transmitted flux is suppressed at small distance from bright galaxies as these reside in overdense regions. We argue that this anti-correlation can be used to constrain large-scale density modes.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

While the evolution of superbubbles driven by clustered supernovae has been studied by numerous authors, the resulting radial momentum yield is uncertain by as much as an order of magnitude depending on the computational methods and assumed properties of the surrounding interstellar medium (ISM). In this work, we study the origin of these discrepancies, and seek to determine the correct momentum budget for a homogeneous ISM. We carry out 3D hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of clustered supernova explosions, using a Lagrangian method and checking for convergence with respect to resolution. We find that the terminal momentum of a shell driven by clustered supernovae is dictated primarily by the mixing rate across the contact discontinuity between the hot and cold phases, and that this energy mixing rate is dominated by numerical diffusion even at the highest resolution we can complete, 0.03 $M_\odot$. Magnetic fields also reduce the mixing rate, so that MHD simulations produce higher momentum yields than HD ones at equal resolution. As a result, we obtain only a lower limit on the momentum yield from clustered supernovae. Combining this with our previous 1D results, which provide an upper limit because they allow almost no mixing across the contact discontinuity, we conclude that the momentum yield per supernova from clustered supernovae in a homogeneous ISM is bounded between $2\times 10^5$ and $3\times 10^6$ $M_\odot$ km s$^{-1}$. A converged value for the simple homogeneous ISM remains elusive.

We use high-resolution cosmological hydrodynamical simulations of Milky Way-sized galaxies with varying supernovae feedback strengths and merger histories to investigate the formation of their gaseous halos and especially their hot ($>10^6$~K) X-ray luminous coronae. Our simulations predict the presence of significant hot gas in the halos as early as $z=3-4$, well before the halos ought to be able to sustain hot mode accretion in the conventional picture. The nascent coronae grow inside-out and initially do so primarily as a result of outflows from the central galaxies powered by merger-induced shock heating and strong supernovae feedback, both of which are elemental features of today's successful galaxy formation models. Furthermore, the outflows and the forming coronae also accelerate the transition from cold to hot mode accretion by contributing to the conditions for sustaining stable accretion shocks. They also disrupt the filamentary streams funneling cold gas onto the central galaxies by causing their mouths to fray into a broad delta, detach from the galaxies, and be pushed away to larger radii. And even though at early times the filaments repeatedly re-form, the hot gas and the outflows act to weaken the filaments and accelerate their ultimate disruption. Although galactic outflows are generally thought of as ejective feedback, we find that their action on the filaments suggests a preventive role as well.

We present further results of a survey for absorption line systems in the spectra of four high redshift quasars (5.79 $\le$ z$_{\textrm{em}}$ $\le$ 6.13) obtained with the ESO Very Large Telescope X-Shooter. We identify 36 $\textrm{CIV}$ and 7 $\textrm{SiIV}$ systems with a $\ge$ 5$\sigma$ significance. The highest redshift $\textrm{CIV}$ and $\textrm{SiIV}$ absorbers identified in this work are at z = 5.80738 $\pm$ 0.00017 and z = 5.77495 $\pm$ 0.00038, respectively. We compute the comoving mass density of $\textrm{SiIV}$ ($\Omega_{\textrm{SiIV}}$) and find that it evolves from $\Omega_{\textrm{SiIV}}$ = 4.3$^{+2.1}_{-2.1}$ $\times$10$^{-9}$ at <z> = 5.05 to $\Omega_{\textrm{SiIV}}$ = 1.4$^{+0.6}_{-0.4}$ $\times$10$^{-9}$ at <z> = 5.66. We also measure $\Omega_{\textrm{CIV}}$ = 1.6$^{+0.4}_{-0.1}$ $\times$10$^{-8}$ at <z> = 4.77 and $\Omega_{\textrm{CIV}}$ = 3.4$^{+1.6}_{-1.1}$ $\times$10$^{-9}$ at <z> = 5.66. We classify our $\textrm{CIV}$ absorber population by the presence of associated $\textit{low}$ and/or $\textit{high ionisation}$ systems and compute their velocity width ($\Delta$v$_{90}$). We find that all $\textrm{CIV}$ systems with $\Delta$v$_{90}$ > 200 kms$^{-1}$ have associated $\textit{low ionisation}$ systems. We investigate two such systems, separated by 550 physical kpc along a line of sight, and find it likely that they are both tracing a multi-phase medium where hot and cold gas is mixing at the interface between the CGM and IGM. We further discuss the \textrm{MgII} systems presented in a previous work and we identify 5 $\textrm{SiII}$, 10 $\textrm{AlII}$, 12 $\textrm{FeII}$, 1 $\textrm{CII}$, 7 $\textrm{MgI}$ and 1 $\textrm{CaII}$ associated transitions. We compute the respective comoving mass densities in the redshift range 2 to 6, as allowed by the wavelength coverage.

The tentative detection by the EDGES experiment of a global 21-cm absorption trough centered at redshift 17 opens up the opportunity to study the birth of the first luminous sources, the intensity of radiation backgrounds at cosmic dawn, the thermal and ionization history of the young intergalactic medium. Here, we focus on the astrophysical implications of the Lyman-alpha photon field needed to couple the spin temperature to the kinetic temperature of the gas at these early epochs. Under the basic assumption that the 21-cm signal is activated by extremely metal-poor stellar systems, we show that the EDGES results are consistent with an extrapolation of the declining galaxy UV luminosity density measured at 4<z<9 by deep HST observations. A substantially enhanced star formation rate density or new exotic sources of UV photons are not required at the redshifts of the EDGES signal. The amount of ionizing radiation produced by the same stellar systems that induce Lyman-alpha coupling is significant, of order 0.5 Lyman-continuum photons per H-atom per 100 Myr. To keep hydrogen largely neutral and delay the reionization process consistently with recent Planck CMB results, mean escape fractions of f_esc < 20% are required at z>15.

The Slug Nebula is one of the largest and most luminous Lyman-alpha (LyA) nebulae discovered to date, extending over 450 kiloparsecs (kpc) around the bright quasar UM287 at z=2.283. Characterized by high surface brightnesses over intergalactic scales, its LyA emission may either trace high-density ionized gas ("clumps") or large column densities of neutral material. To distinguish between these two possibilities, information from a non-resonant line such as Halpha is crucial. Therefore, we analyzed a deep MOSFIRE observation of one of the brightest LyA emitting regions in the Slug Nebula with the goal of detecting associated Halpha emission. We also obtained a deep, moderate resolution LyA spectrum of the nearby brightest region of the Slug. We detected an Halpha flux of F_(Halpha)= 2.62 +/- 0.47 x 10^-17 erg/cm^2/s (SB_(Halpha)=2.70 +/- 0.48 x 10^-18 erg/cm^2/s/sq") at the expected spatial and spectral location. Combining the Halpha detection with its corresponding LyA flux (determined from the narrow-band imaging) we calculate a flux ratio of F_(LyA_/F_(Halpha)= 5.5 +/- 1.1. The presence of a skyline at the location of the Halpha emission decreases the signal to noise ratio of the detection and our ability to put stringent constraints on the Halpha kinematics. Our measurements argue for the origin of the LyA emission being recombination radiation, suggesting the presence of high-density ionized gas. Finally, our high-resolution spectroscopic study of the LyA emission does not show evidence of a rotating disk pattern and suggest a more complex origin for at least some parts of the Slug Nebula.

Recent analyses of the Pal 5 and GD-1 tidal streams suggest that the inner dark matter halo of the Milky Way is close to spherical, in tension with predictions from collisionless N-body simulations of cosmological structure formation. We use the Eris simulation to test whether the combination of dissipative physics and hierarchical structure formation can produce Milky Way-like galaxies whose dark matter halos match the tidal stream constraints from the GD-1 and Pal 5 clusters. We use a dynamical model of the simulated Eris galaxy to generate many realizations of the GD-1 and Pal 5 tidal streams, marginalize over observational uncertainties in the cluster galactocentric positions and velocities, and compare with the observational constraints. We find that the total density and potential of Eris contributed by baryons and dark matter satisfies constraints from the existing Milky Way stellar stream data, as the baryons both round and redistribute the dark matter during the dissipative formation of the galaxy, and provide a centrally-concentrated mass distribution that rounds the inner potential. The Eris dark matter halo or a spherical Navarro-Frenk-White dark matter work comparably well in modeling the stream data. In contrast, the equivalent dark matter-only ErisDark simulation produces a prolate halo that cannot reproduce the observed stream data. The on-going Gaia mission will provide decisive tests of the consistency between $\Lambda$CDM and Milky Way streams, and should distinguish between models like Eris and more spherical halos.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

Recent observations suggest that hydrogen reionization ends late ($z \simeq 6$) and proceeds quickly. We present here a new model of the meta-galactic UV/X-ray background (UVB) that is consistent with this. It adopts the most recent determinations of the ionizing emissivity due to stars and AGN, as well as of the HI absorber column density distribution. Another major improvement is a new treatment of the intergalactic medium (IGM) opacity for ionizing photons that is able to consistently capture the transition from a neutral to an ionized IGM. Previous synthesis models of the UVB, when used in simulations, yield reionization and thermal histories that are inconsistent with the assumed ionizing emissivities. With our new treatment, this discrepancy is fully resolved. In our fiducial model, galaxies leaking $\lesssim 18\%$ of their Lyman continuum emission drive HI reionization, while AGN drive HeII reionization (completing at $z \simeq 6.2$ and $2.8$ respectively). Due to the limited time available for cooling between HI and HeII reionization, higher IGM temperatures are predicted for late reionization scenarios. In our fiducial model, the predicted temperatures agree well with observational constraints at $z \lesssim 4$, while being slightly high compared to (somewhat uncertain) data above that. Models with a larger contribution of AGN are instead disfavoured by the temperature data, as well as by measurements of the HI and HeII Lyman-$\alpha$ forest opacities. We also present "equivalent-equilibrium" ionization/heating rates that mimic our fiducial UVB model for use in simulation codes that assume ionization equilibrium.

While low mass, star forming galaxies are often considered as the primary driver of reionization, their actual contribution to the cosmic ultraviolet background is still uncertain, mostly because the escape fraction of ionizing photons is only poorly constrained. Theoretical studies have shown that efficient supernova feedback is a necessary condition to create paths through which ionizing radiation can escape into the intergalactic medium. We investigate the possibility that accreting supermassive black holes in early dwarf galaxies may provide additional feedback and enhance the leakage of ionizing radiation. We use a series of high resolution cosmological radiation hydrodynamics simulations where we isolate the different sources of feedback. We find that supernova feedback prevents the growth of the black hole, thus quenching its associated feedback. Even in cases where the black hole can grow, the structure of the interstellar medium is strongly dominated by supernova feedback. We conclude that, in the dwarf galaxy regime, supermassive black holes do not appear to play a significant role in enhancing the escape fraction and in contributing to the early UV background.

The Milky Way dark matter halo is formed from the accretion of smaller subhalos. These sub-units also harbor stars---typically old and metal-poor---that are deposited in the Galactic inner regions by disruption events. In this Letter, we show that the dark matter and metal-poor stars in the Solar neighborhood share similar kinematics due to their common origin. Using the high-resolution Eris simulation, which traces the evolution of both the dark matter and baryons in a realistic Milky-Way analog galaxy, we demonstrate that metal-poor stars are indeed effective tracers for the local, virialized dark matter velocity distribution. The local dark matter velocities can therefore be inferred from observations of the stellar halo made by the Sloan Digital Sky Survey within 4 kpc of the Sun. This empirical distribution differs from the Standard Halo Model in important ways and suggests that the bounds on the spin-independent scattering cross section may be weakened for dark matter masses below $\sim$10 GeV. Data from Gaia will allow us to further refine the expected distribution for the smooth dark matter component, and to test for the presence of local substructure.

The reionization of cosmic hydrogen marks a critical juncture in the history of structure formation in the universe. Here we present a new formulation of the standard reionization equation for the evolution of the volume-averaged HII fraction that is more consistent with the accepted conceptual model of inhomogeneous intergalactic absorption. The revised equation retains the basic terminology and simplicity of the classic calculation but explicitly accounts for the presence of the optically thick "Lyman-limit systems" that are known to determine the mean free path of ionizing radiation after overlap. Integration of this equation provides a better characterization of the timing of reionization by smoothly linking the pre-overlap with the post-overlap phases of such process. We confirm the validity of the quasi-instantaneous approximation as predictor of reionization completion/maintenance, and discuss new insights on the sources of cosmic reionization using the improved formalism. A constant emission rate into the intergalactic medium (IGM) of 3 Lyman continuum (LyC) photons per atom per Gyr leads to a reionization history that is consistent with a number of observational constraints on the ionization state of the z=5-9 universe and with the reduced Thomson scattering optical depth recently reported by the Planck Collaboration. While star-forming galaxies can dominate the reionization process if the luminosity-weighted fraction of LyC photons that escape into the IGM, f_esc, exceeds 15% (for a faint magnitude cut-off of the galaxy UV luminosity function of M_lim=-13 and a LyC photon yield per unit 1500 AA luminosity of xi_ion=10^{25.3} Hz/erg, simple models where the product of the two unknowns f_esc xi_ion is not evolving with redshift fail to reproduce the changing neutrality of the IGM observed at these epochs.

We present the results of a survey for intervening MgII absorbers in the redshift range z $\simeq$2-6 in the foreground of four high redshift quasar spectra, 5.79$\le z_{em}\le$6.133, obtained with the ESO VLT X-shooter. We identify 24 absorbers at $\ge 5\sigma$ significance in the equivalent width range 0.117$\le W_{2796}\le 3.655\unicode{xC5}$ with the highest redshift absorber at $z=4.89031\pm4\times10^{-5}$. For weak ($W_{2796}<0.3\unicode{xC5}$) systems, we measure an incidence rate $dN/dz$=1.35$\pm$0.58 at <z>=2.34 and find that it almost doubles to $dN/dz$=2.58$\pm$0.67 by <z>=4.81. Weak absorbers exceeds the number expected from an exponential fit to stronger systems ($W_{2796}>0.3\unicode{xC5}$). We find that there must be significant evolution in the absorption halo properties of MgII absorbers with $W_{2796} >0.1\unicode{xC5}$ by <z>=4.77 and/or that they are associated with galaxies with luminosities beyond the limits of the current luminosity function at z $\sim$5. We find that the incidence rate of strong MgII absorbers ($W_{2796}>1.0\unicode{xC5}$) can be explained if they are associated with galaxies with $L\ge0.29L_{\ast}$ and/or their covering fraction increases. If they continue to only be associated with galaxies with $L\ge0.50L_{\ast}$ then their physical cross section ($\sigma_{phys}$) increases from 0.015 Mpc$^2$ at z=2.3 to 0.041 Mpc$^2$ at <z>=4.77. We measure $\Omega_{MgII}$=2.1$^{+6.3}_{-0.6}\times10^{-8}$, 1.9$^{+2.9}_{-0.2}\times10^{-8}$, 3.9$^{+7.1}_{-2.4}\times10^{-7}$ at <z>=2.48, 3.41, 4.77, respectively. At <z>=4.77, $\Omega_{MgII}$ exceeds the value expected from $\Omega_{HI}$ estimated from the global metallicity of DLAs at z $\simeq$4.85 by a factor of $\sim$44 suggesting that either MgII absorbers trace both ionised and neutral gas and/or are more metal rich than the average DLA at this redshift.

Using the Eris zoom-in cosmological simulation of assembly of a Milky Way analog, we study chemical enrichment of stars due to accretion of metal-enriched gas from the interstellar medium during the Galaxy's development. We consider metal-poor and old stars in both Galactic halo and bulge and make use of stellar orbits, gas density and metallicity distributions in Eris. Assuming spherically symmetric Bondi-Hoyle accretion, we find that halo and bulge stars accrete metals at the rate of about 10^-24 solar mass per year and 10^-22 solar mass per year, respectively, at redshifts z < 3, but this accretion rate increases hundred-fold to about 10^-20 solar mass per year at higher redshifts due to increased gas density. Bulge and halo stars accrete similar amounts of metals at high redshifts as kinematically distinct bulge and halo are not yet developed at these redshifts and both sets of stars encounter similar metal distribution in the ISM on average. Accretion alone can enrich main-sequence stars up to [Fe/H] -2 in extreme cases.Median enrichment level due to accretion in these stars is about [Fe/H]~-6 to -5.Because accretion mostly takes place at high redshifts, it is alpha-enriched to [alpha/Fe]~0.5. We find that accretive metal enrichment is significant enough to affect the predicted metallicity distribution function of halo stars at [Fe/H] < -5.This suggests that attempts to infer the natal chemical environment of the most metal-poor stars from their observed enrichment today can be hindered due to metal accretion. Peculiar enrichment patterns such as those predicted to arise from pair-instability supernovae could help in disentangling natal and accreted metal content of stars.

We compute the expected X-ray diffuse background and radiative feedback on the intergalactic medium (IGM) from X-ray binaries prior and during the epoch of reionization. The cosmic evolution of compact binaries is followed using a population synthesis technique that treats separately neutron stars and black hole binaries in different spectral states and is calibrated to reproduce the observed X-ray properties of galaxies at z<4. Together with an updated empirical determination of the cosmic history of star formation, recent modeling of the stellar mass-metallicity relation, and a scheme for absorption by the IGM that accounts for the presence of ionized HII bubbles during the epoch of reionization, our detailed calculations provide refined predictions of the X-ray volume emissivity and filtered radiation background from "normal" galaxies at z>6. Radiative transfer effects modulate the background spectrum, which shows a characteristic peak between 1 and 2 keV. While the filtering of X-ray radiation through the IGM slightly increases the mean excess energy per photoionization, it also weakens the radiation intensity below 1 keV, lowering the mean photoionization and heating rates. Numerical integration of the rate and energy equations shows that the contribution of X-ray binaries to the ionization of the bulk IGM is negligible, with the electron fraction never exceeding 1%. Direct HeI photoionizations are the main source of IGM heating, and the temperature of the largely neutral medium in between HII cavities increases above the temperature of the cosmic microwave background (CMB) only at z<10, when the volume filling factor of HII bubbles is already >0.1. Therefore, in this scenario, it is only at relatively late epochs that the bulk of neutral intergalactic hydrogen may be observable in 21-cm emission against the CMB.

We present results from the "Ponos" simulation suite on the early evolution of a massive, $M_{\rm vir}(z=0)=1.2\times 10^{13}$ M$_{\odot}$ galaxy. At $z\gtrsim6$, before feedback from a central supermassive black hole becomes dominant, the main galaxy has a stellar mass $\sim 2\times 10^{9}$ M$_{\odot}$ and a star formation rate $\sim 20$ M$_{\odot}$ yr$^{-1}$. The galaxy sits near the expected main sequence of star-forming galaxies at those redshifts, and resembles moderately star-forming systems observed at $z>5$. The high specific star formation rate results in vigorous heating and stirring of the gas by supernovae feedback, and the galaxy develops a thick and turbulent disc, with gas velocity dispersion $\sim 40$ km s$^{-1}$, rotation to dispersion ratio $\sim 2$, and with a significant amount of gas at $\sim 10^5$ K. The Toomre parameter always exceeds the critical value for gravito-turbulence, $Q\sim 1.5-2$, mainly due to the contribution of warm/hot gas inside the disc. Without feedback, a nearly gravito-turbulent regime establishes with similar gas velocity dispersion and lower $Q$. We propose that the "hot and turbulent" disc regime seen in our simulations, unlike the "cold and turbulent" gravito-turbulent regime of massive clumpy disc galaxies at $z\sim 1-2$, is a fundamental characterisation of main sequence galaxies at $z\gtrsim 6$, as they can sustain star formation rates comparable to those of low-mass starbursts at $z=0$. This results in no sustained coherent gas inflows through the disc, and in fluctuating and anisotropic mass transport, possibly postponing the assembly of the bulge and causing the initial feeding of the central black hole to be highly intermittent.

We report the detection of extended Halpha emission from the tip of the HI disk of the nearby edge-on galaxy UGC 7321, observed with the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope. The Halpha surface brightness fades rapidly where the HI column density drops below N(HI) = 10^19 cm^-2 , consistent with fluorescence arising at the ionisation front from gas that is photoionized by the extragalactic ultraviolet background (UVB). The surface brightness measured at this location is (1.2 +/- 0.5)x10^-19 erg/s/cm^2/arcsec^2, where the error is mostly systematic and results from the proximity of the signal to the edge of the MUSE field of view, and from the presence of a sky line next to the redshifted Halpha wavelength. By combining the Halpha and the HI 21 cm maps with a radiative transfer calculation of an exponential disk illuminated by the UVB, we derive a value for the HI photoionization rate of Gamma ~ (6-8)x10^-14 1/s . This value is consistent with transmission statistics of the Lyalpha forest and with recent models of a UVB which is dominated by quasars.

We confirm that the object DDO216-A1 is a substantial globular cluster at the center of Local Group galaxy DDO216 (the Pegasus dwarf irregular), using Hubble Space Telescope ACS imaging. By fitting isochrones, we find the cluster metallicity to be [M/H] = -1.6 +/-0.2, for reddening E(B-V) = 0.16 +/-0.02; the best-fit age is 12.3 +/-0.8 Gyr. There are ~30 RR Lyrae variables in the cluster; the magnitude of the fundamental mode pulsators gives a distance modulus of 24.77 +/-0.08 - identical to the host galaxy. The ratio of overtone to fundamental mode variables and their mean periods make DDO216-A1 an Oosterhoff Type I cluster. We find an I-band central surface brightness 20.85 +/-0.17 F814W mag per square arcsecond, a half-light radius of 3.1 arcsec (13.4 pc), and an absolute magnitude M814 = -7.90 +/-0.16 (approximately 10^5 solar masses). King models fit to the cluster give the core radius and concentration index, r_c = 2.1" +/-0.9" and c = 1.24 +/-0.39. The cluster is an "extended" cluster somewhat typical of some dwarf galaxies and the outer halo of the Milky Way. The cluster is projected <30 pc south of the center of DDO216, unusually central compared to most dwarf galaxy globular clusters. Analytical models of dynamical friction and tidal destruction suggest that it probably formed at a larger distance, up to ~1 kpc, and migrated inward. DDO216 has an unexceptional cluster specific frequency, S_N = 10. DDO216 is the lowest-luminosity Local Group galaxy to host a 10^5 solar mass globular cluster, and the only transition-type (dSph/dIrr) in the Local Group with a globular.

Young stars typically form in star clusters, so the supernovae (SNe) they produce are clustered in space and time. This clustering of SNe may alter the momentum per SN deposited in the interstellar medium (ISM) by affecting the local ISM density, which in turn affects the cooling rate. We study the effect of multiple SNe using idealized 1D hydrodynamic simulations which explore a large parameter space of the number of SNe, and the background gas density and metallicity. The results are provided as a table and an analytic fitting formula. We find that for clusters with up to ~100 SNe the asymptotic momentum scales super-linearly with the number of SNe, resulting in a momentum per SN that can be an order of magnitude larger than for a single SN, with a maximum efficiency for clusters with 10-100 SNe. We argue that additional physical processes not included in our simulations -- self-gravity, breakout from a galactic disk, and galactic shear -- can slightly reduce the momentum enhancement from clustering, but the average momentum per SN still remains a factor of 4 larger than the isolated SN value when averaged over a realistic cluster mass function for a star-forming galaxy. We conclude with a discussion of the possible role of mixing between hot and cold gas, induced by multi-dimensional instabilities or preexisting density variations, as a limiting factor in the buildup of momentum by clustered SNe, and suggest future numerical experiments to explore these effects.

We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretion and is heating the surroundings via thermal AGN feedback. A large central region in the ErisBH disk becomes bar-unstable after z~1.4, but a clear bar-like structure starts to grow significantly only after z~0.4, possibly triggered by the interaction with a massive satellite. At z~0.1 the bar reaches its maximum radial extent of l~2.2 kpc. As the bar grows, it becomes prone to buckling instability, which we quantify based on the anisotropy of the stellar velocity dispersion. The actual buckling event is observable at z~0.1, resulting in the formation of a boxy-peanut bulge clearly discernible in the edge-on view of the galaxy at z=0. The bar in ErisBH does not dissolve during the formation of the bulge but remains strongly non-axisymmetric down to the resolution limit of ~100 pc at z=0. During its early growth, the bar exerts a strong torque on the gas within its extent and drives gas inflows that enhance the nuclear star formation on sub-kpc scales. Later on the infalling gas is nearly all consumed into stars and, to a lesser extent, accreted onto the central black hole, leaving behind a gas-depleted region within the central ~2 kpc. Observations would more likely identify a prominent, large-scale bar at the stage when the galactic central region has already been quenched. Bar-driven quenching may play an important role in disk-dominated galaxies at all redshift. [Abridged]

We study the dependence of fragmentation in massive gas-rich galaxy disks at $z > 1$ on feedback model and hydrodynamical method, employing the GASOLINE2 SPH code and the lagrangian mesh-less code GIZMO in finite mass mode. We compare non-cosmological galaxy disk runs with standard blastwave supernovae (SN)feedback, which introduces delayed cooling in order to drive winds, and runs with the new superbubble SN feedback, which produces winds naturally by modelling the detailed physics of SN-driven bubbles and leads to efficient self-regulation of star formation. We find that, with blastwave feedback, massive star forming clumps form in comparable number and with very similar masses in GASOLINE2 and GIZMO. The typical masses are in the range $10^7-10^8 M_{\odot}$, lower than in most previous works, while giant clumps with masses above $10^9 M_{\odot}$ are exceedingly rare. With superbubble feedback, instead, massive bound star forming clumps do not form because galaxies never undergo a phase of violent disk instability. Only sporadic, unbound star forming overdensities lasting only a few tens of Myr can arise that are triggered by perturbations of massive satellite companions. We conclude that there is a severe tension between explaining massive star forming clumps observed at $z > 1$ primarily as the result of disk fragmentation driven by gravitational instability and the prevailing view of feedback-regulated galaxy formation. The link between disk stability and star formation efficiency should thus be regarded as a key testing ground for galaxy formation theory.

We present initial results from the "Ponos" zoom-in numerical simulations of dark matter substructures in massive ellipticals. Two very highly resolved dark matter halos with $M_{\rm vir}=1.2\times 10^{13}$ $M_{\odot}$ and $M_{\rm vir}=6.5\times 10^{12}$ $M_{\odot}$ and different ("violent" vs. "quiescent") assembly histories have been simulated down to $z=0$ in a $\Lambda$CDM cosmology with a total of 921,651,914 and 408,377,544 particles, respectively. Within the virial radius, the total mass fraction in self-bound $M_{\rm sub}>10^6$ $M_{\odot}$ subhalos at the present epoch is 15% for the violent host and 16.5% for the quiescent one. At $z=0.7$, these fractions increase to 19 and 33%, respectively, as more recently accreted satellites are less prone to tidal destruction. In projection, the average fraction of surface mass density in substructure at a distance of $R/R_{\rm vir}=0.02$ ($\sim 5-10$ kpc) from the two halo centers ranges from 0.6% to $\gtrsim 2$%, significantly higher than measured in simulations of Milky Way-sized halos. The contribution of subhalos with $M_{\rm sub} < 10^9$ $M_{\odot}$ to the projected mass fraction is between one fifth and one third of the total, with the smallest share found in the quiescent host. We assess the impact of baryonic effects via twin, lower-resolution hydrodynamical simulations that include metallicity-dependent gas cooling, star formation, and a delayed-radiative-cooling scheme for supernova feedback. Baryonic contraction produces a super-isothermal total density profile and increases the number of massive subhalos in the inner regions of the main host. The host density profiles and projected subhalo mass fractions appear to be broadly consistent with observations of gravitational lenses.

We use a set of high-resolution N-body simulations of the Galactic disk to study its interactions with the population of satellites predicted cosmologically. One simulation illustrates that multiple passages of massive satellites with different velocities through the disk generate a wobble, having the appearance of rings in face-on projections of the stellar disk. They also produce flares in the disk outer parts and gradually heat the disk through bending waves. A different numerical experiment shows that an individual satellite as massive as the Sagittarius dwarf galaxy passing through the disk will drive coupled horizontal and vertical oscillations of stars in underdense regions, with small significant associated heating. This experiment shows that vertical excursions of stars in these low-density regions can exceed 1 kpc in the Solar neighborhood, resembling the coherent vertical oscillations recently detected locally. They can also induce non-zero vertical streaming motions as large as 10-20 km s$^{-1}$, consistent with recent observations in the Galactic disk. This phenomenon appears as a local ring, with modest associated disk heating.

We simulate the tidal disruption of a collisionless N-body globular star cluster in a total of 300 different orbits selected to have galactocentric radii between 10 and 30 kpc in four dark matter halos: (a) a spherical halo with no subhalos, (b) a spherical halo with subhalos, (c) a realistic halo with no subhalos, and (d) a realistic halo with subhalos. This allows us to isolate and study how the halo's (lack of) dynamical symmetry and substructures affect the dispersal of tidal debris. The realistic halos are constructed from the snapshot of the Via Lactea II simulation at redshift zero. We find that the overall halo's lack of dynamical symmetry disperses tidal debris to make the streams fluffier, consistent with previous studies of tidal debris of dwarf galaxies in larger orbits than ours in this study. On the other hand, subhalos in realistic potentials can locally enhance the densities along streams, making streams denser than their counterparts in smooth potentials. We show that many long and thin streams can survive in a realistic and lumpy halo for a Hubble time. This suggests that upcoming stellar surveys will likely uncover more thin streams which may contain density gaps that have been shown to be promising probes for dark matter substructures.

We investigate the properties of halo gas using three cosmological `zoom-in' simulations of realistic Milky Way-galaxy analogs with varying sub-grid physics. In all three cases, the mass of hot ($T > 10^6$ K) halo gas is $\sim 1\%$ of the host's virial mass. The X-ray luminosity of two of the runs is consistent with observations of the Milky Way, while the third simulation is X-ray bright and resembles more closely a very massive, star-forming spiral. Hot halos extend to 140 kpc from the galactic center and are surrounded by a bubble of warm-hot ($T = 10^5 - 10^6$K) gas that extends to the virial radius. Simulated halos agree well outside 20-30 kpc with the $\beta$-model of Miller & Bregman (2013) based on OVII absorption and OVIII emission measurements. Warm--hot and hot gas contribute up to $80\%$ of the total gas reservoir, and contain nearly the same amount of baryons as the stellar component. The mass of warm-hot and hot components falls into the range estimated for $L^*$ galaxies. With key observational constraints on the density of the Milky Way corona being satisfied, we show that concealing of the ubiquitous warm-hot baryons, along with the ejection of just $20-30 \%$ of the diffuse gas out of the potential wells by supernova-driven outflows, can solve the "missing baryon problem". The recovered baryon fraction within 3 virial radii is $90\%$ of the universal value. With a characteristic density of $\sim 10^{-4}$ cm$^{-3}$ at $50-80$ kpc, diffuse coronae meet the requirement for fast and complete ram--pressure stripping of the gas reservoirs in dwarf galaxy satellites.

The rapid assembly of the massive black holes that power the luminous quasars observed at $z \sim 6-7$ remains a puzzle. Various direct collapse models have been proposed to head-start black hole growth from initial seeds with masses $\sim 10^5\,\rm M_\odot$, which can then reach a billion solar mass while accreting at the Eddington limit. Here we propose an alternative scenario based on radiatively inefficient super-critical accretion of stellar-mass holes embedded in the gaseous circum-nuclear discs (CNDs) expected to exist in the cores of high redshift galaxies. Our sub-pc resolution hydrodynamical simulations show that stellar-mass holes orbiting within the central 100 pc of the CND bind to very high density gas clumps that arise from the fragmentation of the surrounding gas. Owing to the large reservoir of dense cold gas available, a stellar-mass black hole allowed to grow at super-Eddington rates according to the "slim disc" solution can increase its mass by 3 orders of magnitudes within a few million years. These findings are supported by simulations run with two different hydro codes, RAMSES based on the Adaptive Mesh Refinement technique and GIZMO based on a new Lagrangian Godunov-type method, and with similar, but not identical, sub-grid recipes for star formation, supernova feedback, black hole accretion and feedback. The low radiative efficiency of super-critical accretion flows are instrumental to the rapid mass growth of our black holes, as they imply modest radiative heating of the surrounding nuclear environment.

We assess a model of late cosmic reionization in which the ionizing background radiation arises entirely from high redshift quasars and other active galactic nuclei (AGNs). The low optical depth to Thomson scattering reported by the Planck Collaboration pushes the redshift of instantaneous reionization down to z=8.8^{+1.7}_{-1.4} and greatly reduces the need for significant Lyman-continuum emission at very early times. We show that, if recent claims of a numerous population of faint AGNs at z=4-6 are upheld, and the high inferred AGN comoving emissivity at these epochs persists to higher, z~10, redshifts, then active galaxies may drive the reionization of hydrogen and helium with little contribution from normal star-forming galaxies. We discuss an AGN-dominated scenario that satisfies a number of observational constraints: the HI photoionization rate is relatively flat over the range 2<z<5, hydrogen gets fully reionized by z=5.7, and the integrated Thomson scattering optical depth is tau=0.056, in agreement with measurements based on the Lya opacity of the intergalactic medium (IGM) and cosmic microwave background (CMB) polarization. It is a prediction of the model that helium gets doubly reionized before redshift 4, the heat input from helium reionization dominates the thermal balance of the IGM after hydrogen reionization, and z>5 AGNs provide a significant fraction of the unresolved X-ray background at 2 keV. Singly- and doubly-ionized helium contribute about 13% to tau, and the HeIII volume fraction is already 50% when hydrogen becomes fully reionized.

We compare cosmological hydrodynamical simulations combined with the homogeneous metagalactic UV background (UVB) of Haardt & Madau (2012) (HM2012) to observations of the Lyman-alpha forest that are sensitive to the thermal and ionization state of the intergalactic medium (IGM). The transition from optically thick to thin photoheating predicted by the simple one-zone, radiative transfer model implemented by HM2012 predicts a thermal history that is in remarkably good agreement with the observed rise of the IGM temperature at z~3 if we account for the expected evolution of the volume filling factor of HeIII. Our simulations indicate that there may be, however, some tension between the observed peak in the temperature evolution and the rather slow evolution of the HeII opacities suggested by recent Hubble Space Telescope/COS measurements. The HM2012 UVB also underpredicts the metagalactic hydrogen photoionization rate required by our simulations to match the observed opacity of the forest at z>4 and z<2.

We present a new zoom-in hydrodynamical simulation, "Erisbh", which follows the cosmological evolution and feedback effects of a supermassive black hole at the center of a Milky Way-type galaxy. ErisBH shares the same initial conditions, resolution, recipes of gas cooling, star formation and feedback, as the close Milky Way-analog "Eris", but it also includes prescriptions for the formation, growth and feedback of supermassive black holes. We find that the galaxy's central black hole grows mainly through mergers with other black holes coming from infalling satellite galaxies. The growth by gas accretion is minimal because very little gas reaches the sub-kiloparsec scales. The final black hole is, at z=0, about 2.6 million solar masses and it sits closely to the position of SgrA* on the MBH-MBulge and MBH-sigma planes, in a location consistent with what observed for pseudobulges. Given the limited growth due to gas accretion, we argue that the mass of the central black hole should be above 10^5 solar masses already at z~8. The effect of AGN feedback on the host galaxy is limited to the very central few hundreds of parsecs. Despite being weak, AGN feedback seems to be responsible for the limited growth of the central bulge with respect to the original Eris, which results in a significantly flatter rotation curve in the inner few kiloparsecs. Moreover, the disk of ErisBH is more prone to instabilities, as its bulge is smaller and its disk larger then Eris. As a result, the disk of ErisBH undergoes a stronger dynamical evolution relative to Eris and around z=0.3 a weak bar grows into a strong bar of a few disk scale lengths in size. The bar triggers a burst of star formation in the inner few hundred parsecs, provides a modest amount of new fuel to the central black hole, and causes the bulge of ErisBH to have, by z=0, a box/peanut morphology.(Abridged)

Tidal debris from infalling satellites can leave observable structure in the phase-space distribution of the Galactic halo. Such substructure can be manifest in the spatial and/or velocity distributions of the stars in the halo. This paper focuses on a class of substructure that is purely kinematic in nature, with no accompanying spatial features. To study its properties, we use a simulated stellar halo created by dynamically populating the Via Lactea II high-resolution N-body simulation with stars. A significant fraction of the stars in the inner halo of Via Lactea share a common speed and metallicity, despite the fact that they are spatially diffuse. We argue that this kinematic substructure is a generic feature of tidal debris from older mergers and may explain the detection of radial-velocity substructure in the inner halo made by the Sloan Extension for Galactic Understanding and Exploration. The GAIA satellite, which will provide the proper motions of an unprecedented number of stars, should further characterize the kinematic substructure in the inner halo. Our study of the Via Lactea simulation suggests that the stellar halo can be used to map the speed distribution of the local dark-matter halo, which has important consequences for dark-matter direct-detection experiments.

We simulate tidal streams in the presence and absence of substructures inside the zero redshift snapshot of the Via Lactea II (VL-2) simulation. A halo finder is used to remove and isolate the subhalos found inside the high resolution dark matter halo of VL-2, and the potentials for both the main halo and all the subhalos are constructed individually using the self-consistent field (SCF) method. This allows us to make direct comparison of tidal streams between a smooth halo and a lumpy halo without assuming idealized profiles or triaxial fits. We simulate the kinematics of a star cluster starting with the same orbital position but two different velocities. Although these two orbits are only moderately eccentric and have similar apo- and pericentric distances, we find that the two streams have very different morphologies. We conclude that our model of the potential of VL-2 can provide insights about tidal streams that have not been explored by previous studies using idealized or axisymmetric models.

We investigate the production sites and the enrichment history of $r$-process elements in the Galaxy, as traced by the [Eu/Fe] ratio, using the high resolution, cosmological zoom-in simulation `Eris'. At $z=0$, Eris represents a close analog to the Milky Way, making it the ideal laboratory to understand the chemical evolution of our Galaxy. Eris formally traces the production of oxygen and iron due to Type-Ia and Type-II supernovae. We include in post-processing the production of $r$-process elements from compact binary mergers. Unlike previous studies, we find that the nucleosynthetic products from compact binary mergers can be incorporated into stars of very low metallicity and at early times, even with a minimum delay time of 100 Myr. This conclusion is relatively insensitive to modest variations in the merger rate, minimum delay time, and the delay time distribution. By implementing a first-order prescription for metal-mixing, we can further improve the agreement between our model and the data for the chemical evolution of both [$\alpha$/Fe] and [Eu/Fe].We argue that compact binary mergers could be the dominant source of $r$-process nucleosynthesis in the Galaxy.

We report the discovery of an extremely long ($\sim$110 Mpc/$h$) and dark ($\tau_{\rm eff} \gtrsim 7$) Ly$\alpha$ trough extending down to $z \simeq 5.5$ towards the $z_{\rm em} \simeq 6.0$ quasar ULAS J0148+0600. We use these new data in combination with Ly$\alpha$ forest measurements from 42 quasars at $4.5 \le z_{\rm em} \le 6.4$ to conduct an updated analysis of the line-of-sight variance in the intergalactic Ly$\alpha$ opacity over $4 \le z \le 6$. We find that the scatter in transmission among lines of sight near $z \sim 6$ significantly exceeds theoretical expectations for either a uniform ultraviolet background (UVB) or simple fluctuating UVB models in which the mean free path to ionizing photons is spatially invariant. The data, particularly near $z \simeq 5.6$-5.8, instead require fluctuations in the volume-weighted hydrogen neutral fraction that are a factor of 3 or more beyond those expected from density variations alone. We argue that these fluctuations are most likely driven by large-scale variations in the mean free path, consistent with expectations for the final stages of inhomogeneous hydrogen reionization. Even by $z \simeq 5.6$, however, a large fraction of the data are consistent with a uniform UVB, and by $z \sim 5$ the data are fully consistent with opacity fluctuations arising solely from the density field. This suggests that while reionization may be ongoing at $z \sim 6$, it has fully completed by $z \sim 5$.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

We analyze the formation and evolution of the stellar components in "Eris", a 120 pc-resolution cosmological hydrodynamic simulation of a late-type spiral galaxy. The simulation includes the effects of a uniform UV background, a delayed-radiative-cooling scheme for supernova feedback, and a star formation recipe based on a high gas density threshold. It allows a detailed study of the relative contributions of "in-situ" (within the main host) and "ex-situ" (within satellite galaxies) star formation to each major Galactic component in a close Milky Way analog. We investigate these two star-formation channels as a function of galactocentric distance, along different lines of sight above and along the disk plane, and as a function of cosmic time. We find that: 1) approximately 70 percent of today's stars formed in-situ; 2) more than two thirds of the ex-situ stars formed within satellites after infall; 3) the majority of ex-situ stars are found today in the disk and in the bulge; 4) the stellar halo is dominated by ex-situ stars, whereas in-situ stars dominate the mass profile at distances < 5 kpc from the center at high latitudes; and 5) approximately 25% of the inner, r < 20 kpc, halo is composed of in-situ stars that have been displaced from their original birth sites during Eris' early assembly history.

We use detailed nucleosynthesis calculations and a realistic prescription for the environment of the first stars to explore the first episodes of chemical enrichment that occurred during the dark ages. Based on these calculations, we propose a novel explanation for the increased prevalence of carbon-enhanced metal-poor (CEMP) stars with decreasing Fe abundance: The observed chemistry for the most metal-poor Galactic halo stars is the result of an intimate link between the explosions of the first stars and their host minihalo's ability to retain its gas. Specifically, high-energy supernovae produce a near solar ratio of C/Fe, but are effective in evacuating the gas from their host minihalo, thereby suppressing the formation of a second generation of stars. On the other hand, minihalos that host low-energy supernovae are able to retain their gas and form a second stellar generation but, as a result, the second stars are born with a supersolar ratio of C/Fe. Our models are able to accurately reproduce the observed distributions of [C/Fe] and [Fe/H], as well as the fraction of CEMP stars relative to non-CEMP stars as a function of [Fe/H] without any free parameters. We propose that the present lack of chemical evidence for very massive stars (>140 Msun), that ended their lives as a highly energetic pair-instability supernova, does not imply that such stars were rare or did not exist; the chemical products of these very massive first stars may have escaped from their host minihalo, and were never incorporated into subsequent generations of stars. Finally, our models suggest that the most Fe-poor stars currently known may have seen the enrichment from a small multiple of metal-free stars, and need not have been exclusively enriched by a solitary first star. These calculations support the idea that some of the surviving dwarf satellite galaxies of the Milky Way are relics of the first galaxies.

On dwarf galaxy scales, the different shapes of the galaxy stellar mass function and the dark halo mass function require a star-formation efficiency (SFE) in these systems that is currently more than 1 dex lower than that of Milky Way-size halos. Here, we argue that this trend may actually be reversed at high redshift. Specifically, by combining the resolved star-formation histories of nearby isolated dwarfs with the simulated mass-growth rates of dark matter halos, we show that the assembly of these systems occurs in two phases: (1) an early, fast halo accretion phase with a rapidly deepening potential well, characterized by a high SFE; and (2) a late slow halo accretion phase where, perhaps as a consequence of reionization, the SFE is low. Nearby dwarfs have more old stars than predicted by assuming a constant or decreasing SFE with redshift, a behavior that appears to deviate qualitatively from the trends seen amongst more massive systems. Taken at face value, the data suggest that, at sufficiently early epochs, dwarf galaxy halos above the atomic cooling mass limit can be among the most efficient sites of star formation in the universe.

We present results from a fully cosmological, very high-resolution, LCDM "zoom-in" simulation of a group of seven field dwarf galaxies with present-day virial masses in the range M_vir=4.4e8-3.6e10 Msun. The simulation includes a blastwave scheme for supernova feedback, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, a scheme for the turbulent diffusion of metals and thermal energy, and a uniform UV background that modifies the ionization and excitation state of the gas. The properties of the simulated dwarfs are strongly modulated by the depth of the gravitational potential well. All three halos with M_vir < 1e9 Msun are devoid of stars, as they never reach the density threshold for star formation of 100 atoms/cc. The other four, M_vir > 1e9 Msun dwarfs have blue colors, low star formation efficiencies, high cold gas to stellar mass ratios, and low stellar metallicities. Their bursty star formation histories are characterized by peak specific star formation rates in excess of 50-100 1/Gyr, far outside the realm of normal, more massive galaxies, and in agreement with observations of extreme emission-line starbursting dwarfs by the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey. Metal-enriched galactic outflows produce sub-solar effective yields and pollute with heavy elements a Mpc-size region of the intergalactic medium, but are not sufficient to completely quench star formation activity and are not ubiquitous in our dwarfs. Within the limited size of the sample, our simulations appear to simultaneously reproduce the observed stellar mass and cold gas content, resolved star formation histories, stellar kinematics, and metallicities of field dwarfs in the Local Volume.

Over the past two decades, an avalanche of data from multiwavelength imaging and spectroscopic surveys has revolutionized our view of galaxy formation and evolution. Here we review the range of complementary techniques and theoretical tools that allow astronomers to map the cosmic history of star formation, heavy element production, and reionization of the Universe from the cosmic "dark ages" to the present epoch. A consistent picture is emerging, whereby the star-formation rate density peaked approximately 3.5 Gyr after the Big Bang, at z~1.9, and declined exponentially at later times, with an e-folding timescale of 3.9 Gyr. Half of the stellar mass observed today was formed before a redshift z = 1.3. About 25% formed before the peak of the cosmic star-formation rate density, and another 25% formed after z = 0.7. Less than ~1% of today's stars formed during the epoch of reionization. Under the assumption of a universal initial mass function, the global stellar mass density inferred at any epoch matches reasonably well the time integral of all the preceding star-formation activity. The comoving rates of star formation and central black hole accretion follow a similar rise and fall, offering evidence for co-evolution of black holes and their host galaxies. The rise of the mean metallicity of the Universe to about 0.001 solar by z = 6, one Gyr after the Big Bang, appears to have been accompanied by the production of fewer than ten hydrogen Lyman-continuum photons per baryon, a rather tight budget for cosmological reionization.

We present more results from a fully cosmological LCDM simulation of a group of isolated dwarf galaxies that has been shown to reproduce the observed stellar mass and cold gas content, resolved star formation histories, and metallicities of dwarfs in the Local Volume. Here we investigate the energetics and timetable of the cusp-core transformation. As suggested by previous work, supernova-driven gas outflows remove dark matter (DM) cusps and create kpc-size cores in all systems having a stellar mass above 1e6 Msun. The "DM core mass removal efficiency" -- dark mass ejected per unit stellar mass -- ranges today from a few to a dozen, and increases with decreasing host mass. Because dwarfs form the bulk of their stars prior to redshift 1 and the amount of work required for DM heating and core formation scales approximately as Mvir^{5/3}, the unbinding of the DM cusp starts early and the formation of cored profiles is not as energetically onerous as previously claimed. DM particles in the cusp typically migrate to 2-3 core radii after absorbing a few percent of the energy released by supernovae. The present-day slopes of the inner dark matter mass profiles, Gamma=d log M/dlog R=2.5-3, of the simulated "Bashful" and "Doc" dwarfs are similar to those measured in the luminous Fornax and Sculptor dwarf spheroidals. None of the simulated galaxies has a circular velocity profile exceeding 20 km/s in the inner 1 kpc, implying that supernova feedback is key to solve the "too-big-to-fail" problem for Milky Way subhalos.

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

Observations of diffuse Galactic gamma ray emission (DGE) by the Fermi Large Area Telescope (LAT) allow a detailed study of cosmic rays and the interstellar medium. However, diffuse emission models of the inner Galaxy underpredict the Fermi-LAT data at energies above a few GeV and hint at possible non-astrophysical sources including dark matter (DM) annihilations or decays. We present a study of the possible emission components from DM using the high-resolution Via Lactea II N-body simulation of a Milky Way-sized DM halo. We generate full-sky maps of DM annihilation and decay signals that include modeling of the adiabatic contraction of the host density profile, Sommerfeld enhanced DM annihilations, $p$-wave annihilations, and decaying DM. We compare our results with the DGE models produced by the Fermi-LAT team over different sky regions, including the Galactic center, high Galactic latitudes, and the Galactic anti-center. This work provides possible templates to fit the observational data that includes the contribution of the subhalo population to DM gamma-ray emission, with the significance depending on the annihilation/decay channels and the Galactic regions being considered.

We consider super-critical accretion with angular momentum onto stellar-mass black holes as a possible mechanism for growing billion-solar-mass holes from light seeds at early times. We use the radiatively-inefficient "slim disk" solution -- advective, optically thick flows that generalize the standard geometrically thin disk model -- to show how mildly super-Eddington intermittent accretion may significantly ease the problem of assembling the first massive black holes when the Universe was less than 0.8 Gyr old. Because of the low radiative efficiencies of slim disks around non-rotating as well as rapidly rotating holes, the mass e-folding timescale in this regime is nearly independent of the spin parameter. The conditions that may lead to super-critical growth in the early Universe are briefly discussed.

We present an analysis of the effects of dissipational baryonic physics on the local dark matter (DM) distribution at the location of the Sun, with an emphasis on the consequences for direct detection experiments. Our work is based on a comparative analysis of two cosmological simulations of a Milky Way halo: Eris, a full hydrodynamic simulation, and ErisDark, its DM-only realization. We find that two distinct processes lead in Eris to a 30% enhancement of DM in the disk plane at the location of the Sun: the accretion and disruption of satellites resulting in a DM component with net angular momentum and the contraction of baryons pulling DM into the disk plane without forcing it to co-rotate. Owing to its particularly quiescent merger history for dark halos of Milky Way mass, the co-rotating dark disk in Eris is less massive than what has been suggested by previous work, contributing only 9% of the local DM density. The speed distribution in Eris is broadened and shifted to higher speeds compared to its DM-only twin simulation ErisDark. At high speeds $f(v)$ falls more steeply in Eris than in ErisDark or the SHM, easing the tension between recent results from the CDMS-II and XENON100 experiments. The non-Maxwellian aspects of $f(v)$ are still present, but much less pronounced in Eris than in DM-only runs. The weak dark disk increases the time-averaged scattering rate by only a few percent at low recoil energies. On the high velocity tail, however, the increase in typical speeds due to baryonic contraction results in strongly enhanced mean scattering rates compared to ErisDark, although they are still suppressed compared to the SHM. Similar trends are seen regarding the amplitude of the annual modulation, while the modulated fraction is increased compared to the SHM and decreased compared to ErisDark.

We examine recent measurements of the z~2.5 intergalactic medium (IGM) which constrain the HI frequency distribution f(N) and the mean free path MFP to ionizing radiation. We argue that line-blending and the clustering of strong absorption-line systems have led previous authors to systematically overestimate the effective Lyman limit opacity, yielding too small of a MFP for the IGM. We further show that recently published measurements of f(N) at log NHI ~ 16 lie in strong disagreement, implying underestimated uncertainty from sample variance and/or systematics like line-saturation. Allowing for a larger uncertainty in the f(N) measurements, we provide a new cubic Hermite spline model for f(N) which reasonably satisfies all of the observational constraints under the assumption of randomly distributed absorption systems. We caution, however, that this formalism is invalid in light of absorber clustering and use a toy model to estimate the effects. Future work must properly account for the non-Poissonian nature of the IGM.

Simulations of structure formation in the Universe predict that galaxies are embedded in a "cosmic web", where the majority of baryons reside as rarefied and highly ionized gas. This material has been studied for decades in absorption against background sources, but the sparseness of these inherently one-dimensional probes preclude direct constraints on the three-dimensional morphology of the underlying web. Here we report observations of a cosmic web filament in Lyman-alpha emission, discovered during a survey for cosmic gas fluorescently "illuminated" by bright quasars at z=2.3. With a projected size of approximately 460 physical kpc, the Lyman-alpha emission surrounding the radio-quiet quasar UM287 extends well beyond the virial radius of any plausible associated dark matter halo. The estimated cold gas mass of the nebula from the observed emission is at least ten times larger than what is typically found by cosmological simulations, suggesting that a population of intergalactic gas clumps with sub-kpc sizes may be missing within current numerical models.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

If "seed" central black holes were common in the subgalactic building blocks that merged to form present-day massive galaxies, then relic intermediate-mass black holes (IMBHs) should be present in the Galactic bulge and halo. We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with black holes, and assess the size, properties, and detectability of the leftover population. The method assigns a black hole to the most tightly bound central particle of each subhalo at infall according to an extrapolation of the M_BH-sigma_* relation, and self-consistently follows the accretion and disruption of Milky Way progenitor dwarfs and their holes in a cosmological "live" host from high redshift to today. We show that, depending on the minimum stellar velocity dispersion, sigma_m, below which central black holes are assumed to be increasingly rare, as many as ~2000 (sigma_m=3 km/s) or as few as ~70 (sigma_m=12 km/s) IMBHs may be left wandering in the halo of the Milky Way today. The fraction of IMBHs kicked out of their host by gravitational recoil is < 20%. We identify two main Galactic subpopulations, "naked" IMBHs, whose host subhalos were totally destroyed after infall, and "clothed" IMBHs residing in dark matter satellites that survived tidal stripping. Naked IMBHs typically constitute 40-50% of the total and are more centrally concentrated. We show that, in the sigma_m=12 km/s scenario, the clusters of tightly bound stars that should accompany naked IMBHs would be fainter than m_V=16 mag, spatially resolvable, and have proper motions of 0.1-10 milliarcsec per year. Their detection may provide an observational tool to constrain the formation history of massive black holes in the early Universe.

We present a cosmological hydrodynamic simulation of the formation of dwarf galaxies at redshifts z>~2.5 using a physically-motivated model for H2-regulated star formation. Our simulation, performed using the Enzo code and reaching a peak resolution of 109 proper parsecs at z=2.5, extends the results of Kuhlen et al. (2012) to significantly lower redshifts. We show that a star formation prescription regulated by the local H2 abundance leads to the suppression of star formation in dwarf galaxy halos with M_h <~ 10^10 Msun and to a large population of gas-rich "dark galaxies" at z=2.5 with low star formation efficiencies and gas depletion timescales >20 Gyr. The fraction of dark galaxies is 60% at M_h ~ 10^10 Msun and increases rapidly with decreasing halo mass. Dark galaxies form late and their gaseous disks never reach the surface densities, > ~5700 Msun / pc^2 (Z/10^-3 Zsun)^(-0.88), that are required to build a substantial molecular fraction. Despite this large population of dark galaxies, we show that our H2-regulated simulation is consistent with both the observed luminosity function of galaxies and the cosmological mass density of neutral gas at z>~2.5. Moreover, our results provide a theoretical explanation for the recent detection in fluorescent Ly-alpha emission of gaseous systems at high redshift with little or no associated star formation. We further propose that H2-regulation may offer a fresh solution to a number of outstanding "dwarf galaxy problems" in LambdaCDM. In particular, H2-regulation leads galaxy formation to become effectively stochastic on mass scales of M_h ~ 10^10 Msun, and thus these massive dwarfs are not "too big to fail".

We discuss a novel approach to "weighing" the Milky Way dark matter halo, one that combines the latest samples of halo stars selected from the Sloan Digital Sky Survey (SDSS) with state-of-the-art numerical simulations of Milky Way analogs. The fully cosmological runs employed in the present study include "Eris", one of the highest-resolution hydrodynamical simulations of the formation of a M_vir=8e11 M_sun late-type spiral, and the dark-matter only M_vir=1.7e12 M_sun "Via Lactea II" simulation. Eris provides an excellent laboratory for creating mock SDSS samples of tracer halo stars, and we successfully compare their density, velocity anisotropy, and radial velocity dispersion profiles with the observational data. Most mock SDSS realizations show the same "cold veil" recently observed in the distant stellar halo of the Milky Way, with tracers as cold as sigma_los ~ 50 km/s between 100-150 kpc. Controlled experiments based on the integration of the spherical Jeans equation as well as a particle tagging technique applied to Via Lactea II show that a "heavy" M_vir 2e12 M_sun realistic host produces a poor fit to the kinematic SDSS data. We argue that these results offer added evidence for a "light," centrally-concentrated Milky Way halo.

We analyze the present-day structure and assembly history of a high resolution hydrodynamic simulation of the formation of a Milky Way (MW)-like disk galaxy, from the "Eris" simulation suite, dissecting it into cohorts of stars formed at different epochs of cosmic history. At z=0, stars with t_form < 2 Gyr mainly occupy the stellar spheroid, with the oldest (earliest forming) stars having more centrally concentrated profiles. The younger age cohorts populate disks of progressively longer radial scale length and shorter vertical scale height. At a given radius, the vertical density profiles and velocity dispersions of stars vary smoothly as a function of age, and the superposition of old, vertically-extended and young, vertically-compact cohorts gives rise to a double-exponential profile like that observed in the MW. Turning to formation history, we find that the trends of spatial structure and kinematics with stellar age are largely imprinted at birth, or immediately thereafter. Stars that form during the active merger phase at z>3 are quickly scattered into rounded, kinematically hot configurations. The oldest disk cohorts form in structures that are radially compact and relatively thick, while subsequent cohorts form in progressively larger, thinner, colder configurations from gas with increasing levels of rotational support. The disk thus forms "inside-out" in a radial sense and "upside-down" in a vertical sense. Secular heating and radial migration influence the final state of each age cohort, but the changes they produce are small compared to the trends established at formation. The predicted correlations of stellar age with spatial and kinematic structure are in good qualitative agreement with the correlations observed for mono-abundance stellar populations in the MW.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We present the first science results from our Hubble Space Telescope Survey for Lyman limit absorption systems (LLS) using the low dispersion spectroscopic modes of the Advanced Camera for Surveys and the Wide Field Camera 3. Through an analysis of 71 quasars, we determine the incidence frequency of LLS per unit redshift and per unit path length, l(z) and l(x) respectively, over the redshift range 1 < z< 2.6, and find a weighted mean of l(x)=0.29 +/-0.05 for 2.0 < z < 2.5 through a joint analysis of our sample and that of Ribaudo et al. (2011). Through stacked spectrum analysis, we determine a median (mean) value of the mean free path to ionizing radiation at z=2.4 of lambda_mfp = 243(252)h^(-1) Mpc, with an error on the mean value of +/- 43h^(-1) Mpc. We also re-evaluate the estimates of lambda_mfp from Prochaska et al. (2009) and place constraints on the evolution of lambda_mfp with redshift, including an estimate of the "breakthrough" redshift of z = 1.6. Consistent with results at higher z, we find that a significant fraction of the opacity for absorption of ionizing photons comes from systems with N_HI <= 10^{17.5} cm^(-2) with a value for the total Lyman opacity of tau_lyman = 0.40 +/- 0.15. Finally, we determine that at minimum, a 5-parameter (4 power-law) model is needed to describe the column density distribution function f(N_HI, X) at z \sim 2.4, find that f(N_HI,X) undergoes no significant change in shape between z \sim 2.4 and z \sim 3.7, and provide our best fit model for f(N_HI,X).

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift 3, using the "Eris" suite of cosmological "zoom-in" simulations. The reference run adopts a blastwave scheme for supernova feedback that produces galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. Synthetic spectra through the multiphase CGM produce interstellar absorption line strengths of Lya, CII, CIV, SiII, and SiIV as a function of galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high-redshift by Steidel et al. (2010). Only about one third of all the gas within R_vir is outflowing. The fraction of sightlines within one virial radius that intercept optically thick material is 27%, in agreement with recent observations by Rudie et al. (2012). Such optically thick absorption is shown to trace inflowing "cold" streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar, give origin to strong (log N > 13) CII absorption with a covering factor of 22% within R_vir and 10% within 2 R_vir. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large OVI halo, with a typical column density log N>14 and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. (2011). Our zoom-in simulations of this single system appear to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings. (Abridged)

We show that the position of the central dark matter density peak may be expected to differ from the dynamical center of the Galaxy by several hundred parsec. In Eris, a high resolution cosmological hydrodynamics simulation of a realistic Milky-Way-analog disk galaxy, this offset is 300 - 400 pc (~3 gravitational softening lengths) after z=1. In its dissipationless dark-matter-only twin simulation ErisDark, as well as in the Via Lactea II and GHalo simulations, the offset remains below one softening length for most of its evolution. The growth of the DM offset coincides with a flattening of the central DM density profile in Eris inwards of ~1 kpc, and the direction from the dynamical center to the point of maximum DM density is correlated with the orientation of the stellar bar, suggesting a bar-halo interaction as a possible explanation. A dark matter density offset of several hundred parsec greatly affects expectations of the dark matter annihilation signals from the Galactic Center. It may also support a dark matter annihilation interpretation of recent reports by Weniger (2012) and Su & Finkbeiner (2012) of highly significant 130 GeV gamma-ray line emission from a region 1.5 degrees (~200 parsec projected) away from Sgr A* in the Galactic plane.

A galaxy halo may contain a large number of intermediate mass black holes (IMBHs) with masses in the range of 10^{2-6} solar mass. We propose to directly detect these IMBHs by observing multiply imaged QSO-galaxy or galaxy-galaxy strong lens systems in the submillimeter bands with high angular resolution. The silhouette of an IMBH in the lensing galaxy halo would appear as either a monopole-like or a dipole-like variation at the scale of the Einstein radius against the Einstein ring of the dust-emitting region surrounding the QSO. We use a particle tagging technique to dynamically populate a Milky Way-sized dark matter halo with black holes, and show that the surface mass density and number density of IMBHs have power-law dependences on the distance from the center of the host halo if smoothed on a scale of ~ 1 kpc. Most of the black holes orbiting close to the center are freely roaming as they have lost their dark matter hosts during infall due to tidal stripping. Next generation submillimeter telescopes with high angular resolution (< 0.3 mas) will be capable of directly mapping such off-nuclear freely roaming IMBHs with a mass of ~ 10^6 solar mass in a lensing galaxy that harbours a O(10^9) solar mass supermassive black hole in its nucleus.

We investigate the formation and evolution of the pseudobulge in "Eris", a high-resolution N-body + smoothed particle hydrodynamics (SPH) cosmological simulation that successfully reproduces a Milky Way-like massive late-type spiral in a cold dark matter (LCDM) universe. At the present epoch, Eris has a virial mass Mvir=8x10^11 Msun, a photometric stellar mass M*=3.2x10^10 Msun, a bulge-to-total ratio B/T = 0.26, and a weak nuclear bar. We find that the bulk of the pseudobulge forms quickly at high redshift via a combination of non-axisymmetric disk instabilities and tidal interactions or mergers both occurring on dynamical timescales, not through slow secular processes at lower redshift. Its subsequent evolution is not strictly secular either, and is closely intertwined with the evolution of the stellar bar. In fact, the structure that we recognize as a pseudobulge today evolves from a stellar bar that formed at high redshift, was destroyed by minor mergers at z~3, reformed shortly after, and weakened again following a steady gas inflow at z<1. The gradual dissolution of the bar ensues at z=1 and continues until the present without increasing the stellar velocity dispersion in the inner regions. In this scenario the pseudobulge is not a separate component from the inner disk in terms of formation path, rather it is the first step in the inside-out formation of the baryonic disk, in agreement with the fact that pseudobulges of massive spiral galaxies have typically a dominant old stellar population. If our simulations do indeed reproduce the formation mechanisms of massive spirals, then the progenitors of late-type galaxies should have strong bars and small photometric pseudobulges at high redshift.

We present a detailed study of the metal-enriched circum-galactic medium of a massive galaxy at z=3 using "ErisMC", a new cosmological hydrodynamic "zoom-in" simulations of a disk galaxy with mass comparable to the Milky Way. The run adopts a blastwave scheme for supernova feedback that generates galactic outflows without explicit wind particles, a star formation recipe based on a high gas density threshold, and high temperature metal cooling. ErisMC's main progenitor at z=3 resembles a "Lyman break" galaxy of mass M_vir=2.4e11 M_sun, virial radius R_vir=48 kpc, and star formation rate 18 M_sun/yr, and its metal-enriched CGM extends to 200 (physical) kpc from its center. Approximately 41, 9, and 50 percent of all gas-phase metals at z=3 are locked in a hot (T> 3e5 K), warm (3e5 >T> 3e4 K), and cold (T< 3e4 K) medium, respectively. We identify three sources of heavy elements: 1) the main host, responsible for 60% of all the metals found within 3R_vir; 2) its satellite progenitors, responsible for 28% of all the metals within 3R_vir, and for only 5% of those beyond 3R_vir; and nearby dwarfs, which give origin to 12% of all the metals within 3R_vir and 95% of those beyond 3R_vir. Late (z<5) galactic "superwinds" account for only 9% of all the metals observed beyond 2R_vir, the bulk having been released at redshifts 5< z < 8 by early star formation and outflows. In the CGM, lower overdensities are typically enriched by `older', colder metals. Heavy elements are accreted onto Eris along filaments via low-metallicity cold inflows, and are ejected hot via galactic outflows at a few hundred km/s. The outflow mass-loading factor is of order unity for the main halo, but can exceed 10 for nearby dwarfs. We stress that our "zoom-in" simulation focuses on the CGM of a single massive system and cannot describe the enrichment history of the IGM as a whole. (abridged)

We describe cosmological galaxy formation simulations with the adaptive mesh refinement code Enzo that incorporate a star formation prescription regulated by the local abundance of molecular hydrogen. We show that this H2-regulated prescription leads to a suppression of star formation in low mass halos (M_h < ~10^10 M_sun) at z>4, alleviating some of the dwarf galaxy problems faced by theoretical galaxy formation models. H2 regulation modifies the efficiency of star formation of cold gas directly, rather than indirectly reducing the cold gas content with "supernova feedback". We determine the local H2 abundance in our most refined grid cells (76 proper parsec in size at z=4) by applying the model of Krumholz, McKee, & Tumlinson, which is based on idealized 1D radiative transfer calculations of H2 formation-dissociation balance in ~100 pc atomic--molecular complexes. Our H2-regulated simulations are able to reproduce the empirical (albeit lower z) Kennicutt-Schmidt relation, including the low Sigma_gas cutoff due to the transition from atomic to molecular phase and the metallicity dependence thereof, without the use of an explicit density threshold in our star formation prescription. We compare the evolution of the luminosity function, stellar mass density, and star formation rate density from our simulations to recent observational determinations of the same at z=4-8 and find reasonable agreement between the two.

We find that the relative contribution of satellite galaxies accreted at high redshift to the stellar population of the Milky Way's smooth halo increases with distance, becoming observable relative to the classical smooth halo about 15 kpc from the Galactic center. In particular, we determine line-of-sight-averaged [Fe/H] and [alpha/Fe] in the metal-poor main-sequence turnoff (MPMSTO) population along every Sloan Extension for Galactic Understanding and Exploration (SEGUE) spectroscopic line of sight. Restricting our sample to those lines of sight along which we do not detect elements of cold halo substructure (ECHOS), we compile the largest spectroscopic sample of stars in the smooth component of the halo ever observed in situ beyond 10 kpc. We find significant spatial autocorrelation in [Fe/H] in the MPMSTO population in the distant half of our sample beyond about 15 kpc from the Galactic center. Inside of 15 kpc however, we find no significant spatial autocorrelation in [Fe/H]. At the same time, we perform SEGUE-like observations of N-body simulations of Milky Way analog formation. While we find that halos formed entirely by accreted satellite galaxies provide a poor match to our observations of the halo within 15 kpc of the Galactic center, we do observe spatial autocorrelation in [Fe/H] in the simulations at larger distances. This observation is an example of statistical chemical tagging and indicates that spatial autocorrelation in metallicity is a generic feature of stellar halos formed from accreted satellite galaxies.

We present improved synthesis models of the evolving spectrum of the UV/X-ray diffuse background, updating and extending our previous results. Five new main components are added to our radiative transfer code CUBA: (1) the sawtooth modulation of the background intensity from resonant line absorption in the Lyman series of cosmic hydrogen and helium; (2) the X-ray emission from obscured and unobscured quasars; (3) a piecewise parameterization of the distribution in redshift and column density of intergalactic absorbers that fits recent measurements of the mean free path of 1 ryd photons; (4) an accurate treatment of the photoionization structure of absorbers; and (5) the UV emission from star-forming galaxies at all redshifts. We provide tables of the predicted HI and HeII photoionization and photoheating rates for use, e.g., in cosmological hydrodynamics simulations of the Lya forest, and a new metallicity-dependent calibration to the UV luminosity density-star formation rate density relation. A "minimal cosmic reionization model" is also presented in which the galaxy UV emissivity traces recent determinations of the cosmic history of star formation, the luminosity-weighted escape fraction of hydrogen-ionizing radiation increases rapidly with lookback time, the clumping factor of the high-redshift intergalactic medium evolves following the results of recent hydrodynamic simulations, and Population III stars and miniquasars make a negligible contribution to the metagalactic flux. The model provides a good fit to the hydrogen-ionization rates inferred from flux decrement and proximity effect measurements, predicts that cosmological HII (HeIII) regions overlap at redshift 6.7 (2.8), and yields an optical depth to Thomson scattering that is in agreement with WMAP results. (Abridged)

We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with stars. The method is calibrated using the observed luminosity function of Milky Way satellites and the concentration of their stellar populations, and self-consistently follows the accretion and disruption of progenitor dwarfs and the build-up of the stellar halo in a cosmological "live host". Simple prescriptions for assigning stellar populations to collisionless particles are able to reproduce many properties of the observed Milky Way halo and its surviving dwarf satellites, like velocity dispersions, sizes, brightness profiles, metallicities, and spatial distribution. Our model predicts the existence of approximately 1,850 subhalos harboring "extremely faint" satellites (with mass-to-light ratios >5,000) lying beyond the Sloan Digital Sky Survey detection threshold. Of these, about 20 are "first galaxies", i.e. satellites that formed a stellar mass above 10 Msun before redshift 9. The ten most luminous satellites (L> 1e6 Lsun) in the simulation are hosted by subhalos with peak circular velocities today in the range V_max=10-40 km/s that have shed between 80% and 99% of their dark mass after being accreted at redshifts 1.7< z <4.6. The satellite maximum circular velocity and stellar line-of-sight velocity dispersion today follow the relation V_max=2.2 sigma_los. We apply a standard mass estimation algorithm based on Jeans modelling of the line-of-sight velocity dispersion profiles to the simulated dwarf spheroidals, and test the accuracy of this technique. The inner (within 300 pc) mass-luminosity relation for currently detectable satellites is nearly flat in our model, in qualitative agreement with the "common mass scale" found in Milky Way dwarfs. We do, however, predict a weak, but significant positive correlation for these objects: M_300 ~L^{0.088 \pm 0.024}.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

Simulations of the formation of late-type spiral galaxies in a cold dark matter LCDM universe have traditionally failed to yield realistic candidates. Here we report a new cosmological N-body/SPH simulation of extreme dynamic range in which a close analog of a Milky Way disk galaxy arises naturally. Termed Eris, the simulation follows the assembly of a galaxy halo of mass Mvir=7.9x10^11 Msun with a total of N=18.6 million particles (gas + dark matter + stars) within the final virial radius, and a force resolution of 120 pc. It includes radiative cooling, heating from a cosmic UV field and supernova explosions, a star formation recipe based on a high gas density threshold (nSF=5 atoms cm^-3 rather than the canonical nSF=0.1 atoms cm^-3), and neglects AGN feedback. At the present epoch, the simulated galaxy has an extended rotationally-supported disk with a radial scale length Rd=2.5 kpc, a gently falling rotation curve with circular velocity at 2.2 disk scale lenghts of V2.2=214 km/s, a bulge-to-disk ratio B/D=0.35, and a baryonic mass fraction that is 30% below the cosmic value. The disk is thin, is forming stars in the region of the Sigma_SFR - Sigma_HI plane occupied by spiral galaxies, and falls on the photometric Tully-Fisher and the stellar mass-halo virial mass relations. Hot (T>3x10^5 K), X-ray luminous halo gas makes only 26% of the universal baryon fraction and follows a flattened density profile proportional to r^-1.13 out to r=100 kpc. Eris appears then to be the first cosmological hydrodynamic simulation in which the galaxy structural properties, the mass budget in the various components, and the scaling relations between mass and luminosity are all consistent with a host of observational constraints. (Abridged)

We find a bright (L_{UV}=2.5 L*_{z=6}) Lyman alpha emitter at redshift z=5.719 at a projected distance of 79 physical kpc from a strong triply ionized carbon (Civ) absorption system at redshift z=5.7238 previously reported in the spectrum of the z_{em} = 6.309 QSO SDSS J1030+0524. This is the highest redshift galaxy-absorber pair detected to-date, supporting the idea that galaxy-wide outflows were already in place at the end of the epoch of reionization. The proximity of this object makes it the most likely source of metals, consistent with models of outflows at lower redshift where significant observational evidence relates metal absorption systems with galaxies hosting outflows. In a typical outflow scenario, a wind of 200 km/s, active since the universe was only 0.6 Gyr old (z ~8.4), could eject metals out to 79 kpc at z=5.719. Although the origin of metals in the intergalactic medium (IGM) is still under debate, our results are consistent with predictions from cosmological simulations which reproduce the evolution of the cosmic density of Civ, from z ~ 6 to the present day based on outflow-driven enrichment of the IGM. We also report two more Lyman alpha emitters in this field, at z=5.973\pm 0.002 and z=5.676\pm 0.002 respectively, the former confirming the original identification by Stiavelli et al. Our results suggest that the colour cut typically used to identify i-dropouts (i_{775}-z_{850}>1.3) misses a non-negligible fraction of blue galaxies with faint UV continuum at z \geq 5.7.

We present improved synthesis models of the evolving spectrum of the UV/X-ray diffuse background, updating and extending our previous results. A "minimal cosmic reionization model" is presented in which the galaxy UV emissivity traces recent determinations of the cosmic history of star formation, the luminosity-weighted escape fraction of hydrogen-ionizing radiation increases rapidly with redshift, from 1% at z=3 to 50% at z=8.5, the clumping factor of the intergalactic medium is C=3, and Population III stars make a negligible contribution to the metagalactic flux. The model provides a good fit to the hydrogen-ionization rates inferred from flux decrement measurements, and yields an optical depth to Thomson scattering, tau_e=0.085, that is in agreement with WMAP results.

In the near future, measurements of metal absorption features in the intergalactic medium (IGM) will become an important constraint on models of the formation and evolution of the earliest galaxies, the properties of the first stars, and the reionization and enrichment of the IGM. The first measurement of a metal abundance in the IGM at a redshift approaching the epoch of reionization already offers intriguing hints. Between z=5.8 and 4.7 (a 0.3 Gyr interval only 1 Gyr after the big bang), the measured density of CIV absorbers in the IGM increased by a factor of ~ 3.5 (Ryan-Weber et al. 2009; Becker, Rauch & Sargent 2009). If these values prove to be accurate, they pose two puzzles: (1) The total amount of CIV at z=5.8 implies too little star formation to reionize the IGM by z=6 or to match the WMAP electron scattering optical depth (tau). (2) The rapid growth from z = 6-5 is faster than the buildup of stellar mass or the increase in the star formation rate density over the same interval. We show that a delay of ~ 0.4-0.7 Gyr between the instantaneous production of ionizing photons and the later production of metal absorption features (added to the delay due to stellar lifetimes) can provide the full explanation for both puzzles. We calculate the delay in metal production due to finite stellar lifetimes alone and find that it is too short to explain the rapid CIV density increase. The additional delay could naturally be explained as the result of ~ 200 km/s outflows carrying carbon to distances of ~ 100 kpc, the typical distance between galaxies and CIV absorbers in enrichment simulations, and the typical outflow or absorption region scale observed at z ~ 2-3.

Tidal stellar disruptions have traditionally been discussed as a probe of the single, massive black holes (MBHs) that are dormant in the nuclei of galaxies. In Chen et al. (2009), we used numerical scattering experiments to show that three-body interactions between bound stars in a stellar cusp and a non-evolving "hard" MBH binary will also produce a burst of tidal disruptions, caused by a combination of the secular "Kozai effect" and by close resonant encounters with the secondary hole. Here we derive basic analytical scalings of the stellar disruption rates with the system parameters, assess the relative importance of the Kozai and resonant encounter mechanisms as a function of time, discuss the impact of general relativistic (GR) and extended stellar cusp effects, and develop a hybrid model to self-consistently follow the shrinking of an MBH binary in a stellar background, including slingshot ejections and tidal disruptions. In the case of a fiducial binary with primary hole mass M_1=10^7\msun and mass ratio q=M_2/M_1=1/81, embedded in an isothermal cusp, we derive a stellar disruption rate \dot{N_*} ~ 0.2/yr lasting ~ 3X10^5 yr. This rate is 3 orders of magnitude larger than the corresponding value for a single MBH fed by two-body relaxation, confirming our previous findings. For q<<0.01, the Kozai/chaotic effect could be quenched due to GR/cusp effects by an order of magnitude, but even in this case the stellar-disruption rate is still two orders of magnitude larger than that given by standard relaxation processes around a single MBH. Our results suggest that >~10% of the tidal-disruption events may originate in MBH binaries.

The asymmetric emission of gravitational waves produced during the coalescence of a massive black hole (MBH) binary imparts a velocity "kick" to the system that can displace the hole from the center of its host. Here we study the trajectories and observability of MBHs recoiling in three (one major, two minor) gas-rich galaxy merger remnants that were previously simulated at high resolution, and in which the pairing of the MBHs had been shown to be successful. We run new simulations of MBHs recoiling in the major merger remnant with Mach numbers in the range 1<M<6, and use simulation data to construct a semi-analytical model for the orbital evolution of MBHs in gas-rich systems. We show that: 1) in major merger remnants the energy deposited by the moving hole into the rotationally supported, turbulent medium makes a negligible contribution to the thermodynamics of the gas. This contribution becomes significant in minor merger remnants, potentially allowing for an electromagnetic signature of MBH recoil; 2) in major merger remnants, the combination of both deeper central potential well and drag from high-density gas confines even MBHs with kick velocities as high as 1200 km/s within 1 kpc from the host's center; 3) kinematically offset nuclei may be observable for timescales of a few Myr in major merger remnants in the case of recoil velocities in the range 700-1,000 km/s; 4) in minor mergers remnants the effect of gas drag is weaker, and MBHs with recoil speeds in the range 300-600 km/s will wander through the host halo for longer timescales. When accounting for the probability distribution of kick velocities, however, we find that the likelihood of observing recoiling MBHs in gas-rich galaxy mergers is very low, typically below 10^-5 - 10^-6.

We present a study of the ability of the Fermi Gamma-ray Space Telescope to detect dark-matter annihilation signals from the Galactic subhalos predicted by the Via Lactea II N-body simulation. We implement an improved formalism for estimating the boost factor needed to account for the effect of dark-matter clumping on scales below the resolution of the simulation, and we incorporate a detailed Monte Carlo simulation of the response of the Fermi-LAT telescope, including a simulation of its all-sky observing mode integrated over a ten year mission. We find that for WIMP masses up to about 150 GeV in standard supersymmetric models with velocity-averaged cross section 3*10^-26 cm^3 s^-1, a few subhalos could be detectable with >5 standard deviations significance and would likely deviate significantly from the appearance of a point source.

The velocity distribution function of dark matter particles is expected to show significant departures from a Maxwell-Boltzmann distribution. This can have profound effects on the predicted dark matter - nucleon scattering rates in direct detection experiments, especially for dark matter models in which the scattering is sensitive to the high velocity tail of the distribution, such as inelastic dark matter (iDM) or light (few GeV) dark matter (LDM), and for experiments that require high energy recoil events, such as many directionally sensitive experiments. Here we determine the velocity distribution functions from two of the highest resolution numerical simulations of Galactic dark matter structure (Via Lactea II and GHALO), and study the effects for these scenarios. For directional detection, we find that the observed departures from Maxwell-Boltzmann increase the contrast of the signal and change the typical direction of incoming DM particles. For iDM, the expected signals at direct detection experiments are changed dramatically: the annual modulation can be enhanced by more than a factor two, and the relative rates of DAMA compared to CDMS can change by an order of magnitude, while those compared to CRESST can change by a factor of two. The spectrum of the signal can also change dramatically, with many features arising due to substructure. For LDM the spectral effects are smaller, but changes do arise that improve the compatibility with existing experiments. We find that the phase of the modulation can depend upon energy, which would help discriminate against background should it be found.

For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as low angular momentum baryons and DM sink to the center of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different DM particle candidate. This Letter presents new hydrodynamical simulations in a Lambda$CDM framework where analogues of dwarf galaxies, bulgeless and with a shallow central DM profile, are formed. This is achieved by resolving the inhomogeneous interstellar medium, resulting in strong outflows from supernovae explosions which remove low angular momentum gas. This inhibits the formation of bulges and decreases the dark-matter density to less than half within the central kiloparsec. Realistic dwarf galaxies are thus shown to be a natural outcome of galaxy formation in the CDM scenario.

The coalescence of a massive black hole (MBH) binary leads to the gravitational-wave recoil of the system and its ejection from the galaxy core. We have carried out N-body simulations of the motion of a MBH = 3.7x10^6 Msun MBH remnant in the Via Lactea I simulation, a Milky Way sized dark matter halo. The black hole receives a recoil velocity of Vkick = 80, 120, 200, 300, and 400 km/s at redshift 1.5, and its orbit is followed for over 1 Gyr within a live host halo, subject only to gravity and dynamical friction against the dark matter background. We show that, owing to asphericities in the dark matter potential, the orbit of the MBH is hightly non-radial, resulting in a significantly increased decay timescale compared to a spherical halo. The simulations are used to construct a semi-analytic model of the motion of the MBH in a time-varying triaxial Navarro-Frenk-White dark matter halo plus a spherical stellar bulge, where the dynamical friction force is calculated directly from the velocity dispersion tensor. Such a model should offer a realistic picture of the dynamics of kicked MBHs in situations where gas drag, friction by disk stars, and the flattening of the central cusp by the returning black hole are all negligible effects. We find that MBHs ejected with initial recoil velocities Vkick > 500 km/s do not return to the host center within Hubble time. In a Milky Way-sized galaxy, a recoiling hole carrying a gaseous disk of initial mass ~MBH may shine as a quasar for a substantial fraction of its wandering phase. The long decay timescales of kicked MBHs predicted by this study may thus be favorable to the detection of off-nuclear quasar activity.

While the connection between high-redshift star formation and the local universe has recently been used to understand the observed population of faint dwarf galaxies in the Milky Way (MW) halo, we explore how well these nearby objects can probe the epoch of first light. We construct a detailed, physically motivated model for the MW satellites based on the state-of-the-art Via Lactea II dark-matter simulations. Our model incorporates molecular hydrogen (H_2) cooling in low-mass systems and inhomogeneous photo-heating feedback during the internal reionization of our own galaxy. We find that the existence of MW satellites fainter than M_V ~ -5 is strong evidence for H_2 cooling in low-mass halos, while satellites with -5 > M_V > -9 were affected by hydrogen cooling and photoheating feedback. The age of stars in very low-luminosity systems and the minimum luminosity of these satellites are key predictions of our model. Most of the stars populating the brightest MW satellites could have formed after the epoch of reionization. Our models also predict a significantly larger dispersion in M_300 values than observed and a number of luminous satellites with M_300 as low as 10^6 solar masses.

We present new calculations of the evolving UV background out to the epoch of cosmological reionization and make predictions for the amount of GeV gamma-ray attenuation by electron-positron pair production. Our results are based on recent semi-analytic models of galaxy formation, which provide predictions of the dust-extinguished UV radiation field due to starlight, and empirical estimates of the contribution due to quasars. We account for the reprocessing of ionizing photons by the intergalactic medium. We test whether our models can reproduce estimates of the ionizing background at high redshift from flux decrement analysis and proximity effect measurements from quasar spectra, and identify a range of models that can satisfy these constraints. Pair-production against soft diffuse photons leads to a spectral cutoff feature for gamma rays observed between 10 and 100 GeV. This cutoff varies with redshift and the assumed star formation and quasar evolution models. We find only negligible amounts of absorption for gamma rays observed below 10 GeV for any emission redshift. With observations of high-redshift sources in sufficient numbers by the Fermi Gamma-ray Space Telescope and new ground-based instruments it should be possible to constrain the extragalactic background light in the UV and optical portion of the spectrum.

The unambiguous detection of Galactic dark matter annihilation would unravel one of the most outstanding puzzles in particle physics and cosmology. Recent observations have motivated models in which the annihilation rate is boosted by the Sommerfeld effect, a non-perturbative enhancement arising from a long range attractive force. Here we apply the Sommerfeld correction to Via Lactea II, a high resolution N-body simulation of a Milky-Way-size galaxy, to investigate the phase-space structure of the Galactic halo. We show that the annihilation luminosity from kinematically cold substructure can be enhanced by orders of magnitude relative to previous calculations, leading to the prediction of gamma-ray fluxes from up to hundreds of dark clumps that should be detectable by the Fermi satellite.

"Hard" massive black hole (MBH) binaries embedded in steep stellar cusps can shrink via three-body slingshot interactions. We show that this process will inevitably be accompanied by a burst of stellar tidal disruptions, at a rate that can be several orders of magnitude larger than that appropriate for a single MBH. Our numerical scattering experiments reveal that: 1) a significant fraction of stars initially bound to the primary hole are scattered into its tidal disruption loss cone by gravitational interactions with the secondary hole, an enhancement effect that is more pronounced for very unequal-mass binaries; 2) about 25% (40%) of all strongly interacting stars are tidally disrupted by a MBH binary of mass ratio q=1/81 (q=1/243) and eccentricity 0.1; and 3) two mechanisms dominate the fueling of the tidal disruption loss cone, a Kozai non-resonant interaction that causes the secular evolution of the stellar angular momentum in the field of the binary, and the effect of close encounters with the secondary hole that change the stellar orbital parameters in a chaotic way. For a hard MBH binary of 10^7 solar masses and mass ratio 0.01, embedded in an isothermal stellar cusp of velocity dispersion sigma*=100 km/s, the tidal disruption rate can be as large as 1/yr. This is 4 orders of magnitude higher than estimated for a single MBH fed by two-body relaxation. When applied to the case of a putative intermediate-mass black hole inspiraling onto Sgr A*, our results predict tidal disruption rates ~0.05-0.1/yr.

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee outlines some of the outstanding questions regarding the assembly history of Massive Black Holes in the nuclei of galaxies and the revolutionary contributions anticipated in this field from low-frequency gravitational wave astronomy.

We present the results of the largest survey to date for intergalactic metals at redshifts z > 5, using near-IR spectra of nine QSOs with emission redshifts z(em) > 5.7. We find, for the first time, a change in the comoving mass density of C IV ions as we look back to redshifts z > 5. At a mean <z> = 5.76, we deduce Omega(C IV)=(4.4+/-2.6)x10^(-9) which implies a drop by a factor of about 3.5 compared to the value at z < 4.7, after accounting for the differing sensitivities of different surveys. The observed number of C IV doublets is also lower by a similar factor, compared to expectations for a non-evolving column density distribution of absorbers. These results point to a rapid build-up of intergalactic C IV over a period of only 300 Myr; such a build-up could reflect the accumulation of metals associated with the rising levels of star formation activity from z = 9 indicated by galaxy counts, and/or an increasing degree of ionisation of the intergalactic medium (IGM), following the overlap of ionisation fronts from star-forming regions. If the value of Omega(C IV) we derive is typical of the IGM at large, it would imply a metallicity Z(IGM) = 10^(-4) Z(Sun). The early-type stars responsible for synthesising these metals would have emitted only about one Lyman continuum photon per baryon prior to z = 5.8; such a background is insufficient to keep the IGM ionised and we speculate on possible factors which could make up the required shortfall.

One of the next decade's most exciting prospects is to explore the cosmic "dark ages," during which the first stars in the Universe formed, with the 21 cm line of neutral hydrogen. At z>6, this light redshifts into the low-frequency radio regime and can be observed with new generations of wide-field arrays. These experiments have the potential to observe enormous volumes and to improve both cosmological and astrophysical constraints dramatically. In this white paper, we describe how the next decade will see the first steps toward unlocking this potential, allowing us to answer two key questions about the fundamental physics of our Universe and the intergalactic medium: Does the standard cosmological model describe the Universe during the "dark ages?" How did the intergalactic medium evolve during this important time, ending with the reionization of hydrogen?

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The cosmic dark ages and the epoch of reionization, during which the first generations of stars and galaxies formed, are among the most compelling frontiers of extragalactic astrophysics and cosmology. Here we describe an exciting new probe of these eras: the 21 cm line of neutral hydrogen, which will allow us to map the neutral intergalactic medium throughout the era of first galaxy formation. In this white paper, we describe how these studies can answer two key questions about galaxy formation: What were the properties of high-redshift galaxies? How did they affect the Universe around them?

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

Cosmic ultraviolet background radiation between 3 and 4 Ryd is reprocessed by resonant line absorption in the Lyman series of intergalactic He II. This process results in a sawtooth modulation of the radiation spectrum from the He II Lya frequency to the Lyman limit. The size of this modulation is a sensitive probe of the epoch of helium reionization and of the sources that keep the intergalactic medium (IGM) highly ionized. For large absorption opacities, the background intensity will peak at frequencies just above each resonance, go to zero at resonance, and fluctuate greatly just below resonance. The He II sawtooth modulation may be one of the missing ingredients needed in the modelling of the abundances of metal ions like C III and Si IV observed in the IGM at redshift 3.

The final inspiral and coalescence of a black hole binary can produce highly beamed gravitational wave radiation. To conserve linear momentum, the black hole remnant can recoil with "kick" velocity as high as 4000 km/s. We present two sets of full N-body simulations of recoiling massive black holes (MBH) in high-resolution, non-axisymmetric potentials. The host to the first set of simulations is the main halo of the Via Lactea I simulation (Diemand et al. 2007). The nature of the resulting orbits is investigated through a numerical model where orbits are integrated assuming an evolving, triaxial NFW potential, and dynamical friction is calculated directly from the velocity dispersion along the major axes of the main halo of Via Lactea I. By comparing the triaxial case to a spherical model, we find that the wandering time spent by the MBH is significantly increased due to the asphericity of the halo. For kicks larger than 200 km/s, the remnant MBH does not return to the inner 200 pc within 1 Gyr, a timescale an order of magnitude larger than the upper limit of the estimated QSO lifetime. The second set of simulations is run using the outcome of a high-resolution gas-rich merger (Mayer et al. 2007) as host potential. In this case, a recoil velocity of 500 km/s cannot remove the MBH from the nuclear region.

We use the recently completed one billion particle Via Lactea II LambdaCDM simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo can not be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.

We present two models for the cosmological UV background light, and calculate the opacity of GeV gamma--rays out to redshift 9. The contributors to the background include 2 possible quasar emissivities, and output from star--forming galaxies as determined by recent a semi--analytic model (SAM) of structure formation. The SAM used in this work is based upon a hierarchical build-up of structure in a $\Lambda$CDM universe and is highly successful in reproducing a variety of observational parameters. Above 1 Rydberg energy, ionizing radiation is subject to reprocessing by the IGM, which we treat using our radiative transfer code, CUBA. The two models for quasar emissivity differing above z = 2.3 are chosen to match the ionization rates observed using flux decrement analysis and the higher values of the line-of-sight proximity effect. We also investigate the possibility of a flat star formation rate density at z $>5$. We conclude that observations of gamma--rays from 10 to 100 GeV by Fermi (GLAST) and the next generation of ground based experiments should confirm a strongly evolving opacity from $1<$ z $<4$. Observation of attenuation in the spectra of gamma--ray bursts at higher redshift could constrain emission of UV radiation at these early times, either from a flat or increasing star-formation density or an unobserved population of sources.

The discovery of hypervelocity binary stars (HVBs) in the Galactic halo would provide definite evidence of the existence of a massive black hole companion to Sgr A*. Here we use an hybrid approach to compute the rate of ejection and the total number of HVBs produced by a hypothetical intermediate-mass black hole (IMBH, M_2<10^5\msun) orbiting Sgr A*. Depending on the mass of M_2 and on the properties of binary stars in the central parsec of the Milky Way, we show that the number of undisrupted HVBs expected to be expelled from the Galactic Center before binary black hole coalescence ranges from zero to a few dozens at most. Therefore, the non-detection of stellar binaries in a complete survey of hypervelocity stars would not rule out the occurrence of an IMBH-Sgr A* in-spiralling event within the last few 10^8 years.

It is a clear unique prediction of the cold dark matter paradigm of cosmological structure formation that galaxies form hierarchically and are embedded in massive, extended dark halos teeming with self-bound substructure or "subhalos". The amount and spatial distribution of subhalos around their host provide unique information and clues on the galaxy assembly process and the nature of the dark matter. Here we present results from the Via Lactea INCITE simulation, a one billion particle, one million cpu-hour simulation of the formation and evolution of a Galactic dark matter halo and its substructure population.

Our recently completed one billion particle Via Lactea II simulation of a Milky Way-sized dark matter halo resolves over 50,000 gravitationally bound clumps orbiting today within the virialized region of the main host. About 2,300 of these subhalos have one or more "progenitors" above 1e6 Msun at redshift 11, i.e. massive enough for their gas to have cooled via excitation of H2 and fragmented prior to the epoch of cosmic reionization. We count 4,500 such progenitors: if these were able to convert a fraction of their gas content into very metal-poor stars with a Salpeter initial mass function (IMF), they would be shining today with a visual magnitude M_V=6.7 per solar mass in stars. Assuming a universal baryon fraction, we show that mean star formation efficiencies as low as 0.1% in progenitors below 1e8 Msun would overproduce the abundance of the faint Galatic dwarf spheroidals observed by the Sloan Digital Sky Survey. Star formation at first light must have occurred either with an IMF lacking stars below 0.9 Msun, or was intrinsically very inefficient in small dark matter halos. If the latter, our results may be viewed as another hint of a minimum scale in galaxy formation.

During the final phases of inspiral, a massive black hole (MBH) binary experiences a recoil due to the asymmetric emission of gravitational waves. We use recent results from numerical relativity simulations together with models of the assembly and growth of MBHs in hierarchical cosmologies, to study the dynamics, statistics, and observability of recoling MBHs. We find that, at redshift z<3, kicked non-rotating holes are typically found between 1 and 30 kpc from their galaxy centers, while rapidly rotating ones are typically between 10 and a few hundred kpc. A recoiling hole that carries an accretion disk may shine as an "off-nuclear AGN" while it moves away from the center of its host galaxy. We predict that, depending on the hole spin distribution and the duration of their active phase, a population of off-nuclear AGN may already be detectable at low and intermediate redshifts in present deep Hubble Space Telescope observations. The James Webb Space Telescope may discover tens of wandering AGN per square degree, most of them moving within their host halos on unbound trajectories.

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.

We develop a formalism for studying the dynamics of massive black hole binaries embedded in gravitationally-bound stellar cusps, and study the binary orbital decay by three-body interactions, the impact of stellar slingshots on the density profile of the inner cusp, and the properties of the ejected hypervelocity stars (HVSs). We find that the scattering of bound stars shrinks the binary orbit and increases its eccentricity more effectively than that of unbound ambient stars. Binaries with initial eccentricities e>0.3 and/or unequal-mass companions (M_2/M_1<0.1) can decay by three-body interactions to the gravitational wave emission regime in less than a Hubble time. The stellar cusp is significantly eroded, and cores as shallow as \rho\propto r^-0.7 may develop from a pre-existing singular isothermal density profile. A population of HVSs is ejected in the host galaxy halo, with a total mass ~M_2. We scale our results to the scattering of stars bound to Sgr A*, the massive black hole in the Galactic Center, by an inspiraling companion of intermediate mass. Depending on binary mass ratio, eccentricity, and initial slope of the stellar cusp, a core of radius ~0.1 pc typically forms in 1-10 Myr. On this timescale about 500-2500 HVSs are expelled with speeds sufficiently large to escape the gravitational potential of the Milky Way.

We present quantitative predictions for the detectability of individual Galactic dark matter subhalos in gamma-rays from dark matter pair annihilations in their centers. Our method is based on a hybrid approach, employing the highest resolution numerical simulations available (including the recently completed one billion particle Via Lactea II simulation) as well as analytical models for the extrapolation beyond the simulations' resolution limit. We include a self-consistent treatment of subhalo boost factors, motivated by our numerical results, and a realistic treatment of the expected backgrounds that individual subhalos must outshine. We show that for reasonable values of the dark matter particle physics parameters (M_X ~ 50 - 500 GeV and <sigma*v> ~ 10^-26 - 10^-25 cm^3/s) GLAST may very well discover a few, even up to several dozen, such subhalos, at 5 sigma significance, and some at more than 20 sigma. We predict that the majority of luminous sources would be resolved with GLAST's expected angular resolution. For most observer locations the angular distribution of detectable subhalos is consistent with a uniform distribution across the sky. The brightest subhalos tend to be massive (median Vmax of 24 km/s) and therefore likely hosts of dwarf galaxies, but many subhalos with Vmax as low as 5 km/s are also visible. Typically detectable subhalos are 20 - 40 kpc from the observer, and only a small fraction are closer than 10 kpc. The total number of observable subhalos has not yet converged in our simulations, and we estimate that we may be missing up to 3/4 of all detectable subhalos.

In cold dark matter cosmological models, structures form and grow by merging of smaller units. Numerical simulations have shown that such merging is incomplete; the inner cores of halos survive and orbit as "subhalos" within their hosts. Here we report a simulation that resolves such substructure even in the very inner regions of the Galactic halo. We find hundreds of very concentrated dark matter clumps surviving near the solar circle, as well as numerous cold streams. The simulation reveals the fractal nature of dark matter clustering: Isolated halos and subhalos contain the same relative amount of substructure and both have cuspy inner density profiles. The inner mass and phase-space densities of subhalos match those of recently discovered faint, dark matter-dominated dwarf satellite galaxies and the overall amount of substructure can explain the anomalous flux ratios seen in strong gravitational lenses. Subhalos boost gamma-ray production from dark matter annihilation, by factors of 4-15, relative to smooth galactic models. Local cosmic ray production is also enhanced, typically by a factor 1.4, but by more than a factor of ten in one percent of locations lying sufficiently close to a large subhalo. These estimates assume that gravitational effects of baryons on dark matter substructure are small.

We present a simple and efficient method to set up spherical structure models for N-body simulations with a multimass technique. This technique reduces by a substantial factor the computer run time needed in order to resolve a given scale as compared to single-mass models. It therefore allows to resolve smaller scales in N-body simulations for a given computer run time. Here, we present several models with an effective resolution of up to 1.68 x 10^9 particles within their virial radius which are stable over cosmologically relevant time-scales. As an application, we confirm the theoretical prediction by Dehnen (2005) that in mergers of collisonless structures like dark matter haloes always the cusp of the steepest progenitor is preserved. We model each merger progenitor with an effective number of particles of approximately 10^8 particles. We also find that in a core-core merger the central density approximately doubles whereas in the cusp-cusp case the central density only increases by approximately 50%. This may suggest that the central region of flat structures are better protected and get less energy input through the merger process.

We use the ``Via Lactea'' simulation to study the co-evolution of a Milky Way-size LambdaCDM halo and its subhalo population. While most of the host halo mass is accreted over the first 6 Gyr in a series of major mergers, the physical mass distribution [not M_vir(z)] remains practically constant since z=1. The same is true in a large sample of LambdaCDM galaxy halos. Subhalo mass loss peaks between the turnaround and virialization epochs of a given mass shell, and declines afterwards. 97% of the z=1 subhalos have a surviving bound remnant at the present epoch. The retained mass fraction is larger for initially lighter subhalos: satellites with maximum circular velocities Vmax=10 km/s at z=1 have today about 40% of their mass back then. At the first pericenter passage a larger average mass fraction is lost than during each following orbit. Tides remove mass in substructure from the outside in, leading to higher concentrations compared to field halos of the same mass. This effect, combined with the earlier formation epoch of the inner satellites, results in strongly increasing subhalo concentrations towards the Galactic center. We present individual evolutionary tracks and present-day properties of the likely hosts of the dwarf satellites around the Milky Way. The formation histories of ``field halos'' that lie today beyond the Via Lactea host are found to strongly depend on the density of their environment. This is caused by tidal mass loss that affects many field halos on eccentric orbits.

The Via Lactea simulation of the dark matter halo of the Milky Way predicts the existence of many thousands of bound subhalos distributed approximately with equal mass per decade of mass. Here we show that: a) a similar steeply rising subhalo mass function is also present at redshift 0.5 in an elliptical-sized halo simulated with comparable resolution in a different cosmology. Compared to Via Lactea, this run produces nearly a factor of two more subhalos with large circular velocities; b) the fraction of Via Lactea mass brought in by subhalos that have a surviving bound remnant today with present-day peak circular velocity Vmax>2 km/s (>10 km/s) is 45% (30%); c) because of tidal mass loss, the number of subhalos surviving today that reached a peak circular velocity of >10 km/s throughout their lifetime exceeds half a thousand, five times larger than their present-day abundance and more than twenty times larger than the number of known satellites of the Milky Way; e) unless the circular velocity profiles of Galactic satellites peak at values significantly higher that expected from the stellar line-of-sight velocity dispersion, only about one in five subhalos with Vmax>20 km/s today must be housing a luminous dwarf; f) small dark matter clumps appear to be relatively inefficient at forming stars even well beyond the virial radius; g) the observed Milky Way satellites appear to follow the overall dark matter distribution of Via Lactea, while the largest simulated subhalos today are found preferentially at larger radii; h) subhalos have central densities that increase with Vmax and reach 0.1-0.3 Msun/pc3 comparable to the central densities inferred in dwarf spheroidals with core radii >250 pc.

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

Supermassive black holes (SMBHs) are a ubiquitous component of the nuclei of galaxies. It is normally assumed that, following the merger of two massive galaxies, a SMBH binary will form, shrink due to stellar or gas dynamical processes and ultimately coalesce by emitting a burst of gravitational waves. However, so far it has not been possible to show how two SMBHs bind during a galaxy merger with gas due to the difficulty of modeling a wide range of spatial scales. Here we report hydrodynamical simulations that track the formation of a SMBH binary down to scales of a few light years following the collision between two spiral galaxies. A massive, turbulent nuclear gaseous disk arises as a result of the galaxy merger. The black holes form an eccentric binary in the disk in less than a million years as a result of the gravitational drag from the gas rather than from the stars.

These lectures summarize our basic understanding of the emergence of cosmic structures, the epoch of first light, the survival of early substructure in present-day galaxy halos.

Hypervelocity stars (HVSs) ejected by the massive black hole at the Galactic center have unique kinematic properties compared to other halo stars. Their trajectories will deviate from being exactly radial because of the asymmetry of the Milky Way potential produced by the flattened disk and the triaxial dark matter halo, causing a change of angular momentum that can be much larger than the initial small value at injection. We study the kinematics of HVSs and propose an estimator of dark halo triaxiality that is determined only by instantaneous position and velocity vectors of HVSs at large Galactocentric distances (r>~50kpc). We show that, in the case of a substantially triaxial halo, the distribution of deflection angles (the angle between the stellar position and velocity vector) for HVSs on bound orbits is spread uniformly over the range 10--180deg. Future astrometric and deep wide-field surveys should measure the positions and velocities of a significant number of HVSs, and provide useful constraints on the shape of the Galactic dark matter halo.

We present a study of the shapes, orientations, and alignments of Galactic dark matter subhalos in the ``Via Lactea'' simulation of a Milky Way-size LCDM host halo. Whereas isolated dark matter halos tend to be prolate, subhalos are predominantly triaxial. Overall subhalos are more spherical than the host halo, with minor to major and intermediate to major axis ratios of 0.68 and 0.83, respectively. Like isolated halos, subhalos tend to be less spherical in their central regions. The principal axis ratios are independent of subhalo mass, when the shapes are measured within a physical scale like r_Vmax, the radius of the peak of the circular velocity curve. Subhalos tend to be slightly more spherical closer to the host halo center. The spatial distribution of the subhalos traces the prolate shape of the host halo when they are selected by the largest V_max they ever had, i.e. before they experienced strong tidal mass loss. The subhalos' orientation is not random: the major axis tends to align with the direction towards the host halo center. This alignment disappears for halos beyond 3 r_200 and is more pronounced when the shapes are measured in the outer regions of the subhalos. The radial alignment is preserved during a subhalo's orbit and they become elongated during pericenter passage, indicating that the alignment is likely caused by the host halo's tidal forces. These tidal interactions with the host halo act to make subhalos rounder over time.

Hypervelocity stars (HVSs) are a natural consequence of the presence of a massive nuclear black hole (Sgr A*) in the Galactic Center. Here we use the Brown et al. sample of unbound and bound HVSs together with numerical simulations of the propagation of HVSs in the Milky Way halo to constrain three plausible ejection mechanisms: 1) the scattering of stars bound to Sgr A* by an inspiraling intermediate-mass black hole (IMBH); 2) the disruption of stellar binaries in the tidal field of Sgr A*; and 3) the two-body scattering of stars off a cluster of stellar-mass black holes orbiting Sgr A*. We compare the predicted radial and velocity distributions of HVSs with the limited-statistics dataset currently available, and show that the IMBH model appears to produce a spectrum of ejection velocities that is too flat. Future astrometric and deep wide-field surveys of HVSs should shed unambiguous light on the stellar ejection mechanism and probe the Milky Way potential on scales as large as 200 kpc.

We discuss the possibility of GLAST detecting gamma-rays from the annihilation of neutralino dark matter in the Galactic halo. We have used "Via Lactea", currently the highest resolution simulation of Galactic cold dark matter substructure, to quantify the contribution of subhalos to the annihilation signal. We present a simulated allsky map of the expected gamma-ray counts from dark matter annihilation, assuming standard values of particle mass and cross section. In this case GLAST should be able to detect the Galactic center and several individual subhalos.

We study the long-term evolution of massive black hole binaries (MBHBs) at the centers of galaxies using detailed scattering experiments to solve the full three-body problem. Ambient stars drawn from a isotropic Maxwellian distribution unbound to the binary are ejected by the gravitational slingshot. We construct a minimal, hybrid model for the depletion of the loss cone and the orbital decay of the binary, and show that secondary slingshots - stars returning on small impact parameter orbits to have a second super-elastic scattering with the MBHB - may considerably help the shrinking of the pair in the case of large binary mass ratios. In the absence of loss-cone refilling by two-body relaxation or other processes, the mass ejected before the stalling of a MBHB is half the binary reduced mass. About 50% of the ejected stars are expelled ejected in a "burst" lasting ~1E4 yrs M_6^1/4, where M_6 is the binary mass in units of 1E6 Msun. The loss cone is completely emptied in a few bulge crossing timescales, 1E7 yrs M_6^1/4. Even in the absence of two-body relaxation or gas dynamical processes, unequal mass and/or eccentric binaries with M_6 >0.1 can shrink to the gravitational wave emission regime in less than a Hubble time, and are therefore "safe" targets for the planned Laser Interferometer Space Antenna (LISA).

The astrophysical processes that led to the formation of the first seed black holes and to their growth into the supermassive variety that powers bright quasars at redshift 6 are poorly understood. In standard LCDM hierarchical cosmologies, the earliest massive holes (MBHs) likely formed at redshift z~15 at the centers of low-mass (M>5e5 solar masses) dark matter ``minihalos'', and produced hard radiation by accretion. FUV/X-ray photons from such ``miniquasars'' may have permeated the universe more uniformly than EUV radiation, reduced gas clumping, and changed the chemistry of primordial gas. The role of accreting seed black holes in determining the thermal and ionization state of the intergalactic medium depends on the amount of cold and dense gas that forms and gets retained in protogalaxies after the formation of the first stars. The highest resolution N-body simulation to date of Galactic substructure shows that subhalos below the atomic cooling mass were very inefficient at forming stars.

We present initial results from ``Via Lactea'', the highest resolution simulation to date of Galactic CDM substructure. It follows the formation of a Milky Way-size halo with Mvir=1.8x10^12 Msun in a WMAP 3-year cosmology, using 234 million particles. Over 10,000 subhalos can be identified at z=0: Their cumulative mass function is well-fit by N(>Msub)= 0.0064 (Msub/Mvir)^(-1) down to Msun=4x10^6 Msun. The total mass fraction in subhalos is 5.3%, while the fraction of surface mass density in substructure within a projected distance of 10 kpc from the halo center is 0.3%. Because of the significant contribution from the smallest resolved subhalos, these fractions have not converged yet. Sub-substructure is apparent in all the larger satellites, and a few dark matter lumps are resolved even in the solar vicinity. The number of dark satellites with peak circular velocities above 10 km/s (5 km/s) is 124 (812): of these, 5 (26) are found within 0.1 Rvir, a region that appeared practically smooth in previous simulations. The neutralino self-annihilation gamma-ray emission from dark matter clumps is approximately constant per subhalo mass decade. Therefore, while in our run the contribution of substructure to the gamma-ray luminosity of the Galactic halo amounts to only 40% of the total spherically-averaged smooth signal, we expect this fraction to grow significantly as resolution is increased further. An all-sky map of the expected annihilation gamma-ray flux reaching a fiducial observer at 8 kpc from the Galactic center shows that at the current resolution a small number of subhalos start to be bright enough to be visible against the background from the smooth density field surrounding the observer.

We use full three-body scattering experiments to study the ejection of hypervelocity stars (HVSs) by massive black hole binaries (MBHBs) at the center of galaxies. Ambient stars drawn from a Maxwellian distribution unbound to the binary are expelled by the gravitational slingshot. Accurate measurements of thermally averaged hardening, mass ejection, and eccentricity growth rates (H, J, and K) for MBHBs in a fixed stellar background are obtained by numerical orbit integration from initial conditions determined by Monte Carlo techniques. Three-body interactions create a subpopulation of HVSs on nearly radial, corotating orbits, with a spatial distribution that is initially highly flattened in the inspiral plane of the MBHB, but becomes more isotropic with decreasing binary separation. The degree of anisotropy is smaller for unequal mass binaries and larger for stars with higher kick velocities. Eccentric MBHBs produce a more prominent tail of high-velocity stars and break planar symmetry, ejecting HVSs along a broad jet perpendicular to the semimajor axis. The jet two-sidedness decreases with increasing binary mass ratio, while the jet opening-angle increases with decreasing kick velocity and orbital separation. The detection of a numerous population of HVSs in the halo of the Milky Way by the next generation of large astrometric surveys like GAIA may provide a unique signature of the history, nature, and environment of the MBH at the Galactic center.

We report the discovery of a strong CIV 1548,1550 absorption system at z_abs = 5.7238 in the near-infrared spectrum (J-band) of the z_em = 6.28 QSO SDSS J1030+0524. These observations, obtained with the Infrared Spectrometer And Array Camera (ISAAC) on the European Southern Observatory Very Large Telescope (ESO VLT), demonstrate that, with modern instrumentation, QSO absorption line spectroscopy can be successfully extended to near-infrared wavelengths to probe the intergalactic medium near the end of the reionization epoch. Although the statistics of this pilot study are limited, the mass density of triply ionized carbon implied by our data is comparable to the values of Omega_CIV reported at lower redshifts. Neither the column density distribution of CIV absorbers nor its integral show significant redshift evolution over a period of time which stretches from 1 to 4.5 Gyr after the big bang, suggesting that a large fraction of intergalactic metals may already have been in place at redshifts above 6. Alternatively, the strong CIV system we have detected may be associated with outflowing, highly-ionized, gas from a foreground massive galaxy; deep imaging and spectroscopy of galaxies near the QSO sightline should be able to distinguish between these two possibilities.

Earth-mass ``microhalos'' may be the first objects to virialize in the early universe. Their ability to survive the hierarchical clustering process as substructure in the larger halos that form subsequently has implications for dark matter detection experiments. We present a large N-body simulation of early substructure in a supersymmetric cold dark matter (SUSY-CDM) scenario characterized by an exponential cutoff in the power spectrum at M_c=10^-6 Msun. The simulation resolves a 0.014 Msun parent ``SUSY'' halo at z=75 with 14 million particles. On these scales the effective index of the power spectrum approaches -3, and a range of mass scales collapse almost simultaneously. Compared to a z=0 galaxy cluster substructure within our SUSY host is less evident both in phase-space and in physical space, and it is less resistant against tidal disruption. As the universe expands by a factor of 1.3, we find that between 20 and 40 percent of well-resolved SUSY substructure is destroyed, compared to only about 1 percent in the low-redshift cluster. Nevertheless SUSY substructure is just as abundant as in z=0 galaxy clusters, i.e. the normalized mass and circular velocity functions are very similar. The dark matter self-annihilation gamma-ray luminosity from bound subhalos and other deviations from a smooth spherical configuration is at least comparable to the spherically-averaged signal in the SUSY host, and at least three times larger than the spherically-averaged signal in the cluster host. Such components must be taken into account when estimating the total cosmological extragalactic gamma-ray annihilation background.(ABRIDGED)

The question of how much light the first stars produced is fundamental to models of the Universe's development. But observations have so far failed to agree: is the answer a lot, or not very much at all?

We show that the hierarchical assembly of cold dark matter (CDM) haloes preserves the memory of the initial conditions. Using N-body cosmological simulations, we demonstrate that the present-day spatial distribution and kinematics of objects that formed within early (z > 10) protogalactic systems (old stars, satellite galaxies, globular clusters, massive black holes, etc.) depends mostly on the rarity of the peak of the primordial density field which they originally belonged to. Only for objects forming at lower redshifts the exact formation site within the progenitor halo (e.g. whether near the center or in an extended disk) becomes important. In present-day haloes, material from the rarer early peaks is more centrally concentrated and falls off more steeply with radius compared to the overall mass distribution, it has a lower velocity dispersion, moves on more radial orbits, and has a more elongated shape. Population II stars that formed within protogalactic haloes collapsing from >2.5 sigma fluctuations would follow today a r^{-3.5} density profile with a half-light radius of 17 kpc and a velocity anisotropy that increases from isotropic in the inner regions to nearly radial at the halo edge. This agrees well with the radial velocity dispersion profile of Galaxy halo stars from Battaglia et al. (2005) and with the anisotropic orbits of nearby halo stars.

We present the results of the first multi-scale N-Body+SPH simulations of merging galaxies containing central supermassive black holes (SMBHs) and having a spatial resolution of only a few parsecs. Strong gas inflows associated with equal-mass mergers produce non-axisymmetric nuclear disks with masses of order $10^9 M_{\odot}$, resolved by about $10^6$ SPH particles. Such disks have sizes of several hundred parsecs but most of their mass is concentrated within less than $50$ pc. We find that a close SMBH pair forms after the merger. The separation of the two SMBHs then shrinks further owing to dynamical friction against the predominantly gaseous background. The orbits of the SMBHs decay down to the minimum resolvable scale in a few million years for an ambient gas temperature and density typical of a region undergoing a starburst. These results suggest the initial conditions necessary for the eventual coalescence of the two holes arise naturally from the merging of two equal-mass galaxies whose structure and orbits are consistent with the predictions of the $\Lambda$CDM model. Our findings have important implications for planned gravitational wave detection experiments such as {\it LISA}.

We study the long-term evolution of massive black hole binaries (MBHBs) at the centers of galaxies using detailed full three-body scattering experiments. Stars, drawn from a distribution unbound to the binary, are ejected by the gravitational slingshot. We quantify the effect of secondary slingshots -- stars returning on small impact parameter orbits to have a second super-elastic scattering with the MBHB -- on binary separation. Even in the absence of two-body relaxation or gas dynamical processes, very unequal mass binaries of mass M=10^7 solar masses can shrink to the gravitational wave emission regime in less than a Hubble time, and are therefore a target for the planned Laser Interferometer Space Antenna (LISA). Three-body interactions create a subpopulation of hypervelocity stars on nearly radial, corotating orbits, with a spatial distribution that is initially highly flattened in the inspiral plane of the MBHB, but becomes more isotropic with decreasing binary separation. The mass ejected is ~0.7 times the binary reduced mass, and most of the stars are ejected in an initial burst lasting much less than a bulge crossing time.

We use numerical hydrodynamical simulations of early structure formation in a LCDM universe to investigate the spin temperature and 21cm brightness of the diffuse intergalactic medium (IGM) prior to the epoch of cosmic reionization, at z<20. In the absence of any radiative heating, collisions between neutral hydrogen atoms can efficiently decouple the spin temperature from the CMB only in dense minihalos and filaments. Shock heated gas shines in emission, surrounded by cooler gas visible in absorption. In the case of a warm, mostly neutral IGM, produced here by X-ray emission from an early miniquasar, the 21cm signal is strongly enhanced. Even slightly overdense filaments now shine in emission against the CMB, possibly allowing future radio arrays to probe the distribution of neutral hydrogen before reionization.

We consider the production of $^6$Li in spallation reactions by cosmic rays in order to explain the observed abundance in halo metal-poor stars. We show that heating of ambient gas by cosmic rays is an inevitable consequence of this process, and estimate the energy input required to reproduce the observed abundance of $^6$Li/H$\sim 10^{-11}$ to be of order a few hundred eV per particle. We draw attention to the possibility that this could explain the excess entropy in gas in galaxy groups and clusters. The evolution of $^6$Li and the accompanying heating of gas is calculated for structures collapsing at the present epoch with injection of cosmic rays at high redshift. We determine the energy required to explain the abundance of $^6$Li at $z \sim 2$ corresponding to the formation epoch of halo metal-poor stars, and also an increased entropy level of $\sim 300$ keV cm$^2$ necessary to explain X-ray observations of clusters. The energy budget for this process is consistent with the expected energy output of radio-loud AGNs, and the diffusion length scale of cosmic-ray protons responsible for heating is comparable to the size of regions with excess entropy. We also discuss the constraints imposed by the extragalactic gamma-ray background.

The Milky Way contains several distinct old stellar components that provide a fossil record of its formation. We can understand their spatial distribution and kinematics in a hierarchical formation scenario by associating the proto-galactic fragments envisaged by Searle and Zinn (1978) with the rare peaks able to cool gas in the cold dark matter density field collapsing at redshift z>10. We use hierarchical structure formation simulations to explore the kinematics and spatial distribution of these early star-forming structures in galaxy haloes today. Most of the proto-galaxies rapidly merge, their stellar contents and dark matter becoming smoothly distributed and forming the inner Galactic halo. The metal-poor globular clusters and old halo stars become tracers of this early evolutionary phase, centrally biased and naturally reproducing the observed steep fall off with radius. The most outlying peaks fall in late and survive to the present day as satellite galaxies. The observed radial velocity dispersion profile and the local radial velocity anisotropy of Milky Way halo stars are successfully reproduced in this model. If this epoch of structure formation coincides with a suppression of further cooling into lower sigma peaks then we can reproduce the rarity, kinematics and spatial distribution of satellite galaxies as suggested by Bullock et al. (2000). Reionisation at z=12+/-2 provides a natural solution to the missing satellites problem. Measuring the distribution of globular clusters and halo light on scales from galaxies to clusters could be used to constrain global versus local reionisation models. If reionisation occurs contemporary, our model predicts a constant frequency of blue globulars relative to the host halo mass, except for dwarf galaxies where the average relative frequencies become smaller.

We apply a new approach to quantifying galaxy morphology and identifying galaxy mergers to the rest-frame far-ultraviolet images of 82 z ~ 4 Lyman break galaxies (LBGs) and 55 1.2 < z < 1.8 emission-line galaxies in the GOODS and Ultra Deep Fields. We compare the distributions of the Gini coefficient (G), second-order moment of the brightest 20% of galaxy light (M20), and concentration (C) for high-redshift and low-redshift galaxies with average signal to noise per pixel > 2.5 and Petrosian radii > 0.3 arcsec. Ten of the 82 LBGs have M20 >= -1.1 and possess bright double or multiple nuclei, implying a major-merger fraction of star-forming galaxies ~ 10-25% at M_{FUV} < -20, depending on our incompleteness corrections. Galaxies with bulge-like morphologies (G >= 0.55, M20 < -1.6) make up ~ 30% of the z ~ 4 LBG sample, while the remaining ~ 50% have G and M20 values higher than expected for smooth bulges and disks and may be star-forming disks, minor-mergers or post-mergers. The star-forming z ~ 1.5 galaxy sample has a morphological distribution which is similar to the UDF z ~ 4 LBGs, with an identical fraction of major-merger candidates but fewer spheroids. The observed morphological distributions are roughly consistent with current hierarchical model predictions for the major-merger rates and minor-merger induced starbursts at z ~ 1.5 and ~4. We also examine the rest-frame FUV-NUV and FUV-B colors as a function of morphology and find no strong correlations at either epoch.

We investigate the environmental impact of the first active galactic nuclei that may have formed ~150 Myr after the big bang in low-mass ~10^6 Msun minihaloes. Using Enzo, an adaptive-mesh refinement cosmological hydrodynamics code, we carry out three-dimensional simulations of the radiative feedback from `miniquasars' powered by intermediate-mass black holes. We follow the non-equilibrium multispecies chemistry of primordial gas in the presence of a point source of X-ray radiation, which starts shining in a rare high-sigma peak at z=21 and emits a power-law spectrum in the 0.2-10 keV range. We find that, after one Salpeter time-scale, the miniquasar has heated up the simulation box to a volume-averaged temperature of 2800 K. The mean electron and H2 fractions are now 0.03 and 4e-5: the latter is 20 times larger than the primordial value, and will delay the buildup of a uniform UV photodissociating background. The net effect of the X-rays is to reduce gas clumping in the IGM by as much as a factor of 3. While the suppression of baryonic infall lowers the gas mass fraction at overdensities delta in the range 20-2000, enhanced molecular cooling increases the amount of dense material at delta>2000. In many haloes within the proximity of our miniquasar the H2-boosting effect of X-rays is too weak to overcome heating, and the cold and dense gas mass actually decreases. We find little evidence for an entropy floor in gas at intermediate densities preventing gas contraction and H2 formation. Overall, the radiative feedback from X-rays enhances gas cooling in lower-sigma peaks that are far away from the initial site of star formation, thus decreasing the clustering bias of the early pregalactic population, but does not appear to dramatically reverse or promote the collapse of pregalactic clouds as a whole. (abridged)

We study the detectability of primordial metal-free stars. Cosmological enrichment is a local process that takes place over an extended redshift range. While the duration of this transition depends on several unknown factors, in all cases late-forming metal-free stars are found in ~ 10^8 solar mass objects, which are just large enough to cool, but small enough to not be clustered near areas of previous star formation. We discuss the observational properties of these objects, some of which may have already been detected in ongoing Lyman-alpha surveys. If metal-free stars have masses between 140 and 260 solar masses, they are expected to end their lives as pair-production supernovae (PPSNe). Using the KEPLER code, we calculate a set of light curves for these objects. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia supernovae, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches are able to place stringent limits on the fraction of very massive stars out to redshift 2. The planned Joint Dark Energy Mission satellite will be able to extend these constraints out to redshift 6

Nonrotating, zero metallicity stars with initial masses 140 < M < 260 solar masses are expected to end their lives as pair-production supernovae (PPSNe), in which an electron-positron pair-production instability triggers explosive nuclear burning. Interest in such stars has been rekindled by recent theoretical studies that suggest primordial molecular clouds preferentially form stars with these masses. Since metal enrichment is a local process, the resulting PPSNe could occur over a broad range of redshifts, in pockets of metal-free gas. Using the implicit hydrodynamics code KEPLER, we have calculated a set of PPSN light curves that addresses the theoretical uncertainties and allows us to assess observational strategies for finding these objects at intermediate redshifts. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches may soon be able to limit the contribution of these very massive stars to < 1% of the total star formation rate density out to z=2 which already provides useful constraints for theoretical models. The planned Joint Dark Energy Mission satellite will be able to extend these limits out to z=6.

Theoretical and observational arguments suggest that the intergalactic medium (IGM) might have been polluted with metals produced by early star formation. In this scenario, Lyman-break galaxies (LBGs) at redshift z=3 are likely to be surrounded by old metal bubbles, relics of an era when the characteristic mass of galaxies was small and gas retainment more difficult. We find that pregalactic enrichment of the IGM from 10^8-10^{10} M_sun dwarf galaxies at 6<z<12 can quantitatively explain the high cross-correlation between CIV systems and LBGs observed at z=3. The reason is twofold. First, both LBGs and high-z dwarfs are biased tracer of the mass distribution and form from exceptionally high density fluctuations which are strongly clustered. Second, the action of gravity tends to increase the spatial association between metal bubbles and LBGs. Our analysis shows that, in order to match the abundance of CIV systems observed at z=3, the metal bubbles generated by high-z dwarfs must have comoving sizes of ~100 kpc. We conclude that the observed galaxy-CIV spatial association needs not to be generated by late ``superwinds'' from LBGs.

We discuss the possibility of using Supernovae as tracers of the star formation history of the Universe for the range of stellar masses $\sim 3-30$ M$_\odot$ and possibly beyond. We simulate the observations of 350 SNe, up to $z\sim 15$, made with OWL (100m) telescope.

We analyze the effect of dissipation on the orbital evolution of supermassive black holes (SMBHs) using high-resolution self-consistent gasdynamical simulations of binary equal- and unequal-mass mergers of disk galaxies. The galaxy models are consistent with the LCDM paradigm of structure formation and the simulations include the effects of radiative cooling and star formation. We find that equal-mass mergers always lead to the formation of a close SMBH pair at the center of the remnant with separations limited solely by the adopted force resolution of ~ 100 pc. Instead, the final SMBH separation in unequal-mass mergers depends sensitively on how the central structure of the merging galaxies is modified by dissipation. In the absence of dissipation, the satellite galaxy can be entirely disrupted before the merger is completed leaving its SMBH wandering at a distance too far from the center of the remnant for the formation of a close pair. In contrast, we show that gas cooling facilitates the pairing process by increasing the resilience of the companion galaxy to tidal disruption. Moreover, we demonstrate that merging disk galaxies constructed to obey the M(BH)-sigma relation, move relative to it depending on whether they undergo a dissipational or collisionless merger, regardless of the mass ratio of the merging systems. Collisionless simulations reveal that remnants tend to move away from the mean relation highlighting the role of gas-poor mergers as a possible source of scatter. In dissipational mergers, the interplay between strong gas inflows associated with the formation of massive nuclear disks and the consumption of gas by star formation provides the necessary fuel to the SMBHs and allows the merger remnants to satisfy the relation.

We compute the expected low-frequency gravitational wave signal from coalescing massive black hole (MBH) binaries at the center of galaxies. We follow the merging history of halos and associated holes via cosmological Monte Carlo realizations of the merger hierarchy from early times to the present in a LCDM cosmology. MBHs get incorporated through a series of mergers into larger and larger halos, sink to the centre owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become more massive, and form a binary system. Stellar dynamical processes dominates the orbital evolution of the binary at large separations, while gravitational wave emission takes over at small radii, causing the final coalescence of the system. We discuss the observability of inspiraling MBH binaries by a low-frequency gravitational wave experiment such as the planned Laser Interferometer Space Antenna (LISA), discriminating between resolvable sources and unresolved confusion noise. Over a 3-year observing period LISA should resolve this GWB into discrete sources, detecting ~90 individual events above a S/N=5 confidence level, while expected confusion noise is well below planned LISA capabilities.

We make a critical assessment of models that attribute the recently detected near-infrared background ``excess'' (NIRBE) to the redshifted light from Population III objects. To supply the required 25 keV/baryon at redshift 9, Pop III massive stars must form with an efficiency exceeding 30% in all ``minihalos'' with virial temperatures above a few hundred kelvins: to avoid excessive metal pollution, most of the baryons once in Pop III stars must end up in intermediate-mass black holes (IMBHs). Gas accretion onto such IMBHs must either be inhibited or lead to early miniquasars with steep UV/X-ray spectra, in order not to overproduce the present-day unresolved soft X-ray background. In the latter case (NIRBE dominated by "X-ray quiet miniquasars"), the total mass density of IMBHs at redshift 9 must be >50 times higher than the mass density of supermassive black holes observed today in the nuclei of galaxies. A stellar-dominated NIRBE is less economical energetically: more than 5% of all baryons in the universe must be processed into Population III stars. We survey various aspects of the Population III hypothesis for the NIRBE, and show that the ionizing photon budget required to account for the NIRBE is much larger than that required to explain the high electron scattering optical depth measured by the WMAP satellite.

Massive black hole binaries, with masses in the range 1E3-1E8 Msun, are expected to be the most powerful sources of gravitational radiation at mHz frequencies, and hence are among the primary targets for the planned Laser Interferometer Space Antenna (LISA). We extend and refine our previous analysis (Sesana et al. 2004), detailing the gravitational wave signal expected from a cosmological population of massive black hole binaries. As done in our previous paper, we follow the merger history of dark matter halos, the dynamics of the massive black holes they host, and their growth via gas accretion and binary coalescences in a LCDM cosmology. Stellar dynamical processes dominates the orbital evolution of black hole binaries at large separations, while gravitational wave emission takes over at small radii, causing the final coalescence of the pairs. We show that the GW signal from this population, in a 3 year LISA observation, will be resolved into approx 90 discrete events with S/N>5, among which approx. 35 will be observed above threshold until coalescence. These "merging events" involve relatively massive binaries, M=10E5 Msun, in the redshift range 2<z<6. The remaining approx. 55 events come from higher redshift, less massive binaries (M=5E3 Msun at z>6) and, although their S/N integrated over the duration of the observation can be substantial, the final coalescence phase is at too high frequency to be directly observable by space-based interferometers such as LISA. LISA will be able to detect a fraction approx. 90% of all the coalescences of massive black hole binaries occurring at z=5. The residual confusion noise from unresolved massive black hole binaries is expected to be at least an order of magnitude below the estimated stochastic noise.

I review scenarios for the assembly of supermassive black holes (MBHs) at the center of galaxies that trace their hierarchical build-up far up in the dark matter halo "merger tree". Monte Carlo realizations of the merger hierarchy in a LCDM cosmology, coupled to semi-analytical recipes, are a powerful tool to follow the merger history of halos and the dynamics and growth of the MBHs they host. X-ray photons from miniquasars powered by intermediate-mass "seed" holes may permeate the universe more uniformly than EUV radiation, make the low-density diffuse intergalactic medium warm and weakly ionized prior to the epoch of reionization breakthrough, and set an entropy floor. The spin distribution of MBHs is determined by gas accretion, and is predicted to be heavily skewed towards fast-rotating Kerr holes, to be already in place at early epochs, and not to change significantly below redshift 5. Decaying MBH binaries may shape the innermost central regions of galaxies and should be detected in significant numbers by LISA.

We study the expected distribution of massive black hole (MBH) spins and its evolution with cosmic time in the context of hierarchical galaxy formation theories. Our model uses Monte Carlo realizations of the merger hierarchy in a LCDM cosmology, coupled to semi-analytical recipes, to follow the merger history of dark matter halos, the dynamics of the MBHs they host, and their growth via gas accretion and binary coalescences. The coalescence of comparable mass holes increases the spin of MBHs, while the capture of smaller companions in randomly-oriented orbits acts to spin holes down. We find that, given the distribution of MBH binary mass ratios in hierarchical models, binary coalescences alone do not lead to a systematic spin-up or spin-down of MBHs with time: the spin distribution retains memory of its initial conditions. By contrast, because of the Bardeen-Petterson effect, gas accretion via a thin disk tends to spin holes up even if the direction of the spin axis changes randomly in time. In our models, accretion dominates over black hole captures and efficiently spins holes up. The spin distribution is heavily skewed towards fast-rotating Kerr holes, is already in place at early epochs, and does not change much below redshift 5. If accretion is via a thin disk, about 70% of all MBHs are maximally rotating and have radiative efficiencies approaching 30% (assuming a "standard'' spin-efficiency conversion). Even in the conservative case where accretion is via a geometrically thick disk, about 80% of all MBHs have spin parameters a/m > 0.8 and accretion efficiencies > 12%. Rapidly spinning holes with high radiative efficiencies may satisfy constraints based on comparing the local MBH mass density with the mass density inferred from luminous quasars (Soltan's argument).

We use a sample of 42 supernovae detected with the Advanced Camera for Surveys on-board the Hubble Space Telescope as part of the Great Observatories Origins Deep Survey to measure the rate of core collapse supernovae to z~0.7 and type Ia supernovae to z~1.6. This significantly increases the redshift range where supernova rates have been estimated from observations. The rate of core collapse supernovae can be used as an independent probe of the cosmic star formation rate. Based on the observations of 17 core collapse supernovae, we measure an increase in the core collapse supernova rate by a factor of 1.6 in the range 0.3<z<0.7, and an overall increase by a factor of 7 to z~0.7 in comparison to the local core collapse supernova rate. The increase in the rate in this redshift range in consistent with recent measurements of the star formation rate derived from UV-luminosity densities and IR datasets. Based on 25 type Ia supernovae, we find a SN Ia rate that is a factor 3-5 higher at z~1 compared to earlier estimates at lower redshifts (z<0.5), implying that the type Ia supernova rate traces a higher star formation rate at redshifts z>1 compared to low redshift. At higher redshift (z>1), we find a suggested decrease in the type Ia rate with redshift. This evolution of the Ia rate with redshift is consistent with a type Ia progenitor model where there is a substantial delay between the formation of the progenitor star and the explosion of the supernova. Assuming that the type Ia progenitor stars have initial main sequence masses 3-8 M_Sun, we find that 5-7% of the available progenitors explode as type Ia supernovae.

There is now compelling evidence of a link between long-duration gamma-ray bursts (GRBs) and Type Ib/c supernovae (SNe). These core-collapse explosions are conjectured to radiate an anisotropic, beamed component associated with a decelerating, relativistic outflow and an unbeamed, isotropic component associated with the slowly expanding stellar debris. The anisotropic emission remains at a very low level until the Doppler cone of the beam intersects the observer's line of sight, making off-axis GRB jets directly detectable only at long wavelengths and late times. Circumstellar material, however, will Compton scatter the prompt gamma-ray and afterglow radiation flux and give rise to a reflection echo. We show that the Compton echo of a misaligned GRB carries an X-ray luminosity that may exceed by many orders of magnitude that produced by the underlying subrelativistic SN during the first few weeks. Bright scattering echoes may therefore provide a means for detecting a population of misaligned GRBs associated with nearby Type Ib/c SNe and yield crucial information on the environment surrounding a massive star at the time of its death. The question of whether the interpretation of GRB980425 as an ordinary GRB observed off-axis is consistent with the lack of an X-ray echo is addressed, along with the constraints derived on the possible existence of misaligned GRB jets in SN1993J, SN1994I, SN1999em, and SN2002ap.

We compute the expected gravitational wave signal from coalescing massive black hole (MBH) binaries at the center of galaxies in a hierarchical structure formation scenario in which seed holes of intermediate mass form far up in the dark halo merger tree. The merger history of DM halos and MBHs is followed from z=20 to the present in a LCDM cosmology. MBHs get incorporated through halo mergers into larger and larger structures, sink to the center owing to dynamical friction against the DM background, accrete cold material in the merger remnant, and form MBH binary systems. Stellar dynamical interactions cause the hardening of the binary at large separations, while gravitational wave emission takes over at small radii and leads to the final coalescence of the pair. The integrated emission from inspiraling MBH binaries results in a gravitational wave background (GWB). The characteristic strain spectrum has the standard h_c(f)\propto f^{-2/3} behavior only in the range 1E-9<f<1E-6 Hz. At lower frequencies the orbital decay of MBH binaries is driven by the ejection of background stars, and h_c(f) \propto f. At higher frequencies, f>1E-6 Hz, the strain amplitude is shaped by the convolution of last stable circular orbit emission. We discuss the observability of inspiraling MBH binaries by the planned LISA. Over a 3-year observing period LISA should resolve this GWB into discrete sources, detecting ~60 (~250) individual events above a S/N=5 (S/N=1) confidence level. (Abridged)

We present two new non-parametric methods for quantifying galaxy morphology: the relative distribution of the galaxy pixel flux values (the Gini coefficient or G) and the second-order moment of the brightest 20% of the galaxy's flux (M20). We test the robustness of G and M20 to decreasing signal-to-noise and spatial resolution, and find that both measures are reliable to within 10% at average signal-to-noise per pixel greater than 3 and resolutions better than 1000 pc and 500 pc, respectively. We have measured G and M20, as well as concentration (C), asymmetry (A), and clumpiness (S) in the rest-frame near-ultraviolet/optical wavelengths for 150 bright local "normal" Hubble type galaxies (E-Sd) galaxies and 104 0.05 < z < 0.25 ultra-luminous infrared galaxies (ULIRGs).We find that most local galaxies follow a tight sequence in G-M20-C, where early-types have high G and C and low M20 and late-type spirals have lower G and C and higher M20. The majority of ULIRGs lie above the normal galaxy G-M20 sequence, due to their high G and M20 values. Their high Gini coefficients arise from very bright nuclei, while the high second-order moments are produced by multiple nuclei and bright tidal tails. All of these features are signatures of recent and on-going mergers and interactions. We also find that in combination with A and S, G is more effective than C at distinguishing ULIRGs from the "normal" Hubble-types. Finally, we measure the morphologies of 45 1.7 < z < 3.8 galaxies from HST NICMOS observations of the Hubble Deep Field North. We find that many of the z $\sim$ 2 galaxies possess G and A higher than expected from degraded images of local elliptical and spiral galaxies, and have morphologies more like low-redshift single nucleus ULIRGs.

The coalescence of massive black hole (MBH) binaries following galaxy mergers is one of the main sources of low-frequency gravitational radiation. A higher-order relativistic phenomenon, the recoil as a result of the non-zero net linear momentum carried away by gravitational waves, may have interesting consequences for the demography of MBHs at the centers of galaxies. We study the dynamics of recoiling MBHs and its observational consequences. The ``gravitational rocket'' may: i) deplete MBHs from late-type spirals, dwarf galaxies, and stellar clusters; ii) produce off-nuclear quasars, including unusual radio morphologies during the recoil of a radio-loud source; and iii) give rise to a population of interstellar and intergalactic MBHs.

Motivated by the recent detection by WMAP of a large optical depth to Thomson scattering -- implying a very early reionization epoch -- we assess a scenario where the universe was reionized by `miniquasars' powered by intermediate-mass black holes (IMBHs), the remnants of the first generation of massive stars. Pregalactic IMBHs form within minihalos above the cosmological Jeans mass collapsing at z=24, get incorporated through mergers into larger and larger systems, sink to the center owing to dynamical friction, and accrete cold material. The merger history of dark halos and associated IMBHs is followed by Monte Carlo realizations of the merger hierarchy in a LCDM cosmology. While seed IMBHs that are as rare as the 3.5-sigma peaks of the primordial density field evolve largely in isolation, a significant number of black hole binary systems will form if IMBHs populate the more numerous 3-sigma peaks instead. In the case of rapid binary coalescence a fraction of IMBHs will be displaced from galaxy centers and ejected into the IGM by the `gravitational rocket' effect, rather than accrete and shine as miniquasars. We show that, under a number of plausible assumptions for the amount of gas accreted onto IMBHs and their emission spectrum, miniquasars powered by IMBHs may be responsible for cosmological reionization at z~15. Reionization by miniquasars with a hard spectrum may be more `economical' than stellar reionization, as soft X-rays escape more easily from the dense sites of star formation and travel further than EUV radiation. Energetic photons will make the diffuse IGM warm and weakly ionized prior to the epoch of reionization breakthrough, set an entropy floor, and reduce gas clumping. Future 21 cm observations may detect a preheated, weakly-ionized IGM in emission against the CMB. (abridged)

The integrated colors of distant galaxies provide a means for interpreting the properties of their stellar content. Here, we use rest-frame UV-to-optical colors to constrain the spectral-energy distributions and stellar populations of color-selected, B-dropout galaxies at z ~ 4 in the Great Observatories Origins Deep Survey. We combine the ACS data with ground-based near-infrared images, which extend the coverage of galaxies at z ~ 4 to the rest-frame B-band. We observe a color-magnitude trend in the rest-frame m(UV) - B versus B diagram for the z ~ 4 galaxies that has a fairly well-defined "blue-envelope", and is strikingly similar to that of color-selected, U-dropout galaxies at z ~ 3. We also find that although the co-moving luminosity density at rest-frame UV wavelengths (1600 Angstroms) is roughly comparable at z ~ 3 and z ~ 4, the luminosity density at rest-frame optical wavelengths increases by about one-third from z ~ 4 to z ~ 3. Although the star-formation histories of individual galaxies may involve complex and stochastic events, the evolution in the global luminosity density of the UV-bright galaxy population corresponds to an average star-formation history with a star-formation rate that is constant or increasing over these redshifts. This suggests that the evolution in the luminosity density corresponds to an increase in the stellar-mass density of more than 33%.

We investigate a hierarchical structure formation scenario describing the evolution of a Super Massive Black Holes (SMBHs) population. The seeds of the local SMBHs are assumed to be 'pregalactic' black holes, remnants of the first POPIII stars. As these pregalactic holes become incorporated through a series of mergers into larger and larger halos, they sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, form a binary system, and eventually coalesce. A simple model in which the damage done to a stellar cusps by decaying BH pairs is cumulative is able to reproduce the observed scaling relation between galaxy luminosity and core size. An accretion model connecting quasar activity with major mergers and the observed BH mass-velocity dispersion correlation reproduces remarkably well the observed luminosity function of optically-selected quasars in the redshift range 1<z<5. We finally asses the potential observability of the gravitational wave background generated by the cosmic evolution of SMBH binaries by the planned space-born interferometer LISA.

We study lens statistics in flat, low-density universes with different equations of state $w=p_Q/\rho_Q$ for the dark energy component. Dark energy modifies the distance-redshift relation and the mass function of dark matter halos leading to changes in the lensing optical depth as a function of image separation. Those effects must, however, be distinguished from effects associated with the structure of dark matter halos. Baryonic cooling causes galaxy-mass halos to have different central density profiles than group- and cluster-mass halos, which causes the distribution of normal arcsecond-scale lenses to differ from the distribution of ``wide-separation'' ($\Delta\theta \gtrsim 4\arcsec$) lenses. Fortunately, the various parameters related to cosmology and halo structure have very different effects on the overall image separation distribution: (1) the abundance of wide-separation lenses is exremely sensitive (by orders of magnitude) to the distribution of ``concentration'' parameters for massive halos modeled with the Navarro-Frenk-White profile; (2) the transition between normal and wide-separation lenses depends mainly on the mass scale where baryonic cooling ceases to be efficient; and (3) dark energy has effects at all image separation scales. While current lens samples cannot usefully constrain all of the parameters, ongoing and future imaging surveys should discover hundreds or thousands of lenses and make it possible to disentangle the various effects and constrain all of the parameters simultaneously. (abridged)

We have measured the rest--frame lambda~1500 Ang comoving specific luminosity density of star--forming galaxies at redshift 3.5<z<6.5 from deep images taken with the Hubble Space Telescope (HST)and the Advanced Camera for Surveys (ACS), obtained as part of the Great Observatories Origins Deep Survey (GOODS). We used color selection criteria to construct samples of star--forming galaxies at redshifts z~4, 5 and 6, identified by the signature of the 912 Ang Lyman continuum discontinuity and Lyman-alpha forest blanketing in their rest--frame UV colors (Lyman--break galaxies). The ACS samples cover ~0.09 square degree, and are also relatively deep, reaching between 0.2 and 0.5 L_3^*, depending on the redshift, where $L_3^*$ is the characteristic UV luminosity of Lyman--break galaxies at z~3. The specific luminosity density of Lyman--break galaxies appears to be nearly constant with redshift from z~3 to z~6, although the measure at z~6 remains relatively uncertain, because it depends on the accurate estimate of the faint counts of the z~6 sample. If Lyman--break galaxies are fair tracers of the cosmic star formation activity, our results suggest that at z~6 the universe was already producing stars as vigorously as it did near its maximum several Gyr later, at 1<~z<~3. Thus, the onset of large--scale star formation in the universe is to be sought at around z~6 or higher, namely at less than ~7% of the current cosmic age.

Hubble Space Telescope images of high-redshift galaxies selected via color and photometric redshifts are used to examine the size and axial-ratio distribution of galaxies as a function of redshift at lookback times t > 8 Gyr. These parameters are measured at rest-frame UV wavelengths (1200-2000 Angstroms) on images with a rest-frame resolution of less than 0.8 kpc. Galaxy radii are found to scale with redshift approximately as the inverse of the Hubble parameter H(z). This is in accord with the theoretical expectation that the typical sizes of the luminous parts of galaxies should track the expected evolution in the virial radius of dark-matter halos. The mean ratio of semi-major to semi-minor axis for a bright well-resolved sample of galaxies at z ~ 4 is b/a = 0.65, suggesting that these Lyman break galaxies are not drawn from a spheroidal population. However the median concentration index of this sample is C = 3.5, which is closer to the typical concentration indices, C ~ 4, of nearby elliptical galaxies than to the values, C<2 for local disk galaxies of type Sb and later.

We use numerical simulations of hydrogen reionization by stellar sources in the context of LCDM cosmogonies to investigate the 21 (1+z) cm radio signal expected from the diffuse intergalactic medium (IGM) prior to the epoch of reionization breakthrough at redshift z_ion. Two reionization scenarios are analyzed in details: an `early reionization' case with z_ion~13, consistent with the recent discovery by the WMAP satellite of a large optical depth to Thomson scattering, and a `late reionization' case with z_ion~8. It is a generic prediction of these models that the background of Lyalpha photons produced by the early generation of stars which ultimately ionize the universe will be sufficiently intense to make intergalactic neutral hydrogen visible against the cosmic microwave background during the `grey age', i.e. z_ion<z<20. Depending on the redshift of reionization breakthrough, broad-beam observations at frequencies <150 MHz (below 100 MHz for the `early reionization' scenario) with the next generation of radio telescopes should reveal angular fluctuations in the sky brightness temperature in the range 5-20 mK (1 sigma) on scales below 5 arcmin.

Winds from pregalactic starbursts and 'miniquasars' may pollute the IGM with metals and raise its temperature to a high adiabat, and so inhibit the formation of early galaxies. We compute the thermal history of the IGM when it experiences a period of rapid, homogeneous "preheating" at high redshifts. Measurements of the temperature of the Lyamn-alpha forest at z~3 constrain the redshift and energy of preheating, and rule out models that preheat too late or to too high a temperature. We predict galaxy luminosity functions in preheated universes. The results depend crucially on whether the baryonic smoothing scale in the IGM is computed globally, or in a local, density-dependent fashion. Using a globally averaged smoothing scale, we find that models with excessive preheating produce too few L_* and fainter galaxies, and are therefore inconsistent with observational data. A density-dependent smoothing scale requires more energetic preheating to achieve the same degree of suppression in the faint-end slope. All models, however, appear unable to explain the sharp cut-off in the luminosity function at bright magnitudes. Supernova-driven preheating scenarios tend to raise the mean metallicity of the universe well above the minimum levels observed in the Lyman-alpha clouds. We find that ionizing photon escape fractions must be significantly higher than 10% in order to explain the low inferred HI fraction at z~6. While early preheating causes strong suppression of dwarf galaxy formation we show that it is not able to reproduce the observed abundance of satellite galaxies in the Local Group in detail.

In order to search for metals in the Lyman alpha forest at redshifts z > 4, we have obtained spectra of high S/N and resolution of three QSOs at z > 5.4 discovered by the Sloan Digital Sky Survey. These data allow us to probe to metal enrichment of the intergalactic medium at early times with higher sensitivity than previous studies. We find 16 C IV absorption systems with column densities log N(C IV) = 12.50 - 13.98 over a total redshift path Delta X = 3.29. In the redshift interval z = 4.5-5.0, where our statistics are most reliable, we deduce a comoving mass density of C IV ions Omega(C IV) = (4.3 +/- 2.5) x 10(-8) (90% confidence limits) for absorption systems with log N(C IV) > 13.0 (for an Einstein-de Sitter cosmology with h = 0.65). This value of Omega(C IV) is entirely consistent with those measured at z < 4; we confirm the earlier finding by Songaila (2001) that neither the column density distribution of C IV absorbers nor its integral show significant redshift evolution over a period of time which stretches from 1.25 to 4.5 Gyr after the big bang. This somewhat surprising conclusion may be an indication that the intergalactic medium was enriched in metals at redshifts much greater than 5, perhaps by the sources responsible for its reionization. Alternatively, the C IV systems we see may be associated with outflows from massive star-forming galaxies at later times, while the truly intergalactic metals may reside in regions of the Lyman alpha forest of lower density than those probed up to now.

We investigate a hierarchical structure formation scenario in which galaxy stellar cores are created from the binding energy liberated by shrinking supermassive black hole (SMBH) binaries. The binary orbital decay heats the surrounding stars, eroding a preexisting 1/r^2 stellar cusp. We follow the merger history of dark matter halos and associated SMBHs via cosmological Monte Carlo realizations of the merger hierarchy from early times to the present in a LCDM cosmology. Massive black holes get incorporated through a series of mergers into larger and larger halos, sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, and form a binary system. Stellar dynamical processes drive the binary to harden and eventually coalesce. A simple scheme is applied in which the loss cone is constantly refilled and a constant density core forms due to the ejection of stellar mass. We find that a model in which the effect of the hierarchy of SMBH interactions is cumulative and cores are preserved during galaxy mergers produces at the present epoch a correlation between the `mass deficit' (the mass needed to bring a flat inner density profile to a 1/r^2 cusp) and the mass of the nuclear SMBH, with a normalization and slope comparable to the observed relation. Models in which the mass displaced by the SMBH binary is replenished after every major galaxy merger appear instead to underestimate the mass deficit observed in `core' galaxies.

The development of primordial inhomogeneities into the non-linear regime and the formation of the first astrophysical objects within dark matter halos mark the transition from a simple, neutral, cooling universe -- described by just a few parameters -- to a messy ionized one -- the realm of radiative, hydrodynamic, and star formation processes. The recent measurement by the WMAP satellite of a large optical depth to electron scattering implies that this transition must have begun very early, and that the universe was reionized at redshift z_ion=17\pm 5. It is an early generation of extremely metal-poor massive stars and/or `seed' accreting black holes in subgalactic halos that may have generated the ultraviolet radiation and mechanical energy that reheated and reionized most of the hydrogen in the cosmos. The detailed thermal, ionization, and chemical enrichment history of the universe during the crucial formative stages around z=10-20 depends on the power-spectrum of density fluctuations on small scales, the stellar initial mass function and star formation efficiency, a complex network of poorly understood `feedback' mechanisms, and remains one of the crucial missing links in galaxy formation and evolution studies.

In popular cold dark matter cosmological scenarios, stars may have first appeared in significant numbers around a redshift of 10 or so, as the gas within protogalactic halos with virial temperatures in excess of 20,000 K (corresponding to masses comparable to those of present-day dwarf ellipticals) cooled rapidly due to atomic processes and fragmented. It is this `second generation' of subgalactic stellar systems, aided perhaps by an early population of accreting black holes in their nuclei, which may have generated the ultraviolet radiation and mechanical energy that ended the cosmic ``dark ages'' and reheated and reionized most of the hydrogen in the universe by a redshift of 6. The detailed history of the universe during and soon after these crucial formative stages depends on the power-spectrum of density fluctuations on small scales and on a complex network of poorly understood `feedback' mechanisms, and is one of the missing link in galaxy formation and evolution studies. The astrophysics of the epoch of first light is recorded in the thermal state, ionization degree, and chemical composition of the intergalactic medium, the main repository of baryons at high redshifts.

The development of primordial inhomogeneities into the non-linear regime and the formation of the first bound objects mark the transition from a simple cooling universe -- described by just a few parameters -- to a very messy hot one -- the realm of radiative, hydrodynamic and star formation processes. In popular cold dark matter cosmological scenarios, stars and quasars may have first appeared in significant numbers around a redshift of 10 or so, as the gas within subgalactic halos cooled rapidly due to atomic processes and fragmented. This early generation of sources generated the ultraviolet radiation and kinetic energy that ended the cosmic ``dark ages'' and reheated and reionized most of the hydrogen in the cosmos by a redshift of 6. The detailed thermal history of the intergalactic medium -- the main repository of baryons at high redshift -- during and soon after these crucial formative stages depends on the power-spectrum of density fluctuations on small scales and on a complex network of poorly understood `feedback' mechanisms, and is one of the missing link in galaxy formation and evolution studies.

We assess models for the assembly of supermassive black holes (SMBHs) at the center of galaxies that trace their hierarchical build-up far up in the dark halo `merger tree'. We assume that the first `seed' black holes (BHs) formed in (mini)halos collapsing at z=20 from high-sigma density fluctuations. As these pregalactic holes become incorporated through a series of mergers into larger and larger halos, they sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, form a binary system, and eventually coalesce. The merger history of dark matter halos and associated BHs is followed by cosmological Monte Carlo realizations of the merger hierarchy. A simple model, where quasar activity is driven by major mergers and SMBHs accrete at the Eddington rate a mass that scales with the fifth power of the velocity dispersion, is shown to reproduce the optical LF of quasars in the redshift range 1<z<4. Binary and triple BH interactions are followed in our merger tree. The assumptions underlying our scenario lead to the prediction of a population of massive BHs wandering in galaxy halos and the intergalactic medium at the present epoch, and contributing <10% to the total BH mass density. At all epochs the fraction of binary SMBHs in galaxy nuclei is of order 10%, while the fraction of binary quasars is less than 0.3%

In popular cold dark matter cosmological scenarios, stars may have first appeared in significant numbers around a redshift of 10 or so, as the gas within protogalactic halos with virial temperatures in excess of 20,000 K (corresponding to masses comparable to those of present--day dwarf ellipticals) cooled rapidly due to atomic processes and fragmented. It is this `second generation' of subgalactic stellar systems, aided perhaps by an early population of accreting black holes in their nuclei, which may have generated the ultraviolet radiation and mechanical energy that ended the cosmic ``dark ages'' and reheated and reionized most of the hydrogen in the cosmos by a redshift of 6. The detailed history of the universe during and soon after these crucial formative stages depends on the power--spectrum of density fluctuations on small scales and on a complex network of poorly understood `feedback' mechanisms, and is one of the missing link in galaxy formation and evolution studies. The astrophysics of the epoch of first light is recorded in the thermal state, ionization degree, and chemical composition of the intergalactic medium, the main repository of baryons at high redshifts.

Supernova-driven outflows from early galaxies may have had a large impact on the kinetic and chemical structure of the intergalactic medium (IGM). We use three-dimensional Monte Carlo cosmological realizations of a simple linear peaks model to track the time evolution of such metal-enriched outflows and their feedback on galaxy formation. We find that at most 30% of the IGM by volume is enriched to values above 10^-3 solar in models that only include objects that cool by atomic transitions. The majority of enrichment occurs relatively early (5 < z < 12) and resulting in a mass-averaged cosmological metallicity between 10^-3 and 10^-1.5 solar. The inclusion of Population III objects that cool through H2 line emission has only a minor impact on these results: increasing the mean metallicity and filling factor by at most a factor of 1.4, and moving the dawn of the enrichment epoch to a redshift of approximately 14 at the earliest. Thus enrichment by outflowing galaxies is likely to have been incomplete and inhomogeneous, biased to the areas near the starbursting galaxies themselves. Models with a 10% star formation efficiency can satisfactorily reproduce the nearly constant (2 < z < 5, Z approximately 3.5 x 10^-4 solar) metallicity of the low column density Ly-alpha forest derived by Songaila (2001), an effect of the decreasing efficiency of metal loss from larger galaxies. Finally, we show that IGM enrichment is intimately tied to the ram-pressure stripping of baryons from neighboring perturbations. This results in the suppression of at least 20% of the dwarf galaxies in the mass range 10^8.5 to 10^9.5 solar, in all models with filling factors greater than 2%, and an overall suppression of approximately 50% of dwarf galaxies in the most observationally-favored model.

We present results from three-dimensional numerical simulations of the dynamics of SN-driven bubbles as they propagate through and escape the grasp of subgalactic halos with masses M ~ 10^8 Msun at redshift z=9. The hydrodynamic simulations use a nested grid method to follow the evolution of explosive multi-SN events operating on the characteristic timescale of a few x 10^7 yr, the lifetime of massive stars. The results confirm that, if the star formation efficiency of subgalactic halos is close to 10%, a significant fraction of the halo gas will be lifted out of the potential well (`blow-away'), shock the intergalactic medium, and pollute it with metal-enriched material, a scenario recently advocated by Madau, Ferrara & Rees (2001). The volume filling factor of the ejecta is close to unity. Depending on the stellar distribution, we find that less than 30% of the available SN energy is converted into kinetic energy of the blown away material, the remainder being radiated away. It appears that mechanical feedback is less efficient than expected from simple energetic arguments, as off-nuclear SN explosions drive inward-propagating shocks that tend to collect and pile up cold gas in the central regions of the host halo. Low-mass galaxies at early epochs may survive multiple SN events and continue forming stars.

Recent studies of heavy r-process elements in low [Fe/H] halo stars have suggested that an initial population of metal-free very massive stars (VMSs) may be required to provide early Fe enrichment without coproducing heavy r nuclei. We find similar abundance trends in alpha-elements (which should be copiously produced by VMSs), but not in other elements such as carbon (which should not), in agreement with this hypothesis. We then combine the corresponding level of prompt initial enrichment with models of VMS nucleosynthetic yields and spectra to estimate the corresponding ionizing fluxes. The result suggests that there may have been enough VMS activity to reionize the universe. The unusually hard spectrum of VMSs would imply a different reionization history from canonical models. HeII could have been reionized at high redshift, only to recombine as a subsequent generation of stars formed with a ``normal'' initial mass function.

Recent observations suggest that long-duration gamma-ray bursts (GRBs) and their afterglows are produced by highly relativistic jets emitted in core-collapse explosions. The pre-explosive ambient medium provides a natural test for the most likely progenitors of GRBs. Those stars that shed their envelopes most readily have short jet crossing times and are more likely to produce a GRB. We construct a simple computational scheme to explore the expected contribution of the presupernova ejecta of single Wolf-Rayet (WR) stars to the circumstellar environment. Using detailed stellar tracks for the evolution of massive stars, we discuss the effects that the initial main sequence mass, metallicity, rotation and membership in a binary system have on the ambient medium. We extend the theory of GRB afterglows in winds to consider the effect of the relativistic fireball propagating through the WR ejecta. We make specific predictions for the interaction of the relativistic blast wave with the density bumps that arise when the progenitor star rapidly loses a large fraction of its initial mass or when the ejected wind interacts with the external medium and decelerates. A re-brightening of the afterglow with a spectrum redder than the typical synchrotron spectrum (as seen in GRB 970508, GRB 980326 and GRB 000911) is predicted. We also calculate the luminosity of the reflected echo that arises when circumstellar material Compton-scatters the prompt radiation and examine the spectral signatures expected from the interaction of the GRB afterglow with the ejected medium.

We show how observations of multiply-imaged quasars at high redshift can be used as a probe of dark matter clumps (subhalos with masses ~ 10^9 solar masses) within the virialized extent of more massive lensing halos. A large abundance of such satellites is predicted by numerical simulations of galaxy formation in cold dark matter (CDM) cosmogonies. Small-scale structure within galaxy halos affects the flux ratios of the images without appreciably changing their positions. We use numerical simulations to quantify the effect of dark matter substructure on the distribution of magnifications, and find that the magnification ratio of a typical image pair will deviate significantly from the value predicted by a smooth lensing potential if, near the Einstein radius, only a few percent of the lens surface density is contained in subhalos. The angular size of the continuum source dictates the range of subclump masses that can have a detectable effect: to avoid confusion with gravitational microlensing caused by stars in the lens galaxy, the background source must be larger than the optical continuum-emitting region of a QSO. We also find that substructure will cause distortions to images on milli-arcsecond scales and bias the distribution of QSO magnification ratios -- two other possible methods of detection.

In this paper, I will describe the features of the numerical code CUBA, aimed at the solution of the radiative transfer equation in a cosmological context. CUBA will be soon available for public use at the URL http://pitto.mib.infn.it/~haardt/, allowing for several user-supplied input parameters, such as favourite cosmology, luminosity functions, Type II object evolution, stellar spectra, and many others. I will also present some new results of the UV/X-ray cosmic background as produced by the observed populations of QSOs and star forming galaxies, updating and extending our previous works. The background evolution is complemented with a number of derived quantities such as the ionization and thermal state of the IGM, the HeII opacity, the HI and HeII ionization rates, and the HI, HeII and Compton heating rates.

Recent numerical simulations of the fragmentation of primordial molecular clouds in hierarchical cosmogonies have suggested that the very first stars (the so-called Population III) may have been rather massive. Here we point out that a numerous population of massive black holes (MBHs) -- with masses intermediate between those of stellar and supermassive holes -- may be the endproduct of such an episode of pregalactic star formation. If only one MBH with m > 150 msun formed in each of the `minihalos' collapsing at z=20 from 3-sigma fluctuations, then the mass density of Pop III MBHs would be comparable to that of the supermassive variety observed in the nuclei of galaxies. Since they form in high-sigma rare density peaks, relic MBHs are predicted to cluster in the bulges of present-day galaxies as they become incorporated through a series of mergers into larger and larger systems. Dynamical friction would cause more than 50 (m/150 msun)^{1/2} such objects to sink towards the center. The presence of a small cluster of MBHs in galaxy nuclei may have several interesting consequences associated with tidal captures of ordinary stars (likely followed by disruption), MBH capture by the central supermassive black hole, gravitational wave radiation from such coalescences. Accreting pregalactic MBHs may be detectable as ultra-luminous, off-nuclear X-ray sources.

We assess supernova (SN)-driven pregalactic outflows as a mechanism for distributing the product of stellar nucleosynthesis over large cosmological volumes prior to the reionization epoch. SN ejecta will escape the grasp of halos with virial temperatures T_vir>10^{4.3} K (corresponding to masses M>10^8 h^{-1} M_sun at redshift z=9 when they collapse from 2-sigma fluctuations) if rapid cooling can take place, and a significant fraction of their baryonic mass is converted into stars over a dynamical timescale. We study the evolution of SN-driven bubbles as they blow out from subgalactic halos and propagate into the intergalactic medium (IGM), and show that to lift the halo gas out of the potential well the energy injection must continue at least until blow-away occurs. The collective explosive output of about ten thousands SNe per halo at these early epochs could pollute vast regions of intergalactic space to a mean metallicity <Z>=Omega_Z/Omega_b>0.003 (comparable to the levels observed in the Lyman-alpha forest at z=3) without much perturbing the IGM hydrodynamically, i.e. producing large variations of the baryons relative to the dark matter. The volume filling factor of the ejecta is higher than 20% if the star formation efficiency is of order 10%. When the filling factor of the ejecta becomes significant, enriched material will typically be at a higher adiabat than expected from photoionization.

The statistics of wide-separation (6'' < theta < 15'') gravitational lenses constrain the amount of mass in the cores of dark matter halos on group and cluster mass scales. For a family of halo models with a central cusp $\rho \propto r^{-alpha}$ (1.0 <= alpha <= 1.9), the lack of wide-separation lenses in the large Cosmic Lens All-Sky Survey yields an upper limit on the fraction of the halo mass that is contained within ~4% of the virial radius, f_{core} < 0.023 (95% confidence level, LCDM). This limit offers an important test of the cold dark matter paradigm. While the halo profiles derived from numerical simulations appear to be consistent with this upper limit, larger surveys currently underway such as the 2dF and SDSS should detect wide-separation lenses and thus provide a measurement of the core mass fraction in massive dark matter halos.

The optical/far-IR extragalactic background light (EBL) from both resolved and unresolved extragalactic sources is an indicator of the total luminosity of cosmic structures, as the cumulative emission from young and evolved galactic systems, as well as from active galactic nuclei (AGNs), is recorded in this radiation. This is a brief review of some of the implications of the observed brightness of the night sky for the stellar mass density and average metallicity of the universe today, and of the possible contribution of MACHO progenitors and QSOs to the EBL. Assuming a Salpeter initial mass function with a cutoff below 0.6 M_sun, a lower limit of Omega_{g+s}h^2>0.0015 I_{60} can be derived to the visible (recycled gas + stars) mass density required to generate an EBL at a level of I_{EBL}=60 I_{60} nW/m^2/sr: our latest, best-guess estimate is Omega_{g+s}h^2 ~ 0.0023 I_{60}, which implies a mean metallicity at the present-epoch of y_Z Omega_{g+s}/Omega_b ~ 0.2 Z_sun. If massive dark halos around spiral galaxies are partially composed of faint, old white dwarfs, i.e. if a non-negligible fraction (a few percent) of the nucleosynthetic baryons is locked in the remnants (MACHOs) of intermediate-mass stars forming at very high redshifts, then the bright early phases of such halos should contribute significantly to the observed EBL. Assuming a standard black hole accretion model for quasar activity and using recent observations of the quasar population and new synthesis models for the cosmic X-ray background, we estimate a present mass density of QSO remnants of rho_{BH} ~ 3e5 M_sun/Mpc^3 for a 10% efficiency of accreted mass-to-radiation conversion. The quasar contribution to the brightness of the night sky is I_QSO ~ 2 nW/m^2/sr.

At the earliest epochs of structure formation in cold dark matter (CDM) cosmologies, the smallest nonlinear objects are the numerous small halos that condense with virial temperatures below 10,000 K. Such ``minihalos'' are not yet resolved in large-scale three-dimensional cosmological simulations. Here we employ a semi-analytic method, combined with three-dimensional simulations of individual minihalos, to examine their importance during cosmological reionization. We show that, depending on when reionization takes place, they potentially play an important role as sinks of ionizing radiation. If reionization occurs at sufficiently high redshifts (z_r > 20), the intergalactic medium is heated to 10,000 K and most minihalos never form. On the other hand, if z_r < 20, a significant fraction (>10 percent) of all baryons have already collapsed into minihalos, and are subsequently removed from the halos by photoevaporation as the ionizing background flux builds up. We show that this process can require a significant budget of ionizing photons; exceeding the production by a straightforward extrapolation back in time of known quasar and galaxy populations by a factor of up to 10 and 3, respectively.

We discuss the ultraviolet to near-IR galaxy counts from the deepest imaging surveys, including the northern and southern Hubble Deep Fields. The logarithmic slope of the galaxy number-magnitude relation is flatter than 0.4 in all seven UBVIJHK optical bandpasses at faint magnitudes, i.e. the light from resolved galaxies has converged from the UV to the near-IR. Most of the galaxy contribution to the extragalactic background light (EBL) comes from relatively bright, low-redshift objects (50% at V_AB<21 and 90% at V_AB<25.5). We find a lower limit to the surface brightness of the optical EBL of about 15 nW/m^2/sterad, comparable to the intensity of the far-IR background from COBE data. Diffuse light, lost because of surface brightness selection effects, may add substantially to the EBL.

We discuss implications of a stellar-dominated UV background at high redshifts. The composite spectrum of 29 LBGs evaluated by Steidel etal (2000) at <z>=3.4 can be well fit by a stellar population with ongoing star formation, a Salpeter initial mass function, modest or negligible dust reddening, and no intrinsic HI photoelectric absorption. Fading starbursts in which star formation has ceased for 10^7 yr or more cannot reproduce the observed flux shortward of 1 Ryd.The escape fraction of H-ionizing photons must be close to 100 percent for the observed sample of LBGs. The spectrum of ionizing photons produced by a stellar population with ongoing star formation is similar to that of QSOs between 1 and 3 Ryd, becomes softer between 3 and 4 Ryd and drops sharply shortward of 4 Ryd. A galaxy-dominated UV background appears inconsistent with the observed HeII/HI opacity ratio at z=2.4, but might be able to explain the SiIV/CIV abundances measured at z>3 in absorption spectra. A scenario may be emerging where star-forming galaxies reionized intergalactic hydrogen at $z>6$ and dominate the 1 Ryd metagalactic flux at z>3, with quasi-stellar sources taking over at low redshifts. If the large amplitude of the H-ionizing flux estimated by Steidel etal is correct, hydrodynamical simulations of structure formation in the IGM require a baryon density (to explain the observed Lyalpha opacity in absorption spectra) which is similar to or larger than that favoured by recent CMB experiments, and is inconsistent with standard nucleosynthesis values.

Recent evidence appears to link gamma-ray bursts (GRBs) to star-forming regions in galaxies at cosmological distances. If short-lived massive stars are the progenitors of GRBs, the rate of events per unit cosmological volume should be an unbiased tracer (i.e. unaffected by dust obscuration and surface brightness limits) of the cosmic history of star formation. Here we use realistic estimates for the evolution of the stellar birthrate in galaxies to model the number counts, redshift distribution, and time-delay factors of GRBs. We present luminosity function fits to the BATSE log N-log P relation for different redshift distributions of the bursts. Our results imply about 1-2 GRBs every one million Type II supernovae, and a characteristic `isotropic-equivalent' burst luminosity in the range 3-20 x 1e51 ergs/s (for H_0=65 km/s/Mpc). We compute the rate of multiple imaging of background GRBs due to foreground mass condensations in a Lambda-dominated cold dark matter cosmology, assuming that dark halos approximate singular isothermal spheres on galaxy scales and Navarro-Frenk-White profiles on group/cluster scales, and are distributed in mass according to the Press-Schechter model. We show that the expected sensivity increase of Swift relative to BATSE could result in a few strongly lensed individual bursts detected down to a photon flux of 0.1 phot/s/cm^2 in a 3-year survey. Because of the partial sky coverage, however, it is unlikely that the Swift satellite will observe recurrent events (lensed pairs).

In currently popular cosmological scenarios -- all variants of the cold dark matter (CDM) cosmogony -- some time beyond a redshift of 15, stars within the numerous small halos that condense with virial temperatures ~1e4 K created the first heavy elements; these protogalactic systems, together perhaps with an early population of mini-quasars, generated the ultraviolet radiation and mechanical energy that reheated and reionized the cosmos. The history of the Universe during and soon after these crucial formative stages is recorded in the thermal, ionization, and chemical state of the all-pervading intergalactic medium (IGM), which contains most of the ordinary baryonic material left over from the big bang. Throughout the epoch of structure formation, the IGM becomes clumpy under the influence of gravity, and acts as a source for the gas that gets accreted, cools, and forms stars within subgalactic fragments, and as a sink for the metal enriched material, energy, and radiation which they eject.

Recent observations of high-redshift galaxies and quasars indicate that the hydrogen component of the intergalactic medium (IGM) must have been reionized at some redshift z>6. Prior to complete reionization, sources of ultraviolet radiation will be seen behind intervening gas that is still neutral, and their spectra should show the red damping wing of the Gunn-Peterson trough. While this characteristic feature may, in principle, totally suppress the Lyman-alphaa emission line in the spectra of the first generation of objects in the universe, we show here that the IGM in the vicinity of luminous quasars will be highly photoionized on several Mpc scales due to the source emission of Lyman-continuum photons. If the quasar lifetime is shorter than the expansion and gas recombination timescales, the volume ionized will be proportional to the total number of photons produced above 13.6 eV: the effect of this local photoionization is to greatly reduce the scattering opacity between the redshift of the quasar and the boundary of its HII region. We find that the transmission on the red side of the Lyman-alpha resonance is always greater than 50% for sources radiating a total of 1e69.5 ionizing photons into the IGM. The detection of a strong Lyman-alpha emission line in the spectra of bright QSOs shining for >1e7 yr cannot then be used, by itself, as a constraint on the reionization epoch. The first signs of an object radiating prior to the transition from a neutral to an ionized universe may be best searched for in the spectra of luminous sources with a small escape fraction of Lyman-continuum photons into the IGM, or sources with a short duty cycle.

Distant Type Ia and II supernovae (SNe) can serve as valuable probes of the history of the cosmic expansion and star formation, and provide important information on their progenitor models. At present, however, there are few observational constraints on the abundance of SNe at high redshifts. A major science driver for the Next Generation Space Telescope (NGST) is the study of such very distant supernovae. In this paper we discuss strategies for finding and counting distant SNe by using repeat imaging of super-critical intermediate redshift clusters whose mass distributions are well-constrained via modelling of strongly-lensed features. For a variety of different models for the star formation history and supernova progenitors, we estimate the likelihood of detecting lensed SNe as a function of their redshift. In the case of a survey conducted with HST, we predict a high probability of seeing a supernova in a single return visit with either WFPC-2 or ACS, and a much higher probability of detecting examples with z>1 in the lensed case. Most events would represent magnified SNe II at z~1, and a fraction will be more distant examples. We discuss various ways to classify such events using ground-based infrared photometry. We demonstrate an application of the method using the HST archival data and discuss the case of a possible event found in the rich cluster AC 114 (z=0.31).

In popular cosmological scenarios, some time beyond a redshift of 10, stars within protogalaxies created the first heavy elements; these systems, together perhaps with an early population of quasars, generated the ultraviolet radiation and mechanical energy that reheated and reionized the cosmos. The history of the Universe during and soon after these crucial formative stages is recorded in the all-pervading intergalactic medium (IGM), which contains most of the ordinary baryonic material left over from the big bang. Throughout the epoch of structure formation, the IGM becomes clumpy and acquires peculiar motions under the influence of gravity, and acts as a source for the gas that gets accreted, cools, and forms stars within galaxies, and as a sink for the metal enriched material, energy, and radiation which they eject.

About half a million years after the Big Bang, the ever-fading cosmic blackbody radiation cooled below 3000 K and shifted first into the infrared and then into the radio, and the smooth baryonic plasma that filled the Universe became neutral. The Universe then entered a ``dark age'' which persisted until the first cosmic structures collapsed into gravitationally-bound systems, and evolved into stars, galaxies, and black holes that lit up the Universe again. Some time between redshift of 7 and 15, stars within protogalaxies created the first heavy elements; these systems, together perhaps with an early population of quasars, generated the ultraviolet radiation that reheated and reionized the cosmos. The history of the Universe during and soon after these crucial formative stages is recorded in the all-pervading intergalactic medium (IGM), which is believed to contain most of the ordinary baryonic material left over from the Big Bang. Throughout the epoch of structure formation, the IGM becomes clumpy and acquires peculiar motions under the influence of gravity, and acts as a source for the gas that gets accreted, cools, and forms stars within galaxies, and as a sink for the metal enriched material, energy, and radiation which they eject. Observations of absorption lines in quasar spectra at redshifts up to 5 provide invaluable insight into the chemical composition of the IGM and primordial density fluctuation spectrum of some of the earliest formed cosmological structures, as well as of the ultraviolet background radiation that ionizes them.

Recent observations of gamma-ray bursts (GRBs) have provided growing evidence for collimated outflows and emission, and strengthened the connection between GRBs and supernovae. If massive stars are the progenitors of GRBs, the hard photon pulse will propagate in the pre-burst, dense environment. Circumstellar material will Compton scatter the prompt GRB radiation and give rise to a reflection echo. We calculate luminosities, spectra, and light curves of such Compton echoes in a variety of emission geometries and ambient gas distributions, and show that the delayed hard X-ray flash from a pulse propagating into a red supergiant wind could be detectable by Swift out to z~0.2. Independently of the gamma-ray spectrum of the prompt burst, reflection echoes will typically show a high-energy cutoff between m_ec^2/2 and m_ec^2 because of Compton downscattering. At fixed burst energy per steradian, the luminosity of the reflected echo is proportional to the beaming solid angle, Omega_b, of the prompt pulse, while the number of bright echoes detectable in the sky above a fixed limiting flux increases as Omega_b^{1/2}, i.e. it is smaller in the case of more collimated jets. The lack of an X-ray echo at one month delay from the explosion poses severe constraints on the possible existence of a lateral GRB jet in SN 1987A. The late r-band afterglow observed in GRB990123 is fainter than the optical echo expected in a dense red supergiant environment from a isotropic prompt optical flash. Significant MeV delayed emission may be produced through the bulk Compton (or Compton drag) effect resulting from the interaction of the decelerating fireball with the scattered X-ray radiation.

We consider the effects of rapid pair creation by an intense pulse of gamma-rays propagating ahead of a relativistic shock. Side-scattered photons colliding with the main gamma-ray beam amplify the density of scattering charges. The acceleration rate of the pair-loaded medium is calculated, and its limiting bulk Lorentz factor related to the spectrum and compactness of the photon source. One obtains, as a result, a definite prediction for the relative inertia in baryons and pairs. The deceleration of a relativistic shock in the moving medium, and the resulting synchrotron emissivity, are compared with existing calculations for a static medium. The radiative efficiency is increased dramatically by pair loading. When the initial ambient density exceeds a critical value, the scattering depth traversed by the main gamma-ray pulse rises above unity, and the pulse is broadened. These considerations place significant constraints on burst progenitors: a pre-burst mass loss rate exceeding 10^{-5} M_\odot per year is difficult to reconcile with individual pulses narrower than 10 s, unless the radiative efficiency is low. An anisotropic gamma-ray flux (on an angular scale \Gamma^{-1} or larger) drives a large velocity shear that greatly increases the energy in the seed magnetic field forward of the propagating shock.

We consider the radiative acceleration to relativistic bulk velocities of a cold, optically thin plasma which is exposed to an external source of gamma-rays. The flow is driven by radiative momentum input to the gas, the accelerating force being due to Compton scattering in the relativistic Klein-Nishina limit. The bulk Lorentz factor of the plasma, Gamma, derived as a function of distance from the radiating source, is compared with the corresponding result in the Thomson limit. Depending on the geometry and spectrum of the radiation field, we find that particles are accelerated to the asymptotic Lorentz factor at infinity much more rapidly in the relativistic regime; and the radiation drag is reduced as blueshifted, aberrated photons experience a decreased relativistic cross section and scatter preferentially in the forward direction. The random energy imparted to the plasma by gamma-rays can be converted into bulk motion if the hot particles execute many Larmor orbits before cooling. This `Compton afterburn' may be a supplementary source of momentum if energetic leptons are injected by pair creation, but can be neglected in the case of pure Klein-Nishina scattering. Compton drag by side-scattered radiation is shown to be more important in limiting the bulk Lorentz factor than the finite inertia of the accelerating medium. The processes discussed here may be relevant to a variety of astrophysical situations where luminous compact sources of hard X- and gamma-ray photons are observed, including active galactic nuclei, galactic black hole candidates, and gamma-ray bursts.

We assess the constraints imposed by the observed extragalactic background light (EBL) on the cosmic history of star formation and the stellar mass density today. The logarithmic slope of the galaxy number-magnitude relation from the Southern Hubble Deep Field} imaging survey is flatter than 0.4 in all seven UBVIJHK optical bandpasses, i.e. the light from resolved galaxies has converged from the UV to the near-IR. We find a lower limit to the surface brightness of the optical extragalactic sky of about 15 nW/m^2/sr, comparable to the intensity of the far-IR background from COBE data. Assuming a Salpeter initial mass function with a lower cutoff consistent with observations of M subdwarf disk stars, we set a lower limit of Omega_g+s h^2>0.0013 I_50 to the visible (processed gas + stars) mass density required to generate an EBL at a level of 50 I_50 nW/m^2/sr; our `best-guess' value is Omega_g+s h^2=0.003 I_50. Motivated by the recent microlensing results of the MACHO collaboration, we consider the possibility that massive dark halos around spiral galaxies are composed of faint white dwarfs, and show that only a small fraction (<5%) of the nucleosynthetic baryons can be locked in the remnants of intermediate-mass stars forming at z_F<5, as the bright early phases of such halos would otherwise overproduce the observed EBL.

Ongoing searches for supernovae (SNe) at cosmological distances have recently started to provide large numbers of events with measured redshifts and apparent brightnesses. Compared to quasars or galaxies, Type Ia SNe represent a population of sources with well-known intrinsic properties, and could be used to detect gravitational lensing even in the absence of multiple or highly distorted images. We investigate the lensing effect of background SNe due to mass condensations in three popular cold dark matter cosmologies (LambdaCDM, OCDM, SCDM), and compute lensing frequencies, rates of SN explosions, and distributions of arrival time differences and image separations. If dark halos approximate singular isothermal spheres on galaxy scales and NFW profiles on group/cluster scales, and are distributed in mass according to the Press-Schechter theory, then about one every 12 SNe at z~1 will be magnified by Delta m>0.1 mag (SCDM). The detection rate of SN Ia with magnification Delta m>0.3 is estimated to be of order a few events per year per square degree at maximum B-light and I<25, a hundred time smaller than the total rate expected at these magnitude levels. In the field, events magnified by more than 0.75 mag are 7 times less frequent: about one fifth of them gives rise to observable multiple images. Because of the flat K-correction and wide luminosity function, Type II SNe dominate the number counts at I>25 and have the largest fraction of lensed objects. The optimal survey sensitivity for Type Ia's magnified by Delta m>0.75 mag is I=23. Magnification bias increases their incidence by a factor of 50 in samples with I<22, dropping to a factor of 3.5 at 24 mag. At faint magnitudes the enhancement is larger for SN II.

I review the constraints imposed by the observed extragalactic background light (EBL) on the history of the stellar birthrate in galaxies. At faint magnitudes, the logarithmic slope of the galaxy counts is flatter than 0.4 in all seven UBVIJHK optical bandpasses of the Hubble Deep Field-South imaging survey. The integration of the number counts provides a lower limit to the surface brightness of the optical extragalactic sky of 15 nW/m2/sr, comparable to the intensity of the far-IR background from COBE data. If the initial mass function has a Salpeter slope with a lower mass cutoff consistent with observations of M subdwarf disk stars, a lower limit of Omega_*>0.005 (at Hubble constant 50 km/s/Mpc) is derived for the visible (processed gas + stars) mass density needed to generate an extragalactic background light (EBL) at a level of 50 nW/m2/sr. The current `best-guess' estimate to Omega_* is 0.012, about 16% of the nucleosynthetic baryon density. The contribution of quasar activity to the observed EBL is unlikely to exceed 20%.

The emission of 21-cm radiation from a neutral intergalactic medium (IGM) at high redshift is discussed in connection with the thermal and ionization history of the universe. The physical mechanisms that make such radiation detectable against the cosmic microwave background include Ly_alpha coupling of the hydrogen spin temperature to the kinetic temperature of the gas and preheating of the IGM by the first generation of stars and quasars. Three different signatures are investigated in detail: (a) the fluctuations in the redshifted 21-cm emission induced by the gas density inhomogeneities that develop at early times in cold dark matter (CDM) dominated cosmologies; (b) the sharp absorption feature in the radio sky due to the rapid rise of the Ly_alpha continuum background that marks the birth of the first UV sources in the universe; and (c) the 21-cm emission and absorption shells that are generated on several Mpc scales around the first bright quasars. Future radio observations with projected facilities like the Giant Metrewave Radio Telescope and the Square Kilometer Array may shed light on the power spectrum of density fluctuations at z>5, and map the end of the "dark ages", i.e. the transition from the post-recombination universe to one populated with radiation sources.

We show how the investigation of the redshifted 21-cm radiation can give insight into the development of structures in the early universe (at redshifts z>5). In particular we investigate: the epoch of the first light; the fluctuations in the redshifted 21-cm emission induced by the density inhomogeneities in CDM dominated universes; the emission and absorption shells that are generated around the first bright quasars. Such features can be observed with the next generation radio facilities.

High-resolution hydrodynamics simulations of the Ly-alpha forest in cold dark matter dominated cosmologies appear to predict line widths that are substantially narrower than those observed. Here we point out that Compton heating of the intergalactic gas by the hard X-ray background (XRB), an effect neglected in all previous investigations, may help to resolve this discrepancy. The rate of gain in thermal energy by Compton scattering will dominate over the energy input from hydrogen photoionization if the XRB energy density is 0.2x/<epsilon> times higher than the energy density of the UV background at a given epoch, where x is the hydrogen neutral fraction in units of 1e-6 and <epsilon> is the mean X-ray photon energy in units of m_ec^2. The numerical integration of the time-dependent rate equations shows that the intergalactic medium approaches a temperature of about 1.5e4 K at z>3 in popular models for the redshift evolution of the extragalactic background radiation. The importance of Compton heating can be tested experimentally by measuring the Ly-alpha line-width distribution as a function of redshift, thus the Lyman-alpha forest may provide a useful probe of the evolution of the XRB at high redshifts.

There has been remarkable progress recently in both observational and theoretical studies of galaxy formation and evolution. Largely due to a combination of deep Hubble Space Telescope imaging, Keck spectroscopy, and COBE far-IR background measurements, new constraints have emerged on the emission history of the galaxy population as a whole. While the bulk of the stars present today appears to have formed relatively recently, the existence of a decline in the star-formation density above z=2 remains uncertain. The study of the transition from the cosmic `dark age' to an ionized universe populated with luminous sources can shed new light on the star formation activity at high redshifts, and promises answers to some fundamental questions on the formation of cosmic structures. If stellar sources are responsible for photoionizing the intergalactic medium at z=5, the rate of star formation at this epoch must be comparable or greater than the one inferred from optical observations of galaxies at z=3. A population of quasars at z<2 could make a significant contribution to the extragalactic background light if dust-obscured accretion onto supermassive black holes is an efficient process.

The reionization of the Universe is expected to leave a signal in the form of a sharp step in the spectrum of the sky. If reionization occurs at 5 < z < 20, a feature should appear in the radio sky at 70 - 240 MHz due to redshifted HI 21-cm line emission, accompanied by another feature in the optical/near-IR at 0.7 - 2.6 micron due to hydrogen recombination radiation. The expected amplitude is well above fundamental detection limits, and the sharpness of the feature may make it distinguishable from variations due to terrestrial, galactic and extragalactic foregrounds. Because this is essentially a continuum measurement of a signal which occurs over the whole sky, relatively small telescopes may suffice for detection in the radio. In the optical/near-IR, a space telescope is needed with the lowest possible background conditions, since the experiment will be severely background-limited.

I review recent observational and theoretical progress in our understanding of the cosmic evolution of luminous sources. Largely due to a combination of deep HST imaging, Keck spectroscopy, and COBE far-IR background measurements, new constraints have emerged on the emission history of the galaxy population as a whole. Barring large systematic effects, the global ultraviolet, optical, near- and far-IR photometric properties of galaxies as a function of cosmic time cannot be reproduced by a simple stellar evolution model defined by a constant (comoving) star-formation density and a universal (Salpeter) IMF, and require instead a substantial increase in the stellar birthrate with lookback time. While the bulk of the stars present today appears to have formed relatively recently, the existence of a decline in the star-formation density above z=2 remains uncertain. The history of the transition from the cosmic `dark age' to a ionized universe populated with luminous sources can constrain the star formation activity at high redshifts. If stellar sources are responsible for photoionizing the intergalactic medium at z=5, the rate of star formation at this epoch must be comparable or greater than the one inferred from optical observations of galaxies at z=3. A population of dusty, Type II AGNs at z<2 could make a significant contribution to the FIR background if the accretion efficiency is of order 10%.

Many questions in physical cosmology regarding the thermal and ionization history of the intergalactic medium are now successfully studied with the help of cosmological hydrodynamical simulations. Here we present a numerical method that solves the radiative transfer around point sources within a three dimensional cartesian grid. The method is energy conserving independently of resolution: this ensures the correct propagation speeds of ionization fronts. We describe the details of the algorithm, and compute as first numerical application the ionized region surrounding a mini-quasar in a cosmological density field at z=7.

The combined use of the ESO Very Large Telescope (VLT) UT1 Science Verification (SV) images and of the Hubble Space Telescope (HST) Hubble Deep Field South observations allows us to strengthen the identification as a candidate elliptical galaxy of the Extremely Red Object HDFS 223251-603910 previously identified by us on the basis of NICMOS and Cerro Tololo Interamerican Observatory imaging. The photometry presented here includes VLT data in U, B, V, R, I, a STIS unfiltered image, NICMOS J, H, and K band data, thus combining the 16.5 hours of VLT SV exposures with 101 hours of HST observing. The object is detected in all images except the VLT U band and is one of the reddest known with B-K=9.7+-0.5. We consider a wide range of models with different ages, metallicities, star formation histories and dust content, and conclude that the observed spectral energy distribution agrees best with that of an old elliptical galaxy at redshift just below 2. Alternative possibilities are discussed in light of their likelihood and of the perspective of spectroscopic confirmation.

I review recent observational and theoretical progress in our understanding of the cosmic evolution of luminous sources. Through a combination of deep imaging, Keck spectroscopy, and COBE background measurements, important constraints have emerged on the emission history of the galaxy population as a whole. A simple stellar evolution model, defined by a star-formation density that rises from z=0 to z=1.5, a universal Salpeter IMF, and a moderate amount of dust with A_V=0.23 mag (A_1500=1.2 mag), is able to account for most of the optical-FIR extragalactic background light, and reproduces the global ultraviolet, optical, and near-IR photometric properties of the universe. By contrast, a star-formation density that stayed roughly constant at all epochs appears to overproduce the local K-band luminosity density. While the bulk of the stars present today formed relatively recently, the existence of a decline in the star-formation density above z=2 remains uncertain. If stellar sources are responsible for photoionizing the intergalactic medium at z=5, the rate of star formation at this epoch must be comparable or greater than the one inferred from optical observations of galaxies at z=3. A population of dusty AGNs at z>2 could make a significant contribution to the FIR background if the accretion efficiency is close to 10%.

The history of the transition from a neutral intergalactic medium (IGM) to one that is almost fully ionized can reveal the character of cosmological ionizing sources. We study the evolution of the volume filling factor of HII and HeIII regions in a clumpy IGM, and discuss the implications for rival reionization scenarios of the rapid decline in the space density of radio-loud quasars and of the large population of star-forming galaxies recently observed at z>3. The hydrogen component in a highly inhomogeneous universe is completely reionized when the number of photons emitted above 1 ryd in one recombination time equals the mean number of hydrogen atoms. If stellar sources are responsible for keeping the IGM ionized at z=5, the rate of star formation at this epoch must be comparable or greater than the one inferred from optical observations of galaxies at z=3, and the mean metallicity per baryon in the universe of order 1/500 solar. An early generation of stars in dark matter halos with circular velocities v_circ=50 km/s, possibly one of the main source of UV photons at high-z, could be detectable with the Next Generation Space Telescope. Models in which the quasar emissivity declines rapidly at z>3 predict a late HeII reionization epoch, a feature that could explain the recent detection of patchy HeII Lyman-alpha at z=2.9 by Reimers et al. (1997) and the abrupt change observed by Songaila (1998) at about the same epoch of the SiIV/CIV ratio, but appear unable to provide the required number of hydrogen-ionizing photons at z=5.

We report the serendipitous discovery of an extremely red object in the Hubble Deep Field South (HDFS) Test NICMOS (Near Infrared Camera and Multi Object Spectrograph) field of view. The object is resolved in the NICMOS image and has a light profile very well described by an r^{1/4} law with effective radius r_e=0."20+-0."05 and H_{AB}=21.7+-0.1 magnitudes. In contrast, the galaxy is undetected in the R and I band ground based images taken at the CTIO 4 m Blanco Telescope, giving a lower limit to the color of (R-H)_{AB}>3.9 and (I-H)_{AB}>3.5 at the 95% confidence level. The colors of a range of synthetic galactic spectra are computed, showing that the object is likely to be an ``old'' elliptical galaxy at redshift z>~1.7. Alternatively the colors can be reproduced by an ``old'' elliptical galaxy at somewhat lower redshift (z>~1) with significant amount of dust, or by a younger galaxy at higher redshift. This object represents a very interesting target for future VLT observations.

Through a combination of deep wide-field imaging and multi-object spectroscopy, the Next Generation Space Telescope will be able to chart with unprecedented accuracy the evolution of cosmic structures after the `dark ages' (z<5), when galaxies are thought to assemble and form the bulk of their stars. In this talk I will discuss how NGST observations might determine the history of the stellar birthrate in the universe and reveal the character of cosmological ionizing sources at high redshifts.

The history of the transition from a neutral intergalactic medium (IGM) to one that is almost fully ionized can reveal the character of cosmological ionizing sources and set important constraints on the stellar birthrate at high redshifts. The hydrogen component in a highly inhomogeneous universe is completely reionized when the number of photons emitted above 1 ryd in one recombination time equals the mean number of hydrogen atoms. If stellar sources are responsible for photoionizing the IGM at z=5, the rate of star formation at this epoch must be comparable or greater than the one inferred from optical observations of galaxies at z=3, and the mean metallicity per baryon in the universe greater than 1/500 solar. In hierarchical clustering scenarios, high-z dwarfs (i.e. an early generation of stars in dark matter halos with circular velocities of order 50 km/s) are expected to be one of the main source of UV photons and heavy elements at early epochs. They would be very numerous, >0.2 arcsec^{-2}, and faint, I_AB> 29.5 mag: their detection may have to wait for the Next Generation Space Telescope.

High resolution, deep imaging surveys are instrumental in setting constraints on semi-analytical structure formation models in Cold Dark Matter (CDM) cosmologies. We show here that the lack of unresolved B-band ``dropouts'' with V > 25 mag in the Hubble Deep Field (HDF) appears to be inconsistent with the expected number of quasars if massive black holes form with a constant universal efficiency in all CDM halos. To reconcile the models with the data, a mechanism is needed that suppresses the formation of quasars in halos with circular velocities v(circ) < 50-75 km/s. This feedback naturally arises due to the photoionization heating of the gas by the UV background. We consider several alternative effects that would help reduce the quasar number counts, and find that these can not alone account for the observed lack of detections. If reddening by dust can be neglected at early epochs, consistency with the optical data also requires that the luminous extent of dwarf galaxies at high redshifts be larger than a few percent of the virial radii of their dark matter halos, in order not to overpredict the number of point-like B-band dropouts. Future deep observations in the J and H bands with NICMOS might reveal several z > 5 objects per field or provide even stronger constraints on the models than existing B, V, and I data.

The history of the transition from a neutral intergalactic medium to one that is almost fully ionized can reveal the character of cosmological ionizing sources. In this talk I will discuss the implications for rival reionization scenarios of the rapid decline observed in the space density of quasars and star-forming galaxies at redshifts z>3. The hydrogen component in a highly inhomogeneous universe is completely reionized when the number of ionizing photons emitted in one recombination time equals the mean number of hydrogen atoms. At z=5, the local character of the UV metagalactic flux allows one to define a critical emission rate of hydrogen-ionizing photons per unit comoving volume. Models based on photoionization by bright QSOs and/or young galaxies with star formation rates in excess of 0.3-1 M_sun/yr appear to fail to provide the required number of hydrogen-ionizing photons at these redshifts by large factors. If stellar sources are responsible for keeping the universe ionized at z=5, the rate of star formation per unit comoving volume at this epoch must be comparable or greater than observed at z=3.

Ongoing searches for supernovae (SNe) at cosmological distances have recently started to provide a link between SN Ia statistics and galaxy evolution. We use recent estimates of the global history of star formation to compute the theoretical Type Ia and Type II SN rates as a function of cosmic time from the present epoch to high redshifts. We show that accurate measurements of the frequency of SN events in the range 0<z<1 will be valuable probes of the nature of Type Ia progenitors and the evolution of the stellar birthrate in the universe. The Next Generation Space Telescope should detect of order 20 Type II SNe per 4'x 4' field per year in the interval 1<z<4.

We discuss the deep galaxy counts from the Hubble Deep Field (HDF) imaging survey. At faint magnitudes, the slope of the differential number-magnitude relation is flatter than 0.2 in all four HDF bandpasses. In the ultraviolet, a fluctuation analysis shows that the flattening observed below U=26 mag is not due to incompleteness and is more pronounced than in the other bands, consistent with the idea that a redshift limit has been reached in the galaxy distribution. A reddening trend of about 0.5 magnitude is observed at faint fluxes in the colour-magnitude diagram, (U-V)_{eff} vs. V. We interpret these results as the effect of intergalactic attenuation on distant galaxies. At flux levels of AB=27 mag and in agreement with the fluctuation analysis and the colour-magnitude relation, about 7% of the sources in U, 30% in B and 35% in V are Lyman-break ``dropouts'', i.e. candidate star-forming galaxies at z>2. By integrating the number counts to the limits of the HDF survey we find that the mean surface brightness of the extragalactic sky is dominated by galaxies that are relatively bright and are known to have <z>=0.6. To AB=29 mag, the integrated light from resolved galaxies in the I-band is 2.1\pm 0.4 \times 10^{-20} erg/cm^2/s/Hz/sr, and its spectrum is well described by a broken power-law. We discuss the predictions for the counts, colours, and luminosity densities from standard low-q_0 pure-luminosity-evolution models without dust obscuration, and find that they are unable to reproduce all the observed properties of faint field galaxies.

Ongoing searches for supernovae (SNe) at cosmological distances have recently started to provide a link between SN Ia statistics and galaxy evolution. In this talk I will use recent estimates of the global history of star formation to compute the theoretical Type Ia and Type II SN rates as a function of cosmic time. I will compare the predicted values with the rates observed in the range 0<z<0.4, and show how accurate measurements of the frequency of SN events at intermediate redshifts are valuable probes of the evolution of the stellar birthrate in the universe and the nature of Type Ia progenitors. The Next Generation Space Telescope should detect of order 20 Type II SNe per field per year in the interval 1<z<4.

I review a technique for interpreting faint galaxy data which traces the evolution with cosmic time of the galaxy luminosity density, as determined from several deep spectroscopic samples and the HDF. The method relies on the rest frame UV and near-IR continua of galaxies as indicators, for a given IMF and dust content, of their instantaneous SFR and total stellar mass, and offers the prospect of addressing in a coherent framework an important set of subjects: cosmic star formation history, dust in primeval galaxies, shape of the IMF, stellar mass-to-light ratios of present-day galaxies, extragalactic background light, Type II supernovae and heavy element enrichment history of the universe. The global spectrophotometric properties of field galaxies are well fit by a simple stellar evolution model, defined by a time-dependent SFR per unit comoving volume, a universal IMF which is relatively rich in massive stars, and a modest amount of dust reddening. The model is able to account for the entire background light recorded in the galaxy counts down to the very faint magnitude levels probed by the HDF, and produces visible mass-to-light ratios at the present epoch which are consistent with the values observed in nearby galaxies of various morphological types. The bulk (>60%) of the stars present today formed relatively recently (z<1.5), consistently with the expectations from a broad class of hierarchical clustering cosmologies, and in good agreement with the low level of metal enrichment observed at high redshifts in the damped Lyman-alpha systems. Throughout this review I emphasize how the poorly constrained amount of starlight that was absorbed by dust and reradiated in the far-IR at early epochs represents one of the biggest uncertainties in our understanding of the evolution of luminous matter in the universe.

We develop a method for interpreting faint galaxy data which focuses on the integrated light radiated from the galaxy population as a whole. The emission history of the universe at ultraviolet, optical, and near-infrared wavelengths is modeled from the present epoch to z~4 by tracing the evolution with cosmic time of the galaxy luminosity density, as determined from several deep spectroscopic samples and the Hubble Deep Field (HDF) imaging survey. The global spectrophotometric properties of field galaxies can be well fit by a simple stellar evolution model, defined by a time-dependent star formation rate (SFR) per unit comoving volume and a universal IMF extending from 0.1 to 125 M_sun. In the best-fit models, the global SFR rises sharply, by about an order of magnitude, from a redshift of zero to a peak value at z~1.5, to fall again at higher redshifts. The models are able to account for the entire background light recorded in the galaxy counts down to the very faint magnitude levels probed by the HDF. Since only 20% of the current stellar content of galaxies is produced at z>2, a rather low cosmic metallicity is expected at these early times, in good agreement with the observed enrichment history of the damped Lyman-\alpha systems. A ``monolithic collapse'' model, where half of the present-day stars formed at z>2.5 and were shrouded by dust, can be made consistent with the global history of light, but overpredicts the metal mass density at high redshifts as sampled by QSO absorbers.

I review some recent progress made in our understanding of galaxy evolution and the cosmic history of star formation. The Hubble Deep Field (HDF) imaging survey has achieved the sensitivity to capture the bulk of the extragalactic background light from discrete sources. No evidence is found in the optical number-magnitude relation down to AB=29 mag for a large amount of star formation at high redshifts. The emission history of the universe at ultraviolet, optical, and near-infrared wavelengths can be modeled from the present epoch to z~4 by tracing the evolution with cosmic time of the galaxy luminosity density, as determined from several deep spectroscopic samples and the HDF. The global spectrophotometric properties of field galaxies are well fitted by a simple stellar evolution model, defined by a time-dependent star formation rate (SFR) per unit comoving volume and a universal initial mass function which is relatively rich in massive stars. The SFR density is found to rise sharply, by about an order of magnitude, from a redshift of zero to a peak value at z~1.5, to fall again by a factor of 2 (4) out of a redshift of 3 (4). Since only 10% of the current stellar content of galaxies is produced at z>2.5, a rather low cosmic metallicity is predicted at these early times, in good agreement with the observed enrichment history of the damped Lyman-\alpha systems. The biggest uncertainty is represented by the poorly constrained amount of starlight that was absorbed by dust and reradiated in the IR at early epochs.

We use new studies of the cosmic evolution of star-forming galaxies to estimate the production rate of ionizing photons from hot, massive stars at low and intermediate redshifts. The luminosity function of blue galaxies in the Canada-France Redshift Survey shows appreciable evolution in the redshift interval z=0-1.3, and generates a background intensity at 1 ryd of J_L~ 1.3 x 10^{-21} f_{esc} ergs cm^{-2} s^{-1} Hz^{-1} sr^{-1} at z~0.5, where f_esc is the unknown fraction of stellar Lyman-continuum photons which can escape into the intergalactic space, and we have assumed that the absorption is picket fence-type. We argue that recent upper limits on the H-alpha surface brightness of nearby intergalactic clouds constrain this fraction to be <~ 20%. The background ionizing flux from galaxies can exceed the QSO contribution at z~ 0.5 if f_{esc}>~ 6%. We show that, in the general framework of a diffuse background dominated by QSOs and/or star-forming galaxies, the cosmological baryon density associated with photoionized, optically thin gas decreases rapidly with cosmic time. The results of a recent Hubble Space Telescope survey of OVI absorption lines in QSO spectra suggest that most of this evolution may be due to the bulk heating and collisional ionization of the intergalactic medium by supernova events in young galaxy halos.

I review some recent progress made in our understanding of galaxy evolution and the cosmic history of star formation. Like bookends, the results obtained from deep ground-based spectroscopy and from the Hubble Deep Field imaging survey put brackets around the intermediate redshift interval, $1<z<2$, where starbirth probably peaked at a rate 10 times higher than today. The steady decline observed since $z\sim 1$ is largely associated with late-type galaxies. At $z\gtrsim 2.5$, the Lyman-break selected objects may represent the precursors of present-day spheroids, but appear, on average, quite underluminous relative to the expectations of the standard early-and-rapidly forming picture for spheroidal systems. The observed ultraviolet light density accounts for the bulk of the metals seen today in ``normal'' massive galaxies.

The detection of a HeII absorption trough in the spectra of three high-redshift quasars provides new constraints on the spectral shape of the UV extragalactic background. The relative strengths of the observed flux decrements at the rest-frame wavelengths of 1216 and 304 A require a relatively soft radiation field at 4 ryd compared to 1 ryd. If the ionizing metagalactic flux is dominated by the integrated light from QSOs, and the HeII reionization of the universe was completed well before z~3, then the UV/soft X-ray spectrum of individual quasars at these epochs must, once the cosmological ``filtering'' through material along the line of sight is taken into account, satisfy similar constraints on the average. We model the propagation of AGN-like ionizing radiation through the intergalactic medium using CUBA, a numerical code developed in our prior work. We show that, in order to explain the reported HeII absorption, any thermal component responsible for the ``soft X-ray excess'' observed in the spectra of AGNs at z<0.3 must, at z~3, have a typical temperature >80 eV, and a luminosity not exceeding 20% of that of the ``UV bump''.

We investigate the 21-cm signature that may arise from the intergalactic medium (IGM) prior to the epoch of full reionization (z>5). In scenarios in which the IGM is reionized by discrete sources of photoionizing radiation, the neutral gas which has not yet been engulfed by an H II region may easily be preheated to temperatures well above that of the cosmic background radiation (CBR), rendering the IGM invisible in absorption against the CBR. We identify three possible preheating mechanisms: (1) photoelectric heating by soft X-rays from QSOs, (2) photoelectric heating by soft X-rays from early galactic halos, and (3) resonant scattering of the continuum UV radiation from an early generation of stars. In the presence of a sufficiently strong ambient flux of Lyman-alpha photons, the hyperfine transition in the warmed H I will be excited. A beam differencing experiment would detect a patchwork of emission, both in frequency and in angle across the sky. This patchwork could serve as a valuable tool for understanding the epoch, nature, and sources of the reionization of the universe, and their implications for cosmology. We demonstrate that isolated QSOs will produce detectable signals at meter wavelengths within their ``spheres of influence'' over which they warm the IGM. As a result of the redshifted 21-cm radiation emitted by warm H I bubbles, the spectrum of the radio extragalactic background will display frequency structure with velocity widths up to 10,000 km/s. Broad beam observations would reveal corresponding angular fluctuations in the sky intensity with (delta T/T)=1E-3 on degree-scales. Radio measurements near 235 and 150 MHz, as will be possible in the near future using the Giant Metrewave Radio Telescope, may provide the first detection of a neutral IGM at 5<z<10.

The Lyman decrement associated with the cumulative effect of HI in QSO absorption systems along the line of sight provides a distinctive feature for identifying galaxies at z>2.5. The Hubble Deep Field (HDF) observations offer the opportunity to exploit the ubiquitous effect of intergalactic absorption and obtain useful statistical constraints on the redshift distribution of galaxies considerably fainter than current spectroscopic limits. We model the HI cosmic opacity as a function of redshift, including scattering in resonant lines of the Lyman series and Lyman-continuum absorption, and use stellar population synthesis models with a wide variety of ages, metallicities, dust contents, and redshifts, to derive color selection criteria that provide a robust separation between high redshift and low redshift galaxies. From the HDF images we construct a sample of star-forming galaxies at 2<z<4.5. While none of the more than 60 objects in the HDF having known Keck/LRIS spectroscopic redshifts in the range 0<z<1.4 is found to contaminate our high-redshift sample, our color criteria are able to efficiently select the 2.6<z<3.2 galaxies identified by Steidel et al. (1996b). The ultraviolet (and blue) dropout technique opens up the possibility of investigating cosmic star and element formation in the early universe.

The recent discovery by Cowie \etal and Tytler \etal of metals in the \Lya clouds shows that the intergalactic medium (IGM) at high redshift is contaminated by the products of stars, and suggests that ionizing photons from massive star formation may be a significant contributor to the UV background radiation at early epochs. We assess the validity of the stellar photoionization hypothesis. Based on recent computations of metal yields and O-star Lyman continuum (LyC) fluxes, we find that 0.2\% of the rest-mass energy of the metals produced is radiated as LyC. By modeling the transfer of ionizing radiation through the IGM and the rate of chemical enrichment, we demonstrate that the background intensity of photons at 1 Ryd that accompanies the production of metals in the \Lya forest clouds may be significant, approaching $0.5\times 10^{-21}\uvunits$ at $z\approx 3$ if the LyC escape fraction is $\gta 0.25$. Together with quasars, massive stars could then, in principle, provide the hydrogen and helium LyC photons required to ionize the universe at high redshifts. We propose that observations of the \HeII Gunn-Peterson effect and of the metal ionization states of the \Lya forest and Lyman-limit absorbers should show the signature of a stellar spectrum. We also note that the stellar photoionization model fails if a large fraction of the UV radiation emitted from stars cannot escape into the IGM, as suggested by the recent {\it Hopkins Ultraviolet Telescope} observations by Leitherer \etal of low-redshift starburst galaxies, or if most of the metals observed at $z\approx 3$ were produced at much earlier epochs.

We present a detailed calculation of the propagation of AGN-like ionizing radiation through the intergalactic space. We model the ionization state of absorbing clouds, and show that the universe will be more opaque above $4\,$Ryd than previously estimated. Singly ionized helium in \Lya forest clouds and Lyman-limit systems is found to be very efficient in reprocessing soft X-ray, helium-ionizing photons into ultraviolet, hydrogen-ionizing ones. We demonstrate that a significant fraction of the absorbed primary photons (emitted, e.g., by quasar sources) will be reradiated by the photoionized gas through \Lya line emission, two-photon continuum, and recombination continuum radiation. In the light of new data and recent studies, we also reassess the contribution of the QSOs observed in optical surveys to the UV extragalactic background, and find that the stochastic reprocessing of quasar Lyman continuum radiation by hydrogen and helium along the line of sight will significantly affect the amplitude, spectral shape, and fluctuations properties of the metagalactic flux. In a scenario in which QSOs are the primary source of ionizing photons in the universe, the integrated \HI \Lya emission at $z=0$ from photoionized \Lya clouds and Lyman-limit systems is found to be at a level of less than 5\% of current observational limits on the far-UV extragalactic radiation flux. We show that $J_{912}$ increases from $\approx 10^{-23} \uvunits$ at the present epoch to $\approx 5\times 10^{-22}\uvunits$ at $z=2.5$. The attenuated direct quasar emission plus recombination radiation from intergalactic gas appears to provide enough hydrogen-ionizing photons to satisfy the proximity effect at large redshift. The \HeII/\HI ratio in the diffuse intergalactic medium and the \Lya clouds increases from $\approx25$ at